A pressure‐based method with AUSM‐type fluxes for MHD flows at arbitrary Mach numbers

In this paper, we propose an extension of a PISO method, previously developed to solve the Euler equations, and which is here extended to the ideal magnetohydrodynamic (MHD) equations. By following a pressure‐based approach, we make use of the flexibility given by pressure equation for calculating flows at arbitrary Mach numbers. To handle MHD discontinuities, we have adapted the MHD‐Advection Upstream Splitting Method for our pressure‐based formulation. With the purpose of validation, four sets of test cases are presented and discussed. We start with the circularly polarized Alfvén waves that serves as a smooth flow validation. The second case is the 1‐D Riemann problem that is calculated using both 1‐D and 2‐D formulation of the MHD equations. The third and fourth problems are the Orszag–Tang vortex and the supersonic low‐ β cylinder allowing validation of the method in complex 2‐D MHD shock interaction. Copyright © 2013 John Wiley & Sons, Ltd.     

The DRBEM solution of incompressible MHD flow equations

This paper presents a dual reciprocity boundary element method (DRBEM) formulation coupled with an implicit backward difference time integration scheme for the solution of the incompressible magnetohydrodynamic (MHD) flow equations. The governing equations are the coupled system of Navier‐Stokes equations and Maxwell's equations of electromagnetics through Ohm's law. We are concerned with a stream function‐vorticity‐magnetic induction‐current density formulation of the full MHD equations in 2D. The stream function and magnetic induction equations which are poisson‐type, are solved by using DRBEM with the fundamental solution of Laplace equation. In the DRBEM solution of the time‐dependent vorticity and current density equations all the terms apart from the Laplace term are treated as nonhomogeneities. The time derivatives are approximated by an implicit backward difference whereas the convective terms are approximated by radial basis functions. The applications are given for the MHD flow, in a square cavity and in a backward‐facing step. The numerical results for the square cavity problem in the presence of a magnetic field are visualized for several values of Reynolds, Hartmann and magnetic Reynolds numbers. The effect of each parameter is analyzed with the graphs presented in terms of stream function, vorticity, current density and magnetic induction contours. Then, we provide the solution of the step flow problem in terms of velocity field, vorticity, current density and magnetic field for increasing values of Hartmann number. Copyright © 2010 John Wiley & Sons, Ltd.     

An extension of the SIMPLE based discontinuous Galerkin solver to unsteady incompressible flows

In this paper, we present a SIMPLE based algorithm in the context of the discontinuous Galerkin method for unsteady incompressible flows. Time discretization is done fully implicit using backward differentiation formulae (BDF) of varying order from 1 to 4. We show that the original equation for the pressure correction can be modified by using an equivalent operator stemming from the symmetric interior penalty (SIP) method leading to a reduced stencil size. To assess the accuracy as well as the stability and the performance of the scheme, three different test cases are carried out: the Taylor vortex flow, the Orr‐Sommerfeld stability problem for plane Poiseuille flow and the flow past a square cylinder. (1) Simulating the Taylor vortex flow, we verify the temporal accuracy for the different BDF schemes. Using the mixed‐order formulation, a spatial convergence study yields convergence rates of k + 1 and k in the L²‐norm for velocity and pressure, respectively. For the equal‐order formulation, we obtain approximately the same convergence rates, while the absolute error is smaller. (2) The stability of our method is examined by simulating the Orr–Sommerfeld stability problem. Using the mixed‐order formulation and adjusting the penalty parameter of the symmetric interior penalty method for the discretization of the viscous part, we can demonstrate the long‐term stability of the algorithm. Using pressure stabilization the equal‐order formulation is stable without changing the penalty parameter. (3) Finally, the results for the flow past a square cylinder show excellent agreement with numerical reference solutions as well as experiments. Copyright © 2015 John Wiley & Sons, Ltd.     

Time‐domain BEM solution of convection–diffusion‐type MHD equations

