A kinematic subgrid scale model for large-eddy simulation of turbulence generated sound

In the hybrid approach of large-eddy simulation (LES) and Lighthill’s acoustic analogy for turbulence-generated sound, the turbulence source fields are obtained using an LES and the turbulence-generated sound at far fields is calculated from Lighthill’s acoustic analogy. As only the velocity fields at resolved scales are available from the LES, the Lighthill stress tensor, serving as a source term in Lighthill’s acoustic equation, has to be evaluated from the resolved velocity fields. As a result, the contribution from the unresolved velocity fields is missing in the conventional LES. The sound of missing scales is shown to be important and hence needs to be modeled. The present study proposes a kinematic subgrid-scale (SGS) model which recasts the unresolved velocity fields into Lighthill’s stress tensors. A kinematic simulation is used to construct the unresolved velocity fields with the imposed temporal statistics, which is consistent with the random sweeping hypothesis. The kinematic SGS model is used to calculate sound power spectra from isotropic turbulence and yields an improved result: the missing portion of the sound power spectra is approximately recovered in the LES.

Subgrid-scale contributions to Lagrangian time correlations in isotropic turbulence

The application of large-eddy simulation (LES) to particle-laden turbulence raises such a fundamental question as whether the LES with a subgrid scale (SGS) model can correctly predict Lagrangian time correlations (LTCs). Most of the currently existing SGS models are constructed based on the energy budget equations. Therefore, they are able to correctly predict energy spectra, but they may not ensure the correct prediction on the LTCs. Previous researches investigated the effect of the SGS modeling on the Eulerian time correlations. This paper is devoted to study the LTCs in LES. A direct numerical simulation (DNS) and the LES with a spectral eddy viscosity model are performed for isotropic turbulence and the LTCs are calculated using the passive vector method. Both a priori and a posteriori tests are carried out. It is observed that the subgrid-scale contributions to the LTCs cannot be simply ignored and the LES overpredicts the LTCs than the DNS. It is concluded from the straining hypothesis that an accurate prediction of enstrophy spectra is most critical to the prediction of the LTCs.

Large-eddy simulation of turbulent collision of heavy particles in isotropic turbulence

The small-scale motions relevant to the collision of heavy particles represent a general challenge to the conventional large-eddy simulation (LES) of turbulent particle-laden flows. As a first step toward addressing this challenge, we examine the capability of the LES method with an eddy viscosity subgrid scale (SGS) model to predict the collision-related statistics such as the particle radial distribution function at contact, the radial relative velocity at contact, and the collision rate for a wide range of particle Stokes numbers. Data from direct numerical simulation (DNS) are used as a benchmark to evaluate the LES using both a priori and a posteriori tests. It is shown that, without the SGS motions, LES cannot accurately predict the particle-pair statistics for heavy particles with small and intermediate Stokes numbers, and a large relative error in collision rate up to 60% may arise when the particle Stokes number is near St_K=0.5. The errors from the filtering operation and the SGS model are evaluated separately using the filtered-DNS (FDNS) and LES flow fields. The errors increase with the filter width and have nonmonotonic variations with the particle Stokes numbers. It is concluded that the error due to filtering dominates the overall error in LES for most particle Stokes numbers. It is found that the overall collision rate can be reasonably predicted by both FDNS and LES for St_K>3. Our analysis suggests that, for St_K<3, a particle SGS model must include the effects of SGS motions on the turbulent collision of heavy particles. The spectral analysis of the concentration fields of the particles with different Stokes numbers further demonstrates the important effects of the small-scale motions on the preferential concentration of the particles with small Stokes numbers.

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.

A novel method to measure coma aberration of projection system

A novel method to measure coma aberration by pattern displacements at different defocus positions is proposed in this paper. The effect of defocus on coma-induced pattern displacement is analyzed. The measuring principle of the method is described in detail. Using the simulation program PROLITH, the proportionality factors between pattern displacement and coma aberration at different defocus positions are calculated. It is proved that the method is simple to perform and the measurement accuracy of coma can increase approximately by 25% by this novel method. (c) 2006 Elsevier GmbH. All rights reserved.

Strong visible up-conversion emissions and thermometric applications of Er3+-Yb3+ codoped Al2O3 prepared by the sol-gel method

