Entropy dissipation scheme and minimums‐increase‐and‐maximums‐decrease slope limiter

In this paper, we continue to study the entropy dissipation scheme developed in former. We start with a numerical study of the scheme without the entropy dissipation term on the linear advection equation, which shows that the scheme is stable and numerical dissipation and numerical dispersion free for smooth solutions. However, the numerical results for discontinuous solutions show nonlinear instabilities near jump discontinuities. This is because the scheme enforces two related conservation properties in the computation. With this study, we design a so‐called ‘minimums‐increase‐and‐maximums‐decrease’ slope limiter in the reconstruction step of the scheme and delete the entropy dissipation in the linear fields and reduce the entropy dissipation terms in the nonlinear fields. Numerical experiments show improvements of the designed scheme compared with the results presented in former. However, the minimums‐increase‐and‐maximums‐decrease limiter is still not perfect yet, and better slope limiters are still sought. Copyright © 2011 John Wiley & Sons, Ltd.     

A mass‐fraction‐based interface‐capturing method for multi‐component flow

An interface‐capturing method based on mass fraction is developed to solve the Riemann problem in multi‐component compressible flow. Equations of mass fraction with modified form, which is derived from conservative equations of mass, are employed here to capture the interface. By introducing mass fraction into Euler equations system, as well as other conservative coefficients, a quasi‐conservative numerical model is created. Numerical examples show that the mass fraction model performs well not only in multi‐component fluids modeled by simple stiffened gas equation of state (EOS) but also in that modeled by complex Mie–Grüneisen EOS. Moreover, the mass fraction model is applied to Riemann problem with piecewise EOS; the expression of which depends on density. It is found that the mass fraction model can well adapt to the analytic change in piecewise EOS and produce accuracy solutions with fewer unknown quantities, and the model can be easily extended to m‐component fluid mixture by using only m + 4 equations with no additional conditions. Copyright © 2013 John Wiley & Sons, Ltd.     

Time‐linearized time‐harmonic 3‐D Navier–Stokes shock‐capturing schemes

In the present paper, a numerical method for the computation of time‐harmonic flows, using the time‐linearized compressible Reynolds‐averaged Navier–Stokes equations is developed and validated. The method is based on the linearization of the discretized nonlinear equations. The convective fluxes are discretized using an O(Δx) MUSCL scheme with van Leer flux‐vector‐splitting. Unsteady perturbations of the turbulent stresses are linearized using a frozen‐turbulence‐Reynolds‐number hypothesis, to approximate eddy‐viscosity perturbations. The resulting linear system is solved using a pseudo‐time‐marching implicit ADI‐AF (alternating‐directions‐implicit approximate‐factorization) procedure with local pseudo‐time‐steps, corresponding to a matrix‐successive‐underrelaxation procedure. The stability issues associated with the pseudo‐time‐marching solution of the time‐linearized Navier–Stokes equations are discussed. Comparison of computations with measurements and with time‐nonlinear computations for 3‐D shock‐wave oscillation in a square duct, for various back‐pressure fluctuation frequencies (180, 80, 20 and 10 Hz), assesses the shock‐capturing capability of the time‐linearized scheme. Copyright © 2007 John Wiley & Sons, Ltd.     

Interface‐preserving level‐set reinitialization for DG‐FEM

This paper presents a contribution to level‐set reinitialization in the context of discontinuous Galerkin finite element methods. We focus on high‐order polynomials for the discretization and level set geometries, which are comparable to the element size. In contrast to hyperbolic and geometric reinitialization techniques, our method relies on solving a nonlinear elliptic PDE iteratively. We critically compare two different variants of the algorithm experimentally in numerical studies. The results demonstrate that the method is stable for nontrivial test cases and shows high‐order accuracy. Copyright © 2016 John Wiley & Sons, Ltd.     

Depth‐integrated free‐surface flow with a two‐layer non‐hydrostatic formulation

