自然对流换热的高精度非结构网格有限体积法

用非结构网格有限体积法求解自然对流换热时,传统的对流项离散格式难以兼顾数值精度与计算效率,我们发展了一种耦合高精度格式的延迟修正方法,用于对流项的离散.高Re数下方腔驱动流数值计算验证了该方法具有较高的计算精度和较好的稳定性.Boussinesq流体的自然对流换热数值模拟,表明该方法能有效克服高Ra数时数值计算发散,可准确捕捉自然对流换热问题中不同偏心率下的等温线和流线分布特征.     

改进的节块展开方法求解圆柱几何对流扩散方程

为高效求解球床高温气冷堆物理-热工耦合问题,发展改进节块展开法求解圆柱几何下的对流扩散方程.针对圆柱几何和对流扩散方程的特殊性,采用三阶多项式和指数函数作为r向横向积分方程的展开函数,在节块展开法的框架下高效求解对流扩散方程.数值验证表明,改进的节块展开方法具有固有的迎风特性,在使用粗网节块时依然能保持稳定性和较高的计算精度.     

静态溶液中KDP籽晶长大的浮力对流模拟研究

采用有限容积法,对KDP籽晶在静态溶液中的生长过程进行了数值模拟,研究了籽晶长大过程中的形状变化,考察了体过饱和度和籽晶特征尺寸对籽晶表面过饱和度及剪切力的影响,重点分析了籽晶长大过程中的尺寸效应.结果表明,当籽晶的特征长度小于临界尺寸时,籽晶的生长速率随着晶体尺寸的增大而加快;而当籽晶的特征长度大于临界尺寸时,籽晶的生长速率随着晶体尺寸的增大而减慢.     

Evolution of a two-mode squeezed vacuum in the amplitude dissipative channel

For the first time we derive the dissipating result of an initial two-mode squeezed pure vacuum state passing through a two-mode amplitude dissipative channel described by the direct product of two independent single-mode master equations. Although these two master equations do not mix the two modes (there is no coupling between them), since the two-mode squeezed state is simultaneously an entangled state, the final state which emerges from passing this channel is a two-mode mixed density operator. The compact expression of the outcoming state is obtained, which manifestly shows that as time evolves, the squeezing effect decreases.     

Conduction-radiation effect on transient natural convection with thermophoresis

This paper discusses the effect of thermophoretic particle deposition on the transient natural convection laminar flow along a vertical flat surface,which is immersed in an optically dense gray fluid in the presence of thermal radiation.In the analysis,the radiative heat flux term is expressed by adopting the Rosseland diffusion approximation.The governing equations are reduced to a set of parabolic partial differential equations.Then,these equations are solved numerically with a finite-difference scheme in the entire time regime.The asymptotic solutions are also obtained for sufficiently small and large time.The obtained asymptotic solutions are then compared with the numerical solutions,and they are found in excellent agreement.Moreover,the effects of different physical parameters,i.e.,the thermal radiation parameter,the surface temperature parameter,and the thermophoretic parameter,on the transient surface shear stress,the rate of surface heat transfer,and the rate of species concentration,as well as the transient velocity,temperature,and concentration profiles are shown graphically for a fluid(i.e.,air) with the Prandtl number of 0.7 at 20 C and 1.013 × 10 5 Pa.     

脑Willis环的一维血流动力学及氧输运特性的数值研究

