基于菲佐干涉仪与多通道光电倍增管阵列的多普勒频移检测技术

提出了基于菲佐干涉仪和多通道光电倍增管(PMT)阵列探测器组合的多普勒频移检测的方案, 适用于风速测量的直接探测多普勒激光雷达。首先介绍了工作原理, 再根据菲佐干涉仪光谱特征对频移检测用干涉仪进行了优化设计, 优化设计的菲佐干涉仪腔长150mm、平板反射率0.755。对提出的菲佐干涉仪和多通道光电倍增管阵列探测器组合的方案进行了数值模拟, 以分子散射作为背景噪声, 计算了该方法的风速测量误差。模拟结果表明, 设计的基于菲佐干涉仪的直接探测多普勒测风激光雷达, 在30 s的积分时间内、探测高度5 km以下,

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.

Mechanisms of upper-ocean thermal variability in a 1970-1988 simulation and observations

EXTRACT (SEE PDF FOR FULL ABSTRACT): Variations in temperature that occurred in the North Pacific thermocline (250 to 400 meters) during the 1970s and 1980s are described in both a numerical simulation and XBT observations.

Influence of a bow of finite width on bowed string motion: numerical modelling and experimental evidence

An enhanced physical model of the bowed string presented previously [1] is explored. It takes into account: the width of the bow, the angular motion of the string, bow-hair elasticity and string bending stiffness. The results of an analytical investigation of a model system - an infinite string sticking to a bow of finite width and driven on one side of the bow - are compared with experimental results published by Cremer [2] and reinterpreted here. Comparison shows that both the width of the bow and the bow-hair elasticity have a large impact on the reflection and transmission behaviour. In general, bending stiffness plays a minor role. Furthermore, a method of numerical simulation of the stiff string bowed with a bow of finite width is presented along with some preliminary results.

CFD-DEM simulation of propagation of sound waves in fluid particles fluidised medium

In this work, speed of sound in 2 phase mixture has been explored using CFD-DEM (Computational Fluid Dynamcis - Discrete Element Modelling). In this method volume averaged Navier Stokes, continuity and energy equations are solved for fluid. Particles are simulated as individual entities; their behaviour is captured by Newton's laws of motion and classical contact mechanics. Particle-fluid interaction is captured using drag laws given in literature.The speed of sound in a medium depends on physical properties. It has been found experimentally that speed of sound drops significantly in 2 phase mixture of fluidised particles because of its increased density relative to gas while maintaining its compressibility. Due to the high rate of heat transfer within 2 phase medium as given in Roy et al. (1990), it has been assumed that the fluidised gas-particle medium is isothermal.The similar phenomenon has been tried to be captured using CFD-DEM numerical simulation. The disturbance is introduced and fundamental frequency in the medium is noted to measure the speed of sound for e.g. organ pipe. It has been found that speed of sound is in agreement with the relationship given in Roy et al. (1990). Their assumption that the system is isothermal also appears to be valid.

Boundary integral technique for the numerical modelling of the air flow for the Aaberg exhaust system

A boundary integral technique has been developed for the numerical simulation of the air flow for the Aaberg exhaust system. For the steady, ideal, irrotational air flow induced by a jet, the air velocity is an analytical function. The solution of the problem is formulated in the form of a boundary integral equation by seeking the solution of a mixed boundary-value problem of an analytical function based on the Riemann-Hilbert technique. The boundary integral equation is numerically solved by converting it into a system of linear algebraic equations, which are solved by the process of the Gaussian elimination. The air velocity vector at any point in the solution domain is then computed from the air velocity on the boundary of the solution domains.

Review of computational aeroacoustics for application in electronics cooler noise

A review of computational aeroacoustics (CCA) was made for application in electronics cooler noise. Computational aeroacoustics encompasses all numerical methods where the purposes is to predict the noise emissions from a simulated flow. Numerical simulation of the flow inside and around heat sinks and fans can lead to a prediction of the emitted noise while they are still in the design phase. Direct CCA is theoretically the best way to predict flow-based acoustic phenomena numerically. It is typically used only for low-frequency sound prediction. The boundary element method offers low computational cost and does not use a computational grid, but instead use vortex-surface calculations to determine tonal noise. Axial fans are commonly used to increase the airflow and thus the heat transfer over the heat sinks within the computer cases. Very detailed source simulations in the fan and heat sink region coupled with the use of analogy methods could result in excellent simulation results with a reasonable computational effort.

