Stability and consistency of nonhydrostatic free‐surface models using the semi‐implicit θ‐method

The θ‐method is a popular semi‐implicit finite‐difference method for simulating free‐surface flows. Problem stiffness, arising because of the presence of both fast and slow timescale processes, is easily handled by the θ‐method. In most ocean, coastal, and estuary modeling applications, stiffness is caused by fast surface gravity wave timescales imposed on slower timescales of baroclinic variability. The method is well known to be unconditionally stable for shallow water (hydrostatic) models when , where θ is the implicitness parameter. In this paper, we demonstrate that the method is also unconditionally stable for nonhydrostatic models, when for both pressure projection and pressure correction methods. However, the methods result in artificial damping of the barotropic mode. In addition to investigating stability, we also estimate the form of artificial damping induced by both the free surface and nonhydrostatic pressure solution methods. Finally, this analysis may be used to estimate the damping or growth associated with a particular wavenumber and the overall order of accuracy of the discretization. Copyright © 2012 John Wiley & Sons, Ltd.     

Building generalized open boundary conditions for fluid dynamics problems

This paper deals with the design of an efficient open boundary condition (OBC) for fluid dynamics problems. Such problematics arise, for instance, when one solves a local model on a fine grid that is nested in a coarser one of greater extent. Usually, the local solution U^loc is computed from the coarse solution U^ext, thanks to an OBC formulated as , where B_h and B_H are discretizations of the same differential operator (B_h being defined on the fine grid and B_H on the coarse grid). In this paper, we show that such an OBC cannot lead to the exact solution, and we propose a generalized formulation , where g is a correction term. When B_h and B_H are discretizations of a transparent operator, g can be computed analytically, at least for simple equations. Otherwise, we propose to approximate g by a Richardson extrapolation procedure. Numerical test cases on a 1D Laplace equation and on a 1D shallow water system illustrate the improved efficiency of such a generalized OBC compared with usual ones. Copyright © 2012 John Wiley & Sons, Ltd.     

A novel interface‐tracking method based on Lagrangian particles for deformation analysis of a red blood cell in a capillary

A novel efficient interface‐tracking method is developed to gain an insight into the interface in a multiphase or multifluid system, called the modified particle binary level set (MPBLS) method, in which the binary level set function is defined to distinguish the different phases or fluids and further modified by Lagrangian particles scattered along the interface for achieving higher accuracy. The validation of the MPBLS method is carried out first by simulating the free motion of a red blood cell (RBC) in the rotating, shear and Poiseuille flows, respectively. Subsequently, further validations are performed by comparing with the experimental and numerical results published previously. As one of important applications, the MPBLS method is employed to investigate the deformation behaviors of RBCs with different shapes in a capillary. The simulations show that the healthy RBC gradually changes the geometric shape from a biconcave to a steady parachute shape. It is thus guaranteed that the RBC successfully traverses through the smaller capillaries compared with undeformed RBC. However, the unhealthy RBC with the circular or elliptical shape has different deformation behaviors, in which the steady parachute shape is much less concave at the rear and more convex in the front. Copyright © 2011 John Wiley & Sons, Ltd.     

An edge-based pressure stabilization technique for finite elements on arbitrarily anisotropic meshes

In this paper, we analyze a stabilized equal-order finite element approximation for the Stokes equations on anisotropic meshes. In particular, we allow arbitrary anisotropies in a subdomain, for example, along the boundary of the domain, with the only condition that a maximum angle is fulfilled in each element. This discretization is motivated by applications on moving domains as arising, for example, in fluid-structure interaction or multiphase-flow problems. To deal with the anisotropies, we define a modification of the original continuous interior penalty stabilization approach. We show analytically the discrete stability of the method and convergence of order in the energy norm and in the L²-norm of the velocities. We present numerical examples for a linear Stokes problem and for a nonlinear fluid-structure interaction problem, which substantiate the analytical results and show the capabilities of the approach.     

Fast Ewald summation for Stokesian particle suspensions

