Energy conversion in shape memory alloy heat engine - Part II: Simulation

In this paper, a theoretical model proposed in Part I (Zhu et al., 2001a) is used to simulate the behavior of a twin crank NiTi SMA spring based heat engine, which has been experimentally studied by Iwanaga et al. (1988). The simulation results are compared favorably with the measurements. It is found that (1) output torque and heat efficiency decrease as rotation speed increase; (2) both output torque and output power increase with the increase of hot water temperature; (3) at high rotation speed, higher water temperature improves the heat efficiency. On the contrary, at low rotation speed, lower water temperature is more efficient; (4) the effects of initial spring length may not be monotonic as reported. According to the simulation, output torque, output power and heat efficiency increase with the decrease of spring length only in the low rotation speed case. At high rotation speed, the result might be on the contrary.

Elasticity, stability, and ideal strength of beta-SiC in plane-wave-based ab initio calculations

On the basis of the pseudopotential plane-wave method and the local-density-functional theory, this paper studies energetics, stress-strain relation, stability, and ideal strength of beta-SiC under various loading modes, where uniform uniaxial extension and tension and biaxial proportional extension are considered along directions [001] and [111]. The lattice constant, elastic constants, and moduli of equilibrium state are calculated and the results agree well with the experimental data. As the four SI-C bonds along directions [111], [(1) over bar 11], [11(1) over bar] and [111] are not the same under the loading along [111], internal relaxation and the corresponding internal displacements must be considered. We find that, at the beginning of loading, the effect of internal displacement through the shuffle and glide plane diminishes the difference among the four Si-C bonds lengths, but will increase the difference at the subsequent loading, which will result in a crack nucleated on the {111} shuffle plane and a subsequently cleavage fracture. Thus the corresponding theoretical strength is 50.8 GPa, which agrees well with the recent experiment value, 53.4 GPa. However, with the loading along [001], internal relaxation is not important for tetragonal symmetry. Elastic constants during the uniaxial tension along [001] are calculated. Based on the stability analysis with stiffness coefficients, we find that the spinodal and Born instabilities are triggered almost at the same strain, which agrees with the previous molecular-dynamics simulation. During biaxial proportional extension, stress and strength vary proportionally with the biaxial loading ratio at the same longitudinal strain.

Formation of Fivefold Deformation Twins in Nanocrystalline Face-Centered-Cubic Copper Based on Molecular Dynamics Simulations

Fivefold deformation twins were reported recently to be observed in the experiment of the nanocrystalline face-centered-cubic metals and alloys. However, they were not predicted previously based on the molecular dynamics (MD) simulations and the reason was thought to be a uniaxial tension considered in the simulations. In the present investigation, through introducing pretwins in grain regions, using the MD simulations, the authors predict out the fivefold deformation twins in the grain regions of the nanocrystal grain cell, which undergoes a uniaxial tension. It is shown in their simulation results that series of Shockley partial dislocations emitted from grain boundaries provide sequential twining mechanism, which results in fivefold deformation twins. (c) 2006 American Institute of Physics.

Rossby Vortex Simulation on a Paraboloidal Coordinate System Using the Lattice Boltzmann Method

In this paper, we apply our compressible lattice Boltzmann model to a rotating parabolic coordinate system to simulate Rossby vortices emerging in a layer of shallow water flowing zonally in a rotating paraboloidal vessel. By introducing a scaling factor, nonuniform curvilinear mesh can be mapped to a flat uniform mesh and then normal lattice Boltzmann method works. Since the mass per unit area on the two-dimensional (2D) surface varies with the thickness of the water layer, the 2D flow seems to be "compressible" and our compressible model is applied. Simulation solutions meet with the experimental observations qualitatively. Based on this research, quantitative solutions and many natural phenomena simulations in planetary atmospheres, oceans, and magnetized plasma, such as the famous Jovian Giant Red Spot, the Galactic Spiral-vortex, the Gulf Stream, and the Kuroshio Current, etc,, can be expected.

