非结构网格上的B-B单方程模型湍流计算

研究如何在非结构网格上进行Navier Stokes(N-S)方程湍流计算.采用格心有限体积方法离散N-S方程.为了适应非结构网格,计算所用的湍流模型特别选用Baldwin Barth(B-B)单方程模型.此模型由一个单一的具有源项的对流扩散方程组成.为了能在非结构网格上求解B B单方程模型,提出一显式有限体积格式,并直接对带源项的格式进行稳定性分析,得到了相应的时间步长限制条件.最后以平板、RAE 2822翼型、多段翼型绕流等数值算例验证了计算方法的有效性.     

基于MOF界面重构的多物质ALE方法

提出一种基于MOF(Moment-of-Fluid)界面重构的多物质ALE(Arbitrary Lagrangian-Eulerian)方法.流体力学方程组采用相容有限元方法进行空间离散.提出一种新的二维子网格力学模型,用来计算混合网格中的物理量经过一个拉氏步后发生的变化,混合网格内的界面重构采用MOF方法.提出一种精确积分守恒重映方法.给出数值算例,如空气和水的Riemann问题,Dukowicz问题,水中强激波与空气泡相互作用问题等.结果表明,方法具有较高的精度,能够处理物质界面和网格的大变形问题.     

一种基于特征理论模拟凝聚炸药爆轰的单元中心型Lagrange方法

提出一种数值模拟凝聚炸药爆轰问题的单元中心型Lagrange方法.利用有限体积离散爆轰反应流动方程组,基于双曲型偏微分方程组的特征理论获得离散网格节点的速度与压力,获得的网格节点速度与压力用于更新网格节点位置以及计算网格单元边的数值通量.以这种方式获得的网格节点解是一种"真正多维"的理论解,是一维Godunov格式在二维Riemann问题的推广.有限体积离散得到的爆轰反应流动的半离散系统使用一种显-隐Runge-Kutta格式来离散求解：显式格式处理对流项,隐式格式处理化学反应刚性源项.算例表明,提出的单元中心型Lagrange方法能够较好地模拟凝聚炸药的爆轰反应流动.     

隐-显积分因子间断Galerkin方法求解二维辐射扩散方程

提出一种求解二维非平衡辐射扩散方程的数值方法.空间离散上采用加权间断Galerkin有限元方法,其中数值流量的构造采用一种新的加权平均;时间离散上采用隐-显积分因子方法,将扩散系数线性化,然后用积分因子方法求解间断Galerkin方法离散后的非线性常微分方程组.数值试验中在非结构网格上求解了多介质的辐射扩散方程.结果表明:对于强非线性和强耦合的非线性扩散方程组,该方法是一种非常有效的数值算法.     

带源参数的热传导反问题的无网格方法

利用无网格有限点法求带有源参数的一维热传导反问题，推导了相应的离散方程.与其他基于网格的方法相比，有限点法采用移动最小二乘法构造形函数，只需要节点信息，不需要划分网格，用配点法离散求解方程，可以直接施加边界条件，不需要在区域内部求积分，减小了计算量.用有限点法求解热传导反问题具有数值实现简单、计算量小、可以任意布置节点等优点.最后通过算例验证了该方法的有效性.     

三维非结构网格Euler方程的LU-SGS算法及其改进

将LU-SGS隐式时间推进格式运用到非结构网格Euler方程的求解中,并对传统LU-SGS格式进行改进,结合网格重排序,发展了一套效率更高的三维Euler方程求解器.以M6机翼及超临界LANN机翼的跨音速无粘流场为算例,将改进LU-SGS格式与四步龙格-库塔显式格式及传统LU-SGS格式进行了比较.计算结果表明:所有格式的计算结果与实验结果都符合很好;传统LU-SGS格式计算效率为显式格式的3倍多,而改进LU-SGS格式计算效率为显式格式的7倍多.     

粒子输运离散纵标方程基于界面修正的并行计算方法

为了改造粒子输运方程求解的隐式格式,研究设计适应大型并行计算机的并行计算方法,介绍一类求解粒子输运方程离散纵标方程组的基于界面修正的源迭代并行计算方法.应用空间区域分解,在子区域内界面处首先采用迎风显式差分格式进行预估,构造子区域的入射边界条件,然后,在各个子区域内部进行源迭代求解隐式离散纵标方程组.在源迭代过程中,在内界面入射边界处采用隐式格式进行界面修正.数值算例表明该并行计算方法在精度、并行度、简单性诸方面均具有良好的性质.     

非定常Navier-Stokes方程基于两重网格离散的有限元并行算法

基于两重网格离散和区域分解技巧，提出三种求解非定常Navier-Stokes方程的有限元并行算法.算法的基本思想是在每一时间迭代步，在粗网格上采用Oseen迭代法求解非线性问题，在细网格上分别并行求解Oseen、Newton、Stokes线性问题以校正粗网格解.对于空间变量采用有限元离散，时间变量采用向后Euler格式离散.数值实验验证了算法的有效性.     