We present a numerical method for suspensions of spheroids of arbitrary aspect ratio, which sediment under gravity. The method is based on a periodized boundary integral formulation using the Stokes double layer potential. The resulting discrete system is solved iteratively using generalized minimal residual accelerated by the spectral Ewald method, which reduces the computational complexity to , where N is the number of points used to discretize the particle surfaces. We develop predictive error estimates, which can be used to optimize the choice of parameters in the Ewald summation. Numerical tests show that the method is well conditioned and provides good accuracy when validated against reference solutions. Copyright © 2014 John Wiley & Sons, Ltd.     

A correction for Reynolds-averaged-Navier–Stokes turbulence model under the effect of shock waves in hypersonic flows

Studies on the unphysical increase of turbulent quantities for RANS simulation induced by shock waves in hypersonic flows are carried out. Numerical experiments on the hypersonic flow over a blunt body reveal that the phenomenon of unphysical increase of turbulent quantities across the detached shock wave is induced by the strain-rate-based production terms of the k-$\omega$ and k-$\omega$ SST turbulence models, which leads to the over-prediction of aerothermal prediction. While this phenomenon does not occur for Spalart–Allmaras (S–A) turbulence model because of its vorticity-based production term. In order to eliminate this unphysical phenomenon, and to maintain the accuracy of the original models for boundary layer and separation flows, a new correction method for the k-$\omega$ and k-$\omega$ SST models is proposed: by comparing the orders of magnitude between the strain-rate-based and vorticity-based production terms, the vorticity-based production term is used near the shock waves, while the original strain-rate-based production term is still used in other regions. Finally, the correction method is applied to turbulence and transition flows over blunt bodies, and the numerical results show that the correction method effectively eliminates the unphysical increase of turbulent quantities across shock waves and improves the accuracy of aerothermal and transition onset location prediction.     

A high‐order PIC method for advection‐dominated flow with application to shallow water waves

A hybrid Eulerian‐Lagrangian particle‐in‐cell–type numerical method is developed for the solution of advection‐dominated flow problems. Particular attention is given over to the high‐order transfer of flow properties from the particles to the grid. For smooth flows, the method presented is of formal high‐order accuracy in space. The method is applied to solve the nonlinear shallow water equations resulting in a new, and novel, shock capturing shallow water solver. The approach is able to simulate complex shallow water flows, which can contain an arbitrary number of discontinuities. Both trivial and nontrivial bottom topography is considered, and it is shown that the new scheme is inherently well balanced, exactly satisfying the ‐property. The scheme is verified against several one‐dimensional benchmark shallow water problems. These include cases that involve transcritical flow regimes, shock waves, and nontrivial bathymetry. In all the test cases presented, very good results are obtained.     

Axisymmetric boundary integral formulation for a two‐fluid system

A 3D axisymmetric Galerkin boundary integral formulation for potential flow is employed to model two fluids of different densities, one fluid enclosed inside the other. The interface variables are the velocity potential and the normal velocity, and they can be solved for separately, the second linear system being symmetric. The algorithm is validated by comparing with the analytic solutions for a static interior spherical drop over a range of values for the relative densities of exterior and interior fluids and various boundary conditions. For time‐dependent simulations utilizing a level set method for the interface tracking, the accuracy has been checked by comparing against the known oscillation frequency of the sphere. Pinch‐off profiles corresponding to an initial two‐lobe geometry drop and D = 6 are also presented. Published in 2011 by John Wiley & Sons, Ltd.     

A new σ‐transform based Fourier‐Legendre‐Galerkin model for nonlinear water waves

This paper presents a new spectral model for solving the fully nonlinear potential flow problem for water waves in a single horizontal dimension. At the heart of the numerical method is the solution to the Laplace equation which is solved using a variant of the $\sigma$‐transform. The method discretizes the spatial part of the governing equations using the Galerkin method and the temporal part using the classical fourth‐order Runge‐Kutta method. A careful investigation of the numerical method's stability properties is carried out, and it is shown that the method is stable up to a certain threshold steepness when applied to nonlinear monochromatic waves in deep water. Above this threshold artificial damping may be employed to obtain stable solutions. The accuracy of the model is tested for: (i) highly nonlinear progressive wave trains, (ii) solitary wave reflection, and (iii) deep water wave focusing events. In all cases it is demonstrated that the model is capable of obtaining excellent results, essentially up to very near breaking.     