The two‐dimensional convection–diffusion‐type equations are solved by using the boundary element method (BEM) based on the time‐dependent fundamental solution. The emphasis is given on the solution of magnetohydrodynamic (MHD) duct flow problems with arbitrary wall conductivity. The boundary and time integrals in the BEM formulation are computed numerically assuming constant variations of the unknowns on both the boundary elements and the time intervals. Then, the solution is advanced to the steady‐state iteratively. Thus, it is possible to use quite large time increments and stability problems are not encountered. The time‐domain BEM solution procedure is tested on some convection–diffusion problems and the MHD duct flow problem with insulated walls to establish the validity of the approach. The numerical results for these sample problems compare very well to analytical results. Then, the BEM formulation of the MHD duct flow problem with arbitrary wall conductivity is obtained for the first time in such a way that the equations are solved together with the coupled boundary conditions. The use of time‐dependent fundamental solution enables us to obtain numerical solutions for this problem for the Hartmann number values up to 300 and for several values of conductivity parameter. Copyright © 2007 John Wiley & Sons, Ltd.     

Stabilized finite‐element computation of compressible flow with linear and quadratic interpolation functions

The unsteady compressible flow equations are solved using a stabilized finite‐element formulation with C⁰ elements. In 2D, the performance of three‐noded linear and six‐noded quadratic triangular elements is compared. In 3D, the relative performance is evaluated for 6‐noded linear and 18‐noded quadratic wedge elements. Results are compared for the solutions to Euler, laminar, and turbulent flows at different Mach numbers for several flow problems. The finite‐element meshes considered for comparison have same location of nodes for the linear and quadratic interpolations. For the turbulent flow, the Spalart–Allmaras model is used for closure. It is found that the quadratic elements yield better performance than the linear elements. This is attributed to accurate representation of the stabilization terms that involve second‐order derivatives in the formulation. When these terms are dropped from the formulation with quadratic interpolation, the numerical results are similar to those obtained with linear interpolation. The absence of these terms result in added numerical diffusion that seems to give the effect of a relatively reduced Reynolds number. For the same location of nodes, the computations with the linear triangular and wedge elements are approximately 20% and 100% faster than those with quadratic triangular and wedge elements, respectively. However, if the same quadrature rule for numerical integration is used for both interpolations, the computations with quadratic elements are approximately 20% and 45% faster in 2D and 3D, respectively. Copyright © 2014 John Wiley & Sons, Ltd.     

A high‐order flux reconstruction adaptive mesh refinement method for magnetohydrodynamics on unstructured grids

We report our recent development of the high‐order flux reconstruction adaptive mesh refinement (AMR) method for magnetohydrodynamics (MHD). The resulted framework features a shock‐capturing duo of AMR and artificial resistivity (AR), which can robustly capture shocks and rotational and contact discontinuities with a fraction of the cell counts that are usually required. In our previous paper, 36 we have presented a shock‐capturing framework on hydrodynamic problems with artificial diffusivity and AMR. Our AMR approach features a tree‐free, direct‐addressing approach in retrieving data across multiple levels of refinement. In this article, we report an extension to MHD systems that retains the flexibility of using unstructured grids. The challenges due to complex shock structures and divergence‐free constraint of magnetic field are more difficult to deal with than those of hydrodynamic systems. The accuracy of our solver hinges on 2 properties to achieve high‐order accuracy on MHD systems: removing the divergence error thoroughly and resolving discontinuities accurately. A hyperbolic divergence cleaning method with multiple subiterations is used for the first task. This method drives away the divergence error and preserves conservative forms of the governing equations. The subiteration can be accelerated by absorbing a pseudo time step into the wave speed coefficient, therefore enjoys a relaxed CFL condition. The AMR method rallies multiple levels of refined cells around various shock discontinuities, and it coordinates with the AR method to obtain sharp shock profiles. The physically consistent AR method localizes discontinuities and damps the spurious oscillation arising in the curl of the magnetic field. The effectiveness of the AMR and AR combination is demonstrated to be much more powerful than simply adding AR on finer and finer mesh, since the AMR steeply reduces the required amount of AR and confines the added artificial diffusivity and resistivity to a narrower and narrower region. We are able to verify the designed high‐order accuracy in space by using smooth flow test problems on unstructured grids. The efficiency and robustness of this framework are fully demonstrated through a number of two‐dimensional nonsmooth ideal MHD tests. We also successfully demonstrate that the AMR method can help significantly save computational cost for the Orszag‐Tang vortex problem.     

