一类一维的辐射流体力学方程组激波的存在性

该文主要讨论一维空间中一类辐射流体力学方程组的激波.由Rankine-Hugoniot条件及熵条件得此问题可表述为关于辐射流体力学方程组带自由边界的初边值问题.首先通过变量代换,将其自由边界转换为固定边界,然后研究关于此非线性方程组的一个初边值问题解的存在唯一性.为此先构造了此问题的一个近似解,然后分别通过Picard迭代与Newton迭代对此非线性问题构造近似解序列.通过一系列估计与紧性理论得到此近似解序列的收敛性,其极限即为原辐射热力学方程组的一个激波.

辐射流体力学Lagrange方程组一类人为解构造方法

针对一维Lagrange辐射流体力学方程组,基于物理量在Euler空间和Lagrange空间的微分关系,提出了一种人为解构造方法,并构造了一类一维Lagrange辐射流体力学方程组人为解.构造的人为解在整个计算区域光滑可微,质量方程无源项.将构造的人为解应用至二维辐射流体力学Lagrange程序中,从数值误差、收敛阶方面验证了程序的正确性,展示了人为解的可行性和适用性.

二维无界区域中不可压非牛顿流体力学方程组的整体吸引子(英文)

研究一类二维无界区域中的等热双极不可压粘性非牛顿流体力学方程组,通过证明相应的解半群的紧性,得到整体吸引子的存在性.

一维可压流体力学方程组整体经典解

本文研究了一维可压流体力学方程组的Cauchy问题,利用特征线法和分析的技巧,在一些合理的假设下,得到了Cauchy问题整体经典解存在的一个充分必要条件.

考虑中子输运对介质流体力学运动发生影响的辐射流体力学、核子数变化方程组

带电粒子平均自由程远比中子的小,在一级近似下,可以把系统中的带电粒子看成随时都处于局部热平衡状态,而考虑中子在较大范围内所进行的质量、动量、能量的输运.这种输运过程要影响到介质的流体力学运动.本文从玻尔兹曼方程出发,导出了考虑中子输运对介质流体力学运动发生影响的辐射流体力学方程组,核子数变化方程组.

二维三温辐射流体力学拉氏计算中的子网格压力方法

研究对应于电子压、离子压、光子压的扰动隅角力的计算方法,提出二维三温辐射流体力学拉氏计算中的两种子网格压力方法,即基于状态方程和基于几何的方法.数值试验表明,这两种方法均能很好地抑制二维三温辐射流体力学拉氏计算中出现的网格非物理畸变.

单丝Z箍缩的一维单温辐射磁流体力学计算

介绍在正弦电流波形下铝丝的内爆动力学过程的数值模拟计算 ,计算采用一维单温辐射磁流体力学模型 ,能流中除了常规热传导外 ,也包括黑体辐射流项 ,以考虑辐射对内爆动力学过程的影响。计算了用于箍缩预脉冲的低电流 (幅值 1 0～ 5 0kA)和较高电流 (幅值 40 0kA)情况下单丝的内爆结果。

二维三温辐射(磁)流体力学拉氏计算中的人工粘性分裂技术

针对使用边粘性和张量粘性进行二维三温辐射(磁)流体力学拉氏计算时激波处离子温度出现非物理尖峰的问题,提出一种人工粘性分裂技术,将人工粘性耗散能量按照电子压和离子压的比分成两部分,并分别加在电子内能和离子内能上.在二维三温辐射磁流体力学程序MARED上实现该技术,并对Z-pinch的内爆动力学过程进行数值模拟.消除了模拟结果中激波处离子温度的尖峰.

移动网格方法求解二维流体力学方程组

移动网格方法是求解发展方程的一种方法，它根据物理解的变化，随时调整网格的疏密。在流体力学问题中，由于激波、接触间断等变化剧烈的区域在全流场中只占很小一部分，使用移动网格方法可以在不减少计算精度的情况下，较大地减少计算量。

流体力学方程的间断有限元方法

在二维区域三角形网格上应用一阶、二阶和三阶精度间断有限元方法,对流体力学方程和方程组进行了数值模拟.计算结果与差分方法计算结果比较,认为间断有限元方法在求解复杂边界条件和区域问题上有一定的优势.

一种二维三温辐射热传导方程组SFVEM格式的并行预条件子

在三角形剖分下,针对二维三温辐射热传导方程组的SFVEM离散系统,设计了一类基于块对角逆的预条件,给出了条件数的理论估计式.在此基础上,设计了相应的PCG算法.数值实验验证了理论结果的正确性.

