Many engineering and scientific problems require the solution of partial differential equations in complex geometries. Often, these problems involve parametrized geometries, e.g. design optimization, or moving domains, e.g. fluid-structure interaction problems. For such cases, traditional methods based on body-fitted grids require time-consuming mesh generation or re-meshing techniques. Unfitted finite element methods, e.g. CutFEM of AgFEM, are appealing techniques that address these challenges. However, they require ad-hoc integration methods and stabilization techniques to prevent instabilities for small cut cells. Recently, the Shifted Boundary Method (SBM), was introduced to prevent integration over cut cells and small cut-cell instabilities. An extension of the SBM was recently introduced, the Weighted Shifted Boundary Method (WSBM), where the variational form is weighted by the elemental active volume fraction, improving discrete mass/momentum conservation properties in simulations with moving domains. In this work we introduce the Generalized Shifted Boundary Method (GSBM), a geometry-agnostic generalization of the SBM and WSBM formulations that avoids the need of redefinition of integration domains and finite element spaces. The GSBM enables a unified formulation for problems with evolving geometries, supports gradient-based optimization of problems with varying geometries including topological changes, and unifies SBM, WSBM, and optimal-surrogate variants within a single framework. In this work we describe the formulation, and corresponding tests, for three model problems, namely: the Poisson problem, linear elasticity and transient Stokes flow.

The Weighted Shifted Boundary Method (WSBM) was recently introduced as an enhanced Shifted Boundary Method (SBM) for the simulation of flows with moving boundaries. Earlier work of the authors on no-slip boundary conditions for the two-dimensional Stokes flow is extended here to the more challenging case of the three-dimensional incompressible Navier-Stokes equations at low and moderate Reynolds numbers. The SBM is an immersed finite element method that reformulates an infinite-dimensional boundary value problem over a surrogate (approximate) computational domain – to avoid integrating over cut cells – and modifies the original boundary conditions using Taylor expansions – to maintain accuracy. The WSBM weights the SBM’s variational form with the elemental volume fraction of active fluid, drastically reducing spurious pressure oscillations in time that occur when the total volume of active fluid changes abruptly over a time step. The WSBM induces small mass (i.e., volume) conservation errors, which converge quadratically in the case of piecewise-linear finite element interpolations, as the grid is refined. An extensive set of two- and three-dimensional tests demonstrates the robustness and accuracy of the proposed approach.

The Shifted Boundary Method (SBM) belongs to the class of unfitted (or immersed, or embedded) finite element methods, and relies on reformulating the original boundary value problem over a surrogate (approximate) computational domain. The surrogate domain is constructed so as to avoid cut cells and the associated problematic implementation and numerical integration issues. Accuracy is maintained by modifying the original boundary conditions using Taylor expansions: hence the name of the method, that shifts the location and values of the boundary conditions. In this article, we extend the SBM to the simulation of incompressible Stokes flow, by appropriately weighting its variational form with the elemental volume fraction of active fluid. This approach allows to drastically reduce spurious pressure oscillations in time, which are produced if the total volume of active fluid were to change abruptly over a time step. The proposed Weighted SBM (W-SBM) exactly preserves states of hydrostatic equilibrium, and induces small mass and momentum conservation errors, which converge as the grid is refined. This is in analogy to cutFEMs and related unfitted approaches, which rely on an affine representation of cut boundaries. We demonstrate the robustness and accuracy of the proposed method with an extensive suite of two-dimensional tests.

We propose a Lagrangian solid mechanics framework for the simulation of salt tectonics and other large-deformation geomechanics problems at the basin scale. Our approach relies on general elastic-viscoplastic constitutive models to characterize the deformation of geologic strata, in contrast with the majority of published works on the subject, which utilize nonlinear Stokes flow models. By means of multiscale asymptotics, we also show that the inertia term in the momentum balance equation can be safely neglected, if the goal is to track the Earth's crust deformation over long periods of time. Our time integration strategy is a blended transient/quasistatic approach, in that it consists of a constitutive stress update, subject to the constraint that the stresses must satisfy static equilibrium. In addition, we use stabilized finite element methods specifically built for triangular and tetrahedral grids, which can also perform well under incompressibility constraints. Our approach offers computational geologists the following advantages: (1) improved flexibility in the choice of subsurface constitutive models with respect to the nonlinear Stokes flow; (2) improved efficiency over transient dynamics algorithms used in this context in the past, which are forced to resolve seismic events over geologic time scales; and (3) improved robustness in large strain computations over quadrilateral/hexahedral finite elements. We demonstrate the performance of the proposed approach with simulations of passive diapirism.