A multiscale control volume finite element method for advection–diffusion equations

We present a new stabilized method for advection–diffusion equations, which combines a control volume FEM formulation of the governing equations with a novel multiscale approximation of the total flux. The latter incorporates information about the exact solution that cannot be represented on the mesh. To define this flux, we solve the governing equations along suitable mesh segments under the assumption that the flux varies linearly along these segments. This procedure yields second‐order accurate fluxes on the edges of the mesh. Then, we use curl‐conforming elements of the same order to lift these edge fluxes into the mesh elements. In so doing, we obtain a stabilized control volume FEM formulation that is second‐order accurate and does not require mesh‐dependent stabilization parameters. Numerical convergence studies on uniform and nonuniform grids along with several standard advection tests illustrate the computational properties of the new method. Published 2015. This article is a U.S. Government work and is in the public domain in the USA.     

Stabilized finite element formulation for incompressible flow on distorted meshes

Flow computations frequently require unfavourably meshes, as for example highly stretched elements in regions of boundary layers or distorted elements in deforming arbitrary Lagrangian Eulerian meshes. Thus, the performance of a flow solver on such meshes is of great interest. The behaviour of finite elements with residual‐based stabilization for incompressible Newtonian flow on distorted meshes is considered here. We investigate the influence of the stabilization terms on the results obtained on distorted meshes by a number of numerical studies. The effect of different element length definitions within the elemental stabilization parameter is considered. Further, different variants of residual‐based stabilization are compared indicating that dropping the second derivatives from the stabilization operator, i.e. using a streamline upwind Petrov–Galerkin type of formulation yields better results in a variety of cases. A comparison of the performance of linear and quadratic elements reveals further that the inconsistency of linear elements equipped with residual‐based stabilization introduces significant errors on distorted meshes, while quadratic elements are almost unaffected by moderate mesh distortion. Copyright © 2008 John Wiley & Sons, Ltd.     

Discontinuous Galerkin finite element method applied to the coupled unsteady Stokes/Cahn-Hilliard equations

Two-phase flows driven by the interfacial dynamics are studied by tracking implicitly interfaces in the framework of the Cahn-Hilliard theory. The fluid dynamics is described by the Stokes equations with an additional source term in the momentum equation taking into account the capillary forces. A discontinuous Galerkin finite element method is used to solve the coupled Stokes/Cahn-Hilliard equations. The Cahn-Hilliard equation is treated as a system of two coupled equations corresponding to the advection-diffusion equation for the phase field and a nonlinear elliptic equation for the chemical potential. First, the variational formulation of the Cahn-Hilliard equation is presented. A numerical test is achieved showing the optimal order in error bounds. Second, the variational formulation in discontinuous Galerkin finite element approach of the Stokes equations is recalled, in which the same space of approximation is used for the velocity and the pressure with an adequate stabilization technique. The rates of convergence in space and time are evaluated leading to an optimal order in error bounds in space and a second order in time with a backward differentiation formula at the second order. Numerical tests devoted to two-phase flows are provided on ellipsoidal droplet retraction, on the capillary rising of a liquid in a tube, and on the wetting drop over a horizontal solid wall.     

A stabilized finite element method for the Stokes problem based on polynomial pressure projections

A new stabilized finite element method for the Stokes problem is presented. The method is obtained by modification of the mixed variational equation by using local L² polynomial pressure projections. Our stabilization approach is motivated by the inherent inconsistency of equal‐order approximations for the Stokes equations, which leads to an unstable mixed finite element method. Application of pressure projections in conjunction with minimization of the pressure–velocity mismatch eliminates this inconsistency and leads to a stable variational formulation. Unlike other stabilization methods, the present approach does not require specification of a stabilization parameter or calculation of higher‐order derivatives, and always leads to a symmetric linear system. The new method can be implemented at the element level and for affine families of finite elements on simplicial grids it reduces to a simple modification of the weak continuity equation. Numerical results are presented for a variety of equal‐order continuous velocity and pressure elements in two and three dimensions. Copyright © 2004 John Wiley & Sons, Ltd.     

