A simulation study of the mixing, atomic flow and velocity profiles of crystalline materials during sliding

Self-contained Non-Equilibrium Molecular Dynamics (NEMD) simulations using Lennard-Jones potentials were performed to identify the origin and mechanisms of atomic scale interfacial behavior between sliding metals. The mixing sequence and velocity profiles were compared via MD simulations for three cases, viz.: sell-mated, similar and hard-softvcrystal pairs. The results showed shear instability, atomic scale mixing, and generation of eddies at the sliding interface. Vorticity at the interface suggests that atomic flow during sliding is similar to fluid flow under Kelvin-Helmholtz instability and this is supported by velocity profiles from the simulations. The initial step-function velocity profile spreads during sliding. However the velocity profile does not change much at later stages of the simulation and it eventually stops spreading. The steady state friction coefficient during simulation was monitored as a function of sliding velocity. Frictional behavior can be explained on the basis of plastic deformation and adiabatic effects. The mixing layer growth kinetics was also investigated.

Response of a turbulent boundary layer to a step change in surface heat flux

Measurements of both the velocity and the temperature field have been made in the thermal layer that grows inside a turbulent boundary layer which is subjected to a small step change in surface heat flux. Upstream of the step, the wall heat flux is zero and the velocity boundary layer is nearly self-preserving. The thermal-layer measurements are discussed in the context of a self-preserving analysis for the temperature disturbance which grows underneath a thick external turbulent boundary layer. A logarithmic mean temperature profile is established downstream of the step but the budget for the mean-square temperature fluctuations shows that, in the inner region of the thermal layer, the production and dissipation of temperature fluctuations are not quite equal at the furthest downstream measurement station. The measurements for both the mean and the fluctuating temperature field indicate that the relaxation distance for the thermal layer is quite large, of the order of 1000θ0, where θ0 is the momentum thickness of the boundary layer at the step. Statistics of the thermal-layer interface and conditionally sampled measurements with respect to this interface are presented. Measurements of the temperature intermittency factor indicate that the interface is normally distributed with respect to its mean position. Near the step, the passive heat contaminant acts as an effective marker of the organized turbulence structure that has been observed in the wall region of a boundary layer. Accordingly, conditional averages of Reynolds stresses and heat fluxes measured in the heated part of the flow are considerably larger than the conventional averages when the temperature intermittency factor is small.

Capturability analysis of a geometric guidance law in relative velocity space

In this paper, a relative velocity approach is used to analyze the capturability of a geometric guidance law. Point mass models are assumed for both the missile and the target. The speeds of the missile and target are assumed to remain constant throughout the engagement. Lateral acceleration, obtained from the guidance law, is applied to change the path of the missile. The kinematic equations for engagements in the horizontal plane are derived in the relative velocity space. Some analytical results for the capture region are obtained for non-maneuvering and maneuvering targets. For non-maneuvering targets it is enough for the navigation gain to be a constant to intercept the target, while for maneuvering targets a time varying navigation gain is needed for interception. These results are then verified through numerical simulations.

Sound Velocity Anomaly at the Mott Transition: Application to Organic Conductors and V2O3

Close to the Mott transition, lattice degrees of freedom react to the softening of electron degrees of freedom. This results in a change of lattice spacing, a diverging compressibility, and a critical anomaly of the sound velocity. These effects are investigated within a simple model, in the framework of dynamical mean-field theory. The results compare favorably to recent experiments on the layered organic-conductor kappa-(BEDT-TTF)(2)Cu[N(CN)(2)]Cl. We predict that effects of a similar magnitude are expected for V2O3, despite the much larger value of the elastic modulus of this material.

Experimental Study on the Thermal Argon Plasma Generation and Jet Length Change Characteristics at Atmospheric Pressure

