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.

Prediction of mixed-mode cracking direction in random, short-fibre composite-materials

The prediction of cracking direction in composite materials is of significance to the design of composite structures. This paper presents several methods for predicting the cracking direction in the double grooved tension-shear specimen which gives mixed-mode cracking. Five different criteria are used in this analysis: two of them have been used by other investigators and the others are proposed by the present authors. The strain energy density criterion proposed by G.C. Sih is modified to take account of the influence of the anisotropy of the strength on the direction of crack. The two failure criteria of Tsai-Hill and Norris are extended to predict the crack orientation. The stress distributions in the near-notch zone are calculated by using the 8-node quadrilateral isoparametric finite element method. The predictions of all the criteria except one are in good agreement with the experimental measurement. In addition, on the basis of the FEM results, the size of the zone in which the singular term is dominant is estimated.

A Hybrid Random Number Generator Using Single Electron Tunneling Junctions and MOS Transistors

This paper proposes a novel single electron random number generator (RNG). The generator consists of multiple tunneling junctions (MTJ) and a hybrid single electron transistor (SET)/MOS output circuit. It is an oscillator-based RNG. MTJ is used to implement a high-frequency oscillator,which uses the inherent physical randomness in tunneling events of the MTJ to achieve large frequency drift. The hybrid SET and MOS output circuit is used to amplify and buffer the output signal of the MTJ oscillator. The RNG circuit generates high-quality random digital sequences with a simple structure. The operation speed of this circuit is as high as 1GHz. The circuit also has good driven capability and low power dissipation. This novel random number generator is a promising device for future cryptographic systems and communication applications.

Second-order solutions for random interfacial waves under steady uniform currents

In the present research, the study of Song (2004) for random interfacial waves in two-layer fluid is extended to the case of fluids moving at different steady uniform speeds. The equations describing the random displacements of the density interface and the associated velocity potentials in two-layer fluid are solved to the second order, and the wave-wave interactions of the wave components and the interactions between the waves and currents are described. As expected, the extended solutions include those obtained by Song (2004) as one special case where the steady uniform currents of the two fluids are taken as zero, and the solutions reduce to those derived by Sharma and Dean (1979) for random surface waves if the density of the upper fluid and the current of the lower fluid are both taken as zero.

LEVEL SET输运方程的求解方法及其对气-液两相流运动界面数值模拟的影响

用数值模拟方法来研究气-液两相流动与传热现象是当今多相流领域的一个热门课题.由于两相流固有的复杂性,气-液两相流界面迁移现象的数值模拟一直是两相流研究中的一大难点.本文介绍了捕捉气-液两相流相界面运动的水平集方法(Level Set)及其研究进展,介绍了求解Level Set输运方程的3种方法,即一般差分格式、Superbee-TVD格式和Runge-Kutta法-5阶WENO组合格式.结合主流场的求解,分别用这3种方法对4种典型相界面在5种流场中的迁移特性进行了模拟计算,并对计算结果进行了比较和分析.结果表明,Runge-Kutta法-5阶WENO组合格式求解Level Set输运方程的效果最好,在以后的计算中将主要采用这种组合格式来进行气-液相界面输运方程的求解.

Modification of Bertrand's theorem and extended Runge-Lenz vector

It is shown that for a particle with suitable angular moments in the screened Coulomb potential or isotropic harmonic potential, there still exist closed orbits rather than ellipse, characterized by the conserved aphelion and perihelion vectors, i.e. extended Runge-Lenz vector, which implies a higher dynamical symmetry than the geometrical symmetry O-3. The closeness of a planar orbit implies the radial and angular motional frequencies are commensurable.

Mapping Closure Approximation for Reactive Scalars in Random Flows

The Mapping Closure Approximation (MCA) approach is developed to describe the statistics of both conserved and reactive scalars in random flows. The statistics include Probability Density Function (PDF), Conditional Dissipation Rate (CDR) and Conditional Laplacian (CL). The statistical quantities are calculated using the MCA and compared with the results of the Direct Numerical Simulation (DNS). The results obtained from the MCA are in agreement with those from the DNS. It is shown that the MCA approach can predict the statistics of reactive scalars in random flows.

Optimal bounded control of first-passage failure of quasi-integrable Hamiltonian systems with wide-band random excitation

The optimal bounded control of quasi-integrable Hamiltonian systems with wide-band random excitation for minimizing their first-passage failure is investigated. First, a stochastic averaging method for multi-degrees-of-freedom (MDOF) strongly nonlinear quasi-integrable Hamiltonian systems with wide-band stationary random excitations using generalized harmonic functions is proposed. Then, the dynamical programming equations and their associated boundary and final time conditions for the control problems of maximizinig reliability and maximizing mean first-passage time are formulated based on the averaged It$\ddot{\rm o}$ equations by applying the dynamical programming principle. The optimal control law is derived from the dynamical programming equations and control constraints. The relationship between the dynamical programming equations and the backward Kolmogorov equation for the conditional reliability function and the Pontryagin equation for the conditional mean first-passage time of optimally controlled system is discussed. Finally, the conditional reliability function, the conditional probability density and mean of first-passage time of an optimally controlled system are obtained by solving the backward Kolmogorov equation and Pontryagin equation. The application of the proposed procedure and effectiveness of control strategy are illustrated with an example.

Study on current-random wave forces acting on a framework

The forces of random wave plus current acting on a simplified offshore platform (jacket) model have been studied numerically and experimentally. The numerical results are in good agreement with experiments. The mean force can be approximated as a function of equivalent velocity parameter and the root-mean-square force as a function of equivalent significant wave height parameter.