Large-Eddy-Simulation of 3-Dimensional Rayleigh-Taylor Instability in Incompressible Fluids

The 3-dimensiqnal incompressible Rayleigh-Taylor instability is numerically studied through the large-eddy-simulation (LES) approach based on the passive scalar transport model. Both the instantaneous velocity and the passive scalar fields excited by sinu

Simulation of Free-Surface Flow in a Tank Using the Navier-Stokes Model and Unstructured Finite Volume Method

Modelling free-surface flow has very important applications in many engineering areas such as oil transportation and offshore structures. Current research focuses on the modelling of free surface flow in a tank by solving the Navier-Stokes equation. An unstructured finite volume method is used to discretize the governing equations. The free surface is tracked by dynamically adapting the mesh and making it always surface conforming. A mesh-smoothing scheme based on the spring analogy is also implemented to ensure mesh quality throughout the computaiton. Studies are performed on the sloshing response of a liquid in an elastic container subjected to various excitation frequencies. Further investigations are also carried out on the critical frequency that leads to large deformation of the tank walls. Another numerical simulation involves the free-surface flow past as submerged obstacle placed in the tank to show the flow separation and vortices. All these cases demonstrate the capability of this numerical method in modelling complicated practical problems.

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.

Correlative reference model and molecular dynamics simulation of dislocation emission process

A correlative reference model for computer molecular dynamics simulations is proposed. Based on this model, a flexible displacement boundary scheme is introduced and the dislocations emitted from a crack tip can continuously pass through the border of the inner discrete atomic region and pile up at the outer continuum region. The effect of the emitted dislocations within the plastic zone on the inner atomistic region can be clearly demonstrated. The simulations for a molybdinum crystal show that a full dislocation in a bcc crystal is dissociated into three partial dislocations and interaction between the crack and the emitted dislocations results in gradual decrease of the local stress intensity factor.

Molecular dynamics simulation of hydrogen enhancing dislocation emission

The interactive pair potential between Al and H is obtained based on the ab initio calculation and the Chen-Mobius 3D lattice inversion formula. By utilizing the pair potentials calculated, the effects of hydrogen on the dislocation emission from crack tip have been studied. The simulated result shows that hydrogen can reduce the cohesive strength for Al single crystal, and then the critical stress intensity factor for partial dislocation emission decreases from 0.11 MPa root m (C-H = 0) to 0.075 MPa root m (C-H=0.72%) and 0.06 MPa root m (C-H = 1.44%). This indicates thar hydrogen can enhance the dislocation emission. The simulation also shows that atoms of hydrogen can gather and turn into small bubbles, resulting in enhancement of the equilibrium vacancy concentration.

Sensitivity Analysis of Dimensionless Parameters for Physical Simulation of Water-Flooding Reservoir

A numerical approach to optimize dimensionless parameters of water-flooding porous media flows is proposed based on the analysis of the sensitivity factor defined as the variation ration of a target function with respect to the variation of dimensionless parameters. A complete set of scaling criteria for water-flooding reservoir of five-spot well pattern case is derived from the 3-D governing equations, involving the gravitational force, the capillary force and the compressibility of water, oil and rock. By using this approach, we have estimated the influences of each dimensionless parameter on experimental results and thus sorted out the dominant ones with larger sensitivity factors ranging from10-4to10-0 .

Large Eddy Simulation of Complex Turbulent Flows: Physical Aspects and Research Trends

