Numerical investigation of particle velocity distributions in aeolian sand transport

Particle velocity distribution in a blowing sand cloud is a reflection of saltation movement of many particles. Numerical analysis is performed for particle velocity distribution with a discrete particle model. The probability distributions of resultant particle velocity in the impact-entrainment process, particle horizontal and vertical velocities at different heights and the vertical velocity of ascending particles are analyzed. The probability distributions of resultant impact and lift-off velocities of saltating particles can be expressed by a log-normal function, and that of impact angle comply with an exponential function. The probability distribution of particle horizontal and vertical velocities at different heights shows a typical single-peak pattern. In the lower part of saltation layer, the particle horizontal velocity distribution is positively skewed. Further analysis shows that the probability density function of the vertical velocity of ascending particles is similar to the right-hand part of a normal distribution function, and a general equation is acquired for the probability density function of non-dimensional vertical velocity of ascending particles which is independent of diameter of saltating particles, wind strength and height. These distributions in the present numerical analysis are consistent with reported experimental results. The present investigation is important for understanding the saltation state in wind-blown sand movement. (C) 2009 Elsevier B.V. All rights reserved.

Vertical velocity variance and skewness in clear and cloud-topped boundary layers as revealed by Doppler lidar

In this paper, observations by a ground-based vertically pointing Doppler lidar and sonic anemometer are used to investigate the diurnal evolution of boundary-layer turbulence in cloudless, cumulus and stratocumulus conditions. When turbulence is driven primarily by surface heating, such as in cloudless and cumulus-topped boundary layers, both the vertical velocity variance and skewness follow similar profiles, on average, to previous observational studies of turbulence in convective conditions, with a peak skewness of around 0.8 in the upper third of the mixed layer. When the turbulence is driven primarily by cloud-top radiative cooling, such as in the presence of nocturnal stratocumulus, it is found that the skewness is inverted in both sign and height: its minimum value of around −0.9 occurs in the lower third of the mixed layer. The profile of variance is consistent with a cloud-top cooling rate of around 30Wm−2. This is also consistent with the evolution of the thermodynamic profile and the rate of growth of the mixed layer into the stable nocturnal boundary layer from above. In conditions where surface heating occurs simultaneously with cloud-top cooling, the skewness is found to be useful for diagnosing the source of the turbulence, suggesting that long-term Doppler lidar observations would be valuable for evaluating boundary-layer parametrization schemes. Copyright c 2009 Royal Meteorological Society

Secular changes in BMI and obesity risk in Japanese children: Considerations from a morphologic perspective

The increasing prevalence of childhood obesity is a global health issue. Past studies in Japan have reported an increase in both body mass index (BMI) and risk of obesity among children and adolescents. However, changes in body size and proportion in this population over time have also influenced BMI. To date, no study of secular changes in childhood obesity has considered the impact of changes in morphological factors. The current study explored the secular changes in BMI and childhood obesity risk among Japanese children from 1950 to 2000 with consideration of changes in body size and the proportions using The Statistical Report of the School Health Survey (SHS). The age of peak velocity (PV) occurred approximately two years earlier in both genders across this period. While the increments in height, sitting height and sub-ischial leg length relative to height levelled off by 1980, weight gain continued in boys. Between 1980 and 2000, the rate of the upper body weight gain in boys and girls were 0.7-1.3 kg/decade and 0.2-1.0 kg/decade, respectively. After considering body proportions, increments in body weight were small. It could be suggested that the increments in weight and BMI across the 50-year period may be due to a combination of changes including the tempo of growth and body size due to lifestyle factors.

Numerical simulation of idealised three-dimensional downburst wind fields