The Shifted Boundary Method (SBM) belongs to the class of unfitted (or immersed, or embedded) finite element methods and was recently introduced for the Poisson, linear advection/diffusion, Stokes, Navier-Stokes, acoustics, and shallow-water equations. By reformulating the original boundary value problem over a surrogate (approximate) computational domain, the SBM avoids integration over cut cells and the associated problematic issues regarding numerical stability and matrix conditioning. Accuracy is maintained by modifying the original boundary conditions using Taylor expansions. Hence the name of the method, that shifts the location and values of the boundary conditions. In this article, we extend the SBM to the simulation of incompressible Navier-Stokes flows with moving free-surfaces, by appropriately weighting its variational form with the elemental volume fraction of active fluid. This approach prevents spurious pressure oscillations in time, which would otherwise be produced if the total active fluid volume were to change abruptly over a time step. In fact, the proposed weighted SBM method induces small mass (i.e., volume) conservation errors, which converge quadratically in the case of piecewise-linear finite element interpolations, as the grid is refined. We present an extensive set of two- and three-dimensional tests to demonstrate the robustness and accuracy of the method.

The discontinuous Galerkin (DG) method has found widespread application in elliptic problems with rough coefficients, of which the Darcy flow equations are a prototypical example. One of the long-standing issues of DG approximations is the overall computational cost, and many different strategies have been proposed, such as the variational multiscale DG method, the hybridizable DG method, the multiscale DG method, the embedded DG method, and the Enriched Galerkin method. In this work, we propose a mixed dual-scale Galerkin method, in which the degrees-of-freedom of a less computationally expensive coarse-scale approximation are linked to the degrees-of-freedom of a base DG approximation. We show that the proposed approach has always similar or improved accuracy with respect to the base DG method, with a considerable reduction in computational cost. For the specific definition of the coarse-scale space, we consider Raviart–Thomas finite elements for the mass flux and piecewise-linear continuous finite elements for the pressure. We provide a complete analysis of stability and convergence of the proposed method, in addition to a study on its conservation and consistency properties. We also present a battery of numerical tests to verify the results of the analysis, and evaluate a number of possible variations, such as using piecewise-linear continuous finite elements for the coarse-scale mass fluxes.

The numerical simulation of physical phenomena and engineering problems can be affected by numerical errors and various types of uncertainties. Characterizing the former in computational frameworks involving system parameter uncertainties becomes a key issue. In this work, we study the behavior of new variational multiscale (VMS) error estimators for the propagation of parametric uncertainties in a Convection–Diffusion–Reaction (CDR) problem. A sensitivity analysis is performed to assess the performance of the error estimator with respect to the mesh discretization and physical parameters (here, the viscosity value and advection velocity). Three different manufactured analytical solutions are considered as benchmarking tests. Next, the performance of the VMS error estimators is evaluated for the CDR problem with uncertain input parameters. For this purpose, two probabilistic models are constructed for the viscosity and advection direction, and the uncertainties are propagated using a polynomial chaos expansion approach. A convergence analysis is specifically carried out for different configurations, corresponding to regimes where the CDR operator is either smooth or non-smooth. An assessment of the proposed error estimator is finally conducted for the three tests, considering both the viscous- and convection-dominated regimes.

In this work, we define a family of explicit a posteriori error estimators for Finite Volume methods in computational fluid dynamics. The proposed error estimators are inspired by the Variational Multiscale method, originally defined in a Finite Element context. The proposed error estimators are tested in simulations of the incompressible Navier-Stokes equations, the thermally-coupled Navier-Stokes equations, and the fully-coupled compressible large eddy simulation of the HIFiRE Direct Connect Rig Scramjet combustor.

In this article, we develop a dynamic version of the variational multiscale (D-VMS) stabilization for nearly/fully incompressible solid dynamics simulations of viscoelastic materials. The constitutive models considered here are based on Prony series expansions, which are rather common in the practice of finite element simulations, especially in industrial/commercial applications. Our method is based on a mixed formulation, in which the momentum equation is complemented by a pressure equation in rate form. The unknown pressure, displacement, and velocity are approximated with piecewise linear, continuous finite element functions. To prevent spurious oscillations, the pressure equation is augmented with a stabilization operator specifically designed for viscoelastic problems, in that it depends on the viscoelastic dissipation. We demonstrate the robustness, stability, and accuracy properties of the proposed method with extensive numerical tests in the case of linear and finite deformations.