迎风紧致格式求解Hamilton-Jacobi方程

基于Hamilton-Jaeobi(H-J)方程和双曲型守恒律之间的关系,将三阶和五阶迎风紧致格式推广应用于求解H-J方程,建立了高精度的H-J方程求解方法.给出了一维和二维典型数值算例的计算结果,其中包括一个平面激波作用下的Richtmyer-Meshkov界面不稳定性问题.数值试验表明,在解的光滑区域该方法具有高精度,而在导数不连续的不光滑区域也获得了比较好的分辨效果.相比于同阶精度的WENO格式,本方法具有更小的数值耗散,从而有利于多尺度复杂流动的模拟中H-J方程的求解.

NS方程激波计算的摄动有限差分方法

摄动有限差分（PFD）方法从一阶迎风差分格式出发，将差分系数展开为网格步长的幂级数，通过提高修正微分方程的逼近精度来获得更高精度的差分格式。由于格式基于一阶迎风格式，因此具有迎风效应、网格节点少等特点。本文首先通过对Burgers方程的摄动差分格式的推导，将摄动有限差分格式引入时间相关法的计算，并构造了守恒形式的摄动有限差分格式，然后推广到一维Navier-Stokes方程组的计算。数值比较研究表明：本文构造的NS方程摄动有限差分格式具有比一阶迎风较高的精度和分辨率，而且保持了一阶迎风格式的无振荡性质。

Direct Numerical Simulation Of Hypersonic Boundary-Layer Transition Over A Blunt Cone

Direct numerical simulation of transition How over a blunt cone with a freestream Mach number of 6, Reynolds number of 10,000 based on the nose radius, and a 1-deg angle of attack is performed by using a seventh-order weighted essentially nonoscillatory scheme for the convection terms of the Navier-Stokes equations, together with an eighth-order central finite difference scheme for the viscous terms. The wall blow-and-suction perturbations, including random perturbation and multifrequency perturbation, are used to trigger the transition. The maximum amplitude of the wall-normal velocity disturbance is set to 1% of the freestream velocity. The obtained transition locations on the cone surface agree well with each other far both cases. Transition onset is located at about 500 times the nose radius in the leeward section and 750 times the nose radius in the windward section. The frequency spectrum of velocity and pressure fluctuations at different streamwise locations are analyzed and compared with the linear stability theory. The second-mode disturbance wave is deemed to be the dominating disturbance because the growth rate of the second mode is much higher than the first mode. The reason why transition in the leeward section occurs earlier than that in the windward section is analyzed. It is not because of higher local growth rate of disturbance waves in the leeward section, but because the growth start location of the dominating second-mode wave in the leeward section is much earlier than that in the windward section.

Receptivity To Free-Stream Disturbance Waves For Hypersonic Flow Over A Blunt Cone

A high-order shock-fitting finite difference scheme is studied and used to do direction numerical simulation (DNS) of hypersonic unsteady flow over a blunt cone with fast acoustic waves in the free stream, and the receptivity problem in the blunt cone hypersonic boundary layers is studied. The results show that the acoustic waves are the strongest disturbance in the blunt cone hypersonic boundary layers. The wave modes of disturbance in the blunt cone boundary layers are first, second, and third modes which are generated and propagated downstream along the wall. The results also show that as the frequency decreases, the amplitudes of wave modes of disturbance increase, but there is a critical value. When frequency is over the critial value, the amplitudes decrease. Because of the discontinuity of curvature along the blunt cone body, the maximum amplitudes as a function of frequencies are not monotone.

Dns Of Compressible Turbulent Boundary Layer Around A Sharp Cone

Direct numerical simulation of the turbulent boundary layer over a sharp cone with 20 degrees cone angle (or 10 degrees half-cone angle) is performed by using the mixed seventh-order up-wind biased finite difference scheme and sixth-order central difference scheme. The free stream Mach number is 0.7 and free stream unit Reynolds number is 250000/inch. The characteristics of transition and turbulence of the sharp cone boundary layer are compared with those of the flat plate boundary layer. Statistics of fully developed turbulent flow agree well with the experimental and theoretical data for the turbulent flat-plate boundary layer flow. The near wall streak-like structure is shown and the average space between streaks (normalized by the local wall unit) keeps approximately invariable at different streamwise locations. The turbulent energy equation in the cylindrical coordinate is given and turbulent energy budget is studied. The computed results show that the effect of circumferential curvature on turbulence characteristics is not obvious.

Diffusion analogy and shock capturing for solving aerodynamic equations

