对称性及多群中子扩散方程数值解

在多群中子扩散方程解析解的基础上,利用方程及求解域的对称性建立了新的数值求解中子扩散方程的理论模型.该模型显著的优点是适用于各种对称区域(二维、三维区域)尤其是非正方形区域中子扩散方程的求解,它彻底避免了常规节块法应用于非正方形几何时所出现的奇异性问题,且所得的解在求解域内任意点上均满足扩散方程.以二、三维六角形几何为例建立了数学模型,并用基准问题校核了模型的正确性. 关键词： 中子扩散方程 对称群 数值解 解析     

局域浓度调控扩散系数的次氯酸-碘离子-丙二酸系统图灵斑图形成中的反常扩散

通过数值模拟及振幅方程解析解方法,从实空间和倒空间分析了受局域浓度扩散系数调控下次氯酸-碘离子-丙二酸反应扩散系统图灵斑图形成的扩散机理.在零扩散系数调节下,斑图形成为典型的菲克扩散;而在负向正向扩散系数调节下,斑图的形成依赖欠扩散和超扩散.图灵系统的浓度稳态振幅对随机初始条件敏感性随局域浓度扩散调控系数k的增大而增加.     

耐腐蚀聚合药物基体上药物的扩散释放

主要对耐腐蚀聚合药物基体上药物的整个释放过程进行解析和数值研究.数学模型为基于Fick第二法则的边值问题,分为状态Ⅰ和Ⅱ.状态Ⅰ为具有移动性扩散波头的Stefan问题,而当此扩散波头消失时,即开始状态Ⅱ.该研究工作可作为数值研究更加复杂系统的基础之一.     

一维非稳态多层介质导热问题的研究

本文针对一维非稳态多层介质的导热问题,建立了导热的数学模型,通过数值解代替解析解的方法,求得导热过程中的温度分布的情况。首先,建立各层介质热传导的偏微分方程,确定了该方程的定解条件,从而得到了一维非稳态多层介质导热的数学模型。其次,由于偏微分方程的解析解难以求解,所以利用差分法将偏微分方程的解析解用数值解代替,再运用Crank-Nicolson方法将差分后的方程转化为三对角线性方程组,利用追赶法进行求解。最后,将所建立的导热模型运用到实例中,验证了该模型的实用性和可行性。     

一类非线性微分方程的近似解析解

从理论上研究了一类广义扩散方程的求解问题. 利用相似变换和解析拆分技巧给出了求解该类非线性微分方程近似解的一种有效方法, 方程的解可以表示为一个收敛的幂级数. 近似解结果和数值结果非常符合，证明了所提出的方法的准确性和可靠性, 该方法可以用于解决其他科学和工程技术问题. 关键词： 广义扩散方程 非线性边界值问题 解析拆分 近似解析解     

三能级Upper-ladder型系统中的

三能级Upper-ladder型系统中,在旋波、慢变振幅近似下,求解了考虑驱动场相位扩散后的系统密度矩阵运动方程,并给出了这个三能级梯型系统稳态线性解析解.利用对密度矩阵运动方程的稳态线性解析解的数值模拟结果,研究相位扩散对无反转激光增益、色散和粒子数差的影响;利用对密度矩阵运动方程的数值模拟结果,分析相位的扩散对无反...     

分形油藏非牛顿幂律流体低速非达西不稳定渗流的组合数学模型

结合分形理论与渗流理论,对分形油藏非牛顿幂律流体低速非达西不稳定渗流的试井分析问题的数学模型进行了推导.该分形油藏模型由内域为非牛顿幂律流体低速非达西渗流,外域为非牛顿幂律流体达西渗流的同心圆域组成.在考虑井筒储存、表皮效应影响下,建立了该油藏的不稳定渗流有效井径组合数学模型,在3种外边界条件下求出了两个区域内压力在Laplace空间的解析解,应用Stehfest数值反演方法求得井底的无因次压力,分析了井底压力动态特征和参数影响.非牛顿幂律流体的幂律指数、分形参数均对典型曲线产生较大的影响,呈现出与牛顿流体和均质油藏明显不同的特征.这对非均质油藏非牛顿流体的不稳定试井分析及研究其非线性渗流特征均十分重要.     

