An optimal theory for an expansion of flow quantities to capture the flow structures

An optimal theory on how database analysis to capture the flow structures has been developed in this paper, which include the POD method as its special case. By means of the remainder minimization method in the Sobolev space, for more general optimal conditions the new theory has the potential to overcome an inherent limitation of the POD method, i.e., it cannot be used to the situations in which the optimal condition is other than the inner product global one. As an example, using the new theory, the database of a two-dimensional flow over a backward-facing step is analyzed in detail, with velocity and vorticity bases.

Application of the compatibility relationship(1) in seeking a functional relationship among species' equilibrium quantities in complex systems

A new algorithm, representing an important advance in determination of the functional relationship, is first reported here. The algorithm is very useful and convenient for analyzing the incorporation of impurities. To show how the algorithm works, two early and well-known vapor phase epitaxy (VPE) experiments-Ashen's (Ashen, D. J.; Dean, P. J.; Hurle, D. T. J.; Mullin, J. B.; Royle, A.; White, A. M. Gallium Arsenide and Related Compounds, Institute of Physics Conference Series 24, 1974; Institute of Physics: London, 1975; p 229.), involving the doping of silicon and DiLorenzo's (DiLorenzo, J. V. J. Cryst. Growth 1972, 17, 189.), involving the mole fraction effect-are calculated to find the functional relationship between the Si contamination and the partial pressure of HCl. The calculated curves agree with the experimental results. A conclusion that the calculated values are greater than the true values has been drawn.

Multiple Axial Cracks in A Coated Hollow Cylinder due to Thermal Shock

The transient thermal stress problem of an inner-surface-coated hollow cylinder with multiple pre-existing surface cracks contained in the coating is considered. The transient temperature, induced thermal stress, and the crack tip stress intensity factor (SIF) are calculated for the cylinder via finite element method (FEM), which is exposed to convective cooling from the inner surface. As an example, the material pair of a chromium coating and an underlying steel substrate 30CrNi2MoVA is particularly evaluated. Numerical results are obtained for the stress intensity factors as a function of normalized quantities such as time, crack length, convection severity, material constants and crack spacing. (c) 2005 Elsevier Ltd. All rights reserved.

A Continuation Method of Parameter Inversion for Non-Equilibrium Convection-Dispersion Equation

Based on the homotopy mapping, a globally convergent method of parameter inversion for non-equilibrium convection-dispersion equations (CDEs) is developed. Moreover, in order to further improve the computational efficiency of the algorithm, a properly smooth function, which is derived from the sigmoid function, is employed to update the homotopy parameter during iteration. Numerical results show the feature of global convergence and high performance of this method. In addition, even the measurement quantities are heavily contaminated by noises, and a good solution can be found.

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.

Optimized Group Velocity Control Scheme and DNS of Decaying Compressible Turbulence of Relative High Turbulent Mach Number

To overcome the difficulty in the DNS of compressible turbulence at high turbulent Mach number, a new difference scheme called GVC8 is developed. We have succeeded in the direct numerical simulation of decaying compressible turbulence up to turbulent Mach number 0.95. The statistical quantities thus obtained at lower turbulent Mach number agree well with those from previous authors with the same initial conditions, but they are limited to simulate at lower turbulent Mach numbers due to the so-called start-up problem. The energy spectrum and coherent structure of compressible turbulent flow are analysed. The scaling law of compressible turbulence is studied. The computed results indicate that the extended self-similarity holds in decaying compressible turbulence despite the occurrence of shocklets, and compressibility has little effects on relative scaling exponents when turbulent Mach number is not very high.

Chemical kinetics and self-ignition in a model supersonic hydrogen-air combustor

Self-ignition tests of a model scramjet combustor were conducted by using parallel sonic injection of gaseous hydrogen from the base of a blade-like strut into a supersonic vitiated airstream. The range of stagnation pressure and temperature studied varied from 1.0 to 4.5 MPa and from 1300 to 2200 K, respectively. Experimental results show that the self-ignition limit, in terms of either global or local quantities of pressure and temperature, exhibits a nonmonotonic behavior resembling the classical homogeneous explosion limit of the hydrogen-oxygen system. Specifically, for a given temperature, increasing pressure from a low value can render a nonignitable mixture to first become ignitable, then nonignitable again, This correspondence shows that, despite the globally supersonic nonpremixed configuration studied herein, ignition is strongly influenced by the intricate chemical reaction mechanism and thereby exhibits the homogeneous explosion character. Consequently, self-ignition criteria based on a global reaction rate approximating the complex chemistry are inadequate. An auxiliary computational study on counterflow ignition was also conducted to systematically investigate the contamination effects of vitiated air. Results indicate that the net contamination effects for the present experimental data are expected to be substantially smaller than contributions from the individual contamination species because of the counterbalancing influences of the H2O-inhibition and NO-promotion reactions in effecting ignition.

