Scheme Construction with Numerical Flux Residual Correction (NFRC) and Group Velocity Control (GVC)

For simulating multi-scale complex flow fields like turbulent flows, the high order accurate schemes are preferred. In this paper, a scheme construction with numerical flux residual correction (NFRC) is presented. Any order accurate difference approximation can be obtained with the NFRC. To improve the resolution of the shock, the constructed schemes are modified with group velocity control (GVC) and weighted group velocity control (WGVC). The method of scheme construction is simple, and it is used to solve practical problems.

Numerical simulation of rock failure and earthquake process on mesoscopic scale

On the basis of the lattice model of MORA and PLACE, Discrete Element Method, and Molecular Dynamics approach, another kind of numerical model is developed. The model consists of a 2-D set of particles linked by three kinds of interactions and arranged into triangular lattice. After the fracture criterion and rules of changes between linking states are given, the particle positions, velocities and accelerations at every time step are calculated using a finite-difference scheme, and the configuration of particles can be gained step by step. Using this model, realistic fracture simulations of brittle solid (especially under pressure) and simulation of earthquake dynamics are made.

Numerical Solutions for the Transient Flow in the Homogenous Closed Circle Reservoirs

There are many fault block fields in China. A fault block field consists of fault pools. The small fault pools can be viewed as the closed circle reservoirs in some case. In order to know the pressure change of the developed formation and provide the formation data for developing the fault block fields reasonably, the transient flow should be researched. In this paper, we use the automatic mesh generation technology and the finite element method to solve the transient flow problem for the well located in the closed circle reservoir, especially for the well located in an arbitrary position in the closed circle reservoir. The pressure diffusion process is visualized and the well-location factor concept is first proposed in this paper. The typical curves of pressure vs time for the well with different well-location factors are presented. By comparing numerical results with the analytical solutions of the well located in the center of the closed circle reservoir, the numerical method is verified.

Long-Range Stress Redistribution Resulting From Damage in Heterogeneous Media

It has been shown in CA simulations and data analysis of earthquakes that declustered or characteristic large earthquakes may occur with long-range stress redistribution. In order to understand long-range stress redistribution, we propose a linear-elastic but heterogeneous-brittle model. The stress redistribution in the heterogeneous-brittle medium implies a longer-range interaction than that in an elastic medium. Therefore, it is surmised that the longer-range stress redistribution resulting from damage in heterogeneous media may be a plausible mechanism governing main shocks.

Calculations of longitudinal and transverse velocity structure functions using a vortex model of isotropic turbulence

The longitudinal structure function (LSF) and the transverse structure function (TSF) in isotropic turbulence are calculated using a vortex model. The vortex model is composed of the Rankine and Burgers vortices which have the exponential distributions in the vortex Reynolds number and vortex radii. This model exhibits a power law in the inertial range and satisfies the minimal condition of isotropy that the second-order exponent of the LSF in the inertial range is equal to that of the TSF. Also observed are differences between longitudinal and transverse structure functions caused by intermittency. These differences are related to their scaling differences which have been previously observed in experiments and numerical simulations.

Numerical Simulation of Particle Separation in An Oil-Sand Separator

The gathering systems of crude oil are greatly endangered by the fine sand and soil in oil. Up to now , how to separate sand from the viscid oil is still a technical problem for oil production home or abroad. Recently , Institute of Mechanics in Chinese Academy of Sciences has developed a new type of oil-sand separator , which has been applied successfully in oil field in situ. In this paper, the numerical method of vortex-stream function is used to predict the liquid-solid separating course and the efficiency for this oil-sand separator. Results show that the viscosity and particle diameter have much influence on the particle motion. The calculating separating efficiency is compared with that of experiment and indicates that this method can be used to model the complex two-phase flow in the separator.

Direct Numerical Simulation of Supersonic Turbulent Boundary Layer Flow

Direct numerical simulations of a spatially evolving supersonic flat-plate turbulent boundary layer flow with free Mach number M = 2.25 and Reynolds number Re = 365000/in are performed. The transition process from laminar to turbulent flow is obtained by solving the three-dimensional compressible Navier-Stokes, equations, using high-order accurate difference schemes. The obtained statistical results agree well with the experimental and theoretical data. From the numerical results it can be seen that the transition process under the considered conditions is the process which skips the Tolimien-Schlichting instability and the second instability through the instability of high gradient shear layer and becomes of laminar flow breakdown. This means that the transition process is a bypass-type transition process. The spanwise asymmetry of the disturbance locally upstream imposed is important to induce the bypass-type transition. Furthermore, with increasing the time disturbance frequency the transition will delay. When the time disturbance frequency is large enough, the transition will disappear.