An efficient and stable finite element solver of higher order in space and time for nonstationary incompressible flow

In this paper, we present fully implicit continuous Galerkin–Petrov (cGP) and discontinuous Galerkin (dG) time‐stepping schemes for incompressible flow problems which are, in contrast to standard approaches like for instance the Crank–Nicolson scheme, of higher order in time. In particular, we analyze numerically the higher order dG(1) and cGP(2) methods, which are super convergent of third, resp., fourth order in time, whereas for the space discretization, the well‐known LBB‐stable finite element pair of third‐order accuracy is used. The discretized systems of nonlinear equations are treated by using the Newton method, and the associated linear subproblems are solved by means of a monolithic (geometrical) multigrid method with a blockwise Vanka‐like smoother treating all components simultaneously. We perform nonstationary simulations (in 2D) for two benchmarking configurations to analyze the temporal accuracy and efficiency of the presented time discretization schemes w.r.t. CPU and numerical costs. As a first test problem, we consider a classical ‘flow around cylinder’ benchmark. Here, we concentrate on the nonstationary behavior of the flow patterns with periodic oscillations and examine the ability of the different time discretization schemes to capture the dynamics of the flow. As a second test case, we consider the nonstationary ‘flow through a Venturi pipe’. The objective of this simulation is to control the instantaneous and mean flux through this device. Copyright © 2013 John Wiley & Sons, Ltd.     

A stabilized finite element method for transient Navier–Stokes equations based on two local Gauss integrations

On the basis of two local Gauss integrations, a stabilized finite element method for transient Navier–Stokes equations is presented, which is defined by the lowest equal‐order conforming finite element subspace such as (or ) elements. The best algorithmic feature of our method is using two local Gauss integrations to replace projection operator. The diffusion term in these equations is discretized by using finite element method, and the temporal differentiation and advection terms are treated by characteristic schemes. Moreover, we present some numerical simulations to demonstrate the effectiveness, good stability, and accuracy properties of our method. Especially, the rate of convergence study tells us that the stability still keeps well when the Reynolds number is increasing. Copyright © 2011 John Wiley & Sons, Ltd.     

Comparison of POD reduced order strategies for the nonlinear 2D shallow water equations

This paper introduces tensorial calculus techniques in the framework of POD to reduce the computational complexity of the reduced nonlinear terms. The resulting method, named tensorial POD, can be applied to polynomial nonlinearities of any degree p. Such nonlinear terms have an online complexity of , where k is the dimension of POD basis and therefore is independent of full space dimension. However, it is efficient only for quadratic nonlinear terms because for higher nonlinearities, POD model proves to be less time consuming once the POD basis dimension k is increased. Numerical experiments are carried out with a two‐dimensional SWE test problem to compare the performance of tensorial POD, POD, and POD/discrete empirical interpolation method (DEIM). Numerical results show that tensorial POD decreases by 76× the computational cost of the online stage of POD model for configurations using more than 300,000 model variables. The tensorial POD SWE model was only 2 to 8× slower than the POD/DEIM SWE model but the implementation effort is considerably increased. Tensorial calculus was again employed to construct a new algorithm allowing POD/DEIM SWE model to compute its offline stage faster than POD and tensorial POD approaches. Copyright © 2014 John Wiley & Sons, Ltd.     

Fujiwhara interaction of tropical cyclone scale vortices using a weighted residual collocation method

