Computation of Viscous Incompressible Flow using Pressure Correction Method on Unstructured Chimera Grid

In this paper, a pressure correction algorithm for computing incompressible flows is modified and implemented on unstructured Chimera grid. Schwarz method is used to couple the solutions of different sub-domains. A new interpolation to ensure consistency between primary variables and auxiliary variables is proposed. Other important issues such as global mass conservation and order of accuracy in the interpolations are also discussed. Two numerical simulations are successfully performed. They include one steady case, the lid-driven cavity and one unsteady case, the flow around a circular cylinder. The results demonstrate a very good performance of the proposed scheme on unstructured Chimera grids. It prevents the decoupling of pressure field in the overlapping region and requires only little modification to the existing unstructured Navier–Stokes (NS) solver. The numerical experiments show the reliability and potential of this method in applying to practical problems.

Perturbational Finite Volume Method For The Solution of 2-D Navier-Stokes Equations on Unstructured and Structured Colocated Meshes

Based on the first-order upwind and second-order central type of finite volume( UFV and CFV) scheme, upwind and central type of perturbation finite volume ( UPFV and CPFV) schemes of the Navier-Stokes equations were developed. In PFV method, the mass fluxes of across the cell faces of the control volume (CV) were expanded into power series of the grid spacing and the coefficients of the power series were determined by means of the conservation equation itself. The UPFV and CPFV scheme respectively uses the same nodes and expressions as those of the normal first-order upwind and second-order central scheme, which is apt to programming. The results of numerical experiments about the flow in a lid-driven cavity and the problem of transport of a scalar quantity in a known velocity field show that compared to the first-order UFV and second-order CFV schemes, upwind PFV scheme is higher accuracy and resolution, especially better robustness. The numerical computation to flow in a lid-driven cavity shows that the under-relaxation factor can be arbitrarily selected ranging from 0.3 to 0. 8 and convergence perform excellent with Reynolds number variation from 102 to 104.

Development of an efficient particle approach for micro-scale gas flow simulations

The micro-scale gas flows are usually low-speed flows and exhibit rarefied gas effects. It is challenging to simulate these flows because traditional CFD method is unable to capture the rarefied gas effects and the direct simulation Monte Carlo (DSMC) method is very inefficient for low-speed flows. In this study we combine two techniques to improve the efficiency of the DSMC method. The information preservation technique is used to reduce the statistical noise and the cell-size relaxed technique is employed to increase the effective cell size. The new cell-size relaxed IP method is found capable of simulating micro-scale gas flows as shown by the 2D lid-driven cavity flows.

Numerical Investigation of Richtmyer-Meshkov Instability Driven by Cylindrical Shocks

In this paper, a numerical method with high order accuracy and high resolution was developed to simulate the Richtmyer-Meshkov(RM) instability driven by cylindrical shock waves. Compressible Euler equations in cylindrical coordinate were adopted for the cylindrical geometry and a third order accurate group control scheme was adopted to discretize the equations. Moreover, an adaptive grid technique was developed to refine the grid near the moving interface to improve the resolution of numerical solutions. The results of simulation exhibited the evolution process of RM instability, and the effect of Atwood number was studied. The larger the absolute value of Atwood number, the larger the perturbation amplitude. The nonlinear effect manifests more evidently in cylindrical geometry. The shock reflected from the pole center accelerates the interface for the second time, considerably complicating the interface evolution process, and such phenomena of reshock and secondary shock were studied.

A Nanomechanical Device Based on Light-Driven Proton Pumps

In this paper, a hybrid device based on a microcantilever interfaced with bacteriorhodopsin (bR) is constructed. The microcantilever, on which the highly oriented bR film is self-assembled, undergoes controllable and reversible bending when the light-driven proton pump protein, bR, on the microcantilever surface is activated by visible light. Several control experiments are carried out to preclude the influence of heat and photothermal effects. It is shown that the nanomechanical motion is induced by the resulting gradient of protons, which are transported from the KCl solution on the cytoplasmic side of the bR film towards the extracellular side of the bR film. Along with a simple physical interpretation, the microfabricated cantilever interfaced with the organized molecular film of bR can simulate the natural machinery in converting solar energy to mechanical energy.

Oxyhydrogen combustion and detonation driven shock tube

The performance of combustion driver ignited by multi-spark plugs distributed along axial direction has been analysed and tested. An improved ignition method with three circumferential equidistributed ignitors at main diaphragm has been presented, by which the produced incident shock waves have higher repeatability, and better steadiness in the pressure, temperature and velocity fields of flow behind the incident shock, and thus meets the requirements of aerodynamic experiment. The attachment of a damping section at the end of the driver can eliminate the high reflection pressure produced by detonation wave, and the backward detonation driver can be employed to generate high enthalpy and high density test flow. The incident shock wave produced by this method is well repeated and with weak attenuation. The reflection wave caused by the contracted section at the main diaphragm will weaken the unfavorable effect of rarefaction wave behind the detonation wave, which indicates that the forward detonation driver can be applied in the practice. For incident shock wave of identical strength, the initial pressure of the forward detonation driver is about 1 order of magnitude lower than that of backward detonation.

Flow Characteristics of Liquids in Microtubes Driven by a High Pressure

