Numerical Investigation of Toroidal Shock Wave Focusing

In this paper,focusing of a toroidal shock wave propagating from a shock tube of an- nular cross-section into a cylindrical chamber was investigated numerically with the dispersion- controlled scheme. For CFD validation, the numerical code was rst applied to calculate both viscous and inviscid ows at a low Mach number of 1.5, which was compared with the experi- ment results and got better consistency. Then the validated code was used to calculate several cases for high Mach numbers. From the result, several major factors that in uent the ow, such as the Mach number and the viscosity, were analyzed detailedly and along with the high Mach number some unusual ow structure was observed and explained theoretically

Cone Effect in Astronomical Adaptive Optics System Investigated by a Pure Numerical Simulation

It is well-known that cone effect or focus anisoplanatism is produced by the limited distance of a laser guide star (LGS) which is created within the Earth atmosphere and consequently located at a finite distance from the observer. In this paper, the cone effect of the LGS for different vertical profiles of the refractive index structure constant Cn2 is numerically investigated by using a revised computer program of atmospheric propagation of optical wave and an adaptive optics (AO) system including dynamic control process. According to the practice, the overall tilt for the tilt-correction mirror is obtained from a natural star and the aberrated wavefront for phase correction of the deformable mirror is obtained from a LGS in our numerical simulation. It is surprisingly found that the effect of altitude of the LGS on the AO phase compensation effectiveness by using the commonly-available vertical profiles of Cn2 and the lateral wind speed in the atmosphere is relatively weak, and the cone effect for some Cn2 profiles is even negligible. It is found that the cone effect does not have obvious relationship with the turbulence strength, however, it depends on the vertical distribution profile of Cn 2 apparently. On the other hand, the cone effect depends on the vertical distribution of the lateral wind speed as well. In comparison to a longer wavelength, the cone effect becomes more obvious in the case of a shorter wavelength. In all cases concerned in this paper, an AO system by using a sodium guide star has almost same phase compensation effectiveness as that by using the astronomical target itself as a beacon. Effect of dynamic control process in an AO system on the cone effect is studied in this paper for the first time within our knowledge.

Numerical Study on Kerosene Spray in Supersonic Flow

Numerical simulation was conducted to study the kerosene spray characteristics injecting into supersonic cross flow. The verification of the simulation was carried out by experimental Schlieren image, and the agreement was obtained by compared the spray plume pictures. Furthermore, the aerodynamic secondary breakup effect of the supersonic cross flow on the initial droplets was investigated. 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 Simulation of the Flow Field and Concentration Distribution in the Bulk Growth of Silicon Carbide Crystals

The physical vapor transport (PVT) method is being widely used to grow large-size single SiC crystals. The growth process is associated with heat and mass transport in the growth chamber, chemical reactions among multiple species as well as phase change at the crystal/gas interface. The current paper aims at studying and verifying the transport mechanism and growth kinetics model by demonstrating the flow field and species concentration distribution in the growth system. We have developed a coupled model, which takes into account the mass transport and growth kinetics. Numerical simulation is carried out by employing an in-house developed software based on finite volume method. The results calculated are in good agreement with the experimental observation.

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.

Numerical Analysis of the Effects of Free-Stream Thermo-Chemical State on the Test Model Flow Field in High-Enthalpy Tunnel

The effects of the free-stream thermo-chemical state on the test model flow field in the high-enthalpy tunnel are studied numerically. The properties of the free-stream, which is in thermo-chemical non-equilibrium, are determined by calculating the nozzle flow field. A free-stream with total enthalpy equal to the real one in the tunnel while in thermo-chemical equilibrium is constructed artificially to simulate the natural atmosphere condition. The flow fields over the test models (blunt cone and Apollo command capsule model) under both the non-equilibrium and the virtual equilibrium free-stream conditions are calculated. By comparing the properties including pressure, temperature, species concentration and radiation distributions of these two types of flow fields, the effects of the non-equilibrium state of the free-stream in the high-enthalpy shock tunnel are analyzed.

Numerical Study of Heat Transfer Mechanism in Turbulent Supercritical CO2 Channel Flow

Direct numerical simulation (DNS) of supercritical CO2 turbulent channel flow has been performed to investigate the heat transfer mechanism of supercritical fluid. In the present DNS, full compressible Navier-Stokes equations and Peng-Robison state equation are solved. Due to effects of the mean density variation in the wall normal direction, mean velocity in the cooling region becomes high compared with that in the heating region. The mean width between high-and low-speed streaks near the wall decreases in the cooling region, which means that turbulence in the cooling region is enhanced and lots of fine scale eddies are created due to the local high Reynolds number effects. From the turbulent kinetic energy budget, it is found that compressibility effects related with pressure fluctuation and dilatation of velocity fluctuation can be ignored even for supercritical condition. However, the effect of density fluctuation on turbulent kinetic energy cannot be ignored. In the cooling region, low kinematic viscosity and high thermal conductivity in the low speed streaks modify fine scale structure and turbulent transport of temperature, which results in high Nusselt number in the cooling condition of the supercritical CO2.