A closure theory of intermittency of turbulence

The variational approach to the closure problem of turbulence theory, proposed in an earlier article [Phys. Fluids 26, 2098 (1983); 27, 2229 (1984)], is extended to evaluate the flatness factor, which indicates the degree of intermittency of turbulence. Since the flatness factor is related to the fourth moment of a turbulent velocity field, the corresponding higher-order terms in the perturbation solution of the Liouville equation have to be considered. Most closure methods discard these higher-order terms and fail to explain the intermittency phenomenon. The computed flatness factor of the idealized model of infinite isotropic turbulence ranges from 9 to 15 and has the same order of magnitude as the experimental data of real turbulent flows. The intermittency phenomenon does not necessarily negate the Kolmogorov k−5/3 inertial range spectrum. The Kolmogorov k−5/3 law and the high degree of intermittency can coexist as two consistent consequences of the closure theory of turbulence. The Kolmogorov 1941 theory [J. Fluid Mech. 62, 305 (1974)] cannot be disqualified merely because the energy dissipation rate fluctuates.

太阳风动量涨落激发磁层亚暴的机制

本文将太阳风涨落传输能量产生磁层亚暴的机制推广到无碰撞等离子体过程。太阳风的涨落在磁层顶激发压缩阿尔文波,并在磁尾的无碰撞等离子体中传播。尾瓣中满足条件β<<1,而等离子体片中β≥1,其中β为等离子体压力与磁压之比。这样,快磁声波在尾瓣中几乎不衰减,而在等离子体片中很快衰减,将波动能量耗散在等离子体片中使等离子体加热或者粒子加速。这种机制还表明,磁尾等离子体片中的高能粒子可以由太阳风涨落动能耗散而被加速,不一定是直接源于太阳。

A passive scalar field convected by turbulence

Classical statistical mechanics is applied to the study of a passive scalar field convected by isotropic turbulence. A complete set of independent real parameters and dynamic equations are worked out to describe the dynamic state of the passive scalar field. The corresponding Liouville equation is solved by a perturbation method based upon a Langevin–Fokker–Planck model. The closure problem is treated by a variational approach reported in earlier papers. Two integral equations are obtained for two unknown functions: the scalar variance spectrum F(k) and the effective damping coefficient (k). The appearance of the energy spectrum of the velocity field in the two integral equations represents the coupling of the scalar field with the velocity field. As an application of the theory, the two integral equations are solved to derive the inertial-convective-range spectrum, obtaining F(k)=0.61 −1/3 k−5/3. Here is the dissipation rate of the scalar variance and is the dissipation rate of the energy of the velocity field. This theoretical value of the scalar Kolmogorov constant, 0.61, is in good agreement with experiments.

A Multiscale-Domain Decomposition Method for High Reynolds Number Flows

The high Reynolds number flow contains a wide range of length and time scales, and the flow domain can be divided into several sub-domains with different characteristic scales. In some sub-domains, the viscosity dissipation scale can only be considered in a certain direction; in some sub-domains, the viscosity dissipation scales need to be considered in all directions; in some sub-domains, the viscosity dissipation scales are unnecessary to be considered at all. For laminar boundary layer region, the characteristic length scales in the streamwise and normal directions are L and L Re-1/ 2 , respectively. The characteristic length scale and the velocity scale in the outer region of the boundary layer are L and U, respectively. In the neighborhood region of the separated point, the length scale l<

A New Model for Sediment Transport

The problem of predicting sediment transportation by water waves is treated analytically with the rate of wave energy dissipation or wave damping. With resorting to the theory of shallow water waves and the basis of Yamamoto’s Coulomb-damped poroelastic model, the Boussinesq-type equation has been derived over a variation depth bed. For convenience Cnoidal wave is just discussed, The Cnoidal wave with complex wave length and wave velocity, which are as a function of wave frequency, water depth, permeability, Poisson’s ratio and complex elastic moduli of bed soil, is applied to analyse the rate of sediment transportation. Considering the sediment transportation depended on the shear stress near-bed or the horizontal velocity, the conclusion of Yamamoto’s experiment in clay bed has been extended to general situation. It could be figured out that the model should provide a method to avoid the undistinguishable factors during sediment transport processes and relate mass transport with the sediment peculiarities.

“Deborah Numbers”, Coupling Multiple Space and Time Scales and Governing Damage Evolution to Failure

