A finite volume scheme preserving extremum principle for convection–diffusion equations on polygonal meshes

We propose a nonlinear finite volume scheme for convection–diffusion equation on polygonal meshes and prove that the discrete solution of the scheme satisfies the discrete extremum principle. The approximation of diffusive flux is based on an adaptive approach of choosing stencil in the construction of discrete normal flux, and the approximation of convection flux is based on the second‐order upwind method with proper slope limiter. Our scheme is locally conservative and has only cell‐centered unknowns. Numerical results show that our scheme can preserve discrete extremum principle and has almost second‐order accuracy. Copyright © 2017 John Wiley & Sons, Ltd.     

三维多面体网格上扩散方程的保正格式

 针对三维任意(星形)多面体网格, 本文构造了扩散方程的一种单元中心型非线性有限体积格式, 证明了该格式具有保正性. 在该格式设计中, 除引入网格中心量外, 还引入网格节点量和网格面中心量作为中间未知量, 它们将用网格中心未知量线性组合表示, 使得格式仅有网格中心未知量作为基本未知量. 在节点量计算中, 利用网格面上的调和平均点, 设计了一种适用于三维多面体网格的局部显式加权方法. 该格式适用于求解非平面的网格表面和间断扩散系数的问题. 数值例子验证了它对光滑解具有二阶精度和保正性.     

Effect of pH value on the growth morphology of KH2PO4 crystal grown in defined crystallographic direction

KH₂PO₄ single crystals were grown in aqueous solution at different pH values by using “point seeds” with a defined crystallographic direction at 59 degree to the Z axis. Atomic Force Microscope (AFM) was applied to observe the surface morphology of (100) face. It was found that at the same supersaturation, the larger steps appeared at the lower pH value before appearance of 2D nucleus. We found that 2D nucleus was occurred at σ ≤ 0.04 when pH value is <2.8. The occurrence of 2D nucleus was caused by the decreasing step‐edge free energy with the decreasing of pH value in the growth solution. In this paper, we observed the morphologies of (100) faces of KDP crystals which grew in solutions with different pH values. 2D nucleuses appeared on the terrace of growth steps when pH value down to 2.8 and 3.2 at supersaturation of 0.04, while pH value down to 2.4, only 2D nucleation control the growth. Therefore, the pH value can change the growth mechanism of KDP crystals.     

基于贝叶斯参数优化的无模型自适应硅单晶直径控制

直拉硅单晶的生长过程涉及多场多相耦合与复杂的物理化学变化,其中工艺参数的波动是导致晶体直径不均匀的重要原因,如何实现工艺参数的控制以获得理想的、均匀的晶体直径具有重要的研究意义。本文分析现有控制方法存在不稳定以及控制效果不佳的问题后,提出基于贝叶斯参数优化的无模型自适应控制模型来控制硅单晶生长过程中的晶体直径。首先以坩埚上升速度与加热器的功率作为控制输入参数,晶体直径作为输出,搭建无模型自适应控制模型,并分析算法的稳定性。其次将控制模型进行仿真实验,发现硅单晶直径控制模型中不同的超参数设定会影响控制过程的迭代次数以及控制效果。最后,利用贝叶斯优化超参数的取值范围,并进行最终的仿真实验,结果表明,经贝叶斯参数优化后的控制模型计算快、迭代次数少,输出的晶体直径稳定,同时将生长工艺参数控制在实际生产要求范围内。因此,基于贝叶斯参数优化的无模型自适应控制实现了硅单晶直径均匀稳定的有效控制,具有结合工程背景的实际应用前景。     

红黑排序混合算法收敛速度分析

The algorithm of applying the block Gauss elimination to the Red-Black or-dering matrix to reduce the order of the system then solve the reduced system byiterative methods is called Hybrid Red-Black Ordering algorithm.In this paper,we discuss the convergence rate of the hybrid methods combined with JACOBI,CG,GMRES(m).Theoretical analysis shows that without preconditioner thesethree hybrid methods converge about 2 times as fast as the corresponding natural ordering methods.For the case that all the eigenvalues is near the real axis, the GMRES(m) algorithm converges about 3 times faster than the natural ordering GMRES(m).Various numerical experiments are presented.For large scale prob-lem with preconditioners, numerical experiments show that the GMRES(m) hybrid methods converge from about 3 times to even 5 times as fast as the natural order-ing methods and the computing time is reduced to about 1/3 even 1/6 of that of the natural ordering methods.     

非线性抛物组非均匀网格差分解的唯一性和稳定性

１．引言 １．对一维非线性抛物组,在文献山中已构造一般非均匀网格差分格式,其中差分逼近的组合系数对不同的网格点和不同的网格层可以不同,并且运用不动点原理证明了差分解的存在性和收敛性．在非均匀网格差分格式中差分逼近的组合系数为常数的情形,文献［２］证明了具有有界二阶差商的离散向量解的存在性、唯一性和稳定性．本文将对文献［１］中构造的一般非均匀网格差分格式,证明所得到的差分解的唯一性和稳定性． 考虑如下非线性抛物组其中是未知的ｍ－维向量函数是给定的矩阵函数,ｊ（ｘ,ｔ,ｕ,ｐ）。是给定的ｍ－维向量函数…     

带有动力学条件的高维一相Stefan问题

本文我们证明了带有动力学条件的高维一相Ｓｔｅｆａｎ总题局部古典解的存在唯一性。     

A Finite Volume Method Preserving Maximum Principle for the Conjugate Heat Transfer Problems with General Interface Conditions

In this paper, we present a unified finite volume method preserving discrete maximum principle (DMP) for the conjugate heat transfer problems with general interface conditions. We prove the existence of the numerical solution and the DMP-preserving property. Numerical experiments show that the nonlinear iteration numbers of the scheme in [24] increase rapidly when the interfacial coefficients decrease to zero. In contrast, the nonlinear iteration numbers of the unified scheme do not increase when the interfacial coefficients decrease to zero, which reveals that the unified scheme is more robust than the scheme in [24]. The accuracy and DMP-preserving property of the scheme are also veried in the numerical experiments.