The generation, jet length and flow-regime change characteristics of argon plasma issuing into ambient air have been experimentally examined. Different torch structures have been used in the tests. Laminar plasma jets can be generated within a rather wide range of working-gas flow rates, and an unsteady transitional flow state exists between the laminar and turbulent flow regimes. The high-temperature region length of the laminar plasma jet can be over an order longer than that of the turbulent plasma jet and increases with increasing argon flow rate or arc current, while the jet length of the turbulent plasma is less influenced by the generating parameters. The flow field of the plasma jet has very high radial gradients of plasma parameters, and a Reynolds number alone calculated in the ordinary manner may not adequately serve as a criterion for transition. The laminar plasma jet can have a higher velocity than that of an unsteady or turbulent jet. The long laminar plasma jet has good stiffness to withstand the impact of laterally injected cold gas and particulate matter. It could be used as a rather ideal object for fundamental studies and be applied to novel materials processing due to its attractive stable and adjustable properties.

Electromechanical influence of crack velocity at bifurcation for poled ferroelectric materials

The piezoelastodynamic field equations are solved to determine the crack velocity at bifurcation for poled ferroelectric materials where the applied electrical field and mechanical stress can be varied. The underlying physical mechanism, however, may not correspond to that assumed in the analytical model. Bifurcation has been related to the occurrence of a pair of maximum circumferential stress oriented symmetrically about the moving crack path. The velocity at which this behavior prevails has been referred to as the limiting crack speed. Unlike the classical approach, bifurcation will be identified with finite distances ahead of a moving crack. Nucleation of microcracks can thus be modelled in a single formulation. This can be accomplished by using the energy density function where fracture initiation is identified with dominance of dilatation in relation to distortion. Poled ferroelectric materials are selected for this study because the microstructure effects for this class of materials can be readily reflected by the elastic, piezoelectic and dielectric permittivity constants at the macroscopic scale. Existing test data could also shed light on the trend of the analytical predictions. Numerical results are thus computed for PZT-4 and compared with those for PZT-6B in an effort to show whether the branching behavior would be affected by the difference in the material microstructures. A range of crack bifurcation speed upsilon(b) is found for different r/a and E/sigma ratios. Here, r and a stand for the radial distance and half crack length, respectively, while E and a for the electric field and mechanical stress. For PZT-6B with upsilon(b) in the range 100-1700 m/s, the bifurcation angles varied from +/-6degrees to +/-39degrees. This corresponds to E/sigma of -0.072 to 0.024 V m/N. At the same distance r/a = 0.1, PZT-4 gives upsilon(b) values of 1100-2100 m/s; bifurcation angles of +/-15degrees to +/-49degrees; and E/sigma of -0.056 to 0.059 V m/N. In general, the bifurcation angles +/-theta(0) are found to decrease with decreasing crack velocity as the distance r/a is increased. Relatively speaking, the speed upsilon(b) and angles +/-theta(0) for PZT-4 are much greater than those for PZT-6B. This may be attributed to the high electromechanical coupling effect of PZT-4. Using upsilon(b)(0) as a base reference, an equality relation upsilon(b)(-) < upsilon(b)(0) < upsilon(b)(+) can be established. The superscripts -, 0 and + refer, respectively, to negative, zero and positive electric field. This is reminiscent of the enhancement and retardation of crack growth behavior due to change in poling direction. Bifurcation characteristics are found to be somewhat erratic when r/a approaches the range 10(-2)-10(-1) where the kinetic energy densities would fluctuate and then rise as the distance from the moving crack is increased. This is an artifact introduced by the far away condition of non-vanishing particle velocity. A finite kinetic energy density prevails at infinity unless it is made to vanish in the boundary value problem. Future works are recommended to further clarify the physical mechanism(s) associated with bifurcation by means of analysis and experiment. Damage at the microscopic level needs to be addressed since it has been known to affect the macrocrack speeds and bifurcation characteristics. (C) 2002 Published by Elsevier Science Ltd.

Crystal-growth in floating zone with phase-change and thermosolutal convections

Floating zone crystal growth in microgravity environment is investigated numerically by a finite element method for semiconductor growth processing, which involves thermocapillary convection, phase change convection, thermal diffusion and solutal diffusion. The configurations of phase change interfaces and distributions of velocity, temperature and concentration fields are analyzed for typical conditions of pulling rates and segregation coefficients. The influence of phase change convection on the distribution of concentration is studied in detail. The results show that the thermocapillary convection plays an important role in mixing up the melt with dopant. The deformations of phase change interfaces by thermal convection-diffusion and pulling rods make larger variation of concentration field in comparison with the case of plane interfaces.