Convective downburst wind storms generate the peak annual gust wind speed for many parts of the non-cyclonic world at return periods of importance for ultimate limit state design. Despite this there is little clear understanding of how to appropriately design for these wind events given their significant dissimilarities to boundary layer winds upon which most design is based. To enhance the understanding of wind fields associated with these storms a three-dimensional numerical model was developed to simulate a multitude of idealised downburst scenarios and to investigate their near-ground wind characteristics. Stationary and translating downdraft wind events in still and sheared environments were simulated with baseline results showing good agreement with previous numerical work and full-scale observational data. Significant differences are shown in the normalised peak wind speed velocity profiles depending on the environmental wind conditions in the vicinity of the simulated event. When integrated over the height of mid- to high rise structures, all simulated profiles are shown to produce wind loads smaller than an equivalent 10 m height matched open terrain boundary layer profile. This suggests that for these structures the current design approach is conservative from an ultimate loading standpoint. Investigating the influence of topography on the structure of the simulated near-ground downburst wind fields, it is shown that these features amplify wind speeds in a manner similar to that expected for boundary layer winds, but the extent of amplification is reduced. The level of reduction is shown to be dependent on the depth of the simulated downburst outflow.

Large Low-Frequency Fluctuations in the Velocity of a Driven Vortex Lattice in a Single Crystal of 2H-NbSe2 Superconductor

The driven state of a well-ordered flux line lattice in a single crystal of 2H-NbSe2 in the time domain has revealed the presence of substantial fluctuations in velocity, with large and distinct time periods (similar to seconds). A superposition of a periodic drive in the driven vortex lattice causes distinct changes in these fluctuations. We propose that prior to the onset of the peak effect there exists a heretofore unexplored regime of coherent dynamics, with unexpected behavior in velocity fluctuations.

Prediction of flame liftoff height of diffusion/partially premixed jet flames and modeling of mild combustion burners

In this article, a new flame extinction model based on the k/epsilon turbulence time scale concept is proposed to predict the flame liftoff heights over a wide range of coflow temperature and O-2 mass fraction of the coflow. The flame is assumed to be quenched, when the fluid time scale is less than the chemical time scale ( Da < 1). The chemical time scale is derived as a function of temperature, oxidizer mass fraction, fuel dilution, velocity of the jet and fuel type. The present extinction model has been tested for a variety of conditions: ( a) ambient coflow conditions ( 1 atm and 300 K) for propane, methane and hydrogen jet flames, ( b) highly preheated coflow, and ( c) high temperature and low oxidizer concentration coflow. Predicted flame liftoff heights of jet diffusion and partially premixed flames are in excellent agreement with the experimental data for all the simulated conditions and fuels. It is observed that flame stabilization occurs at a point near the stoichiometric mixture fraction surface, where the local flow velocity is equal to the local flame propagation speed. The present method is used to determine the chemical time scale for the conditions existing in the mild/ flameless combustion burners investigated by the authors earlier. This model has successfully predicted the initial premixing of the fuel with combustion products before the combustion reaction initiates. It has been inferred from these numerical simulations that fuel injection is followed by intense premixing with hot combustion products in the primary zone and combustion reaction follows further downstream. Reaction rate contours suggest that reaction takes place over a large volume and the magnitude of the combustion reaction is lower compared to the conventional combustion mode. The appearance of attached flames in the mild combustion burners at low thermal inputs is also predicted, which is due to lower average jet velocity and larger residence times in the near injection zone.

Vibrational dynamics and boson peak in a supercooled polydisperse liquid

Vibrational density of states (VDOS) in a supercooled polydisperse liquid is computed by diagonalizing the Hessian matrix evaluated at the potential energy minima for systems with different values of polydispersity. An increase in polydispersity leads to an increase in the relative population of localized high-frequency modes. At low frequencies, the density of states shows an excess compared to the Debye squared-frequency law, which has been identified with the boson peak. The height of the boson peak increases with polydispersity and shows a rather narrow sensitivity to changes in temperature. While the modes comprising the boson peak appear to be largely delocalized, there is a sharp drop in the participation ratio of the modes that exist just below the boson peak indicative of the quasilocalized nature of the low-frequency vibrations. Study of the difference spectrum at two different polydispersity reveals that the increase in the height of boson peak is due to a population shift from modes with frequencies above the maximum in the VDOS to that below the maximum, indicating an increase in the fraction of the unstable modes in the system. The latter is further supported by the facilitation of the observed dynamics by polydispersity. Since the strength of the liquid increases with polydispersity, the present result provides an evidence that the intensity of boson peak correlates positively with the strength of the liquid, as observed earlier in many experimental systems.

