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.

Full and Model-Reduced Structure-Preserving Simulation of Incompressible Fluids

This thesis outlines the construction of several types of structured integrators for incompressible fluids. We first present a vorticity integrator, which is the Hamiltonian counterpart of the existing Lagrangian-based fluid integrator. We next present a model-reduced variational Eulerian integrator for incompressible fluids, which combines the efficiency gains of dimension reduction, the qualitative robustness to coarse spatial and temporal resolutions of geometric integrators, and the simplicity of homogenized boundary conditions on regular grids to deal with arbitrarily-shaped domains with sub-grid accuracy.

Both these numerical methods involve approximating the Lie group of volume-preserving diffeomorphisms by a finite-dimensional Lie-group and then restricting the resulting variational principle by means of a non-holonomic constraint. Advantages and limitations of this discretization method will be outlined. It will be seen that these derivation techniques are unable to yield symplectic integrators, but that energy conservation is easily obtained, as is a discretized version of Kelvin's circulation theorem.

Finally, we outline the basis of a spectral discrete exterior calculus, which may be a useful element in producing structured numerical methods for fluids in the future.

Spectral narrowing of negatively chirped femtosecond pulse by cross-phase modulation in a single-mode optical fiber

We demonstrate theoretically that the negatively chirped femtosecond laser pulse can be spectrally narrowed by cross-phase modulation. The new view is well Supported by numerical simulation. The negative chirp method in fibers might be useful in all optical wavelength switching applications. (c) 2005 Elsevier B.V. All rights reserved.

Improved uniformity of target illumination by combining a lens array and the technique of spectral dispersion

A scheme using a lens array and the technique of spectral dispersion is presented to improve target illumination uniformity in laser produced plasmas. Detailed two-dimensional simulation shows that a quasi-near-field target pattern, of steeper edges and without side lobes, is achieved with a lens array, while interference stripes inside the pattern are smoothed out by the use of the spectral dispersion technique. Moving the target slightly from the exact focal plane of the principal focusing lens can eliminate middle-scale-length intensity fluctuation further. Numerical results indicate that a well-irradiated laser spot with small nonuniformity and great energy efficiency can be obtained in this scheme. (c) 2007 American Institute of Physics.

Beam smoothing by lens array with spectral dispersion

A scheme of combining technology of lens array (LA) and smoothing by spectral dispersion (SSD) is introduced to improve the irradiation uniformity in laser fusion based on the earlier works on LA. The feasibility of the scheme is also analyzed by numerical simulation. It shows that a focal pattern with flat-top and sharp-edge profile could be obtained, and the irradiation nonuniformity can fall down from 14% with only LA to 3% with both SSD and LA. And this smoothing scheme is depended less on the incidence comparing to other smoothing methods. The preliminary experiment has demonstrated its effectiveness. (c) 2006 Elsevier B.V. All rights reserved.

Detached-eddy simulation of transonic flow past a fan-blade section

Detached-eddy simulation of transonic flow past a thin section of a fan blade has been carried out. The inflow Mach number is 1.03, and a bow shock forms upstream of the blade. The shock (corresponding to an adjacent blade) impinges on the suction-side boundary layer which causes separation and rapid transition to turbulence. The boundary layer later re-attaches near the trailing edge. The pressure-side boundary layer transitions near the leading edge and remains attached. Mean surface pressure shows basic agreement with a steady RANS calculation; strong shock motion in the DES is the major cause of discrepancy. Surface pressure spectra are investigated, and low-frequency two-dimensional disturbances associated with the shock motion are dominant. Removing the two-dimensional component from the spectra, the pressure-side three-dimensional spectra reproduce the spectral shape given by a correlation for flat-plate boundary layer wall-pressure spectra developed by Goody. 1 The suction-side disturbances produce similar high- and intermediate-frequency scalings despite substantially different boundary layer development. Near-wake results show that disturbance kinetic energy peaks at the suction-side inflection point of the mean profile, and that the energy is concentrated at low frequencies relative to the near-trailing edge surface pressure. Copyright © 2009 by the authors.

Hole concentration test of p-type GaN by analyzing the spectral response of p-n(+) structure GaN ultraviolet photodetector

