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.

Transonic Aeroelastic Numerical Simulation in Aeronautical Engineering

A lower-upper symmetric Gauss-Seidel (LU-SGS) subiteration scheme is constructed for time-marching of the fluid equations. The Harten-Lax-van Leer-Einfeldt-Wada (HLLEW) scheme is used for the spatial discretization. The same subiteration formulation is applied directly to the structural equations of motion in generalized coordinates. Through subiteration between the fluid and structural equations, a fully implicit aeroelastic solver is obtained for the numerical simulation of fluid/structure interaction. To improve the ability for application to complex configurations, a multiblock grid is used for the flow field calculation and transfinite interpolation (TFI) is employed for the adaptive moving grid deformation. The infinite plate spline (IPS) and the principal of virtual work are utilized for the data transformation between the fluid and structure. The developed code was first validated through the comparison of experimental and computational results for the AGARD 445.6 standard aeroelastic wing. Then, the flutter character of a tail wing with control surface was analyzed. Finally, flutter boundaries of a complex aircraft configuration were predicted.

A Study of Composite Beam with Shape Memory Alloy Arbitrarily Embedded under Thermal and Mechanical Loadings

The constitutive relations and kinematic assumptions on the composite beam with shape memory alloy (SMA) arbitrarily embedded are discussed and the results related to the different kinematic assumptions are compared. As the approach of mechanics of materials is to study the composite beam with the SMA layer embedded, the kinematic assumption is vital. In this paper, we systematically study the kinematic assumptions influence on the composite beam deflection and vibration characteristics. Based on the different kinematic assumptions, the equations of equilibrium/motion are different. Here three widely used kinematic assumptions are presented and the equations of equilibrium/motion are derived accordingly. As the three kinematic assumptions change from the simple to the complex one, the governing equations evolve from the linear to the nonlinear ones. For the nonlinear equations of equilibrium, the numerical solution is obtained by using Galerkin discretization method and Newton-Rhapson iteration method. The analysis on the numerical difficulty of using Galerkin method on the post-buckling analysis is presented. For the post-buckling analysis, finite element method is applied to avoid the difficulty due to the singularity occurred in Galerkin method. The natural frequencies of the composite beam with the nonlinear governing equation, which are obtained by directly linearizing the equations and locally linearizing the equations around each equilibrium, are compared. The influences of the SMA layer thickness and the shift from neutral axis on the deflection, buckling and post-buckling are also investigated. This paper presents a very general way to treat thermo-mechanical properties of the composite beam with SMA arbitrarily embedded. The governing equations for each kinematic assumption consist of a third order and a fourth order differential equation with a total of seven boundary conditions. Some previous studies on the SMA layer either ignore the thermal constraint effect or implicitly assume that the SMA is symmetrically embedded. The composite beam with the SMA layer asymmetrically embedded is studied here, in which symmetric embedding is a special case. Based on the different kinematic assumptions, the results are different depending on the deflection magnitude because of the nonlinear hardening effect due to the (large) deflection. And this difference is systematically compared for both the deflection and the natural frequencies. For simple kinematic assumption, the governing equations are linear and analytical solution is available. But as the deflection increases to the large magnitude, the simple kinematic assumption does not really reflect the structural deflection and the complex one must be used. During the systematic comparison of computational results due to the different kinematic assumptions, the application range of the simple kinematic assumption is also evaluated. Besides the equilibrium study of the composite laminate with SMA embedded, the buckling, post-buckling, free and forced vibrations of the composite beam with the different configurations are also studied and compared.

Forward-Running Detonation Drivers for High-Enthalpy Shock Tunnels

To improve the quality of driving flows generated with detonation-driven shock tunnels operated in the forward-running mode, various detonation drivers with specially designed sections were examined. Four configurations of the specially designed section, three with different converging angles and one with a cavity ring, were simulated by solving the Euler equations implemented with a pseudo kinetic reaction model. From the first three cases, it is observed that the reflection of detonation fronts at the converging wall results in an upstream-traveling shock wave that can increase the flow pressure that has decreased due to expansion waves, which leads to improvement of the driving flow. The configuration with a cavity ring is found to be more promising because the upstream-traveling shock wave appears stronger and the detonation front is less overdriven. Although pressure fluctuations due to shock wave focusing and shock wave reflection are observable in these detonation-drivers, they attenuate very rapidly to an acceptable level as the detonation wave propagates downstream. Based on the numerical observations, a new detonation-driven shock tunnel with a cavity ring is designed and installed for experimental investigation. Experimental results confirm the conclusion drawn from numerical simulations. The generated driving flow in this shock tunnel could maintain uniformity for as long as 4 ms. Feasibility of the proposed detonation driver for high-enthalpy shock tunnels is well demonstrated.

