An Energy Conserving Finite Difference Scheme for Simulation of Collisions

Nonlinear phenomena play an essential role in the sound production process of many musical instruments. A common source of these effects is object collision, the numerical simulation of which is known to give rise to stability
issues. This paper presents a method to construct numerical schemes that conserve the total energy in simulations of one-mass systems involving collisions, with no conditions imposed on any of the physical or numerical parameters.
This facilitates the adaptation of numerical models to experimental data, and allows a more free parameter adjustment in sound synthesis explorations. The energy preservedness of the proposed method is tested and demonstrated though several examples, including a bouncing ball and a non-linear oscillator, and implications regarding the wider applicability are discussed.

An efficient finite difference scheme for three-dimensional, non-equilibrium flows using operator splitting

An efficient finite difference scheme is presented for the inviscid terms of the three-dimensional, compressible flow equations for chemical non-equilibrium gases. This scheme represents an extension and an improvement of one proposed by the author, and includes operator splitting.

A TVD finite difference scheme with non-uniform meshes and without upstream weighting

We present a finite difference scheme, with the TVD (total variation diminishing) property, for scalar conservation laws. The scheme applies to non-uniform meshes, allowing for variable mesh spacing, and is without upstream weighting. When applied to systems of conservation laws, no scalar decomposition is required, nor are any artificial tuning parameters, and this leads to an efficient, robust algorithm.

A finite difference scheme for inviscid flows with non-eequilibrium chemistry and internal energy

A finite difference scheme is presented for the inviscid terms of the equations of compressible fluid dynamics with general non-equilibrium chemistry and internal energy.

A finite difference scheme for steady, supersonic, two-dimensional, compressible flow of real gases

A finite difference scheme is presented for the solution of the two-dimensional equations of steady, supersonic, compressible flow of real gases. The scheme incorparates numerical characteristic decomposition, is shock-capturing by design and incorporates space-marching as a result of the assumption that the flow is wholly supersonic in at least one space dimension. Results are shown for problems involving oblique hydraulic jumps and reflection from a wall.

FTCS finite difference scheme GPGPU parallel computing for the heat conduction equation = Programación en paralelo GPGPU del método en diferencias finitas FTCS para la ecuación del calor

En el presente artículo se muestran las ventajas de la programación en paralelo resolviendo numéricamente la ecuación del calor en dos dimensiones a través del método de diferencias finitas explícito centrado en el espacio FTCS. De las conclusiones de este trabajo se pone de manifiesto la importancia de la programación en paralelo para tratar problemas grandes, en los que se requiere un elevado número de cálculos, para los cuales la programación secuencial resulta impracticable por el elevado tiempo de ejecución. En la primera sección se describe brevemente los conceptos básicos de programación en paralelo. Seguidamente se resume el método de diferencias finitas explícito centrado en el espacio FTCS aplicado a la ecuación parabólica del calor. Seguidamente se describe el problema de condiciones de contorno y valores iniciales específico al que se va a aplicar el método de diferencias finitas FTCS, proporcionando pseudocódigos de una implementación secuencial y dos implementaciones en paralelo. Finalmente tras la discusión de los resultados se presentan algunas conclusiones. In this paper the advantages of parallel computing are shown by solving the heat conduction equation in two dimensions with the forward in time central in space (FTCS) finite difference method. Two different levels of parallelization are consider and compared with traditional serial procedures. We show in this work the importance of parallel computing when dealing with large problems that are impractical or impossible to solve them with a serial computing procedure. In the first section a summary of parallel computing approach is presented. Subsequently, the forward in time central in space (FTCS) finite difference method for the heat conduction equation is outline, describing how the heat flow equation is derived in two dimensions and the particularities of the finite difference numerical technique considered. Then, a specific initial boundary value problem is solved by the FTCS finite difference method and serial and parallel pseudo codes are provided. Finally after results are discussed some conclusions are presented.

