面向杆系结构的离散元异步长计算方法

为了提升离散元法处理连续介质问题的计算效率, 提出一种基于重叠颗粒的离散元分区异步长计算方法来处理连续介质问题。该方法采用颗粒分割将连续介质划分为若干子分区, 各分区采用velocity-Verlet积分格式求解运动方程。相邻分区通过重叠颗粒构成局部耦合区域, 边界数据传递过程中不涉及插值和截断过程。数值算例表明, 该方法在保证高精度的同时有效地降低了计算时间。在实际应用中可以有针对性地将连续介质划分为不同尺度颗粒的分区, 根据问题规模及分区颗粒尺度特性采用不同时间步长, 节省存储空间且大幅提升计算效率。     

An adaptive implicit–explicit scheme for the DNS and LES of compressible flows on unstructured grids

An adaptive implicit–explicit scheme for Direct Numerical Simulation (DNS) and Large-Eddy Simulation (LES) of compressible turbulent flows on unstructured grids is developed. The method uses a node-based finite-volume discretization with Summation-by-Parts (SBP) property, which, in conjunction with Simultaneous Approximation Terms (SAT) for imposing boundary conditions, leads to a linearly stable semi-discrete scheme. The solution is marched in time using an Implicit–Explicit Runge–Kutta (IMEX-RK) time-advancement scheme. A novel adaptive algorithm for splitting the system into implicit and explicit sets is developed. The method is validated using several canonical laminar and turbulent flows. Load balance for the new scheme is achieved by a dual-constraint, domain decomposition algorithm. The scalability and computational efficiency of the method is investigated, and memory savings compared with a fully implicit method is demonstrated. A notable reduction of computational costs compared to both fully implicit and fully explicit schemes is observed.     

二维抛物型问题的Strang型交替分段区域分裂格式

给出求解二维抛物型方程的Strang型的交替分段区域分裂格式。交替分段思想可以将区域分为一些不重叠的子区域，Strang型算子分裂技巧通过将高维问题的求解分解为几个低维问题的求解来降低其求解的复杂度。方法是无条件稳定的，理论分析了截断误差。数值算例说明格式的有效性及时空的二阶精度.     

A new compact scheme for parallel computing using domain decomposition

Direct numerical simulation (DNS) of complex flows require solving the problem on parallel machines using high accuracy schemes. Compact schemes provide very high spectral resolution, while satisfying the physical dispersion relation numerically. However, as shown here, compact schemes also display bias in the direction of convection – often producing numerical instability near the inflow and severely damping the solution, always near the outflow. This does not allow its use for parallel computing using domain decomposition and solving the problem in parallel in different sub-domains. To avoid this, in all reported parallel computations with compact schemes the full domain is treated integrally, while using parallel Thomas algorithm (PTA) or parallel diagonal dominant (PDD) algorithm in different processors with resultant latencies and inefficiencies. For domain decomposition methods using compact scheme in each sub-domain independently, a new class of compact schemes is proposed and specific strategies are developed to remove remaining problems of parallel computing. This is calibrated here for parallel computing by solving one-dimensional wave equation by domain decomposition method. We also provide the error norm with respect to the wavelength of the propagated wave-packet. Next, the advantage of the new compact scheme, on a parallel framework, has been shown by solving three-dimensional unsteady Navier–Stokes equations for flow past a cone-cylinder configuration at a Mach number of 4.Additionally, a test case is conducted on the advection of a vortex for a subsonic case to provide an estimate for the error and parallel efficiency of the method using the proposed compact scheme in multiple processors.     

An efficient multi-scale Poisson solver for the incompressible Navier–Stokes equations with immersed boundaries

The iterative-multi-scale-finite-volume (IMSFV) procedure is applied as an efficient solver for the pressure Poisson equation arising in numerical methods for the simulation of incompressible flows with the immersed-interface method (IIM). Motivated by the requirements of the specific IIM implementation, a modified version of the IMSFV algorithm is presented to allow the solution of problems, in which the varying coefficient of the elliptic equation (e.g. the permeability of the medium in the context of the simulation of flows in porous media) varies over several orders of magnitude or even becomes zero within the integration domain. Furthermore, a strategy is proposed to incorporate the iterative procedure needed by the IIM to converge out constraints at immersed boundaries into the iterative IMSFV cycle. No significant deterioration of performance of the IMSFV method is observed with respect to cases, in which no iterative improvement of the boundary conditions is considered.     

Ghost Fluid方法与双介质可压缩流动计算