Part I. A study of the velocity structure of the earth by the use of core phases. Part II. The 1971 San Fernando earthquake series focal mechanisms and tectonics

The initial objective of Part I was to determine the nature of upper mantle discontinuities, the average velocities through the mantle, and differences between mantle structure under continents and oceans by the use of P'dP', the seismic core phase P'P' (PKPPKP) that reflects at depth d in the mantle. In order to accomplish this, it was found necessary to also investigate core phases themselves and their inferences on core structure. P'dP' at both single stations and at the LASA array in Montana indicates that the following zones are candidates for discontinuities with varying degrees of confidence: 800-950 km, weak; 630-670 km, strongest; 500-600 km, strong but interpretation in doubt; 350-415 km, fair; 280-300 km, strong, varying in depth; 100-200 km, strong, varying in depth, may be the bottom of the low-velocity zone. It is estimated that a single station cannot easily discriminate between asymmetric P'P' and P'dP' for lead times of about 30 sec from the main P'P' phase, but the LASA array reduces this uncertainty range to less than 10 sec. The problems of scatter of P'P' main-phase times, mainly due to asymmetric P'P', incorrect identification of the branch, and lack of the proper velocity structure at the velocity point, are avoided and the analysis shows that one-way travel of P waves through oceanic mantle is delayed by 0.65 to 0.95 sec relative to United States mid-continental mantle.

A new P-wave velocity core model is constructed from observed times, dt/dΔ's, and relative amplitudes of P'; the observed times of SKS, SKKS, and PKiKP; and a new mantle-velocity determination by Jordan and Anderson. The new core model is smooth except for a discontinuity at the inner-core boundary determined to be at a radius of 1215 km. Short-period amplitude data do not require the inner core Q to be significantly lower than that of the outer core. Several lines of evidence show that most, if not all, of the arrivals preceding the DF branch of P' at distances shorter than 143° are due to scattering as proposed by Haddon and not due to spherically symmetric discontinuities just above the inner core as previously believed. Calculation of the travel-time distribution of scattered phases and comparison with published data show that the strongest scattering takes place at or near the core-mantle boundary close to the seismic station.

In Part II, the largest events in the San Fernando earthquake series, initiated by the main shock at 14 00 41.8 GMT on February 9, 1971, were chosen for analysis from the first three months of activity, 87 events in all. The initial rupture location coincides with the lower, northernmost edge of the main north-dipping thrust fault and the aftershock distribution. The best focal mechanism fit to the main shock P-wave first motions constrains the fault plane parameters to: strike, N 67° (± 6°) W; dip, 52° (± 3°) NE; rake, 72° (67°-95°) left lateral. Focal mechanisms of the aftershocks clearly outline a downstep of the western edge of the main thrust fault surface along a northeast-trending flexure. Faulting on this downstep is left-lateral strike-slip and dominates the strain release of the aftershock series, which indicates that the downstep limited the main event rupture on the west. The main thrust fault surface dips at about 35° to the northeast at shallow depths and probably steepens to 50° below a depth of 8 km. This steep dip at depth is a characteristic of other thrust faults in the Transverse Ranges and indicates the presence at depth of laterally-varying vertical forces that are probably due to buckling or overriding that causes some upward redirection of a dominant north-south horizontal compression. Two sets of events exhibit normal dip-slip motion with shallow hypocenters and correlate with areas of ground subsidence deduced from gravity data. Several lines of evidence indicate that a horizontal compressional stress in a north or north-northwest direction was added to the stresses in the aftershock area 12 days after the main shock. After this change, events were contained in bursts along the downstep and sequencing within the bursts provides evidence for an earthquake-triggering phenomenon that propagates with speeds of 5 to 15 km/day. Seismicity before the San Fernando series and the mapped structure of the area suggest that the downstep of the main fault surface is not a localized discontinuity but is part of a zone of weakness extending from Point Dume, near Malibu, to Palmdale on the San Andreas fault. This zone is interpreted as a decoupling boundary between crustal blocks that permits them to deform separately in the prevalent crustal-shortening mode of the Transverse Ranges region.