波纹板片上降膜吸收过程的传热传质数值仿真

建立了在连续的竖直波纹板上降膜吸收过程的传热传质数学模型,并进行了无因次化处理和仿真数值计算。给出了溶液液膜内流场、温度场、质量分数分布的计算结果,比较了交叉双尺度波纹板片和当量平板上的吸收传热和传质系数随不同溶液喷淋密度的变化结果,显示出前者具有更好的性能,并讨论了波纹节距对传热和传质系数的影响关系。     

解二维时间相关光子扩散方程的交替方向隐式法

研究了应用交替方向隐式法求解非均匀生物组织中的二维时间相关光子扩散方程,并同解析解以及Monte Carlo模拟的结果作了比较,数值结果表明,该方法具有较好的稳定性和较高的精度.     

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

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

A Monte Carlo model for seeded atomic flows in the transition regime

A simple model for the numerical determination of separation effects in seeded atomic gas flows is presented. The model is based on the known possibility to provide a statistically convergent estimate of the exact solution for a linear transport equation using the test particle Monte Carlo method. Accordingly, the flow field of the main gas is preliminary calculated and as a second step the linear transport equations obtained by fixing the target distribution in the collision term of the Boltzmann equation for both main and minority components are solved. Both solutions are based on appropriately devised test particle Monte Carlo methods. The second step, the critical one in evaluating the separation effects, is exact and thereby completely free of numerical diffusion. The model is described in details and illustrated by 2D test cases of atomic separation in shock fronts.     

Thermal Diffusion Process Estimation for Fabrication of Graded Plastic Optical Fibre

Thermal diffusion of dopants is investigated in the process of generating the graded-index profile of plastic optical fibres. Because the diffusion coefficient in high polymers has been shown to depend strongly on dopant concentration, it is allowed in this work to vary with the radial coordinate of the multistep-core fibre. A novel multi-layer model is presented for solving the diffusion equation with the variable diffusion coefficient. It is solved by the finite difference method. The solution determines the dopant diffusion profile in the fibre. It is verified against a solution from the literature and two cases of fibres with diffused profiles. The application is demonstrated on two examples of graded-index plastic optical fibres, one originally with a two-step and the other with four-step core. The results indicate that closer to the core-cladding interface, the computed diffused profile with variable diffusion coefficient D is closer to target profile than the profile obtained with constant D for the same time of thermal process.     

Diffusion in multiphase systems with nonstationary boundaries I. Solution of diffusion in the region with boundaries movig in one direction

In this work the problem of diffusion in a multiphase system whose both boundaries move in one direction according to the parabolic law with different velocities is analysed. The mathematical problem is solved exactly by means of thermal potentials of a double layer. The solution of the diffusion equation in the proximity to the boundary was derived and the concentration gradients on these boundaries were calculated. The numerical procedure of determining the diffusion characteristics from experimental concentration gradients on the phase boundaries was presented. As the zero approximation the result of calculations according to Vasileff and Smoluchowski, that can lead to considerable differences in the determined diffusion coefficients, was used. theory numerical calculations     

Diffusion in multiphase systems with nonstationary boundaries

In this work the problem of diffusion in the definite diffusion field whose one boundary is fixed and the other — interphase boundary — moves in time according to the parabolic law. The mathematical problem is solved exactly by means of thermal potentials of a double layer. The solution of the diffusion equation in the proximity to the boundary was derived and the concentration gradients on these boundaries were calculated. The numerical procedure of determining the diffusion characteristics from experimental concentration gradients on the phase boundaries was presented. As the zero approximation the result of calculations according to Wagner’s solution was used. Theory. Numerical calculations.     

自驱动Janus微球近壁运动特性实验与数值模拟研究

自驱动Janus微球是形状规则但表面构成不同的特殊活性颗粒. 针对微米级Pt-SiO₂型Janus 微球近壁面自驱动现象, 实验测得了微球的自驱动速度V_Janus, 并观察到微球运动过程中与垂直方向存在一偏转仰角ψ, 且ψ角随H₂O₂溶液浓度的增大呈减小趋势. 在此基础上, 建立自驱动Janus微球的数值模型, 通过模拟得到了微球在不同浓度H₂O₂溶液中的偏转仰角ψ及距底面的高度δ, 模拟与实验一致. 利用这些数据进一步讨论了壁面效应对微球旋转特征时间τ_R的影响. 这一工作对于理解Janus 微球的运动机理及发展相关应用具有重要意义.     

