Streamline Upwind Petrov-Galerkin Research Articles

Stabilized isogeometric formulation of the multi-network poroelasticity and transport model (MPET[formula omitted]) for subcutaneous injection of monoclonal antibodies

The MPET2 model couples the multi-network poroelastic theory (MPET) with solute transport equations and provides predictions of the material deformation, fluid dynamics, and solute transport in different compartments of a deformable multiple-porosity medium. MPET2 offers a comprehensive framework for understanding complex porous media across multiple disciplines. Examples of its applications include studying rock formations, soil mechanics and subsurface reservoirs, investigating biological tissues, modeling groundwater flow and contaminant transport, and optimizing the design of porous materials. Despite the wide range of applications of the model, its numerical discretization has received little attention. Here we propose a stabilized formulation of the MPET2 model. To address the unique challenges posed by the discretization of the MPET2 model, we use multiple techniques including the Fluid Pressure Laplacian stabilization, Streamline Upwind Petrov–Galerkin stabilization, and discontinuity capturing. Our spatial discretization is based on Isogeometric Analysis with higher-order continuity basis functions. The fully discretized governing equations are solved simultaneously with a monolithic algorithm. We perform a convergence study of the proposed formulation. Then, we conduct a series of simulations of subcutaneous injection of monoclonal antibodies under different injection conditions. Our simulations show that the stabilized MPET2 formulation can provide oscillation-free solutions for tissue deformation, fluid flow in the interstitial tissue, blood vessels, and lymphatic vessels, drug absorption in blood vessels and lymphatic vessels, as well as drug transport in each compartment. We also study the effects of different injection conditions on drug absorption, showing the potential of the proposed model and algorithm in the future optimization of injection strategy.

Simulation of Drag Crisis in Flow Past a Circular Cylinder using 2D Computations

The unsteady incompressible 2D Navier-Stokes equations are solved using a finite element method with Streamline-Upwind/Petrov-Galerkin (SUPG) and Pressure-Stabilizing/Petrov- Galerkin (PSPG) stabilization terms for high Reynolds number (105 to I06) flow past a cylinder to stimulate drag-crisis. The mean drag coefficient for various Re are compared with experimental data. It is observed that the low at Re =106 separates much later, from the cylinder surface as compared to that at Re =105. The wake at Re =106 is narrower than that at Re = 105. The time-averaged flow field shows a standing bubble just downstream of the flow separation point for the Re = 105 flow. The bubble is not observed for the flow at Re = 106. The energy spectrum of the flow exhibits the characteristics of 2D isotropic turbulence.

Comparing different stabilization strategies for reduced order modeling of viscoelastic fluid flow problems

This paper presents a computational framework for model order reduction of viscoelastic fluid flows, with a particular focus on the stabilization of the reduced order model (ROM) and the efficient approximation of non-linear terms appearing in the governing equations. We compare three different stabilization approaches, two of which are based on offline stabilization only and one that uses both offline and online stabilization. Two main techniques for the treatment of non-linear terms are examined and compared: the discrete empirical interpolation method, and a reduced quadrature method based on a sparse recovery technique. Numerical experiments are conducted on two benchmark flows such as the flow past a sphere and the 4:1 contraction flow problem. We show that significant speedups can be achieved using offline streamline upwind Petrov–Galerkin stabilization together with hyper-reduction. We also discuss the potential accuracy losses of the ROM that could result from mesh-related issues of the full-order model (FOM) and propose possible remedies.

QPDE: Quantum Neural Network Based Stabilization Parameter Prediction for Numerical Solvers for Partial Differential Equations

We propose a Quantum Neural Network (QNN) for predicting stabilization parameter for solving Singularly Perturbed Partial Differential Equations (SPDE) using the Streamline Upwind Petrov Galerkin (SUPG) stabilization technique. SPDE-Q-Net, a QNN, is proposed for approximating an optimal value of the stabilization parameter for SUPG for 2-dimensional convection-diffusion problems. Our motivation for this work stems from the recent progress made in quantum computing and the striking similarities observed between neural networks and quantum circuits. Just like how weight parameters are adjusted in traditional neural networks, the parameters of the quantum circuit, specifically the qubits’ degrees of freedom, can be fine-tuned to learn a nonlinear function. The performance of SPDE-Q-Net is found to be at par with SPDE-Net, a traditional neural network-based technique for stabilization parameter prediction in terms of the numerical error in the solution. Also, SPDE-Q-Net is found to be faster than SPDE-Net, which projects the future benefits which can be earned from the speed-up capabilities of quantum computing.

Space–time computational flow analysis: Unconventional methods and first-ever solutions

