高精度有限差分法模拟Kelvin-Helmholtz不稳定性

 WENO有限差分格式有较高的分辨精度,适合复杂流场的计算,在国际上被广泛采用。本文利用WENO有限差分格式求解2维守恒型欧拉方程,实现了对无粘流体中Kelvin-Helmholtz不稳定性的数值模拟。速度剪切方向采用周期边界条件;扰动增长方向采用嵌边出流边界条件,一个不稳定波长分布64个网格。数值模拟给出的扰动幅值线性增长率与线性稳定性分析给出的结果很好符合,显示了该格式的有效性和精度。数值模拟给出了清晰的密度等值线,表明该方法还具有较好的界面变形捕捉能力。     

高精度有限差分法模拟Kelvin-Helmholtz不稳定性

WENO有限差分格式有较高的分辨精度,适合复杂流场的计算,在国际上被广泛采用。本文利用WENO有限差分格式求解2维守恒型欧拉方程,实现了对无粘流体中Kelvin-Helmholtz不稳定性的数值模拟。速度剪切方向采用周期边界条件;扰动增长方向采用嵌边出流边界条件,一个不稳定波长分布64个网格。数值模拟给出的扰动幅值线性增长率与线性稳定性分析给出的结果很好符合,显示了该格式的有效性和精度。数值模拟给出了清晰的密度等值线,表明该方法还具有较好的界面变形捕捉能力。     

A mass-conserved multiphase lattice Boltzmann method based on high-order difference

The Z–S–C multiphase lattice Boltzmann model [Zheng, Shu, and Chew(ZSC), J. Comput. Phys. 218, 353(2006)]is favored due to its good stability, high efficiency, and large density ratio. However, in terms of mass conservation, this model is not satisfactory during the simulation computations. In this paper, a mass correction is introduced into the ZSC model to make up the mass leakage, while a high-order difference is used to calculate the gradient of the order parameter to improve the accuracy. To verify the improved model, several three-dimensional multiphase flow simulations are carried out,including a bubble in a stationary flow, the merging of two bubbles, and the bubble rising under buoyancy. The numerical simulations show that the results from the present model are in good agreement with those from previous experiments and simulations. The present model not only retains the good properties of the original ZSC model, but also achieves the mass conservation and higher accuracy.     

A shock-detecting sensor for filtering of high-order compact finite difference schemes

A new shock-detecting sensor for properly switching between a second-order and a higher-order filter is developed and assessed. The sensor is designed based on an order analysis. The nonlinear filter with the proposed sensor ensures damping of the high-frequency waves in smooth regions and at the same time removes the Gibbs oscillations around the discontinuities when using high-order compact finite difference schemes. In addition, a suitable scaling is proposed to have dissipation proportional to the shock strength and also to minimize the effects of the second-order filter on the very small scales. Several numerical experiments are carried out and the accuracy of the nonlinear filter with the proposed sensor is examined. In addition, some comparisons with other filters and sensors are made.     

Stochastic finite difference lattice Boltzmann method for steady incompressible viscous flows

With the advent of state-of-the-art computers and their rapid availability, the time is ripe for the development of efficient uncertainty quantification (UQ) methods to reduce the complexity of numerical models used to simulate complicated systems with incomplete knowledge and data. The spectral stochastic finite element method (SSFEM) which is one of the widely used UQ methods, regards uncertainty as generating a new dimension and the solution as dependent on this dimension. A convergent expansion along the new dimension is then sought in terms of the polynomial chaos system, and the coefficients in this representation are determined through a Galerkin approach. This approach provides an accurate representation even when only a small number of terms are used in the spectral expansion; consequently, saving in computational resource can be realized compared to the Monte Carlo (MC) scheme. Recent development of a finite difference lattice Boltzmann method (FDLBM) that provides a convenient algorithm for setting the boundary condition allows the flow of Newtonian and non-Newtonian fluids, with and without external body forces to be simulated with ease. Also, the inherent compressibility effect in the conventional lattice Boltzmann method, which might produce significant errors in some incompressible flow simulations, is eliminated. As such, the FDLBM together with an efficient UQ method can be used to treat incompressible flows with built in uncertainty, such as blood flow in stenosed arteries. The objective of this paper is to develop a stochastic numerical solver for steady incompressible viscous flows by combining the FDLBM with a SSFEM. Validation against MC solutions of channel/Couette, driven cavity, and sudden expansion flows are carried out.     

A second-order accurate immersed boundary method for fully resolved simulations of particle-laden flows

