天气预报有限差分格点模式并行计算研究

有限差分格点模式是目前广泛采用的中尺度数值天气预报模式之一。本文分析了有限差分格点模式典型代表MM5并行计算的特点,对其多重网格嵌套、数据划分、通信、优化等做了深入研究,最后给出了MM5模式在银河巨型计算机上的测试结果。

不同亚格子模型在亚声速槽道流大涡模拟中的应用对比

湍流边界层流动是一种广泛存在于飞行器内部和外部的流动现象,是基础理论和模型验证的重要研究对象.能够捕捉大部分流动细节且计算量适中的大涡模拟(large-eddy simulation, LES)方法在湍流数值模拟中得到了越来越广泛的应用.文章基于格心有限差分方法,使用4阶紧致中心格式离散N-S方程无黏项,分别应用5种不同的亚格子(subgrid-scale, SGS)模型,即隐式, SM(Smagorinsky model), DSM(dynamic Smagorinsky model), WALE(wall-adapting local eddy-viscosity model)和CSM(coherent structures model),对Re=3 000, Ma=0.5的等温壁面槽道流动进行了大涡模拟研究.与实验值和直接数值模拟(direct numerical simulation, DNS)结果对比后发现,流场平均温度、平均密度等热力学量以及平均流向速度对亚格子模型不敏感,不适宜作为判断模型优劣的判据.亚格子模型在壁面附近的耗散越大,壁面摩擦速度以及阻力系数就越小.对于与速度相关的脉动量来说,不同模型得到的结果在壁面和脉动峰值附近误差比较大,中心线附近较小;显式模型结果在流向速度峰值处均高于参考值,而在展向和壁面法向速度脉动峰值处则均偏低.考虑显式的4种模型在壁面附近的涡黏系数分布, DSM和CSM曲线满足涡黏系数与无量纲壁面距离3次方成正比的分布规律, SM曲线斜率偏小而WALE曲线斜率偏大.

PML在三维柱坐标系下的应用及近场-远场转换

在三维柱坐标系下推导了麦克斯韦方程组在完全匹配层(PML)吸收边界条件中的时域有限差分(FDTD)格式并编程实现.然后用一种简便的线性插值法给出了在柱坐标系下由近场到远场转换的FDTD计算方法.最后用这种方法对脉冲天线的远区辐射场进行了数值模拟.

CONSERVATIVE DIFFERENCE SCHEME BASED ON NUMERICAL ANALYSIS FOR NONLINEAR SCHRDINGER EQUATION WITH WAVE OPERATOR

A new conservative finite difference scheme is presented based on the numerical analysis for an initialboundary value problem of a class of Schroedinger equation with the wave operator. The scheme can be linear and implicit or explicit based on the parameter choice. The initial value after discretization has second-order accuracy that is consistent with the scheme accuracy. The existence and the uniqueness of the difference solution are proved. Based on the priori estimates and an inequality about norms, the stability and the convergence of difference solutions with the second-order are proved in the energy norm. Experimental results demonstrate the efficiency of the new scheme.

A-high-order Accuraqcy Implicit Difference Scheme for Solving the Equation of Parabolic Type

In this paper,a implicit difference scheme is proposed for solving the equation of one_dimension parabolic type by undetermined paameters.The stability condition is r=αΔt/Δx 2 1/2 and the truncation error is o(Δt 4+Δx 4) It can be easily solved by double sweeping method.

An Explicit Difference Scheme with High Accuracy and Branching Stability for Solving Parabolic Partial Differential Equation

This paper presents an explicit difference scheme with accuracy and branching stability for solving onedimensional parabolic type equation by the method of undetermined parameters and its truncation error is O(△t4+△x4). The stability condition is r=a△t/△x2<1/2.

A Family of High_order Accuracy Explicit Difference Schemes for Solving 2-D Parabolic Partial Differential Equation