一种求解二维辐射流体力学方程组的显隐式格式

构造一种求解二维辐射流体力学方程组的有限体积方法。相较于Euler方程组,辐射流体力学方程组的数值格式设计更为困难。不仅辐射压力与辐射能量的强非线性增加了数值计算的难度,而且求解强激波问题也是一大难点。与此同时,物质量以声速传播,辐射量以光速传播也增加了该系统求解的难度。为了克服这些难点,我们使用MUSCL-Hancock方法求解模型的双曲部分,利用TR/BDF2的变换格式求解其扩散部分。给出了灰非平衡扩散模型关于时间的渐近分析。最后,通过数值实验展现该格式的表现。

基于变分原理的二维热传导方程差分格式

研究二维热传导方程的差分数值模拟 .用变分原理在不规则结构网格上建立热流通量形式的差分格式 .将热流通量作为未知函数求泛函极值 ,并与温度函数联立求解 .克服通常九点格式用插值方法计算网格边界上的热传导系数和网格结点上的温度所引入的误差 .

通过多孔介质的二维磁流体力学方程组的全局正则性

本文的目的是研究以不断吹/吸方式通过孔道流体的二维磁流体动力学方程(MHD)弱解的全局正则性.为此,我们将利用抛物正则化过程和在流动模型中流体的的Darcy定律.

辐射扩散方程的自适应坐标变换方法

借鉴自适应坐标变换流体力学方法和其他一些有益工作，研究辐射扩散方程的自适应坐标变换方法，并考虑辐射扩散和流体力学耦合的辐射流体力学问题。对辐射扩散方程和一般二维对流扩散方程，从积分和微分形式出发，导出自适应坐标变换下的等价形式，对数值离散格式进行初步讨论。对二维三温辐射扩散方程，采用时空有限体积方法，进行数值离散；利用局部函数近似离散扩散算子。利用合理的进程分裂，简化耦合问题的计算。

重新初始化的LS方法跟踪二维可压缩多介质流界面运动

用非耦合求解方法计算LevelSet函数方程与流体力学方程组 ,应用重新初始化的LevelSet函数确保距离函数性质 ,流体力学方程组采用二阶精度多介质流波传播差分格式计算 ,重新初始化方程采用五阶WENO格式计算。并给出了二维可压缩多介质流界面运动的计算结果。

黏性依赖于密度的一维可压缩黏性辐射反应气体方程组的Cauchy问题

本文考虑一维可压缩黏性辐射反应气体方程组,对一类依赖于密度的、退化的黏性系数,得到了其Cauchy问题大初值非真空整体解的存在性.问题的关键在于如何得到密度和温度的正的上下界估计.

衬套内爆ALE方法二维MHD数值模拟

早期关于二维磁流体力学方程组的描述都是在Lagrangian或Eulerian坐标内进行的。在Lagrangian描述中，计算网格固定在物体上随物体一起运动，网格点与物质点在物体的变形过程中始终保持一致，因此能准确描述物体的移动界面。但对于大变形问题，物质的扭曲将导致计算网格的畸形而使计算无法进行下去。在Eulerian描述中，网格固定在空间中，因而计算网格在物体的变形过程中保持不变，因此很容易处理物质的扭曲。但对于运动界面，需要引入非常复杂的数学映射，将可能导致较大的误差。

Kr喷气箍缩等离子体的数值研究

探讨了Z箍缩等离子体辐射磁流体模拟二维三温物理模型及其数值计算方法,针对"强光一号"加速器上Kr喷气实验的具体条件,利用辐射磁流体力学程序Zpinch2D DG模拟了Kr喷气等离子体聚爆过程,分析了X射线辐射的特点,给出了等离子体密度和温度的演化图像以及喷气箍缩中一些带有普遍性的结论。

在不规则四边形网格上逼近扩散算子的五点及九点差分格式的测试

在Lagrange坐标下使用四边形网格进行二维辐射流体力学数值计算的难点之一是需要构造在不规则四边形网格上仍能较好地逼近扩散算子的差分格式 .本文就五点差分格式和目前常用的九点差分格式进行了比较全面的数值测试和理论分析 .结果表明五点格式仅在均匀矩形网格上具有二阶逼近精度 ,九点格式仅在均匀平行四边形网格上具有二阶逼近精度 ,这两种格式在一般的不规则四边形网格上通常都是不相容的 .尽管九点格式优于五点格式 ,但它对不规则网格的适应性远不如人们以前所想象的那么好 .由此可见 ,为了进一步改进二维辐射流体力学的数值计算 。