Equivalency and Locality in Nano-Scale Measurement

Two principal problems of equivalency and locality in nano-scale measurement are considered in this paper. The conventional measurements of force and displacement are always closely related to the equivalency problem between the measuremental results by experimental system and the real physical status of the sample, and the locality of the mechanical quantities to be measured. There are some noticeable contradictions in nano-scale measurements induced by the two problems. In this paper, by utilizing a coupled molecular-continuum method, we illustrate the important effects of the two principal problems in atomic force microscopy (AFM) measurements on nano-scale. Our calculations and analysis of these typical mechanical measurement problems suggest that in nano-meter scale measurements, the two principal problems must be carefully dealt with. The coupled molecular-continuum method used in this paper is very effective in solving these problems on nano-scale.

A Microtube Viscometer with a Thermostat

The viscometer presented in this paper is suit-able for measuring the viscosity of liquids in micro-litre quantities. It consists of a micro-flow experimental system with a thermostat. Using the measurements of the flow rates and pressure drops of a liquid passing through a microtube, the liquid's viscosity can be calculated from on Hagen-Poiseuille theory. After calibration, the viscometer was used to measure viscosities of deionized water and ethyl alcohol at temperatures ranging from 0 to 40 "C. For both test liquids, the relative deviation of the measured values from those quoted in the literature (obtained using other viscometers) was less than 2.6o/o. The relative uncertainty of the experimental system was reduced to +-l.8% using the relative measuring method. Due to the micro-scale of the test section, only a micro-litre quantity of liquid is needed for a test, this is a potential advantage for measurement of bio-liquid viscosities.

Thermocapillary Convection in Floating Zones

This paper provides an overview of ongoing studies in the area of thermocapillary convection driven by a surface tension gradient parallel to the free surface in a floating zone. Here, research interests are focused around the onset of oscillatory thermocapillary convection, also known as the transition from quasisteady convection to oscillatory convection. The onset of oscillation depends on a set of critical parameters, and the margin relationship can be represented by a complex function of the critical parameters. The experimental results indicate that the velocity deviation of an oscillatory flow has the same order of magnitude as that of an average flow, and the deviations of other quantities, such as temperature and free surface radii fluctuations, are much smaller when compared with their normal counterparts. Therefore, the onset of oscillation should be a result of the dynamic process in a fluid, and the problem is a strongly nonlinear one. In the past few decades, several theoretical models have been introduced to tackle the problem using analytical methods, linear instability analysis methods, energy instability methods, and unsteady 3D numerical methods. The last of the above mentioned methods is known to be the most suitable for a thorough analysis of strong nonlinear processes, which generally leads to a better comparison with the experimental results. The transition from oscillatory thermocapillary convection to turbulence falls under the studies of chaotic behavior in a new system, which opens a fascinating new frontier in nonlinear science, a hot research area drawing many recent works. This paper reviews theoretical models and analysis, and also experimental research, on thermocapillary connection in floating zones. It cites 93 references.

Energy Principle And Nonlinear Electric-Mechanical Behavior Of Ferroelectric Ceramics

Many experimental observations have shown that a single domain in a ferroelectric material switches by progressive movement of domain walls, driven by a combination of electric field and stress. The mechanism of the domain switch involves the following steps: initially, the domain has a uniform spontaneous polarization; new domains with the reverse polarization direction nucleate, mainly at the surface, and grow though the crystal thickness; the new domain expands sideways as a new domain continues to form; finally, the domain switch coalesces to complete the polarization reversal. According to this mechanism, the volume fraction of the domain switching is introduced in the constitutive law of the ferroelectric material and used to study the nonlinear constitutive behavior of a ferroelectric body in this paper. The principle of stationary total potential energy is put forward in which the basic unknown quantities are the displacement u(i), electric displacement D-i and volume fraction rho(I) of the domain switching for the variant I. The mechanical field equation and a new domain switching criterion are obtained from the principle of stationary total potential energy. The domain switching criterion proposed in this paper is an expansion and development of the energy criterion established by Hwang et al. [ 1]. Based on the domain switching criterion, a set of linear algebraic equations for determining the volume fraction rho(I) of domain switching is obtained, in which the coefficients of the linear algebraic equations only contain the unknown strain and electric fields. If the volume fraction rho(I) of domain switching for each domain is prescribed, the unknown displacement and electric potential can be obtained based on the conventional finite element procedure. It is assumed that a domain switches if the reduction in potential energy exceeds a critical energy barrier. According to the experimental results, the energy barrier will strengthen when the volume fraction of the domain switching increases. The external mechanical and electric loads are increased step by step. The volume fraction rho(I) of domain switching for each element obtained from the last loading step is used as input to the constitutive equations. Then the strain and electric fields are calculated based on the conventional finite element procedure. The finite element analysis is carried out on the specimens subjected to uniaxial coupling stress and electric field. Numerical results and available experimental data are compared and discussed. The present theoretic prediction agrees reasonably with the experimental results.