A new method to test the hole concentration of p-type GaN is proposed, which is carried out by analyzing the spectral response of p-n(+) structure GaN ultraviolet photodetector. It is shown that the spectral response of the photodetector changes considerably with reversed bias. It is found that the difference between photodetector's quantum efficiency at two wavelengths, i.e. 250 and 361 nm, varies remarkably with increasing reversed bias. According to the simulation calculation, the most characteristic change occurs at a reversed voltage under which the p-GaN layer starts to be completely depleted. Based on this effect the carrier concentration of p-GaN can be derived.

Simulation of light coupling enhancement and localization of transmission field via subwavelength metallic gratings

Surface plasmons(SPs) generated in nano metallic gratings on medium layer can greatly enhance the transmission field through the metallic gratings. The enhancement effect is achieved from lambda = 500 nm to near-infrared domain. The enhancement rate is about 110 % at the wavelength of about 6 10 nm and about 180 % at lambda = 700 nm and 740 nm where most kinds of thin film solar cells have a high spectral response. These structures should provide a promising way to increase the coupling efficiency of thin film solar cells and optical detectors of different wavelength response.

Computer simulation of a-Si : H/mu c-Si : H diphasic silicon solar cells

Based on our experimental research on diphasic silicon films, the parameters such as absorption coefficient, mobility lifetime product and bandgap were estimated by means of effective-medium theory. And then computer simulation of a-Si: H/mu c-Si: H diphasic thin film solar cells was performed. It was shown that the more crystalline fraction in the diphasic silicon films, the higher short circuit density, the lower open-circuit voltage and the lower efficiency. From the spectral response, we can see that the response in long wave region was improved significantly with increasing crystalline fraction in the silicon films. Taking Lambertian back refraction into account, the diphasic silicon films with 40%-50% crystalline fraction was considered to be the best intrinsic layer for the bottom solar cell in micromorph tandem.

Numerical simulation of nc-Si : H/c-Si heterojunction solar cells

AMPS simulator, which was developed by Pennsylvania State University, has been used to simulate photovoltaic performances of nc-Si:H/c-Si solar cells. It is shown that interface states are essential factors prominently influencing open circuit voltages (V-OC) and fill factors (FF) of these structured solar cells. Short circuit current density (J(SC)) or spectral response seems more sensitive to the thickness of intrinsic a-Si:H buffer layers inserted into n(+)-nc-Si:H layer and p-c-Si substrates. Impacts of bandgap offset on solar cell performances have also been analyzed. As DeltaE(C) increases, degradation of VOC and FF owing to interface states are dramatically recovered. This implies that the interface state cannot merely be regarded as carrier recombination centres, and impacts of interfacial layer on devices need further investigation. Theoretical maximum efficiency of up to 31.17% (AM1.5,100mW/cm(2), 0.40-1.1mum) has been obtained with BSF structure, idealized light-trapping effect(R-F=0, R-B=1) and no interface states.

Analysis and numerical simulation of eigenmode characteristics for semiconductor lasers with an equilateral triangle micro-resonator

The eigenmodes confined in the equilateral triangle resonator (ETR) are analyzed by deriving the eigenvalues and the mode field distributions and by the finite difference time domain (FDTD) technique. The analytical results show that the one-period-length for the mode light rays inside the ETR is the perimeter of the ETR, and the number of transverse modes is limited by the condition of total internal reflection. In addition, the sum of the longitudinal mode index and the transverse mode index should be an even number, which limits the number of confined modes again. Based on the FDTD technique and the Pade approximation, we calculate the mode resonant frequencies and the quality factors from the local maximum and the width of the spectral distribution of the intensity The numerical results of mode frequencies agree very well with the analytical results, and the quality factor of the fundamental mode is usually higher than that of the higher order transverse modes. The results show that the ETR is suitable to realize single-made operation as semiconductor microcavity lasers.

Simulation of a monolithically integrated CMOS bioamplifier for EEG recordings