应用带有Isobaric修正的GhostFluid方法配合LevelSet方法计算可压缩双介质无粘流动.该方法可以消除计算流体界面时所产生的数值跳动和耗散,且编程上比界面跟踪法简单.应用WENO格式数值求解欧拉方程和LevelSet方程,对由刚性气体状态方程所支配的一二维双介质流动进行数值计算,得到了分辨率较高的计算结果.     

一种求解可压缩Navier-Stokes方程的隐式迭代时间推进方法

针对有壁面边界的可压缩流动问题,提出与基于非等距网格的高精度紧致型差分格式相结合的简化隐式迭代时间推进法,建立求解可压缩Navier-Stokes方程的直接数值模拟方法,提高了计算效率.应用该方法,直接数值模拟两种有壁面边界的二维可压缩流动问题,即可压缩平板边界层流动和可压缩槽道流动.     

A hierarchical fracture model for the iterative multiscale finite volume method

An iterative multiscale finite volume (i-MSFV) method is devised for the simulation of multiphase flow in fractured porous media in the context of a hierarchical fracture modeling framework. Motivated by the small pressure change inside highly conductive fractures, the fully coupled system is split into smaller systems, which are then sequentially solved. This splitting technique results in only one additional degree of freedom for each connected fracture network appearing in the matrix system. It can be interpreted as an agglomeration of highly connected cells; similar as in algebraic multigrid methods. For the solution of the resulting algebraic system, an i-MSFV method is introduced. In addition to the local basis and correction functions, which were previously developed in this framework, local fracture functions are introduced to accurately capture the fractures at the coarse scale. In this multiscale approach there exists one fracture function per network and local domain, and in the coarse scale problem there appears only one additional degree of freedom per connected fracture network. Numerical results are presented for validation and verification of this new iterative multiscale approach for fractured porous media, and to investigate its computational efficiency. Finally, it is demonstrated that the new method is an effective multiscale approach for simulations of realistic multiphase flows in fractured heterogeneous porous media.     

A Fast High Order Iterative Solver for the Electromagnetic Scattering by Open Cavities Filled with the Inhomogeneous Media

The scattering of the open cavity filled with the inhomogeneous media is studied. The problem is discretized with a fourth order finite difference scheme and the immersed interface method, resulting in a linear system of equations with the high order accurate solutions in the whole computational domain. To solve the system of equations, we design an efficient iterative solver, which is based on the fast Fourier transformation, and provides an ideal preconditioner for Krylov subspace method. Numerical experiments demonstrate the capability of the proposed fast high order iterative solver.     

Multiscale finite-volume method for parabolic problems arising from compressible multiphase flow in porous media

The multiscale finite-volume (MSFV) method was originally developed for the solution of heterogeneous elliptic problems with reduced computational cost. Recently, some extensions of this method for parabolic problems have been proposed. These extensions proved effective for many cases, however, they are neither general nor completely satisfactory. For instance, they are not suitable for correctly capturing the transient behavior described by the parabolic pressure equation. In this paper, we present a general multiscale finite-volume method for parabolic problems arising, for example, from compressible multiphase flow in porous media. Opposed to previous methods, here, the basis and correction functions are solutions of full parabolic governing equations in localized domains. At the same time, to enhance the computational efficiency of the scheme, the basis functions are kept pressure independent and do not have to be recalculated as pressure evolves. This general approach requires no additional assumptions and its good efficiency and high accuracy is demonstrated for various challenging test cases. Finally, to improve the quality of the results and also to extend the scheme for highly anisotropic heterogeneous problems, it is combined with the iterative MSFV (i-MSFV) method for parabolic problems. As one iterates, the i-MSFV solutions of compressible multiphase problems (parabolic problems) converge to the corresponding fine-scale reference solutions in the same way as demonstrated recently for incompressible cases (elliptic problems). Therefore, the proposed MSFV method can also be regarded as an efficient linear solver for parabolic problems and studies of its efficiency are presented for many test cases.     

An implicit parallel multigrid computing scheme to solve coupled thermal-solute phase-field equations for dendrite evolution

