稀薄流非线性模型方程离散速度坐标法有限差分解

从一般非线性Bo ltzm ann方程出发,发展并实现了一套适于大范围K nudsen数稀薄流问题数值模拟的统一算法。采用BGK模型和Shakov模型近似碰撞项,进而引入两个二速度无量纲简化分布函数,通过关于分子速度第三分量取矩积分,将三速度单一模型方程变换为二速度微分方程组。基于G auss-H erm ite积分公式和正交多项式G auss积分公式,借助离散速度坐标法消除简化模型方程对分子速度空间的连续依赖性,从相空间到物理空间得到一组带源项双曲守恒离散方程,并给出其显式和隐式二阶迎风TVD有限差分解。以二维圆柱A r气体超声速绕流算例,验证了数值算法的有效性,比较分析了漫反射和镜面反射两种气体分子壁面反射模型的计算结果。     

杜芬方程的1／3纯亚谐解及其过渡过程

通过谐波平衡法和数值积分法研究了杜芬方程的1／3纯亚谐解,对参数变化的过渡过程的敏感性和扰动初始值的过滤进行了研究,同时,还对谐波平衡法和数值积分法的精度进行了比较,也考虑了亚谐波成分的渐近稳定性。     

非线性抛物型方程有限元法数值积分的有效性

Abstract. The effect of numerical integration in finite element methods applied to a class of nonlinear parabolic equations is considered and some sufficient conditions on the quadrature scheme to ensure that the order of convergence is unaltered in the presence of numerical integration are given. Optimal Lz and H1 estimates for the error and its time derivative are established.     

载流直螺线管磁场的三维数值解

本文利用严格的解析和数值方法讨论了载流直螺线管产生的三维磁场 ,并在计算机上进行了高精度的数值计算。     

块结构自适应网格上任意区域和任意界面的数值积分

研究块结构自适应网格上计算任意界面上和任意区域内的数值积分方法,其中任意界面和任意区域通过一个水平集函数表示.首先介绍在一致网格上任意界面上和任意区域内的数值积分方法.然后,将该方法推广到块结构的自适应网格上.数值算例表明,自适应网格方法有二阶精度.同一致网格方法相比,自适应网格方法显著地减少了计算机存储量的需求.     

有限元实现含应变非强化区的双面循环塑性模型

在有限单元方法的框架下,由广义径向返回法来实现含有应变非强化区的循双面环塑性本构模型.在小变形的假设下,本文详细描述了有限元方法实现率无关的双面循环塑性本构模型的完整算法.采用向后Euler隐式算法进行离散积分,其最终可化简为一个非线性方程,并且由Peg-asus方法来求解.对径向返回法进行线性化,可以得到切向刚度算子.最后,文中通过两个算例表征了离散积分的精确性和材料模型的正确性.     

一个新的高精度二重数值积分公式

本文将文献[1]的结果推广到二维情形,得到了相应的计算二重积分的高精度求积公式及其复合型公式。它具有辛普生公式的优点,但代数精确度比辛苦生公式提高2阶。 本文数值试验表明,这个新的高精度数值积分公式是十分有效的。     

直流圆形线圈径向磁感应强度研究

基于毕奥-萨伐尔定律，分别考虑及不考虑直流通电圆形线圈的径向宽度和纵向厚度，推导出两种情况下该线圈沿径向方向的磁感应强度分布公式.然后，利用Matlab软件对获得的分布公式进行数值积分，得到理论上的磁感应强度分布数据.为了考察理论值的正确性，通过设计自制圆形线圈，进行实际测量.结果表明，考虑径向宽度和纵向厚度的影响时，磁感应强度实际测量值与理论值的相对误差较小.通过理论计算与实验研究相结合考察直流通电圆形线圈沿径向方向的磁感应强度分布，既能消除电磁感应法测交变磁场实验中磁场分布的系统误差，还可加强学生对电磁理论的理解与认识，提高电磁学课程的教学效果.