A stabilized mixed formulation for unsteady Brinkman equation based on the method of horizontal lines

In this paper, we present a stabilized mixed formulation for unsteady Brinkman equation. The formulation is systematically derived based on the variational multiscale formalism and the method of horizontal lines. The derivation does not need the assumption that the fine‐scale variables do not depend on the time, which is the case with the conventional derivation of multiscale stabilized formulations for transient mixed problems. An expression for the stabilization parameter is obtained in terms of a bubble function, and appropriate bubble functions for various finite elements are also presented. Under the proposed formulation, equal‐order interpolation for the velocity and pressure (which is computationally the most convenient) is stable. Representative numerical results are presented to illustrate the performance of the proposed formulation. Spatial and temporal convergence studies are also performed, and the proposed formulation performed well. Copyright © 2011 John Wiley & Sons, Ltd.     

Fourth‐order compact formulation of Navier–Stokes equations and driven cavity flow at high Reynolds numbers

A new fourth‐order compact formulation for the steady 2‐D incompressible Navier–Stokes equations is presented. The formulation is in the same form of the Navier–Stokes equations such that any numerical method that solve the Navier–Stokes equations can easily be applied to this fourth‐order compact formulation. In particular, in this work the formulation is solved with an efficient numerical method that requires the solution of tridiagonal systems using a fine grid mesh of 601 × 601. Using this formulation, the steady 2‐D incompressible flow in a driven cavity is solved up to Reynolds number with Re = 20 000 fourth‐order spatial accuracy. Detailed solutions are presented. Copyright © 2005 John Wiley & Sons, Ltd.     

A stabilized finite element formulation for high‐speed inviscid compressible flows using finite calculus

This paper aims at the development of a new stabilization formulation based on the finite calculus (FIC) scheme for solving the Euler equations using the Galerkin FEM on unstructured triangular grids. The FIC method is based on expressing the balance of fluxes in a space–time domain of finite size. It is used to prevent the creation of instabilities typically present in numerical solutions due to the high convective terms and sharp gradients. Two stabilization terms, respectively called streamline term and transverse term, are added via the FIC formulation to the original conservative equations in the space–time domain. An explicit fourth‐order Runge–Kutta scheme is implemented to advance the solution in time. The presented numerical test examples for inviscid flows prove the ability of the proposed stabilization technique for providing appropriate solutions especially near shock waves. Although the derived methodology delivers precise results with a nearly coarse mesh, a mesh refinement technique is coupled to the solution process for obtaining a suitable mesh particularly in the high‐gradient zones. Copyright © 2014 John Wiley & Sons, Ltd.     

A curvature‐free multiphase flow solver via surface stress‐based formulation

In this note, we show the link between the classical continuous surface stress and continuous surface force approaches together with special finite element method techniques toward a fully implicit level set method. Based on a modified surface stress formulation, neither normals nor curvature has to be explicitly calculated. The method is space‐dimension independent. Prototypical numerical tests of benchmarking character for a rising 2D bubble are provided for validating the accuracy of this new approach. We show additionally that the explicit redistancing can be avoided using a nonlinear PDE so that a fully implicit and even monolithic formulation of the corresponding multiphase problem gets feasible.     

Comparison of spectral and finite element methods applied to the study of the core‐annular flow in an undulating tube