Improving the resolution of the shock is one of the most important subjects in computational aerodynamics. In this paper the behaviour of the solutions near the shock is discussed and the reason of the oscillation production is investigated heuristically. According to the differential approximation of the difference scheme the so-called diffusion analogy equation and the diffusion analogy coefficient are defined. Four methods for improving the resolution of the shock are presented using the concept of diffusion analogy.

Characteristics of a CW HF chemical-laser calculated from a simplified two-dimensional model

A two-dimensional simplified model of an HF chemical laser is introduced. Using an implicit finite difference scheme, the solution of two adjacent parallel streams with diffusion mixing and chemical reaction is generated. A contour of mixing and reaction boundary is obtained without presupposition. The distribution of the HF(v) concentrations, gas temperature and the optical small signal gain (alpha sub V, J) on the flowing plane (X, Y) are presented. Compared with the solution solved directly from a set of Navier-Stokes equations, the results of these two methods agree with each other qualitatively. The influences of the different velocity, temperature (T sub 0) and composition of the two streams on the small signal gain after the nozzle exit are investigated. It is interesting that for larger J with a fixed v, the peaks of alpha sub v-T sub 0 profiles move towards higher T sub 0. The computing method is simple and only a short computing time is needed.

粘性可压混合层时间稳定性对称紧致差分求解

基于可压扰动方程组的一阶改型，将高精度对称紧致格式引入边值法数值线性稳定分析。对所获非线性离散特征值问题给出了一个通用形式二阶迭代局部算法，实现了时间模式和空间模式的统一求解，并将扰动特征及其特征函数同时得到。据此分析了可压平面自由混合层时间稳定性，涉及二维/三维扰动波、粘性/无粘扰动波、第一/第二模态、特征函数、伪特征值谱等。研究表明，压缩性效应和粘性效应对最不稳定扰动波数和增长率呈相似的减抑作用；在Mc = 1附近，从高波数段开始，粘性效应可强化二维不稳定扰动波由第一模态向第二模态的过渡。

若干流体力学问题的数学性质研究

该文利用高智的扩散抛物化方程组理论及流体力学基本方程组的特征次特征理论,流体大小尺度（LSS）方程组理论以及摄动有限差分（PFD）方法,研究若干流体力学问题的数学性质.该文得到的主要结论有：1.利用湍流大小尺度（LSS）方程组推导出湍流大小尺度涡量（LSSV）方程组,并证明两个关于湍流大小尺度涡量的命题,从而得到湍流封闭大小尺度涡量（CLSSV）方程组,并对已有的近程相互作用命题进行推广.2.根据扩散抛物化方程组理论和流体力学层次结构方程组的特征和次特征方法,研究了抛物化稳定性方程组（PSE）的特征和次特征以及消除PSE的剩余椭圆特性的问题.3.利用摄动有限差分（PFD）方法得到对流扩散反应方程的变步长摄动有限差分格式,是等步长摄动有限差分格式的推广.

A highly efficient beam propagation method for modeling step-index waveguides with tilt interface

Based on a new finite-difference scheme and Runge-Kutta method together with transparent boundary conditions (TBCs), a novel beam propagation method to model step-index waveguides with tilt interfaces is presented. The modified scheme provides an precies description of the tilt interface of the nonrectangular waveguide structure, showing a much better efficiency and accuracy comparing with the previously presented formulas.

A High-Efficiency Fiber-to-Waveguide Coupler with Low Polarization Dependence Using a Diluted Waveguide in InP Substrate with a 1.55μm Wavelength

We propose a fiber-to-waveguide coupler for side-illuminated p-i-n photodiodes to obtain high responsivity and low polarization dependence that is grown on InP substrate and is suitable for surface hybrid integration in low cost modules. The fiber-to-waveguide coupler is based on a diluted waveguide,which is composed of ten periods of undoped 120nm InP/80nm InGaAsP （1.05μm bandgap） multiple layers. Using the semi-vectorial three dimensional beam propagation method （BPM） with the central difference scheme,the coupling efficiency of fiber-to-waveguide under different conditions is simulated and studied,and the optimized conditions for fiber-to-waveguide coupling are obtained. For TE-like and TM-like modes,the calculated maximum coupling efficiency is higher than 94% and 92% ,respectively. The calculated polarization dependence is less than 0. ldB,showing good polarization independence.

Effects of Magnetohydrodynamic Interaction-Zone Position on Shock-Wave/Boundary-Layer Interaction

