Impact of spatial resolution and spatial difference accuracy on the performance of Arakawa A-D grids

This paper alms at illustrating the impact of spatial difference scheme and spatial resolution on the performance of Arakawa A-D grids in physical space. Linear shallow water equations are discretized and forecasted on Arakawa A-D grids for 120-minute using the ordinary second-order (M and fourth-order (C4) finite difference schemes with the grid spacing being 100 km, 10 km and I km, respectively. Then the forecasted results are compared with the exact solution, the result indicates that when the grid spacing is I kin, the inertial gravity wave can be simulated on any grid with the same results from C2 scheme or C4 scheme, namely the impact of variable configuration is neglectable; while the inertial gravity wave is simulated with lengthened grid spacing, the effects of different variable configurations are different. However, whether for C2 scheme or for C4 scheme, the RMS is minimal (maximal) on C (D) grid. At the same time it is also shown that when the difference accuracy increases from C2 scheme to C4 scheme, the resulted forecasts do not uniformly decrease, which is validated by the change of the group A velocity relative error from C2 scheme to C4 scheme. Therefore, the impact of the grid spacing is more important than that of the difference accuracy on the performance of Arakawa A-D grid.

Uniformly convergent finite element and finite difference methods for singularly perturbed ordinary differential equations

This thesis is concerned with uniformly convergent finite element and finite difference methods for numerically solving singularly perturbed two-point boundary value problems. We examine the following four problems: (i) high order problem of reaction-diffusion type; (ii) high order problem of convection-diffusion type; (iii) second order interior turning point problem; (iv) semilinear reaction-diffusion problem. Firstly, we consider high order problems of reaction-diffusion type and convection-diffusion type. Under suitable hypotheses, the coercivity of the associated bilinear forms is proved and representation results for the solutions of such problems are given. It is shown that, on an equidistant mesh, polynomial schemes cannot achieve a high order of convergence which is uniform in the perturbation parameter. Piecewise polynomial Galerkin finite element methods are then constructed on a Shishkin mesh. High order convergence results, which are uniform in the perturbation parameter, are obtained in various norms. Secondly, we investigate linear second order problems with interior turning points. Piecewise linear Galerkin finite element methods are generated on various piecewise equidistant meshes designed for such problems. These methods are shown to be convergent, uniformly in the singular perturbation parameter, in a weighted energy norm and the usual L2 norm. Finally, we deal with a semilinear reaction-diffusion problem. Asymptotic properties of solutions to this problem are discussed and analysed. Two simple finite difference schemes on Shishkin meshes are applied to the problem. They are proved to be uniformly convergent of second order and fourth order respectively. Existence and uniqueness of a solution to both schemes are investigated. Numerical results for the above methods are presented.

Energy conserving schemes for the simulation of musical instrument contact dynamics

Collisions are an innate part of the function of many musical instruments. Due to the nonlinear nature of contact forces, special care has to be taken in the construction of numerical schemes for simulation and sound synthesis. Finite difference schemes and other time-stepping algorithms used for musical instrument modelling purposes are normally arrived at by discretising a Newtonian description of the system. However because impact forces are non-analytic functions of the phase space variables, algorithm stability can rarely be established this way. This paper presents a systematic approach to deriving energy conserving schemes for frictionless impact modelling. The proposed numerical formulations follow from discretising Hamilton׳s equations of motion, generally leading to an implicit system of nonlinear equations that can be solved with Newton׳s method. The approach is first outlined for point mass collisions and then extended to distributed settings, such as vibrating strings and beams colliding with rigid obstacles. Stability and other relevant properties of the proposed approach are discussed and further demonstrated with simulation examples. The methodology is exemplified through a case study on tanpura string vibration, with the results confirming the main findings of previous studies on the role of the bridge in sound generation with this type of string instrument.

Binaural Reproduction of Finite Difference Simulations Using Spherical Array Processing

Due to its efficiency and simplicity, the finite-difference time-domain method is becoming a popular choice for solving wideband, transient problems in various fields of acoustics. So far, the issue of extracting a binaural response from finite difference simulations has only been discussed in the context of embedding a listener geometry in the grid. In this paper, we propose and study a method for binaural response rendering based on a spatial decomposition of the sound field. The finite difference grid is locally sampled using a volumetric array of receivers, from which a plane wave density function is computed and integrated with free-field head related transfer functions, in the spherical harmonics domain. The volumetric array is studied in terms of numerical robustness and spatial aliasing. Analytic formulas that predict the performance of the array are developed, facilitating spatial resolution analysis and numerical binaural response analysis for a number of finite difference schemes. Particular emphasis is placed on the effects of numerical dispersion on array processing and on the resulting binaural responses. Our method is compared to a binaural simulation based on the image method. Results indicate good spatial and temporal agreement between the two methods.

Stochastic regularization method for backward Cauchy problem in Banach spaces

We consider a stochastic regularization method for solving the backward Cauchy problem in Banach spaces. An order of convergence is obtained on sourcewise representative elements.

A second-order accurate numerical approximation for the Riesz space fractional advection-dispersion equation

In this paper, we consider a space Riesz fractional advection-dispersion equation. The equation is obtained from the standard advection-diffusion equation by replacing the ¯rst-order and second-order space derivatives by the Riesz fractional derivatives of order β 1 Є (0; 1) and β2 Є(1; 2], respectively. Riesz fractional advection and dispersion terms are approximated by using two fractional centered difference schemes, respectively. A new weighted Riesz fractional ¯nite difference approximation scheme is proposed. When the weighting factor Ѳ = 1/2, a second- order accurate numerical approximation scheme for the Riesz fractional advection-dispersion equation is obtained. Stability, consistency and convergence of the numerical approximation scheme are discussed. A numerical example is given to show that the numerical results are in good agreement with our theoretical analysis.

A novel numerical approximation for the space fractional advection-dispersion equation

In this paper, we consider a space fractional advection–dispersion equation. The equation is obtained from the standard advection–diffusion equation by replacing the first- and second-order space derivatives by the Riesz fractional derivatives of order β1 ∈ (0, 1) and β2 ∈ (1, 2], respectively. The fractional advection and dispersion terms are approximated by using two fractional centred difference schemes. A new weighted Riesz fractional finite-difference approximation scheme is proposed. When the weighting factor θ equals 12, a second-order accuracy scheme is obtained. The stability, consistency and convergence of the numerical approximation scheme are discussed. A numerical example is given to show that the numerical results are in good agreement with our theoretical analysis.

On the propagation of a multi-dimensional shock of arbitrary strength

In this paper the kinematics of a curved shock of arbitrary strength has been discussed using the theory of generalised functions. This is the extension of Moslov’s work where he has considered isentropic flow even across the shock. The condition for a nontrivial jump in the flow variables gives the shock manifold equation (sme). An equation for the rate of change of shock strength along the shock rays (defined as the characteristics of the sme) has been obtained. This exact result is then compared with the approximate result of shock dynamics derived by Whitham. The comparison shows that the approximate equations of shock dynamics deviate considerably from the exact equations derived here. In the last section we have derived the conservation form of our shock dynamic equations. These conservation forms would be very useful in numerical computations as it would allow us to derive difference schemes for which it would not be necessary to fit the shock-shock explicitly.

Direct numerical simulation of turbulent spots

A group of high-order finite-difference schemes for incompressible flow was implemented to simulate the evolution of turbulent spots in channel flows. The long-time accuracy of these schemes was tested by comparing the evolution of small disturbances to a plane channel flow against the growth rate predicted by linear theory. When the perturbation is the unstable eigenfunction at a Reynolds number of 7500, the solution grows only if there are a comparatively large number of (equispaced) grid points across the channel. Fifth-order upwind biasing of convection terms is found to be worse than second-order central differencing. But, for a decaying mode at a Reynolds number of 1000, about a fourth of the points suffice to obtain the correct decay rate. We show that this is due to the comparatively high gradients in the unstable eigenfunction near the walls. So, high-wave-number dissipation of the high-order upwind biasing degrades the solution especially. But for a well-resolved calculation, the weak dissipation does not degrade solutions even over the very long times (O(100)) computed in these tests. Some new solutions of spot evolution in Couette flows with pressure gradients are presented. The approach to self-similarity at long times can be seen readily in contour plots.

双曲守恒型方程的二阶摄动有限差分格式

对双曲守恒型方程，将其一阶迎风格式空间差商的常系数摄动展开为时间步长和空间步长的幂级数，通过确定幂级数系数而获得二阶精度的摄动有限差分（PFD）格式。进而从双曲守恒型方程的通量分裂型一阶迎风格式出发，通过娄似的摄动展开方法，获得空间精度为二阶的通量分裂形式的摄动有限差分（FPFD）格式。这两类格式保留了一阶守恒迎风格式的简洁结构形式，使用三节点即可达到二阶精度，又避免了三点二阶格式的非物理数值振荡。并将这两类格式推广应用到双曲守恒型方程组，最后通过模型方程和一维激波管流动的数值算例验证了格式的高精度、高分辨率性质。

Direct Numerical Simulation of Supersonic Turbulent Boundary Layer Flow

Direct numerical simulations of a spatially evolving supersonic flat-plate turbulent boundary layer flow with free Mach number M = 2.25 and Reynolds number Re = 365000/in are performed. The transition process from laminar to turbulent flow is obtained by solving the three-dimensional compressible Navier-Stokes, equations, using high-order accurate difference schemes. The obtained statistical results agree well with the experimental and theoretical data. From the numerical results it can be seen that the transition process under the considered conditions is the process which skips the Tolimien-Schlichting instability and the second instability through the instability of high gradient shear layer and becomes of laminar flow breakdown. This means that the transition process is a bypass-type transition process. The spanwise asymmetry of the disturbance locally upstream imposed is important to induce the bypass-type transition. Furthermore, with increasing the time disturbance frequency the transition will delay. When the time disturbance frequency is large enough, the transition will disappear.

Numerical-simulation of low supersonic-flow around a sphere

The controlled equations defined in a physical plane are changed into those in a computational plane with coordinate transformations suitable for different Mach number M(infinity). The computational area is limited in the body surface and in the vicinities of detached shock wave and sonic line. Thus the area can be greatly cut down when the shock wave moves away from the body surface as M(infinity) --> 1. Highly accurate, total variation diminishing (TVD) finite-difference schemes are used to calculate the low supersonic flowfield around a sphere. The stand-off distance, location of sonic line, etc. are well comparable with experimental data. The long pending problem concerning a flow passing a sphere at 1.3 greater-than-or-equal-to M(infinity) > 1 has been settled, and some new results on M(infinity) = 1.05 have been presented.

粘性可压混合层时间稳定性对称紧致差分求解

基于可压扰动方程组的一阶改型，将高精度对称紧致格式引入边值法数值线性稳定分析。对所获非线性离散特征值问题给出了一个通用形式二阶迭代局部算法，实现了时间模式和空间模式的统一求解，并将扰动特征及其特征函数同时得到。据此分析了可压平面自由混合层时间稳定性，涉及二维/三维扰动波、粘性/无粘扰动波、第一/第二模态、特征函数、伪特征值谱等。研究表明，压缩性效应和粘性效应对最不稳定扰动波数和增长率呈相似的减抑作用；在Mc = 1附近，从高波数段开始，粘性效应可强化二维不稳定扰动波由第一模态向第二模态的过渡。

槽道湍流的直接数值模拟

通过直接数值模拟(DNS)研究槽道湍流的性质和机理。包含五个部分：1）湍流直接数值模拟的差分方法研究。2）求解不可压N－S方程的高效算法和不可压槽道湍流的直接数值模拟。3）可压缩槽道湍流的直接数值模拟和压缩性机理分析。4）“二维湍流”的机理分析。5）槽道湍流的标度律分析。1．针对壁湍流计算网格变化剧烈的特点，构造了基于非等距网格的的迎风紧致格式。该方法直接针对计算网格构造格式中的系数，克服了传统方法采用 Jacobian 变换因网格变化剧烈而带来的误差。针对湍流场的多尺度特性分析了差分格式的精度、网格尺度与数值模拟能分辨的最小尺度的关系，给出不同差分格式对计算网格步长的限制。同时分析了计算中混淆误差的来源和控制方法，指出了迎风型紧致格式能很好地控制混淆误差。2．将上述格式与三阶精度的Adams半隐格式相结合，构造了不可压槽道湍流直接数值模拟的高效算法。该算法利用基于交错网格的离散形式的压力Poisson方程求解压力项，避免了压力边界条件处理的困难。利用FFT对方程中的隐式部分进行解耦，解耦后的方程采用追赶法（LU分解法）求解，大大减少了计算量。为了检验该方法，进行了三维不可压槽道湍流的直接数值模拟，得到了Re＝2800的充分发展不可压槽道湍流，并对该湍流场进行了统计分析。包括脉动速度偏斜因子在内的各阶统计量与实验结果及Kim等人的计算结果吻合十分理想，说明本方法是行之有效的。3．进行了三维充分发展的可压缩槽道湍流的直接数值模拟。得到了 Re=3300，Ma=0.8的充分发展可压槽道湍流的数据库。流场的统计特征（如等效平均速度分布，“半局部”尺度无量纲化的脉动速度均方根）和他人的数值计算结果吻合。得到了可压槽道湍流的各阶统计量，其中脉动速度的偏斜因子和平坦因子等高阶统计量尚未见其他文献报道。同时还分析了压缩性效应对壁湍流影响的机理，指出近壁处的压力－膨胀项将部分湍流脉动的动能转换成内能，使得可压湍流近壁速度条带结构更加平整。4．模拟了二维不可压槽道流动的饱和态（所谓“二维湍流”），分析了“二维槽道湍流”的非线性行为特征。分析了流场中的上抛－下扫和间歇现象，研究了“二维湍流”与三维湍流的区别。指出“二维湍流”反映了三维湍流的部分特征，同时指出了展向扰动对于湍流核心区发展的重要性。5．首次对可压缩槽道湍流及“二维槽道湍流”标度律进行了分析，得出了以下结论：a）槽道湍流中，在槽道中心线附近较宽的区域，存在标度律。b）该区域流场存在扩展自相似性（ESS）。c）在Mach数不是很高时，压缩性对标度指数影响不大。本文结果同SL标度律的理论值吻合较好，有效支持了该理论。对“二维槽道湍流”也有相似的结论，但与三维湍流不同的是，“二维槽道湍流”存在标度律的区域更宽，近壁处的标度指数比中心处有所升高。