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.

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.

The interaction of two pairs of solitary waves in a two-fluid system

This paper deals with the interaction of solitary waves in a two-fluid system which consistsof two superimposed incompressible inviscid fluids with a free surface and a horizontal rigidbottom. Under the assumption of shallow water wave, we first derive the basic equationssuitable for the model considered, a generalized form of the Boussinesq equations, then usingthe PLK method and the reductive perturbation method, obtain the second-order approximatesolution for the head-on collision between two pairs of interface and surface solitary waves,and give their maximum amplitudes during the collision and the nonuniform phase shiftsafter the collision which lead to the distortion of the wave profiles.

Variational approach to the closure problem of turbulence theory

A new method is proposed to solve the closure problem of turbulence theory and to drive the Kolmogorov law in an Eulerian framework. Instead of using complex Fourier components of velocity field as modal parameters, a complete set of independent real parameters and dynamic equations are worked out to describe the dynamic states of a turbulence. Classical statistical mechanics is used to study the statistical behavior of the turbulence. An approximate stationary solution of the Liouville equation is obtained by a perturbation method based on a Langevin-Fokker-Planck (LFP) model. The dynamic damping coefficient eta of the LFP model is treated as an optimum control parameter to minimize the error of the perturbation solution; this leads to a convergent integral equation for eta to replace the divergent response equation of Kraichnan's direct-interaction (DI) approximation, thereby solving the closure problem without appealing to a Lagrangian formulation. The Kolmogorov constant Ko is evaluated numerically, obtaining Ko = 1.2, which is compatible with the experimental data given by Gibson and Schwartz, (1963).

The asymptotic near-tip solution for mode-iii crack in steady growth in power hardening media

Two local solutions, one perpendicular and one parallel to the direction of solar gravitational field, are discussed. The influence of gravity on the gas-dynamical process driven by the piston is discussed in terms of characteristic theory, and the flow field is given quantitatively. For a typical piston trajectory similar to the one for an eruptive prominence, the velocity of the shock front which locates ahead the transient front is nearly constant or slightly accelerated, and the width of the compressed flow region may be kept nearly constant or increased linearly, depending on the velocity distribution of the piston. Based on these results, the major features of the transient may be explained. Some of the fine structure of the transient is also shown, which may be compared in detail with observations.

On the magnetospheric structure of pulsars

The magnetospheric structure of a pulsar is discussed for a non-force-free magnetic field. The local solution to the axisymmetric equations of the pulsar is obtained by the method of expanding in the polar angle. Particular attention is given to the solutions near the polar axis and the equator. Near the pulsar surface, the magnetic field energy density is found to be larger than the other energy components; the gravitational potential and the kinetic energy are relatively larger far away from the pulsar surface. It is shown that these relations influence the mass distribution in the pulsar magnetosphere. The results also show that the plasma rotation may be nonrigid and, hence, that a corotational region with a closed magnetic field may not exist.

PICES Press, Vol. 10, No. 2, July 2002

International symposium on North Pacific transitional areas [pp. 1-4] [pdf, 0.8 Mb] PICES Volunteer Observing Ship (VOS) Workshop [pp. 5-7] [pdf, 0.3 Mb] Joint meeting on Causes of marine mortality of salmon [pp. 8-9] [pdf, 0.3 Mb] The state of the western North Pacific in the second half of 2001 [pp. 10-11] [pdf, 0.5 Mb] State of the eastern North Pacific in spring 2002 [pp. 12-13] [pdf. 0.4 Mb] The status of the Bering Sea in the second half of 2001 [pp. 14-15] [pdf. 0.3 Mb] PICES Workshop on “Perturbation analysis” on subarctic Pacific gyre ecosystem models [pp. 16-17] [pdf. 0.4 Mb] Status and future plans for SOLAS-Japan [pp. 18-20] [pdf. 0.5 Mb] China-Korea Joint Ocean Research Center: A bridge across the Yellow Sea to connect Chinese and Korean oceanographic institutes and scientists [pp. 21-22] [pdf. 0.3 Mb] Persistent changes in the California Current ecosystem [pp. 23-24] [pdf. 0.2 Mb] The Hokusei Maru: 53 years of research in the Pacific [pp. 25-28] [pdf. 0.5 Mb] First meeting of the CLIVAR Pacific Panel [pp. 29-30] [pdf. 0.3 Mb] Call for contributions to the North Pacific Ecosystem Status Report [p. 31] [pdf. 0.2 Mb] PICES announcements [p. 32] [pdf. 0.2 Mb]

PICES-GLOBEC International Program on Climate Change and Carrying Capacity

Introduction [pdf, 0.27 MB] Methods [pdf, 0.15 MB] Results and discussion [pdf, 2.1 MB] Conclusions [pdf, 0.12 MB] Appendix A: Data gathering review, results and balancing [pdf, 0.3 MB] Appendix B: Data tables [pdf, 0.35 MB] Appendix C: BASS Workshop on the "Development of a conceptual model of the subarctic Pacific Basin ecosystems" [pdf, 0.16 MB] Appendix D: BASS/MODEL Workshop on "Higher trohic level modeling" [pdf, 0.24 MB] Appendix E: BASS/MODEL Workshop to review ecosystem models for the subarctic Pacific gyres [pdf, 4.39 MB] Appendix F: BASS/MODEL Workshop on "Perturbation analysis" on subarctic Pacific gyre ecosystem models using ECOPATH/ECOSIM" [pdf, 0.37 MB] Appendix G: Proposal for a BASS Workshop on "Linkages between open and coastal systems" [pdf, 0.15 MB] References [pdf, 0.14 MB] (97 page document)