The Er3+-Yb3+ codoped Al2O3 has been prepared by the sol-gel method using the aluminium isopropoxide [Al(OC3H7)(3)]-derived Al2O3 sols with addition of the erbium nitrate [Er(NO3)(3) center dot 5H(2)O] and ytterbium nitrate [Yb(NO3)(3) center dot 5H(2)O]. The phase structure, including only two crystalline types of doped Al2O3 phases, theta and gamma, was obtained for the 1 mol% Er3+ and 5 mol% Yb3+ codoped Al2O3 at the sintering temperature of 1,273 K. By a 978 nm semiconductor laser diodes excitation, the visible up-conversion emissions centered at about 523, 545, and 660 nm were obtained. The temperature dependence of the green up-conversion emissions was studied over a wide temperature range of 300-825 K, and the reasonable agreement between the calculated temperature by the fluorescence intensity ratio (FIR) theory and the measured temperature proved that Er3+-Yb3+ codoped Al2O3 plays an important role in the application of high temperature sensor.

A new method of modeling cladding band structure of air-guiding photonic crystal fibers

Cladding band structure of air-guiding photonic crystal fibers with high air-filling fraction is calculated in terms of fiber shape variation. The fundamental photonic band gap dependence on structure parameters, air-filling fraction and spacing, is also investigated. The numerical results show that the band gap edges shift toward longer wavelength as the air-filling fraction is increased, whereas the relative band gap width increases linearly. For a fixed air-filling fraction, the band gap edges with respect to spacing keep constant. With this method, the simulation results agree well with the reported data. © 2007 Elsevier B.V. All rights reserved.

Gamma-irradiation effects on Ce-doped YAG crystals grown by Cz and TGT method

Gamma-rays radiation effects on Ce:YAG crystals grown by Czochralski (Cz) and temperature gradient techniques (TGT) have been studied by means of optical absorption and luminescence spectra. Valence of Ce3+ ion changes during the gamma-ray irradiation process and this result indicates Ce4+ ion may exist in both Cz-Ce:YAG and TGT-Ce:YAG crystals. Thermally stimulated luminescence measurements reveal intense thermoluminescence peaks in gamma-irradiated Ce:YAG crystals and trap parameters were calculated by general-order kinetics expression. (C) 2006 Elsevier B.V. All rights reserved.

Estimation of the detection probability for Yangtze finless porpoises (Neophocaena phocaenoides asiaeorientalis) with a passive acoustic method

Yangtze finless porpoises were surveyed by using simultaneous visual and acoustical methods from 6 November to 13 December 2006. Two research vessels towed stereo acoustic data loggers, which were used to store the intensity and sound source direction of the high frequency sonar signals produced by finless porpoises at detection ranges up to 300 m on each side of the vessel. Simple stereo beam forming allowed the separation of distinct biosonar sound source, which enabled us to count the number of vocalizing porpoises. Acoustically, 204 porpoises were detected from one vessel and 199 from the other vessel in the same section of the Yangtze River. Visually, 163 and 162 porpoises were detected from two vessels within 300 m of the vessel track. The calculated detection probability using acoustic method was approximately twice that for visual detection for each vessel. The difference in detection probabilities between the two methods was caused by the large number of single individuals that were missed by visual observers. However, the sizes of large groups were underestimated by using the acoustic methods. Acoustic and visual observations complemented each other in the accurate detection of porpoises. The use of simple, relatively inexpensive acoustic monitoring systems should enhance population surveys of free-ranging, echolocating odontocetes. (C) 2008 Acoustical Society of America.

Estimated detection distance of a baiji's (Chinese river dolphin, Lipotes vexillifer) whistles using a passive acoustic survey method

Source levels and phonation intervals of whistles produced by a free-ranging baiji (Chinese river dolphin) were measured in the seminatural reserve of Shishou in Hubei, China. A total of 43 whistles were recorded over 12 recording sessions. The mean dominant frequency (the frequency at the highest energy) was 5.7 kHz (s.d.=0.67). The calculated source level was 143.2 dB rms re 1 mu Pa (s.d.=5.8). Most phonation intervals were shorter than 460 s, and the average interval was 205 s (s.d. = 254). Theoretical detection range of baiji's whistle was 6600 m at the present study site, but it could reduce a couple of hundred meters in practical noisy situation in the Yangtze River. Sporadic phonation (205 s interval on average) with relatively faint signal of baiji was considered, to be difficult to be detected by a towing hydrophone system. Stationed monitoring or slow speed towing of hydrophones along the river current is recommended. (c) 2006 Acoustical Society of America.