Two different spatial levels are involved concerning damage accumulation to eventual failure. nucleation and growth rates of microdamage nN* and V*. It is found that the trans-scale length ratio c*/L does not directly affect the process. Instead, two independent dimensionless numbers: the trans-scale one * * ( V*)including the * **5 * N c V including mesoscopic parameters only, play the key role in the process of damage accumulation to failure. The above implies that there are three time scales involved in the process: the macroscopic imposed time scale tim = /a and two meso-scopic time scales, nucleation and growth of damage, (* *4) N N t =1 n c and tV=c*/V*. Clearly, the dimensionless number De*=tV/tim refers to the ratio of microdamage growth time scale over the macroscopically imposed time scale. So, analogous to the definition of Deborah number as the ratio of relaxation time over external one in rheology. Let De be the imposed Deborah number while De represents the competition and coupling between the microdamage growth and the macroscopically imposed wave loading. In stress-wave induced tensile failure (spallation) De* < 1, this means that microdamage has enough time to grow during the macroscopic wave loading. Thus, the microdamage growth appears to be the predominate mechanism governing the failure. Moreover, the dimensionless number D* = tV/tN characterizes the ratio of two intrinsic mesoscopic time scales: growth over nucleation. Similarly let D be the “intrinsic Deborah number”. Both time scales are relevant to intrinsic relaxation rather than imposed one. Furthermore, the intrinsic Deborah number D* implies a certain characteristic damage. In particular, it is derived that D* is a proper indicator of macroscopic critical damage to damage localization, like D* ∼ (10–3~10–2) in spallation. More importantly, we found that this small intrinsic Deborah number D* indicates the energy partition of microdamage dissipation over bulk plastic work. This explains why spallation can not be formulated by macroscopic energy criterion and must be treated by multi-scale analysis.

Subsonic interface crack with crack face contact

A steady-state subsonic interface crack propagating between an elastic solid and a rigid substrate with crack face contact is studied. Two cases with respective to the contact length are considered, i.e., semi-infinite and finite crack face contact. Different from a stationary or an open subsonic interface crack, stress singularity at the crack tip in the present paper is found to be non-oscillatory. Furthermore, in the semi-infinite contact case, the singularity of the stress field near the crack tip is less than 1/2. In the finite contact case, no singularity exists near the crack tip, but less than 1/2 singularity does at the end of the contact zone. In both cases, the singularity depends on the linear contact coefficient and the crack speed. Asymptotic solutions near the crack tip are given and analyzed. In order to satisfy the contact conditions, reasonable region of the linear contact coefficient is found. In addition, the solution predicts a non-zero-energy dissipation rate due to crack face contact.

Adiabatic shear localization: occurrence, theories, and applications

Nanoscale periodic corrugation to dimple transition due to "beat" in a bulk metallic glass

We report an intriguing observation that the interaction of brittle nanoscale periodic corrugations (NPCs) can lead to the formation of ductile dimples on the dynamic fracture surface of a tough Vit 1 bulk metallic glass (BMG) under high-velocity plate impact. A “beat” phenomenon due to superposition of simple harmonic vibrations, approximately characterizing NPCs, is proposed to explain this unusual brittle-to-ductile transition. The present results agree well with our previously revealed energy dissipation mechanism in the fracture of BMGs.

Work and energy equations and the principle of generalized effective stress for unsaturated soils

The effective stress principle has been efficiently applied to saturated soils in the soil mechanics and geotechnical engineering practice; however, its applicability to unsaturated soils is still under debate. The appropriate selection of stress state variables is essential for the construction of constitutive models for unsaturated soils. Owing to the complexity of unsaturated soils, it is difficult to determine the deformation and strength behaviors of unsaturated soils uniquely with the previous single-effective-stress variable theory and two-effective-stress-variable theory in all the situations. In this paper, based on the porous media theory, the specific expression of work is proposed, and the effective stress of unsaturated soils conjugated with the displacement of the soil skeleton is further derived. In the derived work and energy balance equations, the energy dissipation in unsaturated soils is taken into account. According to the derived work and energy balance equations, all of the three generalized stresses and the conjugated strains have effects on the deformation of unsaturated soils. For considering these effects, a principle of generalized effective stress to describe the behaviors of unsaturated soils is proposed. The proposed principle of generalized effective stress may reduce to the previous effective stress theory of single-stress variable or the two-stress variables under certain conditions. This principle provides a helpful reference for the development of constitutive models for unsaturated soils.

