短脉冲激光与介质相互作用的扩散模型

推导了扩散近似方程,通过半无限大均匀介质计算,用扩散理论分析解验证了数值方法的有效性.模拟了光在非均匀介质内的传输过程,给出了介质内光通量随时间变化的空间分布.结果表明,该基于扩散模型的数值方法能够模拟短脉冲光在强散射介质中的传播过程以及漫散射光的时间变化特性,并且借助于光通量空间分布能够准确模拟非均匀介质内内含物的位置.     

Wavelength selection in Rayleigh-Bénard convection

Rayleigh-Bénard convection is considered in the presence of a weakly nonuniform temperature distribution and weakly nonuniform layer height with their gradients perpendicular to the (parallel) convection rolls. In the infinite Prandtl number limit the phase equation that describes the slow spatial and temporal evolution of the local wave number is derived. Evaluating the equation near threshold and for stress-free boundary conditions we find nonuniversal wavenumber selection, forced phase diffusion with moving patterns, and other measurable effects.     

Studies on aluminum powder combustion in detonation environment

The combustion mechanism of aluminum particles in a detonation environment characterized by high temperature (in unit 10³ K), high pressure (in unit GPa), and high-speed motion (in units km/s) was studied, and a combustion model of the aluminum particles in detonation environment was established. Based on this model, a combustion control equation for aluminum particles in detonation environment was obtained. It can be seen from the control equation that the burning time of aluminum particle is mainly affected by the particle size, system temperature, and diffusion coefficient. The calculation result shows that a higher system temperature, larger diffusion coefficient, and smaller particle size lead to a faster burn rate and shorter burning time for aluminum particles. After considering the particle size distribution characteristics of aluminum powder, the application of the combustion control equation was extended from single aluminum particles to nonuniform aluminum powder, and the calculated time corresponding to the peak burn rate of aluminum powder was in good agreement with the experimental electrical conductivity results. This equation can quantitatively describe the combustion behavior of aluminum powder in different detonation environments and provides technical means for quantitative calculation of the aluminum powder combustion process in detonation environment.     

Diffusion near nonequilibrium steady state

The diffusion equation describing the response of a system in a steady nonequilibrium state to a nonuniform external perturbation is derived from the Boltzmann equation. For transient processes, the diffusion equation governs the slow (final) stage of the process, irrespective to the form of the initial distribution. The initial velocity distribution becomes the local steady-state distribution during the first fast stage of the process. This fast initial relaxation affects only the initial condition for the diffusion equation. Some principal features of the nonequilibrium diffusion constant are discussed. As an illustration, the diffusion constant of hot electrons is calculated in the electron temperature approximation.     

Diffusion in inhomogeneous media

The problem is to establish the correct diffusion equation in a medium that is inhomogeneous and whose temperature also varies in space. As a special model we study particles whose phase space distribution obeys Kramers' equation with a generalized collision operator. In the usual limit of strong collisions a diffusion equation is obtained. This equation contains additional drift terms, which depend on the form of the collision operator. They cannot be expressed as a mobility and a diffusion coefficient, unless the decay law of the velocity happens to be linear. Conclusion: no universal form of the diffusion equation exists, but each system has to be studied individually.Dedicated to Professor Harry Thomas on the occasion of his 60th birthday     

Diffusion in flows with a linearly nonuniform velocity profile

Methods of group analysis are applied to determine the Green’s function in the infinite space of the test particle diffusion equation for flows with a linearly nonuniform velocity profile. Using this function, the solutions of the Cauchy problem are obtained for certain particular initial conditions.     

Large energy behavior of the velocity distribution for the hard-sphere gas

