单晶硅化学气相沉积反应器流场初探

对单晶硅化学气相沉积(CVD)反应器在沉积过程中的流场进行了初步分析。通过数值求解三维层流Navier-Stokes方程，研究了反应器内浮力效应所引起的流场对称性破坏。结果表明，由于存在浮力效应，轴对称几何体中也会发生非轴对称流场分布，从而影响单晶硅的均匀生长。

LOCAL BENDING OF THIN FILM ON VISCOUS LAYER

Effects of deposition layer position and number/density on local bending of a thin film are systematically investigated. Because the deposition layer interacts with the thin film at the interface and there is an offset between the thin film neutral surface and the interface, the deposition layer generates not only axial stress but also bending moment. The bending moment induces an instant out-of-plane deflection of the thin film, which may or may not cause the so-called local bending. The deposition layer is modeled as a local stressor, whose location and density are demonstrated to be vital to the occurrence of local bending. The thin film rests on a viscous layer, which is governed by the Navier-Stokes equation and behaves like an elastic foundation to exert transverse forces on the thin film. The unknown feature of the axial constraint force makes the governing equation highly nonlinear even for the small deflection case. The constraint force and film transverse deflection are solved iteratively through the governing equation and the displacement constraint equation of immovable edges. This research shows that in some special cases, the deposition density increase does not necessarily reduce the local bending. By comparing the thin film deflections of different deposition numbers and positions, we also present the guideline of strengthening or suppressing the local bending.

零攻角NADTPS重构高超声速流场气动力/热研究

长时间、远距离高超声速飞行要求飞行器设计兼顾减阻和防热特点,论文提出无烧蚀自适应减阻防热新概念(NADTPS)并展开相关研究。NADTPS基于支杆-钝头体结构重构流场,利用侧向射流重点防护支杆头部。本文通过数值求解二维轴对称Navier-Stokes方程,对零攻角下NADTPS重构流场进行模拟,研究流场基本特征、气动力/热特性,分析新概念的减阻和防热机理。

超声速平板绕流的全流域阻力特性研究

超声速流动中摩擦阻力的准确预测一直是气动研究中的难点。本文以零攻角的平板绕流为对象，结合传统的CFD方法求解Navier-Stokes方程和直接模拟Monte-Carlo（DSMC）方法模拟粒子的微观运动，研究平板在超声速来流下的全流域阻力特性。研究发现，在超声速流中会产生激波-边界层相互干扰作用并随着来流马赫数的增加而加强，同时平板前缘存在非平衡效应显著的区域；努森数的增大使激波和边界层均增厚，粘性作用增强，非平衡区的范围增大，非平衡强度增强，而激波强度减弱；最后通过拟合连续流条件下的摩阻系数和修正自由分子流结果，结合桥函数，获得了全流域范围内均有效的平板阻力系数表达式。

支杆-钝头体带攻角飞行流场和“军刺”挡板作用机理研究

本文通过三维Navier-Stokes方程的数值求解，研究了支杆-钝头体结构在10°攻角M∞=6.0飞行条件下的流场结构和特点，指出其气动力特性恶化的原因，提出采用“军刺”挡板改善流场和气动力特性，并通过对比两种不同挡板作用下的流场和气动力特性变化分析其作用机理，发现“军刺”挡板结构分割流场抑制三维效应形成的周向流动，迎风面形成稳定的回流区和剪切层结构，将迎风面锥激波推离轴线，降低钝头体肩部流动结构相互作用强度，并在一定程度上缓解背风面流动干扰，明显改善支杆-钝头体带攻角飞行时的气动力特性。

用双流体模型模拟超声速气固两相流动

采用了双流体模型对JPL(Jet Propulsion Laboratory)喷管中气固两相流动以及超声速射流两相流动进行了数值模拟,并研究了可压缩两相流动中气相与颗粒的相互作用规律.自主开发的一般曲线坐标系下二维轴对称可压缩双流体程序Solve2D,对气相求解Navier-Stokes(N-S)方程组,采用k-ε湍流模型,颗粒相求解Euler方程组.对JPL喷管内的两相流场和湍流两相射流流场进行了数值模拟,研究了不同颗粒质量百分数以及不同颗粒直径时的气固两相流场的流动规律.

气泡/液滴运动的Level Set方法模拟研究

利用Level Set方法，结合投影法求解了描述气泡/液滴运动的Navier-Stokes方程。对地面常重力场中不同大小的空气泡在高黏度糖浆溶液中的自由上升运动现象，数值模拟结果与实验观测结果符合甚好，表明该方法能够计算大密度比和黏度比$(>1000:1)$情况下的气液两相流动。而对等密度液滴的热毛细迁移现象的数值模拟结果同样能够与实验结果相一致，表明该方法同样适于研究具有Marangoni效应的两相流动现象，特别是在空间微重力环境中的气液两相传热现象中的局部流动与传热问题。

