The application of B-P constitutive equations in finite element analysis of high velocity impact

We present in this paper the application of B-P constitutive equations in finite element analysis of high velocity impact. The impact process carries out in so quick time that the heat-conducting can be neglected and meanwhile, the functions of temperature in equations need to be replaced by functions of plastic work. The material constants in the revised equations can be determined by comparison of the one-dimensional calculations with the experiments of Hopkinson bar. It can be seen from the comparison of the calculation with the experiment of a tungsten alloy projectile impacting a three-layer plate that the B-P constitutive equations in that the functions of temperature were replaced by the functions of plastic work can be used to analysis of high velocity impact.

Choice of calibration equations in calibrating microwave test fixtures

The error theory of linear equation system has been applied to the calibration procedure of microwave network analyser in this article. A new explanation for the choice of the linear calibration equations is proposed and a general principle for choosing calibration equations is presented. The method can also be used to predict the occurrence of the problem of frequency limitation at some periodic frequencies. This principle is employed to the thru-short-delay (TSD) method and the solution using the chosen equations gives the most accurate results. A good agreement between the theory and the experiment has been obtained.

H4 exponential finite-difference scheme for convective diffusion-equations

Perturbations are applied to the convective coefficients and source term of a convection-diffusion equation so that second-order corrections may be applied to a second-order exponential scheme. The basic Structure of the equations in the resulting fourth-order scheme is identical to that for the second order. Furthermore, the calculations are quite simple as the second-order corrections may be obtained in a single pass using a second-order scheme. For one to three dimensions, the fourth-order exponential scheme is unconditionally stable. As examples, the method is applied to Burgers' and other fluid mechanics problems. Compared with schemes normally used, the accuracies are found to be good and the method is applicable to regions with large gradients.

Higher-order Boussinesq-type equations for interfacial waves in a two-fluid system

Interfacial waves propagating along the interface between a three-dimensional two-fluid system with a rigid upper boundary and an uneven bottom are considered. There is a light fluid layer overlying a heavier one in the system, and a small density difference exists between the two layers. A set of higher-order Boussinesq-type equations in terms of the depth-averaged velocities accounting for stronger nonlinearity are derived. When the small parameter measuring frequency dispersion keeping up to lower-order and full nonlinearity are considered, the equations include the Choi and Camassa's results (1999). The enhanced equations in terms of the depth-averaged velocities are obtained by applying the enhancement technique introduced by Madsen et al. (1991) and Schaffer and Madsen (1995a). It is noted that the equations derived from the present study include, as special cases, those obtained by Madsen and Schaffer (1998). By comparison with the dispersion relation of the linear Stokes waves, we found that the dispersion relation is more improved than Choi and Camassa's (1999) results, and the applicable scope of water depth is deeper.

Conjugate Model for Heat and Mass Transfer of Porous Wall in the High Temperature Gas Flow

Heat and mass transfer of a porous permeable wall in a high temperature gas dynamical flow is considered. Numerical simulation is conducted on the ground of the conjugate mathematical model which includes filtration and heat transfer equations in a porous body and boundary layer equations on its surface. Such an approach enables one to take into account complex interaction between heat and mass transfer in the gasdynamical flow and in the structure subjected to this flow. The main attention is given to the impact of the intraporous heat transfer intensity on the transpiration cooling efficiency.

Flood Routing Models in Confluent and Dividing Channels

By introducing a water depth connecting formula, the hydraulic equations in the dividing channel system were coupled and the relation of discharge distribution between the branches of the dividing channels can be yielded. In this manner, a numerical model for the confluent channels was established to study the variation of backwater effects with the parameters in the channel junction. The meeting of flood peaks in the mainstream and tributary can be analyzed with this model. The flood peak meeting is found to be a major factor for the extremely high water level in the mainstream during the 1998 Yangtze River flood. Subsequently the variations of discharge distribution and water level with channel parameters between each branch in this system were studied as well. As a result, flood evolution caused by Jingjiang River shortcut and sediment deposition in the entrance of dividing channels of the Yangtze River may be qualitatively elucidated. It is suggested to be an effective measure for flood mitigation to enhance regulation capability of reservoirs available upstream of the tributaries and harness branch entrance channels.

A Constrained Particle Dynamics for Continuum-Particle Hybrid Method in Micro- and Nano-Fluidics

A hybrid method of continuum and particle dynamics is developed for micro- and nano-fluidics, where fluids are described by a molecular dynamics (MD) in one domain and by the Navier-Stokes (NS) equations in another domain. In order to ensure the continuity of momentum flux, the continuum and molecular dynamics in the overlap domain are coupled through a constrained particle dynamics. The constrained particle dynamics is constructed with a virtual damping force and a virtual added mass force. The sudden-start Couette flows with either non-slip or slip boundary condition are used to test the hybrid method. It is shown that the results obtained are quantitatively in agreement with the analytical solutions under the non-slip boundary conditions and the full MD simulations under the slip boundary conditions.

A Cip/Multi-Moment Finite Volume Method For Shallow Water Equations With Source Terms

A novel finite volume method has been presented to solve the shallow water equations. In addition to the volume-integrated average (VIA) for each mesh cell, the surface-integrated average (SIA) is also treated as the model variable and is independently predicted. The numerical reconstruction is conducted based on both the VIA and the SIA. Different approaches are used to update VIA and SIA separately. The SIA is updated by a semi-Lagrangian scheme in terms of the Riemann invariants of the shallow water equations, while the VIA is computed by a flux-based finite volume formulation and is thus exactly conserved. Numerical oscillation can be effectively avoided through the use of a non-oscillatory interpolation function. The numerical formulations for both SIA and VIA moments maintain exactly the balance between the fluxes and the source terms. 1D and 2D numerical formulations are validated with numerical experiments. Copyright (c) 2007 John Wiley & Sons, Ltd.

