A New Ignition Criterion Based Upon a Critical Thermal Energy Density

It has long been known that various ignition criteria of energetic materials have been limited in applicability to small regions. In order to explore the physical nature of ignition, we calculated how much thermal energy per unit mass of energetic materials was absorbed under different external stimuli. Hence, data of several typical sensitivity tests were analyzed by order of magnitude estimation. Then a new concept on critical thermal energy density was formulated. Meanwhile, the chemical nature of ignition was probed into by chemical kinetics.

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.

Marangoni-Benard instability with the exchange of evaporation at liquid-vapour interface

A new two-sided model rather than the one-sided model in previous works is put forward. The linear instability analysis is performed on the Marangoni-Benard convection in the two-layer system with an evaporation interface. We define a new evaporation Biot number which is different from that in the one-sided model, and obtain the curves of critical Marangoni number versus wavenumber. The influence of evaporation velocity and Biot number on the system is discussed and a new phenomenon uninterpreted before is now explained from our numerical results.

Optimized Group Velocity Control Scheme and DNS of Decaying Compressible Turbulence of Relative High Turbulent Mach Number

To overcome the difficulty in the DNS of compressible turbulence at high turbulent Mach number, a new difference scheme called GVC8 is developed. We have succeeded in the direct numerical simulation of decaying compressible turbulence up to turbulent Mach number 0.95. The statistical quantities thus obtained at lower turbulent Mach number agree well with those from previous authors with the same initial conditions, but they are limited to simulate at lower turbulent Mach numbers due to the so-called start-up problem. The energy spectrum and coherent structure of compressible turbulent flow are analysed. The scaling law of compressible turbulence is studied. The computed results indicate that the extended self-similarity holds in decaying compressible turbulence despite the occurrence of shocklets, and compressibility has little effects on relative scaling exponents when turbulent Mach number is not very high.

Calculations of longitudinal and transverse velocity structure functions using a vortex model of isotropic turbulence

The longitudinal structure function (LSF) and the transverse structure function (TSF) in isotropic turbulence are calculated using a vortex model. The vortex model is composed of the Rankine and Burgers vortices which have the exponential distributions in the vortex Reynolds number and vortex radii. This model exhibits a power law in the inertial range and satisfies the minimal condition of isotropy that the second-order exponent of the LSF in the inertial range is equal to that of the TSF. Also observed are differences between longitudinal and transverse structure functions caused by intermittency. These differences are related to their scaling differences which have been previously observed in experiments and numerical simulations.

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.

Modeling Nonlinear Peeling of Ductile Thin Films-Critical Assessment of Analytical Bending Models Using FE Simulations

Three analytical double-parameter criteria based on a bending model and a two-dimensional finite element analysis model are presented for the modeling of ductile thin film undergoing a nonlinear peeling process. The bending model is based on different governing parameters: (1) the interfacial fracture toughness and the separation strength, (2) the interfacial fracture toughness and the crack tip slope angle, and (3) the interfacial fracture toughness and the critical Mises effective strain of the delaminated thin film at the crack tip. Thin film nonlinear peeling under steady-state condition is solved with the different governing parameters. In addition, the peeling test problem is simulated by using the elastic-plastic finite element analysis model. A critical assessment of the three analytical bending models is made by comparison of the bending model solutions with the finite element analysis model solutions. Furthermore, through analyses and comparisons for solutions based on both the bending model and the finite element analysis model, some connections between the bending model and the finite element analysis model are developed. Moreover, in the present research, the effect of different selections for cohesive zone shape on the ductile film peeling solutions is discussed.

Influences of slope gradient on soil erosion

The main factors influencing soil erosion include the net rain excess, the water depth, the velocity, the shear stress of overland flows, and the erosion-resisting capacity of soil. The laws of these factors varying with the slope gradient were investigated by using the kinematic wave theory. Furthermore, the critical slope gradient of erosion was driven. The analysis shows that the critical slope gradient of soil erosion is dependent on grain size, soil bulk density, surface roughness, runoff length, net rain excess, and the friction coefficient of soil, etc. The critical slope gradient has been estimated theoretically with its range between 41.5 degrees similar to 50 degrees.

Experimental Investigation of the Velocity Effect on Adhesion Forces with an Atomic Force Microscope Atomic Force Microscope

Capillary forces are significantly dominant in adhesive forces measured with an atomic force microscope (AFM) in ambient air, which are always thought to be dependent on water film thickness, relative humidity, and the free energy of water film. We study the nature of the pull-off force on a variety of surfaces as a function of tip velocity. It is found that the capillary forces are of relatively strong dependence on tip velocity. The present experiment is expected to provide a better understanding of the work mechanism of AFM in ambient air.

Critical Sensitivity and Trans-Scale Fluctuations in Catastrophic Rupture

Rupture in the heterogeneous crust appears to be a catastrophe transition. Catastrophic rupture sensitively depends on the details of heterogeneity and stress transfer on multiple scales. These are difficult to identify and deal with. As a result, the threshold of earthquake-like rupture presents uncertainty. This may be the root of the difficulty of earthquake prediction. Based on a coupled pattern mapping model, we represent critical sensitivity and trans-scale fluctuations associated with catastrophic rupture. Critical sensitivity means that a system may become significantly sensitive near catastrophe transition. Trans-scale fluctuations mean that the level of stress fluctuations increases strongly and the spatial scale of stress and damage fluctuations evolves from the mesoscopic heterogeneity scale to the macroscopic scale as the catastrophe regime is approached. The underlying mechanism behind critical sensitivity and trans-scale fluctuations is the coupling effect between heterogeneity and dynamical nonlinearity. Such features may provide clues for prediction of catastrophic rupture, like material failure and great earthquakes. Critical sensitivity may be the physical mechanism underlying a promising earthquake forecasting method, the load-unload response ratio (LURR).