三维波流数值水槽的构建方法

为了模拟波浪与水流的联合作用,基于Reynolds平均的Navier-Stokes方法,构建了一个三维波流数值水槽. 使用该模型对不同出流边界条件进行了对比分析. 结果表明:引入合理的入流和出流边界条件,可以保证在计算稳定后,水槽平均水面基本保持不变. 垂向时均流速分布的计算结果与实验数据吻合良好. 这表明:该波流数值水槽具有较好的波流特性,可用于计算波浪与水流、以及与三维固定结构之间相互作用.

可压缩湍流直接数值模拟

本书是以作者及其研究小组多年来的研究工作为基础撰写的，主要是可压缩湍流的直接数值模拟。全书共分两部分。第一部分为数值模拟方法：针对可压缩湍流的多尺度非定常非稳定流动特征，所构造的高精度差分格式（高于三阶精度），分析了各种高精度差分格式的数值行为特征和模拟小尺度物理结构、分辨非定常间断面的能力，重点分析了正确模拟多尺度物理结构及其非线性干扰对数值方法的要求。第二部分是物理问题的直接数值模拟：根据可压缩湍流的流动物理特征，所建立的求解可压缩Navier-Stokes方程的高精度数值模拟方法，详细描述了不稳定流和可压缩湍流的直接数值模拟，以及在此基础上所发现的一些新的物理现象。直接数值模拟主要针对典型可压缩湍流，包括均匀各向同性湍流、槽道湍流、平面混合流、射流、平板边界层湍流和钝楔、钝锥边界层湍流等，重点是可压缩湍流从层流转捩到湍流的直接数值模拟及其流动机理的研究。书中还介绍了我课组所研制的高精度、高效并行算法的软件。

Effects of buoyancy on Langmuir circulation

Based on the Navier-Stokes equation, an equation describing the Langmuir circulation is derived by a perturbation method when the influences of Coriolis force and buoyancy force are both considered. The approach used in the analysis is similar to the works carried out by Craik and Leibovich [J. Fluid Mech. 73 (1976) 401], Leibovich [J. Fluid Mech. 79 (1977) 715] and Huang [J. Fluid Mech. 91 (1979) 191]. Potential applications of the equation proposed are discussed in the area of Antarctic circumpolar current.

浮力对Langmuir环流的影响

Langmuir环流是发生在上层海洋中纵向的螺旋状涡旋运动，其轴向与风向基本一致。它是一常见的重要海洋现象，与海洋研究的许多重要问题，如海－气交换、海洋混合、海洋环境、海洋遥感、海洋生态和海洋灾害等有着十分密切的关系。 本文在同时考虑科氏力与浮力作用下，从Navier-Stokes方程出发，应用摄动理论，导出一个包含科氏力与浮力影响的Langmuir环流控制方程组。所采用的方法类似于Craik and Leibovich（1976）、Leibovich（1977）及Huang（1979）等人研究Langmuir环流时所采用的研究方法。此外，本文还讨论了所导出的方程组在南极绕极流区域的潜在应用。

波动方程模拟与偏移的实用化方法研究

Along with the widespread and in-depth applications in petroleum prospecting and development, the seismic modeling and migration technologies are proposed with a higher requirement by oil industrial, and the related practical demand is getting more and more urgent. Based on theories of modeling and migration methods for wave equation, both related with velocity model, I thoroughly research and develop some methods for the goal of highly effective and practical in this dissertation. In the first part, this dissertation probes into the layout designing by wave equations modeling, focusing on the target-oriented layout designing method guided by wave equation modeling in complicated structure areas. It is implemented by using the fourth order staggered grid finite difference (FD) method in velocity-stress 2D acoustic wave equations plus perfectly matched layer (PML) absorbing boundary condition. To design target-oriented layout: (a) match the synthetic record on the surface with events of subsurface structures by analyzing the snapshots of theoretical model; (b) determine the shot-gather distance by tracking the events of target areas and measuring the receiving range when it reaches the surface; (c) restrict the range of valid shot-gather distance by drawing seismic windows in single shot records; (d) choose the best trace distance by comparing the resolution of prospecting targets from the simulated records with different trace distance. Eventually, we obtained the observation system parameters, which achieve the design requirements. In the second part, this dissertation presents the practical method to improve the 3D Fourier Finite Difference (FFD) migration, and carefully analyzes all the factors which influence 3D FFD migration’s efficiency. In which, one of the most important parameters of migration is the extrapolating step. This dissertation presents an efficient 3D FFD migration algorithm, which use FFD propagator to extrapolate wavefields over big layers, and use Born-Kirchhoff interpolator to image wavefields over small layers between the big ones. Finally, I show the effectiveness of this hybrid migration method by comparing migration results from 3D SEG/EAGE model with different methods.

The effects of random surface waves on the steady Ekman current solutions