Avalanche behavior and statistical properties in a microcrack coalescence process

Fracture owing to the coalescence of numerous microcracks can be described by a simple statistical model, where a coalescence event stochastically occurs as the number density of nucleated microcracks increases. Both numerical simulation and statistical analysis reveal that a microcrack coalescence process may display avalanche behavior and that the final failure is catastrophic. The cumulative distribution of coalescence events in the vicinity of critical fracture follows a power law and the fracture profile has self-affine fractal characteristic. Some macromechanical quantities may be traced back and extracted from the mesoscopic process based on the statistical analysis of coalescence events.

Compact difference approximation with consistent boundary condition

For simulating multi-scale complex flow fields it should be noted that all the physical quantities we are interested in must be simulated well. With limitation of the computer resources it is preferred to use high order accurate difference schemes. Because of their high accuracy and small stencil of grid points computational fluid dynamics (CFD) workers pay more attention to compact schemes recently. For simulating the complex flow fields the treatment of boundary conditions at the far field boundary points and near far field boundary points is very important. According to authors' experience and published results some aspects of boundary condition treatment for far field boundary are presented, and the emphasis is on treatment of boundary conditions for the upwind compact schemes. The consistent treatment of boundary conditions at the near boundary points is also discussed. At the end of the paper are given some numerical examples. The computed results with presented method are satisfactory.

Energy principle of ferroelectric ceramics and single domain mechanical model

Many physical experiments have shown that the domain switching in a ferroelectric material is a complicated evolution process of the domain wall with the variation of stress and electric field. According to this mechanism, the volume fraction of the domain switching is introduced in the constitutive law of ferroelectric ceramic and used to study the nonlinear constitutive behavior of ferroelectric body in this paper. The principle of stationary total energy is put forward in which the basic unknown quantities are the displacement u (i) , electric displacement D (i) and volume fraction rho (I) of the domain switching for the variant I. Mechanical field equation and a new domain switching criterion are obtained from the principle of stationary total energy. The domain switching criterion proposed in this paper is an expansion and development of the energy criterion. On the basis of the domain switching criterion, a set of linear algebraic equations for the volume fraction rho (I) of domain switching is obtained, in which the coefficients of the linear algebraic equations only contain the unknown strain and electric fields. Then a single domain mechanical model is proposed in this paper. The poled ferroelectric specimen is considered as a transversely isotropic single domain. By using the partial experimental results, the hardening relation between the driving force of domain switching and the volume fraction of domain switching can be calibrated. Then the electromechanical response can be calculated on the basis of the calibrated hardening relation. The results involve the electric butterfly shaped curves of axial strain versus axial electric field, the hysteresis loops of electric displacement versus electric filed and the evolution process of the domain switching in the ferroelectric specimens under uniaxial coupled stress and electric field loading. The present theoretic prediction agrees reasonably with the experimental results given by Lynch.

Transition of oscillatory floating half zone convection from Earth's gravity to microgravity

Oscillatory features of floating half zone convection were experimentally studied by using the drop shaft facility of Japan Microgravity Center which supported microgravity period of 10 s. Coordinated measurements including free surface deformation and oscillation, temperature and flow pattern in both 1-g and micro-g environment were obtained. The oscillatory frequency and amplitude in micro-g condition were lower and larger than the ones in l-g condition, respectively. The results gave, at first time, the oscillatory features such as free surface wave in micro-g, coordinated measurements of more than two physical quantities in the micro-g, and transition of thermocapillary oscillatory convection from I-g to micro-g.