Numerical Simulation and Experimental Investigation about Internal and External Flows

In this paper, TASCflow3D is used to solve inner and outer 3D viscous incompressible turbulent flow (R-e = 5.6 X 10(6)) around axisymmetric body with duct. The governing equation is a RANS equation with standard k-epsilon turbulence model. The discrete method used is a finite volume method based on the finite element approach. In this method, the description of geometry is very flexible and at the same time important conservative properties are retained. The multi-block and algebraic multi-grid techniques are used for the convergence acceleration. Agreement between experimental results and calculation is good. It indicates that this novel approach can be used to simulate complex flow such as the interaction between rotor and stator or propulsion systems containing tip clearance and cavitation.

Numerical Simulation of Heat Transfer and Kinetics in the Bulk Growth of Silicon Carbide

The growth process of 2-inch silicon carbide (SiC) single crystals by the physical vapor transport method (or modified Lely method) has been modeled and simulated. The comprehensive process model incorporates the calculations of radio frequency (RF) induction heating, heat and mass transfer and growth kinetics. The transport equations for electromagnetic field, heat transfer, and species transport are solved using a finite volume-based numerical scheme called MASTRAPP (Multizone Adaptive Scheme for Transport and Phase Change Process). Temperature distribution for a 2-inch growth system is calculated, and the effects of induction heating frequency and current on the temperature distribution and growth rate are investigated. The predicted results have been compared with the experimental data.

Numerical Simulation of R-M Instability

In order to capture shock waves and contact discontinuities in the field and easy to program with parallel computation a new algorithm is developed to solve the N-S equations for simulation of R-M instability problems. The method with group velocity control is used to suppress numerical oscillations, and an adaptive non-uniform mesh is used to get fine resolution. Numerical results for cylindrical shock-cylindrical interface interaction with a shock Mach number Ms=1.2 and Atwood number A=0.818, 0.961, 0.980 (the interior density of the interface/outer density p(1)/p(2) = 10, 50, 100, respectively), and for the planar shock-spherical interface interaction with Ms=1.2 and p(1)/p(2) = 14.28are presented. The effect of Atwood number and multi-mode initial perturbation on the R-M instability are studied. Multi-collisions of the reflected shock with the interface is a main reason of nonlinear development of the interface instability and formation of the spike-bubble structures In simulation with double mode perturbation vortex merging and second instability are found. After second instability the small vortex structures near the interface produced. It is important factor for turbulent mixing.

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

Numerical Research on Liquid Spray in Supersonic Cross Flow

Numerical simulation was conducted to characterize the kerosene spray injecting into supersonic cross flow, especially focusing on the aerodynamic secondary breakup effect of the supersonic cross flow on the initial droplets. It was revealed that the initial parent drops were broken up into small drops whose diameter is about O(10) micrometers soon after they entered into the supersonic cross flow. During the appropriate range of initial drop size, the parent droplets would be broken up into small drops with the same magnitude diameter no matter how large the initial drops SMD was.

Numerical Solution for Polymer Transient Flows in a Circle Bounded Composite Formation

A new numerical model for transient flows of polymer solution in a circular bounded composite formation is presented in this paper. Typical curves of the wellbore transient pressure are yielded by FEM. The effects of non-Newtonian power-law index, mobility and boundary distance have been considered. It is found that for the mobility ratio larger than 1, which is favorable for the polymer flooding, the pressure derivative curve in log-log form rises up without any hollow. On the other hand, if the pressure derivative curve has a hollow and then is raised up, we say that the polymer flooding fails. Finally, the new model has been extended to more complicated boundary case.