A simulation model of a floating half zone

A simulation model of a floating half zone was suggested by steady numerical simulation and experiment respectively, in the previous papers [Q.S. Chen, W.R. Hu, Int. J. Mass Heat Transfer 40 (1997) 757; J.H. Han, Y. Ar, R. Zhou, W.R. Hue, Int. J. Mass Heat Transfer 40 (1997) 2671]. In the present paper, the simulation model is studied by using the method of unsteady and three-dimensional numerical simulation, and the transient process from steady convection to oscillatory convection is especially analyzed. Comparison of onsets of oscillation for both simulation model and the usual model were obtained, and the results show that the critical Marangoni number of the simulation model is obviously smaller than that of the usual model for the same slender liquid bridge. This implies that the usual model of a floating half zone gives a lower estimation on the onset of oscillation for floating zone convection.

Two-dimensional Cartesian air-gap element (CAGE) for dynamic finite-element modeling of electrical machines with a flat air gap

This paper extends the air-gap element (AGE) to enable the modeling of flat air gaps. AGE is a macroelement originally proposed by Abdel-Razek et al.for modeling annular air gaps in electrical machines. The paper presents the theory of the new macroelement and explains its implementation within a time-stepped finite-element (FE) code. It validates the solution produced by the new macroelement by comparing it with that obtained by using an FE mesh with a discretized air gap. It then applies the model to determine the open-circuit electromotive force of an axial-flux permanent-magnet machine and compares the results with measurements.

Seismic Responses of a Hemispherical Alluvial Valley to SV Waves: A Three-Dimensional Analytical Approximation

An analytical solution to the three-dimensional scattering and diffraction of plane SV-waves by a saturated hemispherical alluvial valley in elastic half-space is obtained by using Fourier-Bessel series expansion technique. The hemispherical alluvial valley with saturated soil deposits is simulated with Biot's dynamic theory for saturated porous media. The following conclusions based on numerical results can be drawn: (1) there are a significant differences in the seismic response simulation between the previous single-phase models and the present two-phase model; (2) the normalized displacements on the free surface of the alluvial valley depend mainly on the incident wave angles, the dimensionless frequency of the incident SV waves and the porosity of sediments; (3) with the increase of the incident angle, the displacement distributions become more complicated; and the displacements on the free surface of the alluvial valley increase as the porosity of sediments increases.

High Aspect Ratio (L/D) Riser Viv Prediction Using Wake Oscillator Model

A two-dimensional (2-D) vortex-induced vibration (VIV) prediction model for high aspect ratio (LID) riser subjected to uniform and sheared flow is studied in this paper. The nonlinear structure equations are considered. The near wake dynamics describing the fluctuating nature of vortex shedding is modeled using classical van der Pol equation. A new approach was applied to calibrate the empirical parameters in the wake oscillator model. Compared the predicted results with the experimental data and computational fluid dynamic (CFD) results. Good agreements are observed. It can be concluded that the present model can be used as simple computational tool in predicting some aspects of VIV of long flexible structures. (C) 2008 Elsevier Ltd. All rights reserved.

Pull-In Instability Of Micro-Switch Actuators: Model Review

The existing three widely used pull-in theoretical models (i.e., one-dimensional lumped model, linear supposition model and planar model) are compared with the nonlinear beam mode in this paper by considering both cantilever and fixed-fixed type micro and nano-switches. It is found that the error of the pull-in parameters between one-dimensional lumped model and the nonlinear beam model is large because the denominator of the electrostatic force is minimal when the electrostatic force is computed at the maximum deflection along the beam. Since both the linear superposition model and the slender planar model consider the variation of electrostatic force with the beam's deflection, these two models not only are of the same type but also own little error of the pull-in parameters with the nonlinear beam model, the error brought by these two models attributes to that the boundary conditions are not completely satisfied when computing the numerical integration of the deflection.

Two-Dimensional Kinetics Of Beta(2)-Integrin And Icam-1 Bindings Between Neutrophils And Melanoma Cells In A Shear Flow