Research on the influence of fabrication parameters on the performance of buried ion-exchanged glass waveguides using the variational method

The variational method is proposed to analyze the influence of the fabrication parameters on the performance of buried K+-Na+ ion-exchanged Er3+-Yb3+ ions co-doped glass waveguide. The unknown parameters of the Hermite-Gaussian functions as the trial field distribution are determined based on the scalar variational principle. It is demonstrated that the results calculated in this paper agree with those measured in the experiment. The mode dimensions, the effective refractive index, and the overlap factor as the functions of the fabrication parameters are investigated. These results of the variational analysis are useful for the design and optimization of Er3+-Yb3+ ions co-doped waveguides.

CALCULATED ELECTRONIC AND MAGNETIC-STRUCTURES OF THE NEW TERNARY RARE-EARTH IRON NITRIDE ND2FE17N3

The electronic and magnetic structures of Nd2Fe17 and Nd2Fe17N3 have been calculated using the first-principle, spin-polarized orthogonalized linear combination of atomic orbitals method. Comparative studies of the two materials reveal important effects of the nitrogen atoms (at 9e site) on the electronic and magnetic structures. Results are presented for the total density of states, site-projected partial density of states and the spin magnetic moments on four nonequivalent Fe sites. The highest magnetic moments are found to be located on the 6c site for Nd2Fe17 and on the 9d site for Nd2Fe17N3, in agreement with the neutron and Mossbauer experiments. The variation trends of the magnetic moments on different Fe sites are discussed in terms of the separation between Fe and N atoms. Compared with Nd2Fe17, an increase in the exchange splitting of the Fe d band is found in Nd2Fe17N3, which accounts for its higher Curie temperature as observed in experiments. The calculated results show that the nitrogen atoms are charge acceptors in these compounds.

CALCULATION OF DIELECTRIC-CONSTANTS FOR SEMICONDUCTORS - AN APPLICATION OF ABINITIO PSEUDOPOTENTIAL METHOD

We have used ab initio pseudopotential method to generate basis wavefunctions and eigen energies to carry out first principle calculations of the static macroscopic dielectric constant for GaAs and GaP. The resulted converged random phase approximation (RPA) value is 12.55 and 10.71, in excellent agreement to the experimental value of 12.4 and 10.86, respectively. The inclusion of the exchange correlation contribution makes the calculated result slightly worsen. A convergence test with respect to the number of k points in Brillouin zone (BZ) integration was carried out. Sixty irreducible BZ k points were used to achieve the converged results. Integration with only 10 special k points increased the RPA value by 15%.

The scattering matrix method for quantum waveguides

A scattering matrix method for investigating the electron transport in quantum waveguides is presented. By dividing the structure into a number of transverse slices, the global scattering matrix is obtained by the composition of the individual scattering matrices associated with each interface. Complicated geometries and inhomogeneous external potentials are included in the formulation. It is shown that the proposed scattering matrix method possesses many advantages over the traditional mode-matching and transfer matrix methods, especially in treating the electron wave propagation in complicated geometries. Justification for the method is provided by the unitarity of the calculated scattering matrix, and the consistency of the results with those obtained by the recursive Green's function method.

A Numerical Method for Calculating Transmission Coefficients Across Arbitrary Potential Barriers with High Accuracy

We report a new method for calculating transmission coefficients across arbitrary potential barriers based on the Runge-Kutta method. A numerical solution of the Schrodinger equation is calculated using the Runge-Kutta method,and a new model is established to analyze the numerical results to find the transmission coefficient. This technique is applied to various cases, such as parabolic potential barrier and double-barrier structures. Transmission probability with high precision is obtained and discussed. The tunnelling current density through a MOS structure is also explored and the result coincides with the Fowler-Nordheim model,which indicates the applicability of our method.