On the usage of NURBS as interface representation in free‐surface flows

When simulating free‐surface flows using the finite element method, there are many cases where the governing equations require information which must be derived from the available discretized geometry. Examples are curvature or normal vectors. The accurate computation of this information directly from the finite element mesh often requires a high degree of refinement—which is not necessarily required to obtain an accurate flow solution. As a remedy and an option to be able to use coarser meshes, the representation of the free surface using non‐uniform rational B‐splines (NURBS) curves or surfaces is investigated in this work. The advantages of a NURBS parameterization in comparison with the standard approach are discussed. In addition, it is explored how the pressure jump resulting from surface tension effects can be handled using doubled interface nodes. Numerical examples include the computation of surface tension in a two‐phase flow as well as the computation of normal vectors as a basis for mesh deformation methods. For these examples, the improvement of the numerical solution compared with the standard approaches on identical meshes is shown. Copyright © 2011 John Wiley & Sons, Ltd.     

Time‐dependent finite element simulations of a shear‐thinning viscoelastic fluid with application to blood flow

A new finite element method is developed to simulate time‐dependent viscoelastic shear‐thinning flows characterized by the generalized Oldroyd‐B model. The focus of the algorithm is improved stability through a free‐energy dissipative scheme by using low‐order piecewise‐constant finite element approximations for stress. The algorithm is further modified by incorporating a pressure‐projection method, a DG‐upwinding scheme, a symmetric interior penalty DG method to solve the elliptic pressure‐update equation and a geometric multigrid preconditioner. The improved stability and cost to accuracy is compared when using higher order discontinuous bilinear approximation, where in addition, we consider the influence of a slope limiter for these elements. The algorithm is applied to the 2D start‐up‐driven cavity problem, and the stability of the free energy is illustrated and compared between element choices. An application of the model to modelling blood in small arterioles and channels is considered by simulating pulsatile blood flow through a stenotic arteriole. The individual influences of viscoelasticity and shear‐thinning within the generalized Oldroyd‐B model are investigated by comparing results to the Newtonian, generalized Newtonian and Oldroyd‐B models. Copyright © 2014 John Wiley & Sons, Ltd.     

Development of a hybrid particle‐mesh method for two‐phase flow simulations

A hybrid particle‐mesh method was developed for efficient and accurate simulations of two‐phase flows. In this method, the main component of the flow is solved using the constrained interpolated profile/multi‐moment finite volumemethod; the two‐phase interface is rendered using the finite volume particle (FVP) method. The effect of surface tension is evaluated using the continuum surface force model. Numerical particles in the FVP method are distributed only on the surface of the liquid in simulating the interface between liquid and gas; these particles are used to determine the density of each mesh grid. An artificial term was also introduced to mitigate particle clustering in the direction of maximum compression and sparse discretization errors in the stretched direction. This enables accurate interface tracking without diminishing numerical efficiency. Two benchmark simulations are used to demonstrate the validity of the method developed and its numerical stability. Copyright © 2016 John Wiley & Sons, Ltd.     

High‐order strand grid methods for low‐speed and incompressible flows

A novel high‐order finite volume scheme using flux correction methods in conjunction with structured finite differences is extended to low Mach and incompressible flows on strand grids. Flux correction achieves a high order by explicitly canceling low‐order truncation error terms across finite volume faces and is applied in unstructured layers of the strand grid. The layers are then coupled together using a source term containing summation‐by‐parts finite differences in the strand direction. A preconditioner is employed to extend the method to low speed and incompressible flows. We further extend the method to turbulent flows with the Spalart–Allmaras model. Laminar flow test cases indicate improvements in accuracy and convergence using the high‐order preconditioned method, while turbulent body‐of‐revolution flow results show improvements in only some cases, perhaps because of dominant errors arising from the turbulence model itself. Copyright © 2016 John Wiley & Sons, Ltd.     

Retracted and replaced: A flow‐condition‐based interpolation finite element procedure for triangular grids

