三峡永久船闸高边坡开挖三维离散元数值模拟:3-D Simulation of the Excavation of High Steep Slope of Three-Gorges Permanent Lock by Distinct Element Method

采用面一面接触的三维离散元刚性块体模型,从实测节理面中取出其中的三组,按照其倾向、倾角和节理间距将三峡永久船闸未开挖的区域划分为10~5个离散单元,通过施加力边界条件,给出了与实测初始地应力场接近的数值模拟结果；然后,分4步模拟了永久船闸的开挖过程。计算结果表明：开挖过程会引起节理面出现张开趋势,个别岩体还会沿着节理面滑移。岩体位移的不对称现象较为自然地说明了由节理引起的岩体各向异性特征。

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.

Direct numerical simulation of transition and turbulence in compressible mixing layer

The three-dimensional compressible Navier-Stokes equations are approximated by a fifth order upwind compact and a sixth order symmetrical compact difference relations combined with three-stage Ronge-Kutta method. The computed results are presented for convective Mach number Mc = 0.8 and Re = 200 with initial data which have equal and opposite oblique waves. From the computed results we can see the variation of coherent structures with time integration and full process of instability, formation of Lambda-vortices, double horseshoe vortices and mushroom structures. The large structures break into small and smaller vortex structures. Finally, the movement of small structure becomes dominant, and flow field turns into turbulence. It is noted that production of small vortex structures is combined with turning of symmetrical structures to unsymmetrical ones. It is shown in the present computation that the flow field turns into turbulence directly from initial instability and there is not vortex pairing in process of transition. It means that for large convective Mach number the transition mechanism for compressible mixing layer differs from that in incompressible mixing layer.

Effect of Microtopography, Slope Length and Gradient, and Vegetative Cover on Overland Flow Through Simulation

Overland flow on a hillslope is significantly influenced by its microtopography, slope length and gradient, and vegetative cover. A 1D kinematic wave model in conjunction with a revised form of the Green-Ampt infiltration equation was employed to evaluate the effect of these surface conditions. The effect of these conditions was treated through the resistance parameter in the kinematic wave model. The resistance in this paper was considered to be made up of grain resistance, form resistance, and wave resistance. It was found that irregular slopes with microtopography eroded more easily than did regular slopes. The effect of the slope gradient on flow velocity and flow shear stress could be negative or positive. With increasing slope gradient, the flow velocity and shear stress first increased to a peak value, then decreased again, suggesting that there exists a critical slope gradient for flow velocity and shear stress. The vegetative cover was found to protect soil from erosion primarily by enhancing erosion-resisting capacity rather than by decreasing the eroding capability of overland flow.

Large-Eddy-Simulation of 3-Dimensional Rayleigh-Taylor Instability in Incompressible Fluids

The 3-dimensiqnal incompressible Rayleigh-Taylor instability is numerically studied through the large-eddy-simulation (LES) approach based on the passive scalar transport model. Both the instantaneous velocity and the passive scalar fields excited by sinu

A Quantitative Analysis of the Effect of Laser Transformation Hardening on Crack Driving Force in Steels

A mechanical model of a laser transformation hardening specimen with a crack in the middle of the hardened layer is developed to quantify the effects of the residual stress and hardness gradient on crack driving force in terms of J-integral. It is assumed

Simulation of Free-Surface Flow in a Tank Using the Navier-Stokes Model and Unstructured Finite Volume Method

Modelling free-surface flow has very important applications in many engineering areas such as oil transportation and offshore structures. Current research focuses on the modelling of free surface flow in a tank by solving the Navier-Stokes equation. An unstructured finite volume method is used to discretize the governing equations. The free surface is tracked by dynamically adapting the mesh and making it always surface conforming. A mesh-smoothing scheme based on the spring analogy is also implemented to ensure mesh quality throughout the computaiton. Studies are performed on the sloshing response of a liquid in an elastic container subjected to various excitation frequencies. Further investigations are also carried out on the critical frequency that leads to large deformation of the tank walls. Another numerical simulation involves the free-surface flow past as submerged obstacle placed in the tank to show the flow separation and vortices. All these cases demonstrate the capability of this numerical method in modelling complicated practical problems.

Numerical Simulation of Transient Thermal Field and Residual Stress in Laser Melting Process