Laser-induced changes on the complex refractive indices of phase-change thin film

The refractive indices of crystalline phase-change films are usually obtained by thermal-induced crystallization. However, this is not accurate, because the crystallization of phase-change film in rewritable optical disks is laser induced. In this study, we use the initializer to crystallize the phase-change films. The dependence of the refractive index n and the extinction coefficient k of the phase-change films on the initialization conditions are studied. Remarkable changes of the refractive indices (especially k) are found when the initialization laser power density is 6.63 mW/mum(2) and the initialization velocity is 4.0 m/s. At the same time, the structure changes of the phase-change films are also studied. This dependence is explained by the structure change of the films. These results are significant in improving the accuracy of optical design and the thermal simulation of phase-change optical disks, as well as in the study of phase-change optical disks at shorter wavelengths. (C) 2003 Society of Photo-Optical Instrumentation Engineers.

Methane hydrate growth velocity in seepage system of Qiongdongnan Basin, South China Sea

The occurrences of diapirs, gas-filled zones and gas plumes in seawater in Qiongdongnan Basin of South China Sea indicate that there may exist seepage system gas-hydrate reservoirs. Assuming there has a methane venting zone of 1500 m in diameter, and the methane flux is 1000 kmol/a, and the temperature of methane hydrate-bearing sediments ranges from 3 degrees C to 20 degrees C, then according to the hydrate film growth theory, by numerical simulation, this paper computes the temperatures and velocities in 0 mbsf, 100 mbsf, 200 mbsf, 425 mbsf over discrete length, and gives the change charts. The results show that the cementation velocity in sediments matrix of methane hydrate is about 0.2 nm/s, and the seepage system will evolve into diffusion system over probably 35000 years. Meanwhile, the methane hydrate growth velocity in leakage system is 20 similar to 40 times faster than in diffusion system.

The change of characteristics of tidal wave in the Changjiang River mouth area since the post-glacial transgression maximum

During the period of the post-glacial transgression maximum (PGTM), there was a huge trumpet estuary in the modern Changjiang River Delta area. The location and the shape of the Paleo-Changjiang River Estuary (PCRE) were much different from those of the present Changjiang River Estuary. The study on the change of characteristics of tidal wave in the Changjiang River mouth area since the PGTM can help to understand better the dynamic development of the Changjiang River Delta. The course curves of tidal level and tidal current velocity during a single tidal cycle for 35 points are calculated, and characteristics of tidal waves in the PCRE and its adjacent area are compared with those of tidal waves in the modern Changjiang River mouth area. The results show that the tidal waves within the PCRE and in its adjacent area during the period of the PGTM belonged to standing wave or a mixture of standing wave and progressive wave. Since then, the tidal wave in the Changjiang River mouth become gradually to be progressive wave with the PCRE being filled and the Changjiang River mouth shifting southeastwards.

Climate change impacts on sea–air fluxes of CO2 in three Arctic seas: a sensitivity study using Earth observation