A monolithically integrated CMOS bioamplifier is presented in this paper for EEG recording applications. The capacitive-coupled circuit input structure is utilized to eliminate the large and random DC offsets existing in the electrode-tissue interface. Diode-connected NMOS transistors with negative voltage between gate and source are candidates for large resistors necessary to the bioamplifier. A passive BEF (Band Eliminator Filter) can reduce 50 Hz noise disturbance strength by more than 60 dB. A novel analysis approach is given to help determine the noise power spectral density. Simulation results show that the two-stage CMOS bioamplifier in a closed-loop capacitive feedback configuration,provides an AC in-band gain of 39.6 dB, a DC gain of zero, and an input-referred noise of 87 nVrms integrated from 0.01 Hz to 100 Hz.

Simulation of a Monolithic Integrated CMOS Preamplifier for Neural Recordings

A monolithic integrated CMOS preamplifier is presented for neural recording applications. Two AC-coupied capacitors are used to eliminate the large and random DC offsets existing in the electrode-electrolyte interface. Diode-connected nMOS transistors with a negative voltage between the gate and source are candidates for the large resistors necessary for the preamplifier. A novel analysis is given to determine the noise power spectral density. Simulation results show that the two-stage CMOS preamplifier in a closed-loop capacitive feedback configuration provides an AC in-band gain of 38.8dB,a DC gain of 0,and an input-referred noise of 277nVmax, integrated from 0. 1Hz to 1kHz. The preamplifier can eliminate the DC offset voltage and has low input-referred noise by novel circuit configuration and theoretical analysis.

Radiation-use efficiency and the harvest index of winter wheat at different nitrogen levels and their relationships to canopy spectral reflectance

Radiation-use efficiency (RUE, g/MJ) and the harvest index (HI, unitless) are two helpful characteristics in interpreting crop response to environmental and climatic changes. They are also increasingly important for accurate crop yield simulation, but they are affected by various environmental factors. In this study, the RUE and HI of winter wheat and their relationships to canopy spectral reflectance were investigated based on the massive field measurements of five nitrogen (N) treatments. Crop production can be separated into light interception and RUE. The results indicated that during a long period of slow growth from emergence to regreening, the effect of N on crop production mainly showed up in an increased light interception by the canopy. During the period of rapid growth from regreening to maturity, it was present in both light interception and RUE. The temporal variations of RUEAPAR (aboveground biomass produced per unit of photosynthetically active radiation absorbed by the canopy) during the period from regreening to maturity had different patterns corresponding to the N deficiency, N adequacy and N-excess conditions. Moreover, significant relationships were found between the RUEAPAR and the accumulative normalised difference vegetation index (NDVI) in the integrated season (R-2 = 0.68), between the HI and the accumulative NDVI after anthesis (R-2 = 0.89), and between the RUEgrain (ratio of grain yield to the total amount of photosynthetically active radiation absorbed by the canopy) and the accumulative NDVI of the whole season (R-2 = 0.89) and that after anthesis (R-2 = 0.94). It suggested that canopy spectral reflectance has the potential to reveal the spatial information of the RUEAPAR, HI and RUEgrain. It is hoped that this information will be useful in improving the accuracy of crop yield simulation in large areas.

Numerical simulation and validation of ocean waves measured by an Along-Track Interferometric Synthetic Aperture Radar

A new nonlinear integral transform of ocean wave spectra into Along-Track Interferometric Synthetic Aperture Radar (ATI-SAR) image spectra is described. ATI-SAR phase image spectra are calculated for various sea states and radar configurations based on the nonlinear integral transform. The numerical simulations show that the slant range to velocity ratio (R/V), significant wave height to ocean wavelength ratio (H-s/lambda), the baseline (2B) and incident angle (theta) affect ATI-SAR imaging. The ATI-SAR imaging theory is validated by means of Two X-band, HH-polarized ATI-SAR phase images of ocean waves and eight C-band, HH-polarized ATI-SAR phase image spectra of ocean waves. It is shown that ATI-SAR phase image spectra are in agreement with those calculated by forward mapping in situ directional wave spectra collected simultaneously with available ATI-SAR observations. ATI-SAR spectral correlation coefficients between observed and simulated are greater than 0.6 and are not sensitive to the degree of nonlinearity. However, the ATI-SAR phase image spectral turns towards the range direction, even if the real ocean wave direction is 30 degrees. It is also shown that the ATI-SAR imaging mechanism is significantly affected by the degree of velocity bunching nonlinearity, especially for high values of R/V and H-s/lambda.