Complex chemistry DNS of n-heptane spray autoignition at high pressure and intermediate temperature conditions

Direct Numerical Simulations (DNS) of turbulent n-heptane sprays autoigniting at high pressure (P=24bar) and intermediate air temperature (Tair=1000K) have been performed to investigate the physical mechanisms present under conditions where low-temperature chemistry is expected to be important. The initial turbulence in the carrier gas, the global equivalence ratio in the spray region, and the initial droplet size distribution of the spray were varied. Results show that spray ignition exhibits a spotty nature, with several kernels developing independently in those regions where the mixture fraction is close to its most reactive value ξMR (as determined from homogeneous reactor calculations) and the scalar dissipation rate is low. Turbulence reduces the ignition delay time as it promotes mixing between air and the fuel vapor, eventually resulting in lower values of scalar dissipation. High values of the global equivalence ratio are responsible for a larger number of ignition kernels, due to the higher probability of finding regions where ξ=ξMR. Spray polydispersity results in the occurrence of ignition over a wider range of mixture fraction values. This is a consequence of the inhomogeneities in the mixing field that characterize these sprays, where poorly mixed rich spots are seen to alternate with leaner ones which are well-mixed. The DNS simulations presented in this work have also been used to assess the applicability of the Conditional Moment Closure (CMC) method to the simulation of spray combustion. CMC is found to be a valid method for capturing spray autoignition, although care should be taken in the modelling of the unclosed terms appearing in the CMC equations. © 2013 The Combustion Institute.

Application of scalar dissipation rate model to siemens DLE combustors

The standard design process for the Siemens Industrial Turbomachinery, Lincoln, Dry Low Emissions combustion systems has adopted the Eddy Dissipation Model with Finite Rate Chemistry for reacting computational fluid dynamics simulations. The major drawbacks of this model have been the over-prediction of temperature and lack of species data limiting the applicability of the model. A novel combustion model referred to as the Scalar Dissipation Rate Model has been developed recently based on a flamelet type assumption. Previous attempts to adopt the flamelet philosophy with alternative closure models have failed, with the prediction of unphysical phenomenon. The Scalar Dissipation Rate Model (SDRM) was developed from a physical understanding of scalar dissipation rate, signifying the rate of mixing of hot and cold fluids at scales relevant to sustain combustion, in flames and was validated using direct numerical simulations data and experimental measurements. This paper reports on the first industrial application of the SDRM to SITL DLE combustion system. Previous applications have considered ideally premixed laboratory scale flames. The industrial application differs significantly in the complexity of the geometry, unmixedness and operating pressures. The model was implemented into ANSYS-CFX using their inbuilt command language. Simulations were run transiently using Scale Adaptive Simulation turbulence model, which switches between Large Eddy Simulation and Unsteady Reynolds Averaged Navier Stokes using a blending function. The model was validated in a research SITL DLE combustion system prior to being applied to the actual industrial geometry at real operating conditions. This system consists of the SGT-100 burner with a glass square-sectioned combustor allowing for detailed diagnostics. This paper shows the successful validation of the SDRM against time averaged temperature and velocity within measurement errors. The successful validation allowed application of the SDRM to the SGT-100 twin shaft at the relevant full load conditions. Limited validation data was available due to the complexity of measurement in the real geometry. Comparison of surface temperatures and combustor exit temperature profiles showed an improvement compared to EDM/FRC model. Furthermore, no unphysical phenomena were predicted. This paper presents the successful application of the SDRM to the industrial combustion system. The model shows a marked improvement in the prediction of temperature over the EDM/FRC model previously used. This is of significant importance in the future applications of combustion CFD for understanding of hardware mechanical integrity, combustion emissions and dynamics of the flame. Copyright © 2012 by ASME.

Statistics of reaction progress variable and mixture fraction gradients from DNS of turbulent partially premixed flames

