Fast spectral solutions of the double-gyre problem in a turbulent flow regime

Several semi-analytical models are considered for a double-gyre problem in a turbulent flow regime for which a reference fully numerical eddy-resolving solution is obtained. The semi-analytical models correspond to solving the depth-averaged Navier–Stokes equations using the spectral Galerkin approach. The robustness of the linear and Smagorinsky eddy-viscosity models for turbulent diffusion approximation is investigated. To capture essential properties of the double-gyre configuration, such as the integral kinetic energy, the integral angular momentum, and the jet mean-flow distribution, an improved semi-analytical model is suggested that is inspired by the idea of scale decomposition between the jet and the surrounding flow.     

Direct solution of Navier–Stokes equations by radial basis functions

The pressure–velocity formulation of the Navier–Stokes (N–S) equation is solved using the radial basis functions (RBF) collocation method. The non-linear collocated equations are solved using the Levenberg–Marquardt method. The primary novelty of this approach is that the N–S equation is solved directly, instead of using an iterative algorithm for the primitive variables. Two flow situations are considered: Couette flow with and without pressure gradient, and 2D laminar flow in a duct with and without flow obstruction. The approach is validated by comparing the Couette flow results with the analytical solution and the 2D results with those obtained using the well-validated CFD-ACE™ commercial package.     

Block mesh refinement for incompressible flows in curvilinear domains

A method for the solution of the Navier–Stokes equation for the prediction of flows inside domains of arbitrary shaped bounds by the use of Cartesian grids with block-refinement in space is presented. In order to avoid the complexity of the body fitted numerical grid generation procedure, we use a saw tooth method for the curvilinear geometry approximation. By using block-nested refinement, we achieved the desired geometry Cartesian approximation in order to find an accurate solution of the N–S equations. The method is applied to incompressible laminar flows and is based on a cell-centred approximation. We present the numerical simulation of the flow field for two geometries, driven cavity and stenosed tubes. The utility of the algorithm is tested by comparing the convergence characteristics and accuracy to those of the standard single grid algorithm. The Cartesian block refinement algorithm can be used in any complex curvilinear geometry simulation, to accomplish a reduction in memory requirements and the computational time effort.     

High-order approximation of Pearson diffusion processes

This paper focuses on Pearson diffusions and the spectral high-order approximation of their related Fokker–Planck equations. The Pearson diffusions is a class of diffusions defined by linear drift and quadratic squared diffusion coefficient. They are widely used in the physical and chemical sciences, engineering, rheology, environmental sciences and financial mathematics. In recent years diffusion models have been studied analytically and numerically primarily through the solution of stochastic differential equations. Analytical solutions have been derived for some of the Pearson diffusions, including the Ornstein–Uhlenbeck, Cox–Ingersoll–Ross and Jacobi processes. However, analytical investigations and computations for diffusions with so-called heavy-tailed ergodic distributions are more difficult to perform. The novelty of this research is the development of an accurate and efficient numerical method to solve the Fokker–Planck equations associated with Pearson diffusions with different boundary conditions. Comparisons between the numerical predictions and available time-dependent and equilibrium analytical solutions are made. The solution of the Fokker–Planck equation is approximated using a reduced basis spectral method. The advantage of this approach is that many models for pricing options in financial mathematics cannot be expressed in terms of a stochastic partial differential equation and therefore one has to resort to solving Fokker–Planck type equations.     

Hybrid LES-RANS Modeling of Complex Turbulent Flows