The initial-value problem for the Boltzmann-Lorentz equation for hard spheres at zero temperature is shown to be ill defined, the general solution depending on an arbitrary function. The uniqueness of the solution can be obtained by imposing the conservation of the number of particles (Carleman's type of condition does not suffice). The linearized Boltzmann equation for hard spheres is then analyzed, as it occurs in Enskog's method for calculating transport coefficients. It is demonstrated that in the case of viscosity and diffusion it is necessary to add supplementary conditions to obtain the uniqueness of the solution. The nonuniform character of Enskog's expansion and violation of positivity in the large velocity region are exhibited.     

多维梯度折射率介质内辐射传递的扩散近似

推导多维梯度折射率介质内稳态辐射传递的扩散近似方程.使用有限元法对扩散近似进行离散和求解,利用两个二维半透明介质的稳态辐射传递问题验证该扩散近似的精度及适用性.算例考虑介质为均匀折射率及梯度折射率两种情况.利用扩散近似分别求解辐射平衡时的边界热流、介质内温度场分布,并与辐射传递方程的求解结果进行对比分析.结果表明:介质折射率变化、散射特性、光学厚度及散射反照率均直接影响扩散近似的精度;在光学厚及强散射条件下,该扩散近似可以作为一种快速算法应用于梯度折射率介质稳态辐射传递的求解.     

Multigrid method based on the transformation-free HOC scheme on nonuniform grids for 2D convection diffusion problems

In this paper, a multigrid method based on the high order compact (HOC) difference scheme on nonuniform grids, which has been proposed by Kalita et al. [J.C. Kalita, A.K. Dass, D.C. Dalal, A transformation-free HOC scheme for steady convection–diffusion on non-uniform grids, Int. J. Numer. Methods Fluids 44 (2004) 33–53], is proposed to solve the two-dimensional (2D) convection diffusion equation. The HOC scheme is not involved in any grid transformation to map the nonuniform grids to uniform grids, consequently, the multigrid method is brand-new for solving the discrete system arising from the difference equation on nonuniform grids. The corresponding multigrid projection and interpolation operators are constructed by the area ratio. Some boundary layer and local singularity problems are used to demonstrate the superiority of the present method. Numerical results show that the multigrid method with the HOC scheme on nonuniform grids almost gets as equally efficient convergence rate as on uniform grids and the computed solution on nonuniform grids retains fourth order accuracy while on uniform grids just gets very poor solution for very steep boundary layer or high local singularity problems. The present method is also applied to solve the 2D incompressible Navier–Stokes equations using the stream function–vorticity formulation and the numerical solutions of the lid-driven cavity flow problem are obtained and compared with solutions available in the literature.     

Effective relaxation time of the impurity concentration in a medium with a slightly nonuniform diffusion coefficient

Based on the previously introduced integral criterion and a new approach to estimating the time characteristics of diffusion, the relaxation time of impurity concentration in a medium with a slightly nonuniform diffusion coefficient is determined. The diffusion coefficient profile that cuts the relaxation time is found.     

Self-Consistent Model and Numerical Analysis of VCSEL''s Laser Diodes

A self-consistent model of vertical-cavity surface-emitting lasers (VCSEL's) is presented in this paper, in which the carrier diffusion rate equation, the photon density rate equation, and the thermal conduction equation are considered simultaneously. The nonuniform heat flux density distribution in the active region due to the current-spreading effect is taken into account. The effects of temperature on gain and transparency carrier density are also included in this model. It is nonlinearly and self-consistently solved in the Matlab environment. The transient and lateral distribution characteristics of the carrier and photon densities and the junction temperature are investigated. Some interesting results are gotten and analyzed.     

Correlation transfer and diffusion of ultrasound-modulated multiply scattered light

We develop a temporal correlation transfer equation (CTE) and a temporal correlation diffusion equation (CDE) for ultrasound-modulated multiply scattered light. These equations can be applied to an optically scattering medium with embedded optically scattering and absorbing objects to calculate the power spectrum of light modulated by a nonuniform ultrasound field. We present an analytical solution based on the CDE and Monte Carlo simulation results for light modulated by a cylinder of ultrasound in an optically scattering slab. We further validate with experimental measurements the numerical calculations for an actual ultrasound field. The CTE and CDE are valid for moderate ultrasound pressures and on a length scale comparable with the optical transport mean-free path. These equations should be applicable to a wide spectrum of conditions for ultrasound-modulated optical tomography of soft biological tissues.     

Inverse Problems of the Near-Field Radiothermometry

It is shown that the one-dimensional problem of near-field microwave radiometry in the case of a smoothly inhomogeneous medium can be reduced to an integral equation similar to the well-known solution of the radiation transfer equation. The kernel of this integral equation depends on only one parameter determined experimentally. This makes it possible to determine the subsurface temperature profiles by using a set of near-field antennas of arbitrary design for which the sounding-depth range is overlapped by the effective depth of formation of the radiation received by these antennas. In the case of near-field measurements, the previously obtained joint solution of the radiation transfer equation and the thermal diffusion equation can also be modified and applied for monitoring the subsurface temperature profile using the temporal dependence of the received signal.     

双次曝光积分效应实现杂质浓度分布均匀化

 激光诱导扩散中,当入射激光光强为高斯分布甚至均匀分布时，微小扩散区的温度分布不均匀。由于扩散系数是温度的函数，必将导致扩散后杂质浓度分布的均匀性较差，无法制作出高性能的p-n结。提出采用多次激光诱导扩散的积分效应来实现杂质浓度分布的均匀化整形。对于InP衬底的CO_{2激光诱导Zn扩散,利用温度闭环测控系统测得的基片表面热斑温度场分布，分析计算了两次激光诱导扩散重叠区域的浓度分布积分效应。在此基础上模拟计算出，用双次曝光积分效应做杂质浓度分布的均匀化整形时，基片上两次激光照射位置的最佳间隔为20 μm。这为改进激光诱导扩散工艺，用多次曝光实现面均匀的杂质浓度分布奠定了理论基础。}     

奇异退化扩散反应方程的高阶紧致差分及网格自适应方法

利用余项修正法建立奇异退化扩散反应方程非均匀网格上的高阶紧致差分式,其时间具有二阶精度,空间具有三阶至四阶精度. 利用等分布原理建立时间和空间的网格自适应方法.最后通过具有精确解的数值算例验证方法的可靠性和精确性,并研究一维爆破问题.     

声致发光气泡内水蒸气的影响

提出了一个单气泡声致发光的简单计算模型.这个模型是在均匀压强近似下，考虑质量和温度在气泡内的非均匀分布，同时考虑了水蒸气在气泡壁上的凝结与蒸发以及水蒸气在气泡内相对惰性气体的质量扩散.通过Saha方程估算气体电离密度，利用电子与离子、电子与中性粒子的轫致辐射，电子与离子的复合辐射公式估算气泡的辐射强度.不考虑化学反应，计算了不同水温时的气泡发光强度，发现当水温在0 ℃时轫致辐射发光模型比较符合实验结果，水温升高时，如水温为20 ℃或以上，轫致辐射发光模型的计算与实验结果出现数量级差别.考虑化学反应，轫致辐射发光模型的计算则总是比实验结果低2个数量级.     

Conductivity of structurally anisotropic 2D composites

The character of electric conduction in a two-dimensional structurally anisotropic composite containing elliptic inclusions with a large ratio of the semiaxis. Dielectric and ideally conducting inclusions were considered. It was demonstrated that the approximation linear in concentration of inclusions fails beginning from rather small concentrations. The region of intermediate concentrations was studied using approximate methods based on the analogy between electric conductivity and the diffusion of a single particle in a nonuniform medium. The conductivity of the model system under study was analyzed over a wide concentration range with the use of the effective medium theory; the results are in close agreement with those obtained by qualitative methods. The behavior of conductivity near the percolation threshold (metal-to-dielectric phase transition) was examined using the similarity. The conductivity of thin composite films with inclusions in the form of carbon nanotubes was examined, as an example.     

Scattering of long-wavelength acoustic phonons by isotopic impurities

We consider the problem of the relaxation of an arbitrary initial distribution function of a gas of long-wave acoustic phonons scattered by isotopic impurities embedded in a crystalline medium with cubic symmetry. The spectral decomposition of the collision integral of the suitable Boltzmann-Peierls equation is obtained. The spectrum of the collision operator is purely discrete and in addition to the eigenvalue 0 consists of three other eigenvalues. Explicit analytic expressions for these eigenvalues are obtained. Within the Chapman-Enskog approximation we derive the diffusion equation for the density of phonons and obtain the explicit expression for the diffusion coefficient. The dependency of the eigenvalues of the collision operator and the diffusion coefficient on the elastic constants of the medium is studied.     

Use of the two-moment boundary condition in the problem of rarefied gas slip over a solid cylindrical surface

The slip velocity of a rarefied gas nonuniform in temperature and mass velocity is calculated for gas slip over the surface of a right circular cylinder. The calculation uses the two-moment boundary condition in an approximation linear in Knudsen number. Corrections to the slip velocity that are due to the interface curvature, volume temperature stresses, and nonuniform temperature distribution in the Knudsen layer are studied as func-tions of the accommodation coefficients in the first two moments of the distribution function. The Bhatnagar-Gross-Krook model of the Boltzmann kinetic equation is employed as the basic equation for the gas state.