Support vector machine for evaluating seismic-liquefaction potential using shear wave velocity

The use of the shear wave velocity data as a field index for evaluating the liquefaction potential of sands is receiving increased attention because both shear wave velocity and liquefaction resistance are similarly influenced by many of the same factors such as void ratio, state of stress, stress history and geologic age. In this paper, the potential of support vector machine (SVM) based classification approach has been used to assess the liquefaction potential from actual shear wave velocity data. In this approach, an approximate implementation of a structural risk minimization (SRM) induction principle is done, which aims at minimizing a bound on the generalization error of a model rather than minimizing only the mean square error over the data set. Here SVM has been used as a classification tool to predict liquefaction potential of a soil based on shear wave velocity. The dataset consists the information of soil characteristics such as effective vertical stress (sigma'(v0)), soil type, shear wave velocity (V-s) and earthquake parameters such as peak horizontal acceleration (a(max)) and earthquake magnitude (M). Out of the available 186 datasets, 130 are considered for training and remaining 56 are used for testing the model. The study indicated that SVM can successfully model the complex relationship between seismic parameters, soil parameters and the liquefaction potential. In the model based on soil characteristics, the input parameters used are sigma'(v0), soil type. V-s, a(max) and M. In the other model based on shear wave velocity alone uses V-s, a(max) and M as input parameters. In this paper, it has been demonstrated that Vs alone can be used to predict the liquefaction potential of a soil using a support vector machine model. (C) 2010 Elsevier B.V. All rights reserved.

Dark Matter Halo Properties as Deduced from the Observed H I Scale Height Data, in The Low-Frequency Radio Universe

We use the HΙ scale height data along with the HΙ rotation curve as constraints to probe the shape and density profile of the dark matter halos of M31 (Andromeda) and the superthin, low surface brightness (LSB) galaxy UGC 07321. We model the galaxy as a two component system of gravitationally-coupled stars and gas subjected to the force field of a dark matter halo. For M31, we get a flattened halo which is required to match the outer galactic HΙ scale height data, with our best-fit axis ratio (0.4) lying at the most oblate end of the distributions obtained from cosmological simulations. For UGC 07321, our best-fit halo core radius is only slightly larger than the stellar disc scale length, indicating that the halo is important even at small radii in this LSB galaxy. The high value of the gas velocity dispersion required to match the scale height data can explain the low star-formation rate of this galaxy.

Universality of scaling laws in correlation between velocity and shear stress in turbulent boundary layers

In this study, we analyse simultaneous measurements (at 50 Hz) of velocity at several heights and shear stress at the surface made during the Utah field campaign for the presence of ranges of scales, where distinct scale-to-scale interactions between velocity and shear stress can be identified. We find that our results are similar to those obtained in a previous study [Venugopal et al., 2003] (contrary to the claim in V2003, that the scaling relations might be dependent on Reynolds number) where wind tunnel measurements of velocity and shear stress were analysed. We use a wavelet-based scale-to-scale cross-correlation to detect three ranges of scales of interaction between velocity and shear stress, namely, (a) inertial subrange, where the correlation is negligible; (b) energy production range, where the correlation follows a logarithmic law; and (c) for scales larger than the boundary layer height, the correlation reaches a plateau.

On the height of fall of dropping mass in SASW measurements for asphaltic road pavements

A number of spectral analysis of surface waves (SASW) tests were performed on asphaltic road pavements by dropping a metallic 6.5 kg sphere, from a height (H) ranging from 1 to 3 m. Various combinations of source to first receiver distance (S) and receiver spacing (X) were employed. By increasing the height of the fall of the dropping mass, the maximum wavelength (lambda(max)), up to which the shear wave velocity profile can be predicted with the usage of the SASW measurements, was found to increase continuously. The height of fall of the dropping mass also seems to affect the admissible range of the wavelength for given combinations of X and S. Irrespective of different chosen combinations of S, X and H, a unique combined dispersion curve was generated in all the cases for a given pavement site as long as the threshold minimum value of the coherence function is greater than 0.90.

Wallner lines, crack velocity and mechanisms of crack nucleation and growth in a brittle bulk metallic glass

Mode I fracture experiments were conducted on brittle bulk metallic glass (BMG) samples and the fracture surface features were analyzed in detail to understand the underlying physical processes. Wollner lines, which result from the interaction between the propagating crack front and shear waves emanating from a secondary source, were observed on the fracture surface and geometric analysis of them indicates that the maximum crack velocity is similar to 800 m s(-1), which corresponds to similar to 0.32 times the shear wave speed. Fractography reveals that the sharp crack nucleation at the notch tip occurs at the mid-section of the specimens with the observation of flat and half-penny-shaped cracks. On this basis, we conclude that the crack initiation in brittle BMGs is stress-controlled and occurs through hydrostatic stress-assisted cavity nucleation ahead of the notch tip. High magnification scanning electron and atomic force microscopies of the dynamic crack growth regions reveal highly organized, nanoscale periodic patterns with a spacing of similar to 79 nm. Juxtaposition of the crack velocity with this spacing suggests that the crack takes similar to 10(-10) s for peak-to-peak propagation. This, and the estimated adiabatic temperature rise ahead of the propagating crack tip that suggests local softening, is utilized to critically discuss possible causes for the nanocorrugation formation. Taylor's fluid meniscus instability is unequivocally ruled out. Then, two other possible mechanisms, viz. (a) crack tip blunting and resharpening through nanovoid nucleation and growth ahead of the crack tip and eventual coalescence, and (b) dynamic oscillation of the crack in a thin slab of softened zone ahead of the crack-tip, are critically discussed. (C) 2014 Acta Materialia Inc. Published by Elsevier Ltd. All rights reserved.

Air microjet system for non-contact force application and the actuation of micro-structures

We demonstrate a non-contact technique to apply calibrated and localized forces in the micro-Newton to milli-Newton range using an air microjet. An electromagnetically actuated diaphragm controlled by a signal generator is used to generate the air microjet. With a nozzle diameter of 150 mu m, the microjet diameter was maintained to a maximum of 1 mm at a distance of 5 mm from the nozzle. The force generated by the microjet was measured using a commercial force sensor to determine the velocity profile of the jet. Axial flow velocities of up to 25 m s(-1) were obtained at distances as long as 6 mm. The microjet exerted a force up to 1 mu N on a poly dimethyl siloxane (PDMS) micropillar (50 mu m in diameter, 157 mu m in height) and 415 mu N on a PDMS membrane (3 mm in diameter, 28 mu m thick). We also demonstrate that from a distance of 6 mm our microjet can exert a peak pressure of 187 Pa with a total force of about 84 mu N on a flat surface with 8 V operating voltage. Out of the cleanroom fabrication and robust design make this system cost effective and durable.

DPIV and Interferometry Technique for Velocity Field and Concentration Fields Measurement in Crystal Growth

The property of crystal depends seriously on the solution concentration distribution near the growth surface of a crystal. However, the concentration distributions are affected by the diffusion and convection of the solution. In the present experiment, the two methods of optical measurement are used to obtained velocity field and concentration field of NaClO₃ solution. The convection patterns in sodium chlorate (NaClO₃) crystal growth are measured by Digital Particle image Velocimetry (DPIV) technology. The 2-dimentional velocity distributions in the solution of NaClO3 are obtained from experiments. And concentration field are obtained by a Mach-Zehnder interferometer with a phase shift servo system. Interference patterns were recorded directly by a computer via a CCD camera. The evolution of velocity field and concentration field from dissolution to crystallization are visualized clearly. The structures of velocity fields were compared with that of concentration field.