A flow‐condition‐based interpolation finite element scheme is presented for use of triangular grids in the solution of the incompressible Navier–Stokes equations. The method provides spatially isotropic discretizations for low and high Reynolds number flows. Various example solutions are given to illustrate the capabilities of the procedure. This article and been retracted and replaced. See retraction and replacement notice DOI: 10.1002/fld.1247 . Copyright © 2005 John Wiley & Sons, Ltd.     

Finite‐volume multi‐stage schemes for shallow‐water flow simulations

A finite‐volume multi‐stage (FMUSTA) scheme is proposed for simulating the free‐surface shallow‐water flows with the hydraulic shocks. On the basis of the multi‐stage (MUSTA) method, the original Riemann problem is transformed to an independent MUSTA mesh. The local Lax–Friedrichs scheme is then adopted for solving the solution of the Riemann problem at the cell interface on the MUSTA mesh. The resulting first‐order monotonic FMUSTA scheme, which does not require the use of the eigenstructure and the special treatment of entropy fixes, has the generality as well as simplicity. In order to achieve the high‐resolution property, the monotonic upstream schemes for conservation laws (MUSCL) method are used. For modeling shallow‐water flows with source terms, the surface gradient method (SGM) is adopted. The proposed schemes are verified using the simulations of six shallow‐water problems, including the 1D idealized dam breaking, the steady transcritical flow over a hump, the 2D oblique hydraulic jump, the circular dam breaking and two dam‐break experiments. The simulated results by the proposed schemes are in satisfactory agreement with the exact solutions and experimental data. It is demonstrated that the proposed FMUSTA schemes have superior overall numerical accuracy among the schemes tested such as the commonly adopted Roe and HLL schemes. Copyright © 2007 John Wiley & Sons, Ltd.     

Developing a unified FVE‐ALE approach to solve unsteady fluid flow with moving boundaries

In this study, an arbitrary Lagrangian–Eulerian (ALE) approach is incorporated with a mixed finite‐volume–element (FVE) method to establish a novel moving boundary method for simulating unsteady incompressible flow on non‐stationary meshes. The method collects the advantages of both finite‐volume and finite‐element (FE) methods as well as the ALE approach in a unified algorithm. In this regard, the convection terms are treated at the cell faces using a physical‐influence upwinding scheme, while the diffusion terms are treated using bilinear FE shape functions. On the other hand, the performance of ALE approach is improved by using the Laplace method to improve the hybrid grids, involving triangular and quadrilateral elements, either partially or entirely. The use of hybrid FE grids facilitates this achievement. To show the robustness of the unified algorithm, we examine both the first‐ and the second‐order temporal stencils. The accuracy and performance of the extended method are evaluated via simulating the unsteady flow fields around a fixed cylinder, a transversely oscillating cylinder, and in a channel with an indented wall. The numerical results presented demonstrate significant accuracy benefits for the new hybrid method on coarse meshes and where large time steps are taken. Of importance, the current method yields the second‐order temporal accuracy when the second‐order stencil is used to discretize the unsteady terms. Copyright © 2009 John Wiley & Sons, Ltd.     

Pressure‐stabilized maximum‐entropy methods for incompressible Stokes

We present a parameter‐free stable maximum‐entropy method for incompressible Stokes flow. Derived from a least‐biased optimization inspired by information theory, the meshfree maximum‐entropy method appears as an interesting alternative to classical approximation schemes like the finite element method. Especially compared with other meshfree methods, e.g. the moving least‐squares method, it allows for a straightforward imposition of boundary conditions. However, no Eulerian approach has yet been presented for real incompressible flow, encountering the convective and pressure instabilities. In this paper, we exclusively address the pressure instabilities caused by the mixed velocity‐pressure formulation of incompressible Stokes flow. In a preparatory discussion, existing stable and stabilized methods are investigated and evaluated. This is used to develop different approaches towards a stable maximum‐entropy formulation. We show results for two analytical tests, including a presentation of the convergence behavior. As a typical benchmark problem, results are also shown for the leaky lid‐driven cavity. The already presented information‐flux method for convection‐dominated problems in mind, we see this as the last step towards a maximum‐entropy method capable of simulating full incompressible flow problems. Copyright © 2015 John Wiley & Sons, Ltd.     