Governing Equations and Numerical Solutions of Tension Leg Platform with Finite Amplitude Motion

It is demonstrated that when tension leg platform (TLP) moves with finite amplitude in waves, the inertia force, the drag force and the buoyancy acting on the platform are nonlinear functions of the response of TLP. The tensions of the tethers are also nonlinear functions of the displacement of TLP. Then the displacement, the velocity and the acceleration of TLP should be taken into account when loads are calculated. In addition, equations of motions should be set up on the instantaneous position. A theoretical model for analyzing the nonlinear behavior of a TLP with finite displacement is developed, in which multifold nonlinearities are taken into account, i.e., finite displacement, coupling of the six degrees of freedom, instantaneous position, instantaneous wet surface, free surface effects and viscous drag force. Based on the theoretical model, the comprehensive nonlinear differential equations are deduced. Then the nonlinear dynamic analysis of ISSC TLP in regular waves is performed in the time domain. The degenerative linear solution of the proposed nonlinear model is verified with existing published one. Furthermore, numerical results are presented, which illustrate that nonlinearities exert a significant influence on the dynamic responses of the TLP.

Experimental Investigation and Numerical Simulation on the Effect of Fissure Water Pressure in Vertical Sliding Surface

This paper studies the effect of fissure water pressure in different fractures on the critical angle of landslide by laboratory investigation and numerical simulation in order to understand the mechanisms of fissure water pressure on landslide stability. Laboratory observations show that the effect of fissure water pressure on the critical angle of landslide is little when the distance between water-holding fracture and slope toe is three times greater than the depth of fissure water. These experimental results are also simulated by a three-dimensional face-to-face contact discrete element method. This method has included the fissure water pressure and can accurately calculate the critical angle of jointed slope when fissure water pressure in vertical sliding surface exists. Numerical results are in good agreement with experimental observations. It is revealed that the location of water-holding structural surface is important to landslide stability. The ratio of the distance between water-holding fissure and slope toe to the depth of fissure water is a key parameter to justify the effect of fissure water pressure on the critical angle of landslide.

Switching performance of epitaxially grown normally-off 4H-SiC JFET

Static and dynamic behavior of the epitaxially grown dual gate trench 4H-SiC junction field effect transistor (JFET) is investigated. Typical on-state resistance Ron was 6-10mΩcm2 at VGS = 2.5V and the breakdown voltage between the range of 1.5-1.8kV was realized at VGS = -5V for normally-off like JFETs. It was found that the turn-on energy delivers the biggest part of the switching losses. The dependence of switching losses from gate resistor is nearly linear, suggesting that changing the gate resistor, a way similar to Si-IGBT technology, can easily control di/dt and dv/dt. Turn-on losses at 200°C are lower compared to those at 25°C, which indicates the influence of the high internal p-type gate layer resistance. Inductive switching numerical analysis suggested the strong influence of channel doping conditions on the turn-on switching performance. The fast switching normally-off JFET devices require heavily doped narrow JFET channel design. © (2009) Trans Tech Publications, Switzerland.

Weak Zone Related Seismic Cycles in Progressive Failure Leading to Collapse in Brittle Crust

Until quite recently our understanding of the basic mechanical process responsible for earthquakes and faulting was not well known. It can be argued that this was partly a consequence of the complex nature of fracture in crust and in part because evidence of brittle phenomena in the natural laboratory of the earth is often obliterated or obscured by other geological processes. While it is well understood that the spatial and temporal complexity of earthquakes and the fault structures emerge from geometrical and material built-in heterogeneities, one important open question is how the shearing becomes localized into a band of intense fractures. Here the authors address these questions through a numerical approach of a tectonic plate by considering rockmass heterogeneity both in microscopic scale and in mesoscopic scale. Numerical simulations of the progressive failure leading to collapse under long-range slow driving forces in the far-field show earthquake-like rupture behavior. $En Echelon$ crack-arrays are reproduced in the numerical simulation. It is demonstrated that the underlying fracturing induced acoustic emissions (or seismic events) display self-organized criticality------from disorder to order. The seismic cycles and the geometric structures of the fracture faces, which are found greatly depending on the material heterogeneity (especially on the macroscopic scale), agree with that observed experimentally in real brittle materials. It is concluded that in order to predict a main shock, one must have extremely detailed knowledge on very minor features of the earth's crust far from the place where the earthquake originated. If correct, the model proposed here seemingly provides an explanation as to why earthquakes to date are not predicted so successfully. The reason is not that the authors do not understand earthquake mechanisms very well but that they still know little about our earth's crust.