自适应间断有限元方法求解三维欧拉方程

将龙格库塔间断有限元方法(RDDG)与自适应方法相结合,求解三维欧拉方程.区域剖分采用非结构四面体网格,依据数值解的变化采用自适应技术对网格进行局部加密或粗化,减少总体网格数目,提高计算效率.给出四种自适应策略并分析不同自适应策略的优缺点.数值算例表明方法的有效性.     

半解析对偶棱边元及其在波导不连续性问题中的应用

给出电磁波导的对偶变量变分原理,并采用对偶棱边元对波导的横截面进行半解析离散. 将波导中沿纵向均匀的区段视为子结构,运用基于Riccati方程的精细积分算法求出其出口刚度阵,然后与不均匀区段的常规有限元网格拼装即可对波导不连续性问题进行求解. 半解析对偶棱边元的采用可以在最大程度上对有限元网格进行缩减,并且能够在不增加计算量的前提下任意增加子结构的长度,从而可以将截断求解区域的人工边界设置在距离不均匀区段充分远的地方,极大地减少了近似边界条件所带来的误差. 数值算例证明这种方法具有很高的精度与效率. 关键词： 波导的不连续性 半解析辛分析 对偶棱边元 精细积分     

Viscoelastic theory for nematic interfaces

A complete macroscopic theory for compressible nematic-viscous fluid interfaces is developed and used to characterize the interfacial elastic, viscous, and viscoelastic material properties. The derived expression for the interfacial stress tensor includes elastic and viscous components. Surface gradients of the interfacial elastic stress tensor generates tangential Marangoni forces as well as normal forces. The latter may be present even in planar surfaces, implying that in principle static planar interfaces may accommodate pressure jumps. The asymmetric interfacial viscous stress tensor takes into account the surface nematic ordering and is given in terms of the interfacial rate of deformation and interfacial Jaumann derivative. The material function that describes the anisotropic viscoelasticity is the dynamic interfacial tension, which includes the interfacial tension and dilational viscosities. Viscous dissipation due to interfacial compressibility is described by the anisotropic dilational viscosity, and it is shown to describe the Boussinesq surface fluid appropriate for Newtonian interfaces when the director is homeotropic. Three characteristic interfacial shear viscosities are defined according to whether the surface orientation is along the velocity direction, the velocity gradient, or the unit normal. In the last case the expression reduces to the interfacial shear viscosity of the Boussinesq surface fluid. The theory provides a theoretical framework to study interfacial stability, thin liquid film stability and hydrodynamics, and any other interfacial rheology phenomena.     

高温高压流体在线毛细管黏度测量

提出了一种适用于高温高压条件下的在线毛细管黏度测量方法,设计并搭建了相应的实验系统。通过与高温高压循环系统连接,可实现压力（约30MPa）、温度（约773K）范围内的流体黏度在线测量。利用该实验系统,在线测量了压力分别为12、23、25、28MPa下,温度303．46－383．28K范围内纯水的动力黏度。较之NIST文...     

An Analytical Study of Hydrodynamic Instability in the Flame. 1. Viscous Gas in the Flame Zone

The problem of the hydrodynamic stability of slow combustion is analytically solved with consideration given to the viscosity of the gas in the flame zone, the temperature dependence of the viscosity, and the dependences of the flame speed on the front curvature according to the Markstein model and on the pressure. The viscous forces in the flame zone alone cannot ensure the stability of the flame at any values of the Reynolds number. These forces act only as amplifiers of the stabilizing factor according to the Markstein model or in the case of a negative dependence of the flame velocity on the pressure. This property of internal friction forces is the more pronounced, the stronger the viscosity increases with the temperature. Thermal expansion is not only a destabilizing factor, leading to an increase in viscosity and other transport coefficients, but also produces a stabilizing effect.     

假塑性流体纳米压印中影响填充度的因素

作为新一代的半导体加工工艺, 直接金属纳米压印以其步骤简单、成本低等显著优点得到迅速的发展. 然而目前纳米压印中所采用的转移介质在流动状态下为牛顿流体, 牛顿流体的黏度是一个常量, 而假塑性流体具有黏度随着剪切速率的增大而逐渐减小的趋势, 更适用于纳米压印. 综合假塑性流体的剪切稀化特性以及直接金属图形转移的优点, 将不同大小的金属纳米粒子分散在基液中制成假塑性金属纳米流体并将其作为转移介质用于纳米压印中. 基于假塑性流体的Carreau流变模型利用COMSOL软件仿真分析金属纳米粒子假塑性流体参数集对图形压印转移的影响, 完成假塑性流体与牛顿流体分别作为转移介质实现图形转移的对比分析. 同时还得到了压印过程中影响填充度的各个因素, 如流体黏度、施加压强、掩模板移动速度等. 研究工作为金属纳米粒子假塑性流体制备以及纳米压印流程的设计提供了理论基础. 关键词： 纳米压印 假塑性流体 填充度     