In the current paper, we have primarily addressed one powerful simulation tool developed during the last decades-Large Eddy Simulation (LES), which is most suitable for unsteady three-dimensional complex turbulent flows in industry and natural environment. The main point in LES is that the large-scale motion is resolved while the small-scale motion is modeled or, in geophysical terminology, parameterized. With a view to devising a subgrid-scale(SGS) model of high quality, we have highlighted analyzing physical aspects in scale interaction and-energy transfer such as dissipation, backscatter, local and non-local interaction, anisotropy and resolution requirement. They are the factors responsible for where the advantages and disadvantages in existing SGS models come from. A case study on LES of turbulence in vegetative canopy is presented to illustrate that LES model is more based on physical arguments. Then, varieties of challenging complex turbulent flows in both industry and geophysical fields in the near future-are presented. In conclusion; we may say with confidence that new century shall see the flourish in the research of turbulence with the aid of LES combined with other approaches.

Simulation of the Load-Unload Response Ratio and Critical Sensitivity in the Lattice Solid Model

The Load-Unload Response Ratio (LURR) method is an intermediate-term earthquake prediction approach that has shown considerable promise. It involves calculating the ratio of a specified energy release measure during loading and unloading where loading and unloading periods are determined from the earth tide induced perturbations in the Coulomb Failure Stress on optimally oriented faults. In the lead-up to large earthquakes, high LURR values are frequently observed a few months or years prior to the event. These signals may have a similar origin to the observed accelerating seismic moment release (AMR) prior to many large earthquakes or may be due to critical sensitivity of the crust when a large earthquake is imminent. As a first step towards studying the underlying physical mechanism for the LURR observations, numerical studies are conducted using the particle based lattice solid model (LSM) to determine whether LURR observations can be reproduced. The model is initialized as a heterogeneous 2-D block made up of random-sized particles bonded by elastic-brittle links. The system is subjected to uniaxial compression from rigid driving plates on the upper and lower edges of the model. Experiments are conducted using both strain and stress control to load the plates. A sinusoidal stress perturbation is added to the gradual compressional loading to simulate loading and unloading cycles and LURR is calculated. The results reproduce signals similar to those observed in earthquake prediction practice with a high LURR value followed by a sudden drop prior to macroscopic failure of the sample. The results suggest that LURR provides a good predictor for catastrophic failure in elastic-brittle systems and motivate further research to study the underlying physical mechanisms and statistical properties of high LURR values. The results provide encouragement for earthquake prediction research and the use of advanced simulation models to probe the physics of earthquakes.

The Mach reflection of a detonation based on soot track measurements

This paper presents a series of soot tracks formed by gaseous detonation waves diffracting around wedges with different wedge angles. These cellular structure patterns describe the Mach-reflection processes of a detonation and reveal some unique characteristics. They can be used to analyze the relationship between the trajectory angle of the triple point, wedge angle, and initial pressure in Mach reflection. Compared to the Mach-reflected one-dimensional shock wave in nonreactive air, all these unique characteristics for a Mach-reflected detonation should be attributed to the transverse-wave structure of the detonation front; meanwhile, the precursor shock wave and transverse wave influence the Mach-reflected detonation, respectively. The experimental results support the recently published numerical simulation of this complex phenomenon.

Molecular dynamics simulation on dislocation emission process

A correlative reference model for a computer simulation of molecular dynamics is proposed in this paper. Based on this model, a flexible displacement boundary scheme is naturally introduced and the dislocations emitted from a crack tip are presumed to continuously pass through the border of an inner discrete atomic region to pile up at an outer continuum region. The simulations for a Mo crystal show that the interaction between a crack and emitted dislocations results in the decrease in local stress intensity factor gradually.

3D thermo-electro-mechanical simulations of gas sensors based on SOI membranes

3D thermo-electro-mechanical device simulations are presented of a novel fully CMOS-compatible MOSFET gas sensor operating in a SOI membrane. A comprehensive stress analysis of a Si-SiO2-based multilayer membrane has been performed to ensure a high degree of mechanical reliability at a high operating temperature (e.g. up to 400°C). Moreover, optimisation of the layout dimensions of the SOI membrane, in particular the aspect ratio between the membrane length and membrane thickness, has been carried out to find the best trade-off between minimal device power consumption and acceptable mechanical stress.