基于二维广义横向积分综合一维半解析的三维中子扩散节块法

目前, 轻水堆堆芯计算广泛采用基于横向积分技术的中子扩散方程节块展开法, 该方法需要对横向泄漏进行多项式近似而不严格, 而且堆芯设计还需要额外的精细功率重构模块用于获得组件内各栅元的功率分布。本文提出两维广义横向积分方法, 直接采用源展开以及表面流耦合方法, 可以避免上述两个不足。由于反应堆堆芯不均匀性更多发生在径向, 因此采用径向基于广义横向积分方法结合轴向常规节块法的综合方法进行三维中子扩散计算。通过基准问题的数值计算, 验证了该方法对于堆芯扩散计算的正确性和可行性。     

基于格子Boltzmann方法的超声速非平衡流模拟

基于D1Q4可压缩格子Boltzmann模型,按照流通矢量分裂方法的思路,采用坐标旋转技术构造求解三维带化学反应Navier-Stokes方程对流通量求解器.结合有限体积法求解三维化学非平衡流Navier-Stokes方程,采用时间算子分裂算法解决化学反应刚性问题,数值模拟超声速化学非平衡流的三个经典算例.数值结果表明:在高马赫数下,采用D1Q4可压缩格子Boltzmann模型构造的三维对流通量求解器数值模拟中没有出现非物理解,同时在超声速化学非平衡流场中正确分辨激波、燃烧波等物理现象,精度和分辨率均较高,验证了本文构造的三维对流通量求解器的可靠性,拓宽了D1Q4可压缩格子Boltzmann模型的应用范围,为计算超声速化学非平衡流提供一种新方法.     

Analysis of first stage ignition delay times of dimethyl ether in a laminar flow reactor

The combustion chemistry of the first stage ignition and chemistry/flow interactions are studied for dimethyl ether (DME) with a mathematical analysis of two systems: a plug flow reactor study is used to reduce the reaction chemistry systematically. A skeletal reaction mechanism for the low temperature chemistry of DME until the onset of ignition is derived on the basis of the detailed DME mechanism of the Lawrence Livermore National Laboratory – see Curran, Fischer and Dryer, Int. J. Chem. Kinetics, Vol. 32 (2000). It is shown that reasonably good results for ignition delay times can be reached using a simple system of three ordinary differential equations and that the resulting analytical solution depends only on two reaction rates and the initial fuel concentration. The stepwise reduction of the system based on assumptions yields an understanding on why these reactions are so important. Furthermore, the validation of the assumptions yields insight into the influence of the fuel and the oxygen concentration on the temperature during the induction phase. To investigate the influence of chemistry/flow interactions, a 2D model with a laminar Hagen–Poiseuille flow and 2D-polynomial profiles for the radial species concentration is considered. For the 2D model, it is found that only the diffusion coefficients and the reactor radius need to be taken into consideration additionally to describe the system sufficiently. Also, the coupling of flow and chemistry is clarified in the mathematical analysis. The insight obtained from the comparison of the 2D model and the plug flow model is used to establish an average velocity for the conversion of ignition locations to ignition delay times in a laminar flow reactor. Finally, the 2D analytical solution is compared against new experimental data, obtained in such a laminar flow reactor for an undiluted DME/air mixture with an equivalence ratio of φ = 0.835 and a temperature range of 555 to 585 K at atmospheric pressure.     

两相湍流色噪声扩维法PDF模型及数值模拟

由颗粒运动的朗之万方程出发,对流体脉动速度采用扩维方法,得到两个不同层次的PDF输运方程.通过对颗粒运动方程求解和高斯分布假设,解决PDF方程的封闭问题,获得颗粒二阶矩模型,然后将颗粒应力方程简化成代数方程,建立代数应力模型.将对流扩散方程的有限分析法运用到求解两相流模型中,对壁面两相射流进行数值模拟,对比分析数值结果与实验结果.