A family of high_order accuracy explicit difference schemes for solving 2_dimension parabolic P.D.E. are constructed. Th e stability condition is r=Δt/Δx 2=Δt/Δy 2【1/2 and the truncation err or is O(Δt 3+Δx 4).

LARGE-EDDY SIMULATION OF TWO-PHASE REACTING FLOW IN MODEL COMBUSTOR

The gas-droplet two-phase reacting flow in a model combustor with the V-gutter flame holder is studied by an Eulerian-Lagrangian large-eddy simulation (LES) approach. The k-equation subgrid-scale model is used to simulate the subgrid eddy viscosity, and the eddy-break-up (EBU) combustion subgrid-scale model is used to determine the chemical reaction rate. A two-step turbulent combustion subgrid-scale model is employed for calculating carbon monoxide CO concentration, and the NO subgrid-scale pollutant formation model for the evaluation of the rate of NO formation. The heat flux model is applied to the prediction of radiant heat transfer. The gas phase is solved with the SIMPLE algorithm and a hybrid scheme in the staggered grid system. The liquid phase equations are solved in a Lagrangian frame in reference of the particle-source-in-cell (PSIC) algorithm. From simulation results, the exchange of mass, moment and energy between gas and particle fields for the reacting flow in the afterburner with a V-gutter flame holder can be obtained. By the comparison of experimental and simulation results, profile temperature and pollutant of the outlet are quite in agreement with experimental data. Results show that the LES approach for predicting the two-phase instantaneous reacting flow and pollutant emissions in the afterburner is feasible.

Multicomponent seismic forward modeling of gas hydrates beneath the seafloor

We investigated the effect of microscopic distribution modes of hydrates in porous sediments, and the saturation of hydrates and free gas on the elastic properties of saturated sediments. We simulated the propagation of seismic waves in gas hydrate-bearing sediments beneath the seafloor, and obtained the common receiver gathers of compressional waves（P-waves） and shear waves（S-waves）. The numerical results suggest that the interface between sediments containing gas hydrates and free gas produces a large-amplitude bottomsimulating reflector. The analysis of multicomponent common receiver data suggests that ocean-bottom seismometers receive the converted waves of upgoing P- and S-waves, which increases the complexity of the wavefield record.

HIGH RESOLUTION POSITIVITY-PRESERVING DIFFERENCE SCHEMES FOR TWO DIMENSIONAL EULER EQUATIONS

A class of high resolution positivity preserving Boltzmann type difference schemes for one and two dimensional Euler equations is studied. First, the relation between Boltzmann and Euler equations is analyzed. By using a kind of special interpolation, the high resolution Boltzmann type difference scheme is constructed. Finally, numerical tests show that the schemes are effective and useful.

NUMERICAL SIMULATIONS OF EVOLUTIONS FROM WALL PULSE TO TURBULENT COHE RENT STRUCTURES

The high order compact d if ference method is developed for solving the perturbation equations based on Navi er Stokes equations, and is used in studying complex evolution processes from w all negative pulse to the turbulent coherent structure in the channel flow. Th is method contains three dimensional coupling difference scheme with high accur acy and high resolution, and the high order time splitting methods. Compared with the general spectral method, the method can be used to research turbule nt coherent structure under more general boundary conditions and in flow domains . In this paper, the generation and evolution of the turbulent coherent structur es ind uced by wall pulse in the channel flow are simulated, and the basic characterist ics and rules of the turbulent coherent structure are shown. Computational r esults indicate that a wall negative pulse is more convenient than the resonant three wave model.

Stability of finite difference numerical simulations of acoustic logging-while-drilling with different perfectly matched layer schemes