We applied coincident Earth observation data collected during 2008 and 2009 from multiple sensors (RA2, AATSR and MERIS, mounted on the European Space Agency satellite Envisat) to characterise environmental conditions and integrated sea-air fluxes of CO2 in three Arctic seas (Greenland, Barents, Kara). We assessed net CO2 sink sensitivity due to changes in temperature, salinity and sea ice duration arising from future climate scenarios. During the study period the Greenland and Barents seas were net sinks for atmospheric CO2, with integrated sea-air fluxes of -36 +/- 14 and -11 +/- 5 Tg C yr(-1), respectively, and the Kara Sea was a weak net CO2 source with an integrated sea-air flux of +2.2 +/- 1.4 TgC yr(-1). The combined integrated CO2 sea-air flux from all three was -45 +/- 18 TgC yr(-1). In a sensitivity analysis we varied temperature, salinity and sea ice duration. Variations in temperature and salinity led to modification of the transfer velocity, solubility and partial pressure of CO2 taking into account the resultant variations in alkalinity and dissolved organic carbon (DOC). Our results showed that warming had a strong positive effect on the annual integrated sea-air flux of CO2 (i.e. reducing the sink), freshening had a strong negative effect and reduced sea ice duration had a small but measurable positive effect. In the climate change scenario examined, the effects of warming in just over a decade of climate change up to 2020 outweighed the combined effects of freshening and reduced sea ice duration. Collectively these effects gave an integrated sea-air flux change of +4.0 TgC in the Greenland Sea, +6.0 Tg C in the Barents Sea and +1.7 Tg C in the Kara Sea, reducing the Greenland and Barents sinks by 11% and 53 %, respectively, and increasing the weak Kara Sea source by 81 %. Overall, the regional integrated flux changed by +11.7 Tg C, which is a 26% reduction in the regional sink. In terms of CO2 sink strength, we conclude that the Barents Sea is the most susceptible of the three regions to the climate changes examined. Our results imply that the region will cease to be a net CO2 sink in the 2050s.

Stellar jitter from variable gravitational redshift: implications for radial velocity confirmation of habitable exoplanets

A variation of gravitational redshift, arising from stellar radius fluctuations, will introduce astrophysical noise into radial velocity measurements by shifting the centroid of the observed spectral lines. Shifting the centroid does not necessarily introduce line asymmetries. This is fundamentally different from other types of stellar jitter so far identified, which do result from line asymmetries. Furthermore, only a very small change in stellar radius, ˜0.01 per cent, is necessary to generate a gravitational redshift variation large enough to mask or mimic an Earth-twin. We explore possible mechanisms for stellar radius fluctuations in low-mass stars. Convective inhibition due to varying magnetic field strengths and the Wilson depression of starspots are both found to induce substantial gravitational redshift variations. Finally, we investigate a possible method for monitoring/correcting this newly identified potential source of jitter and comment on its impact for future exoplanet searches.

The effects of an acute bout of whole-body vibration on pulse wave velocity in individuals with SCI

Cardiovascular disease is a leading cause of mortality in the spinal cord injured (SCI) population. Reduced arterial compliance is a cardiovascular risk factor and whole body vibration (WBV) has be en shown to improve arterial compliance in able-bodied individuals. The study investigated the effect of an acute session ofWBV on arterial compliance as measured by pulse wave velocity (PWV). On separate days, arm, leg and aortic PWV were measured pre- and post- a 45 minute session of passive stance (PS) and WBV. The WBV was intermittent with a set frequency of 45Hz and amplitude of O.6mm. There was no condition by time effect when comparing PWV after WBV and PS. Following WBV, aortic (928.6±127.7 vs. 901.1±96.6cm/sec), leg (1035.2±113.8 vs.l099.8±114.2cm/sec) and arm PWV (1118.9±119.8 vs. 1181.1±124.4cm/s) did not change. As such, WBV did not reduce arterial compliance, however future research with protocol modifications is recommended.

Numerical Simulation of Electroosmotic Flow with Step Change in Zeta Potential

Electroosmotic flow is a convenient mechanism for transporting polar fluid in a microfluidic device. The flow is generated through the application of an external electric field that acts on the free charges that exists in a thin Debye layer at the channel walls. The charge on the wall is due to the chemistry of the solid-fluid interface, and it can vary along the channel, e.g. due to modification of the wall. This investigation focuses on the simulation of the electroosmotic flow (EOF) profile in a cylindrical microchannel with step change in zeta potential. The modified Navier-Stoke equation governing the velocity field and a non-linear two-dimensional Poisson-Boltzmann equation governing the electrical double-layer (EDL) field distribution are solved numerically using finite control-volume method. Continuities of flow rate and electric current are enforced resulting in a non-uniform electrical field and pressure gradient distribution along the channel. The resulting parabolic velocity distribution at the junction of the step change in zeta potential, which is more typical of a pressure-driven velocity flow profile, is obtained.