Flows with Inertia in a Three-Dimensional Random Fiber Network

A fiber web is modeled as a three-dimensional random cylindrical fiber network. Nonlinear behavior of fluid flowing through the fiber network is numerically simulated by using the lattice Boltzmann (LB) method. A nonlinear relationship between the friction factor and the modified Reynolds number is clearly observed and analyzed by using the Fochheimer equation, which includes the quadratic term of velocity. We obtain a transition from linear to nonlinear region when the Reynolds numbers are sufficiently high, reflecting the inertial effect of the flows. The simulated permeability of such fiber network has relatively good agreement with the experimental results and finite element simulations.

Influence of stochastic mesoscopic structure on macroscopic mechanical behavior of brittle material

A newly developed numerical code, MFPA(2D) (Material Failure Process Analysis), is applied to study the influence of stochastic mesoscopic structure on macroscopic mechanical behavior of rock-like materials. A set of uniaxial compression tests has been numerically studied with numerical specimens containing pre-existing crack-like flaw. The numerical results reveal the influence of random mesoscopic structure on failure process of brittle material, which indicates that the variation of failure mode is strongly sensitive to the local disorder feature of the specimen. And the patterns of the crack evolution in the specimens are very different from each other due to the random mesoscopic structure in material. The results give a good explanation for various kinds of fracture modes and peak strength variation observed in laboratory studies with specimens made from the same rock block being statistically homogenous in macro scale. In addition, the evolution of crack is more complicated in heterogeneous cases than in homogeneous cases.

Direct Numerical Simulation Of Hypersonic Boundary-Layer Transition Over A Blunt Cone

Direct numerical simulation of transition How over a blunt cone with a freestream Mach number of 6, Reynolds number of 10,000 based on the nose radius, and a 1-deg angle of attack is performed by using a seventh-order weighted essentially nonoscillatory scheme for the convection terms of the Navier-Stokes equations, together with an eighth-order central finite difference scheme for the viscous terms. The wall blow-and-suction perturbations, including random perturbation and multifrequency perturbation, are used to trigger the transition. The maximum amplitude of the wall-normal velocity disturbance is set to 1% of the freestream velocity. The obtained transition locations on the cone surface agree well with each other far both cases. Transition onset is located at about 500 times the nose radius in the leeward section and 750 times the nose radius in the windward section. The frequency spectrum of velocity and pressure fluctuations at different streamwise locations are analyzed and compared with the linear stability theory. The second-mode disturbance wave is deemed to be the dominating disturbance because the growth rate of the second mode is much higher than the first mode. The reason why transition in the leeward section occurs earlier than that in the windward section is analyzed. It is not because of higher local growth rate of disturbance waves in the leeward section, but because the growth start location of the dominating second-mode wave in the leeward section is much earlier than that in the windward section.

Computer Simulation Of Random Sphere Packing In An Arbitrarily Shaped Container

Most simulations of random sphere packing concern a cubic or cylindric container with periodic boundary, containers of other shapes are rarely studied. In this paper, a new relaxation algorithm with pre-expanding procedure for random sphere packing in an arbitrarily shaped container is presented. Boundaries of the container are simulated by overlapping spheres which covers the boundary surface of the container. We find 0.4 similar to 0.6 of the overlap rate is a proper value for boundary spheres. The algorithm begins with a random distribution of small internal spheres. Then the expansion and relaxation procedures are performed alternately to increase the packing density. The pre-expanding procedure stops when the packing density of internal spheres reaches a preset value. Following the pre-expanding procedure, the relaxation and shrinking iterations are carried out alternately to reduce the overlaps of internal spheres. The pre-expanding procedure avoids the overflow problem and gives a uniform distribution of initial spheres. Efficiency of the algorithm is increased with the cubic cell background system and double link data structure. Examples show the packing results agree well with both computational and experimental results. Packing density about 0.63 is obtained by the algorithm for random sphere packing in containers of various shapes.

Extended Stoney'S Formula For A Film-Substrate Bilayer With The Effect Of Interfacial Slip

The curvature-stress relation is studied for a film-substrate bilayer with the effect of interfacial slip and compared with that of an ideal interface without interfacial slip. The interfacial slip together with the dimensions, elastic and interfacial properties of the film and substrate layers can cause a significant deviation of curvature-stress relation from that with an ideal interface. The interfacial slip also results in the so-called free edge effect that the stress, constraint force, and curvature vary dramatically around the free edges. The constant curvature as predicted by Stoney's formula and the Timoshenko model of an ideal interface is no longer valid for a bilayer with a nonideal interface. The models with the assumption of an ideal interface can also lead to an erroneous evaluation on the true stress state inside a bilayer with a nonideal interface. The extended Stoney's formula incorporating the effects of both the layer dimensions and interfacial slip is presented.

First-passage failure of quasi-integrable Hamiltonian systems

The first-passage failure of quasi-integrable Hamiltonian si-stems (multidegree-of-freedom integrable Hamiltonian systems subject to light dampings and weakly random excitations) is investigated. The motion equations of such a system are first reduced to a set of averaged Ito stochastic differential equations by using the stochastic averaging method for quasi-integrable Hamiltonian systems. Then, a backward Kolmogorov equation governing the conditional reliability function and a set of generalized Pontryagin equations governing the conditional moments of first-passage time are established. Finally, the conditional reliability function, and the conditional probability density and moments of first-passage time are obtained by solving these equations with suitable initial and boundary conditions. Two examples are given to illustrate the proposed procedure and the results from digital simulation are obtained to verify the effectiveness of the procedure.