Classical alternating direction (AD) methods for parabolic equations, based on some standard implicit time stepping procedure such as Crank-Nicolson, can have errors associated with the AD perturbations that are much larger than the errors associated with the underlying time stepping procedure . We plan to show that minor modifications in the AD procedures can virtually eliminate the perturbation errors at an minor additional computational cost. A mixed finite element method is applied in the spactial variables. Similar to the finite difference and finite element methods in spactial variables, we plan to have the same accuracy in time. A convergence analysis can also be shown .

The present study developes a mixed finite element formulation for the analysis of piezoelectric circular and annular plates. This formulation combines the conventional displacement-electric potential type variational principle and the piezoelectric Reissner`s principle with a weighting factor which represents ratio of weights imposed on the above two variational principles, and which can be adjusted by comparing with experiment results. With this formulation, stresses and electric displacements, like displacements and electric potential, are primary variables and are continuous across elements and element interfaces. Also, all displacement, stress, electric displacement, and electric potential boundary conditions can be easily and exactly imposed. Static deformations and vibration frequencies of some typical piezoelectric circular and annular plates are then obtained with the present approach and are compared with those by other methods. Based on experiment results in the literature, it is found that better results could be obtained in general by the present mixed finite element formulation than the others when 1 is chosen as the weighting factor.

In my dissertation, I developed mixed hp-finite element methods for linear elasticity with weakly imposed symmetry, which is based on Arnold-Falk-Winther's stable mixed finite elements. I have proved the h-stability of my method for meshes with arbitrary variable orders. In order to show the h-stability, I need an upper limit of the highest order of meshes, which can be an arbitrary nonnegative integer. / text

In the proposed project¡Athe three - dimensional vibration of circular and annular plates is analyzed by a mixed finite element¡C Stresses¡Aas well as displacements¡A are primary variables in the mixed finite element formulation¡Atherefore¡Aall the stress and displacement boundary conditions can be imposed exactly¡CMeanwhile¡Athe proposed finite element is a modification of axisymmetric finite element which is based on three ¡V dimensional elasticity¡Aso general results of both axisymmetric and unaxisymmetric vibration of circular and annular plates can be obtained¡C Results of the present project will be compared to those by conventional displacement ¡V type finite element¡ARitz method and series method to show the difference among these theories¡CEspecially¡Athe effect of satisfying the stress boundary conditions on the unaxisymmetric vibration analyses can be demonstrated¡Awhich is not available in the literature up to date¡C

Swelling of a polymer gel is a kinetic process coupling mass transport and mechanical deformation. A comparison between a nonlinear theory for polymer gels and the classical theory of linear poroelasticity is presented. It is shown that the two theories are consistent within the linear regime under the condition of a small perturbation from an isotropically swollen state of the gel. The relationships between the material properties in the linear theory and those in the nonlinear theory are established by a linearization procedure. Both linear and nonlinear solutions are presented for swelling kinetics of substrate-constrained and freestanding hydrogel layers. A new procedure is suggested to fit the experimental data with the nonlinear theory. A nonlinear, transient finite element formulation is presented for initial boundary value problems associated with swelling and deformation of hydrogels, based on nonlinear continuum theories for hydrogels with compressible and incompressible constituents. The incompressible instantaneous response of the aggregate imposes a constraint to the finite element discretization in order to satisfy the LBB condition for numerical stability of the mixed method. Three problems of practical interests are considered: constrained swelling, flat-punch indentation, and fracture of hydrogels. Constrained swelling may lead to instantaneous surface instability. Indentation relaxation of hydrogels is simulated beyond the linear regime under plane strain conditions, and is compared with two elastic limits for the instantaneous and equilibrium states. The effects of Poisson’s ratio and loading rate are discussed. On the study of hydrogel fracture, a method for calculating the transient energy release rate for crack growth in hydrogels, based on a modified path-independent J-integral, is presented. The transient energy release rate takes into account the energy dissipation due to diffusion. Numerical simulations are performed for a stationary center crack loaded in mode I, with both immersed and non-immersed chemical boundary conditions. Both sharp crack and blunted notch crack models are analyzed over a wide range of applied remote tensile strains. Comparisons to linear elastic fracture mechanics are presented. A critical condition is proposed for crack growth in hydrogels based on the transient energy release rate. The applicability of this growth condition for simulating concomitant crack propagation and solvent diffusion in hydrogels is discussed. / text

The Stefan-Maxwell equations are a system of nonlinear partial differential equations that describe the diffusion of multiple chemical species in a container. These equations are of particular interest for their applications to biology and chemical engineering. The nonlinearity and coupled nature of the equations involving many variables make finding solutions difficult, so numerical methods are often used. In the engineering literature the system is inverted to write fluxes as functions of the species gradient before any numerical method is applied. In this thesis it is shown that employing a mixed finite element method makes the inversion unnecessary, allowing the numerical solution of Stefan-Maxwell equations in their primitive form. The plan of the thesis is as follows, first a mixed variational formulation will be derived for the Stefan-Maxwell equations. The nonlinearity will be dealt with through a linearization. Conditions for well-posedness of the linearized formulation are then determined. Next, the linearized variational formulation is approximated using mixed finite element methods. The finite element methods will then be shown to converge to an approximate solution. A priori error estimates are obtained between the solution to the approximate problem and the exact solution. The convergence order is then verified through an analytic test case and compared to standard methods. Finally, the solution is computed for another test case involving the diffusion of three species and compared to other methods.

The Stefan-Maxwell equations are a system of nonlinear partial differential equations that describe the diffusion of multiple chemical species in a container. These equations are of particular interest for their applications to biology and chemical engineering. The nonlinearity and coupled nature of the equations involving many variables make finding solutions difficult, so numerical methods are often used. In the engineering literature the system is inverted to write fluxes as functions of the species gradient before any numerical method is applied. In this thesis it is shown that employing a mixed finite element method makes the inversion unnecessary, allowing the numerical solution of Stefan-Maxwell equations in their primitive form. The plan of the thesis is as follows, first a mixed variational formulation will be derived for the Stefan-Maxwell equations. The nonlinearity will be dealt with through a linearization. Conditions for well-posedness of the linearized formulation are then determined. Next, the linearized variational formulation is approximated using mixed finite element methods. The finite element methods will then be shown to converge to an approximate solution. A priori error estimates are obtained between the solution to the approximate problem and the exact solution. The convergence order is then verified through an analytic test case and compared to standard methods. Finally, the solution is computed for another test case involving the diffusion of three species and compared to other methods.

In this dissertation, we study the mixed finite element method for the linear plane elasticity problem and iterative solvers for the resulting discrete system. We use the Arnold-Winther Element in the mixed finite element discretization. An overlapping Schwarz preconditioner and a multigrid preconditioner for the discrete system are developed and analyzed. We start by introducing the mixed formulation (stress-displacement formulation) for the linear plane elasticity problem and its discretization. A detailed analysis of the Arnold-Winther Element is given. The finite element discretization of the mixed formulation leads to a symmetric indefinite linear system. Next, we study efficient iterative solvers for the symmetric indefinite linear system which arises from the mixed finite element discretization of the linear plane elasticity problem. The preconditioned Minimum Residual Method is considered. It is shown that the problem of constructing a preconditioner for the indefinite linear system can be reduced to the problem of constructing a preconditioner for the H(div) problem in the Arnold-Winther finite element space. Our main work involves developing an overlapping Schwarz preconditioner and a multigrid preconditioner for the H(div) problem. We give condition number estimates for the preconditioned systems together with supporting numerical results.

The miscible displacement equations provide the mathematical model for simulating the displacement of a mixture of oil and miscible fluid in underground reservoirs during the Enhance Oil Recovery(EOR) process. In this thesis, I propose a stable numerical scheme combining a mixed finite element method and space-time discontinuous Galerkin method for solving miscible displacement equations under low regularity assumption. Convergence of the discrete solution is investigated using a compactness theorem for functions that are discontinuous in space and time. Numerical experiments illustrate that the rate of convergence is improved by using a high order time stepping method. For petroleum engineers, it is essential to compute finely detailed fluid profiles in order to design efficient recovery procedure thereby increase production in the EOR process. The method I propose takes advantage of both high order time approximation and discontinuous Galerkin method in space and is capable of providing accurate numerical solutions to assist in increasing the production rate of the miscible displacement oil recovery process.

In this dissertation, we study the mixed finite element method for the linear plane elasticity problem and iterative solvers for the resulting discrete system. We use the Arnold-Winther Element in the mixed finite element discretization. An overlapping Schwarz preconditioner and a multigrid preconditioner for the discrete system are developed and analyzed. We start by introducing the mixed formulation (stress-displacement formulation) for the linear plane elasticity problem and its discretization. A detailed analysis of the Arnold-Winther Element is given. The finite element discretization of the mixed formulation leads to a symmetric indefinite linear system. Next, we study efficient iterative solvers for the symmetric indefinite linear system which arises from the mixed finite element discretization of the linear plane elasticity problem. The preconditioned Minimum Residual Method is considered. It is shown that the problem of constructing a preconditioner for the indefinite linear system can be reduced to the problem of constructing a preconditioner for the H(div) problem in the Arnold-Winther finite element space. Our main work involves developing an overlapping Schwarz preconditioner and a multigrid preconditioner for the H(div) problem. We give condition number estimates for the preconditioned systems together with supporting numerical results.