An implicit, second-order space and time discretization scheme together with a parallel multigrid method involving a strip grid domain partitioning has been developed to solve fully coupled, nonlinear phase field equations involving solute and heat transport for multiple solidifying dendrites. The computational algorithm has been shown to be stable and monotonously convergent, and allowed time marching steps that were 3–4 orders of magnitude larger than those employed in similar explicit approaches, resulting in an increase of 3–4 orders of magnitude in computing efficiency. Full solute and thermal coupling was achieved for metallic alloys with a realistic, high Lewis number of >10⁴. The parallel multigrid computing scheme is shown to provide a scalable methodology that allowed the efficient use of distributed supercomputing resource to simulate the evolution of tens of complex shaped 2D dendrites in a computational domain containing tens or even hundreds of millions of grid points. The simulations have provided insight into the dynamic interplay of many growing dendrites in a more realistic fully coupled thermal-solute condition, capturing for the first time fine scale features such as dendrite splitting.     

Dark and kink soliton solutions of the generalized ZK–BBM equation by iterative scheme

The generalized ZK–BBM equation is solved using iterative scheme of the Adomian decomposition method (ADM) and variational iteration method (VIM). A dark and a kink soliton solutions of the generalized ZK–BBM equation are obtained under initial conditions. The convergence analysis of the ADM and VIM solution shows that these solutions are convergent. The comparison of the ADM and VIM solutions with the exact solution shows that the solutions of the generalized ZK–BBM equation by the iterative methods are almost exact. The absolute errors show that the accuracy and efficiency of the ADM and VIM depend on the problem and its domain. It is found that the iterative scheme of Adomian decomposition method and variational iteration method are quite efficient for the soliton solution of the generalized ZK–BBM equation.     

An adaptive GRP scheme for compressible fluid flows

This paper presents a second-order accurate adaptive generalized Riemann problem (GRP) scheme for one and two dimensional compressible fluid flows. The current scheme consists of two independent parts: Mesh redistribution and PDE evolution. The first part is an iterative procedure. In each iteration, mesh points are first redistributed, and then a conservative interpolation formula is used to calculate the cell-averages and the slopes of conservative variables on the resulting new mesh. The second part is to evolve the compressible fluid flows on a fixed nonuniform mesh with the Eulerian GRP scheme, which is directly extended to two-dimensional arbitrary quadrilateral meshes. Several numerical examples show that the current adaptive GRP scheme does not only improve the resolution as well as accuracy of numerical solutions with a few mesh points, but also reduces possible errors or oscillations effectively.     

Uniformly stable numerical schemes for the Boltzmann equation preserving the compressible Navier–Stokes asymptotics

In this paper, we develop a numerical method to solve Boltzmann like equations of kinetic theory which is able to capture the compressible Navier–Stokes dynamics at small Knudsen numbers. Our approach is based on the micro/macro decomposition technique, which applies to general collision operators. This decomposition is performed in all the phase space and leads to an equivalent formulation of the Boltzmann (or BGK) equation that couples a kinetic equation with macroscopic ones. This new formulation is then discretized with a semi-implicit time scheme combined with a staggered grid space discretization. Finally, several numerical tests are presented in order to illustrate the efficiency of our approach. Incidentally, we also introduce in this paper a modification of a standard splitting method that allows to preserve the compressible Navier–Stokes asymptotics in the case of the simplified BGK model. Up to our knowledge, this property is not known for general collision operators.     

An adaptive MHD method for global space weather simulations

A 3D parallel adaptive mesh refinement (AMR) scheme is described for solving the partial-differential equations governing ideal magnetohydrodynamic (MHD) flows. This new algorithm adopts a cell-centered upwind finite-volume discretization procedure and uses limited solution reconstruction, approximate Riemann solvers, and explicit multi-stage time stepping to solve the MHD equations in divergence form, providing a combination of high solution accuracy and computational robustness across a large range in the plasma β (β is the ratio of thermal and magnetic pressures). The data structure naturally lends itself to domain decomposition, thereby enabling efficient and scalable implementations on massively parallel supercomputers. Numerical results for MHD simulations of magnetospheric plasma flows are described to demonstrate the validity and capabilities of the approach for space weather applications     

Pores resolving simulation of Darcy flows

A theoretical formulation and corresponding numerical solutions are presented for microscopic fluid flows in porous media with the domain sufficiently large to reproduce integral Darcy scale effects. Pore space geometry and topology influence flow through media, but the difficulty of observing the configurations of real pore spaces limits understanding of their effects. A rigorous direct numerical simulation (DNS) of percolating flows is a formidable task due to intricacies of internal boundaries of the pore space. Representing the grain size distribution by means of repelling body forces in the equations governing fluid motion greatly simplifies computational efforts. An accurate representation of pore-scale geometry requires that within the solid the repelling forces attenuate flow to stagnation in a short time compared to the characteristic time scale of the pore-scale flow. In the computational model this is achieved by adopting an implicit immersed-boundary method with the attenuation time scale smaller than the time step of an explicit fluid model. A series of numerical simulations of the flow through randomly generated media of different porosities show that computational experiments can be equivalent to physical experiments with the added advantage of nearly complete observability. Besides obtaining macroscopic measures of permeability and tortuosity, numerical experiments can shed light on the effect of the pore space structure on bulk properties of Darcy scale flows.     