Cell adhesion, mediated by specific receptor-ligand interactions, plays an important role in biological processes such as tumor metastasis and inflammatory cascade. For example, interactions between beta(2)-integrin ( lymphocyte function-associated antigen-1 and/or Mac-1) on polymorphonuclear neutrophils (PMNs) and ICAM-1 on melanoma cells initiate the bindings of melanoma cells to PMNs within the tumor microenvironment in blood flow, which in turn activate PMN-melanoma cell aggregation in a near-wall region of the vascular endothelium, therefore enhancing subsequent extravasation of melanoma cells in the microcirculations. Kinetics of integrin-ligand bindings in a shear flow is the determinant of such a process, which has not been well understood. In the present study, interactions of PMNs with WM9 melanoma cells were investigated to quantify the kinetics of beta(2)-integrin and ICAM-1 bindings using a cone-plate viscometer that generates a linear shear flow combined with a two-color flow cytometry technique. Aggregation fractions exhibited a transition phase where it first increased before 60 s and then decreased with shear durations. Melanoma-PMN aggregation was also found to be inversely correlated with the shear rate. A previously developed probabilistic model was modified to predict the time dependence of aggregation fractions at different shear rates and medium viscosities. Kinetic parameters of beta(2)-integrin and ICAM-1 bindings were obtained by individual or global fittings, which were comparable to respectively published values. These findings provide new quantitative understanding of the biophysical basis of leukocyte-tumor cell interactions mediated by specific receptor-ligand interactions under shear flow conditions.

Effective elastic moduli of inhomogeneous solids by embedded cell model

An embedded cell model is presented to obtain the effective elastic moduli for three-dimensional two-phase composites which is an exact analytic formula without any simplified approximation and can be expressed in an explicit form. For the different cells such as spherical inclusions and cracks surrounded by sphere and oblate ellipsoidal matrix, the effective elastic moduli are evaluated and the results are compared with those from various micromechanics models. These results show that the present model is direct, simple and efficient to deal with three-dimensional tyro-phase composites.

A non-equilibrium sediment transport model for rill erosion

Sediment transport in rill flows exhibits the characteristics of non-equilibrium transport, and the sediment transport rate of rill flow gradually recovers along the flow direction by erosion. By employing the concept of partial equilibrium sediment transport from open channel hydraulics, a dynamic model of rill erosion on hillslopes was developed. In the model, a parameter, called the restoration coefficient of sediment transport capacity, was used to express the recovery process of sediment transport rate, which was analysed by dimensional analysis and determined from laboratory experimental data. The values of soil loss simulated by the model were in agreement with observed values. The model results showed that the length and gradient of the hillslope and rainfall intensity had different influences on rill erosion. Copyright (c) 2006 John Wiley & Sons, Ltd.

Three-dimensional crack problem analysis using boundary element method with finite-part integrals

A set of hypersingular integral equations of a three-dimensional finite elastic solid with an embedded planar crack subjected to arbitrary loads is derived. Then a new numerical method for these equations is proposed by using the boundary element method combined with the finite-part integral method. According to the analytical theory of the hypersingular integral equations of planar crack problems, the square root models of the displacement discontinuities in elements near the crack front are applied, and thus the stress intensity factors can be directly calculated from these. Finally, the stress intensity factor solutions to several typical planar crack problems in a finite body are evaluated.

An improvement on the theoretical model for the large-scale coherent structures in the outer region of a turbulent boundary layer

According to the experimental results, there exist large-scale coherent structures in the outer region of a turbulent boundary layer, which have been studied by many authors.As experimental results, Antonia (1990) showed the phase- aver aged streamlines and isovorticity lines of the large-scale coherent structures in a turbulent boundary layer for different Reynolds numbers. Based on the hydrodynamic stability theory, the 2-D theoretical model for the large-scale structures was proposed by Luo and Zhou, in which the eddy viscosity was defined as a complex function of the position in the normal direction. The theoretical results showed in ref. were in agreement with those in ref. However, there were two problems in the results. One is that in the experimental results, there were divergent focuses between two saddle points in the streamlines, but in the theoretical results, there were centers. The other is that the stretched parts of the isovorticity lines appear at the location of centers in the theoretical results, while in the experimental results they located somewhere between the focuses and saddle points. The reason is that the computations were based on a 2-D model.

