An asymptotic solution of hydrodynamic equations at the mid-plane of a plume in the earths mantle

From the partial differential equations of hydrodynamics governing the movements in the Earth's mantle of a Newtonian fluid with a pressure- and temperature-dependent viscosity, considering the bilateral symmetry of velocity and temperature distributions at the mid-plane of the plume, an analytical solution of the governing equations near the mid-plane of the plume was found by the method of asymptotic analysis. The vertical distribution of the upward velocity, viscosity and temperature at the mid-plane, and the temperature excess at the centre of the plume above the ambient mantle temperature were then calculated for two sets of Newtonian rheological parameters. The results obtained show that the temperature at the mid-plane and the temperature excess are nearly independent of the rheological parameters. The upward velocity at the mid-plane, however, is strongly dependent on the rheological parameters.

The thermal and mechanical structure of a two-dimensional plume in the earths mantle

On the condition that the distribution of velocity and temperature at the mid-plane of a mantle plume has been obtained (pages 213–218, this issue), the problem of determining the lateral structure of the plume at a given depth is reduced to solving an eigenvalue problem of a set of ordinary differential equations with five unknown functions, with an eigenvalue being related to the thermal thickness of the plume at this depth. The lateral profiles of upward velocity, temperature and viscosity in the plume and the thickness of the plume at various depths are calculated for two sets of Newtonian rheological parameters. The calculations show that the precondition for the existence of the plume, δT/L 1 (L = the height of the plume, δT = lateral distance from the mid-plane), can be satisfied, except for the starting region of the plume or near the base of the lithosphere. At the lateral distance, δT, the upward velocity decreases to 0.1 – 50% of its maximum value at different depths. It is believed that this model may provide an approach for a quantitative description of the detailed structure of a mantle plume.

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).

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<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.

Modified Simplified Rate Equation Model of Flowing Chemical Oxygen-Iodine Laser and its Application

A modified simplified rate equation (RE) model of flowing chemical oxygen-iodine laser (COIL), which is adapted to both the condition of homogeneous broadening and inhomogeneous broadening being of importance and the condition of inhomogeneous broadening being predominant, is presented for performance analyses of a COIL. By using the Voigt profile function and the gain-equal-loss approximation, a gain expression has been deduced from the rate equations of upper and lower level laser species. This gain expression is adapted to the conditions of very low gas pressure up to quite high pressure and can deal with the condition of lasing frequency being not equal to the central one of spectral profile. The expressions of output power and extraction efficiency in a flowing COIL can be obtained by solving the coupling equations of the deduced gain expression and the energy equation which expresses the complete transformation of the energy stored in singlet delta state oxygen into laser energy. By using these expressions, the RotoCOIL experiment is simulated, and obtained results agree well with experiment data. Effects of various adjustable parameters on the performances of COIL are also presented.

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.

Basic equations for ferroelectrics

A set of new formula of energy functions for ferroelectrics was proposed, and then the new basic equations were derived in this paper. The finite element formulation based on the new basic equations was improved to avoid the equivalent nodal load produced by remnant polarization. With regard to the fundamentals of mathematics and physics, the new energy functions and basic equations are reasonable for the material element of ferroelectrics in finite element analysis.

Occurrence of spanning of a submarine pipeline with initial embedment

For better understanding the mechanism of the occurrence of pipeline span for a pipeline with initial embedment, physical and numerical methods are adopted in this study. Experimental observations show that there often exist three characteristic phases in the process of the partially embedded pipeline being suspended: (a) local scour around pipe; (b) onset of soil erosion beneath pipe; and (c) complete suspension of pipe. The effects of local scour on the onset of soil erosion beneath the pipe are much less than those of soil seepage failure induced by the pressure drop. Based on the above observations and analyses, the mechanism of the occurrence of pipeline spanning is analyzed numerically in view of soil seepage failure. In the numerical analyses, the current-induced pressure along the soil surface in the vicinity of the pipe (i.e. the pressure drop) is firstly obtained by solving the N-S equations, thereafter the seepage flow in the soil is calculated with the obtained pressure drop as the boundary conditions along the soil surface. Numerical results indicate that the seepage failure (or piping) may occur at the exit of the seepage path when the pressure gradient gets larger than the critical value. The numerical treatment provides a practical tool for evaluating the potentials for the occurrence of pipe span due to the soil seepage failure.

Energy Functions And Basic Equations For Ferroelectrics

Energy functions (or characteristic functions) and basic equations for ferroelectrics in use today are given by those for ordinary dielectrics in the physical and mechanical communications. Based on these basic equations and energy functions, the finite element computation of the nonlinear behavior of the ferroelectrics has been carried out by several research groups. However, it is difficult to process the finite element computation further after domain switching, and the computation results are remarkably deviating from the experimental results. For the crack problem, the iterative solution of the finite element calculation could not converge and the solutions for fields near the crack tip oscillate. In order to finish the calculation smoothly, the finite element formulation should be modified to neglect the equivalent nodal load produced by spontaneous polarization gradient. Meanwhile, certain energy functions for ferroelectrics in use today are not compatible with the constitutive equations of ferroelectrics and need to be modified. This paper proposes a set of new formulae of the energy functions for ferroelectrics. With regard to the new formulae of the energy functions, the new basic equations for ferroelectrics are derived and can reasonably explain the question in the current finite element analysis for ferroelectrics.

Numerical Investigation On The Flowfield Of "Swallowtail" Cavity For Supersonic Mixing Enhancement