Three-Dimensional Network Model And Its Parameter Calibration

A material model, whose framework is parallel spring-bundles oriented in 3-D space, is proposed. Based on a discussion of the discrete schemes and optimum discretization of the solid angles, a 3-D network cell consisted of one-dimensional components is developed with its geometrical and physical parameters calibrated. It is proved that the 3-D network model is able to exactly simulate materials with arbitrary Poisson ratio from 0 to 1/2, breaking through the limit that the previous models in the literature are only suitable for materials with Poisson ratio from 0 to 1/3. A simplified model is also proposed to realize high computation accuracy within low computation cost. Examples demonstrate that the 3-D network model has particular superiority in the simulation of short-fiber reinforced composites.

环形激波和爆轰波会聚过程的气体动力学特性

首先应用修正的激波CCW关系式分析了激波和爆轰波会聚过程及其化学反应的影响, 得到了初始Mach数相同条件下激波聚焦效果比爆轰波强的定性结论. 然后求解Euler方程和基元化学反应模型对环形激波和爆轰波的会聚过程进行了数值模拟, 分析了会聚过程的气体动力学特性、波系结构的差异及其产生原因. 通过对波绕射、会聚和聚焦点附近状态参数变化的对比分析, 揭示了爆轰波化学反应在绕射和会聚不同阶段所起的作用, 验证了理论分析结果.

红壤农田水热动态耦合模式和观测研究

建立了红壤农田水热动态耦合模式,分析了模式中温度变化与水分运动分层的物理原因,说明了气候状况对地表面能量交换的 影响,给出了净辐射和蒸散量的计算方法,提出了有限差分计算中具有二阶精度的Euler隐式格式,介绍了红壤站的气候概况和野外观 测情况.最后利用本模式对红壤花生地陆气水热交换过程进行了数值模拟,模拟值与实测值吻合较好.

Wave dynamic processes induced by a supersonic projectile discharging from a shock tube

A numerical study on wave dynamic processes occurring in muzzle blast flows, which are created by a supersonic projectile released from the open-end of a shock tube into ambient air, is described in this paper. The Euler equations, assuming axisymmetric flows, are solved by using a dispersion-controlled scheme implemented with moving boundary conditions. Three test cases are simulated for examining friction effects on the muzzle flow. From numerical simulations, the wave dynamic processes, including two blast waves, two jet flows, the bow shock wave and their interactions in the muzzle blasts, are demonstrated and discussed in detail. The study shows that the major wave dynamic processes developing in the muzzle flow remain similar when the friction varies, but some wave processes, such as shock-shock interactions, shock-jet interactions and the contact surface instability, get more intensive, which result in more complex muzzle blast flows.

Three-Dimensional Propagation of the Transmitted Shock Wave in a Square Cross-Sectional Chamber

An investigation into the three-dimensional propagation of the transmitted shock wave in a square cross-section chamber was described in this paper, and the work was carried out numerically by solving the Euler equations with a dispersion-controlled scheme. Computational images were constructed from the density distribution of the transmitted shock wave discharging from the open end of the square shock tube and compared directly with holographic interferograms available for CFD validation. Two cases of the transmitted shock wave propagating at different Mach numbers in the same geometry were simulated. A special shock reflection system near the corner of the square cross-section chamber was observed, consisting of four shock waves: the transmitted shock wave, two reflection shock waves and a Mach stem. A contact surface may appear in the four-shock system when the transmitted shock wave becomes stronger. Both the secondary shock wave and the primary vortex loop are three-dimensional in the present case due to the non-uniform flow expansion behind the transmitted shock.

An Extension of 2D Janbu’s Generalized Procedure of Slices for 3D Slope Stability Analysis II—Numerical Method and Applications