Numerical Simulation of R-M Instability

In order to capture shock waves and contact discontinuities in the field and easy to program with parallel computation a new algorithm is developed to solve the N-S equations for simulation of R-M instability problems. The method with group velocity control is used to suppress numerical oscillations, and an adaptive non-uniform mesh is used to get fine resolution. Numerical results for cylindrical shock-cylindrical interface interaction with a shock Mach number Ms=1.2 and Atwood number A=0.818, 0.961, 0.980 (the interior density of the interface/outer density p(1)/p(2) = 10, 50, 100, respectively), and for the planar shock-spherical interface interaction with Ms=1.2 and p(1)/p(2) = 14.28are presented. The effect of Atwood number and multi-mode initial perturbation on the R-M instability are studied. Multi-collisions of the reflected shock with the interface is a main reason of nonlinear development of the interface instability and formation of the spike-bubble structures In simulation with double mode perturbation vortex merging and second instability are found. After second instability the small vortex structures near the interface produced. It is important factor for turbulent mixing.

Perturbational finite volume scheme for the one-dimensional Navier-Stokes equations

Starting from the second-order finite volume scheme,though numerical value perturbation of the cell facial fluxes, the perturbational finite volume (PFV) scheme of the Navier-Stokes (NS) equations for compressible flow is developed in this paper. The central PFV scheme is used to compute the one-dimensional NS equations with shock wave.Numerical results show that the PFV scheme can obtain essentially non-oscillatory solution.

About the Stabilization of a Nonlinear Perturbed Difference Equation

This paper investigates the local asymptotic stabilization of a very general class of instable autonomous nonlinear difference equations which are subject to perturbed dynamics which can have a different order than that of the nominal difference equation. In the general case, the controller consists of two combined parts, namely, the feedback nominal controller which stabilizes the nominal (i.e., perturbation-free) difference equation plus an incremental controller which completes the stabilization in the presence of perturbed or unmodeled dynamics in the uncontrolled difference equation. A stabilization variant consists of using a single controller to stabilize both the nominal difference equation and also the perturbed one under a small-type characterization of the perturbed dynamics. The study is based on Banach fixed point principle, and it is also valid with slight modification for the stabilization of unstable oscillatory solutions.

若干流体力学问题的数学性质研究

该文利用高智的扩散抛物化方程组理论及流体力学基本方程组的特征次特征理论,流体大小尺度（LSS）方程组理论以及摄动有限差分（PFD）方法,研究若干流体力学问题的数学性质.该文得到的主要结论有：1.利用湍流大小尺度（LSS）方程组推导出湍流大小尺度涡量（LSSV）方程组,并证明两个关于湍流大小尺度涡量的命题,从而得到湍流封闭大小尺度涡量（CLSSV）方程组,并对已有的近程相互作用命题进行推广.2.根据扩散抛物化方程组理论和流体力学层次结构方程组的特征和次特征方法,研究了抛物化稳定性方程组（PSE）的特征和次特征以及消除PSE的剩余椭圆特性的问题.3.利用摄动有限差分（PFD）方法得到对流扩散反应方程的变步长摄动有限差分格式,是等步长摄动有限差分格式的推广.

Is General Relativity a simplified theory?

IARD 8th Biennial Conference on Classical and Quantum Relativistic Dynamics of Particles and Fields - Galileo Galilei Inst Theoret Phys (GGI), Florence, ITALY - MAY 29-JUN 01, 2012. Edited by:Horowitz, LP

Pumped double quantum dot with spin-orbit coupling

We study driven by an external electric field quantum orbital and spin dynamics of electron in a one-dimensional double quantum dot with spin-orbit coupling. Two types of external perturbation are considered: a periodic field at the Zeeman frequency and a single half-period pulse. Spin-orbit coupling leads to a nontrivial evolution in the spin and orbital channels and to a strongly spin-dependent probability density distribution. Both the interdot tunneling and the driven motion contribute into the spin evolution. These results can be important for the design of the spin manipulation schemes in semiconductor nanostructures.

Formation Mechanism Of Lamellar Chips During Machining Of Bulk Metallic Glass

The uniqThe unique lamellar chips formed in turning–machining of a Vit 1 bulk metallic glass (BMG) are found to be due to repeated shearband formation in the primary shear zone (PSZ). A coupled thermomechanical orthogonal cutting model, taking into account force, free volume and energy balance in the PSZ, is developed to quantitatively characterize lamellar chip formation. Its onset criterion is revealed through a linear perturbation analysis. Lamellar chip formation is understood as a self-sustained limit-cycle phenomenon: there is autonomous feedback in stress, free volume and temperature in the PSZ. The underlying mechanism is the symmetry breaking of free volume flow and source, rather than thermal instability. These results are fundamentally useful for machining BMGs and even for understanding the physical nature of inhomogeneous flow in BMGs.ue lamellar chips formed in turning–machining of a Vit 1 bulk metallic glass (BMG) are found to be due to repeated shearband.