Direct Numerical Simulation Of Three-Dimensional Richtmyer-Meshkov Instability

Direct numerical simulation (DNS) is used to study flow characteristics after interaction of a planar shock with a spherical media interface in each side of which the density is different. This interfacial instability is known as the Richtmyer-Meshkov (R-M) instability. The compressible Navier-Stoke equations are discretized with group velocity control (GVC) modified fourth order accurate compact difference scheme. Three-dimensional numerical simulations are performed for R-M instability installed passing a shock through a spherical interface. Based on numerical results the characteristics of 3D R-M instability are analysed. The evaluation for distortion of the interface, the deformation of the incident shock wave and effects of refraction, reflection and diffraction are presented. The effects of the interfacial instability on produced vorticity and mixing is discussed.

可压缩流气动声场的数值模拟

气动声学是一门流动力学和声学之间的交叉学科，主要研究流动及其与物体相互作用产生噪声的机理。动用计算技术研究气动声学问题的手段称为计算气动声学。本文的目的是，基于高精度数值算法的研究，分别运用Lighthill比拟理论、Kirchhoff积分和直接数值模拟等方法，针对翼型绕流、激波-涡干扰和轴对称射流，研究了物面非定常脉动压力、涡脱落、激波-涡干扰以及涡对并等产生噪声的机理。首先针对声场与主流场在能级和特征尺度等方面的差异，从空间离散角度分析了几种差分格式，表明迎风紧致格式/对称紧致格式有较小的数值色散、耗散和各向异性误差，因而适用于气动噪声的计算。以Runge-Kutta格式为例，对时间离散带来的误差进行了分析。指出对声波计算来说，仅考虑格式稳定性是不够的，时间步长还受到允许色散误差和耗散误差的限制。基于保色戎关系的思想，构造了优化Runge-Kutta格式。处例显示优化Runge-Kutta格式相对于经典格式有更高的计算效率。采用3阶迎风紧致格式和3阶Runge-Kutta格式数值模拟了NACA0012翼型的可压缩非定常绕流流场，并将此流场作为近场声源，运用声学比拟理论对偶极子声和四极子声进行研究。结果指出，主流速度对远场声压有决定性影响，在来流马赫数较大时，四极子噪声和偶极子噪声具有相同量级，不能被忽略，表明了可压缩效应对声场的影响。采用5阶迎风紧致格式和4阶Runge-Kutta格式求解非定常可压缩Navier-Stokes方程，对激波-单涡/双涡干扰导致的声场进行了直接数值模拟。详细研究了激波-涡干扰产生噪声的机理，指出噪声的产生及其性质和激波变形密切相关。研究了近场噪声衰减和传播距离r的关系，发现噪声衰减大致和r~(4/5)而不是r~(1/2)成反比关系，提出这种差异是由流场的非线性效应引起的。构造了Kirchhoff积分和非定常流动计算相结合的算法。采用5阶迎风紧致格式和3阶Runge-Kutta格式对亚声速轴对称射流进行直接数值模拟。将射流流场作为近场声源，结合Kirchhoff方法求解远场 气动噪声。数值结果表明远场噪声具有方向性，噪声声压在离开对称轴20°处达到最大值。随着传播距离增大，噪声方向性逐渐减弱。

可压平面混合层稳定性分析及数值模拟

