粘弹性二阶流体混合层流场拟序结构的数值研究

本文用拟谱方法对随时间发展的二维粘弹性二阶流体混合层流场进行了直接数据值模拟,给出在高雷诺数和低Deborah数下大涡的卷起、配对和合并等过程,通过与相同雷诺数下牛顿流体的比较,揭示了弱粘弹性对混合层中大涡拟序结构演变的影响.     

直圆管突扩通道内宾汉流体湍流流场的数值研究

本文依据牛顿流体中建立的标准ｋ＿ε湍流模型这一基本思想，考虑宾汉流体的本构方程，建立了适用于求解宾汉流体湍流流动的控制方程·采用压力修正算法，实现了宾汉流体速度场与压力场的关联·在理论研究基础上，对直圆管突扩通道内宾汉流体湍流流动进行了数值研究，并探讨了直圆管突扩通道内宾汉流体湍流流动机理·     

金属泡沫微流道热沉内流体流动与传热特性的数值研究

为改善高能量密度电子设备的冷却效率,提出了在微流道热沉内填充金属泡沫的新型热沉结构,并数值研究了金属泡沫的孔隙率、孔密度、材质(铜、镍及铝)、流体工质(水、乙二醇及纳米流体)等相关参数对微流道流动与换热特性的影响.研究结果表明:金属泡沫可以显著地强化微流道热沉的换热特性;添加金属泡沫后微流道热沉的换热性能可提高2倍以上;采用纳米流体与金属泡沫相结合的双重强化换热手段可以进一步地增强微流道热沉的冷却能力;在层流流动状态下金属泡沫微流道热沉可以对发热量为200 W/cm2的电子设备进行有效地冷却,表明其在高功率密度电子设备热管理领域具有广阔的应用前景.     

近临界流体微通道内流动稳定性和换热特性研究

微尺度条件下的化工、医药、传热与能源利用等系统的研究已经成为极具潜力和挑战性的课题.相应条件下流体流动和换热的分析必须考虑尺度效应所带来的系列问题.该研究采用了数值模拟方法对近临界二氧化碳流体在微尺度通道内的流动稳定性和换热特性进行了探索.研究发现,在近临界区域内由于流体较强的膨胀特性和较低的热扩散特性,在微尺度几何条件下会产生瞬态不稳定的漩涡流动.该种条件下微尺度对流换热和混合效率都得到了大幅提高.进一步,研究针对微尺度局部稳定性演化进行了机理分析并应用了参数估计,总结获得了微通道内近临界流体瞬态换热和混合的基本特性.     

三维磁微极流体方程弱解的正则性

本文给出了磁微极流体方程弱解的一个新的正则性准则:如果u满足uz ∈Lq(0,T;Lp(R3)),其中p≥3且满足3/p+2/q≤1,那么弱解(u,ω,b)在(0,T)是光滑解.     

Bingham流体数值模拟的双罚函数—正交投影的隐式求解方法

本文对Bingham流体给出了双罚函数逼近和正交投影的隐式求解方法.这个方法把Bingham流体处理为承受不等式应力约束的Newton流体的逼近解.有效地模拟了可以出现流动或不流动的“刚性核”的Bingham流体的流动问题.     

服从OldroydB型微分模型的粘弹性流体问题的数值解法

本文对服从OldroydB型微分模型的粘弹性流体问题给出了一种数值逼近算法．该算法对压力方程采用标准混合有限元方法,对速度方程采用并行非重叠区域分解方法和特征线法．这种并行算法在子区域上用Galerkin方法,通过积分平均方法显式地给出内边界的数值流．在本文最后还给出了该算法的最优L^2。一误差估计．     

热传导的微极固体和流体介质界面上波的传播

研究微极广义热弹性固体半空间和热传导微极流体半空间界面上波的传播.讨论微极广义热弹性固体半空间和热传导微极流体半空间之间平面界面上,斜向入射平面波的反射和透射现象.假设入射波穿过微极广义热弹性固体,射向平面界面后传播.得到了封闭形式的、不同反射和透射波的波幅比,它们是入射角、频率的函数,并为介质的弹性性质所影响.对一些特定的类型,显示出微极和热松弛对波幅比的影响.还从本文的研究中推演出一些早期工作的结果.     

微极流体向受热面的MHD驻点流动

分析了有均匀横向磁场作用时,导电微极流体垂直冲击受热面时形成的二维驻点流动问题．应用适当的相似转换,将连续、动量、角动量及热量的控制方程,及其相应的边界条件,简化为无量纲形式．然后,利用以有限差分离散化为基础的算法,求解简化了的自相似非线性方程．用Richardson外推法,进一步求精其结果．以图表形式表示磁场参数、微极性参数、Prandtl数对流动和温度场的影响,说明了其解的重要特性．研究表明,随着磁场参数的增大,速度和热边界层厚度变小了．与Newton流体相比较,微极流体的剪应力和传热率出现明显的减少,这对聚合物生产过程中流体的流动和热量控制是有益的．     

化学反应对竖直平板边界磁流体动力学微极流体滑流的影响

对半无限竖直平板为边界的多孔介质，研究了传热、传质对微极流体不稳定滑流的影响，其化学反应是一级均匀的。均匀磁场垂直作用于可以吸收微极流体的多孔表面，吸引速度随着时间而变化。自由流动的速度随着微小扰动而呈指数增大或减小。采用近似方法获得了微极流体的速度、微转动、温度、浓度的表达式，还得到了在不同流体特征和流动条件下，壁面的摩擦系数、耦合应力系数、传热率和传质率。     

不可压缩流问题低次元稳定有限体积数值方法研究