Thermocapillary Convection in Floating Zones

This paper provides an overview of ongoing studies in the area of thermocapillary convection driven by a surface tension gradient parallel to the free surface in a floating zone. Here, research interests are focused around the onset of oscillatory thermocapillary convection, also known as the transition from quasisteady convection to oscillatory convection. The onset of oscillation depends on a set of critical parameters, and the margin relationship can be represented by a complex function of the critical parameters. The experimental results indicate that the velocity deviation of an oscillatory flow has the same order of magnitude as that of an average flow, and the deviations of other quantities, such as temperature and free surface radii fluctuations, are much smaller when compared with their normal counterparts. Therefore, the onset of oscillation should be a result of the dynamic process in a fluid, and the problem is a strongly nonlinear one. In the past few decades, several theoretical models have been introduced to tackle the problem using analytical methods, linear instability analysis methods, energy instability methods, and unsteady 3D numerical methods. The last of the above mentioned methods is known to be the most suitable for a thorough analysis of strong nonlinear processes, which generally leads to a better comparison with the experimental results. The transition from oscillatory thermocapillary convection to turbulence falls under the studies of chaotic behavior in a new system, which opens a fascinating new frontier in nonlinear science, a hot research area drawing many recent works. This paper reviews theoretical models and analysis, and also experimental research, on thermocapillary connection in floating zones. It cites 93 references.

The Hydrothermal Wave Of Large-Prandtl-Number Fluid In A Shallow Cavity

The hydrothermal wave was investigated numerically for large-Prandtl-number fluid (Pr = 105.6) in a shallow cavity with different heated sidewalls. The traveling wave appears and propagates in the direction opposite to the surface flow (upstream) in the case of zero gravity when the applied temperature difference grows and over the critical value. The phase relationships of the disturbed velocity, temperature and pressure demonstrate that the traveling wave is driven by the disturbed temperature, which is named hydrothermal wave. The hydrothermal wave is so weak that the oscillatory flow field and temperature distribution can hardly be observed in the liquid layer. The exciting mechanism of hydrothermal wave is analyzed and discussed in the present paper.

Influence Of Free Surface Curvature Of A Liquid Layer On The Critical Marangoni Convection

The Pearson instability was suggested to discuss the onset of Marangoni convection in a liquid layer of large Prandtl number under an applied temperature difference perpendicular to the free surface in the microgravity environment. In this case, the temperature distribution on the curved free surface is nonuniform, and the thermocapillary convection is induced and coupled with the Marangoni convection. In the present paper the effect of volume ratio of the liquid layer on the critical Marangoni convection and the corresponding spatial variation of the convection structure in zero-gravity condition were numerically investigated by two-dimensional model. (C) 2008 Elsevier Ltd. All rights reserved.

The Head-On Colliding Process Of Binary Liquid Droplets At Low Velocity: High-Speed Photography Experiments And Modeling

The experimental and theoretical studies are reported in this paper for the head-on collisions of a liquid droplet with another of the same fluid resting on a solid substrate. The droplet on the hydrophobic polydimethylsiloxane (PDMS) substrate remains in a shape of an approximately spherical segment and is isometric to an incoming droplet. The colliding process of the binary droplets was recorded with high-speed photography. Head-on collisions saw four different types of response in our experiments: complete rebound, coalescence, partial rebound With conglutination, and coalescence accompanied by conglutination. For a complete rebound, both droplets exhibited remarkable elasticity and the contact time of the two colliding droplets was found to be in the range of 10-20 ms. With both droplets approximately considered as elastic bodies, Hertz contact theory was introduced to estimate the contact time for the complete rebound case. The estimated result Was found to be on the same order of magnitude as the experimental data, which indicates that the present model is reasonable. (C) 2008 Elsevier Inc. All rights reserved.

Experimental Investigation Of Particle Velocity Distributions In Windblown Sand Movement

With the PDPA (Phase Doppler Particle Analyzer) measurement technology, the probability distributions of particle impact and lift-off velocities on bed surface and the particle velocity distributions at different heights are detected in a wind tunnel. The results show that the probability distribution of impact and lift-off velocities of sand grains can be expressed by a log-normal function, and that of impact and lift-off angles complies with an exponential function. The mean impact angle is between 28 degrees and 39 degrees, and the mean lift-off angle ranges from 30 degrees to 44 degrees. The mean lift-off velocity is 0.81-0.9 times the mean impact velocity. The proportion of backward-impacting particles is 0.05-0.11, and that of backward-entrained particles ranges from 0.04 to 0.13. The probability distribution of particle horizontal velocity at 4 mm height is positive skew, the horizontal velocity of particles at 20 mm height varies widely, and the variation of the particle horizontal velocity at 80 mm height is less than that at 20 mm height. The probability distribution of particle vertical velocity at different heights can be described as a normal function.