Search for: All records

Abstract The numerical analysis of stochastic parabolic partial differential equations of the form$$\begin{aligned} du + A(u)\, dt = f \,dt + g \, dW, \end{aligned}$$ $\begin{array}{c}du+A\left(u\right)dt=fdt+gdW,\end{array}$is surveyed, whereAis a nonlinear partial operator andWa Brownian motion. This manuscript unifies much of the theory developed over the last decade into a cohesive framework which integrates techniques for the approximation of deterministic partial differential equations with methods for the approximation of stochastic ordinary differential equations. The manuscript is intended to be accessible to audiences versed in either of these disciplines, and examples are presented to illustrate the applicability of the theory. 

Abstract The modeling of coupled fluid transport and deformation in a porous medium is essential to predict the various geomechanical process such as CO2 sequestration, hydraulic fracturing, and so on. Current applications of interest, for instance, that include fracturing or damage of the solid phase, require a nonlinear description of the large deformations that can occur. This paper presents a variational energy‐based continuum mechanics framework to model large‐deformation poroelasticity. The approach begins from the total free energy density that is additively composed of the free energy of the components. A variational procedure then provides the balance of momentum, fluid transport balance, and pressure relations. A numerical approach based on finite elements is applied to analyze the behavior of saturated and unsaturated porous media using a nonlinear constitutive model for the solid skeleton. Examples studied include the Terzaghi and Mandel problems; a gas–liquid phase‐changing fluid; multiple immiscible gases; and unsaturated systems where we model injection of fluid into soil. The proposed variational approach can potentially have advantages for numerical methods as well as for combining with data‐driven models in a Bayesian framework. 

In this work, we apply a phase-field modeling framework to elucidate the interplay between nucleation and kinetics in the dynamic evolution of twinning interfaces. The key feature of this phase-field approach is the ability to transparently and explicitly specify nucleation and kinetic behavior in the model, in contrast to other regularized interface models. We use this to study 2 distinct problems where it is essential to explicitly specify the kinetic and nucleation behavior governing twin evolution. First, we study twinning interfaces in 2-d. When these interfaces are driven to move, we find that significant levels of twin nucleation occur ahead of the moving interface. Essentially, the finite interface velocity and the relaxation time of the stresses ahead of the interface allows for nucleation to occur before the interface is able to propagate to that point. Second, we study the growth of needle twins in antiplane elasticity. We show that both nucleation and anisotropic kinetics are essential to obtain predictions of needle twins. While standard regularized interface approaches do not permit the transparent specification of anisotropic kinetics, this is readily possible with the phase-field approach that we have used here. 

Accretion and ablation, i.e., the addition and removal of mass at the surface, are important in a wide range of physical processes, including solidification, growth of biological tissues, environmental processes, and additive manufacturing. The description of accretion requires the addition of new continuum particles to the body, and is therefore challenging for standard continuum formulations for solids that require a reference configuration. Recent work has proposed an Eulerian approach to this problem, enabling side-stepping of the issue of constructing the reference configuration. However, this raises the complementary challenge of determining the stress response of the solid, which typically requires the deformation gradient that is not immediately available in the Eulerian formulation. To resolve this, the approach introduced the elastic deformation as an additional kinematic descriptor of the added material, and its evolution has been shown to be governed by a transport equation. In this work, the method of characteristics is applied to solve concrete simplified problems motivated by biomechanics and manufacturing. Specifically, (1) for a problem with both ablation and accretion in a fixed domain and (2) for a problem with a time-varying domain, the closed-form solution is obtained in the Eulerian framework using the method of characteristics without explicit construction of the reference configuration.