The mixed-type integral equations for transient elastodynamic plane crack problems

In the present paper, by use of the boundary integral equation method and the techniques of Green fundamental solution and singularity analysis, the dynamic infinite plane crack problem is investigated. For the first time, the problem is reduced to solving a system of mixed-typed integral equations in Laplace transform domain. The equations consist of ordinary boundary integral equations along the outer boundary and Cauchy singular integral equations along the crack line. The equations obtained are strictly proved to be equivalent with the dual integral equations obtained by Sih in the special case of dynamic Griffith crack problem. The mixed-type integral equations can be solved by combining the numerical method of singular integral equation with the ordinary boundary element method. Further use the numerical method for Laplace transform, several typical examples are calculated and their dynamic stress intensity factors are obtained. The results show that the method proposed is successful and can be used to solve more complicated problems.

Domain decomposition hybrid method for numerical simulation of bluff body flows:theoretical model and application

The discrete vortex method is not capable of precisely predicting the bluff body flow separation and the fine structure of flow field in the vicinity of the body surface. In order to make a theoretical improvement over the method and to reduce the difficulty in finite-difference solution of N-S equations at high Reynolds number, in the present paper, we suggest a new numerical simulation model and a theoretical method for domain decomposition hybrid combination of finite-difference method and vortex method. Specifically, the full flow. field is decomposed into two domains. In the region of O(R) near the body surface (R is the characteristic dimension of body), we use the finite-difference method to solve the N-S equations and in the exterior domain, we take the Lagrange-Euler vortex method. The connection and coupling conditions for flow in the two domains are established. The specific numerical scheme of this theoretical model is given. As a preliminary application, some numerical simulations for flows at Re=100 and Re-1000 about a circular cylinder are made, and compared with the finite-difference solution of N-S equations for full flow field and experimental results, and the stability of the solution against the change of the interface between the two domains is examined. The results show that the method of the present paper has the advantage of finite-difference solution for N-S equations in precisely predicting the fine structure of flow field, as well as the advantage of vortex method in efficiently computing the global characteristics of the separated flow. It saves computer time and reduces the amount of computation, as compared with pure N-S equation solution. The present method can be used for numerical simulation of bluff body flow at high Reynolds number and would exhibit even greater merit in that case.

Thermophoretic deposition of particles onto cold surfaces of bodies in two-dimensional and axisymmetric flows

The problem of thermophoretic deposition of small particles onto cold surfaces is studied in two-dimensional and axisymmetric flow fields. The particle concentration equation is solved numerically together with the momentum and energy equations in the laminar boundary layer with variable density effect included. It is shown explicitly to what extent the particle concentration and deposition rate at the wall are influenced by the density variation effect for external flow past bodies. The general numerical procedure is given for two-dimensional and axisymmetric cases and is illustrated with examples of thermophoretic deposition of particles in flows past a cold cylinder and a sphere.

Depth dependent elastic strain in ZnO epilayer: combined Rutherford backscattering/channeling and X-ray diffraction

A ZnO layer was grown by metalorganic chemical vapor deposition (MOCVD) on a sapphire (0 0 0 1) substrate. The perpendicular and parallel elastic strain of the ZnO epilayer, e(perpendicular to) = 0.19%, e(parallel to) = -0.29%, respectively, were derived by using the combination of Rutherford backscattering (RBS)/channeling and X-ray diffraction (XRD). The ratio vertical bar e(parallel to)/ e(perpendicular to)vertical bar = 1.5 indicates that ZnO layer is much stiffer in the a-axis direction than in the c-axis direction. By using RBS/C, the depth dependent elastic strain was deduced. The strain is higher at the depth close to the interface and decreases towards the surface. The negative tetragonal distortion was explained by considering the lattice mismatch and thermal mismatch in ZnO thin film. (c) 2004 Elsevier B.V. All rights reserved.

PIXE analysis of Fe content in Fe-implanted GaN film

The diluted magnetic semiconductors (DMSs) were achieved by the ion implantation. Fe+ ions (250 keV) were implanted into n-type GaN at room temperature with doses ranging from 8 X 10(15) cm(-2) to 8 X 10(16) cm(-2) and subsequently rapidly annealed at 800 degrees C for 5 m in N-2 ambient. PIXE was employed to determine the Fe-implanted content. The magnetic property was measured by the Quantum Design MPMS SQUID magnetometer. No secondary phases or clusters are detected within the sensitivity of XRD. Apparent ferromagnetic hysteresis loops measured at 10 K were presented. The relationships between the Fe-implanted content and the ferromagnetic property are discussed. (c) 2006 Elsevier B.V. All rights reserved.

Properties of proton-implanted p-type Si: supports for the models explaining a novel p-n junction in Si

Proton-implanted and annealed p-type Si wafers were investigated by using both transmission electron microscopy and spreading resistivity probe. The novel pn junction [Li et al., Mat. Res. Sec. Symp, Proc. 396 (1996) 745], as obtained by using n-type Si subjected to the process as this work, was not observed in the p-type Si wafers in this work. A drop of superficial resistivity in the sample was found and is explained by the proposed models interpreting the novel pn junction. (C) 2000 Elsevier Science B.V. All rights reserved.