以激光熔凝表面强韧化处理为背景,应用空间弹塑性有限单元和高精度数值算法同时考虑材料组织性能的变化模拟工件的温度场及残余应力,研究激光熔凝加工中瞬时温度场及残余应力数值模拟,同时考虑相变潜热及相变塑性的影响,用算例验证了模型的正确性,给出了不同时刻温度场分布及残余应力分布。

Numerical Simulation of Gaseous Detonation of H_2-O_2 Mixture with Detailed Chemical Reaction Model

采用高精度的ENO格式和基于基元化学反应的真实化学反应模型求解氢氧混合气体一维爆轰波的精细结构。采用直接起爆方法得到稳定传播的爆轰波,计算的爆轰波阵面参数和实验相当符合。对爆轰波反应区化学反应的研究表明,参与反应的不同组分具有不同类型的变化特征。网格尺寸影响的研究表明,计算结果的精度随着网格尺寸的增加而增加,并能保持较好的收敛性。移动网格研究结果表明,网格运动速度和爆轰速度接近时,两者的相互作用对计算结果产生一定影响。

Direct Numerical Simulation of Passive Scalar in Decaying Compressible Turbulence

The passive scalars in the decaying compressible turbulence with the initial Reynolds number (defined by Taylor scale and RMS velocity) Re=72, the initial turbulent Mach numbers (defined by RMS velocity and mean sound speed) Mt=0.2-0.9, and the Schmidt numbers of passive scalar Sc=2-10 are numerically simulated by using a 7th order upwind difference scheme and 8th order group velocity control scheme. The computed results are validated with different numerical methods and different mesh sizes. The Batchelor scaling with k(-1) range is found in scalar spectra. The passive scalar spectra decay faster with the increasing turbulent Mach number. The extended self-similarity (ESS) is found in the passive scalar of compressible turbulence.

Numerical simulation of standing solitons and their interaction

Standing soliton was studied by numerical simulation of ifs governing equation, a cubic Schrodiger equation with a complex conjugate term, which was derived by Miles and was accepted. The value of linear damping in Miles equation was studied. Calculations showed that linear damping effects strongly on the formation of a standing soliton and Laedke and Spatschek stable condition is only a necessary condition, but not a sufficient one. The interaction of two standing solitons was simulated. Simulations showed that the interaction pattern depends on system parameters. Calculations for the different initial condition and its development indicated that a stable standing soliton can be fanned only for proper initial disturbance, otherwise the disturbance will disappear or develop into several solitons.

Numerical Simulation of the Flow Field and Concentration Distribution in the Bulk Growth of Silicon Carbide Crystals

The physical vapor transport (PVT) method is being widely used to grow large-size single SiC crystals. The growth process is associated with heat and mass transport in the growth chamber, chemical reactions among multiple species as well as phase change at the crystal/gas interface. The current paper aims at studying and verifying the transport mechanism and growth kinetics model by demonstrating the flow field and species concentration distribution in the growth system. We have developed a coupled model, which takes into account the mass transport and growth kinetics. Numerical simulation is carried out by employing an in-house developed software based on finite volume method. The results calculated are in good agreement with the experimental observation.

Effects of Subgrid-Scale Modeling on Time Correlations in Large Eddy Simulation

The effects of the unresolved subgrid-scale (SGS) motions on the energy balance of the resolved scales in large eddy simulation (LES) have been investigated actively because modeling the energy transfer between the resolved and unresolved scales is crucial to constructing accurate SGS models. But the subgrid scales not only modify the energy balance, they also contribute to temporal decorrelation of the resolved scales. The importance of this effect in applications including the predictability problem and the evaluation of sound radiation by turbulent flows motivates the present study of the effect of SGS modeling on turbulent time correlations. This paper compares the two-point, two-time Eulerian velocity correlation in isotropic homogeneous turbulence evaluated by direct numerical simulation (DNS) with the correlations evaluated by LES using a standard spectral eddy viscosity. It proves convenient to express the two-point correlations in terms of spatial Fourier decomposition of the velocity field. The LES fields are more coherent than the DNS fields: their time correlations decay more slowly at all resolved scales of motion and both their integral scales and microscales are larger than those of the DNS field. Filtering alone is not responsible for this effect: in the Fourier representation, the time correlations of the filtered DNS field are identical to those of the DNS field itself. The possibility of modeling the decorrelating effects of the unresolved scales of motion by including a random force in the model is briefly discussed. The results could have applications to the problem of computing sound sources in isotropic homogeneous turbulence by LES