The paper describes a state-of-the-art hybrid LES-URANS method for the simulation of complex internal and external turbulent flows. Relying on a unified LES-URANS approach with a soft interface the methodology is designed for wall-bounded non-equilibrium flows. The unsteady Reynolds-averaged Navier-Stokes (URANS) mode within the hybrid approach is taken into account by an explicit algebraic Reynolds stress model (EARSM), which guarantees an appropriate representation of the anisotropic near-wall turbulence. All non-closed terms in the transport equation of the turbulent kinetic energy are modeled by enhanced formulations using the EARSM (production and diffusion term) and the splitting of the dissipation rate into a homogeneous and an inhomogeneous contribution. The former is expressed analytically by a Taylor series expansion of the homogeneous lateral Taylor microscale in the vicinity of the wall guaranteeing the correct asymptotic behavior. The latter is incorporated into the diffusion term. The interface location between the large-eddy simulation (LES) mode and the URANS mode is determined automatically on-the-fly based on the modeled turbulent kinetic energy and the distance to the wall. For transitional (external) flows an additional dynamic transition criterion is applied which determines the laminar and the turbulent flow regimes and is combined with the existing interface criterion. An internal flow over a periodic arrangement of hills and an external flow past a SD7003 airfoil with a laminar separation bubble are taken into account for a detailed evaluation of the method. (© 2014 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

A 3D non-linear k–ε turbulent model for prediction of flow and mass transport in channel with vegetation

The results from a 3D non-linear k–ε turbulence model with vegetation are presented to investigate the flow structure, the velocity distribution and mass transport process in a straight compound open channel and a curved open channel. The 3D numerical model for calculating flow is set up in non-orthogonal curvilinear coordinates in order to calculate the complex boundary channel. The finite volume method is used to disperse the governing equations and the SIMPLEC algorithm is applied to acquire the coupling of velocity and pressure. The non-linear k–ε turbulent model has good useful value because of taking into account the anisotropy and not increasing the computational time. The water level of this model is determined from 2D Poisson equation derived from 2D depth-averaged momentum equations. For concentration simulation, an expression for dispersion through vegetation is derived in the present work for the mixing due to flow over vegetation. The simulated results are in good agreement with available experimental data, which indicates that the developed 3D model can predict the flow structure and mass transport in the open channel with vegetation.     

PSE在边界层中预测转捩位置的应用

层流到湍流的转捩是自然界和各项工程实践中广泛存在的现象,层流和湍流的性质大不相同．因此,预测转捩位置是流体力学中的重要理论和实际问题．针对不可压缩边界层,入口加入展向等幅值型和展向波包型两类扰动,展向等幅值型扰动是由一个二维波(2-D)和两个三维波(3-D)组成,使用抛物化稳定性方程(PSE)的方法来研究扰动的演化和预测转捩位置,并且与数值模拟的结果相比较．结果表明,PSE可以研究扰动的演化和预测转捩位置,同时其计算比数值模拟快得多．     

Numerical simulation of thermal problems coupled with magnetohydrodynamic effects in aluminium cell

A system of partial differential equations describing the thermal behavior of aluminium cell coupled with magnetohydrodynamic effects is numerically solved. The thermal model is considered as a two-phases Stefan problem which consists of a non-linear convection–diffusion heat equation with Joule effect as a source. The magnetohydrodynamic fields are governed by Navier–Stokes and by static Maxwell equations. A pseudo-evolutionary scheme (Chernoff) is used to obtain the stationary solution giving the temperature and the frozen layer profile for the simulation of the ledges in the cell. A numerical approximation using a finite element method is formulated to obtain the fluid velocity, electrical potential, magnetic induction and temperature. An iterative algorithm and 3-D numerical results are presented.     

Depth-averaged modeling of free surface flows in open channels with emerged and submerged vegetation

The resistance induced by vegetation on the flow in a watercourse should be considered in projects of watercourse management and river restoration. Depth-averaged numerical model is an efficient tool to study this problem. In this study, a depth-averaged model using the finite volume method on a staggered curvilinear grid and the SIMPLEC algorithm for numerical solution is developed for simulating the hydrodynamics of free surface flows in watercourses with vegetation. For the model formulation the vegetation resistance is treated as a momentum sink and represented by a Manning type equation, and turbulence is parameterized by the k–ε equations. An analytical equation is derived to represent the resistance induced by submerged vegetation by an equivalent Manning roughness coefficient. Numerical simulation is carried out for the flow in an open channel with a 180° bend, and the flow in a curved open channel partly covered by emerged vegetation, as well as the flow in a straight trapezoidal channel with submerged vegetation. The agreement between the computed results and the measured data is generally good, showing that the resistance due to emerged or submerged vegetation can be represented accurately by the Manning roughness equation. The computed results demonstrate that the depth-averaged modeling is a reasonable and efficient tool to study flows in watercourses with vegetations.     

Coloured noise for dispersion of contaminants in shallow waters

In this article, we explore the application of a set of stochastic differential equations called particle model in simulating the advection and diffusion of pollutants in shallow waters. The Fokker–Planck equation associated with this set of stochastic differential equations is interpreted as an advection–diffusion equation. This enables us to derive an underlying particle model that is exactly consistent with the advection–diffusion equation. Still, neither the advection–diffusion equation nor the related traditional particle model accurately takes into account the short-term spreading behaviour of particles. To improve the behaviour of the model shortly after the deployment of contaminants, a particle model forced by a coloured noise process is developed in this article. The use of coloured noise as a driving force unlike Brownian motion, enables to us to take into account the short-term correlated turbulent fluid flow velocity of the particles. Furthermore, it is shown that for long-term simulations of the dispersion of particles, both the particle due to Brownian motion and the particle model due to coloured noise are consistent with the advection–diffusion equation.     

Simulation of planar ionization wave front propagation on an unstructured adaptive grid

The paper deals with the numerical solution of a basic 2D model of the propagation of an ionization wave. The system of equations describing this propagation consists of a coupled set of reaction–diffusion-convection equations and a Poissons equation. The transport equations are solved by a finite volume method on an unstructured triangular adaptive grid. The upwind scheme and the diamond scheme are used for the discretization of the convection and diffusion fluxes, respectively. The Poisson equation is also discretized by the diamond scheme. Numerical results are presented. We deal in more detail with numerical tests of the grid adaptation technique and its influence on the numerical results. An original behavior is observed. The grid refinement is not sufficient to obtain accurate results for this particular phenomenon. Using a second order scheme for convection is necessary.     

Numerical simulations of water–gas flow in heterogeneous porous media with discontinuous capillary pressures by the concept of global pressure

We present an approach and numerical results for a new formulation modeling immiscible compressible two-phase flow in heterogeneous porous media with discontinuous capillary pressures. The main feature of this model is the introduction of a new global pressure, and it is fully equivalent to the original equations. The resulting equations are written in a fractional flow formulation and lead to a coupled degenerate system which consists of a nonlinear parabolic (the global pressure) equation and a nonlinear diffusion–convection one (the saturation equation) with nonlinear transmission conditions at the interfaces that separate different media. The resulting system is discretized using a vertex-centred finite volume method combined with pressure and flux interface conditions for the treatment of heterogeneities. An implicit Euler approach is used for time discretization. A Godunov-type method is used to treat the convection terms, and the diffusion terms are discretized by piecewise linear conforming finite elements. We present numerical simulations for three one-dimensional benchmark tests to demonstrate the ability of the method to approximate solutions of water–gas equations efficiently and accurately in nuclear underground waste disposal situations.     

A fourth-order compact ADI method for solving two-dimensional unsteady convection–diffusion problems

In this article, an exponential high-order compact (EHOC) alternating direction implicit (ADI) method, in which the Crank–Nicolson scheme is used for the time discretization and an exponential fourth-order compact difference formula for the steady-state 1D convection–diffusion problem is used for the spatial discretization, is presented for the solution of the unsteady 2D convection–diffusion problems. The method is temporally second-order accurate and spatially fourth order accurate, which requires only a regular five-point 2D stencil similar to that in the standard second-order methods. The resulting EHOC ADI scheme in each ADI solution step corresponds to a strictly diagonally dominant tridiagonal matrix equation which can be inverted by simple tridiagonal Gaussian decomposition and may also be solved by application of the one-dimensional tridiagonal Thomas algorithm with a considerable saving in computing time. The unconditionally stable character of the method was verified by means of the discrete Fourier (or von Neumann) analysis. Numerical examples are given to demonstrate the performance of the method proposed and to compare mostly it with the high order ADI method of Karaa and Zhang and the spatial third-order compact scheme of Note and Tan.     

A computational meshless method for the generalized Burger’s–Huxley equation

A numerical solution of the generalized Burger’s–Huxley equation, based on collocation method using Radial basis functions (RBFs), called Kansa’s approach is presented. The numerical results are compared with the exact solution, Adomian decomposition method (ADM) and Variational iteration method (VIM). Highly accurate and efficient results are obtained by RBFs method. Excellent agreement with the exact solution is observed while better (or same) accuracy is obtained than other numerical schemes cited in this work.     

Velocity and concentration profiles in uniform sediment-laden flow

This paper presents a theoretical model for computing the velocity and sediment concentration profiles in a uniform sediment-laden flow carrying all fine, medium and coarse sediments. The proposed model essentially includes the effect of sediment concentration in total turbulent shear stress and eddy diffusivity in addition to the modified mixing length derived by Umeyama and Gerritsen [J. Hydr. Engng., ASCE, 118 (2) (1992) 229–245] applied to Hunt’s diffusion equation. Numerical solution of coupled differential equations for velocity and sediment concentration is carried out. The theoretical results show quite good agreement with the experimental data.     

Weighted quadrature rules for finite element methods

We discuss the numerical integration of polynomials times non-polynomial weighting functions in two dimensions arising from multiscale finite element computations. The proposed quadrature rules are significantly more accurate than standard quadratures and are better suited to existing finite element codes than formulas computed by symbolic integration. We validate this approach by introducing the new quadrature formulas into a multiscale finite element method for the two-dimensional reaction–diffusion equation.     

Parameter estimation in convection dominated nonlinear convection–diffusion problems by the relaxation method and the adjoint equation

The development of numerical methods for strongly nonlinear convection–diffusion problems with dominant convection is an ongoing topic in numerical analysis. For inverse problems in this setting, there is a need of fast and accurate solvers. Here, we present operator splitting with a Riemann solver for the convective part and a relaxation method for the diffusive part, as a means to achieve this goal. Combined with the adjoint equation method this allows us to solve inverse problems within reasonable time frames and with modest computing power. As an example, the dual-well experiment is considered and the adjoint method is compared with a conjugate gradient algorithm and a Levenberg–Marquardt type of iteration method.     

Direct numerical simulation of the laminar–turbulent transition at hypersonic flow speeds on a supercomputer

A method for direct numerical simulation of three-dimensional unsteady disturbances leading to a laminar–turbulent transition at hypersonic flow speeds is proposed. The simulation relies on solving the full three-dimensional unsteady Navier–Stokes equations. The computational technique is intended for multiprocessor supercomputers and is based on a fully implicit monotone approximation scheme and the Newton–Raphson method for solving systems of nonlinear difference equations. This approach is used to study the development of three-dimensional unstable disturbances in a flat-plate and compression-corner boundary layers in early laminar–turbulent transition stages at the free-stream Mach number M = 5.37. The three-dimensional disturbance field is visualized in order to reveal and discuss features of the instability development at the linear and nonlinear stages. The distribution of the skin friction coefficient is used to detect laminar and transient flow regimes and determine the onset of the laminar–turbulent transition.     

Cross-relaxation dynamics on anomalous saturation processes in low pressure supersonic cw HF chemical laser amplifier

It is assumed that both translational and rotational nonequilibrium cross-relaxations play a role simultaneoulsy in low pressure supersonic cw HF chemical laser amplifier. For two-type models of gas flow medium with laminar and turbulent flow diffusion mixing, the expressions of saturated gain spectrum are derived respectively, and the numerical calculations are performed as well. The numerical results show that turbulent flow diffusion mixing model is in the best agreement with the experimental result. Project supported by the National Natural Science Foundation of China (Grant No. 19474036) and by the Laboratory of Hih-temperature Gas Dynamics. Insititute of Mechanics, Chinese Academy of Sciences (Grant No. KJ951-E-202).     

Numerical solution of 2D and 3D backward facing step flows

The work deals with numerical modelling of flow through 2-dimensional (2D) and 3-dimensional (3D) backward facing step. In laminar case, we apply several higher order upwind and central discretizations and compare numerical results with measurements. The turbulent regime is considered in 2D as well as in 3D and influence of secondary flow is observed. Different modifications of low-Re two equation turbulence models and an explicit algebraic Reynolds stress model (EARSM) are considered. (© 2005 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim)