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.

Experimental Study and Numerical Simulation of Cellular Structures and Mach Reflection of Gaseous Detonation Wave

In this paper the Deflagration to Detonation Transition (DDT) process of gaseous H-2-O-2 mixture and Mach reflection of gaseous detonation wave on a wedge have been conducted experimentally. The cellular pattern of DDT process and Mach reflection were obtained from experiments with wedge angle theta = 10(0) similar to 40(0) and initial pressure of gaseous mixture 16kPa similar to 26.7kPa. The 2-D numerical simulations of DDT process and Mach reflection of detonation wave were performed by using the simplified ZND model and improved space-time conservation element and solution element (CE/SE) method. The numerical cellular structures were compared with the cellular patterns of soot track. Compared results were shown that it is satisfactory. The characteristic comparisons on Mach reflection of air shock wave and detonation wave were carried also out and their differences were given.

Transonic Aeroelastic Numerical Simulation in Aeronautical Engineering

A lower-upper symmetric Gauss-Seidel (LU-SGS) subiteration scheme is constructed for time-marching of the fluid equations. The Harten-Lax-van Leer-Einfeldt-Wada (HLLEW) scheme is used for the spatial discretization. The same subiteration formulation is applied directly to the structural equations of motion in generalized coordinates. Through subiteration between the fluid and structural equations, a fully implicit aeroelastic solver is obtained for the numerical simulation of fluid/structure interaction. To improve the ability for application to complex configurations, a multiblock grid is used for the flow field calculation and transfinite interpolation (TFI) is employed for the adaptive moving grid deformation. The infinite plate spline (IPS) and the principal of virtual work are utilized for the data transformation between the fluid and structure. The developed code was first validated through the comparison of experimental and computational results for the AGARD 445.6 standard aeroelastic wing. Then, the flutter character of a tail wing with control surface was analyzed. Finally, flutter boundaries of a complex aircraft configuration were predicted.

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 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 simulation of a single fluid particle rising in non-Newtonian fluids

In this work, a level set method is developed for simulating the motion of a fluid particle rising in non-Newtonian fluids described by generalized Newtonian as well as viscoelastic model fluids. As the shear-thinning model we use a Carreau-Yasuda model, and the viscoelastic effect can be modeled with Oldroyd-B constitutive equations. The control volume formulation with the SIMPLEC algorithm incorporated is used to solve the governing equations on a staggered Eulerian grid. The level set method is implemented to compute the motion of a bubble in a Newtonian fluid as one of typical examples for validation, and the computational results are in good agreement with the reported experimental data．The level set method is also applied for simulating a Newtonian drop rising in Carreau-Yasuda and Oldroyd-B fluids．Numerical results including noticeably negative wake behind the drop and viscosity field are obtained, and compare satisfactorily with the known literature data.

Direct numerical simulation of hypersonic boundary layer transition over a blunt cone with a small angle of attack

The direct numerical simulation of boundary layer transition over a 5° half-cone-angle blunt cone is performed. The free-stream Mach number is 6 and the angle of attack is 1°. Random wall blow-and-suction perturbations are used to trigger the transition. Different from the authors’ previous work [Li et al., AIAA J. 46, 2899(2008)], the whole boundary layer flow over the cone is simulated (while in the author’s previous work, only two 45° regions around the leeward and the windward sections are simulated). The transition location on the cone surface is determined through the rapid increase in skin fraction coefficient (Cf). The transition line on the cone surface shows a nonmonotonic curve and the transition is delayed in the range of 0° ≤ θ ≤ 30° (θ = 0° is the leeward section). The mechanism of the delayed transition is studied by using joint frequency spectrum analysis and linear stability theory (LST). It is shown that the growth rates of unstable waves of the second mode are suppressed in the range of 20° ≤ θ ≤ 30°, which leads to the delayed transition location. Very low frequency waves VLFWs� are found in the time series recorded just before the transition location, and the periodic times of VLFWs are about one order larger than those of ordinary Mack second mode waves. Band-pass filter is used to analyze the low frequency waves, and they are deemed as the effect of large scale nonlinear perturbations triggered by LST waves when they are strong enough.The direct numerical simulation of boundary layer transition over a 5° half-cone-angle blunt cone is performed. The free-stream Mach number is 6 and the angle of attack is 1°. Random wall blow-and-suction perturbations are used to trigger the transition. Different from the authors’ previous work [ Li et al., AIAA J. 46, 2899 (2008) ], the whole boundary layer flow over the cone is simulated (while in the author’s previous work, only two 45° regions around the leeward and the windward sections are simulated). The transition location on the cone surface is determined through the rapid increase in skin fraction coefficient (Cf). The transition line on the cone surface shows a nonmonotonic curve and the transition is delayed in the range of 20° ≤ θ ≤ 30° (θ = 0° is the leeward section). The mechanism of the delayed transition is studied by using joint frequency spectrum analysis and linear stability theory (LST). It is shown that the growth rates of unstable waves of the second mode are suppressed in the range of 20° ≤ θ ≤ 30°, which leads to the delayed transition location. Very low frequency waves (VLFWs) are found in the time series recorded just before the transition location, and the periodic times of VLFWs are about one order larger than those of ordinary Mack second mode waves. Band-pass filter is used to analyze the low frequency waves, and they are deemed as the effect of large scale nonlinear perturbations triggered by LST waves when they are strong enough.

An Improved CE/SE Scheme for Numerical Simulation of Gaseous and Two-Phase Detonations

A new structure of solution elements and conservation elements based on rectangular mesh was pro- posed and an improved space-time conservation element and solution element (CE/SE) scheme with sec- ond-order accuracy was constructed. Furthermore, the application of improved CE/SE scheme was extended to detonation simulation. Three models were used for chemical reaction in gaseous detonation. And a two-fluid model was used for two-phase (gas–droplet) detonation. Shock reflections were simu- lated by the improved CE/SE scheme and the numerical results were compared with those obtained by other different numerical schemes. Gaseous and gas–droplet planar detonations were simulated and the numerical results were carefully compared with the experimental data and theoretical results based on C–J theory. Mach reflection of a cellular detonation was also simulated, and the numerical cellular pat- terns were compared with experimental ones. Comparisons show that the improved CE/SE scheme is clear in physical concept, easy to be implemented and high accurate for above-mentioned problems.

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.

Numerical Simulation of Liquefaction Deformation of Sand Layer around Bucket Foundations under Dynamic Loadings

Bucket Foundations under Dynamic Loadings The liquefaction deformation of sand layer around a bucket foundation is simulated under equivalent dynamic ice-induced loadings. A simplified numerical model is presented by taking the bucket-soil interaction into consideration. The development of vertical and horizontal liquefaction deformations are computed under equivalent dynamic ice-induced loadings. Firstly, the numerical model and results are proved to be reliable by comparing them with the centrifuge testing results. Secondly, the factors and the development characteristics of liquefaction deformation are analyzed. Finally, the following numerical simulation results are obtained: the liquefaction deformation of sand layer increases with the increase of loading amplitude and with the decrease of loading frequency and sand skeleton’s strength. The maximum vertical deformation is located on the sand layer surface and 1/4 times of the bucket’s height apart from the bucket’s side wall (loading boundary). The maximum horizontal deformation occurs at the loading boundary. When the dynamic loadings is applied for more than 5 hours, the vertical deformation on the sand layer surface reaches 3 times that at the bottom, and the horizontal deformation at 2.0 times of the bucket height apart from the loading boundary is 3.3% of which on the loading boundary.

Numerical simulation of wavefronts diffracted by apertures with circular symmetry

The numerical simulation of the wavefronts diffracted by apertures with circular symmetry is realized by a numerical method. It is based on the angular spectrum of plane waves, which ignored the vector nature of light. The on-axial irradiance distributions of plane wavefront and Gauss wavefront diffracted by the circular aperture have been calculated along the propagation direction. Comparisons of the simulation results with the analytical results and the experimental results tell us that it is a feasible method to calculate the diffraction of apertures. (c) 2006 Published by Elsevier GmbH.

Study of the surface plasma transmission properties of a Fabry-Perot resonator by numerical simulation

It is demonstrated with powerful evidence that the extraordinary transmission of a metallic grating is undoubtedly due to the excitation of standing surface plasma waves in the Fabry-Perot like resonator. This is the first time that the strong standing waves set up in the groove of a sub-wavelength double-layer grating (SWDG) for the surface plasma waves have been reported. Moreover, about 90% transmission is gained with an SWDG, more easily fabricated than ordinary metallic gratings, in the first peak of transmission spectrum.

Photon iterative numerical technique for steady-state simulation of gain-clamped semiconductor optical amplifiers

The photon iterative numerical technique, which chooses the outputs of the amplified spontaneous emission spectrum and lasing mode as iteration variables to solve the rate equations, is proposed and applied to analyse the steady behaviour of conventional semiconductor optical amplifiers (SOAs) and gain-clamped semiconductor optical amplifiers (GCSOAs). Numerical results show that the photon iterative method is a much faster and more efficient algorithm than the conventional approach, which chooses the carrier density distribution of the SOAs as the iterative variable. It is also found that the photon iterative method has almost the same computing efficiency for conventional SOAs and GCSOAs.

Symmetry analysis and numerical simulation of mode characteristics for equilateral-polygonal optical microresonators

Mode characteristics for two-dimensional equilateral-polygonal microresonators are investigated based on symmetry analysis and finite-difference time-domain numerical simulation. The symmetries of the resonators can be described by the point group C-Nv, accordingly, the confined modes in these resonators can be classified into irreducible representations of the point group C-Nv. Compared with circular resonators, the modes in equilateral-polygonal resonators have different characteristics due to the break of symmetries, such as the split of double-degenerate modes, high field intensity in the center region, and anomalous traveling-wave modes, which should be considered in the designs of the polygonal resonator microlasers or optical add-drop filters.

Numerical simulation of nc-Si : H/c-Si heterojunction solar cells

AMPS simulator, which was developed by Pennsylvania State University, has been used to simulate photovoltaic performances of nc-Si:H/c-Si solar cells. It is shown that interface states are essential factors prominently influencing open circuit voltages (V-OC) and fill factors (FF) of these structured solar cells. Short circuit current density (J(SC)) or spectral response seems more sensitive to the thickness of intrinsic a-Si:H buffer layers inserted into n(+)-nc-Si:H layer and p-c-Si substrates. Impacts of bandgap offset on solar cell performances have also been analyzed. As DeltaE(C) increases, degradation of VOC and FF owing to interface states are dramatically recovered. This implies that the interface state cannot merely be regarded as carrier recombination centres, and impacts of interfacial layer on devices need further investigation. Theoretical maximum efficiency of up to 31.17% (AM1.5,100mW/cm(2), 0.40-1.1mum) has been obtained with BSF structure, idealized light-trapping effect(R-F=0, R-B=1) and no interface states.