The flow characteristics of liquids in microtubes driven by a high pressure ranging from 1 MPa to 30 MPa are studied in this paper. The diameter of the microtube is from 3 μm to 10 μm and liquids composed of simple small molecules are chosen as the working fluids. The Reynolds number ranges from 0. 1 to 24. The behavior of isopropanol and carbon tetrachloride under high pressure is found different from the prediction from conventional Hagen-Poiseuille (HP) equation. The normalized friction coefficient C* increases significantly with the pressure. From an analysis of the microtube deformation, liquid compressibility, viscous heating and wall slip, it may be seen that the viscosity at high pressure plays an important role here. An exponential function of viscosity vs pressure is introduced into the HP equation to counteract the difference between experimental and theoretical values. However, this difference is not so marked for di-water.

Upwind compact method and direct numerical-simulation of driven flow in a square cavity

A high-order accurate finite-difference scheme, the upwind compact method, is proposed. The 2-D unsteady incompressible Navier-Stokes equations are solved in primitive variables. The nonlinear convection terms in the governing equations are approximated by using upwind biased compact difference, and other spatial derivative terms are discretized by using the fourth-order compact difference. The upwind compact method is used to solve the driven flow in a square cavity. Solutions are obtained for Reynolds numbers as high as 10000. When Re less than or equal to 5000, the results agree well with those in literature. When Re = 7500 and Re = 10000, there is no convergence to a steady laminar solution, and the flow becomes unsteady and periodic.

Failure process analysis of cement under explosive loading with a generalized driven nonlinear threshold model

The microstructural heterogeneity and stress fluctuation play important roles in the failure process of brittle materials. In this paper, a generalized driven nonlinear threshold model with stress fluctuation is presented to study the effects of microstructural heterogeneity on continuum damage evolution. As an illustration, the failure process of cement material under explosive loading is analyzed using the model. The result agrees well with the experimental one, which proves the efficiency of the model.

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.

Influence of Contact Angle and Tube Size on Capillary-Driven Flow Under Microgravity

The influence of contact angle and tube radius on the capillary-driven flow for circular cylindrical tubes is studied systematically by microgravity experiments using the drop tower. Experimental results show that the velocity of the capillary flow decreases monotonically with an increase in the contact angle. However, the time-evolution of the velocity of the capillary flow is different for different sized tubes. At the beginning of the microgravity period, the capillary flow in a thinner tube moves faster than that in a thicker tube, and then the latter overtakes the former. Therefore, there is an intersection between the curves of meniscus velocity vs microgravity time for two differently sized tubes. In addition, for two given sized tubes this intersection is delayed when the contact angle increases. The experimental results are analyzed theoretically and also supported by numerical computations.

Electron dynamics and harmonics emission spectra due to electron oscillation driven by intense laser pulses

The dynamics and harmonics emission spectra due to electron oscillation driven by intense laser pulses have been investigated considering a single electron model. The spectral and angular distributions of the harmonics radiation are numerically analyzed and demonstrate significantly different characteristics from those of the low-intensity field case. Higher-order harmonic radiation is possible for a sufficiently intense driving laser pulse. A complex shifting and broadening structure of the spectrum is observed and analyzed for different polarization. For a realistic pulsed photon beam, the spectrum of the radiation is redshifted for backward radiation and blueshifted for forward radiation, and spectral broadening is noticed. This is due to the changes in the longitudinal velocity of the electron during the laser pulse. These effects are much more pronounced at higher laser intensities giving rise to even higher-order harmonics that eventually leads to a continuous spectrum. Numerical simulations have further shown that broadening of the high harmonic radiation can be limited by increasing the laser pulse width. The complex shifting and broadening of the spectra can be employed to characterize the ultrashort and ultraintense laser pulses and to study the ultrafast dynamics of the electrons. (c) 2006 American Institute of Physics.

Energetic protons emitted from coulomb explosion dynamics of large-sized hydrogen clusters driven by an ultrashort intense laser pulse

Protons with very high kinetic energy of about 10keV and the saturation effect of proton energy for laser intensity have been observed in the interaction of an ultrashort intense laser pulse with large-sized hydrogen clusters. Including the cluster-size distribution as well as the laser-intensity distribution on the focus spot, the theoretical calculations based on a simplified Coulomb explosion model have been compared with our experimental measurements, which are in good agreement with each other.

Electron quantum path tuning and isolated attosecond pulse emission driven by a waveform-controlled multi-cycle laser field

We investigate high-order harmonic emission and isolated attosecond pulse (IAP) generation in atoms driven by a two-colour multi-cycle laser field consisting of an 800 nm pulse and an infrared laser pulse at an arbitrary wavelength. With moderate laser intensity, an IAP of similar to 220 as can be generated in helium atoms by using two-colour laser pulses of 35 fs/800 nm and 46 fs/1150 nm. The discussion based on the three-step semiclassical model, and time-frequency analysis shows a clear picture of the high-order harmonic generation in the waveform-controlled laser field which is of benefit to the generation of XUV IAP and attosecond electron pulses. When the propagation effect is included, the duration of the IAP can be shorter than 200 as, when the driving laser pulses are focused 1 mm before the gas medium with a length between 1.5 mm and 2 mm.

Laser-driven coherent betatron oscillation in a laser-wakefield cavity

The origin of beam disparity in emittance and betatron oscillation orbits, in and out of the polarization plane of the drive laser of laser-plasma accelerators, is explained in terms of betatron oscillations driven by the laser field. As trapped electrons accelerate, they move forward and interact with the laser pulse. For the bubble regime, a simple model is presented to describe this interaction in terms of a harmonic oscillator with a driving force from the laser and a restoring force from the plasma wake field. The resulting beam oscillations in the polarization plane, with period approximately the wavelength of the driving laser, increase emittance in that plane and cause microbunching of the beam. These effects are observed directly in 3D particle-in-cell simulations.