The Deforming-Spatial-Domain/Stabilized Space–Time (DSD/SST) method was introduced in 1990 as a moving-mesh method for computational analysis of flows with moving boundaries and interfaces (MBI), which is a wide class of problems that includes fluid–particle and fluid–structure interactions and free-surface and multi-fluid flows. The method was inspired by Thomas Hughes’s 1987 work on space–time finite element methods for elastodynamics. The original DSD/SST method is now called “ST-SUPS”, reflecting its stabilization components, which are the Streamline-Upwind/Petrov–Galerkin (SUPG) method, pioneered by Hughes, and the Pressure-Stabilizing/Petrov–Galerkin (PSPG) method, inspired by Hughes’s work on a Stokes-flow Petrov–Galerkin formulation allowing equal-order interpolations for velocity and pressure. Hughes’s work on the residual-based variational multiscale (RBVMS) method inspired the ST-VMS method, which is the VMS version of the DSD/SST. A number of special methods were introduced in connection with the core methods ST-SUPS and ST-VMS. Hughes’s work on isogeometric analysis (IGA) inspired one of those special methods, the “ST-IGA”, with IGA basis functions not only in space but also in time. The core and special ST methods enabled first-ever solutions in some of the most challenging classes of MBI problems, including particle-laden flows, spacecraft parachute fluid–structure interactions, and car and tire aerodynamics. We provide an overview of the ST methods inspired by Hughes’s work and highlight some of the first-ever solutions in these three classes of MBI problems.

Stabilization in 3‐D FEM and solution of the MHD equations

In this study, we consider the numerical solution of convection‐diffusion typed equations defined in 3‐D domain using the finite element method (FEM) with the stabilized version in order to fill the gap in literature on this area in the sense of different applications. The formulation processed is performed step by step. For this purpose, we started with FEM formulation of the 3‐D Laplace equation. Then, the obtained formulation is extended to the convection‐diffusion equation with standard Galerkin FEM. In order to solve the stability problems for the convection dominated case, the most popular stabilized formulation known as the streamline upwind Petrov‐Galerkin (SUPG) type stabilization method is considered. As a further step, obtained stabilized formulation is used in the solutions of the MHD equations by transforming the coupled system of differential equations to the decoupled convection‐diffusion type equations defined on different 3‐D domains. For the comparison purpose, MHD equations are also formulated with finite volume method (FVM). Obtained numerical results are displayed in terms of a table and figures as the 2‐D slices of the 3‐D plots in order to visualize the effect and accuracy of the proposed numerical scheme over the different test problems.

Adaptive Finite Element Simulation of Double-Diffusive Convection

Double-diffusive convection plays an important role in many physical phenomena of practical importance. However, the numerical simulation of these phenomena is challenging since fine meshes are often required to capture the flow physics. Hence, several different numerical methods have been employed in the past. This work reports the development and application of an adaptive finite element method for the simulation of these phenomena, thereby avoiding the need for the use of very fine meshes over the whole domain. The weak formulation of the conservation equations for mass, momentum, energy and species concentration is used. The Boussinesq approximation relates the density of the fluid to the temperature and/or the species concentration. A second-order backward difference method is used for time discretization and the Galerkin method is employed for spatial discretization. Both adaptive time step and grid refinement techniques are employed, and the code is parallelized using MPI. Three different stabilization methods of the convective-diffusion equations are compared; namely, the streamline upwind Petrov–Galerkin (SUPG) method, and two modified methods aimed at diminishing spurious oscillations that include an artificial diffusion term. This diffusion term may be either isotropic or orthogonal to the streamlines. The addition of artificial isotropic diffusion to the SUPG method provides enhanced stability. The method is applied to double-diffusive finger convection in a sucrose-salt aqueous mixture and a stratified salt solution heated from below. The method accurately reproduces the experimentally observed temporal evolution of the salt fingers in the former case and the location of the interfaces between convective and non-convective zones in the latter.

SUPG formulation augmented with YZβ shock‐capturing for computing shallow‐water equations

AbstractWe demonstrate that the streamline‐upwind/Petrov–Galerkin (SUPG) formulation enhanced with YZβ discontinuity‐capturing, that is, the SUPG‐YZβ formulation, is an efficient and robust method for computing 2D shallow‐water equations (SWEs). The SUPG‐stabilized semi‐discrete formulation is discretized in time by employing the backward Euler time‐integration scheme. The nonlinear equation systems arising from the space and time discretizations are handled using the Newton–Raphson (N–R) method at each time step. The resulting linear equation systems are solved directly at each nonlinear iteration. Two challenging test problems are provided to examine the performance of the proposed formulation and techniques. To that end, we consider a full dam‐break and a partial dam‐break problem. We develop the solvers in the FEniCS environment. Test computations reveal that the SUPG‐YZβ formulation successfully eliminates spurious oscillations that cannot be captured with the SUPG‐stabilized formulation alone in narrow regions where steep gradients occur.