A Numerical Study for Turbulent Flow and Thermal Influence Over Inhomogenous Canopy of Roughness Elements

A large-eddy simulation with transitional structure function(TSF) subgrid model we previously proposed was performed to investigate the turbulent flow with thermal influence over an inhomogeneous canopy, which was represented as alternative large and small roughness elements. The aerodynamic and thermodynamic effects of the presence of a layer of large roughness elements were modelled by adding a drag term to the three-dimensional Navier-Stokes equations and a heat source/sink term to the scalar equation, respectively. The layer of small roughness elements was simply treated using the method as described in paper (Moeng 1984, J. Atmos Sci. 41, 2052-2062) for homogeneous rough surface. The horizontally averaged statistics such as mean vertical profiles of wind velocity, air temperature, et al., are in reasonable agreement with Gao et al.(1989, Boundary layer meteorol. 47, 349-377) field observation (homogeneous canopy). Not surprisingly, the calculated instantaneous velocity and temperature fields show that the roughness elements considerably changed the turbulent structure within the canopy. The adjustment of the mean vertical profiles of velocity and temperature was studied, which was found qualitatively comparable with Belcher et al. (2003, J Fluid Mech. 488, 369-398)'s theoretical results. The urban heat island(UHI) was investigated imposing heat source in the region of large roughness elements. An elevated inversion layer, a phenomenon often observed in the urban area (Sang et al., J Wind Eng. Ind. Aesodyn. 87, 243-258)'s was successfully simulated above the canopy. The cool island(CI) was also investigated imposing heat sink to simply model the evaporation of plant canopy. An inversion layer was found very stable and robust within the canopy.

Effects of surface roughness on airframe noise

The generation of sound by turbulent boundary layer flow at low Mach number over a rough wall is investigated by applying the theoretical model which describes the scattering of the turbulence near field into sound by roughness elements. Attention is focused on the numerical method to approximately quantify the absolute level of the roughness noise radiated to far field. Empirical models for the source statistics are obtained by scaling smooth-wall data through increased skin friction velocity and boundary layer thickness for the rough surface. Numerical integration is performed to determine the roughness noise, and it reproduces the spectral characteristics of the available empirical formula and experimental data. Experiments are conducted to measure the radiated sound from two rough plates in an open jet by four 1/2'' free-field condenser microphones. The measured noise spectra of the rough plates are above that of a smooth plate in 1-2.5 kHz frequency and exhibits encouraging agreement with the predicted spectra. Also, a phased microphone array is utilized to localize the sound source, and it confirms that the rough plates generate higher source strengthes in this frequency range. A parametric study illustrates that the roughness height and roughness density significantly affect the far-field radiated roughness noise with the roughness height having the dominant effect. The estimates of the roughness noise for a Boeing 757 sized aircraft wing show that in high frequency region the sound radiated from surface roughness may exceed that from the trailing edge, and higher overall sound pressure levels for the roughness noise are also observed.

Numerical Investigation on Laser Transformation Hardening with Different Temporal Pulse Shapes

A 3-D numerical model for pulsed laser transformation hardening (LTH) is developed using the finite element method. In this model, laser spatial and temporal intensity distribution, temperature-dependent thermophysical properties of material, and multi-phase transformations are considered. The influence of laser temporal pulse shape on connectivity of hardened zone, maximum surface temperature of material and hardening depth is numerically investigated at different pulse energy levels. Results indicate that these hardening parameters are strongly dependent on the temporal pulse shape. For the rectangular temporal pulse shape, the temperature field obtained from this model is in excellent agreement with analytical solution, and the predicted hardening depth is favorably compared with experimental one. It should be pointed out that appropriate temporal pulse shape should be selected according to pulse energy level in order to achieve desirable hardening quality under certain laser spatial intensity distribution.