Magnetostatic relations including fluctuation fields - the local expansion

Using the approach of local expansion, we analyze the magnetostatic relations in the case of conventional turbulence. The turbulent relations are obtained consisten tly for themomentum equation and induction equation of both the average and fluctuation relations.In comparison with the magnetostatic relations as discussed usually, turbulent fluctuationfields produce forces, one of which 1/(4π)(α1×B0)×B0 may have parallel and perpendicular components in the direction of magnetic field, the other of which 1/(4π)K×B0 is introduced by the boundary value of turbulence and is perpendicular to the magnetic field. In the case of 2-dimensional configuration of magnetic field, the basic equation will be reduced into a second-order elliptic equation, which includes some linear and nonlinear terms introduced by turbulent fluctuation fields. Turbulent fields may change the configuration of magnetic field and even shear it non-uniformly. The study on the influence of turbulent fields is significant since they are observed in many astrophysical environments.

The influence of fluctuation fields on the force-free field .2. The local expansion

In Paper I (Hu, 1982), we discussed the the influence of fluctuation fields on the force-free field for the case of conventional turbulence and demonstrated the general relationships. In the present paper, by using the approach of local expansion, the equation of average force-free field is obtained as (1+b)×B 0=(#x002B;a)B 0#x002B;a (1)×B 0#x002B;K. The average coefficientsa,a (1),b, andK show the influence of the fluctuation fields in small scale on the configurations of magnetic field in large scale. As the average magnetic field is no longer parallel to the average electric current, the average configurations of force-free fields are more general and complex than the usual ones. From the view point of physics, the energy and momentum of the turbulent structures should have influence on the equilibrium of the average fields. Several examples are discussed, and they show the basic features of the fluctuation fields and the influence of fluctuation fields on the average configurations of magnetic fields. The astrophysical environments are often in the turbulent state, the results of the present paper may be applied to the turbulent plasma where the magnetic field is strong.

The magnetostatic relations including fluctuation fields

The influences of the fluctuation fields are important in many astrophysical environments as shown by the observations, and can not be neglected. On the basis of the first-order smoothing approximation, in the present paper, we demonstrate the magnetostatic equations for both the cases of the conventional turbulence aud the random waves, and discuss the consistent conditions of the equations. In the static problem, the fluctuation Lorentz force(▽×δB)×δB influences the large-scale configurations of magnetic field. To study this influence in detail is quite necessary for the explanations of the observation features, especially for the astrophysical environments where the magnetic fields, including the fluctuation fields, are the dominant factors in the equilibrium of momentum and energy.

Real-Time Measurement on Deformation Fields of Hole-Excavated Samples under Impact Tension

In this paper, the real-time deformation fields are observed in two different kinds of hole-excavated dog-bone samples loaded by an SHTB, including single hole sample and dual holes sample with the aperture size of 0.8mm. The testing system consists of a high-speed camera, a He-Ne laser, a frame grabber and a synchronization device with the controlling accuracy of I microsecond. Both the single hole expanding process and the interaction of the two holes are recorded with the time interval of 10 mu s. The observed images on the sample surface are analyzed by newly developed software based on digital correlation theory and a modified image processing method. The 2-D displacement fields in plane are obtained with a resolution of 50 mu m and an accuracy of 0.5 mu m. Experimental results obtained in this paper are proofed, by compared with FEM numerical simulations.

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<numbers RDxi between different cells in Navier-Stokes (NS) equations computations for high Reynolds number flows, an idea of solving the conservation equations for discrete cells was proposed and named the discrete fluid dynamics (DFD) algorithm. Analysis shows that the basic conservative equations for discrete cells are the Euler equations, NS- and diffusion parabolized (DP) NS equations. In this paper, a new multiscale-domain decomposition method is developed for the high Reynolds number flow. First, the whole domain is decomposed to different sub-domains with the different characteristic scales. Then the different dominant equation of all sub-domains is defined according to the diffusion parabolized (DP) theory of viscous flow. Finally these different equations are solved simultaneously in whole computational region. For numerical tests of high Reynolds numerical flows, two-dimensional supersonic flows over rearward and frontward steps as well as an interaction flow between shock wave and boundary layer were solved numerically. The pressure distributions and local coefficients of skin friction on the wall are given. The numerical results obtained by the multiscale-domain decomposition algorithm are well agreement with those by NS equations. Comparing with the usual method of solving the Navier-Stokes equations in the whole flow, under the same numerical accuracy, the present multiscale domain decomposition method decreases CPU consuming about 20% and reflects the physical mechanism of practical flow more accurately.

Numerical-Value Perturbation Method with Application of Solving Convective-Diffusion Integral Equation

The convective--diffusion equation is of primary importance in such fields as fluid dynamics and heat transfer hi the numerical methods solving the convective-diffusion equation, the finite volume method can use conveniently diversified grids (structured and unstructured grids) and is suitable for very complex geometry The disadvantage of FV methods compared to the finite difference method is that FV-methods of order higher than second are more difficult to develop in three-dimensional cases. The second-order central scheme (2cs) offers a good compromise among accuracy, simplicity and efficiency, however, it will produce oscillatory solutions when the grid Reynolds numbers are large and then very fine grids are required to obtain accurate solution. The simplest first-order upwind (IUW) scheme satisfies the convective boundedness criteria, however. Its numerical diffusion is large. The power-law scheme, QMCK and second-order upwind (2UW) schemes are also often used in some commercial codes. Their numerical accurate are roughly consistent with that of ZCS. Therefore, it is meaningful to offer higher-accurate three point FV scheme. In this paper, the numerical-value perturbational method suggested by Zhi Gao is used to develop an upwind and mixed FV scheme using any higher-order interpolation and second-order integration approximations, which is called perturbational finite volume (PFV) scheme. The PFV scheme uses the least nodes similar to the standard three-point schemes, namely, the number of the nodes needed equals to unity plus the face-number of the control volume. For instanc6, in the two-dimensional (2-D) case, only four nodes for the triangle grids and five nodes for the Cartesian grids are utilized, respectively. The PFV scheme is applied on a number of 1-D problems, 2~Dand 3-D flow model equations. Comparing with other standard three-point schemes, The PFV scheme has much smaller numerical diffusion than the first-order upwind (IUW) scheme, its numerical accuracy are also higher than the second-order central scheme (2CS), the power-law scheme (PLS), the QUICK scheme and the second-order upwind(ZUW) scheme.

Characteristics of flow fields induced by interfacial waves in two-layer fluid

When designing deep ocean structures, it is necessary to estimate the effects of internal waves on the platform and auxiliary parts such as tension leg, riser and mooring lines. Up to now, only a few studies are concerned with the internal wave velocity fields. By using the most representative two-layer model, we have analyzed the behavior of velocity field induced by interfacial wave in the present paper. We find that there may exist velocity shear of fluid particles in the upper and lower layers so that any structures in the ocean are subjected to shear force nearby the interface. In the meantime, the magnitude of velocity for long internal wave appears spatially uniform in the respective layer although they still decay exponentially. Finally, the temporal variation for Stokes and solitary waves are shown to be of periodical and pulse type.

Floating zone convection

The microgravity research, as a branch of the advanced sciences and a spe- cialized field of high technology, has been made in China since the late 1980's. The research group investigating microgravity fluid physics consisted of our col- leagues and the authors in the Institute of Mechanics of the Chinese Academy of Sciences (CAS), and we pay special attention to the floating zone convection as our first research priority. Now, the research group has expanded and is a part of the National Microgravity Laboratory of the CAS, and the research fields have been extended to include more subjects related to microgravity science. Howev- er, the floating zone convection is still an important topic that greatly holds our research interests.

Smaller Deborah number inducing more serrated plastic flow of metallic glass

Spherical nano-indentations of Cu46Zr54 bulk metallic glass (BMG) model systems were performed using molecular dynamics (MD) computer simulations, focusing specifically on the physical origin of serrated plastic flow. The results demonstrate that there is a direct correlation between macroscopic flow serration and underlying irreversible rearrangement of atoms, which is strongly dependent on the loading (strain) rate and the temperature. The serrated plastic flow is, therefore, determined by the magnitude of such irreversible rearrangement that is inhomogeneous temporally. A dimensionless Deborah number is introduced to characterize the effects of strain rate and temperature on serrations. Our simulations are shown to compare favorably with the available experimental observations.

A universal Biot number determining the susceptibility of ceramics to quenching

A universal Biot number of ceramics, which not only determines the susceptibility of the ceramics to quenching but also indicates the duration that the ceramics fail during thermal shock, is theoretically obtained. The present analysis shows that the thermal shock failure of the ceramics with a Biot number greater than this universal value is a very rapid process that just occurs in the initial regime of the heat conduction of the ceramics. This universal Biot number provides a guide to the selection of the ceramics applying to the thermostructural engineering including thermal shock.

Large-eddy simulation of turbulent collision of heavy particles in isotropic turbulence

The small-scale motions relevant to the collision of heavy particles represent a general challenge to the conventional large-eddy simulation (LES) of turbulent particle-laden flows. As a first step toward addressing this challenge, we examine the capability of the LES method with an eddy viscosity subgrid scale (SGS) model to predict the collision-related statistics such as the particle radial distribution function at contact, the radial relative velocity at contact, and the collision rate for a wide range of particle Stokes numbers. Data from direct numerical simulation (DNS) are used as a benchmark to evaluate the LES using both a priori and a posteriori tests. It is shown that, without the SGS motions, LES cannot accurately predict the particle-pair statistics for heavy particles with small and intermediate Stokes numbers, and a large relative error in collision rate up to 60% may arise when the particle Stokes number is near St_K=0.5. The errors from the filtering operation and the SGS model are evaluated separately using the filtered-DNS (FDNS) and LES flow fields. The errors increase with the filter width and have nonmonotonic variations with the particle Stokes numbers. It is concluded that the error due to filtering dominates the overall error in LES for most particle Stokes numbers. It is found that the overall collision rate can be reasonably predicted by both FDNS and LES for St_K>3. Our analysis suggests that, for St_K<3, a particle SGS model must include the effects of SGS motions on the turbulent collision of heavy particles. The spectral analysis of the concentration fields of the particles with different Stokes numbers further demonstrates the important effects of the small-scale motions on the preferential concentration of the particles with small Stokes numbers.

Beam splitter entangler for light fields

We propose an experimentally feasible scheme to generate various types of entangled states of light fields by using beam splitters and single-photon detectors. Two beams of light fields are incident on two beam splitters respectively with each beam being asymmetrically split into two parts in which one part is supposed to be so weak that it contains at most one photon. We let the two weak output modes interfere at a third beam splitter. A conditional joint measurement on both weak output modes may result in an entanglement between the other two output modes. The conditions for the maximal entanglement are discussed based on the concurrence. Several specific examples are also examined.

Topology of toroidal helical fields in non-circular cross-sectional tokamaks

The ordinary differential magnetic field line equations are solved numerically; the tokamak magnetic structure is studied on Hefei Tokamak-7 Upgrade (HT-7U) when the equilibrium field with a monotonic q-profile is perturbed by a helical magnetic field. We find that a single mode (m, n) helical perturbation can cause the formation of islands on rational surfaces with q = m/n and q = (m +/- 1, +/- 2, +/- 3,...)/n due to the toroidicity and plasma shape (i.e. elongation and triangularity), while there are many undestroyed magnetic surfaces called Kolmogorov-Arnold-Moser (KAM) barriers on irrational surfaces. The islands on the same rational surface do not have the same size. When the ratio between the perturbing magnetic field B-r(r) and the toroidal magnetic field amplitude B(phi)0 is large enough, the magnetic island chains on different rational surfaces will overlap and chaotic orbits appear in the overlapping area, and the magnetic field becomes stochastic. It is remarkable that the stochastic layer appears first in the plasma edge region.

Ion acceleration in laser-excited plasma wake fields

The dynamics of the plasma ions in the wake fields of short, ultraintense laser pulses in underdense plasmas are investigated analytically and numerically. Owing to the large ion-to-electron mass ratio, the motion of plasma ions in-such wake fields has often been assumed to be neglectable. It is shown that when the laser intensity exceeds 10(20) W/cm(2), the ion motion can no longer be ignored. In this case, ion momentum peaks appear behind the laser pulse, which correspond with the ion density peaks. The laser-excited wake field appears to be effective for ion acceleration, in particular to ions with high-charge numbers. The dependence of ion acceleration on the laser intensity, pulse width, and background plasma density is discussed. (c) 2006 Optical Society of America.