Thesis (Ph.D)--Computing, Georgia Institute of Technology, 2009. / Committee Chair: Vuduc, Richard; Committee Member: Biros, George; Committee Member: Davatzikos, Christos; Committee Member: Tannenbaum, Allen; Committee Member: Zhou, Hao Min. Part of the SMARTech Electronic Thesis and Dissertation Collection.

published_or_final_version / Civil Engineering / Master / Master of Philosophy

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

Not available / text

Not available / text

Not available / text

Not available / text

Not available / text

Not available / text

We present an algorithm for fully automatic hp-adaptivity for finite element approximations of elliptic and Maxwell boundary value problems in three dimensions. The algorithm automatically generates a sequence of coarse grids, and a corresponding sequence of fine grids, such that the energy norm of error decreases exponentially with respect to the number of degrees of freedom in either sequence. At each step, we employ a discrete optimization algorithm to determine the refinements for the current coarse grid such that the projection-based interpolation error for the current fine grid solution decreases with an optimal rate with respect to the number of degrees of freedom added by the refinement. The refinements are restricted only by the requirement that the resulting mesh is at most 1-irregular, but they may be anisotropic in both element size h and order of approximation p. While we cannot prove that our method converges at all, we present numerical evidence of exponential convergence for a diverse suite of model problems from acoustic and electromagnetic scattering. In particular we show that our method is well suited to the automatic resolution of exterior problems truncated by the introduction of a perfectly matched layer. To enable and accelerate the solution of these problems on commodity hardware, we include a detailed account of three critical aspects of our implementations, namely an efficient implementations of sum factorization, several interfaces to the direct multi-frontal solver MUMPS, and some fast direct solvers for the computation of a sequence of nested projections. / text