Finite‐element/level‐set/operator‐splitting (FELSOS) approach for computing two‐fluid unsteady flows with free moving interfaces

The present work is devoted to the study on unsteady flows of two immiscible viscous fluids separated by free moving interface. Our goal is to elaborate a unified strategy for numerical modelling of two‐fluid interfacial flows, having in mind possible interface topology changes (like merger or break‐up) and realistically wide ranges for physical parameters of the problem. The proposed computational approach essentially relies on three basic components: the finite element method for spatial approximation, the operator‐splitting for temporal discretization and the level‐set method for interface representation. We show that the finite element implementation of the level‐set approach brings some additional benefits as compared to the standard, finite difference level‐set realizations. In particular, the use of finite elements permits to localize the interface precisely, without introducing any artificial parameters like the interface thickness; it also allows to maintain the second‐order accuracy of the interface normal, curvature and mass conservation. The operator‐splitting makes it possible to separate all major difficulties of the problem and enables us to implement the equal‐order interpolation for the velocity and pressure. Diverse numerical examples including simulations of bubble dynamics, bifurcating jet flow and Rayleigh–Taylor instability are presented to validate the computational method. Copyright © 2004 John Wiley & Sons, Ltd.     

Simplex space‐time meshes in two‐phase flow simulations

In this paper, we present the numerical solution of 2‐phase flow problems of engineering significance with a space‐time finite element method that allows for local temporal refinement. Arbitrary temporal refinement is applied to preselected regions of the mesh and is governed by a quantity that is part of the solution process, namely, the interface position in 2‐phase flow. Because of local effects such as surface tension, jumps in material properties, etc, the interface can in general be considered a region that requires high flexibility and high resolution, both in space and in time. The new method, which leads to tetrahedral (for 2D problems) and pentatope (for 3D problems) meshes, offers an efficient yet accurate approach to the underlying 2‐phase flow problems.     

Implementation of high‐order compact finite‐difference method to parabolized Navier–Stokes schemes

The numerical solution to the parabolized Navier–Stokes (PNS) and globally iterated PNS (IPNS) equations for accurate computation of hypersonic axisymmetric flowfields is obtained by using the fourth‐order compact finite‐difference method. The PNS and IPNS equations in the general curvilinear coordinates are solved by using the implicit finite‐difference algorithm of Beam and Warming type with a high‐order compact accuracy. A shock‐fitting procedure is utilized in both compact PNS and IPNS schemes to obtain accurate solutions in the vicinity of the shock. The main advantage of the present formulation is that the basic flow variables and their first and second derivatives are simultaneously computed with the fourth‐order accuracy. The computations are carried out for a benchmark case: hypersonic axisymmetric flow over a blunt cone at Mach 8. A sensitivity study is performed for the basic flowfield, including profiles and their derivatives obtained from the fourth‐order compact PNS and IPNS solutions, and the effects of grid size and numerical dissipation term used are discussed. The present results for the flowfield variables and also their derivatives are compared with those of other basic flow models to demonstrate the accuracy and efficiency of the proposed method. The present work represents the first known application of a high‐order compact finite‐difference method to the PNS schemes, which are computationally more efficient than Navier–Stokes solutions. Copyright © 2008 John Wiley & Sons, Ltd.     

A spectral/hp least‐squares finite element analysis of the Carreau–Yasuda fluids

