An efficient parallel algorithm with application to computational fluid dynamics

When solving time-dependent partial differential equations on parallel computers using the nonoverlapping domain decomposition method, one often needs numerical boundary conditions on the boundaries between subdomains. These numerical boundary conditions can significantly affect the stability and accuracy of the final algorithm.In this paper, a stability and accuracy analysis of the existing methods for generating numerical boundary conditions will be presented, and a new approach based on explicit predictors and implicit correctors will be used to solve convection-diffusion equations on parallel computers, with application to aerospace engineering for the solution of Euler equations in computational fluid dynamics simulations. Both theoretical analyses and numerical results demonstrate significant improvement in stability and accuracy by using the new approach.     

求解任意波数的三维Helmholtz方程

提出通过Adomian分解法求解任意波数的三维Helmholtz方程。通过Adomian分解法可以把三维Helmholtz微分方程转换成递归代数公式,并进一步把其边界条件转换成适用符号计算的简单代数公式。利用边界条件可以很容易得到方程的解析解表达式。Adomian分解法的主要特点在于计算简单快速,并且不需要进行线性化或离散化。最后通过数值计算以验证Adomian分解法求解任意波数下三维Helmholtz方程的有效性。数值计算结果表明：Adomian分解法的计算结果非常接近精确解,并且该方法在大波数情况下还具有良好的收敛性。     

全波场地震波数值模拟中的边界处理实践研究

在用有限差分法或有限元法模拟无界区域中的波动时,需要对计算区域的边界做特殊处理,以消除由于把地震波的传播设定在有限区域而产生的边界反射。为了这一目的,人们研究出了多种人工边界处理方法,完全匹配层（PML）吸收边界条件就是理想的方法之一,现已被广泛应用。本文将PML吸收边界条件应用于全波场地震波的数值模拟,数值计算实验表明,对qP波,匹配层的厚度为5个网格间距即可达到要求,而对qSV波与qSH波,为达到理想的吸收效果,匹配层的厚度应当增大,当厚度为13个网格间距时达到了理想的吸收效果。     

Development of PML absorbing boundary conditions for computational aeroacoustics: A progress review

Recent advances in the development of perfectly matched layer (PML) as absorbing boundary conditions for computational aeroacoustics are reviewed. The PML methodology is presented as a complex change of variables. In this context, the importance of a proper space-time transformation in the PML technique for Euler equations is emphasized. A unified approach for the derivation of PML equations is offered that involves three essential steps. The three-step approach is illustrated in details for the PML of linear and non-linear Euler equations. Numerical examples are also given that include non-reflecting boundary conditions for a ducted channel flow and mixing layer roll-up vortices.     

Numerical Simulation of Pulse Detonation Engine Phenomena

This computational study examines transient, reactive compressible flow phenomena associated with the pulse detonation wave engine. The PDWE is an intermittent combustion engine that relies on unsteady detonation wave propagation for combustion and compression elements of the propulsive cycle. The present computations focus on high order numerical simulations of the generic PDWE configuration with simplified reaction kinetics, so that rapid, straightforward estimates of engine performance may be made. Both one- and two-dimensional simulations of the high speed reactive flow phenomena are performed and compared to determine the applicability of 1D simulations for performance characterization. Examination of the effects of the combustion reaction mechanism and the use of a pressure relaxation length for 1D simulations is made. Characteristic engine performance parameters, in addition to engine noise estimates within and external to the detonation tube, are presented.     

Perfectly matched layers for time-harmonic elastodynamics of unbounded domains: theory and finite-element implementation

One approach to the numerical solution of a wave equation on an unbounded domain uses a bounded domain surrounded by an absorbing boundary or layer that absorbs waves propagating outwards from the bounded domain. A perfectly matched layer (PML) is an unphysical absorbing layer model for linear wave equations that absorbs, almost perfectly, outgoing waves of all non-tangential angles-of-incidence and of all non-zero frequencies. This paper develops the PML concept for time-harmonic elastodynamics in Cartesian coordinates, utilising insights obtained with electromagnetics PMLs, and presents a novel displacement-based, symmetric finite-element implementation of the PML for time-harmonic plane-strain or three-dimensional motion. The PML concept is illustrated through the example of a one-dimensional rod on elastic foundation and through the anti-plane motion of a two-dimensional continuum. The concept is explored in detail through analytical and numerical results from a PML model of the semi-infinite rod on elastic foundation, and through numerical results for the anti-plane motion of a semi-infinite layer on a rigid base. Numerical results are presented for the classical soil–structure interaction problems of a rigid strip-footing on a (i) half-plane, (ii) layer on a half-plane, and (iii) layer on a rigid base. The analytical and numerical results obtained for these canonical problems demonstrate the high accuracy achievable by PML models even with small bounded domains.     

