Numerical solution of the incompressible Navier-Stokes equations with an upwind compact difference scheme

A new finite difference method for the discretization of the incompressible Navier-Stokes equations is presented. The scheme is constructed on a staggered-mesh grid system. The convection terms are discretized with a fifth-order-accurate upwind compact difference approximation, the viscous terms are discretized with a sixth-order symmetrical compact difference approximation, the continuity equation and the pressure gradient in the momentum equations are discretized with a fourth-order difference approximation on a cell-centered mesh. Time advancement uses a three-stage Runge-Kutta method. The Poisson equation for computing the pressure is solved with preconditioning. Accuracy analysis shows that the new method has high resolving efficiency. Validation of the method by computation of Taylor's vortex array is presented.

Dynamic Characteristics of Spherically Converging Detonation Waves

The spherically converging detonation wave was numerically investigated by solving the one-dimensional multi-component Euler equations in spherical coordinates with a dispersion-controlled dissipative scheme. Finite rate and detailed chemical reaction models were used and numerical solutions were obtained for both a spherical by converging detonation in a stoichiometric hydrogen-oxygen mixture and a spherically focusing shock in air. The results showed that the post-shock pressure approximately arises to the same amplitude in vicinity of the focal point for the two cases, but the post-shock temperature level mainly depends on chemical reactions and molecular dissociations of a gas mixture. While the chemical reaction heat plays an important role in the early stage of detonation wave propagation, gas dissociations dramatically affect the post-shock flow states near the focal point. The maximum pressure and temperature, non-dimensionalized by their initial value, are approximately scaled to the propagation radius over the initial detonation diameter. The post-shock pressure is proportional to the initial pressure of the detonable mixture, and the post-shock temperature is also increased with the initial pressure, but in a much lower rate than that of the post-shock pressure.

Elasticity, stability, and ideal strength of beta-SiC in plane-wave-based ab initio calculations

On the basis of the pseudopotential plane-wave method and the local-density-functional theory, this paper studies energetics, stress-strain relation, stability, and ideal strength of beta-SiC under various loading modes, where uniform uniaxial extension and tension and biaxial proportional extension are considered along directions [001] and [111]. The lattice constant, elastic constants, and moduli of equilibrium state are calculated and the results agree well with the experimental data. As the four SI-C bonds along directions [111], [(1) over bar 11], [11(1) over bar] and [111] are not the same under the loading along [111], internal relaxation and the corresponding internal displacements must be considered. We find that, at the beginning of loading, the effect of internal displacement through the shuffle and glide plane diminishes the difference among the four Si-C bonds lengths, but will increase the difference at the subsequent loading, which will result in a crack nucleated on the {111} shuffle plane and a subsequently cleavage fracture. Thus the corresponding theoretical strength is 50.8 GPa, which agrees well with the recent experiment value, 53.4 GPa. However, with the loading along [001], internal relaxation is not important for tetragonal symmetry. Elastic constants during the uniaxial tension along [001] are calculated. Based on the stability analysis with stiffness coefficients, we find that the spinodal and Born instabilities are triggered almost at the same strain, which agrees with the previous molecular-dynamics simulation. During biaxial proportional extension, stress and strength vary proportionally with the biaxial loading ratio at the same longitudinal strain.

Capillary Effect on Vertically Excited Surface Wave in Circular Cylindrical Vessel

In a vertically oscillating circular cylindrical container, singular perturbation theory of two-time scale expansions was developed in inviscid fluids to investigate the motion of single free surface standing wave including the effect of surface tension.

Interferometric Measurement of Density in Nonstationary Shock Wave Reflection Flow and Comparison with CFD

We present density measurements from the application of interferometry and Fourier transform fringe analysis to the problem of nonstationary shock wave reflection over a semicircular cylinder and compare our experimental measurements to theoretical results from a CFD simulation of the same problem. The experimental results demonstrate our ability to resolve detailed structure in this complex shock wave reflection problem, allowing visualization of multiple shocks in the vicinity of the triple point, plus visualization of the shear layer and an associated vortical structure. Comparison between CFD and experiment show significant discrepancies with experiment producing a double Mach Reflection when CFD predicts a transitional Mach reflection.

Comparison of the Dynamic Response of a Saturated Soil Between Two Models

By comparing the dynamic responses of saturated soil to Biot's and Yamamoto's models, the properties of the two models have be pointed out. First of all, an analysis has been made for energy loss of each model from the basic equations. Then the damping of elastic waves in coarse sand and fine sand with loading frequency and soil's parameters have been calculated and the representation of viscous friction and Coulomb friction in the two models has been concluded. Finally, the variations of loading wave damping and stress phase angles with water depth and soil's parameters have been obtained as loading waves range in ocean waves.

The cellular structure of a two-dimensional H2/O2/Ar detonation wave

In this paper, the cellular structure of a two-dimensional detonation wave in a low pressure H2/O2/Ar mixture calculated with a detailed chemical reaction model, high order scheme and high resolution grids is investigated. The regular cellular structure is produced about 1 ms after introducing perturbations in the reaction zone of a steady one-dimensional detonation wave. It is found from the present resolution study that the discrepancies concerning the structure type arising from the coarser grid employed can be resolved using a sufficiently fine grid size of 0.05 mm and below and shows a double-Mach-like strong-type configuration. During the structure evolution process, the structure configuration does not change much in the periods before and after the triple point collision. Through the triple point collision, three regular collision processes are observed and are followed by a quick change to the double-Mach-like configuration. The simulated structure tracks show that there are three different tracks associated with different triple points or the kink on the transverse wave. Comparisons with previous work and experiments indicate the presence of a strong structure for an ordinary detonation.