A least‐squares finite element model with spectral/hp approximations was developed for steady, two‐dimensional flows of non‐Newtonian fluids obeying the Carreau–Yasuda constitutive model. The finite element model consists of velocity, pressure, and stress fields as independent variables (hence, called a mixed model). Least‐squares models offer an alternative variational setting to the conventional weak‐form Galerkin models for the Navier–Stokes equations, and no compatibility conditions on the approximation spaces used for the velocity, pressure, and stress fields are necessary when the polynomial order (p) used is sufficiently high (say, p > 3, as determined numerically). Also, the use of the spectral/hp elements in conjunction with the least‐squares formulation with high p alleviates various forms of locking, which often appear in low‐order least‐squares finite element models for incompressible viscous fluids, and accurate results can be obtained with exponential convergence. To verify and validate, benchmark problems of Kovasznay flow, backward‐facing step flow, and lid‐driven square cavity flow are used. Then the effect of different parameters of the Carreau–Yasuda constitutive model on the flow characteristics is studied parametrically. Copyright © 2016 John Wiley & Sons, Ltd.     

An improved unstructured WENO method for compressible multi‐fluids flow

An improved high‐order accurate WENO finite volume method based on unstructured grids for compressible multi‐fluids flow is proposed in this paper. The third‐order accuracy WENO finite volume method based on triangle cell is used to discretize the governing equations. To have higher order of accuracy, the P1 polynomial is reconstructed firstly. After that, the P2 polynomial is reconstructed from the combination of the P1. The reconstructed coefficients are calculated by analytical form of inverse matrix rather than the numerical inversion. This greatly improved the efficiency and the robustness. Four examples are presented to examine this algorithm. Numerical results show that there is no spurious oscillation of velocity and pressure across the interface and high‐order accurate result can be achieved. Copyright © 2016 John Wiley & Sons, Ltd.     

A depth‐integrated non‐hydrostatic finite element model for wave propagation

This paper presents a two‐dimensional finite element model for simulating dynamic propagation of weakly dispersive waves. Shallow water equations including extra non‐hydrostatic pressure terms and a depth‐integrated vertical momentum equation are solved with linear distributions assumed in the vertical direction for the non‐hydrostatic pressure and the vertical velocity. The model is developed based on the platform of a finite element model, CCHE2D. A physically bounded upwind scheme for the advection term discretization is developed, and the quasi second‐order differential operators of this scheme result in no oscillation and little numerical diffusion. The depth‐integrated non‐hydrostatic wave model is solved semi‐implicitly: the provisional flow velocity is first implicitly solved using the shallow water equations; the non‐hydrostatic pressure, which is implicitly obtained by ensuring a divergence‐free velocity field, is used to correct the provisional velocity, and finally the depth‐integrated continuity equation is explicitly solved to satisfy global mass conservation. The developed wave model is verified by an analytical solution and validated by laboratory experiments, and the computed results show that the wave model can properly handle linear and nonlinear dispersive waves, wave shoaling, diffraction, refraction and focusing. Copyright © 2013 John Wiley & Sons, Ltd.     

A finite element thermally coupled flow formulation for phase‐change problems

A finite element, thermally coupled incompressible flow formulation considering phase‐change effects is presented. This formulation accounts for natural convection, temperature‐dependent material properties and isothermal and non‐isothermal phase‐change models. In this context, the full Navier–Stokes equations are solved using a generalized streamline operator (GSO) technique. The highly non‐linear phase‐change effects are treated with a temperature‐based algorithm, which provides stability and convergence of the numerical solution. The Boussinesq approximation is used in order to consider the temperature‐dependent density variation. Furthermore, the numerical solution of the coupled problem is approached with a staggered incremental‐iterative solution scheme, such that the convergence criteria are written in terms of the residual vectors. Finally, this formulation is used for the solutions of solidification and melting problems validating some numerical results with other existing solutions obtained with different methodologies. Copyright © 2000 John Wiley & Sons, Ltd.     

Application of a preconditioned high‐order accurate artificial compressibility‐based incompressible flow solver in wide range of Reynolds numbers