Outflow Boundary Conditions for the Fourier Transformed One-Dimensional Vlasov–Poisson System

In order to facilitate numerical simulations of plasma phenomena where kinetic processes are important, we have studied the technique of Fourier transforming the Vlasov equation analytically in the velocity space, and solving the resulting equation numerically. Special attention has been paid to the boundary conditions of the Fourier transformed system. By using outgoing wave boundary conditions in the Fourier transformed space, small-scale information in velocity space is carried outside the computational domain and is lost. Thereby the so-called recurrence phenomenon is reduced. This method is an alternative to using numerical dissipation or smoothing operators in velocity space. Different high-order methods are used for computing derivatives as well as for the time-stepping, leading to an over-all fourth-order method.     

Curvature and entropy based wall boundary condition for the high order spectral volume Euler solver

A curvature and entropy based wall boundary condition is implemented in the high order spectral volume (SV) context. This method borrows ideas from the “curvature-corrected symmetry technique” developed by (Dadone A, Grossman B. Surface Boundary Conditions for Compressible Flows. AIAA J 1994; 32(2): 285–93), for a low order structured grid Euler solver. After numerically obtaining the curvature, the right state (by convention, the left state is inside the computational domain and the right state lies outside of the computational domain) face pressure values are obtained by solving a linearised system of equations. This is unlike that of the lower order finite volume and difference simulations, wherein the right state face values are trivial to obtain. The right state face density values are then obtained by enforcing entropy conservation. Accuracy studies show that simulations performed by employing the new boundary conditions deliver much more accurate results than the ones which employ traditional boundary conditions, while at the same time asymptotically reaching the desired order of accuracy. Numerical results for two-dimensional inviscid flows around the NACA0012 airfoil and over a bump with the new boundary condition showed dramatic improvements over those with the conventional approach. In all cases and orders, spurious entropy productions with the new boundary treatment are significantly reduced. In general, the numerical results are very promising and indicate that the approach has a great potential for 3D high order simulations.     

Mapped Hybrid Central-WENO Finite Difference Scheme for Detonation Waves Simulations

In this study, we employ the fifth order hybrid Central-WENO conservative finite difference scheme (Hybrid) in the simulation of detonation waves. The Hybrid scheme is used to keep the solutions parts displaying high gradients and discontinuities always captured by the WENO-Z scheme in an essentially non-oscillatory manner while the smooth parts are highly resolved by an efficient and accurate central finite difference scheme and to speedup the computation of the overall scheme. To detect the smooth and discontinuous parts of the solutions, a high order multi-resolution algorithm by Harten is used. A tangent domain mapping is used to cluster grid points near the detonation front in order to enhance the grid resolution within half reaction zone that drives the development of complex nonlinear wave structures behind the front. We conduct several numerical comparisons among the WENO-Z scheme with a uniformly spaced grid, the WENO-Z scheme and the Hybrid scheme with the domain mapping in simulations of classical stable and unstable detonation waves. One- and two-dimensional numerical examples show that the increased grid resolution in the half reaction zone by the Mapped WENO-Z scheme and the Mapped Hybrid scheme allows a significant increased efficiency and accuracy when compares with the solution obtained with a highly resolved one computed by the WENO-Z scheme with a uniformly spaced grid. Results of three-dimensional simulations of stable, slightly unstable and highly unstable detonation waves computed by the Mapped Hybrid scheme are also presented.     

Boussinesq Systems of Bona-Smith Type on Plane Domains: Theory and Numerical Analysis

We consider a class of Boussinesq systems of Bona-Smith type in two space dimensions approximating surface wave flows modelled by the three-dimensional Euler equations. We show that various initial-boundary-value problems for these systems, posed on a bounded plane domain are well posed locally in time. In the case of reflective boundary conditions, the systems are discretized by a modified Galerkin method which is proved to converge in L ² at an optimal rate. Numerical experiments are presented with the aim of simulating two-dimensional surface waves in realistic (plane) domains with a variety of initial and boundary conditions, and comparing numerical solutions of Bona-Smith systems with analogous solutions of the BBM-BBM system.     