A third-order weighted essentially nonoscillatory and non-free-parameter difference scheme magnetohydrodynamic solver has been established to investigate the mechanisms of magnetohydrodynamics controlling separation induced by an oblique shock wave impinging on a flat plate. The effects of magnetohydrodynamic interaction-zone location on the separation point, reattachment point, separation-bubble size, and boundary-layer velocity profiles are analyzed. The results show that there exists a best location for the magnetohydrodynamic zone to be applied, where the separation point is delayed the farthest, and the separation bubble is decreased up to about 50% in size compared to the case without magnetohydrodynamic control, which demonstrated the promising of magnetohydrodynamics suppressing the separation induced by shock-wave/boundary-layer interactions.

地震波模拟：非规则网格离散的数值解法、吸收边界及其实际应用

This dissertation presents a series of irregular-grid based numerical technique for modeling seismic wave propagation in heterogeneous media. The study involves the generation of the irregular numerical mesh corresponding to the irregular grid scheme, the discretized version of motion equations under the unstructured mesh, and irregular-grid absorbing boundary conditions. The resulting numerical technique has been used in generating the synthetic data sets on the realistic complex geologic models that can examine the migration schemes. The motion equation discretization and modeling are based on Grid Method. The key idea is to use the integral equilibrium principle to replace the operator at each grid in Finite Difference scheme and variational formulation in Finite Element Method. The irregular grids of complex geologic model is generated by the Paving Method, which allow varying grid spacing according to meshing constraints. The grids have great quality at domain boundaries and contain equal quantities of nodes at interfaces, which avoids the interpolation of parameters and variables. The irregular grid absorbing boundary conditions is developed by extending the Perfectly Matched Layer method to the rotated local coordinates. The splitted PML equations of the first-order system is derived by using integral equilibrium principle. The proposed scheme can build PML boundary of arbitrary geometry in the computational domain, avoiding the special treatment at corners in a standard PML method and saving considerable memory and computation cost. The numerical implementation demonstrates the desired qualities of irregular grid based modeling technique. In particular, (1) smaller memory requirements and computational time are needed by changing the grid spacing according to local velocity; (2) Arbitrary surfaces and interface topographies are described accurately, thus removing the artificial reflection resulting from the stair approximation of the curved or dipping interfaces; (3) computational domain is significantly reduced by flexibly building the curved artificial boundaries using the irregular-grid absorbing boundary conditions. The proposed irregular grid approach is apply to reverse time migration as the extrapolation algorithm. It can discretize the smoothed velocity model by irregular grid of variable scale, which contributes to reduce the computation cost. The topography. It can also handle data set of arbitrary topography and no field correction is needed.

离散裂缝诱导各向异性粘弹性地震波模拟与特征分析

The real media always attenuate and distort seismic waves as they propagate in the earth. This behavior can be modeled with a viscoelastic and anisotropic wave equation. The real media can be described as fractured media. In this thesis, we present a high-order staggered grid finite-difference scheme for 2-D viscoelastic wave propagation in a medium containing a large number of small finite length fractures. We use the effective medium approach to compute the anisotropic parameters in each grid cell. By comparing our synthetic seismogram by staggered-grid finite-difference with that by complex-ray parameter ray tracing method, we conclude that the high-order staggered-grid finite-difference technique can effectively used to simulate seismic propagation in viscoelastic-anisotropic media. Synthetic seismograms demonstrate that strong attenuation and significant frequency dispersion due to viscosity are important factors of reducing amplitude and delaying arrival time varying with incidence angle or offset. On the other hand, the amount of scattered energy not only provides an indicator of orientation of fracture sets, but can also provide information about the fracture spacing. Analysis of synthetic seismograms from dry- and fluid-filled fractures indicates that dry-filled fractures show more significant scattering on seismic wavefields than fluid-filled ones, and offset-variations in P-wave amplitude are observable. We also analyze seismic response of an anticlinal trap model that includes a gas-filled fractured reservoir with high attenuation, which attenuates and distorts the so-called bright spot.

Fully developed laminar flow and heat transfer in smooth trapezoidal microchannel

Numerical analysis of fully developed laminar slip flow and heat transfer in trapezoidal micro-channels has been studied with uniform wall heat flux boundary conditions. Through coordinate transformation, the governing equations are transformed from physical plane to computational domain, and the resulting equations are solved by a finite-difference scheme. The influences of velocity slip and temperature jump on friction coefficient and Nusselt number are investigated in detail. The calculation also shows that the aspect ratio and base angle have significant effect on flow and heat transfer in trapezoidal micro-channel. (c) 2005 Elsevier Ltd. All rights reserved.