This paper presents the derivation of a depth‐integrated wave propagation and runup model from a system of governing equations for two‐layer non‐hydrostatic flows. The governing equations are transformed into an equivalent, depth‐integrated system, which separately describes the flux‐dominated and dispersion‐dominated processes. The depth‐integrated system reproduces the linear dispersion relation within a 5 error for water depth parameter up to kd = 11, while allowing direct implementation of a momentum conservation scheme to model wave breaking and a moving‐waterline technique for runup calculation. A staggered finite‐difference scheme discretizes the governing equations in the horizontal dimension and the Keller box scheme reconstructs the non‐hydrostatic terms in the vertical direction. An semi‐implicit scheme integrates the depth‐integrated flow in time with the non‐hydrostatic pressure determined from a Poisson‐type equation. The model is verified with solitary wave propagation in a channel of uniform depth and validated with previous laboratory experiments for wave transformation over a submerged bar, a plane beach, and fringing reefs. Copyright © 2011 John Wiley & Sons, Ltd.     

A high‐order fully explicit flux‐form semi‐Lagrangian shallow‐water model

A new approach is proposed for constructing a fully explicit third‐order mass‐conservative semi‐Lagrangian scheme for simulating the shallow‐water equations on an equiangular cubed‐sphere grid. State variables are staggered with velocity components stored pointwise at nodal points and mass variables stored as element averages. In order to advance the state variables in time, we first apply an explicit multi‐step time‐stepping scheme to update the velocity components and then use a semi‐Lagrangian advection scheme to update the height field and tracer variables. This procedure is chosen to ensure consistency between dry air mass and tracers, which is particularly important in many atmospheric chemistry applications. The resulting scheme is shown to be competitive with many existing numerical methods on a suite of standard test cases and demonstrates slightly improved performance over other high‐order finite‐volume models. Copyright © 2014 John Wiley & Sons, Ltd.     

Computer‐algebra multiple‐timescale method for spatially periodic thin‐film viscous‐flow problems

A flexible, fully automated, computer‐algebra algorithm is developed for solving a class of non‐linear partial‐differential evolution equations arising frequently in the modeling of two‐dimensional transient free‐surface viscous thin‐film flows. The method, which is formulated for solving spatially periodic problems, is based upon an explicit multiple‐timescale asymptotic approximation of the thin‐film thickness. It admits the resolution of diverse physical phenomena by employing a finite geometric progression of increasingly slow timescales. The method is implemented on a challenging test problem comprising the evolution of an annular film of viscous liquid, with a free surface, adhering to the exterior of a horizontal rotating circular cylinder; as a model for numerous industrially motivated coating flows, this benchmark problem has been analyzed in diverse numerical and theoretical studies, against whose results those of the present method are compared. The explicit algebraic form of the solution admits a study of large‐time evolutionary dynamics that lies beyond the reach of considerably more expensive conventional numerical solvers, thereby shedding new light on the hitherto‐undiscovered explicit dependence of large‐time evolutionary fluid dynamics in terms of independent parameters describing gravitational and capillary effects. The results obtained from the new computer‐algebra procedure are demonstrated to be in good agreement with those obtained from a bespoke efficient numerical integration method that is spectrally accurate in space and 8th‐order (Runge–Kutta) in time. Newly discovered mechanisms describing the decay of free‐surface wave modes, from arbitrary initial conditions to the steady state, are presented. Copyright © 2010 John Wiley & Sons, Ltd.     

Multi‐size‐mesh multi‐time‐step algorithm for noise computation on curvilinear meshes