An immersed boundary method (IBM) with second-order spatial accuracy is presented for fully resolved simulations of incompressible viscous flows laden with rigid particles. The method is based on the computationally efficient direct-forcing method of Uhlmann [M. Uhlmann, An immersed boundary method with direct forcing for simulation of particulate flows, J. Comput. Phys. 209 (2005) 448–476] that is embedded in a finite-volume/pressure-correction method. The IBM consists of two grids: a fixed uniform Eulerian grid for the fluid phase and a uniform Lagrangian grid attached to and moving with the particles. A regularized delta function is used to communicate between the two grids and proved to be effective in suppressing grid locking. Without significant loss of efficiency, the original method is improved by: (1) a better approximation of the no-slip/no-penetration (ns/np) condition on the surface of the particles by a multidirect forcing scheme, (2) a correction for the excess in the effective particle diameter by a slight retraction of the Lagrangian grid from the surface towards the interior of the particles with a fraction of the Eulerian grid spacing, and (3) an enhancement of the numerical stability for particle–fluid mass density ratios near unity by a direct account of the inertia of the fluid contained within the particles. The new IBM contains two new parameters: the number of iterations N_s of the multidirect forcing scheme and the retraction distance r_d. The effect of N_s and r_d on the accuracy is investigated for five different flows. The results show that r_d has a strong influence on the effective particle diameter and little influence on the error in the ns/np condition, while exactly the opposite holds for N_s. A novel finding of this study is the demonstration that r_d has a strong influence on the order of grid convergence. It is found that for spheres the choice of r_d = 0.3Δx yields second-order accuracy compared to first-order accuracy of the original method that corresponds to r_d = 0. Finally, N_s = 2 appears optimal for reducing the error in the ns/np condition and maintaining the computational efficiency of the method.     

A front-tracking method for computation of interfacial flows with soluble surfactants

A finite-difference/front-tracking method is developed for computations of interfacial flows with soluble surfactants. The method is designed to solve the evolution equations of the interfacial and bulk surfactant concentrations together with the incompressible Navier–Stokes equations using a non-linear equation of state that relates interfacial surface tension to surfactant concentration at the interface. The method is validated for simple test cases and the computational results are found to be in a good agreement with the analytical solutions. The method is then applied to study the cleavage of drop by surfactant—a problem proposed as a model for cytokinesis [H.P. Greenspan, On the dynamics of cell cleavage, J. Theor. Biol. 65(1) (1977) 79; H.P. Greenspan, On fluid-mechanical simulations of cell division and movement, J. Theor. Biol., 70(1) (1978) 125]. Finally the method is used to model the effects of soluble surfactants on the motion of buoyancy-driven bubbles in a circular tube and the results are found to be in a good agreement with available experimental data.     

Pressure based finite volume method for calculation of compressible viscous gas flows

A pressure based, iterative finite volume method is developed for calculation of compressible, viscous, heat conductive gas flows at all speeds. The method does not need the use of under-relaxation coefficient in order to ensure a convergence of the iterative process. The method is derived from a general form of system of equations describing the motion of compressible, viscous gas. An emphasis is done on the calculation of gaseous microfluidic problems. A fast transient process of gas wave propagation in a two-dimensional microchannel is used as a benchmark problem. The results obtained by using the new method are compared with the numerical solution obtained by using SIMPLE (iterative) and PISO (non-iterative) methods. It is shown that the new iterative method is faster than SIMPLE. For the considered problem the new method is slightly faster than PISO as well. Calculated are also some typical microfluidic subsonic and supersonic flows, and the Rayleigh–Bénard convection of a rarefied gas in continuum limit. The numerical results are compared with other analytical and numerical solutions.     

A finite element variational multiscale method for incompressible flows based on two local gauss integrations

In this article, we present a finite element variational multiscale (VMS) method for incompressible flows based on two local Gauss integrations, and compare it with common VMS method which is defined by a low order finite element space L_h$L_h$ on the same grid as X_h$X_h$ for the velocity deformation tensor and a stabilization parameter α$\alpha$. The best algorithmic feature of our method is using two local Gauss integrations to replace projection operator. We theoretically discuss the relationship between our method and common VMS method for the Taylor–Hood elements, and show that the nonlinear system derived from our method by finite element discretization is much smaller than that of common VMS method computationally.     

Geometric conservation law and applications to high-order finite difference schemes with stationary grids