In acoustic logging-while-drilling （ALWD） finite difference in time domain （FDTD） simulations, large drill collar occupies, most of the fluid-filled borehole and divides the borehole fluid into two thin fluid columns （radius -27 mm）. Fine grids and large computational models are required to model the thin fluid region between the tool and the formation. As a result, small time step and more iterations are needed, which increases the cumulative numerical error. Furthermore, due to high impedance contrast between the drill collar and fluid in the borehole （the difference is 〉30 times）, the stability and efficiency of the perfectly matched layer （PML） scheme is critical to simulate complicated wave modes accurately. In this paper, we compared four different PML implementations in a staggered grid finite difference in time domain （FDTD） in the ALWD simulation, including field-splitting PML （SPML）, multiaxial PML（M- PML）, non-splitting PML （NPML）, and complex frequency-shifted PML （CFS-PML）. The comparison indicated that NPML and CFS-PML can absorb the guided wave reflection from the computational boundaries more efficiently than SPML and M-PML. For large simulation time, SPML, M-PML, and NPML are numerically unstable. However, the stability of M-PML can be improved further to some extent. Based on the analysis, we proposed that the CFS-PML method is used in FDTD to eliminate the numerical instability and to improve the efficiency of absorption in the PML layers for LWD modeling. The optimal values of CFS-PML parameters in the LWD simulation were investigated based on thousands of 3D simulations. For typical LWD cases, the best maximum value of the quadratic damping profile was obtained using one do. The optimal parameter space for the maximum value of the linear frequency-shifted factor （a0） and the scaling factor （β0） depended on the thickness of the PML layer. For typical formations, if the PML thickness is 10 grid points, the global error can be reduced to 〈1% using the optimal PML parameters, and the error will decrease as the PML thickness increases.

True-amplitude wavefield separation using staggered-grid interpolation in the wavenumber domain

Wavefield separation of multicomponent seismic data to image subsurface structures can be realized in either the space domain or the wavenumber domain. However, as the particle velocity components used in the wavenumber-domain wavefield separation are not defined at the same grid point with the staggered-grid finite-difference method for elastic wavefield simulation, we propose the wavenumber-domain interpolation method to estimate the required values at the common grid points prior to the wavenumber-domain true-amplitude wavefield separation. Moreover, numerical experiments show that the wavenumber-domain interpolation method has high interpolation accuracy and the trueamplitude wavefield separation method shows good amplitude preservation. The application of the proposed methodology to elastic reverse-time migration can obtain good amplitudepreserved images even in the case of some velocity error.

The Stability Research for the Finite Difference Scheme of a Nonlinear Reaction-diffusion Equation

In the article, the fully discrete finite difference scheme for a type of nonlinear reaction-diffusion equation is established. Then the new function space is introduced and the stability problem for the finite difference scheme is discussed by means of variational approximation method in this function space. The approach used is of a simple characteristic in gaining the stability condition of the scheme.

Second-order difference scheme for a nonlinear model of wood drying process

A numerical simulation for a model of wood drying process is considered. The model is given by a couple of nonlinear differential equations. One is a nonlinear parabolic equation and the other one is a nonlinear ordinary equation. A difference scheme is derived by the method of reduction of order. First, a new variable is introduced and the original problem is rewritten into a system of the first-order differential equations. Secondly, a difference scheme is constructed for the later problem. The solvability, stability and convergence of the difference scheme are proved by the energy method. The convergence order of the difference scheme is secondorder both in time and in space. A prior error estimate is put forward. The new variable is put aside to reduce the computational cost. A numerical example testifies the theoretical result.

各向同性湍流通过正激波的演化特征研究

