Shape effects on the yield stress and deformation of silicon nanowires: A molecular dynamics simulation

The tension and compression of single-crystalline silicon nanowires (SiNWs) with different cross-sectional shapes are studied systematically using molecular dynamics simulation. The shape effects on the yield stresses are characterized. For the same surface to volume ratio, the circular cross-sectional SiNWs are stronger than the square cross-sectional ones under tensile loading, but reverse happens in compressive loading. With the atoms colored by least-squares atomic local shear strain, the deformation processes reveal that the failure modes of incipient yielding are dependent on the loading directions. The SiNWs under tensile loading slip in {111} surfaces, while the compressive loading leads the SiNWs to slip in the {110} surfaces. The present results are expected to contribute to the design of the silicon devices in nanosystems.

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.

Shape Optimization for Hypersonic Arc-Wing Missiles

National Natural Science Foundation of China (NO.90916013)

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.

Self-focusing, Raman, and modulation instability of shaped relativistic laser pulse in plasmas

The interaction of shaped laser pulses with plasmas is studied in a strict theoretical framework without adopting the slow-varying envelope approximation (SVEA). Any physical quantities involved in the interaction are denoted as a summation of different real quantities of respective phases. The relationships among the phases of those real quantities and their moduli are strictly analyzed. Such strict analyses lead to a more exact equation set for the three-dimensional envelope of the laser pulse, which is not based on SVEA. Based on this equation set, self-focusing, Raman, and modulation instabilities could be discussed in a unified framework. The solutions of this equation set for the laser envelope reveal many possible multicolor laser modes in plasmas. The energy and the shape of a pulse determine its propagation through plasmas in a multicolor mode or in a monochromic mode. A global growth rate is introduced to measure the speed of the transition from the monochromic mode in vacuum to a possible mode in plasmas. (c) 2006 American Institute of Physics.

Quantum interference effects of two coherent population trapping states on the atomic spectral lines of a F-e=0 -> F-g=1 transition

We investigate the fluorescence spectrum in a nearly degenerate atomic system of a F-e = 0 -> F-g = 1 transition by analytically solving Schrodinger equations. An ultranarrow fluorescence spectral line in between the two coherent population trapping windows has been found. Our analytic solutions clearly show the origin of the ultranarrow spectral line. Due to quantum interference effects between two coherent population trapping states, the width and intensity of the central spectral line can be controlled by an external magnetic field. Such an effect may be used to detect a magnetic field.

Coherent spectroscopy of the Rb-87 Fg=1 -> F-e=0 closed transition

This paper has observed linewidth narrowing of dark states in rubidium cell by using the Hanle configuration. The reduction of the coherent resonance width under the transition of Rb-87 F-g = 1 -> F-e = 0 is observed and the qualitative explanation about its mechanism is presented. Multiple subnatural width dips are obtained with a linearly polarized laser beam for the transition of Rb-87 F-g = 0, 1, 2. The feature of negative and positive slope, namely dispersionlike feature, is observed in the transmitted light.