Analytical Interaction of the Acoustic Wave and Turbulent Flame

A modified resonance model of a weakly turbulent flame in a high-frequency acoustic wave is derived analytically. Under the mechanism of Darrieus-Landau instability, the amplitude of flame wrinkles, which is as functions of the expansion coefficient and the perturbation wave number, increases greatly independent of the 'stationary' turbulence. The high perturbation wave number makes the resonance easier to be triggered but weakened with respect to the extra acoustic wave. In a closed burning chamber with the acoustic wave induced by the flame itself, the high perturbation wave number is to restrain the resonance for a realistic flame.

Characteristics and Sub-Characteristics of Two-Dimensional Parabolized Stability Equations

特征分析表明：对原始扰动量的抛物化稳定性方程组(PSE),它在亚超音速区分别具有椭圆和抛物特性,给出PSE特征对马赫数的依赖关系,阐明PSE仅把信息对流-扩散传播特性抛物化,而保留了信息对流-扰动传播特性,因此PSE应称为扩散抛物化稳定性方程(DPSE)。

Comparison of existing poroelastic models for wave damping in a porous seabed

Mechanism of wave-seabed interaction has been extensively studied by coastal geotechnical engineers in recent years. Numerous poro-elastic models have been proposed to investigate the mechanism of wave propagation on a seabed in the past. The existing poro-elastic models include drained model, consolidation model, Coulomb-damping model, and full dynamic model. However, to date, the difference between the existing models is unclear. In this paper, the fully dynamic poro-elastic model for the wave-seabed interaction will be derived first. Then, the existing models will be reduced from the proposed fully dynamic model. Based on the numerical comparisons, the applicable range of each model is also clarified for the engineering practice.

Perturbational Finite Volume Method For The Solution of 2-D Navier-Stokes Equations on Unstructured and Structured Colocated Meshes

Based on the first-order upwind and second-order central type of finite volume( UFV and CFV) scheme, upwind and central type of perturbation finite volume ( UPFV and CPFV) schemes of the Navier-Stokes equations were developed. In PFV method, the mass fluxes of across the cell faces of the control volume (CV) were expanded into power series of the grid spacing and the coefficients of the power series were determined by means of the conservation equation itself. The UPFV and CPFV scheme respectively uses the same nodes and expressions as those of the normal first-order upwind and second-order central scheme, which is apt to programming. The results of numerical experiments about the flow in a lid-driven cavity and the problem of transport of a scalar quantity in a known velocity field show that compared to the first-order UFV and second-order CFV schemes, upwind PFV scheme is higher accuracy and resolution, especially better robustness. The numerical computation to flow in a lid-driven cavity shows that the under-relaxation factor can be arbitrarily selected ranging from 0.3 to 0. 8 and convergence perform excellent with Reynolds number variation from 102 to 104.

Analysis of Wave-Induced Effective Stresses of Seabed

波浪作用下海床的稳定性分析是海洋工程地质评价的重要内容。海床的稳定性可通过计算分析其随时间变化的有效应力场来评估。建议了一个周期载荷作用下土体的本构模型,并用于计算波浪作用下海床的应力与变形。采用Biot固结理论和有限单元法,分析了海床的动态应力场与孔隙水压力场。波浪作用下两种渗透系数时有效应力的动态变化过程结果对比,反映了渗透消散作用对海床有效应力变化的影响。

Investigation of Dynamic Fracture Toughness under Compressive-Shear Combinede Stress Wave Loading

提出亚微秒单脉冲应力波载荷作用下II型裂纹的平板冲击实验技术。加载率为dK/dt-10~8MPa·m~{”/d}·s~{-1}。实验中由锰铜应力片和弹性波理论分别测定和计算了压应力；通过微观分析确定了动态裂纹的平均扩展长度；引进等效应力强度因子,用动态断裂理论确定了60号钢的动态断裂韧性K_{Id}和K_{IId}；建立了亚微秒冲击载荷作用下确定材料动态断裂韧性的方法。

Stability and Liqueafaction Analysis of Seabed to Wave Loading

波浪作用下海床的稳定性与液化分析是海底管线、防波堤和海洋平台设计中必须仔细考虑的问题。推荐了一个循环荷载作用下土体的弹塑性实用本构模型,并给出了一种粉土的模型参数。该模型直接根据初始应力状态和循环应力的大小与作用时间计算土体的塑性应变增量,在有限元计算中不需要引入弹塑性矩阵。采用Biot理论和有限单元法,计算了粉土海床在波浪作用下的孔隙水压力和有效应力的变化过程,并对海床的稳定性和液化进行了分析。计算结果与波浪槽实验反映的规律是相符的。

Stress Wave Propagation in a Gradient Elastic Medium

The gradient elastic constitutive equation incorporating the second gradient of the strains is used to determine the monochromatic elastic plane wave propagation in a gradient infinite medium and thin rod. The equation of motion, together with the internal material length, has been derived. Various dispersion relations have been determined. We present explicit expressions for the relationship between various wave speeds, wavenumber and internal material length.