An official website of the United States government
In this paper, we propose and study first- and second-order (in time) stabilized linear finite element schemes for the incompressible Navier-Stokes (NS) equations. The energy, momentum, and angular momentum conserving (EMAC) formulation has emerged as a promising approach for conserving energy, momentum, and angular momentum of the NS equations, while the exponential scalar auxiliary variable (ESAV) has become a popular technique for designing linear energy-stable numerical schemes. Our method leverages the EMAC formulation and the Taylor-Hood element with grad-div stabilization for spatial discretization. We adopt the implicit-explicit backward differential formulas (BDFs) coupled with a novel stabilized ESAV approach for time stepping. For the solution process, we develop an efficient decoupling technique for the resulting fully-discrete systems so that only one linear Stokes solve is needed at each time step, which is similar to the cost of classic implicit-explicit BDF schemes for the NS equations. Robust optimal error estimates are successfully derived for both velocity and pressure for the two proposed schemes, with Gronwall constants that are particularly independent of the viscosity. Furthermore, it is rigorously shown that the grad-div stabilization term can greatly alleviate the viscosity-dependence of the mesh size constraint, which is required for error estimation when such a term is not present in the schemes. Various numerical experiments are conducted to verify the theoretical results and demonstrate the effectiveness and efficiency of the grad-div and ESAV stabilization strategies and their combination in the proposed numerical schemes, especially for problems with high Reynolds numbers.
more »
« less
Implementation and verification of a user‐defined element (UEL) for coupled thermal‐hydraulic‐mechanical‐chemical (THMC) processes in saturated geological media
Abstract Efficient and accurate modeling of the coupled thermal‐hydraulic‐mechanical‐chemical (THMC) processes in various rock formations is indispensable for designing energy geo‐structures such as underground repositories for high‐level nuclear wastes. This work focuses on developing and verifying an implicit finite element solver for generic coupled THMC problems in geological settings. Starting from the mass, momentum, and energy balance laws, a specialized set of governing equations and a thermoporoelastic constitutive model is derived. This system is then solved by an implicit finite element (FE) scheme. Specifically, the residuals and the Jacobians are scripted in a user‐defined element (UEL) subroutine which is then combined with the general‐purpose FE software Abaqus Standard to solve initial‐boundary value problems. Considering the complexity of the system, the UEL development follows a stepwise manner by first solving the coupled hydraulic‐mechanical (HM) and thermal‐hydraulic‐mechanical (THM) equations before moving on to the full THMC problem. Each implementation step consists of at least one verification test by comparing computed results with closed‐form analytical solutions to ensure that the various coupling effects are correctly realized. To demonstrate the robustness of the algorithm and to validate the UEL, a three‐dimensional case study is performed with reference to the in‐situ heating test of ATLAS at Belgium in 1980s. A hypothetical radionuclide leakage event is then simulated by activating the chemical‐concentration degree of freedom and prescribing a constant high concentration at the heater's surface. The model predicts a limited contaminated regime after six years considering both diffusion and advection effects on species transport.
more »
« less
- Award ID(s):
- 2113474
- PAR ID:
- 10419296
- Publisher / Repository:
- Wiley Blackwell (John Wiley & Sons)
- Date Published:
- Journal Name:
- International Journal for Numerical and Analytical Methods in Geomechanics
- Volume:
- 47
- Issue:
- 11
- ISSN:
- 0363-9061
- Page Range / eLocation ID:
- p. 2153-2190
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
Directed energy deposition (DED) has been widely used for component repair. In the repair process, the surface defects are machined to a groove or slot and then refilled. The sidewall inclination angle of the groove geometry has been recognized to have a considerable impact on the mechanical properties of repaired parts. The objective of this work was to investigate the feasibility of repairing various V-shaped defects with both experiments and modeling. At first, the repair volume was defined by scanning the defective zone. Then, the repair volume was sliced to generate the repair toolpath. After that, the DED process was used to deposit Ti6Al4V powder on the damaged plates with two different slot geometries. Mechanical properties of the repaired parts were evaluated by microstructure analysis and tensile test. Testing of the repaired parts showed excellent bonding between the deposits and base materials with the triangular slot repair. 3D finite element analysis (FEA) models based on sequentially coupled thermo-mechanical field analysis were developed to simulate the corresponding repair process. Thermal histories of the substrate on the repair sample were measured to calibrate the 3D coupled thermo-mechanical model. The temperature measurements showed very good verification with the predicted temperature results. After that, the validated model was used to predict the residual stresses and distortions in the parts. Predicted deformation and stress results can guide the evaluation of the repair quality.more » « less
-
We present a fully-coupled, implicit-in-time framework for solving a thermodynamically-consistent Cahn-Hilliard Navier-Stokes system that models two-phase flows. In this work, we extend the block iterative method presented in Khanwale et al. [Simulating two-phase flows with thermodynamically consistent energy stable Cahn-Hilliard Navier-Stokes equations on parallel adaptive octree based meshes, J. Comput. Phys. (2020)], to a fully-coupled, provably second-order accurate scheme in time, while maintaining energy-stability. The new method requires fewer matrix assemblies in each Newton iteration resulting in faster solution time. The method is based on a fully-implicit Crank-Nicolson scheme in time and a pressure stabilization for an equal order Galerkin formulation. That is, we use a conforming continuous Galerkin (cG) finite element method in space equipped with a residual-based variational multiscale (RBVMS) procedure to stabilize the pressure. We deploy this approach on a massively parallel numerical implementation using parallel octree-based adaptive meshes. We present comprehensive numerical experiments showing detailed comparisons with results from the literature for canonical cases, including the single bubble rise, Rayleigh-Taylor instability, and lid-driven cavity flow problems. We analyze in detail the scaling of our numerical implementation.more » « less
-
Trabecular bone, a solid that has a heterogeneous porous structure, demonstrates nonlinear stress–strain relationship, even within the small strain region, when subject to stresses. It also exhibits different responses when subject to tension and compression. This study presents the development of an implicit constitutive relation between the stress and the linearized strain specifically tailored for trabecular bone-like materials. The structure of the constitutive relation requires the solution of the balance of linear momentum and the constitutive relations simultaneously, and in view of this, a two-field mixed finite element model capable of solving general boundary value problems governed by a system of coupled equations is proposed. We investigate the effects of nonlinearity and heterogeneity in a dogbone-shaped sample. Our study is able to capture the significant nonlinear characteristics of the response of the trabecular bone undergoing small strains in experiments, in both tension and compression, very well.more » « less
-
Abstract In this article we formulate a stable computational nonlocal poromechanics model for dynamic analysis of saturated porous media. As a novelty, the stabilization formulation eliminates zero‐energy modes associated with the original multiphase correspondence constitutive models in the coupled nonlocal poromechanics model. The two‐phase stabilization scheme is formulated based on an energy method that incorporates inhomogeneous solid deformation and fluid flow. In this method, the nonlocal formulations of skeleton strain energy and fluid flow dissipation energy equate to their local formulations. The stable coupled nonlocal poromechanics model is solved for dynamic analysis by an implicit time integration scheme. As a new contribution, we validate the coupled stabilization formulation by comparing numerical results with analytical and finite element solutions for one‐dimensional and two‐dimensional dynamic problems in saturated porous media. Numerical examples of dynamic strain localization in saturated porous media are presented to demonstrate the efficacy of the stable coupled poromechanics framework for localized failure under dynamic loads.more » « less