Willis环是大脑侧枝循环的重要组成部分, 研究其血流动力学特性以及氧输运规律对脑缺血疾病的认知和预防有着非常重要的作用. 该文旨在利用一维血流动力学模型模拟整个Willis环的流量变化和压力分布, 并建立动脉内氧输运的一维模型以模拟Willis 环内氧分压的变化规律, 为研究脑组织内血液流动和氧输运打下基础. 首先, 基于弹性圆管内的一维非线性流动方程和状态方程建立血流动力学模型, 在一维对流扩散方程的基础上, 考虑由管腔向壁面的扩散和壁面细胞的新陈代谢消耗推导出氧输运特性方程. 通过 Lax-Wendroff两步法对血流动力学方程进行离散, 而在进行对流扩散方程的离散时, 则运用迎风格式. 通过数值计算得到了正常情况下Willis环各个血管任意位置的流量、压力和氧分压的变化曲线, 正常情况下各个位置的氧分压处于稳定的平衡状态. 最后, 还通过此模型进一步模拟了右侧颈内动脉狭窄对各个血管内流动的影响. 当狭窄程度达到80%时, 中脑动脉的流量和压力会明显下降, 造成其供应区域的血流减少. 同时, Willis环右侧血管内的氧分压会大大降低, 而左侧血管的氧分压会出现上升趋势, 但幅度要小于右侧血管降低的幅度.     

一种新的求解对流扩散边值问题的无网格方法

基于配点法和楔形基函数,提出了一种新的求解对流扩散边值问题的无网格方法。通过一维和二维的问题验证了该数值方法的可行性;并根据数值算例和分析,可以看到该数值方法能达到满意的收敛效果。该数值方法的隐格式形式能够有效地消除对流占优问题的数值振荡现象,是一种真正的无网格方法。     

流过半无限竖直可渗透壁面时的磁流体动力学自由对流

研究不可压缩粘性导电流体,流过半无限竖直可渗透平板时,将其偏微分形式的流动和传热的基本控制方程,应用适当的相似变换,简化为非线性的常微分方程组．对两种抽吸参数:大的和小的抽吸参数,采用摄动法得到变换后方程的近似解．数值结果表明,随着磁场参数和抽吸参数的增大,任意点的速度场在减小;磁场参数的影响,引起热边界层厚度的增大;速度和温度场随着热汇参数的增大而减小．     

考虑Soret和Dufour效应的方腔内双扩散自然对流格子Boltzmann模拟

采用格子Boltzmann方法,考虑Soret和Dufour效应,对内置高浓度发热圆的方腔内部双扩散自然对流现象进行数值模拟.高浓度发热圆位于方腔中心,四周壁面均为低温低浓度.在该模型中,用三个独立的LBGK方程分别模拟速度场、温度场和浓度场,并通过Boussinesq近似将它们耦合起来.分析Soret数和Dufour数对方腔内部双扩散自然对流的影响,得到流线图、等温线图、等浓度线图、发热圆表面平均Nusselt数和平均Sherwood数.结果表明：Soret和Dufour效应对方腔内双扩散自然对流影响明显,不能忽略.     

不同驱动力对Cz硅生长中熔体对流影响的数值模拟研究

为了确定Cz单晶硅生长各种驱动力对熔体对流及固/液界面形状的影响,利用CGSim软件,对典型的Cz单晶硅生长中的熔体对流进行数值模拟.研究了重力、表面张力、平流力、晶转、埚转和氩气剪切力等各种驱动力的大小对熔体对流涡胞、涡胞强度、界面形状、温度分布的影响.结果表明:各种驱动力对熔体对流的影响大小依次为:浮力＞表面张力＞晶转力＞氩气剪切力＞埚转力＞平流力;浮力和表面张力使熔体产生一沿坩埚壁上升、从固/液界面附近下降的涡胞,晶转力和氩气剪切力使熔体产生与前面反方向的涡胞,而埚转力产生多个不同流向的对流涡胞,使熔体混合更加均匀,熔体凝固引起的平流力对熔体对流影响不大;增大埚转,熔体中涡胞数量更多、对流换热更充分、温度梯度更小、熔体内的最高温度更低,有利于减少石英坩埚氧的熔解,但界面更向下凹;增大晶转,熔体内的最高温度无明显变化;固定埚转Ωc=-10 r/min,晶转存在一临界值Ωs(C)=60～ 80 r/min,当Ωs＜Ωs(C)时,增大晶转,固/液界面更向上凹,当Ωs＞Ωs(C)时,增大晶转,固/液界面更向下凹.