MODELING EVOLUTION OF INTERNAL WAVES IN SOUTH CHINA SEA

Internal waves are an important factor in the design of drill operations and production in deep water, because the waves have very large amplitude and may induce large horizontal velocity. How the internal waves occur and propagate over benthal terrain is of great concern for ocean engineers. In the present paper, we have formulated a mathematical model of internal wave propagation in a two-layer deep water, which involves the effects of friction, dissipation and shoaling, and is capable of manifesting the variation of the amplitude and the velocity pattern. After calibration by field data measured at the Continental Slope in the Northern South China Sea, we have applied the model to the South China Sea, investigating the westward propagation of internal waves from the Luzon Strait, where internal waves originate due to the interaction of benthal ridge and tides. We find that the internal wave induced velocity profile is obviously characterized by the opposite flow below and above the pycnocline, which results in a strong shear, threatening safety of ocean structures, such as mooring system of oil platform, risers, etc. When internal waves propagate westwards, the amplitude attenuates due to the effects of friction and dissipation. The preliminary results show that the amplitude is likely to become half of its initial value at Luzon Strait when the internal waves propagate about 400 kilometers westwards.

Intense laser beam guiding in self-induced electron cavitation channel in underdense plasmas

In underdense plasmas, the transverse ponderomotive force of an intense laser beam with Gaussian transverse profile expels electrons radially, and it can lead to an electron cavitation. An improved cavitation model with charge conservation constraint is applied to the determination of the width of the electron cavity. The envelope equation for laser spot size derived by using source-dependent expansion method is extended to including the electron cavity. The condition for self-guiding is given and illuminated by an effective potential for the laser spot size. The effects of the laser power, plasma density and energy dissipation on the self-guiding condition are discussed.

Influence of temperature and salinity fluctuations on propagation behaviour of partially coherent beams in oceanic turbulence

A theoretical study of the behaviour of partially coherent beams propagating through oceanic turbulence has been performed. Based on the previously developed knowledge of beam spreading of a partially coherent beam in the atmosphere and the spatial power spectrum of the refractive index of ocean water, we study the normalized root-mean-square width of a partially coherent beam on propagation through oceanic turbulence and its turbulence distance which may be a measure of turbulence resistance. Our analysis indicates that the behaviour of partially coherent beams on propagation may be described by the rate of dissipation of the mean-squared temperature chi(T) and that of salinity chi(S). In terms of a quantity w that defines the contributions of the temperature and salinity distributions to the distribution of the refractive index, chi(S) could be written as a function of chi(T) and w. Therefore, the behaviour of partially coherent beams on propagation can be characterized only by chi(T) for a given w. The results are shown for curved surfaces, from which one can see that partially coherent beams exhibit robust turbulence resistance when the water volume has a smaller chi(T).

Wavefront measurement of double-clad pulsed ytterbium fiber laser

Because of high efficiencies, compact structure, and excellent heat dissipation, high-power fiber lasers are extremely useful for applications such as cutting, welding, precision drilling, trimming, sensing, optical transmitter, material processing, micromachining, and so on. However, the wavefront of the double clad fiber laser doped with ytterbium is still unknown. In this paper, wavefront of a fiber laser is measured and the traditional Hartmann-shack wavefront sensing method is adopted. We measured a double clad fiber laser doped with ytterbium which produces pulse wave output at infrared wavelength. The wavefront shape and contour are reconstructed and the result shows that wavefront is slightly focused and not an ideal plane wavefront. Wavefront measurement of fiber laser will be useful to improving the lasers' performance and developing the coherent technique for its applications.

千瓦级双包层光纤激光器冷却方案设计理论和实验研究

尽管双包层光纤激光器的散热性能好于传统的固体激光器的散热性能,光纤激光器中的热沉积仍然是限制提高其输出功率的重要因素.以双端抽运的400W双包层光纤激光器为实例,定量分析了光纤内的热沉积分布.根据所建立的散热模型,为了确保千瓦级双包层光纤激光器安全稳定的运行,抽运端附近的对流换热系数应大于2.8×10-2W·cm-2K-1.据此设计出高功率双包层光纤激光器抽运端冷却装置并成功应用在激光系统中,获得了千瓦级的激光输出.