可压平面混合层是包含复杂多时空尺度运动的非定常流体力学部问题，具有深刻的理论意义和广泛的应用背景。针对该问题所涉及内容的多面性，本文的目的是，基于高精度、高分辨率数值算法的构造、发展和数值行为分析，采用线性稳定性分析和直接数值模拟方法。从理论和计算两方面集中研究压缩性效应、粘性效应、初值效应以及燃烧反应放热效应等对可压平面混合层早期稳定性行为和大尺度拟序涡结构非线性演化的影响。以混合层已有研究成果的分析和综述为开端，论文主体共包括四部分：第一部分是可压平面混合层时间/空间模式数值线性稳定性分析。实现了高精度对称紧致差分格式（SCD）对可压粘性扰动线性稳定性边值问题的求解，对导出的线性和非线性离散特征值问题，提出了两个高效局部解法。研究涉及二维/三维扰动波、无粘/粘性扰动波、特征函数和特征值谱、第一/第二模态、超声速快/慢模态、速度比和密度比等。验证了对流Mach数Mc为一个合理的压缩性参数。指出压缩性效应和粘性效应对最不稳定扰动波的波数（频率）和增长率呈相拟的抑制作用，且时间模式稳定性分析结果在许多方面是可信的。从随机和线性扰动场出发，采用高精度五阶迎风紧致和六阶对称紧致混合差分算法（UCD5/SCD6）对可压平面混合层的稳定性特征进行了直接数值模拟，揭示了初始主导线性扰动与一些实际涡结构非线性作用形态间的内在关联，印证了线性稳定性分析方法的合理性和有效性。第二部分是高精度迎风紧致差分格式（UCD）时空全离散数值行为分析。导出了其一维/二维一般色散表达式。研究表明，UCD格式在高波数区具有内在的全离散耗散和色散特性；其数值群速度的快/慢特征可因CFL数不同而改变；在稳定CFL数下简单附加人工粘性可强化UCD格式在高波数区的耗散量；提高时间精度可放宽稳定CFL数限制；UCD格式的二维全离散色散介质中存在三个不同性质的数值波，其全离散稳定性由数值声波主控。第三部分实现了高精度UCD5/SCD6差分算法对空间发展可压平面混合层的直接数值模拟。通过亚谐扰动波的个数和扰动频率的控制，捕捉到了基频涡的饱和、一次和二次对并等现象，显示了大尺度涡结构与入中初始扰动方式之间的内在联系。利用参数Mc观察了压缩性效应对大尺度涡空间演化及其相互作用的影响。第四部分实现了高精度UCD5/SCD6差分算法对非预混扩散火焰化学反应平面混合层的直接数值模拟。研究指出，放热效应可抑制和延迟涡的形成，使基频涡卷拉伸甚至丧失，混合层Reynolds 应力ρu'v'和流向速度波动关联项u'v'下降，以致涡结构与外流动量交换和标量输运减少，脉动输运能力被削弱，从而混合效率、产物生成率和混合层增长率下降，放热主要通过膨胀效应和斜压效应来抑制大尺度涡的演化。

Composite scheme using localized relaxation with non-standard finite difference method for hyperbolic conservation laws

Non-standard finite difference methods (NSFDM) introduced by Mickens [Non-standard Finite Difference Models of Differential Equations, World Scientific, Singapore, 1994] are interesting alternatives to the traditional finite difference and finite volume methods. When applied to linear hyperbolic conservation laws, these methods reproduce exact solutions. In this paper, the NSFDM is first extended to hyperbolic systems of conservation laws, by a novel utilization of the decoupled equations using characteristic variables. In the second part of this paper, the NSFDM is studied for its efficacy in application to nonlinear scalar hyperbolic conservation laws. The original NSFDMs introduced by Mickens (1994) were not in conservation form, which is an important feature in capturing discontinuities at the right locations. Mickens [Construction and analysis of a non-standard finite difference scheme for the Burgers–Fisher equations, Journal of Sound and Vibration 257 (4) (2002) 791–797] recently introduced a NSFDM in conservative form. This method captures the shock waves exactly, without any numerical dissipation. In this paper, this algorithm is tested for the case of expansion waves with sonic points and is found to generate unphysical expansion shocks. As a remedy to this defect, we use the strategy of composite schemes [R. Liska, B. Wendroff, Composite schemes for conservation laws, SIAM Journal of Numerical Analysis 35 (6) (1998) 2250–2271] in which the accurate NSFDM is used as the basic scheme and localized relaxation NSFDM is used as the supporting scheme which acts like a filter. Relaxation schemes introduced by Jin and Xin [The relaxation schemes for systems of conservation laws in arbitrary space dimensions, Communications in Pure and Applied Mathematics 48 (1995) 235–276] are based on relaxation systems which replace the nonlinear hyperbolic conservation laws by a semi-linear system with a stiff relaxation term. The relaxation parameter (λ) is chosen locally on the three point stencil of grid which makes the proposed method more efficient. This composite scheme overcomes the problem of unphysical expansion shocks and captures the shock waves with an accuracy better than the upwind relaxation scheme, as demonstrated by the test cases, together with comparisons with popular numerical methods like Roe scheme and ENO schemes.

