求解声波动方程的隐格式有限体积法

为研究声传播问题,提出一种声波动方程的隐格式有限体积法,该方法将格点型有限体积法与Newmark格式相结合.模拟平面波的传播过程,对比分析隐格式有限体积法和文献中显格式有限体积法的精度、稳定性及计算消耗等方面的性能.数值结果表明:当λ/Δx≥10时,两种算法均能得到满足精度要求的解;采用无条件稳定的隐格式算法,当满足ω0Δt≤0.3时,预测声压的相对峰值误差1%;当采用相同时间、空间步长时,隐格式算法精度高于显格式算法;隐格式算法对吸收边界的处理精度高于显格式算法,但对全反射边界的处理精度低于显格式算法;两种算法内存消耗比较接近,显格式算法的CPU耗时较少.     

Global Avalanche Characteristics of Boolean Functions by Concatenation

In order to measure the correlation propeties of two Boolean functions,the global avalanche characteristics of Boolean functions constructed by concatenation are discussed,i.e.,f_1‖f_2and f_1‖f_2‖f_3‖f_4.Firstly,for the function f = f_1‖f_2,the cross-correlation function of f_1,f_2 in the special condition are studied.In this case,f,f_1,f_2 must be in desired form.By computing their sum-of-squares indicators,the crosscorrelation function between f_1,f_2 is obtained.Secondly,for the function g = f_1‖f_2‖f_3‖f_4,by analyzing the relation among their auto-correlation functions,their sum-of-squares indicators are investigated.Based on them,the sum-of-squares indicators of functions obtained by Canteaut et al.are investigated.The results show that the correlation property of g is good when the correlation properties of Boolean functions f_1,f_2,f_3,f_4 are good.     

Research and Application of Dynamic Equation for Full Frontal Impact

Full frontal impact theory needs researching and exploring to satisfy the primary safety design of occupant restraint system,avoiding the increasingly "engineering"trend in order to develop and design safety vehicle. After occupant restraint system is simulated by using linear elastic stiffness k,the occupant-vehicle frontal rigid barrier impact model is established. Dynamic equation of dummy chest coupling vehicle is built for full frontal impact based on ordinary vehicle deceleration by Hooke law,and the equation is solved by comparing coefficient and satisfying boundary qualifications. While relative vehicle characteristic parameters are kept unchanging,the actual vehicle deceleration is fitted to the simplified equivalent square wave( ESW),tipped equivalent square wave( TESW) and equivalent dual trapezoids wave( EDTW). Phase angle  and amplitude A of dynamic equations based on ESW,TESW and EDTW are calculated and deduced. The results show that: the dynamic equation of dummy chest coupling vehicle can be well utilized to instruct the primary safety design of full frontal impact for objective vehicle to satisfy chest deceleration demands and the equation based on TESW is best for this design.     

基于解析解的四阶PDE曲面裁剪方法

该文提出了一种可以广泛应用于四阶PDE曲面的裁剪方法。利用PDE的参数域内的曲线在曲面上投影,得到所求裁剪曲面的边界曲线,然后通过边界曲线的导矢与曲面在边界曲线处的法向量得到边界曲线处的跨界导矢,最后以求得的裁剪曲面的边界曲线以及裁剪曲面在边界曲线处的跨界导矢为PDE曲面的边界条件,用四阶的PDE曲面方程求得PDE裁剪曲面。     

Method of calculating the collision integral and solution of the Boltzmann kinetic equation for simple gases,gas mixtures and gases with rotational degrees of freedom

The conservative method of calculating the Boltzmann collision integral for simple gases, gas mixtures and gases with rotational degrees of freedom of molecules is presented. In all cases the common approach based on the projection technique of summing up the contributions in the collision integral is used. The method is applied for solving the Boltzmann kinetic equation for two fundamental rarefied gas flow problems: the heat transfer problem and the problem of shock wave structure. A comparison with experimental and numerical data of other authors is reported. It is shown that the considered numerical method allows one to solve the Boltzmann equation for real gases with high accuracy.     

A new high accuracy method for two-dimensional biharmonic equation with nonlinear third derivative terms: application to Navier–Stokes equations of motion

In this paper, we propose a new compact fourth-order accurate method for solving the two-dimensional fourth-order elliptic boundary value problem with third-order nonlinear derivative terms. We use only 9-point single computational cell in the scheme. The proposed method is then employed to solve Navier–Stokes equations of motion in terms of streamfunction–velocity formulation, and the lid-driven square cavity problem. We describe the derivation of the method in details and also discuss how our streamfunction–velocity formulation is able to handle boundary conditions in terms of normal derivatives. Numerical results show that the proposed method enables us to obtain oscillation-free high accuracy solution.     

Compact difference method for solving the fractional reaction–subdiffusion equation with Neumann boundary value condition

In this paper, we derive a high-order compact finite difference scheme for solving the reaction–subdiffusion equation with Neumann boundary value condition. The L1 method is used to approximate the temporal Caputo derivative, and the compact difference operator is applied for spatial discretization. We prove that the compact finite difference method is unconditionally stable and convergent with order O(τ^2?α+h⁴) in L₂ norm, where τ, α, and h are the temporal step size, the order of time fractional derivative and the spatial step size, respectively. Finally, some numerical experiments are carried out to show the effectiveness of the proposed difference scheme.     

Numerical solution of nonlinear system of Klein–Gordon equations by cubic B-spline collocation method

A technique to approximate the solutions of nonlinear Klein–Gordon equation and Klein–Gordon-Schrödinger equations is presented separately. The approach is based on collocation of cubic B-spline functions. The above-mentioned equations are decomposed into a system of partial differential equations, which are further converted to an amenable system of ODEs. The obtained system has been solved by SSP-RK54 scheme. Numerical solutions are presented for five examples, to show the accuracy and usefulness of proposed approach. The approximate solutions of both the equations are computed without using any transformation and linearization. The technique can be applied with ease to solve linear and nonlinear PDEs and also reduces the computational work.     

Continuous-time orbit problems are decidable in polynomial-time

We place the continuous-time orbit problem in P$\mathsf{P}$, sharpening the decidability result shown by Hainry [7].