Statistically planar turbulent partially premixed flames for different initial intensities of decaying turbulence have been simulated for global equivalence ratios = 0.7 and 1.0 using three-dimensional, simplified chemistry-based direct numerical simulations (DNS). The simulation parameters are chosen such that the flames represent the thin reaction zones regime combustion. A random bimodal distribution of equivalence ratio is introduced in the unburned gas ahead of the flame to account for the mixture inhomogeneity. The results suggest that the probability density functions (PDFs) of the mixture fraction gradient magnitude |Δξ| (i.e., P(|Δξ|)) can be reasonably approximated using a log-normal distribution. However, this presumed PDF distribution captures only the qualitative nature of the PDF of the reaction progress variable gradient magnitude |Δc| (i.e., P(|Δc|)). It has been found that a bivariate log-normal distribution does not sufficiently capture the quantitative behavior of the joint PDF of |Δξ| and |Δc| (i.e., P(|Δξ|, |Δc|)), and the agreement with the DNS data has been found to be poor in certain regions of the flame brush, particularly toward the burned gas side of the flame brush. Moreover, the variables |Δξ| and |Δc| show appreciable correlation toward the burned gas side of the flame brush. These findings are corroborated further using a DNS data of a lifted jet flame to study the flame geometry dependence of these statistics. © 2013 Copyright Taylor and Francis Group, LLC.

Numerical study on hydrate bearing sediments during gas production

The fully coupled methane hydrate model developed in Cambridge was adopted in this numerical study on gas production trial at the Eastern Nankai Trough, Japan 2013. Based on the latest experimental data of hydrate soil core samples, the clay parameters at Eastern Nankai site were successfully calibrated. With updated clay parameters and site geometry, a 50 days gas production trail was numerically simulated in FLAC2D. The geomechanical behaviour of hydrate bearing sediments under 3 different depressurization strategies were explored and discussed. The results from both axisymmetrical and plane-strain models suggest, the slope of the seabed only affects mechanical properties while no significant impact on the dissociation, temperature and pore pressure. For mechanical deformation after PT recovery, there are large settlements above the perforation zone and small uplift underneath the production zone. To validate the fully coupled model, numerical simulation with finer mesh in the hydrate production zone was carried out. The simulation results suggest good agreement between our model and JOE's results on history matching of gas and water production during trial. Parameter sensitivity of gas production is also investigated and concluded the sea water salinity is a dominant factor for gas production.

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 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 investigations on dissolution and recrystallization processes of GaSb/InSb/GaSb under microgravity and terrestrial conditions

The effects of gravity and crystal orientation on the dissolution of GaSb into InSb melt and the recrystallization of InGaSb were investigated under microgravity condition using a Chinese recoverable satellite and under normal gravity condition on earth. To investigate the effect of gravity on the solid/liquid interface and compositional profiles. a numerical simulation was carried out. The InSb crystal melted at 525 degrees C and then a part of GaSb dissolved into the InSb melt during heating to 706 degrees C and this process led to the formation of InGaSb solution. InGaSb solidified during the cooling process. The experimental and calculation results clearly show that the shape of the solid/liquid interface and compositional profiles in the solution were significantly affected by gravity. Under microgravity, as the Ga compositional profiles were uniform in the radial direction. the interfaces were almost parallel. On the contrary, for normal gravity condition, as large amounts of Ga moved up in the upper region due to buoyancy, the dissolved zone broadened towards gravitational direction. Also. during the cooling process, needle crystals of InGaSb started appearing and the value of x of InxGa1-xSb crystals increased with the decrease of temperature. The GaSb with the (111)B plane dissolved into the InSb melt much more than that of the (111)A plane. (C) 2000 Elsevier Science B.V. All rights reserved.

Calculation of valence subband structures for strained GaInP/AlGaInP quantum wells without axial approximation

Usually in the calculation of valence subband structure for III-V direct bandgap material, axial approximation had been used in the Luttinger-Kohn model to simplify the computational efforts. In this letter, the valence subband structure for the GaInP/AlGaInP strained and lattice-matched quantum wells was calculated without axial approximation, on the basis of 6x6 Luttinger-Kohn Hamiltonian including strain and spin-orbit splitting effects. The numerical simulation results were presented with help of the finite-difference methods. The calculation results with/without axial approximation were compared and the effect of axial approximation on the valence subband structure was discussed in detail. The results indicated that there was a strong warping in the GaInP valence band, and axial approximation can lead to an error when k was not equal to zero, especially for compressively strained and lattice-matched GaInP/AlGaInP quantum wells.