浅水中污染物扩散的自适应有限元分析法

介绍浅水中污染物扩散分析中的有限元法.分析包括两个部分:1)流场速度、水面高度的计算;2)根据扩散模型计算污染物浓度场.联合使用了自适应网格技术以期提高解的精度,同时减少计算时间和计算机内存的消耗.通过几个有已知解的实例验证了有限元公式和计算机程序.最后,使用这种联合方法分析泰国Chao Phraya河附近海湾中的污染物扩散.     

A numerical experimental study of inverse preconditioning for the parallel iterative solution to 3D finite element flow equations

Integration of the subsurface flow equation by finite elements (FE) in space and finite differences (FD) in time requires the repeated solution to sparse symmetric positive definite systems of linear equations. Iterative techniques based on preconditioned conjugate gradients (PCG) are one of the most attractive tool to solve the problem on sequential computers. A present challenge is to make PCG attractive in a parallel computing environment as well. To this aim a key factor is the development of an efficient parallel preconditioner. FSAI (factorized sparse approximate inverse) and enlarged FSAI relying on the approximate inverse of the coefficient matrix appears to be a most promising parallel preconditioner. In the present paper PCG using FSAI, diagonal and pARMS (parallel algebraic recursive multilevel solvers) preconditioners is implemented on the IBM SP4/512 and CLX/768 supercomputers with up to 32 processors to solve underground flow problems of a large size. The results show that FSAI may allow for a parallel relative efficiency larger than 50% on the largest problems with p=32 processors. Moreover, FSAI turns out to be significantly less expensive and more robust than pARMS. Finally, it is shown that for p in the upper range may be much improved if PCG–FSAI is implemented on CLX.     

Modeling of radiative heat transfer in an axisymmetric cylindrical enclosure with participating medium

Radiative heat transfer in an axisymmetric enclosure with absorbing, emitting, and scattering medium is studied here by using the different methods such as MDOM, FVM, and MFVM with emphasis on the treatment of angular derivative term, which appears in curvilinear coordinates due to angular redistribution. After final discretization equation for MFVM is introduced by using the step scheme and directional weights, the present approach is validated by applying it to three different benchmarking problems with absorbing, emitting, and scattering medium. All of the results presented here support its accuracy as well as moderate efficiency. Finally, the present approaches are applied to a truncated cone-shaped enclosure as a body-fitted geometry case.     

Finite element modal solutions of planar photonic crystal fibers with rectangular air-holes

Modal solutions of planar photonic crystal waveguides with rectangular air-holes are presented by using a rigorous full-vectorial finite element-based approach. The effective indices, mode field profiles, spot-sizes, power confinements, modal hybridness, beat lengths and group velocity dispersions are shown for the fundamental and higher order modes of the quasi-TE and TM polarizations.     

A complete boundary integral formulation for compressible Navier–Stokes equations

A complete boundary integral formulation for compressible Navier–Stokes equations with time discretization by operator splitting is developed using the fundamental solutions of the Helmholtz operator equation with different order. The numerical results for wall pressure and wall skin friction of two‐dimensional compressible laminar viscous flow around airfoils are in good agreement with field numerical methods. Copyright © 2004 John Wiley & Sons, Ltd.     

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.     

空位在金刚石近(001)表面扩散的分子动力学模拟

用分子动力学方法模拟了空位在金刚石近(001)表面的扩散过程,研究了温度对空位扩散的影响.结果表明,当温度为1000K左右时,位于近表面第二层上的空位开始向表面运动;当温度在1400—2000K时,空位完全扩散到表面.这与实验结果和其他计算结果符合得很好.同时发现,温度为1400—1800K时,空位的扩散经历了两次迁移运动,其分别对应了均方位移图中的两个极大值.在不施加任何约束的条件下得到了空位的动态扩散路径,空位在金刚石近(001)表面的扩散势垒约为042eV.并探讨了一定温度下空位数目增多及其不同排列 关键词： 金刚石 空位 扩散 分子动力学     

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.