Large Eddy Simulation of Complex Turbulent Flows: Physical Aspects and Research Trends

In the current paper, we have primarily addressed one powerful simulation tool developed during the last decades-Large Eddy Simulation (LES), which is most suitable for unsteady three-dimensional complex turbulent flows in industry and natural environment. The main point in LES is that the large-scale motion is resolved while the small-scale motion is modeled or, in geophysical terminology, parameterized. With a view to devising a subgrid-scale(SGS) model of high quality, we have highlighted analyzing physical aspects in scale interaction and-energy transfer such as dissipation, backscatter, local and non-local interaction, anisotropy and resolution requirement. They are the factors responsible for where the advantages and disadvantages in existing SGS models come from. A case study on LES of turbulence in vegetative canopy is presented to illustrate that LES model is more based on physical arguments. Then, varieties of challenging complex turbulent flows in both industry and geophysical fields in the near future-are presented. In conclusion; we may say with confidence that new century shall see the flourish in the research of turbulence with the aid of LES combined with other approaches.

Extended self-similar scaling law of multi-scale eddy structure in wall turbulence

The longitudinal fluctuating velocity of a turbulent boundary layer was measured in a water channel at a moderate Reynolds number. The extended self-similar scaling law of structure function proposed by Benzi was verified. The longitudinal fluctuating velocity, in the turbulent boundary layer was decomposed into many multi-scale eddy structures by wavelet transform. The extended self-similar scaling law of structure function for each scale eddy velocity was investigated. The conclusions are I) The statistical properties of turbulence could be self-similar not only at high Reynolds number, but also at moderate and low Reynolds number, and they could be characterized by the same set of scaling exponents xi (1)(n) = n/3 and xi (2)(n) = n/3 of the fully developed regime. 2) The range of scales where the extended self-similarity valid is much larger than the inertial range and extends far deep into the dissipation range,vith the same set of scaling exponents. 3) The extended selfsimilarity is applicable not only for homogeneous turbulence, but also for shear turbulence such as turbulent boundary layers.

冲击载荷下饱和砂土中超孔隙水压力的建立与消散过程

对落锤冲击造成饱和砂土中的冲击压力进行了测量。通过测量得到了冲击时饱和砂土中超孔隙水压力的建立和消散过程，分析了多次冲击对超孔隙水压力建立的影响，得到了一次冲击就能使饱和砂土达到完全液化的冲击强度临界值。

Assessment of Microscale Test Methods of Peeling and Splitting along Surface of Thin-Film/Substrates

Peel test methods are assessed through being applied to a peeling analysis of the ductile film/ceramic substrate system. Through computing the fracture work of the system using the either beam bend model (BB model) or the general plane analysis model (GPA model), surprisingly, a big difference between both model results is found. Although the BB model can capture the plastic dissipation phenomenon for the ductile film case as the GPA model can, it is much sensitive to the choice of the peeling criterion parameters, and it overestimates the plastic bending effect unable to capture crack tip constraint plasticity. In view of the difficulty of measuring interfacial toughness using peel test method when film is the ductile material, a new test method, split test, is recommended and analyzed using the GPA model. The prediction is applied to a wedge-loaded experiment for Al-alloy double-cantilever beam in literature.

Numerical Investigation of Toroidal Shock Wave Focusing in a Cylindrical Chamber

In this paper, focusing of a toroidal shock wave propagating from an annular shock tube into a cylindrical chamber was investigated numerically with the dispersion controlled dissipation (DCD) scheme. The first case for an incident Mach number of 1.5 was conducted and compared with experiments for validation. Then, several cases were calculated for higher incident Mach numbers varying from 2.0 to 5.0, and complicated flow structures were observed. The numerical study was mainly focused on two aspects: focusing process and flow structures. The process, including diffraction, focusing, and reflection, is displayed to reveal the focusing mechanism, and the flow structures at different incident. Mach numbers are used to demonstrate shock reflection styles and focusing characteristics.

Toughness Of Ni/Al2O3 Interfaces As Dependent On Micron-Scale Plasticity And Atomistic-Scale Separation

Ceramic/metal interfaces were studied that fail by atomistic separation accompanied by plastic dissipation in the metal. The macroscopic toughness of the specific Ni alloy/Al2O3 interface considered is typically on the order of ten times the atomistic work of separation in mode I and even higher if combinations of mode I and mode II act on the interface. Inputs to the computational model of interface toughness are: (i) strain gradient plasticity applied to the Ni alloy with a length parameter determined by an indentation test, and (ii) a potential characterizing mixed mode separation of the interface fit to atomistic results. The roles of the several length parameters in the strain gradient plasticity are determined for indentation and crack growth. One of the parameters is shown to be of dominant importance, thus establishing that indentation can be used to measure the relevant length parameter. Recent results for separation of Ni/Al2O3 interfaces computed by atomistic methods are reviewed, including a set of results computed for mixed mode separation. An approximate potential fit to these results is characterized by the work of separation, the peak separation stress for normal separation and the traction-displacement relation in pure shearing of the interface. With these inputs, the model for steady-state crack growth is used to compute the toughness of the interface under mode I and under the full range of mode mix. The effect of interface strength and the work of separation on macroscopic toughness is computed. Fundamental implications for plasticity-enhanced toughness emerge.