Peristaltic Flow of Shear Thinning Fluid via Temperature-Dependent Viscosity and Thermal Conductivity

In this paper Williamson fluid is taken into account to study its peristaltic flow with heat effects. The study is carried out in a wave frame of reference for symmetric channel. Analysis of heat transfer is accomplished by accounting the effects of non-constant thermal conductivity and viscosity and viscous dissipation. Modeling of fundamental equations is followed by the construction of closed form solutions for pressure gradient, stream function and temperature while assuming Reynold's number to be very low and wavelength to be very long. Double perturbation technique is employed, considering Weissenberg number and variable fluid property parameter to be very small. The effects of emerging parameters on pumping, trapping, axial pressure gradient, heat transfer coefficient, pressure rise, velocity profile and temperature are analyzed through the graphical representation. A direct relation is observed between temperature and thermal conductivity whereas the indirect proportionality with viscosity. The heat transfer coefficient is lower for a fluid with variable thermal conductivity and variable viscosity as compared to the fluid with constant thermal conductivity and constant viscosity.     

Ion larmour radius effect on rf ponderomotive forces and induced poloidal flow in tokamak plasmas

Analytical approximations are used to clarify the effect of Larmour radius on rf ponderomotive forces and on poloidal flows induced by them in tokamak plasmas. The electromagnetic force is expressed as a sum of a gradient part and of a wave momentum transfer force, which is proportional to wave dissipation. The first part, called the gradient electromagnetic stress force, is combined with fluid dynamic (Reynolds) stress force, and gyroviscosity is included into viscosity force to model finite ion Larmour radius effects in the momentum response to the rf fields in plasmas. The expressions for the relative magnitude of different forces for kinetic Alfven waves and fast waves are derived.     

Gravitational collapse of dissipative fluid as a source of gravitational waves

Gravitational collapse of cylindrical anisotropic fluid has been considered in analogy with the work of Misner and Sharp. Using Darmois matching conditions, the interior cylindrical dissipative fluid (in the form of shear viscosity and heat flux) is matched to an exterior vacuum Einstein–Rosen space–time. It is found that on the bounding 3-surface the radial pressure of the anisotropic perfect fluid is linearly related to the shear viscosity and the heat flux of the dissipative fluid on the boundary. This non-zero radial pressure on the bounding surface may be considered as the source of gravitational waves outside the collapsing matter distribution.     

The selection of layer thicknesses to control acoustic radiation force profiles in layered resonators

Ultrasonic standing waves can be used to generate radiation forces on particles within a fluid. A number of authors have derived detailed representations of these forces but these are most commonly applied using an approximation to the energy distribution based upon an idealized standing wave within a mode based upon rigid boundaries. An electro-acoustic model of the acoustic energy distribution within a standing wave with arbitrary thickness boundaries has been expanded to model the radiation force on an example particle within the acoustic field. This is used to examine the force profile on a particle at resonances other than those predicted with rigid boundaries, and with pressure nodes at different positions. A simple analytical method for predicting modal conditions for combinations of frequencies and layer thickness characteristics is presented, which predicts that resonances can exist that will produce a pressure node at arbitrary positions in the fluid layer of such a system. This can be used to design resonators that will drive particles to positions other than the center of the fluid layer, including the fluid/solid boundary of the layer, with significant potential applications in sensing systems. Further, the model also predicts conditions for multiple subwavelength resonances within the fluid layer of a single resonator, each resonance having different nodal planes for particle concentration.     

Thermodynamics of the Bulk Viscous Cosmological Fluid in Presence of the Particle Creation Pressure

In the present work, we consider FRW metric and investigate some cosmological quantities in presence of bulk viscosity and particle creation pressure. The obtained results for a viscous cosmological fluid with particle creation show that the Hubble expansion parameter, energy density, bulk viscosity pressure, creation pressure and temperature depend on the particle creation rate and increase with increasing particle creation coefficient. It is found that the bulk viscosity and particle creation pressure seem to play important roles in the evolution of the early Universe.     

The effects of viscoelastic fluid on kinesin transport

Kinesins are molecular motors which transport various cargoes in the cytoplasm of cells and are involved in cell division. Previous models for kinesins have only targeted their in vitro motion. Thus, their applicability is limited to kinesin moving in a fluid with low viscosity. However, highly viscoelastic fluids have considerable effects on the movement of kinesin. For example, the high viscosity modifies the relation between the load and the speed of kinesin. While the velocity of kinesin has a nonlinear dependence with respect to the load in environments with low viscosity, highly viscous forces change that behavior. Also, the elastic nature of the fluid changes the velocity of kinesin. The new mechanistic model described in this paper considers the viscoelasticity of the fluid using subdiffusion. The approach is based on a generalized Langevin equation and fractional Brownian motion. Results show that a single kinesin has a maximum velocity when the ratio between the viscosity and elasticity is about 0.5. Additionally, the new model is able to capture the transient dynamics, which allows the prediction of the motion of kinesin under time varying loads.