Information Preservation Modeling of Rayleigh-Benard Transition from Thermal Conduction to Convection

Onset and evolution of the Rayleigh-Benard (R-B) convection are investigated using the Information Preservation (IP) method. The information velocity and temperature are updated using the Octant Flux Splitting (OFS) model developed by Masters & Ye based on the Maxwell transport equation suggested by Sun & Boyd. Statistical noise inherent in particle approaches such as the direct simulation Monte Carlo (DSMC) method is effectively reduced by the IP method, and therefore the evolutions from an initial quiescent fluid to a final steady state are shown clearly. An interesting phenomenon is observed: when the Rayleigh number (Ra) exceeds its critical value, there exists an obvious incubation stage. During the incubation stage, the vortex structure clearly appears and evolves, whereas the Nusselt number (Nu) of the lower plate is close to unity. After the incubation stage, the vortex velocity and Nu rapidly increase, and the flow field quickly reaches a steady, convective state. A relation of Nu to Ra given by IP agrees with those given by DSMC, the classical theory and experimental data.

Curie Transition Of Nc Nickel By Mechanical Spectroscopy And Magnetization Study

Mechanical spectroscopy measurement is performed to study the internal friction of nanocrystalline ( NC) nickel with an average grain size of 23 nm from room temperature to 610 K. An internal friction peak is observed at about 550 K, which corresponds to the Curie transition process of the NC nickel according to the result of magnetization test. Moreover, the fact that the Curie temperature of NC nickel is lower than that of coarse-grained nickel is explained by an analytical model based on the weakening of cohesive energy.

Numerical study on heavy rigid particle motion in a plane wake flow by spectral element method

Experimental particle dispersion patterns in a plane wake flow at a high Reynolds number have been predicted numerically by discrete vortex method (Phys. Fluids A 1992; 4:2244-2251; Int. J. Multiphase Flow 2000; 26:1583-1607). To address the particle motion at a moderate Reynolds number, spectral element method is employed to provide an instantaneous wake flow field for particle dynamics equations, which are solved to make a detail classification of the patterns in relation to the Stokes and Froude numbers. It is found that particle motion features only depend on the Stokes number at a high Froude number and depend on both numbers at a low Froude number. A ratio of the Stokes number to squared Froude number is introduced and threshold values of this parameter are evaluated that delineate the different regions of particle behavior. The parameter describes approximately the gravitational settling velocity divided by the characteristic velocity of wake flow. In order to present effects of particle density but preserve rigid sphere, hollow sphere particle dynamics in the plane wake flow is investigated. The evolution of hollow particle motion patterns for the increase of equivalent particle density corresponds to that of solid particle motion patterns for the decrease of particle size. Although the thresholds change a little, the parameter can still make a good qualitative classification of particle motion patterns as the inner diameter changes.

Nanoscale periodic corrugation to dimple transition due to "beat" in a bulk metallic glass

We report an intriguing observation that the interaction of brittle nanoscale periodic corrugations (NPCs) can lead to the formation of ductile dimples on the dynamic fracture surface of a tough Vit 1 bulk metallic glass (BMG) under high-velocity plate impact. A “beat” phenomenon due to superposition of simple harmonic vibrations, approximately characterizing NPCs, is proposed to explain this unusual brittle-to-ductile transition. The present results agree well with our previously revealed energy dissipation mechanism in the fracture of BMGs.

Stabilizing to disruptive transition of focal adhesion response to mechanical forces

Strong mechanical forces can, obviously, disrupt cell-cell and cell-matrix adhesions, e.g., cyclic uniaxial stretch induces instability of cell adhesion, which then causes the reorientation of cells away from the stretching direction. However, recent experiments also demonstrated the existence of force dependent adhesion growth (rather than dissociation). To provide a quantitative explanation for the two seemingly contradictory phenomena, a microscopic model that includes both integrin-integrin interaction and integrin-ligand interaction is developed at molecular level by treating the focal adhesion as an adhesion cluster. The integrin clustering dynamics and integrin-ligand binding dynamics are then simulated within one unified theoretical frame with Monte Carlo simulation. We find that the focal adhesion will grow when the traction force is higher than a relative small threshold value, and the growth is dominated by the reduction of local chemical potential energy by the traction force. In contrast, the focal adhesion will rupture when the traction force exceeds a second threshold value, and the rupture is dominated by the breaking of integrin-ligand bonds. Consistent with the experiments, these results suggest a force map for various responses of cell adhesion to different scales of mechanical force. PMID: 20542514