Phonon properties of one-dimensional nanocrystalline solids

We study phonon properties of one-dimensional nanocrystalline solids that are associated with a model nanostructured sequence. A real-space renormalization-group approach, connected with a series of renormalization-group transformations, is developed to calculate numerically the local phonon Green's function at an arbitrary site, and then the phonon density of states of these kinds of nanocrystalline chains. Some interesting phonon properties of nanocrystalline chains are obtained that are in qualitative agreement with the experimental results for the optical-absorption spectra of nanostructured solids.

Two-dimensional simulation of inductive plasma sources with self-consistent power deposition

A time averaged two-dimensional fluid model including an electromagnetic module with self-consistent power deposition was developed to simulate the transport of a low pressure radio frequency inductively coupled plasma source. Comparsions with experiment and previous simulation results show, that the fluid model is feasible in a certain range of gas pressure. In addition, the effects of gas pressure and power input have been discussed.

The thermal and mechanical structure of a two-dimensional plume in the earths mantle

On the condition that the distribution of velocity and temperature at the mid-plane of a mantle plume has been obtained (pages 213–218, this issue), the problem of determining the lateral structure of the plume at a given depth is reduced to solving an eigenvalue problem of a set of ordinary differential equations with five unknown functions, with an eigenvalue being related to the thermal thickness of the plume at this depth. The lateral profiles of upward velocity, temperature and viscosity in the plume and the thickness of the plume at various depths are calculated for two sets of Newtonian rheological parameters. The calculations show that the precondition for the existence of the plume, δT/L 1 (L = the height of the plume, δT = lateral distance from the mid-plane), can be satisfied, except for the starting region of the plume or near the base of the lithosphere. At the lateral distance, δT, the upward velocity decreases to 0.1 – 50% of its maximum value at different depths. It is believed that this model may provide an approach for a quantitative description of the detailed structure of a mantle plume.

The wave model of mesothermal plasma near wakes and Korteweg-de Vries equation

The stationary two-dimensional (x, z) near wakes behind a flat-based projectile which moves at a constant mesothermal speed (V∞) along a z-axis in a rarefied, fully ionized, plasma is studied using the wave model previously proposed by one of the authors (VCL). One-fluid theory is used to depict the free expansion of ambient plasma into the vacuum produced behind a fast-moving projectile. This nonstationary, one-dimensional (x, t) flow which is approximated by the K-dV equation can be transformed, through substitution, t=z/V∞, into a stationary two-dimensional (x, z) near wake flow seen by an observer moving with the body velocity (V∞). The initial value problem of the K-dV equation in (x, t) variables is solved by a specially devised numerical method. Comparisons of the present numerical solution for the asymptotically small and large times with available analytical solutions are made and found in satisfactory agreements.

Starting Nozzle Flow Simulation Using K-G Two-Equation Turbulence Model

The starting process of two-dimensional nozzle flows has been simulated with Euler, laminar and k - g two-equation turbulence Navier-Stokes equations. The flow solver is based on a combination of LUSGS subiteration implicit method and five spatial discretized schemes, which are Roe, HLLE, MHLLE upwind schemes and AUSM+, AUSMPW schemes. In the paper, special attention is for the flow differences of the nozzle starting process obtained from different governing equations and different schemes. Two nozzle flows, previously investigated experimentally and numerically by other researchers, are chosen as our examples. The calculated results indicate the carbuncle phenomenon and unphysical oscillations appear more or less near a wall or behind strong shock wave except using HLLE scheme, and these unphysical phenomena become more seriously with the increase of Mach number. Comparing the turbulence calculation, inviscid solution cannot simulate the wall flow separation and the laminar solution shows some different flow characteristics in the regions of flow separation and near wall.