Propagation and transformation of periodic nonlinear shallow-water waves in basins with selected breakwater systems

This paper describes the development of a Boussinesq three-equation model for simulating propagation and transformation of periodic nonlinear waves (cnoidal waves) in an arbitrary shallow-water basin. The Boussinesq equations in terms of depth-averaged horizontal velocities and free-surface elevation are solved numerically in a curvilinear coordinate system. An Euler’s predictor-corrector finite-difference algorithm is applied for numerical computation. The effects of irregular boundary, non-uniform water depth and coastal structures inside a basin are all included in the model simulation. A second-order cnoidal wave solution for the Boussinesq equations is used as an incident wave condition. A set of open boundary conditions is also applied to effectively transmit waves out of the computational domain. Model tests were conducted by simulating waves propagating past an isolated breakwater. The effect of variable depth was examined with modeling waves over an uneven bottom with convex ramp topography. The overall evolution of wave propagation, diffraction and reflection in coupled harbors with various layouts of inner and outer breakwaters was also studied. Data comparisons reveal that the simulated wave heights agree reasonably well with laboratory measurements, especially in the region of inner basin.     

Block iterative methods for the numerical solution of three dimensional mildly non-linear biharmonic problems of first kind

In this article, we discuss two sets of new finite difference methods of order two and four using 19 and 27 grid points, respectively over a cubic domain for solving the three dimensional nonlinear elliptic biharmonic problems of first kind. For both the cases we use block iterative methods and a single computational cell. The numerical solution of (?u/?n) are obtained as by-product of the methods and we do not require fictitious points in order to approximate the boundary conditions. The resulting matrix system is solved by the block iterative method using a tri-diagonal solver. In numerical experiments the proposed methods are compared with the exact solutions both in singular and non-singular cases.     

波导传感器中高效柱面波平面波转换结构的变换光学设计

基于阻抗可调的变换光学理论,设计了波导传感器中二维平面波和柱面波相互转换的光学结构。通过在虚拟空间设置合适的阻抗函数,变换介质边界上的反射得到有效消除,实现了高效的平面波到柱面波转换;转换效率将随着输入波束宽度的增加而提高。四个相似转换结构的组合,可实现高效柱面波-平面波转换,与采用传统的变换光学设计比较,基于阻抗可调变换光学理论的设计,数值模拟表明转换效率可提高约百分之六。     

Application of shifted Jacobi pseudospectral method for solving (in)finite-horizon min–max optimal control problems with uncertainty

The difficulty of solving the min–max optimal control problems (M-MOCPs) with uncertainty using generalised Euler–Lagrange equations is caused by the combination of split boundary conditions, nonlinear differential equations and the manner in which the final time is treated. In this investigation, the shifted Jacobi pseudospectral method (SJPM) as a numerical technique for solving two-point boundary value problems (TPBVPs) in M-MOCPs for several boundary states is proposed. At first, a novel framework of approximate solutions which satisfied the split boundary conditions automatically for various boundary states is presented. Then, by applying the generalised Euler–Lagrange equations and expanding the required approximate solutions as elements of shifted Jacobi polynomials, finding a solution of TPBVPs in nonlinear M-MOCPs with uncertainty is reduced to the solution of a system of algebraic equations. Moreover, the Jacobi polynomials are particularly useful for boundary value problems in unbounded domain, which allow us to solve infinite- as well as finite and free final time problems by domain truncation method. Some numerical examples are given to demonstrate the accuracy and efficiency of the proposed method. A comparative study between the proposed method and other existing methods shows that the SJPM is simple and accurate.     

A Pseudospectral Penalty Scheme for 2D Isotropic Elastic Wave Computations

In this paper, we present a pseudospectral scheme for solving 2D elastic wave equations. We start by analyzing boundary operators leading to the well-posedness of the problem. In addition, equivalent characteristic boundary conditions of common physical boundary conditions are discussed. These theoretical results are further employed to construct a Legendre pseudospectral penalty scheme based on a tensor product formulation for approximating waves on a general curvilinear quadrilateral domain. A stability analysis of the scheme is conducted for the case where a straight-sided quadrilateral element is used. The analysis shows that, by properly setting the penalty parameters, the scheme is stable at the semi-discrete level. Numerical experiments for testing the performance of the scheme are conducted, and the expected p- and h-convergence patterns are observed. Moreover, the numerical computations also show that the scheme is time stable, which makes the scheme suitable for long time simulations. This work is supported by National Science Council grant No. NSC 95-2120-M-001-003.     