The fundamental interaction between tropical cyclones was investigated through a series of water tank experiements by Fujiwhara  6 - 8 . However, a complete understanding of tropical cyclones remains an open research challenge although there have been numerous investigations through measurments with aircrafts/satellites, as well as with numerical simulations. This article presents a computational model for simulating the interaction between cyclones. The proposed numerical method is presented briefly, where the time integration is performed by projecting the discrete system onto a Krylov subspace. The method filters the large scale fluid dynamics using a multiresolution approximation, and the unresolved dynamics are modeled with a Smagorinsky type subgrid scale parameterization scheme. Numerical experiments with Fujiwhara interactions are considered to verify modeling accuracy. An excellent agreement between the present simulation and a reference simulation at has been demonstrated. At , the kinetic energy of cyclones is seen consolidated into larger scales with concurrent enstrophy cascade – suggesting a steady increase of energy containing scales – a phenomena that is typical in two‐dimensional turbulence theory. The primary results of this article suggest a novel avenue for addressing some of the computational challenges of mesoscale atmospheric circulations. Copyright © 2015 John Wiley & Sons, Ltd.     

Hessian‐based model reduction: large‐scale inversion and prediction

Hessian‐based model reduction was previously proposed as an approach in deriving reduced models for the solution of large‐scale linear inverse problems by targeting accuracy in observation outputs. A control‐theoretic view of Hessian‐based model reduction that hinges on the equality between the Hessian and the transient observability gramian of the underlying linear system is presented. The model reduction strategy is applied to a large‐scale ( degrees of freedom) three‐dimensional contaminant transport problem in an urban environment, an application that requires real‐time computation. In addition to the inversion accuracy, the ability of reduced models of varying dimension to make predictions of the contaminant evolution beyond the time horizon of observations is studied. Results indicate that the reduced models have a factor speedup in computing time for the same level of accuracy. Copyright © 2012 John Wiley & Sons, Ltd.     

A computational methodology for two‐dimensional fluid flows

A weighted residual collocation methodology for simulating two‐dimensional shear‐driven and natural convection flows has been presented. Using a dyadic mesh refinement, the methodology generates a basis and a multiresolution scheme to approximate a fluid flow. To extend the benefits of the dyadic mesh refinement approach to the field of computational fluid dynamics, this article has studied an iterative interpolation scheme for the construction and differentiation of a basis function in a two‐dimensional mesh that is a finite collection of rectangular elements. We have verified that, on a given mesh, the discretization error is controlled by the order of the basis function. The potential of this novel technique has been demonstrated with some representative examples of the Poisson equation. We have also verified the technique with a dynamical core of a two‐dimensional flow in primitive variables. An excellent result has been observed—on resolving a shear layer and on the conservation of the potential and the kinetic energies—with respect to previously reported benchmark simulations. In particular, the shear‐driven simulation at CFL = 2.5 (Courant–Friedrichs–Lewy) and (Reynolds number) exhibits a linear speed up of CPU time with an increase of the time step, Δt. For the natural convection flow, the conversion of the potential energy to the kinetic energy and the conservation of total energy is resolved by the proposed method. The computed streamlines and the velocity fields have been demonstrated. Copyright © 2014 John Wiley & Sons, Ltd.     

Numerical simulation of the transient shape of the red blood cell in microcapillary flow

The transient shape of a red blood cell (RBC) in a microcapillary flow is simulated under different initial conditions, including various axis orientations and centroid locations, using the LBM-DLM/FD method, which is derived from the lattice Boltzmann method and the distributed Lagrange multiplier/fictitious domain method. Although the terminal velocity is not sensitive to the initial configuration, the evolution of the velocity and the shape are determined by the initial conditions. The parachute and the slipper shape are the most probable shapes for a deformed RBC in the flow. An RBC with an initial axis orientation of 90 degrees exhibits a more complicated deformation. RBCs have a tendency to move to the centerline of a tube if an offset between the RBC centroid and the centerline exists. Our numerical results are validated by experiments, and some details beyond the experiment are provided.     

A dissipative particle dynamics study of flow in periodically grooved nanochannels