The geometric conservation law (GCL) includes the volume conservation law (VCL) and the surface conservation law (SCL). Though the VCL is widely discussed for time-depending grids, in the cases of stationary grids the SCL also works as a very important role for high-order accurate numerical simulations. The SCL is usually not satisfied on discretized grid meshes because of discretization errors, and the violation of the SCL can lead to numerical instabilities especially when high-order schemes are applied. In order to fulfill the SCL in high-order finite difference schemes, a conservative metric method (CMM) is presented. This method is achieved by computing grid metric derivatives through a conservative form with the same scheme applied for fluxes. The CMM is proven to be a sufficient condition for the SCL, and can ensure the SCL for interior schemes as well as boundary and near boundary schemes. Though the first-level difference operators δ₃ have no effects on the SCL, no extra errors can be introduced as δ₃ = δ₂. The generally used high-order finite difference schemes are categorized as central schemes (CS) and upwind schemes (UPW) based on the difference operator δ₁ which are used to solve the governing equations. The CMM can be applied to CS and is difficult to be satisfied by UPW. Thus, it is critical to select the difference operator δ₁ to reduce the SCL-related errors. Numerical tests based on WCNS-E-5 show that the SCL plays a very important role in ensuring free-stream conservation, suppressing numerical oscillations, and enhancing the robustness of the high-order scheme in complex grids.     

A finite difference method with non-uniform timesteps for fractional diffusion and diffusion-wave equations

An implicit finite difference method with non-uniform timesteps for solving fractional diffusion and diffusion-wave equations in the Caputo form is presented. The non-uniformity of the timesteps allows one to adapt their size to the behaviour of the solution, which leads to large reductions in the computational time required to obtain the numerical solution without loss of accuracy. The stability of the method has been proved recently for the case of diffusion equations; for diffusion-wave equations its stability, although not proven, has been checked through extensive numerical calculations.     

Numerical stability analysis for the explicit high-order finite difference analysis of rotationally symmetric shells

A numerical stability analysis has been formulated to accompany the already developed explicit high-order finite difference analysis of rotationally symmetric shells subjected to time-dependent impulsive loadings. This already developed analysis utilizes a constant nodal point spacing for the spatial finite difference mesh, with the governing field differential equations formulated in terms of the transverse, meridional, and circumferential displacements as the fundamental variables. The remaining quantities which enter into the natural boundary conditions at each edge of the shell are incorporated into the complete system of equations by defining those quantities at each boundary in terms of the displacements. Surface loadings and inertia forces in each of the three displacement directions of the shell have been considered in the governing equations. Ordinary finite difference representations are used for the time derivatives. All loadings and dependent variables in the circumferential direction of the shell are expressed in Fourier series expansions. The complete system of equations is solved implicitly for the first time increment, while explicit relations are used to determine the three primary displacements within the boundary edges of the shell for the second and succeeding time increments. Separate implicit solutions at each boundary are then used to determine the remaining unspecified primary variables on and outside the boundaries. Subsequently, the remaining primary variables within the boundary edges of the shell and all secondary variables are determined explicitly. Numerical stability (or instability) of numerical solutions for given choices of spatial and time increments is determined by evaluation of the eigenvalues of the explicit coefficient matrix and comparing the maximum eigenvalue with the requirements of a stability criterion developed before by the author. Solutions for typical shells and loadings together with results of stability analyses have been included, and comparisons of the stability requirements and solutions with the requirements and solutions based upon ordinary spatial finite difference representations are included.     

Stabilized finite element method for incompressible flows with high Reynolds number

In the following paper, we discuss the exhaustive use and implementation of stabilization finite element methods for the resolution of the 3D time-dependent incompressible Navier–Stokes equations. The proposed method starts by the use of a finite element variational multiscale (VMS) method, which consists in here of a decomposition for both the velocity and the pressure fields into coarse/resolved scales and fine/unresolved scales. This choice of decomposition is shown to be favorable for simulating flows at high Reynolds number. We explore the behaviour and accuracy of the proposed approximation on three test cases. First, the lid-driven square cavity at Reynolds number up to 50,000 is compared with the highly resolved numerical simulations and second, the lid-driven cubic cavity up to Re = 12,000 is compared with the experimental data. Finally, we study the flow over a 2D backward-facing step at Re = 42,000. Results show that the present implementation is able to exhibit good stability and accuracy properties for high Reynolds number flows with unstructured meshes.     

A finite difference method of solving anisotropic scattering problems