Spectral Analysis of High Order Continuous FEM for Hyperbolic PDEs on Triangular Meshes: Influence of Approximation, Stabilization, and Time-Stepping

In this work we study various continuous finite element discretization for two dimensional hyperbolic partial differential equations, varying the polynomial space (Lagrangian on equispaced, Lagrangian on quadrature points (Cubature) and Bernstein), the stabilization techniques (streamline-upwind Petrov–Galerkin, continuous interior penalty, orthogonal subscale stabilization) and the time discretization (Runge–Kutta (RK), strong stability preserving RK and deferred correction). This is an extension of the one dimensional study by Michel et al. (J Sci Comput 89(2):31, 2021. https://doi.org/10.1007/s10915-021-01632-7), whose results do not hold in multi-dimensional frameworks. The study ranks these schemes based on efficiency (most of them are mass-matrix free), stability and dispersion error, providing the best CFL and stabilization coefficients. The challenges in two-dimensions are related to the Fourier analysis. Here, we perform it on two types of periodic triangular meshes varying the angle of the advection, and we combine all the results for a general stability analysis. Furthermore, we introduce additional high order viscosity to stabilize the discontinuities, in order to show how to use these methods for tests of practical interest. All the theoretical results are thoroughly validated numerically both on linear and non-linear problems, and error-CPU time curves are provided. Our final conclusions suggest that Cubature elements combined with SSPRK and OSS stabilization is the most promising combination.

A Variational Multiscale Method for Particle Dispersion Modeling in the Atmosphere

A LES model is proposed to predict the dispersion of particles in the atmosphere in the context of Chemical, Biological, Radiological and Nuclear (CBRN) applications. The code relies on the Finite Element Method (FEM) for both the fluid and the dispersed solid phases. Starting from the Navier-Stokes equations and a general description of the FEM strategy, the Streamline Upwind Petrov-Galerkin (SUPG) method is formulated putting some emphasis on the related assembly matrix and stabilization coefficients. Then, the Variational Multiscale Method (VMS) is presented together with a detailed illustration of its algorithm and hierarchy of computational steps. It is demonstrated that the VMS can be considered as a more general version of the SUPG method. The final part of the work is used to assess the reliability of the implemented predictor/multicorrector solution strategy.

Thermal and Surface Roughness Effects on Performance Analysis of Slider Bearing via Petrov–Galerkin Finite Element Method

In this study, the performance of slider bearing with impacts of temperature and roughness of surface on 1D longitudinal and transverse roughness type is investigated using the Streamline Upwind Petrov–Galerkin (SUPG)-finite element method (FEM). It is considered that the roughness is stochastic and Gaussian in random distribution. It is also considered that viscosity and density depend upon temperature. The domain’s irregularity caused by surface roughness is changed into a regular domain for numerical computation purposes. To determine the capacity of load-carrying and pressure distribution, the continuity, momentum, modified Reynolds, and energy equations are decoupled and solved by the SUPG-FEM. It can be demonstrated that in the case of the longitudinal model of nonparallel slider bearing, the load-bearing bearing capacity and drag force of friction due to the case i of the combined effects are lower than those attributable to the influence of surface roughness. Nevertheless, the thermal case’s impact on the 1D transverse type capacity of load-carrying capability and drag friction force is below that of the combined and surface roughness impact instances. The reverse is true for parallel slider bearings, however, for a 1D transverse type, the load-carrying ability is not particularly impressive.

Improving the accuracy of discretisations of the vector transport equation on the lowest-order quadrilateral Raviart-Thomas finite elements

Within finite element models of fluids, vector-valued fields such as velocity or momentum variables are commonly discretised using the Raviart-Thomas elements. However, when using the lowest-order quadrilateral Raviart-Thomas elements, standard finite element discretisations of the vector transport equation typically have a low order of spatial accuracy. This paper describes two schemes that improve the accuracy of transporting such vector-valued fields on two-dimensional curved manifolds.The first scheme that is presented reconstructs the transported field in a higher-order function space, where the transport equation is then solved. The second scheme applies a mixed finite element formulation to the vector transport equation, simultaneously solving for the transported field and its vorticity. An approach to stabilising this mixed vector-vorticity formulation is presented that uses a Streamline Upwind Petrov-Galerkin (SUPG) method. These schemes are then demonstrated, along with their accuracy properties, through some numerical tests. Two new test cases are used to assess the transport of vector-valued fields on curved manifolds, solving the vector transport equation in isolation. The improvement of the schemes is also shown through two standard test cases for rotating shallow-water models.

Local maximum-entropy approximation based stabilization methods for the convection diffusion problems

Local Maximum Entropy (LME) approximation in the field of elliptic partial differential equations is more advanced compared to convection dominated flow problems. This paper presents a comparative study on the LME based meshfree stabilization methods for convection dominated convection-diffusion problems. LME approximants account for the disposal of the elements and direct imposition of boundary conditions, whereas the stabilization methods are meant to deal with the issues associated with the non-self-adjoint convective term. A new meshfree stabilization scheme, residual free bubbles based multiscale element free Galerkin method is also presented and compared with streamline upwind Petrov Galerkin meshfree method and variational multiscale element free Galerkin Method. The effect of different priors and radius of support which define the LME basis are studied to test their importance in the three methods with respect to standard benchmark problems. Cubic spline prior based LME is found to be enough to produce better stability and accuracy for all the stabilization methods when compared to quartic spline and Gaussian priors. LME based variational multiscale element free Galerkin method is observed to yield better stability and accuracy for all the priors. The convergence rate for the three stabilization methods is observed to be in good agreement with the classical FEM-SUPG.

Uniform stability of the SUPG method for the evolutionary convection–diffusion–reaction equation

To solve the evolutionary convection-diffusion-reaction equation, the implicit Euler method and the streamline upwind Petrov Galerkin (SUPG) method are used for time discretization and space discretization respectively in this paper. Uniform stability is obtained with respect to the temporal and spatial discretization parameters. A new projection is introduced for the stability analysis instead of the common idea of material derivative. The analysis does not require a strict relation between time step and space step, which is usually needed by derivations based on the technique of material derivative. Furthermore, it can handle the case of time-dependent convection and reaction coefficients.

On the performance of a Chimera-FEM implementation to treat moving heat sources and moving boundaries in time-dependent problems

Problems with moving sources and moving inner boundaries in transient regime are of high interest in many research fields and engineering applications. One approach to properly tackle such problems is based on the Chimera method for non-matching grids, where each moving object is defined on a fine mesh that moves across the fixed coarse background. In this way, the high gradients around moving sources or boundaries are captured by the fine mesh without need of globally fine fixed meshes or adaptive refinement. In this work, Chimera is implemented in the framework of the finite element method, following an arbitrary Lagrangian–Eulerian formulation for the moving meshes and a standard Eulerian formulation for the fixed background mesh. Further, in case of convection-dominated problems, the scheme is stabilized using the streamline upwind Petrov–Galerkin method. The coupling between the moving fine meshes and the coarse fixed one is achieved via Dirichlet boundary conditions and a high-order interpolation algorithm. The performance of the proposed methodology in terms of accuracy and stability is assessed by means of numerical tests to be compared with equivalent problems solved using fixed meshes. Further tests serve to highlight the good performance of the proposed Chimera-based finite element method to address both convection-dominated problems with multiple moving boundaries and sources, and three-dimensional arc welding processes. • An algebraic-based Chimera-FEM scheme for unsteady problems is proposed. • Remeshing is avoided for large movements of heat sources and boundaries. • Advective-dominant problems can be addressed with Chimera-FEM and SUPG. • Good accuracy is obtained compared to conventional one-domain approaches. • The scheme is suitable to be used on 3D welding process simulations.

Totally volume integral of fluxes for discontinuous Galerkin method-two dimensional Euler equations

In this paper, the scheme of constructing high-order accurate totally volume discontinuous finite element method for the numerical solution of the 1D Euler equations is extended to 2D Euler equations on Cartesian meshes. In the present work, the boundary integral fluxes are transformed into volume integral by applying divergence theorem to the boundary integral of the Riemann fluxes. Therefore, the totally volume discontinuous finite element is independent on the boundary integral fluxes at the element boundaries as opposed to the classical discontinuous Galerkin method. The accuracy is obtained by applying high-order polynomial approximations within elements using the tensor product of Lagrange polynomial. For temporal integration, strong stability preserving Runge-Kutta method SSPRK (3, 3) is applied. The scheme is stabilized by using Streamline Upwind Petrov Galerkin (SUPG) stabilization technique. For the spatial discretization, the polynomial of order 1and 2 are used, the shape function is constructed for the master (computational) element after applying the coordinate transformation for the physical domain, the transformation for the governing equations is performed to get it in the function of computational Cartesian, then the governing equations are put in conservative form. The numerical results of applying totally volume integral discontinuous Galerkin method for two-dimensional Euler equations presented in this paper show that the scheme is very accurate, fast, and effective even with shock appearance.

Error analysis of a SUPG-stabilized POD-ROM method for convection-diffusion-reaction equations

A reduced order model (ROM) method based on proper orthogonal decomposition (POD) is analyzed for convection-diffusion-reaction equations. The streamline-upwind Petrov–Galerkin (SUPG) stabilization is used in the practically interesting case of dominant convection, both for the full order method (FOM) and the ROM simulations. The asymptotic choice of the stabilization parameter for the SUPG-ROM is done as proposed in the literature. This paper presents a finite element convergence analysis of the SUPG-ROM method for errors in different norms. The constants in the error bounds are uniform with respect to small diffusion coefficients. An approach for a posteriori error estimation is discussed. Numerical studies illustrate the performance of the SUPG-ROM method.

The Finite Element Numerical Study of Transonic Flow of Moist Air with Nonequilibrium Condensation

The finite element numerical simulation of the transonic flow with condensation is studied in this paper. The transformed weak form is numerically discretized by using the classical Galerkin finite element spatial discretization scheme and the backward Euler difference time discretization scheme. Because the coupled Euler equations are convection equations, the results obtained by the classical Galerkin method often have nonphysical numerical oscillations, so the stabilization method is introduced. Another important feature of compressible Euler equation is the existence of shock wave. Because the thickness of shock wave is small and difficult to capture, the shock wave can be better captured by adding artificial viscosity term. Here, the streamline upwind/Petrov–Galerkin (SUPG) finite element method with good stability and the isotropic diffusion method with shock capture ability are used to solve the problem of transonic flow with condensation. In the simulation of steady problems, in order to improve the calculation accuracy of the existing area of the shock wave, the adaptive mesh refinement technology is added, which refines the mesh in the region of the shock wave to improve the resolution of the shock wave. Numerical experiments verify the feasibility and stabilization of the numerical method. Finally, the simulation of transonic flow around NACA-0012 and parallel-jet nozzle A1 has obtained the better numerical results.

Extending FEniCS to work in higher dimensions using tensor product finite elements

We present a method to extend the finite element library FEniCS to solve problems with domains in dimensions above three by constructing tensor product finite elements. This methodology only requires that the high dimensional domain is structured as a Cartesian product of two lower dimensional subdomains. In this study we consider Dirichlet problems for scalar linear partial differential equations, though the methodology can be extended to non-linear problems. The utilization of tensor product finite elements allows us to construct a global system of linear algebraic equations that only relies on the finite element infrastructure of the lower dimensional subdomains contained in FEniCS. We demonstrate the effectiveness of our methodology in four distinctive test cases. The first test case is a Poisson equation posed in a four dimensional domain which is a Cartesian product of two unit squares solved using the classical Galerkin finite element method. The second test case is the wave equation in space–time, where the computational domain is a Cartesian product of a two dimensional space grid and a one dimensional time interval. In this second case we also employ the Galerkin method. The third test case is an advection dominated advection–diffusion equation where the global domain is a Cartesian product of two one dimensional intervals in which the streamline upwind Petrov–Galerkin method is applied to ensure discrete stability. The final test case uses the Galerkin approach to solve a Poisson problem on a Cartesian product of two intervals with a spatially varying, non-separable diffusivity term In all cases, a p=1 basis is used and optimal L2 convergence rates of order hp+1 of the errors are achieved with respect to h refinement.

A cell-based smoothed finite element model for non-Newtonian blood flow

Smoothed Finite Element Method (S-FEM) has drawn increasing attention in the field of computational fluid dynamics (CFD) and the present work seeks to make further contribution to this growing field of S-FEM by simulating for the non-Newtonian blood flow. This investigation took the form of Streamline Upwind Petrov-Galerkin in conjunction with Stabilized Pressure Gradient Projection (SUPG/SPGP) to alleviate the spatial oscillation and instability problems. The validation of the cell-based S-FEM (CS-FEM) combined with SUPG/SPGP was carried out by the blood flow over a backward-facing step. The performances of the presented method were explored by the blood flow in the carotid bifurcation and blood flow in the intracranial segment of internal carotid artery. Impressively, the results exhibit good features of the CS-FEM on solving severely distorted mesh vis-a-vis the standard finite element method (FEM). The presented method could realize the accurate prediction in the mentioned complex blood flows as the same as Finite Volume Method software STAR-CCM+.