In the present study, the preconditioned incompressible Navier‐Stokes equations with the artificial compressibility method formulated in the generalized curvilinear coordinates are numerically solved by using a high‐order compact finite‐difference scheme for accurately and efficiently computing the incompressible flows in a wide range of Reynolds numbers. A fourth‐order compact finite‐difference scheme is utilized to accurately discretize the spatial derivative terms of the governing equations, and the time integration is carried out based on the dual time‐stepping method. The capability of the proposed solution methodology for the computations of the steady and unsteady incompressible viscous flows from very low to high Reynolds numbers is investigated through the simulation of different 2‐dimensional benchmark problems, and the results obtained are compared with the existing analytical, numerical, and experimental data. A sensitivity analysis is also performed to evaluate the effects of the size of the computational domain and other numerical parameters on the accuracy and performance of the solution algorithm. The present solution procedure is also extended to 3 dimensions and applied for computing the incompressible flow over a sphere. Indications are that the application of the preconditioning in the solution algorithm together with the high‐order discretization method in the generalized curvilinear coordinates provides an accurate and robust solution method for simulating the incompressible flows over practical geometries in a wide range of Reynolds numbers including the creeping flows.     

High‐order numerical simulations of flow‐induced noise

In this paper, the flow/acoustics splitting method for predicting flow‐generated noise is further developed by introducing high‐order finite difference schemes. The splitting method consists of dividing the acoustic problem into a viscous incompressible flow part and an inviscid acoustic part. The incompressible flow equations are solved by a second‐order finite volume code EllipSys2D/3D. The acoustic field is obtained by solving a set of acoustic perturbation equations forced by flow quantities. The incompressible pressure and velocity form the input to the acoustic equations. The present work is an extension of our acoustics solver, with the introduction of high‐order schemes for spatial discretization and a Runge–Kutta scheme for time integration. To achieve low dissipation and dispersion errors, either Dispersion‐Relation‐Preserving (DRP) schemes or optimized compact finite difference schemes are used for the spatial discretizations. Applications and validations of the new acoustics solver are presented for benchmark aeroacoustic problems and for flow over an NACA 0012 airfoil. Copyright © 2010 John Wiley & Sons, Ltd.     

Computation of two‐dimensional flows past ram‐air parachutes

Computational results for flow past a two‐dimensional model of a ram‐air parachute with leading edge cut are presented. Both laminar (Re=10⁴) and turbulent (Re=10⁶) flows are computed. A well‐proven stabilized finite element method (FEM), which has been applied to various flow problems earlier, is utilized to solve the incompressible Navier–Stokes equations in the primitive variables formulation. The Baldwin–Lomax model is employed for turbulence closure. Turbulent flow computations past a Clarck‐Y airfoil without a leading edge cut, for α=7.5°, result in an attached flow. The leading edge cut causes the flow to become unsteady and leads to a significant loss in lift and an increase in drag. The flow inside the parafoil cell remains almost stagnant, resulting in a high value of pressure, which is responsible for giving the parafoil its shape. The value of the lift‐to‐drag ratio obtained with the present computations is in good agreement with those reported in the literature. The effect of the size and location of the leading edge cut is studied. It is found that the flow on the upper surface of the parafoil is fairly insensitive to the configuration of the cut. However, the flow quality on the lower surface improves as the leading edge cut becomes smaller. The lift‐to‐drag ratio for various configurations of the leading edge cut varies between 3.4 and 5.8. It is observed that even though the time histories of the aerodynamic coefficients from the laminar and turbulent flow computations are quite different, their time‐averaged values are quite similar. Copyright © 2001 John Wiley & Sons, Ltd.     

An efficient edge‐based level set finite element method for free surface flow problems

We present an efficient technique for the solution of free surface flow problems using level set and a parallel edge‐based finite element method. An unstructured semi‐explicit solution scheme is proposed. A custom data structure, obtained by blending node‐based and edge‐based approaches is presented so to allow a good parallel performance. In addition to standard velocity extrapolation (for the convection of the level set function), an explicit extrapolation of the pressure field is performed in order to impose both the pressure boundary condition and the volume conservation. The latter is also improved with a modification of the divergence free constrain. The method is shown to allow an efficient solution of both simple benchmark cases and complex industrial examples. Copyright © 2012 John Wiley & Sons, Ltd.