Complex Transition of Double-Diffusive Convection in a Rectangular Enclosure with Height-to-Length Ratio Equal to 4: Part I

This is the first part of direct numerical simulation (DNS) of double-diffusive convection in a slim rectangular enclosure with horizontal temperature and concentration gradients. We consider the case with the thermal Rayleigh number of 10^5, the Pradtle number of 1, the Lewis number of 2, the buoyancy ratio of composition to temperature being in the range of [0,1], and height-to-width aspect ration of 4. A new 7th order upwind compact scheme was developed for approximation of convective terms, and a three-stage third-order Runge-Kutta method was employed for time advancement. Our DNS suggests that with the buoyancy ratio increasing form 0 to 1, the flow of transition is a complex series changing fromthe steady to periodic, chaotic, periodic, quasi-periodic, and finally back to periodic. There are two types of periodic flow, one is simple periodic flow with single fundamental frequency (FF), and another is complex periodic flow with multiple FFs. This process is illustrated by using time-velocity histories, Fourier frequency spectrum analysis and the phase-space rajectories.

Direct numerical simulation of hypersonic boundary layer transition over a blunt cone with a small angle of attack

The direct numerical simulation of boundary layer transition over a 5° half-cone-angle blunt cone is performed. The free-stream Mach number is 6 and the angle of attack is 1°. Random wall blow-and-suction perturbations are used to trigger the transition. Different from the authors’ previous work [Li et al., AIAA J. 46, 2899(2008)], the whole boundary layer flow over the cone is simulated (while in the author’s previous work, only two 45° regions around the leeward and the windward sections are simulated). The transition location on the cone surface is determined through the rapid increase in skin fraction coefficient (Cf). The transition line on the cone surface shows a nonmonotonic curve and the transition is delayed in the range of 0° ≤ θ ≤ 30° (θ = 0° is the leeward section). The mechanism of the delayed transition is studied by using joint frequency spectrum analysis and linear stability theory (LST). It is shown that the growth rates of unstable waves of the second mode are suppressed in the range of 20° ≤ θ ≤ 30°, which leads to the delayed transition location. Very low frequency waves VLFWs� are found in the time series recorded just before the transition location, and the periodic times of VLFWs are about one order larger than those of ordinary Mack second mode waves. Band-pass filter is used to analyze the low frequency waves, and they are deemed as the effect of large scale nonlinear perturbations triggered by LST waves when they are strong enough.The direct numerical simulation of boundary layer transition over a 5° half-cone-angle blunt cone is performed. The free-stream Mach number is 6 and the angle of attack is 1°. Random wall blow-and-suction perturbations are used to trigger the transition. Different from the authors’ previous work [ Li et al., AIAA J. 46, 2899 (2008) ], the whole boundary layer flow over the cone is simulated (while in the author’s previous work, only two 45° regions around the leeward and the windward sections are simulated). The transition location on the cone surface is determined through the rapid increase in skin fraction coefficient (Cf). The transition line on the cone surface shows a nonmonotonic curve and the transition is delayed in the range of 20° ≤ θ ≤ 30° (θ = 0° is the leeward section). The mechanism of the delayed transition is studied by using joint frequency spectrum analysis and linear stability theory (LST). It is shown that the growth rates of unstable waves of the second mode are suppressed in the range of 20° ≤ θ ≤ 30°, which leads to the delayed transition location. Very low frequency waves (VLFWs) are found in the time series recorded just before the transition location, and the periodic times of VLFWs are about one order larger than those of ordinary Mack second mode waves. Band-pass filter is used to analyze the low frequency waves, and they are deemed as the effect of large scale nonlinear perturbations triggered by LST waves when they are strong enough.