A "swallowtail" cavity for the supersonic combustor was proposed to serve as an efficient flame holder for scramjets by enhancing the mass exchange between the cavity and the main flow. A numerical study on the "swallowtail" cavity was conducted by solving the three-dimensional Reynolds-averaged Navier-Stokes equations implemented with a k-epsilon turbulence model in a multi-block mesh. Turbulence model and numerical algorithms were validated first, and then test cases were calculated to investigate into the mechanism of cavity flows. Numerical results demonstrated that the certain mass in the supersonic main flow was sucked into the cavity and moved spirally toward the combustor walls. After that, the flow went out of the cavity at its lateral end, and finally was efficiently mixed with the main flow. The comparison between the "swallowtail" cavity and the conventional one showed that the mass exchanged between the cavity and the main flow was enhanced by the lateral flow that was induced due to the pressure gradient inside the cavity and was driven by the three-dimensional vortex ring generated from the "swallowtail" cavity structure.

Self-organized generation of transverse waves in diverging cylindrical detonations

Self-organized generation of transverse waves associated with the transverse wave instabilities at a diverging cylindrical detonation front was numerically studied by solving two-dimensional Euler equations implemented with an improved two-step chemical kinetic model. After solution validation, four mechanisms of the transverse wave generation were identified from numerical simulations, and referred to as the concave front focusing, the kinked front evolution, the wrinkled front evolution and the transverse wave merging, respectively. The propagation of the cylindrical detonation is maintained by the growth of the transverse waves that match the rate of increase in surface area of the detonation front to asymptotically approach a constant average number of transverse waves per unit length along the circumference of the detonation front. This cell bifurcation phenomenon of cellular detonations is discussed in detail to gain better understanding on detonation physics.

A Hybrid Scheme Based on Finite Element/Volume Methods for Two Immiscible Fluid Flows

We have successfully extended our implicit hybrid finite element/volume (FE/FV) solver to flows involving two immiscible fluids. The solver is based on the segregated pressure correction or projection method on staggered unstructured hybrid meshes. An intermediate velocity field is first obtained by solving the momentum equations with the matrix-free implicit cell-centered FV method. The pressure Poisson equation is solved by the node-based Galerkin FE method for an auxiliary variable. The auxiliary variable is used to update the velocity field and the pressure field. The pressure field is carefully updated by taking into account the velocity divergence field. This updating strategy can be rigorously proven to be able to eliminate the unphysical pressure boundary layer and is crucial for the correct temporal convergence rate. Our current staggered-mesh scheme is distinct from other conventional ones in that we store the velocity components at cell centers and the auxiliary variable at vertices. The fluid interface is captured by solving an advection equation for the volume fraction of one of the fluids. The same matrix-free FV method, as the one used for momentum equations, is used to solve the advection equation. We will focus on the interface sharpening strategy to minimize the smearing of the interface over time. We have developed and implemented a global mass conservation algorithm that enforces the conservation of the mass for each fluid.

A CE/SE Scheme for Flows in Porous Media and Its Applications

In this paper, a new computational scheme for solving flows in porous media was proposed. The scheme was based on an improved CE/SE method (the space-time Conservation Element and Solution Element method). We described porous flows by adopting DFB (Brinkman-Forchheimer extended Darcy) equation. The comparison between our computational results and Ghia's confirmed the high accuracy, resolution, and efficiency of our CE/SE scheme. The proposed first-order CE/SE scheme is a new reliable way for numerical simulations of flows in porous media. After investigation of effects of Darcy number on porous flow, it shows that Darcy number has dominant influence on porous flow for the Reynolds number and porosity considered.

Numerical simulation of condensation of bubbles under microgravity conditions by moving mesh method in the double-staggered grid

A numerical 2D method for simulation of two-phase flows including phase change under microgravity conditions is presented in this paper, with a level set method being coupled with the moving mesh method in the double-staggered grid systems. When the grid lines bend very much in a curvilinear grid, great errors may be generated by using the collocated grid or the staggered grid. So the double-staggered grid was adopted in this paper. The level set method is used to track the liquid-vapor interface. The numerical analysis is fulfilled by solving the Navier-Stokes equations using the SIMPLER method, and the surface tension force is modeled by a continuum surface force approximation. A comparison of the numerical results obtained with different numerical strategies shows that the double-staggered grid moving-mesh method presented in this paper is more accurate than that used previously in the collocated grid system. Based on the method presented in this paper, the condensation of a single bubble in the cold water under different level of gravity is simulated. The results show that the condensation process under the normal gravity condition is different from the condensation process under microgravity conditions. The whole condensation time is much longer under the normal gravity than under the microgravity conditions.

A numerical investigation on the flow physics of koi carp (Cyprinus carpio koi) burst-and-coast swimming

Abstract—Burst-and-coast is the most common locomotion type in freely routine swimming of koi carps (Cyprinus carpio koi), which consists of a burst phase and a coast phase in each cycle and mostly leads to a straight-line trajectory. Combining with the tracking experiment, the flow physics of koi carp’s burst-andcoast swimming is investigated using a novel integrated CFD method solving the body-fluid interaction problem. The dynamical equations of a deforming body are formulated. Following that, the loose-coupled equations of the body dynamics and the fluid dynamics are numerically solved with the integrated method. The two burst modes, MT (Multiple Tail-beat) and HT (Half Tail-beat), which have been reported by the experiments, are investigated by numerical simulations in this paper. The body kinematics is predicted and the flow physics is visualized, which are in good agreement with the corresponding experiments. Furthermore, the optimization on the energy cost and several critical control mechanisms in burst-and-coast swimming of koi carps are explored, by varying the parameters in its selfpropelled swimming. In this paper, energetics is measured by the two mechanical quantities, total output power CP and Froude efficiency Fr. Results and discussion show that from the standpoint of mechanical energy, burst-and-coast swimming does not actually save energy comparing with steady swimming at the same average speed, in that frequently changing of speed leads to decrease of efficiency.