Aeroacoustic problems are often multi‐scale and a zonal refinement technique is thus desirable to reduce computational effort while preserving low dissipation and low dispersion errors from the numerical scheme. For that purpose, the multi‐size‐mesh multi‐time‐step algorithm of Tam and Kurbatskii [AIAA Journal, 2000, 38 (8), p. 1331–1339] allows changes by a factor of two between adjacent blocks, accompanied by a doubling in the time step. This local time stepping avoids wasting calculation time, which would result from imposing a unique time step dictated by the smallest grid size for explicit time marching. In the present study, the multi‐size‐mesh multi‐time‐step method is extended to general curvilinear grids by using a suitable coordinate transformation and by performing the necessary interpolations directly in the physical space due to multidimensional interpolations combining order constraints and optimization in the wave number space. A particular attention is paid to the properties of the Adams–Bashforth schemes used for time marching. The optimization of the coefficients by minimizing an error in the wave number space rather than satisfying a formal order is shown to be inefficient for Adams–Bashforth schemes. The accuracy of the extended multi‐size‐mesh multi‐time‐step algorithm is first demonstrated for acoustic propagation on a sinusoidal grid and for a computation of laminar trailing edge noise. In the latter test‐case, the mesh doubling is close to the airfoil and the vortical structures are crossing the doubling interface without affecting the quality of the radiated field. The applicability of the algorithm in three dimensions is eventually demonstrated by computing tonal noise from a moderate Reynolds number flow over an airfoil. Copyright © 2013 John Wiley & Sons, Ltd.     

Sensitivity‐guided decision‐making for wind farm micro‐siting

This paper presents a quantitative risk assessment for design and development of a renewable energy system to support decision‐making among design alternatives. Throughout the decision‐making phases, resources are allocated among exploration and exploitation tasks to manage the uncertainties in design parameters and to adapt designs to new information for enhanced performance. The resource allocation problem is formulated as a sequential decision feedback loop for a quantitative analysis of exploration and exploitation trade‐offs. We support decision‐making by tracking the evolution of uncertainties, the sensitivity of design alternatives to the uncertainties, and the performance, reliability, and robustness of each design. This is achieved by analyzing the uncertainties in the wind resource, the turbine performance and operation, and the models that define the power curve and wake deficiency. Comparison of the performance, reliability, and robustness of aligned and staggered turbine layouts before and after wind assessment experiments aids in improving micro‐siting decisions. The results demonstrate that design decisions can be supported by efficiently allocating resources towards improved estimates of achievable design objectives and by quantitatively assessing the risk in meeting those objectives. Copyright © 2016 John Wiley & Sons, Ltd.     

Positivity‐preserving,flux‐limited finite‐difference and finite‐element methods for reactive transport

A new class of positivity‐preserving, flux‐limited finite‐difference and Petrov–Galerkin (PG) finite‐element methods are devised for reactive transport problems.The methods are similar to classical TVD flux‐limited schemes with the main difference being that the flux‐limiter constraint is designed to preserve positivity for problems involving diffusion and reaction. In the finite‐element formulation, we also consider the effect of numerical quadrature in the lumped and consistent mass matrix forms on the positivity‐preserving property. Analysis of the latter scheme shows that positivity‐preserving solutions of the resulting difference equations can only be guaranteed if the flux‐limited scheme is both implicit and satisfies an additional lower‐bound condition on time‐step size. We show that this condition also applies to standard Galerkin linear finite‐element approximations to the linear diffusion equation. Numerical experiments are provided to demonstrate the behavior of the methods and confirm the theoretical conditions on time‐step size, mesh spacing, and flux limiting for transport problems with and without nonlinear reaction. Copyright © 2003 John Wiley & Sons, Ltd.     

A particle‐in‐cell method for the solution of two‐layer shallow‐water equations

A particle‐in‐cell (PIC) numerical method developed for the study of shallow‐water dynamics, when the moving fluid layer is laterally confined by the intersection of its top and bottom surfaces, is described. The effect of ambient rotation is included for application to geophysical fluids, particularly open‐ocean buoyant vortices in which the underlying density interface outcrops to the surface around the rim of the vortex. Extensions to include the dynamical effect of a second moving layer (baroclinicity) and the presence of a lateral rigid boundary (sidewall) are also described. Although the method was developed for oceanographic investigations, applications to other fluid mechanics problems would be straightforward. Copyright © 2000 John Wiley & Sons, Ltd.     

