Subgrid scale fluid velocity timescales seen by inertial particles in large-eddy simulation of particle-laden turbulence

Large-eddy simulation (LES) has emerged as a promising tool for simulating turbulent flows in general and, in recent years,has also been applied to the particle-laden turbulence with some success (Kassinos et al., 2007). The motion of inertial particles is much more complicated than fluid elements, and therefore, LES of turbulent flow laden with inertial particles encounters new challenges. In the conventional LES, only large-scale eddies are explicitly resolved and the effects of unresolved, small or subgrid scale (SGS) eddies on the large-scale eddies are modeled. The SGS turbulent flow field is not available. The effects of SGS turbulent velocity field on particle motion have been studied by Wang and Squires (1996), Armenio et al. (1999), Yamamoto et al. (2001), Shotorban and Mashayek (2006a,b), Fede and Simonin (2006), Berrouk et al. (2007), Bini and Jones (2008), and Pozorski and Apte (2009), amongst others. One contemporary method to include the effects of SGS eddies on inertial particle motions is to introduce a stochastic differential equation (SDE), that is, a Langevin stochastic equation to model the SGS fluid velocity seen by inertial particles (Fede et al., 2006; Shotorban and Mashayek, 2006a; Shotorban and Mashayek, 2006b; Berrouk et al., 2007; Bini and Jones, 2008; Pozorski and Apte, 2009).However, the accuracy of such a Langevin equation model depends primarily on the prescription of the SGS fluid velocity autocorrelation time seen by an inertial particle or the inertial particle–SGS eddy interaction timescale (denoted by $\delt T_{Lp}$ and a second model constant in the diffusion term which controls the intensity of the random force received by an inertial particle (denoted by C_0, see Eq. (7)). From the theoretical point of view, dTLp differs significantly from the Lagrangian fluid velocity correlation time (Reeks, 1977; Wang and Stock, 1993), and this carries the essential nonlinearity in the statistical modeling of particle motion. dTLp and C0 may depend on the filter width and particle Stokes number even for a given turbulent flow. In previous studies, dTLp is modeled either by the fluid SGS Lagrangian timescale (Fede et al., 2006; Shotorban and Mashayek, 2006b; Pozorski and Apte, 2009; Bini and Jones, 2008) or by a simple extension of the timescale obtained from the full flow field (Berrouk et al., 2007). In this work, we shall study the subtle and on-monotonic dependence of $\delt T_{Lp}$ on the filter width and particle Stokes number using a flow field obtained from Direct Numerical Simulation (DNS). We then propose an empirical closure model for $\delta T_{Lp}$. Finally, the model is validated against LES of particle-laden turbulence in predicting single-particle statistics such as particle kinetic energy. As a first step, we consider the particle motion under the one-way coupling assumption in isotropic turbulent flow and neglect the gravitational settling effect. The one-way coupling assumption is only valid for low particle mass loading.

Large-eddy simulation of flows past a flapping airfoil using immersed boundary method

The numerical simulation of flows past flapping foils at moderate Reynolds numbers presents two challenges to computational fluid dynamics: turbulent flows and moving boundaries. The direct forcing immersed boundary (IB) method has been devel- oped to simulate laminar flows. However, its performance in simulating turbulent flows and transitional flows with moving boundaries has not been fully evaluated. In the present work, we use the IB method to simulate fully developed turbulent channel flows and transitional flows past a stationary/plunging SD7003 airfoil. To suppress the non-physical force oscillations in the plunging case, we use the smoothed discrete delta function for interpolation in the IB method. The results of the present work demonstrate that the IB method can be used to simulate turbulent flows and transitional flows with moving boundaries.

New Numerical Algorithms in SUPER CE/SE and Their Applications in Explosion Mechanics

The new numerical algorithms in SUPER/CESE and their applications in explosion mechanics are studied. The researched algorithms and models include an improved CE/SE (space-time Conservation Element and Solution Element) method, a local hybrid particle level set method, three chemical reaction models and a two-fluid model. Problems of shock wave reflection over wedges, explosive welding, cellular structure of gaseous detonations and two-phase detonations in the gas-droplet system are simulated by using the above-mentioned algorithms and models. The numerical results reveal that the adopted algorithms have many advantages such as high numerical accuracy, wide application field and good compatibility. The numerical algorithms presented in this paper may be applied to the numerical research of explosion mechanics.

NUMERICAL AND EXPERIMENTAL STUDY ON LIQUID-SOLID FLOW IN A HYDROCYCLONE

Hydrocyclones are widely used in industry, of which the geometrical design using CFD techniques is gaining more popularity in recent years. In this study, the Euler-Euler approach and the Reynolds stress model are applied to simulate the liquid-solid flowfield in a hydrocyclone. The methodology is validated by a good agreement between experimental data and numerical results. Within the research range, the simulation indicates that the liquid-solid separation mainly occurs in the conical segment, and increasing conical height or decreasing cylindrical height helps to improve the grade efficiencies of solid particles. Based on these results, two of the same hydrocyclones are designed and installed in series to establish a liquid-solid separation system. Many experiments are then conducted under different conditions, in which the effects of the water cut and the second hydrocyclone on the separation are investigated. The results also confirm that smaller solid particles are more susceptible to the inlet conditions, and the second hydrocyclone plays a more important role as the water cut reduces.

A numerical model for predicting the jump of lift force on an oscillating cylinder

The fluid force coefficients on a transversely oscillating cylinder are calculated by applying two- dimensional large eddy simulation method. Considering the ‘‘jump’’ phenomenon of the amplitude of lift coefficient is harmful to the security of the submarine slender structures, the characteristics of this ‘‘jump’’ are dissertated concretely. By comparing with experiment results, we establish a numerical model for predicting the jump of lift force on an oscillating cylinder, providing consultation for revising the hydrodynamic parameters and checking the fatigue life scale design of submarine slender cylindrical structures.

The Experimental and Numerical Studies on Gas Production from Hydrate Reservoir by Depressurization

A set of experimental system to study hydrate dissociation in porous media is built and some experiments on hydrate dissociation by depressurization are carried out. A mathematical model is developed to simulate the hydrate dissociation by depressurization in hydrate-bearing porous media. The model can be used to analyze the effects of the flow of multiphase fluids, the kinetic process and endothermic process of hydrate dissociation, ice-water phase equilibrium, the variation of permeability, convection and conduction on the hydrate dissociation, and gas and water productions. The numerical results agree well with the experimental results, which validate our mathematical model. For a 3-D hydrate reservoir of Class 3, the evolutions of pressure, temperature, and saturations are elucidated and the effects of some main parameters on gas and water rates are analyzed. Numerical results show that gas can be produced effectively from hydrate reservoir in the first stage of depressurization. Then, methods such as thermal stimulation or inhibitor injection should be considered due to the energy deficiency of formation energy. The numerical results for 3-D hydrate reservoir of Class 1 show that the overlying gas hydrate zone can apparently enhance gas rate and prolong life span of gas reservoir.

Numerical and Experimental Investigation On The Prevention of Co Deflagration

The exhaust gases from industrial furnaces contain a huge amount of heat and chemical enthalpy. However, it is hard to recover this energy since exhaust gases invariably contain combustible components such as carbon monoxide (CC). If the CO is unexpectedly ignited during the heat recovery process, deflagration or even detonation could occur, with serious consequences such as complete destruction of the equipment. In order to safely utilize the heat energy contained in exhaust gas, danger of its explosion must be fully avoided. The mechanism of gas deflagration and its prevention must therefore be studied. In this paper, we describe a numerical and experimental investigation of the deflagration process in a semi-opened tube. The results show that, upon ignition, a low-pressure wave initially spreads within the tube and then deflagration begins. For the purpose of preventing deflagration, an appropriate amount of nitrogen was injected into the tube at a fixed position. Both simulation and experimental results have shown that the injection of inert gas can successfully interrupt the deflagration process. The peak value of the deflagration pressure can thereby be reduced by around 50%. (C) 2008 Elsevier Ltd. All rights reserved.

Numerical research of heat generation in flashlamp-pumped Nd : glass disk amplifiers

The heat generation in a flashlamp-pumped Nd:glass disk amplifier is studied by the simulation of the whole pumping process, which is based on the ray-tracing method. The results of temperature rise distribution as well as gain distribution are presented. The evolution of heat generation in disk during the pumping process is discussed in detail. Some main factors related with the thermal effect, such as the quantum efficiency, fluorescence lifetime, and pulse duration, are investigated through studying the ratio of the heat generation to energy storage in the gain medium. The influence of each parameter on heat generation is studied carefully, and the results provide ways to decrease the heat generation during the pumping process. (c) 2005 Society of Photo-Optical Instrumentation Engineers.

Simulation of growth process of thin film on non-planar substrate

In (2 + 1) dimension, growth process of thin film on non-planar substrate in Kuramoto-Sivashinsky model is studied with numerical simulation approach. 15 x 15 semi-ellipsoids arranged orderly on the surface of substrate are used to represent initial rough surface. The results show that at the initial stage of growth process, the surface morphology of thin film appears to be grid-structure, and the interface width constantly decreases with the growth time, then reaches minimum. However, the grid-structure becomes ambiguous, and granules of different sizes distribute evenly on the surface of thin film with the increase of growth time. Thereafter, the average size of granules and the interface width gradually increase, and the surface morphology of thin film presents fractal properties. The numerical results of height-height correlation functions of thin film verify the surface morphology of thin film to be fractal for a longer growth time. By fitting of the height-height correlation functions of thin film with different growth times, the growth process is described quantitatively. (c) 2004 Elsevier B.V. All rights reserved.

NUMERICAL ANALYSIS OF A NOVEL 980nm MICROCAVITY SEMICONDUCTOR OPTICAL AMPLIFIER WITH HIGH GAIN AND LOW NOISE

A novel microcavity semiconductor optical amplifier ( MCSOA) was proposed by incorporating top and bottom distributed Bragg reflectors ( DBRs) into the waveguide structure of conventional traveling-wave semiconductor optical amplifiers(TW-SOAs). The incoming( outgoing) light beam incidented onto (escaped from) the waveguide structure at a oblique angle through two optical windows, where the top DBR was etched away, and anti-reflection coating was deposited. The light beams inside the optical cavity were reflected repeatedly between two DBRs and propagated along waveguide in a zigzag optical path. The performance of the MCSOA was systematically investigated by extensive numerical simulation based on a traveling-wave model by taking into account the comprehensive effects of DBRs on both the amplification of signals and the filtering of spontaneous emission( SE). Our results show that the MCSOA is capable of achieving a fiber-to-fiber gain as high as 40dB and a low noise figure is less than 3.5dB.

Numerical study of sub-wavelength plasmonic waveguide

The authors present an analysis of a plasmonic waveguide, simulated using a two-dimensional finite-difference time-domain technique. With the surface structures located on the surface of the metal, the device is able to confine and guide light waves in a sub-wavelength scale. And two waveguides can be placed within 150 nm (similar to 6% of the incident wavelength) that will helpful for the optoelectronic integration. Within the 20 mu m simulation region, it is found that the intensity of the guided light at the interface is roughly two to four times the peak intensity of the incident light, and the propagation length can reach approximately 40 Pm at the wavelength of 2.44 mu m. (c) 2007 Elsevier B.V. All rights reserved.

Numerical analysis for four-wave mixing based wavelength conversion of differential phase-shift keying signals

We numerically investigate the main constrains for high efficiency wavelength conversion of differential phase-shift keying (DPSK) signals based on four-wave mixing (FWM) in highly nonlinear fiber (HNLF). Using multi-tone pump phase modulation techniques, high efficiency wavelength conversion of DPSK signals is achieved with the stimulated Brillouin scattering (SBS) effects effectively suppressed. Our analysis shows that there is a compromise between conversion efficiency and converted idler degradation. By optimizing the pump phase modulation configuration, the converted DPSK idler's degradation can be dramatically decreased through balancing SBS suppression and pump phase modulation degradation. Our simulation results also show that these multi-tone pump phase modulation techniques are more appropriate for the future high bit rate systems.

Numerical and theoretical analysis of the crosstalk in linear optical amplifiers

The dynamic characteristics, including the crosstalk and relaxation oscillation, of linear optical amplifiers (LOAs) are investigated by small-signal analysis under an averaging carrier density approximation and compared with the results of numerical simulation. The good agreement between the numerical simulation and the small-signal analysis indicated the averaging carrier density is an appropriate approximation for analyzing LOAs. Theoretical analyzes also show that the dynamic properties of the vertical laser fields dominate the dynamic performance of LOAs. Based on the small-signal analysis, a concise equation for the crosstalk under high bit rate was derived, which can be applied to measure the differential gain of LOAs.

Experimental and numerical study on AlGaInAs/AlGaAs distributed feedback lasers with GaInP gratings

Ridge-waveguide AlGaInAs/AlGaAs distributed feedback lasers with lattice-matched GaInP gratings were fabricated and their light-current characteristics, spectrum and far-field characteristics were measured. On the basis of our experimental results we analyze the effect of the electron stopper layer on light-current performance using the commercial laser simulation software PICS3D. The simulator is based on the self-consistent solution of drift diffusion equations, the Schrodinger equation, and the photon rate equation. The simulation results suggest that, with the use of a 80 nm-width p-doped Al0.6GaAs electron stopper layer, the slope efficiency can be increased and the threshold current can be reduced by more than 10 mA.

Effects and numerical analysis of argon gas flow on the oxygen concentration in Czochralski silicon single crystal growth

Argon gas, as a protective environment and carrier of latent heat, has an important effect on the temperature distribution in crystals and melts. Numeric simulation is a potent tool for solving engineering problems. In this paper, the relationship between argon gas flow and oxygen concentration in silicon crystals was studied systematically. A flowing stream of argon gas is described by numeric simulation for the first time. Therefore, the results of experiments can be explained, and the optimum argon flow with the lowest oxygen concentration can be achieved. (C) 2002 Elsevier Science B.V. All rights reserved.