A new method of solving radiative transfer problems is described including a comparison of its speed with that of the doubling method, and a discussion of its accuracy and suitability for computations involving variable optical properties. The method uses a discretization in angle to produce a coupled set of first-order differential equations which are integrated between discrete depth points to produce a set of recursion relations for symmetric and anti-symmetric angular sums of the radiation field at alternate depth points. The formulation given here includes depth-dependent anisotropic scattering, absorption, and internal sources, and allows arbitrary combinations of specular and non-Lambertian diffuse reflection at either or both boundaries. The method is shown to be faster than the doubling method when the number of depth points and angular quadrature points is identical. Numerical tests of the method show that it can return accurate emergent intensities even for large optical depths. The method is also shown to conserve flux to machine accuracy in conservative atmospheres. Finally, several checks are made that demonstrate that the new method can compute accurate radiation fields in atmospheres with variable optical properties.     

基于大规模并行计算的三维多群中子扩散方程有限差分方法

三维多群中子扩散方程的精确、高效求解是核动力堆芯设计及燃料管理的基础。应用有限差分方法求解该方程具有简便、精确、成熟的优点;然而,该方法的计算量和存储量均较大,极大地限制了它的计算规模和应用范围。本文基于大规模并行计算,研究三维多群中子扩散方程有限差分方法：采用中心有限差分格式离散中子扩散方程;基于MPI并行编程模型,采用空间区域分解的方式实现大规模并行计算;采用多群多区域耦合PGMRES算法进行并行加速。在集群服务器上开发了ParaFiDi程序,并采用IAEA3D,PHWR等多个基准题对该程序进行验证。数值结果表明,ParaFiDi程序具有较高的计算精度和计算效率。     

基于大规模并行计算的三维多群中子扩散方程有限差分方法

三维多群中子扩散方程的精确、高效求解是核动力堆芯设计及燃料管理的基础。应用有限差分方法求解该方程具有简便、精确、成熟的优点;然而,该方法的计算量和存储量均较大,极大地限制了它的计算规模和应用范围。本文基于大规模并行计算,研究三维多群中子扩散方程有限差分方法:采用中心有限差分格式离散中子扩散方程;基于MPI并行编程模型,采用空间区域分解的方式实现大规模并行计算;采用多群多区域耦合PGMRES算法进行并行加速。在集群服务器上开发了ParaFiDi程序,并采用IAEA3D,PHWR等多个基准题对该程序进行验证。数值结果表明,ParaFiDi程序具有较高的计算精度和计算效率。     

A high-order low-Mach number AMR construction for chemically reacting flows

A high-order projection scheme was developed for the study of chemically reacting flows in the low-Mach number limit. The numerical approach for the momentum transport uses a combination of cell-centered/cell-averaged discretizations to achieve a fourth order formulation for the pressure projection algorithm. This scheme is coupled with a second order in time operator-split stiff approach for the species and energy equations. The code employs a fourth order, block-structured, adaptive mesh refinement approach to address the challenges posed by the large spectrum of spatial scales encountered in reacting flow computations. Results for advection–diffusion-reaction configurations are used to illustrate the performance of the numerical construction.     

A new high-order discontinuous Galerkin spectral finite element method for Lagrangian gas dynamics in two-dimensions

This paper presents a new high-order cell-centered Lagrangian scheme for two-dimensional compressible flow. The scheme uses a fully Lagrangian form of the gas dynamics equations, which is a weakly hyperbolic system of conservation laws. The system of equations is discretized in the Lagrangian space by discontinuous Galerkin method using a spectral basis. The vertex velocities and the numerical fluxes through the cell interfaces are computed consistently in the Eulerian space by virtue of an improved nodal solver. The nodal solver uses the HLLC approximate Riemann solver to compute the velocities of the vertex. The time marching is implemented by a class of TVD Runge–Kutta type methods. A new HWENO (Hermite WENO) reconstruction algorithm is developed and used as limiters for RKDG methods to maintain compactness of RKDG methods. The scheme is conservative for the mass, momentum and total energy. It can maintain high-order accuracy both in space and time, obey the geometrical conservation law, and achieve at least second order accuracy on quadrilateral meshes. Results of some numerical tests are presented to demonstrate the accuracy and the robustness of the scheme.     

A general finite difference time domain method for hybrid dispersive media model

The shift operator finite difference time domain (SO-FDTD) method is developed to deal with a hybrid dispersive media model. First, we prove that the complex permittivity of the hybrid dispersive media model can be described by a rational polynomial fraction with respect to jω. Then, the relation between electric displacement D and electric field strength E is derived in the time domain by introducing z_t as a shift operator. The constitutive relation in the discretized time domain and the recursive formulation of D and E available for FDTD computation are obtained. Finally, the reflection of the hybrid dispersive slab is computed. The computed results are in good agreement with that obtained by analytic method. This illustrates the generalization and the feasibility of the presented scheme.