Difference schemes on non-uniform mesh and their application

High order accurate schemes are needed to simulate the multi-scale complex flow fields to get fine structures in simulation of the complex flows with large gradient of fluid parameters near the wall, and schemes on non-uniform mesh are desirable for many CFD (computational fluid dynamics) workers. The construction methods of difference approximations and several difference approximations on non-uniform mesh are presented. The accuracy of the methods and the influence of stretch ratio of the neighbor mesh increment on accuracy are discussed. Some comments on these methods are given, and comparison of the accuracy of the results obtained by schemes based on both non-uniform mesh and coordinate transformation is made, and some numerical examples with non-uniform mesh are presented.

Upwind compact method and direct numerical-simulation of driven flow in a square cavity

A high-order accurate finite-difference scheme, the upwind compact method, is proposed. The 2-D unsteady incompressible Navier-Stokes equations are solved in primitive variables. The nonlinear convection terms in the governing equations are approximated by using upwind biased compact difference, and other spatial derivative terms are discretized by using the fourth-order compact difference. The upwind compact method is used to solve the driven flow in a square cavity. Solutions are obtained for Reynolds numbers as high as 10000. When Re less than or equal to 5000, the results agree well with those in literature. When Re = 7500 and Re = 10000, there is no convergence to a steady laminar solution, and the flow becomes unsteady and periodic.

Upwind compact scheme with dispersion control for aerodynamic equations

A compact upwind scheme with dispersion control is developed using a dissipation analogy of the dispersion term. The term is important in reducing the unphysical fluctuations in numerical solutions. The scheme depends on three free parameters that may be used to regulate the size of dissipation as well as the size and direction of dispersion. A coefficient to coordinate the dispersion is given. The scheme has high accuracy, the method is simple, and the amount of computation is small. It also has a good capability of capturing shock waves. Numerical experiments are carried out with two-dimensional shock wave reflections and the results are very satisfactory.

High order unconditionally stable difference schemes for the Riesz space-fractional telegraph equation

In this paper, a class of unconditionally stable difference schemes based on the Pad´e approximation is presented for the Riesz space-fractional telegraph equation. Firstly, we introduce a new variable to transform the original dfferential equation to an equivalent differential equation system. Then, we apply a second order fractional central difference scheme to discretise the Riesz space-fractional operator. Finally, we use (1, 1), (2, 2) and (3, 3) Pad´e approximations to give a fully discrete difference scheme for the resulting linear system of ordinary differential equations. Matrix analysis is used to show the unconditional stability of the proposed algorithms. Two examples with known exact solutions are chosen to assess the proposed difference schemes. Numerical results demonstrate that these schemes provide accurate and efficient methods for solving a space-fractional hyperbolic equation.

Unsteady flow and heat transfer between rotating coaxial disks

A study is made on the flow and heat transfer of a viscous fluid confined between two parallel disks. The disks are allowed to rotate with different time dependent angular velocities, and the upper disk is made to approach the lower one with a constant speed. Numerical solutions of the governing parabolic partial differential equations are obtained through a fourth-order accurate compact finite difference scheme. The normal forces and torques that the fluid exerts on the rotating surfaces are obtained at different nondimensional times for different values of the rate of squeezing and disk angular velocities. The temperature distribution and heat transfer are also investigated in the present analysis.