Methods for implementing the stream boundary condition in DSMC computations

In the direct simulation Monte‐Carlo (DSMC) method for simulating rarefied gas flows, the velocities of simulator particles that cross a simulation boundary and enter the simulation space are typically generated using the acceptance–rejection procedure that samples the velocities from a truncated theoretical velocity distribution that excludes low and high velocities. This paper analyses an alternative technique, where the velocities of entering particles are obtained by extending the simulation procedures to a region adjacent to the simulation space, and considering the movement of particles generated within that region during the simulation time step. The alternative method may be considered as a form of acceptance–rejection procedure, and permits the generation of all possible velocities, although the population of high velocities is depleted with respect to the theoretical distribution. Nevertheless, this is an improvement over the standard acceptance–rejection method. Previous implementations of the alternative method gave a number flux lower than the theoretical number required. Two methods for obtaining the correct number flux are presented. For upstream boundaries in high‐speed flows, the alternative method is more computationally efficient than the acceptance–rejection method. However, for downstream boundaries, the alternative method is extremely inefficient. The alternative method, with the correct theoretical number flux, should therefore be used in DSMC computations in favour of the acceptance–rejection method for upstream boundaries in high‐speed flows. Copyright © 2003 John Wiley & Sons, Ltd.     

LES of turbulent heat transfer: proper convection numerical schemes for temperature transport

Large eddy simulations of two basic configurations (decay of isotropic turbulence, and the academic plane channel flow) with heat transfer have been performed comparing several convection numerical schemes, in order to discuss their ability to evaluate temperature fluctuations properly. Results are compared with the available incompressible heat transfer direct numerical simulation data. It is shown that the use of regularizing schemes (such as high order upwind type schemes) for the temperature transport equation in combination with centered schemes for momentum transport equation gives better results than the use of centred schemes for both equations. Copyright © 2004 John Wiley & Sons, Ltd.     

Modelling of a single drop impact onto liquid film using particle method

The process of single liquid drop impact on thin liquid surface is numerically simulated with moving particle semi‐implicit method. The mathematical model involves gravity, viscosity and surface tension. The model is validated by the simulation of the experimental cases. It is found that the dynamic processes after impact are sensitive to the liquid pool depth and the initial drop velocity. In the cases that the initial drop velocity is low, the drop will be merged with the liquid pool and no big splash is seen. If the initial drop velocity is high enough, the dynamic process depends on the liquid depth. If the liquid film is very thin, a bowl‐shaped thin crown is formed immediately after the impact. The total crown subsequently expands outward and breaks into many tiny droplets. When the thickness of the liquid film increases, the direction of the liquid crown becomes normal to the surface and the crown propagates outward. It is also found that the radius of the crown is described by a square function of time: r_C = [c(t ? t₀)]^0.5. When the liquid film is thick enough, a crown and a deep cavity inside it are formed shortly after the impact. The bottom of the cavity is initially oblate and then the base grows downward to form a sharp corner and subsequently the corner moves downward. Copyright © 2004 John Wiley & Sons, Ltd.     

电站锅炉系统性能仿真模型的建立

1引言根据热力学第二定律，燃料中的化学拥将近一半损失在锅炉的燃烧与传热过程中。因此，优化锅炉性能对降低整个火电机组的煤耗率具有明显效果。完整的锅炉系统由磨煤机、送风机、一次风机、弓讯机、炉膛、汽包、过热器、再热器、省煤器、空气预热器等单元组成，这些单元之间的联结错综复杂。本文建立了锅炉系统中通用单元的性能模型，并在此基础上，对锅炉系统的结构特点进行了分析，利用过程系统工程中的序贯模块法建立了整个锅炉系统的性能模型。对单元的划分、回路的切断及断裂流股的收敛与迭代策略提出了自己的观点。最后根据生产厂…     

Champ de vitesse dans un échangeur à vortex en régime turbulentHeat exchanger velocity field with turbulent inlet condition

We consider the isothermal flow through a cylindrical flat chamber, a model of some particular heat exchanger, for which LDV measurements and a numerical simulation have been performed. Experimental results show the establishment of an important vortex zone, the secondary flow extending all along the chamber radius. This observation leads to an expected significant increase of the fluid mixing. Results issued from the numerical simulation appear to be in close agreement with experimental data. Nevertheless, the k–ε model used here must be improved to obtain a better approach near the vortex centre. To cite this article: S. Petitot et al., C. R. Mecanique 330 (2002) 593–599.     

Quantitative Comparison of REAPDOR and TRAPDOR Experiments by Numerical Simulations and Determination of H–Al Distances in Zeolites

Quantitative H–Al distances in acid sites of two zeolites with MFI and IFR framework topology were obtained by numerical simulation of ¹H{²⁷Al} rotational echo adiabatic passage double resonance (REAPDOR) experiments. A ²⁷Al offset-dependent data set yields for each resolved ¹H NMR line a corresponding nuclear electric quadrupole coupling constant of the neighboring ²⁷Al site. This information is used for analyzing a second data set for on-resonance irradiation, where the dipolar evolution time (number of rotor cycles) was varied, to yield the ¹H–²⁷Al dipolar coupling constant. Numerical simulations indicate that the REAPDOR method does not depend significantly on the polar angles, defining the orientation of the electric field gradient tensor of ²⁷Al with respect to the Al–H dipolar vector. In contrast, the transfer of populations in double resonance sequence is sensitive to these angles, and it can be thus used to measure them.     

Optimal design of a convergent-barrel cold spray nozzle by numerical method

A convergent-barrel (CB) cold spray nozzle was designed through numerical simulation. It was found that the main factors influencing significantly particle velocity and temperature include the length and diameter of the barrel section, the nature of the accelerating gas and its pressure and temperature, and the particle size. Particles can achieve a relatively low velocity but a high temperature under the same gas pressure using a CB nozzle compared to a convergent-divergent (CD) nozzle. The experiment results with Cu powder using the designed CB nozzle confirmed that particle deposition can be realized under a lower gas pressure with a CB nozzle.     

Performance and simulation of an on-board guided-wave optical processor for airborne SAR applications

This paper describes the performance and simulation of a compact integrated optical processor for the real-time reconstruction of two-dimensional images in airborne stripmap synthetic aperture radar applications. The functional behavior of the processor is explained in some detail. The design criteria are briefly given. The simulation step allowed the main processor characteristics and properties to be identified. A number of comparisons were obtained in airborne SAR mission scenarios between the predictions of the optical device and those achieved by the modern electronic digital approach, based on the wavefront reconstruction method by matched filtering.     

基于二次插值的非饱和土壤水流问题的特征差分方法及数值模拟

本文针对一类非饱和土壤水流问题,提出了基于二次插值的特征差分格式,得到了严谨的L2模误差估计.并作了数值试验,指明方法的有效性.     

Numerical simulation of the spatiotemporal evolution of a nanoparticle–plasma system

A one-dimensional nanodusty plasma was modeled by self-consistently coupling a plasma model with nanoparticle growth, charging, and transport models. As nanoparticles grow from subnanometer to tens of nm in diameter, the numerical results predict a rich spatiotemporal structure, including four distinct temporal phases: a charge-limited phase, a charge accumulation phase, an early ion drag phase, and a sheath interaction phase.