The response of near-surface current profiles to wind and random surface waves are studied based on the approach of Jenkins [1989. The use of a wave prediction model for driving a near surface current model. Dtsch. Hydrogr. Z. 42,134-149] and Tang et al. [2007. Observation and modeling of surface currents on the Grand Banks: a study of the wave effects on surface currents. J. Geophys. Res. 112, C10025, doi:10.1029/2006JC004028]. Analytic steady solutions are presented for wave-modified Ekman equations resulting from Stokes drift, wind input and wave dissipation for a depth-independent constant eddy viscosity coefficient and one that varies linearly with depth. The parameters involved in the solutions can be determined by the two-dimensional wavenumber spectrum of ocean waves, wind speed, the Coriolis parameter and the densities of air and water, and the solutions reduce to those of Lewis and Belcher [2004. Time-dependent, coupled, Ekman boundary layer solutions incorporating Stokes drift. Dyn. Atmos. Oceans. 37, 313-351] when only the effects of Stokes drift are included. As illustrative examples, for a fully developed wind-generated sea with different wind speeds, wave-modified current profiles are calculated and compared with the classical Ekman theory and Lewis and Belcher's [2004. Time-dependent, coupled, Ekman boundary layer solutions incorporating Stokes drift. Dyn. Atmos. Oceans 37, 313-351] modification by using the Donelan and Pierson [1987. Radar scattering and equilibrium ranges in wind-generated waves with application to scatterometry. J. Geophys. Res. 92, 4971-5029] wavenumber spectrum, the WAM wave model formulation for wind input energy to waves, and wave energy dissipation converted to currents. Illustrative examples for a fully developed sea and the comparisons between observations and the theoretical predictions demonstrate that the effects of the random surface waves on the classical Ekman current are important, as they change qualitatively the nature of the Ekman layer. But the effects of the wind input and wave dissipation on surface current are small, relative to the impact of the Stokes drift. (C) 2008 Elsevier Ltd. All rights reserved.

A Cip/Multi-Moment Finite Volume Method For Shallow Water Equations With Source Terms

A novel finite volume method has been presented to solve the shallow water equations. In addition to the volume-integrated average (VIA) for each mesh cell, the surface-integrated average (SIA) is also treated as the model variable and is independently predicted. The numerical reconstruction is conducted based on both the VIA and the SIA. Different approaches are used to update VIA and SIA separately. The SIA is updated by a semi-Lagrangian scheme in terms of the Riemann invariants of the shallow water equations, while the VIA is computed by a flux-based finite volume formulation and is thus exactly conserved. Numerical oscillation can be effectively avoided through the use of a non-oscillatory interpolation function. The numerical formulations for both SIA and VIA moments maintain exactly the balance between the fluxes and the source terms. 1D and 2D numerical formulations are validated with numerical experiments. Copyright (c) 2007 John Wiley & Sons, Ltd.

3D Finite Volume Scheme for Czochralski Single Crystal Growth

Czochralski (Cz) technique, which is used for growing single crystals, has dominated the production of single crystals for electronic applications. The Cz growth process involves multiple phases, moving interface and three-dimensional behavior. Much has been done to study these phenomena by means of numerical methods as well as experimental observations. A three-dimensional curvilinear finite volume based algorithm has been developed to model the Cz process. A body-fitted transformation based approach is adopted in conjunction with a multizone adaptive grid generation (MAGG) technique to accurately handle the three-dimensional problems of phase-change in irregular geometries with free and moving surfaces. The multizone adaptive model is used to perform a three-dimensional simulation of the Cz growth of silicon single crystals.Since the phase change interface are irregular in shape and they move in response to the solution, accurate treatment of these interfaces is important from numerical accuracy point of view. The multizone adaptive grid generation (MAGG) is the appropriate scheme for this purpose. Another challenge encountered is the moving and periodic boundary conditions, which is essential to the numerical solution of the governing equations. Special treatments are implemented to impose the periodic boundary condition in a particular direction and to determine the internal boundary position and shape varying with the combination of ambient physicochemical transport process and interfacial dynamics. As indicated above that the applications and processes characterized by multi-phase, moving interfaces and irregular shape render the associated physical phenomena three-dimensional and unsteady. Therefore a generalized 3D model rather than a 2D simulation, in which the governing equations are solved in a general non-orthogonal coordinate system, is constructed to describe and capture the features of the growth process. All this has been implemented and validated by using it to model the low pressure Cz growth of silicon. Accuracy of this scheme is demonstrated by agreement of simulation data with available experimental data. Using the quasi-steady state approximation, it is shown that the flow and temperature fields in the melt under certain operating conditions become asymmetric and unsteady even in the absence of extrinsic sources of asymmetry. Asymmetry in the flow and temperature fields, caused by high shear initiated phenomena, affects the interface shape in the azimuthal direction thus results in the thermal stress distribution in the vicinity, which has serious implications from crystal quality point of view.