A dual‐weighted trust‐region adaptive POD 4D‐VAR applied to a finite‐element shallow‐water equations model

We consider a limited‐area finite‐element discretization of the shallow‐water equations model. Our purpose in this paper is to solve an inverse problem for the above model controlling its initial conditions in presence of observations being assimilated in a time interval (window of assimilation). We then attempt to obtain a reduced‐order model of the above inverse problem, based on proper orthogonal decomposition (POD), referred to as POD 4‐D VAR. Different approaches of POD implementation of the reduced inverse problem are compared, including a dual‐weighed method for snapshot selection coupled with a trust‐region POD approach. Numerical results obtained point to an improved accuracy in all metrics tested when dual‐weighing choice of snapshots is combined with POD adaptivity of the trust‐region type. Results of ad‐hoc adaptivity of the POD 4‐D VAR turn out to yield less accurate results than trust‐region POD when compared with high‐fidelity model. Directions of future research are finally outlined. Copyright © 2009 John Wiley & Sons, Ltd.     

Higher‐order‐accurate numerical method for temporal stability simulations of Rayleigh‐Bénard‐Poiseuille flows

For Rayleigh‐Bénard‐Poiseuille flows, thermal stratification resulting from a wall‐normal temperature gradient together with an opposing gravitational field can lead to buoyancy‐driven instability. Moreover, for sufficiently large Reynolds numbers, viscosity‐driven instability can occur. Two higher‐order‐accurate methods based on the full and linearized Navier‐Stokes equations were developed for investigating the temporal stability of such flows. The new methods employ a spectral discretization in the homogeneous directions. In the wall‐normal direction, the convective and viscous terms are discretized with fifth‐order‐accurate biased and fourth‐order‐accurate central compact finite differences. A fourth‐order‐accurate explicit Runge‐Kutta method is employed for time integration. To validate the methods, the primary instability was investigated for different combinations of the Reynolds and Rayleigh number. The results from these primary stability investigations are consistent with linear stability theory results from the literature with respect to both the onset of the instability and the dependence of the temporal growth rate on the wave angle. For the cases with buoyancy‐driven instability, strong linear growth is observed for a broad range of spanwise wavenumbers. The largest growth rates are obtained for a wave angle of 90°. For the cases with viscosity‐driven instability, the linear growth rates are lower and the first mode to experience nonlinear growth is a higher harmonic with half the wavelength of the fundamental.     

High‐order time integrators for front‐tracking finite‐element analysis of viscous free‐surface flows

This paper proposes implicit Runge–Kutta (IRK) time integrators to improve the accuracy of a front‐tracking finite‐element method for viscous free‐surface flow predictions. In the front‐tracking approach, the modeling equations must be solved on a moving domain, which is usually performed using an arbitrary Lagrangian–Eulerian (ALE) frame of reference. One of the main difficulties associated with the ALE formulation is related to the accuracy of the time integration procedure. Indeed, most formulations reported in the literature are limited to second‐order accurate time integrators at best. In this paper, we present a finite‐element ALE formulation in which a consistent evaluation of the mesh velocity and its divergence guarantees satisfaction of the discrete geometrical conservation law. More importantly, it also ensures that the high‐order fixed mesh temporal accuracy of time integrators is preserved on deforming grids. It is combined with the use of a family of L‐stable IRK time integrators for the incompressible Navier–Stokes equations to yield high‐order time‐accurate free‐surface simulations. This is demonstrated in the paper using the method of manufactured solution in space and time as recommended in Verification and Validation. In particular, we report up to fifth‐order accuracy in time. The proposed free‐surface front‐tracking approach is then validated against cases of practical interest such as sloshing in a tank, solitary waves propagation, and coupled interaction between a wave and a submerged cylinder. Copyright © 2015 John Wiley & Sons, Ltd.     

A discontinuous Galerkin‐based sharp‐interface method to simulate three‐dimensional compressible two‐phase flow