激波与湍流相互作用(shock-turbulence interaction,STI)是空气动力学研究中的一个基础问题.基于格心有限差分法(cell-centered finite difference method,CCFDM)求解器Helios,采用五阶加权紧致非线性格式(weighted compact nonlinear scheme,WCNS)对各向同性湍流通过正激波的情形进行直接数值模拟(direct numerical simulation,DNS).对湍流相关物理量进行统计,分析结果表明,在湍流中波后的密度、温度和压力较无湍流情形下略小,而速度则略大,均在波后呈现短暂过冲然后缓慢向理论值逼近的变化趋势;波后流向雷诺应力突降随之快速增长又衰减,呈现非单调变化趋势,线性相互作用分析(linear interaction analysis,LIA)将其归结为波后能量从声模式转移为涡模式方式,与流向不同,横向雷诺应力突增后单调衰减,波后雷诺应力各向异性明显且随下游距离逐渐增强;波后湍动能突增后呈现非单调变化趋势;泰勒微尺度和Kolmogorov尺度过激波后均明显减小,说明波后湍流长度尺度变小,从而对波后网格的分辨率提出了更高的要求;密度、温度和压力过激波后脉动均方根均增加,密度和压力脉动强度减小,温度脉动强度增大.

Finite difference method for dynamic response analysis of anchorage system

Based on some assumptions, the dynamic analysis model of anchorage system is established. The dynamic governing equation is expressed as finite difference format and programmed by using MATLAB language. Compared with theoretical method, the finite difference method has been verified to be feasible by a case study. It is found that under seismic loading, the dynamic response of anchorage system is synchronously fluctuated with the seismic vibration. The change of displacement amplitude of material points is slight, and comparatively speaking, the displacement amplitude of the outside point is a little larger than that of the inside point, which shows amplification effect of surface. While the axial force amplitude transforms considerably from the inside to the outside. It increases first and reaches the peak value in the intersection between the anchoring section and free section, then decreases slowly in the free section. When considering damping effect of anchorage system, the finite difference method can reflect the time attenuation characteristic better, and the calculating result would be safer and more reasonable than the dynamic steady-state theoretical method. What is more, the finite difference method can be applied to the dynamic response analysis of harmonic and seismic random vibration for all kinds of anchor, and hence has a broad application prospect.

Improving the Stability Problem of the Finite Difference Scheme for Reaction-diffusion Equation

This paper deals with the special nonlinear reaction-diffusion equation. The finite difference scheme with incremental unknowns approximating to the differential equation （2.1） is set up by means of introducing incremental unknowns methods. Through the stability analyzing for the scheme, it was shown that the stability conditions of the finite difference schemes with the incremental unknowns are greatly improved when compared with the stability conditions of the corresponding classic difference scheme.

Highly Efficient Lattice Boltzmann Model for Compressible Fluids:Two-Dimensional Case

We present a highly efficient lattice Boltzmann model for simulating compressible flows. This model is based on the combination of an appropriate finite difference scheme, a 16-discrete-velocity model [Kataoka and Tsutahara, Phys. Rev. E 69 （2004） 035701（R）] and reasonable dispersion and dissipation terms. The dispersion term effectively reduces the oscillation at the discontinuity and enhances numerical precision. The dissipation term makes the new model more easily meet with the yon Neumann stability condition. This model works for both high-speed and low-speed flows with arbitrary specific-heat-ratio. With the new model simulation results for the well-known benchmark problems get a high accuracy compared with the analytic or experimental ones. The used benchmark tests include （i） Shock tubes such as the Sod, Lax, Sjogreen, Colella explosion wave, and collision of two strong shocks, （ii） Regular and Mach shock reflections, and （iii） Shock wave reaction on cylindrical bubble problems. With a more realistic equation of state or free-energy functional, the new model has the potential tostudy the complex procedure of shock wave reaction on porous materials.

Uniform Convergence Analysis of Finite Difference Scheme for Singularly Perturbed Delay Differential Equation on an Adaptively Generated Grid

Adaptive grid methods are established as valuable computational technique in approximating effectively the solutions of problems with boundary or interior layers. In this paper,we present the analysis of an upwind scheme for singularly perturbed differential-difference equation on a grid which is formed by equidistributing arc-length monitor function.It is shown that the discrete solution obtained converges uniformly with respect to the perturbation parameter.Numerical experiments illustrate in practice the result of convergence proved theoretically.