A Galerkin/finite element and a pseudo‐spectral method, in conjunction with the primitive (velocity‐pressure) and streamfunction‐vorticity formulations, are tested for solving the two‐phase flow in a tube, which has a periodically varying, circular cross section. Two immiscible, incompressible, Newtonian fluids are arranged so that one of them is around the axis of the tube (core fluid) and the other one surrounds it (annular fluid). The physical and flow parameters are such that the interface between the two fluids remains continuous and single‐valued. This arrangement is usually referred to as Core‐Annular flow. A non‐orthogonal mapping is used to transform the uneven tube shape and the unknown, time dependent interface to fixed, cylindrical surfaces. With both methods and formulations, steady states are calculated first using the Newton–Raphson method. The most dangerous eigenvalues of the related linear stability problem are calculated using the Arnoldi method, and dynamic simulations are carried out using the implicit Euler method. It is shown that with a smooth tube shape the pseudo‐spectral method exhibits exponential convergence, whereas the finite element method exhibits algebraic convergence, albeit of higher order than expected from the relevant theory. Thus the former method, especially when coupled with the streamfunction‐vorticity formulation, is much more efficient. The finite element method becomes more advantageous when the tube shape contains a cusp, in which case the convergence rate of the pseudo‐spectral method deteriorates exhibiting algebraic convergence with the number of the axial spectral modes, whereas the convergence rate of the finite element method remains unaffected. Copyright © 2002 John Wiley & Sons, Ltd.     

A stress‐based least‐squares finite‐element model for incompressible Navier–Stokes equations

In this paper we present a stress‐based least‐squares finite‐element formulation for the solution of the Navier–Stokes equations governing flows of viscous incompressible fluids. Stress components are introduced as independent variables to make the system first order. Continuity equation becomes an algebraic equation and is eliminated from the system with suitable modifications. The h and p convergence are verified using the exact solution of Kovasznay flow. Steady flow past a large circular cylinder in a channel is solved to test mass conservation. Transient flow over a backward‐facing step problem is solved on several meshes. Results are compared with that obtained using vorticity‐based first‐order formulation for both benchmark problems. Copyright © 2007 John Wiley & Sons, Ltd.     

The locally conservative Galerkin (LCG) method for solving the incompressible Navier–Stokes equations

In this paper, the locally conservative Galerkin (LCG) method (Numer. Heat Transfer B Fundam. 2004; 46 :357–370; Int. J. Numer. Methods Eng. 2007) has been extended to solve the incompressible Navier–Stokes equations. A new correction term is also incorporated to make the formulation to give identical results to that of the continuous Galerkin (CG) method. In addition to ensuring element‐by‐element conservation, the method also allows solution of the governing equations over individual elements, independent of the neighbouring elements. This is achieved within the CG framework by breaking the domain into elemental sub‐domains. Although this allows discontinuous trial function field, we have carried out the formulation using the continuous trial function space as the basis. Thus, the changes in the existing CFD codes are kept to a minimum. The edge fluxes, establishing the continuity between neighbouring elements, are calculated via a post‐processing step during the time‐stepping operation. Therefore, the employed formulation needs to be carried out using either a time‐stepping or an equivalent iterative scheme that allows post‐processing of fluxes. The time‐stepping algorithm employed in this paper is based on the characteristic‐based split (CBS) scheme. Both steady‐ and unsteady‐state examples presented show that the element‐by‐element formulation employed is accurate and robust. Copyright © 2007 John Wiley & Sons, Ltd.     

Gauge finite element method for incompressible flows

A finite element method for computing viscous incompressible flows based on the gauge formulation introduced in [Weinan E, Liu J‐G. Gauge method for viscous incompressible flows. Journal of Computational Physics (submitted)] is presented. This formulation replaces the pressure by a gauge variable. This new gauge variable is a numerical tool and differs from the standard gauge variable that arises from decomposing a compressible velocity field. It has the advantage that an additional boundary condition can be assigned to the gauge variable, thus eliminating the issue of a pressure boundary condition associated with the original primitive variable formulation. The computational task is then reduced to solving standard heat and Poisson equations, which are approximated by straightforward, piecewise linear (or higher‐order) finite elements. This method can achieve high‐order accuracy at a cost comparable with that of solving standard heat and Poisson equations. It is naturally adapted to complex geometry and it is much simpler than traditional finite element methods for incompressible flows. Several numerical examples on both structured and unstructured grids are presented. Copyright © 2000 John Wiley & Sons, Ltd.