带矢量ABC底面的共形完全匹配层

共形完全匹配层是一种有耗各向异性媒质组成的凸且光滑的壳体,其底面一般是PEC面或PMC面,但是PEC面或PMC面会对原散射场产生反射;为了减少底面反射,将CPML原有的PEC（或PMC）底面改为矢量ABC吸收边界,并给出了带矢量ABC底面的CPML泛函公式。通过数值算例证明,这种带矢量ABC底面的CPML边界不仅减少了底面反射,而且吸收效果好,计算精度高。     

基于CFD的船舶斜浪三自由度运动仿真研究

针对船舶在斜浪中的多自由度耦合运动问题, 建立了三维粘性流耐波性数值波浪水池, 采用边界条件造波法生成高精度的斜向规则波, 通过同时求解RANS方程和刚体运动学方程, 结合网格整体移动方法和滑移网格方法, 实现了船舶斜浪航行的垂荡、纵摇及横摇三自由度耦合运动数值模拟。给出了DTMB5512船模斜浪中的垂荡、纵摇及横摇的频率响应函数, 与线性切片理论计算结果进行比较, 吻合良好。该方法可为船舶斜浪航行的耐波性预报提供一种有效的途径。     

Estimation of wave responses from subvertical macrofracture systems using a grid characteristic method

This work investigates the formation and propagation of scattered waves forming a response of crack structures on the seismologic record. The initial impulse is a flat wave front propagating into the depth of the medium. The paper considers the periodic structure of the scattered wave response from a system (cluster) of subvertical cracks. The procedures of estimation of geometrical parameters of such crack structures are based on numerical experiments. We use a grid-characteristic method to calculate the triangular grids with the formulation of boundary conditions at the interface of the medium and the cracks, as well as on the borders of the region of integration with the characteristic properties of the governing equations of the hyperbolic type. This numerical method makes it possible to build the most correct numerical algorithms on the boundaries of the integration domain and on the surfaces of the media division (interfaces) and to take into account the domain of the solution’s dependence and the physics of the problem (propagation of disturbances on the characteristic directions). This is why this method is perhaps the most suitable one for the numerical solution of dynamic problems with a distinct wave character in geologically substantially inhomogeneous media, in particular, for the considered problem of seismic waves interaction with fissured structures.     

双螺杆挤出机流场数值模拟中流道进出口边界条件的探讨

在对双螺杆挤出机流场的数值模拟中,流道进出口边界条件的设置一直是一个颇具争议的问题。由于事先无法获得计算域进出口平面上的真实边界条件,研究人员在进行双螺杆挤出机的流场分析时,大都采用放松边界条件。为了考察放松边界条件对双螺杆挤出机流场数值模拟结果的影响,本文采用聚合物流动分析软件POLYFLOW,在流量恒定的前提下对双螺杆挤出机流道进出口给定三种不同分布形式的速度边界条件,对其流场进行了数值模拟。数值计算结果表明,在体积流量恒定的条件下,流道进出口不同分布形式的速度边界条件对流场的影响主要集中在进出口附近区域,但对离进出口边界较远的流场影响很小。一般而言,当计算域所对应的螺杆较长时,可以忽略流道进出口的放松边界条件所引起的误差;当计算域较短时,不宜直接采用放松边界条件,而应根据螺杆的实际构型．在计算域的进出口增加适当长度的发展段。     

High Order Accurate Finite Difference Modeling of Seismo-Acoustic Wave Propagation in a Moving Atmosphere and a Heterogeneous Earth Model Coupled Across a Realistic Topography

We develop a numerical method for simultaneously simulating acoustic waves in a realistic moving atmosphere and seismic waves in a heterogeneous earth model, where the motions are coupled across a realistic topography. We model acoustic wave propagation by solving the linearized Euler equations of compressible fluid mechanics. The seismic waves are modeled by the elastic wave equation in a heterogeneous anisotropic material. The motion is coupled by imposing continuity of normal velocity and normal stresses across the topographic interface. Realistic topography is resolved on a curvilinear grid that follows the interface. The governing equations are discretized using high order accurate finite difference methods that satisfy the principle of summation by parts. We apply the energy method to derive the discrete interface conditions and to show that the coupled discretization is stable. The implementation is verified by numerical experiments, and we demonstrate a simulation of coupled wave propagation in a windy atmosphere and a realistic earth model with non-planar topography.