Design of a Self-Adaptive Down-Hole Layered Steam Injector and Its Application

Layered steam injection, widely used in Liaohe Oilfield at present, is an effective recovery technique to heavy oil reserves. Which makes the steam front-peak push forward uniformly, the amount of steam injection be assigned rationally, and the effect of injection steam be obtained as expected. To maintain a fixed ratio of layered steam injection and solve the problem of nonadjustable hole diameter with the change of layer pressure in the existing injectors, a new method is proposed in this paper to design layered steam injectors based on the dynamic balance theory According to gas-liquid two-phase flow theory and heat transfer theory, the energy equation and the heat conduction equation in boreholes are developed. By analyzing the energy equilibrium of water-steam passing through the injector hole, we find an expression to describe the relation between the cross-sectional area of injector hole and the layer pressure. With this expression, we provide a new set of calculation methods and write the corresponding computer program to design and calculate the main parameters of a steam injector. The actual measurement data show that the theoretically calculated results are accurate, the software runs reliably, and they provide the design of self-adjustable layered steam injectors with the theoretical foundation.

Transient Temperature Distribution in a Sheet Metal by Low-Energy Laser Scanning

Low-energy laser-heating techniques are widely used in engineering applications such as, thinfilm deposition, surface treatment, metal forming and micro-structural pattern formation. In this paper,under the conditions of ignoring the thermo-mechanical coupling, a numerical simulation on the spatialand temporal temperature distribution in a sheet metal produced by the laser beam scanning in virtue of thefinite element method is presented. Both the three-dimensional transient temperature field and thetemperature evolution as a function of heat penetrating depth in the metal sheet are calculated. Thetemperature dependence of material properties was taken into account. It was shown that, after taking thetemperature dependence of the material absorbance effect into consideration, the temperature change ratealong the scanning direction and the temperature maximum were both increased.

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<

Numerical Analyses of Transonic Fluid/Structure interaction in Aircraft Engineering

Through the coupling between aerodynamic and structural governing equations, a fully implicit multiblock aeroelastic solver was developed for transonic fluid/stricture interaction. The Navier-Stokes fluid equations are solved based on LU-SGS (lower-upper symmetric Gauss-Seidel) Time-marching subiteration scheme and HLLEW (Harten-Lax-van Leer-Einfeldt-Wada) spacing discretization scheme and the same subiteration formulation is applied directly to the structural equations of motion in generalized coordinates. Transfinite interpolation (TFI) is used for the grid deformation of blocks neighboring the flexible surfaces. The infinite plate spline (IPS) and the principal of virtual work are utilized for the data transformation between fluid and structure. The developed code was fort validated through the comparison of experimental and computational results for the AGARD 445.6 standard aeroelastic wing. In the subsonic and transonic range, the calculated flutter speeds and frequencies agree well with experimental data, however, in the supersonic range, the present calculation overpredicts the experimental flutter points similar to other computations. Then the flutter character of a complete aircraft configuration is analyzed through the calculation of the change of structural stiffness. Finally, the phenomenon of aileron buzz is simulated for the weakened model of a supersonic transport wing/body model at Mach numbers of 0.98 and l.05. The calculated unsteady flow shows, on the upper surface, the shock wave becomes stronger as the aileron deflects downward, and the flow behaves just contrary on the lower surface of the wing. Corresponding to general theoretical analysis, the flow instability referred to as aileron buzz is induced by a stronger shock alternately moving on the upper and lower surfaces of wing. For the rigid structural model, the flow is stable at all calculated Mach numbers as observed in experiment

Chemical equilibrium analysis and CFD simulation of coal combustion and NOx formation

A full two-fluid model of reacting gas-particle flows with an algebraic unified second-order moment (AUSM) turbulence-chemistry model is used to simulate Beijing coal combustion and NOx formation. The sub-models are the k-epsilon-kp two-phase turbulence model, the EBU-Arrhenius volatile and CO combustion model, the six-flux radiation model, coal devolatilization model and char combustion model. The blocking effect on NOx formation is discussed. In addition, the chemical equilibrium analysis is used to predict NOx concentration at different temperature. Results of CID simulation and chemical equilibrium analysis show that, optimizing air dynamic parameters can delay the NOx formation and decrease NOx emission, but it is effective only in a restricted range. In order to decrease NOx emission near to zero, the re-burning or other chemical methods must be used.

3D Finite Volume Scheme for Czochralski Crystal Growth

A general three-dimensional model is developed for simulation of the growth process of silicon single crystals by Czochralski technique. The numerical scheme is based on the curvilinear non-orthogonal finite volume discretization. Numerical solutions show that the flow and temperature fields in the melt are asymmetric and unsteady for 8’’ silicon growth. The effects of rotation of crystal on the flow structure are studied. The rotation of crystal forms the Ekman layer in which the temperature gradient along solid/melt surface is small.

Numerical Simulation of Pulsed Laser Transformation Hardening