A numerical method for the simulation of compressible two‐phase flows is presented in this paper. The sharp‐interface approach consists of several components: a discontinuous Galerkin solver for compressible fluid flow, a level‐set tracking algorithm to follow the movement of the interface and a coupling of both by a ghost‐fluid approach with use of a local Riemann solver at the interface. There are several novel techniques used: the discontinuous Galerkin scheme allows locally a subcell resolution to enhance the interface resolution and an interior finite volume Total Variation Diminishing (TVD) approximation at the interface. The level‐set equation is solved by the same discontinuous Galerkin scheme. To obtain a very good approximation of the interface curvature, the accuracy of the level‐set field is improved and smoothed by an additional P_NP_M‐reconstruction. The capabilities of the method for the simulation of compressible two‐phase flow are demonstrated for a droplet at equilibrium, an oscillating ellipsoidal droplet, and a shock‐droplet interaction problem at Mach 3. Copyright © 2015 John Wiley & Sons, Ltd.     

The multi‐stage centred‐scheme approach applied to a drift‐flux two‐phase flow model

For two‐phase flow models, upwind schemes are most often difficult do derive, and expensive to use. Centred schemes, on the other hand, are simple, but more dissipative. The recently proposed multi‐stage (MUSTA ) method is aimed at coming close to the accuracy of upwind schemes while retaining the simplicity of centred schemes. So far, the MUSTA approach has been shown to work well for the Euler equations of inviscid, compressible single‐phase flow. In this work, we explore the MUSTA scheme for a more complex system of equations: the drift‐flux model, which describes one‐dimensional two‐phase flow where the motions of the phases are strongly coupled. As the number of stages is increased, the results of the MUSTA scheme approach those of the Roe method. The good results of the MUSTA scheme are dependent on the use of a large‐enough local grid. Hence, the main benefit of the MUSTA scheme is its simplicity, rather than CPU ‐time savings. Copyright © 2006 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.     

A σ‐coordinate non‐hydrostatic model with embedded Boussinesq‐type‐like equations for modeling deep‐water waves

A σ‐coordinate non‐hydrostatic model, combined with the embedded Boussinesq‐type‐like equations, a reference velocity, and an adapted top‐layer control, is developed to study the evolution of deep‐water waves. The advantage of using the Boussinesq‐type‐like equations with the reference velocity is to provide an analytical‐based non‐hydrostatic pressure distribution at the top‐layer and to optimize wave dispersion property. The σ‐based non‐hydrostatic model naturally tackles the so‐called overshooting issue in the case of non‐linear steep waves. Efficiency and accuracy of this non‐hydrostatic model in terms of wave dispersion and nonlinearity are critically examined. Overall results show that the newly developed model using a few layers is capable of resolving the evolution of non‐linear deep‐water wave groups. Copyright © 2009 John Wiley & Sons, Ltd.     

CFD simulation of two‐ and three‐dimensional free‐surface flow

The paper describes the implementation of moving‐mesh and free‐surface capabilities within a 3‐d finite‐volume Reynolds‐averaged‐Navier–Stokes solver, using surface‐conforming multi‐block structured meshes. The free‐surface kinematic condition can be applied in two ways: enforcing zero net mass flux or solving the kinematic equation by a finite‐difference method. The free surface is best defined by intermediate control points rather than the mesh vertices. Application of the dynamic boundary condition to the piezometric pressure at these points provides a hydrostatic restoring force which helps to eliminate any unnatural free‐surface undulations. The implementation of time‐marching methods on moving grids are described in some detail and it is shown that a second‐order scheme must be applied in both scalar‐transport and free‐surface equations if flows driven by free‐surface height variations are to be computed without significant wave attenuation using a modest number of time steps. Computations of five flows of theoretical and practical interest—forced motion in a pump, linear waves in a tank, quasi‐1d flow over a ramp, solitary wave interaction with a submerged obstacle and 3‐d flow about a surface‐penetrating cylinder—are described to illustrate the capabilities of our code and methods. Copyright © 2003 John Wiley & Sons, Ltd.