A consistent-splitting approach to computing stiff steady-state reacting flows with adaptive chemistry

Splitting techniques have been used extensively for computing reacting flows with detailed chemistry. Nevertheless, there are still some open questions with respect to efficiency and the error introduced by splitting. In this paper, the accuracy and effectiveness of split-operator methods for computing steady-state reacting flows are determined. A fully coupled scheme is described together with two splitting schemes: a standard Strang-splitting scheme and a consistent-splitting scheme, all with implicit transport computations. The effect of splitting errors on the convergence and solution accuracy is investigated analytically using a one-dimensional scalar equation. The accuracy with respect to the original discretized equations is tested for an H₂/O₂ burner flame. Finally, consistent splitting is combined with an adaptive chemistry approach to compute three partially premixed laminar methane flames using detailed chemistry (217 reactions). The calculations confirm that HCO radical concentration is an excellent surrogate for heat release rate.     

A coupled pressure-based computational method for incompressible/compressible flows

Pressure-based flow solvers couple continuity and linearized truncated momentum equations to derive a Poisson type pressure correction equation and use the well known SIMPLE algorithm. Momentum equations and the pressure correction equation are typically solved sequentially. In many cases this method results in slow and often difficult convergence. The current paper proposes a novel computational algorithm, solving for pressure and velocity simultaneously within a pressure-correction coupled solution approach using finite volume method on structured and unstructured meshes. The method can be applied to both incompressible and subsonic compressible flows. For subsonic compressible flows, the energy equation is also coupled with flow field and the density of fluid is obtained by equation of state. The procedure eliminates the pressure correction step, the most expensive component of the SIMPLE-like algorithms. The proposed coupled continuity-momentum-energy equation method can be used to simulate steady state or transient flow problems. The method has been tested on several CFD benchmark cases with excellent results showing dramatically improved numerical convergence and significant reduction in computational time.     

孔隙度数值反演计算

由渗流微分方程定解问题和Peaceman方程给出了网格压力、井底压力对网格孔隙度的导数,利用三维渗流方程压强数值解计算井底压力对网格孔隙度的导数;采用共轭梯度法实现孔隙度均匀(或分块均匀)分布油藏模型的反演计算.算例表明,经过8~10次迭代后反演结果与真值的最大相对误差在0.03%以内,反演收敛于真值.     

Study on the unified algorithm for three-dimensional complex problems covering various flow regimes using Boltzmann model equation

The Boltzmann simplified velocity distribution function equation describing the gas transfer phenomena from various flow regimes will be explored and solved numerically in this study. The discrete velocity ordinate method of the gas kinetic theory is studied and applied to simulate the complex multi-scale flows. Based on the uncoupling technique on molecular movement and colliding in the DSMC method, the gas-kinetic finite difference scheme is constructed to directly solve the discrete velocity distribution functions by extending and applying the unsteady time-splitting method from computational fluid dynamics. The Gauss-type discrete velocity numerical quadrature technique for different Mach number flows is developed to evaluate the macroscopic flow parameters in the physical space. As a result, the gas-kinetic numerical algorithm is established to study the three-dimensional complex flows from rarefied transition to continuum regimes. The parallel strategy adapted to the gas-kinetic numerical algorithm is investigated by analyzing the inner parallel degree of the algorithm, and then the HPF parallel processing program is developed. To test the reliability of the present gas-kinetic numerical method, the three-dimensional complex flows around sphere and spacecraft shape with various Knudsen numbers are simulated by HPF parallel computing. The computational results are found in high resolution of the flow fields and good agreement with the theoretical and experimental data. The computing practice has confirmed that the present gas-kinetic algorithm probably provides a promising approach to resolve the hypersonic aerothermodynamic problems with the complete spectrum of flow regimes from the gas-kinetic point of view of solving the Boltzmann model equation. Supported by the National Natural Science Foundation of China (Grant Nos. 90205009 and 10321002) and the National Parallel Computing Center