Displacive To Order-Disorder Two-Step Phase Transition Model For Para-Ferroelectric Transition

Two-step phase transition model, displacive to order-disorder, is proposed. The driving forces for these two transitions are fundamentally different. The displacive phase transition is one type of the structural phase transitions. We clearly define the structural phase transition as the symmetry broking of the unit cell and the electric dipole starts to form in the unit cell. Then the dipole-dipole interaction takes place as soon as the dipoles in unit cells are formed. We believe that the dipole-dipole interaction may cause an order-disorder phase transition following the displacive phase transition. Both structural and order-disorder phase transition can be first-order or second-order or in between. We found that the structural transition temperatures can be lower or equal or higher than the order-disorder transition temperature. The para-ferroelectric phase transition is the combination of the displacive and order-disorder phase transitions. It generates a variety of transition configurations along with confusions. In this paper, we discuss all these configurations using our displacive to order-disorder two-step phase transition model and clarified all the confusions.

Numerical investigation on a transition prediction model

A new transition prediction model is introduced, which couples the intermittency effect into the turbulence transport equations and takes the characteristics of fluid transition into consideration to mimic the exact process of transition. Test cases include a two-dimensional incompressible plate and a two-dimensional NACA0012 airfoil. Performance of this transition model for incompressible flows is studied, with numerical results consistent to experimental data. The requirement of grid resolution for this transition model is also studied.

Space-resolved spectral diagnosis of line-shaped laser-produced magnesium plasmas

Space-resolved spectra of line-shaped laser-produced magnesium plasmas in the normal direction of the target have been obtained using a pinhole crystal spectrograph. These spectra are treated by a spectrum analyzing code for obtaining the true spectra and fine structures of overlapped lines. The spatial distributions of electron temperature and density along the normal direction of the target surface have been obtained with different spectral diagnostic techniques. Especially, the electron density plateaus beyond the critical surface in line-shaped magnesium plasmas have been obtained with a fitting technique applied to the Stark-broadened Ly-alpha wings of hydrogenic ions. The difference of plasma parameters between those obtained by different diagnostic techniques is discussed. Other phenomena, such as plasma satellites, population inversion, etc., which are observed in magnesium plasmas, are also presented.

Spectral effects for an ultrashort pulse train propagating in a two-level atom medium

We investigate the spectra of a femtosecond pulse train propagating in a resonant two-level atom (TLA) medium. it is found that higher spectral components can be produced even for a 2 pi femtosecond pulse train. Furthermore, the spectral effects depend crucially on both the relative shift phi and the delay time tau between the successive pulses of the femtosecond pulse train.

Constrained high-order statistical blind deconvolution of spectral data

A constrained high-order statistical algorithm is proposed to blindly deconvolute the measured spectral data and estimate the response function of the instruments simultaneously. In this algorithm, no prior-knowledge is necessary except a proper length of the unit-impulse response. This length can be easily set to be the width of the narrowest spectral line by observing the measured data. The feasibility of this method has been demonstrated experimentally by the measured Raman and absorption spectral data.

Conical emission by femtosecond pulses with different spectral bandwidths

Conical emission (CE) has been investigated experimentally by laser pulses with different pulse durations and spectral bandwidths. The results show that the overall CE curve will shift as the varying of spectral bandwidth of pump laser pulse. But for pump laser pulses which have same spectral bandwidth but different pulse duration, the CE angles will be same at the spectral region close to the pump wavelength while will be different at the spectral region far away from the pump wavelength. We have also fitted the measured CE angles with X-wave model. The calculated curves and the measured CE curves match reasonably well. The best fits indicate that the group velocity of the filament pulse may be greatly controlled by controlling the spectral bandwidth of pump laser pulse. (C) 2008 Elsevier B.V. All rights reserved.