This paper provides a numerical approach on achieving the limit equilibrium method for 3D slope stability analysis proposed in the theoretical part of the previous paper. Some programming techniques are presented to ensure the maneuverability of the method. Three examples are introduced to illustrate the use of this method. The results are given in detail such as the local factor of safety and local potential sliding direction for a slope. As the method is an extension of 2D Janbu's generalized procedure of slices (GPS), the results obtained by GPS for the longitudinal sections of a slope are also given for comparison with the 3D results. A practical landslide in Yunyang, the Three Gorges, of China, is also analyzed by the present method. Moreover, the proposed method has the advantages and disadvantages of GPS. The problem frequently encountered in calculation process is still about the convergency, especially in analyzing the stability of a cutting corner. Some advice on discretization is given to ensure convergence when the present method is used. However, the problem about convergency still needs to be further explored based on the rigorous theoretical background.

高雷诺数流动的控制方程体系和扩散抛物化Navier-Stokes方程组的意义和用途

在计算机发达的时代,高雷诺(Re)数绕流计算中有无必要使用简化NS方程组,本文讨论这个问题.主要内容如下:(1)高Re数绕流包含3种基本流动:所有方向对流占优流动、所有方向对流扩散竞争流动和部分方向对流占优部分方向对流扩散竞争流动(简称干扰剪切流动),3个基本流动的特征彼此不同且在流场中所占领域大小彼此相差悬殊,NS方程区域很小,它们的最简单控制方程组Euler、Navier-Stokes(NS)和扩散抛物化(DP)NS方程组的数学性质彼此不同,因此利用Euler-DPNS-NS方程组体系分析计算高Re数绕流流动就是一个合乎逻辑的选择,该法与利用单一NS方程组的常用方法可以彼此检验和补充.(2)流体之间以及流体与外界的动量、能量和质量交换,流态从层流到湍流的演化主要发生在干扰剪切流动中,干扰剪切流及其最简单控制方程--DPNS方程组具有基础意义;DPNS方程组笔者在1967年已提出.(3)诸简化NS方程组:DPNS、抛物化(P)NS、薄层(TL)NS、黏性层(VL)NS方程组的发展、相互关系,它们的历史贡献和今后的用途;它们的数学性质均为扩散抛物型,但它们包含的黏性项彼此有所不同;从流体力学角度来看,它们中只有DPNS方程组能够准确描述干扰剪切流动.提出把诸简化NS方程组统一为DPNS方程组的建议.(4)干扰剪切流--DPNS方程组与无干扰剪切流--边界层方程组之间的关系以及进一步研究干扰剪切流的意义.

油水分离技术研究之一--直管和螺旋管的数值模拟

本文采用数值模拟的方法研究油水混合物在直管和螺旋管中的流动状况.计算采用Euler-Euler法和Euler模型:支配油水分离两相流动的基本方程包括连续方程和动量方程,湍流模型采用多相流中混合型k-ε模式,基本方程的离散和求解采用SIMPLE算法.利用Fluent软件,以直管和螺旋管为例进行了计算,获得了初步计算结果.计算表明,本文所用方法可以较好地模拟直管中油水在重力作用下的分离,以及螺旋分离器中油水在重力和离心力共同作用下的分离现象.并可为实验研究提供参考。

Shocked Flows Induced by Supersonic Projectiles Moving in Tubes

A numerical study on shocked flows induced by a supersonic projectile moving in tubes is described in this paper. The dispersion-controlled scheme was adopted to solve the Euler equations implemented with moving boundary conditions. Four test cases were carried out in the present study: the first two cases are for validation of numerical algorithms and verification of moving boundary conditions, and the last two cases are for investigation into wave dynamic processes induced by the projectile moving at Mach numbers of M-p = 2.0 and 2.4, respectively, in a short time duration after the projectile was released from a shock tube into a big chamber. It was found that complex shock phenomena exist in the shocked flow, resulting from shock-wave/projectile interaction, shock-wave focusing, shock-wave reflection and shock-wave/contact-surface interactions, from which turbulence and vortices may be generated. This is a fundamental study on complex shock phenomena, and is also a useful investigation for understanding on shocked flows in the ram accelerator that may provide a highly efficient facility for launching hypersonic projectiles.

脉冲爆轰发动机热射流起爆机理数值分析

应用频散可控耗散差分格式,求解具有化学反应项的Euler方程,探讨了热射流起爆可燃混合气缩短DDT过程的物理机制.数值研究模拟了不同条件下的起爆过程,从氢氧链式反应出发详细分析了氢氧爆轰直接起爆的SWACER（能量释放而形成激波或压缩波的相干放大）机制的建立条件,讨论了热射流起爆存在超临界、临界和亚临界三种直接起爆机制.