A novel pulsed laser surface processing technology is introduced, which can make use of the spatial and temporal profile of laser pulse to obtain ideal hardening parameters. The intensity distribution of laser pulse is spatially and temporally controlled by using laser shape transformation technology. A 3D numerical model including multi-phase transformations is established to explore material microstructure evolution induced by temperature field evolution. The influences of laser spatial-temporal profiles on hardening parameters are investigated. Different from the continuous laser processing technology, results indicate that spatial and temporal profiles are important factors in determining processing quality during pulsed laser processing method.

3D Finite Volume Scheme for Czochralski Single Crystal Growth

Czochralski (Cz) technique, which is used for growing single crystals, has dominated the production of single crystals for electronic applications. The Cz growth process involves multiple phases, moving interface and three-dimensional behavior. Much has been done to study these phenomena by means of numerical methods as well as experimental observations. A three-dimensional curvilinear finite volume based algorithm has been developed to model the Cz process. A body-fitted transformation based approach is adopted in conjunction with a multizone adaptive grid generation (MAGG) technique to accurately handle the three-dimensional problems of phase-change in irregular geometries with free and moving surfaces. The multizone adaptive model is used to perform a three-dimensional simulation of the Cz growth of silicon single crystals.Since the phase change interface are irregular in shape and they move in response to the solution, accurate treatment of these interfaces is important from numerical accuracy point of view. The multizone adaptive grid generation (MAGG) is the appropriate scheme for this purpose. Another challenge encountered is the moving and periodic boundary conditions, which is essential to the numerical solution of the governing equations. Special treatments are implemented to impose the periodic boundary condition in a particular direction and to determine the internal boundary position and shape varying with the combination of ambient physicochemical transport process and interfacial dynamics. As indicated above that the applications and processes characterized by multi-phase, moving interfaces and irregular shape render the associated physical phenomena three-dimensional and unsteady. Therefore a generalized 3D model rather than a 2D simulation, in which the governing equations are solved in a general non-orthogonal coordinate system, is constructed to describe and capture the features of the growth process. All this has been implemented and validated by using it to model the low pressure Cz growth of silicon. Accuracy of this scheme is demonstrated by agreement of simulation data with available experimental data. Using the quasi-steady state approximation, it is shown that the flow and temperature fields in the melt under certain operating conditions become asymmetric and unsteady even in the absence of extrinsic sources of asymmetry. Asymmetry in the flow and temperature fields, caused by high shear initiated phenomena, affects the interface shape in the azimuthal direction thus results in the thermal stress distribution in the vicinity, which has serious implications from crystal quality point of view.

Numeric Calculation for Transient Thermal Stress Field in Wall-Shape Metal Part During Laser Direct Forming

In order to investigate the transient thermal stress field in wall-shape metal part during laser direct forming, a FEM model basing on ANSYS is established, and its algorithm is also dealt with. Calculation results show that while the wall-shape metal part is being deposited, in X direction, the thermal stress in the top layer of the wall-shape metal part is tensile stress and in the inner of the wall-shape metal part is compressive stress. The reason causing above-mentioned thermal stress status in the wall-shape metal part is illustrated, and the influence of the time and the processing parameters on the thermal stress field in wall-shape metal part is also studied. The calculation results are consistent with experimental results in tendency.

A New Model for Sediment Transport

The problem of predicting sediment transportation by water waves is treated analytically with the rate of wave energy dissipation or wave damping. With resorting to the theory of shallow water waves and the basis of Yamamoto’s Coulomb-damped poroelastic model, the Boussinesq-type equation has been derived over a variation depth bed. For convenience Cnoidal wave is just discussed, The Cnoidal wave with complex wave length and wave velocity, which are as a function of wave frequency, water depth, permeability, Poisson’s ratio and complex elastic moduli of bed soil, is applied to analyse the rate of sediment transportation. Considering the sediment transportation depended on the shear stress near-bed or the horizontal velocity, the conclusion of Yamamoto’s experiment in clay bed has been extended to general situation. It could be figured out that the model should provide a method to avoid the undistinguishable factors during sediment transport processes and relate mass transport with the sediment peculiarities.

Numerical Partition of Energy Absorption of Foam-filled Top-Hat and Double-Hat Sections

The interaction effect, i.e., the contribution of each component to the total energy absorption of an axially crushed foam-filled hat section was investigated quantitatively via numerical simulation. The FE results were first verified by experimental work of aluminum foam-filled top-hat and double-hat sections, then the contribution of foam-fillers and that of hat sections to the overall energy absorption were quantitatively obtained, respectively. When foam-filled, increase in energy absorption was found both in hat section component and foam-filler component, whereas the latter contributes predominantly to the interaction effect.

Effect of Pipe Diameter on Flow Pattern Transition in Vertical Oil Water Flows

This paper investiges the effect of pipe diameter on flow pattern transition boundary in oil water vertical flows, and proposes a model to determine the maximum inner diameter (D_{infty s}) of a pipe in which the slug flow would not occur When pipe inner diameter D>D_{infty s}, only bubble flow exists, while D

Simulation of Vortex-dominant and Turbulent Flows in Natural Environment

In the present talk, the simulation of vortex dominant and turbulent flows are primarily addressed. To cope with complicated circumstances in environmental flows we illustrate the strategy of combining simplified physical model and suitable algorithm by a few examples.