Ductile To Brittle Transition In Dynamic Fracture Of Brittle Bulk Metallic Glass

We report an unusual transition from a locally ductile to a pure brittle fracture in the dynamic fracture of brittle Mg65Cu20Gd10 bulk metallic glass. The fractographic evolution from a dimple structure to a periodic corrugation pattern and then to the mirror zone along the crack propagation direction during the dynamic fracture process is discussed within the framework of the meniscus instability of the fracture process zone. This work might provide an important clue in understanding of the energy dissipation mechanism for dynamic crack propagation in brittle glassy materials. (C) 2008 American Institute of Physics.

Dynamic Fracture Instability Of Tough Bulk Metallic Glass

We report the observations of a clear fractographic evolution from vein pattern, dimple structure, and then to periodic corrugation structure, followed by microbranching pattern, along the crack propagation direction in the dynamic fracture of a tough Zr41.2Ti13.8Cu12.5Ni10Be22.5 (Vit.1) bulk metallic glass (BMGs) under high-velocity plate impact. A model based on fracture surface energy dissipation and void growth is proposed to characterize this fracture pattern transition. We find that once the dynamic crack propagation velocity reaches a critical fraction of Rayleigh wave speed, the crack instability occurs; hence, crack microbranching goes ahead. Furthermore, the correlation between the critical velocity of amorphous materials and their intrinsic strength such as Young's modulus is uncovered. The results may shed new insight into dynamic fracture instability for BMGs. (C) 2008 American Institute of Physics.

Instability of the vertically forced surface wave in a circular cylindrical container

The nonlinear amplitude equation, which was derived by Jian Yongjun employing expansion of two-time scales in inviscid fluids in a vertically oscillating circular cylindrical vessel, is modified by introducing a damping term due to the viscous dissipation of this system. Instability of the surface wave is analysed and properties of the solutions of the modified equation are determined together with phase-plane trajectories. A necessary condition of forming a stable surface wave is obtained and unstable regions are illustrated. Research results show that the stable pattern of surface wave will not lose its stability to an infinitesimal disturbance.

On dispersion-controlled principles for non-oscillatory shock-capturing schemes

The role of dispersions in the numerical solutions of hydrodynamic equation systems has been realized for long time. It is only during the last two decades that extensive studies on the dispersion-controlled dissipative (DCD) schemes were reported. The studies have demonstrated that this kind of the schemes is distinct from conventional dissipation-based schemes in which the dispersion term of the modified equation is not considered in scheme construction to avoid nonphysical oscillation occurring in shock wave simulations. The principle of the dispersion controlled aims at removing nonphysical oscillations by making use of dispersion characteristics instead of adding artificial viscosity to dissipate the oscillation as the conventional schemes do. Research progresses on the dispersion controlled principles are reviewed in this paper, including the exploration of the role of dispersions in numerical simulations, the development of the dispersion-controlled principles, efforts devoted to high-order dispersion-controlled dissipative schemes, the extension to both the finite volume and the finite element methods, scheme verification and solution validation, and comments on several aspects of the schemes from author's viewpoint.

阿尔文涨落与地球磁层亚暴

本文讨论了一种地球磁层的亚暴机制。当行星际磁场有大的南向分量时,磁层的位形可由基本闭式转变为开式。磁鞘中的阿尔文波可以携带超过10~(18)尔格/秒的能流传入磁层尾部,并将能量耗散于等离子体片中。等离子体片中的粒子被加热和加速后,注入近地空间,产生环电流和极区亚暴。计算了剪切流场中阿尔文波的传播过程,以及磁层中阿尔文波的耗散。将本文的结算与[4]中的结果合在一起,可以说明当行星际磁场转向南时,容易发生地球磁层亚暴,但这两者并非一一对应的关系,行星际磁场没有南向分量时也可以发生地球磁层亚暴。

Parameter-transformation relations for traveling-wave solutions of kdv-burgers equation

 Burgers suggested that the main properties of free-turbulence in the boundless area without basic flow might be understood with the aid of the following equation, which was much simpler than those of fluid dynamics, 更多还原

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.

Optimization of 2-stage gas gun for fusion refueling pellet injection with consideration of real effects

Analytic expression of pellet acceleration by constant base pressure with consideration of gas-wall friction, heat transfer and viscous dissipation that important for high speed injection is obtained. The process of compression stage is formulated by a set of governing equations and can be numerically integrated. Excellent confirmation with experiments is obtained and the ways to optimum match the compression stage with the launch stage are suggested.

Turbulent passive scalar field of a small prandtl number

The statistical-mechanics theory of the passive scalar field convected by turbulence, developed in an earlier paper [Phys. Fluids 28, 1299 (1985)], is extended to the case of a small molecular Prandtl number. The set of governing integral equations is solved by the equation-error method. The resultant scalar-variance spectrum for the inertial range is F(k)~x−5/3/[1+1.21x1.67(1+0.353x2.32)], where x is the wavenumber scaled by Corrsin's dissipation wavenumber. This result reduces to the − (5)/(3) law in the inertial-convective range. It also approximately reduces to the − (17)/(3) law in the inertial-diffusive range, but the proportionality constant differs from Batchelor's by a factor of 3.6.