The effect of periodic rectangular wall roughness on planar nanochannel flow is investigated by dissipative particle dynamics simulation. The wall protrusion length is varied, and its effect on the flow is examined. Analysis of particle trajectories and average residence time reveals temporary trapping of fluid particles inside the rectangular cavities for a considerable amount of time. This trapping affects the density, velocity, pressure, and temperature distribution inside and close to the cavities. Inside the cavities, low‐velocity regions and regions of high density related to high pressure and high temperature are observed. When compared with that of the channel with flat walls case, lower flow velocities, temperatures, and pressures are observed for grooved channels. The reduction of the above quantities is more pronounced as the protrusion length, that is, the roughness characteristic length, decreases. Finally, the relation of friction factor, f, with the flow Reynolds number is discussed. The model predicts = constant in the range . The results of this work are of direct relevance to the design of nanofluidic devices. Copyright © 2011 John Wiley & Sons, Ltd.     

An augmented velocity–vorticity–pressure formulation for the Brinkman equations

This paper deals with the analysis of a new augmented mixed finite element method in terms of vorticity, velocity, and pressure, for the Brinkman problem with nonstandard boundary conditions. The approach is based on the introduction of Galerkin least‐squares terms arising from the constitutive equation relating the aforementioned unknowns and from the incompressibility condition. We show that the resulting augmented bilinear form is continuous and elliptic, which, thanks to the Lax–Milgram theorem, and besides proving the well‐posedness of the continuous formulation, ensures the solvability and stability of the Galerkin scheme with any finite element subspace of the continuous space. In particular, Raviart–Thomas elements of any order for the velocity field, and piecewise continuous polynomials of degree k + 1 for both the vorticity and the pressure, can be utilized. A priori error estimates and the corresponding rates of convergence are also given here. Next, we derive two reliable and efficient residual‐based a posteriori error estimators for this problem. The ellipticity of the bilinear form together with the local approximation properties of the Clément interpolation operator are the main tools for showing the reliability. In turn, inverse inequalities and the localization technique based on triangle‐bubble and edge‐bubble functions are utilized to show the efficiency. Finally, several numerical results illustrating the good performance of the method, confirming the properties of the estimators and showing the behavior of the associated adaptive algorithms, are reported. Copyright © 2015 John Wiley & Sons, Ltd.     

Convergence of five-point scheme with diagonally dominant characteristic for DLR k−ε model and its application

In this paper, the equivalent equations of DLR k−ε turbulent model in the boundary-fitted curvilinear coordinate are employed. Using the upwind idea that the contribution of the difference coefficients to the main node is positive contribution, and the other nodes are negative contribution or no contribution, new five-point difference schemes with a full diagonally dominant coefficient matrix (5-point-DD difference scheme) are constructed. Finally, taking the $\bar{u}$ equation in the DLR k−ε turbulent model as an example, the mathematical characteristics of the 5-point-DD difference scheme are analyzed, and the uniform boundedness and convergence theorems of the Gauss-Seidel iterative sequence are given. Numerical simulations show that the five-point schemes are strictly diagonally dominant, and the calculated results are in good agreement with the experimental results.     

A consistent and fast weakly compressible smoothed particle hydrodynamics with a new wall boundary condition

A modified weakly compressible smoothed particle hydrodynamics (WCSPH) is presented, which utilizes consistent discretization schemes for spatial derivatives in the flow equations. Here, each SPH particle is considered as a computational point that represents a specific part of the fluid. To overcome non‐physical oscillations that usually arise in standard WCSPH, we modified the mass conservation equation by using a numerical filter. This modification is based on the difference between two discretization schemes used for the term . Furthermore, a new implementation of wall boundary condition in SPH is introduced. This condition is imposed on the pressure of wall boundary particles to ensure that the acceleration of each boundary particle in normal direction to the wall is zero. Thus, no penetration through walls will occur. To examine the performance of the modified method, we solved a series of two‐dimensional incompressible internal flow benchmark problems. By comparing the result with analytical solutions and the results of the standard WCSPH, we show that the use of consistent schemes in conjunction with the proposed numerical filter improves both accuracy and speed of the numerical method. Copyright © 2011 John Wiley & Sons, Ltd.