分析了R^d,d=2,3维不可压缩流Stokes问题低次元稳定有限体积方法,它主要利用局部压力投影方法对两种流行但不满足inf-sup条件的有限元配对(P_1-P_0和P_1-P_1)在有限体积方法的框架下进行稳定;利用有限元与有限体积方法的等价性进行有限体积方法理论分析.结果表明不可压缩流Stokes问题在f∈Hd,d=2,3维不可压缩流Stokes问题低次元稳定有限体积方法,它主要利用局部压力投影方法对两种流行但不满足inf-sup条件的有限元配对(P_1-P_0和P_1-P_1)在有限体积方法的框架下进行稳定;利用有限元与有限体积方法的等价性进行有限体积方法理论分析.结果表明不可压缩流Stokes问题在f∈H¹情况下,本文方法得到的解与稳定有限元方法解之间具有O(h1情况下,本文方法得到的解与稳定有限元方法解之间具有O(h²)阶超收敛阶结果,且稳定有限体积方法取得了与稳定有限元方法相同的收敛速度,与稳定有限元方法比较,稳定有限体积方法计算简单高效,同时保持物理守恒,因此在实际应用中具有很好的潜力。     

A Variational Analysis for the Moving Finite Element Method for Gradient Flows

By using the Onsager principle as an approximation tool, we give a novel derivation for the moving finite element method for gradient flow equations. We show that the discretized problem has the same energy dissipation structure as the continuous one. This enables us to do numerical analysis for the stationary solution of a nonlinear reaction diffusion equation using the approximation theory of free-knot piecewise polynomials. We show that under certain conditions the solution obtained by the moving finite element method converges to a local minimizer of the total energy when time goes to infinity. The global minimizer, once it is detected by the discrete scheme, approximates the continuous stationary solution in optimal order. Numerical examples for a linear diffusion equation and a nonlinear Allen-Cahn equation are given to verify the analytical results.     

周期结构区域中压电问题的双尺度有限元方法

近年来纤维压电复合材料的力电性能预测已发展为一个重要的研究领域．对力电耦合周期结构的复合材料问题,通过引入匹配的边界层得到了电势与位移解的新型双尺度有限元计算方法,建立了电势与位移的双尺度耦合关系,分析了双尺度有限元解的误差.数值算例验证了方法的有效性．     

Numerical simulations of the improved Boussinesq equation

In this study, numerical simulations of the improved Boussinesq equation are obtained using two finite difference schemes and two finite element methods, based on the second‐and third‐order time discretization. The methods are tested on the problems of propagation of a soliton and interaction of two solitons. After the L_∞ error norm is used to measure differences between the exact and numerical solutions, the results obtained by the proposed methods are compared with recently published results. © 2009 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 2010     

求解超声速无粘绕流的杂交型有限元格式

本文从一维双曲型标量方程出发,以一个普通二阶有限元格式及由其单边对角化导出的一阶单调型格式为基础,构造出一种具有单调性的杂交型有限元格式.为了向二维Euler方程组情形推广,所采用的开关函数是基于流场梯度的局部函数,并专门考虑了相邻单元的影响.二维情形的算例表明新格式可以明显抑制激波附近的振荡.     

Numerical study of natural superconvergence in least-squares finite element methods for elliptic problems

Natural superconvergence of the least-squares finite element method is surveyed for the one-and two-dimensional Poisson equation. For two-dimensional problems, both the families of Lagrange elements and Raviart-Thomas elements have been considered on uniform triangular and rectangular meshes. Numerical experiments reveal that many superconvergence properties of the standard Galerkin method are preserved by the least-squares finite element method. The second author was supported in part by the US National Science Foundation under Grant DMS-0612908.     

Numerical results for the Klein‐Gordon equation in de Sitter spacetime

We consider the initial value problem for the Klein‐Gordon equation in de Sitter spacetime. We use the central difference scheme on the temporal discretization. We also discretize the spatial variable using the finite element method with implicit and the Crank‐Nicolson schemes for the numerical solution of the initial value problem. In order to show the accuracy for the results of the solutions, we also examine the finite difference methods. We observe that the numerical results obtained by using these methods are compatible.     

Numerical methods and error analysis for one dimensional elliptic problems containing singularities

The method of auxiliary mapping (MAM), introduced by Babu?ka and Oh, was proven to be very successful in dealing with monotone singularities arising in two‐dimensional problems. In this article, in the framework of the p‐version of FEM, MAM is presented for one‐dimensional elliptic boundary value problems containing singularities. Moreover, in order to show the effectiveness of MAM, a detailed proof of an error estimate is also presented, which gives a sharp error bound of MAM. © 2003 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 19: 399–420, 2003.     

A Two‐Parameter Stabilized Finite Element Method for Incompressible Flows

Based on two‐grid discretizations, a two‐parameter stabilized finite element method for the steady incompressible Navier–Stokes equations at high Reynolds numbers is presented and studied. In this method, a stabilized Navier–Stokes problem is first solved on a coarse grid, and then a correction is calculated on a fine grid by solving a stabilized linear problem. The stabilization term for the nonlinear Navier–Stokes equations on the coarse grid is based on an elliptic projection, which projects higher‐order finite element interpolants of the velocity into a lower‐order finite element interpolation space. For the linear problem on the fine grid, either the same stabilization approach (with a different stabilization parameter) as that for the coarse grid problem or a completely different stabilization approach could be employed. Error bounds for the discrete solutions are estimated. Algorithmic parameter scalings of the method are also derived. The theoretical results show that, with suitable scalings of the algorithmic parameters, this method can yield an optimal convergence rate. Numerical results are provided to verify the theoretical predictions and demonstrate the effectiveness of the proposed method. © 2016 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 33: 425–444, 2017