非均匀剪切流场中液晶聚合物微观结构的无网格模拟

基于宏观流场控制方程与微观分子取向扩散方程耦合的微-宏观双尺度模型,率先采用无网格方法对液晶聚合物在非均匀剪切流场中的微观结构进行了模拟研究.无网格方法精度高、稳定性好的特性保证了模拟结果的可靠性.研究了Deborah数对平板Poiseuille流中液晶聚合物微观结构的影响,预测出非均匀剪切流场中液晶聚合物的一种单一结构和五种复合结构.指出在复合结构的过渡区,分子运动具有不稳定性,可能产生瑕疵. 关键词： 液晶 微观结构 双尺度 无网格     

微尺度型腔内液态金属流动规律模拟研究

修正了传统的Navier-Stoke方程,并利用修正后的方程与Fluent软件对微尺度型腔内Zn-Al合金微流动规律进行了数值模拟.模拟结果表明：液态金属在微尺度管道内的流动规律在管径为0.5mm时出现临界状态,管径越细,速度附面层相对厚度就越大.同时,微管流动存在着宏观流动中没有的“凸进效应”,且随着入口压力增大,流动前沿自由液面的凸进效应减弱.在流动前沿区域和上游区域之间出现负压力梯度区,且管径越细越明显. 关键词： 微管道 微流动 Zn-Al合金     

基于宏观方程数值本构关系的气体动理论加速收敛方法

在跨流域复杂流动问题的模拟中,基于求解速度分布函数演化方程的气体动理论方法的效率问题一直受到工程应用领域关注.研究提升气体动理论方法在定常流动模拟中的计算效率具有重要意义.为了提升定常流动计算收敛速度,本文提出了一种耦合宏观方程数值本构关系的气体动理论加速收敛方法.通过求解Boltzmann模型方程,将应力、热流高阶项的数值解与宏观方程耦合,实现了宏观方程的封闭;另一方面,宏观方程的计算结果被用来更新Boltzmann模型方程的当地平衡态速度分布函数中的宏观物理量,以此构造求解Boltzmann模型方程的全隐式数值格式.通过跨流域方腔流动、超声速圆柱绕流及双圆柱干扰绕流案例的数值模拟,对方法进行了广泛考核.计算结果与常规气体动理论统一算法、直接模拟蒙特卡罗法符合良好,证明该方法很好地描述了稀薄流动中的非线性本构关系,以及激波、强壁面剪切、流动分离等强非平衡特征.进一步,对于低努森数Kn的流动,方法能显著加速收敛过程,提升计算效率;随着努森数Kn增加,气体对流输运效应减弱,方法的加速收敛效果降低.与此同时,如何减少内迭代耗时,进一步提升效率有待更多研究.     

定向凝固宏微观热质传递及双扩散流作用

针对氯化铵水溶液定向凝固实验,从枝晶微尺度到宏观尺度耦合模拟了相变热质传递过程,并考察了双扩散流的影响.研究结果表明,NH4CI-H2O溶液的枝晶微界面浓度变化与非平衡态的谢尔方程趋于一致.两相区内液相率较高,易于流动的产生;热浮升力形成逆时针回流,溶质浮升力形成顺时针回流;热浮升力作用大于溶质浮升力.但双扩散流几乎不影响凝固沿着结晶室高度方向进行的一维定向特性.数值模拟中低估了两相区的渗透率,是造成两相区厚度模拟结果小于实验值的原因.     

非结构化网格下自由分子微电热推力器数值模拟

采用非结构化网格,压力边界条件的DSMC-FVM耦合方法对自由分子微电阻加热推力器微喷管内的流动与换热进行了数值模拟.分析了微喷管的流场特征,证明所采用的DSMC-FVM多尺度方法对微喷管内流场与温度场模拟是合适而且必要的.     

eXtended Pom-Pom黏弹性流体的改进光滑粒子动力学模拟

黏弹性流体广泛存在于自然界和工业生产中,对其复杂流变特性的研究具有重要的学术价值和应用意义.本文提出一种改进的光滑粒子动力学方法,对基于eXtended Pom-Pom模型的黏弹性流动进行了数值模拟.为了提高计算精度,采用一种不含核导数计算的核梯度修正离散格式.为了防止粒子穿透固壁,提出一种增强型的边界处理技术.为了消除张力不稳定性,将人工应力耦合到动量守恒方程中.运用改进光滑粒子动力学方法数值模拟了基于eXtended Pom-Pom模型的黏弹性Poiseuille流和黏弹性液滴撞击固壁问题,通过与解析解或有限差分方法解的比较以及对数值收敛性的评价,验证了改进光滑粒子动力学方法的有效性和优势,并在此基础上,深入分析了Reyonlds数、Weissenberg数、溶剂黏度比、各向异性参数、松弛时间比和分子链臂数等流变参数对流动过程的影响.     

纤维悬浮聚合物熔体描述的均一结构多尺度模型

通过对纤维悬浮聚合物熔体的可逆和不可逆热力学过程的耦合,建立了分子链弹性哑铃模型与悬浮纤维取向描述相耦合的、具有均一 (GENERIC) 结构形式的熔体多尺度模型. 由该均一结构的多尺度模型不仅可以得出熔体不同尺度上的应力贡献,还可为一般多尺度模型方程组的建立提供其结构形式均一化的方法. 关键词： GENERIC结构 纤维取向 黏弹性 聚合物     

基于总能形式的耦合的双分布函数热晶格玻尔兹曼数值方法

双分布函数热晶格玻尔兹曼数值方法在微尺度热流动系统中得到广泛的应用. 本文基于晶格玻尔兹曼平衡分布函数低阶Hermite展开式, 创新性地提出了包含黏性热耗散和压缩功的耦合的双分布函数热晶格玻尔兹曼数值方法, 将能量场内温度的变化以动量源的形式引入晶格波尔兹曼动量演化方程, 实现了能量场与动量场之间的耦合. 研究了考虑黏性热耗散和压缩功的和不考虑的两种热自然对流模型, 重点分析了不同瑞利数和普朗特数下流场内的流动情况以及温度、速度和平均努赛尔数的变化趋势. 本文实验结果与文献结果一致, 验证了本文数值方法的可行性和准确性. 研究结果表明: 随着瑞利数和普朗特数的增大, 方腔内对流传热作用逐渐增强, 边界处形成明显的边界层; 考虑黏性热耗散和压缩功的模型对流作用相对增强, 黏性热耗散和压缩功对自然对流的影响在微尺度流动过程中不能忽略.     

基于CFD/CSD的舵面控制气动弹性数值模拟方法

在动态网格上通过耦合求解流动控制方程和结构动力学方程, 发展了一种舵面控制下飞行器运动响应过程中气动弹性数值模拟研究方法.流动控制方程采用N-S方程, 结构动力学采用线性模态叠加方法, 其中流动控制方程空间离散采用基于非结构网格的有限体积方法, 对流通量采用计算HLLC格式, 非定常时间离散采用基于LU-SGS的双时间步长方法.模拟中, 气动运动和结构变形在双时间步长方法推进过程中采用改进松耦合方法, 气动网格与结构网格之间信息交换采用无限平板样条法实现, 飞行器的运动和变形采用基于重叠网格和Delaunay图映射变形网格相结合的方法进行处理.采用多个考核算例对发展的数值方法进行考核验证, 结果表明该方法可以高效精确模拟舵面开环控制下飞行器运动响应过程中的气动弹性特性.      

生物大分子多尺度理论和计算方法

分子模拟是研究生物大分子的重要手段. 过去二十年来, 人们将分子模拟与实验研究相结合, 揭示出生物大分子结构和动力学方面的诸多重要性质. 传统分子模拟主要采用全原子分子模型或各种粗粒化的分子模型. 在实际应用中, 传统分子模拟方法通常存在精度或效率瓶颈, 一定程度上限制了其应用范围. 近年来, 多尺度分子模型越来越受到人们的关注. 多尺度分子模型基于统计力学原理, 将全原子模型和粗粒化模型相耦合, 有望克服传统分子模拟方法中的精度/效率瓶颈, 进而拓展分子模拟在生物大分子研究中的应用范围. 根据模型之间的耦合方式, 近年来发展起来的多尺度分子模拟方法可归纳为如下四种类型: 混合分辨多尺度模型、并行耦合多尺度模型、单向耦合多尺度模型、以及自学习多尺度模型. 本文将对上述四类多尺度模型做简要介绍, 并讨论其主要优缺点、应用范围以及进一步发展方向.     

Time Periodic Electroosmotic Flow of The Generalized Maxwell Fluids in a Semicircular Microchannel

Analytical solutions are presented using method of separation of variables for the time periodic electroosmotic flow (EOF) of linear viscoelastic fluids in semicircular microchannel. The linear viscoelastic fluids used here are described by the general Maxwell model. The solution involves analytically solving the linearized Poisson-Boltzmann (P -B) equation, together with the Cauchy momentum equation and the general Maxwell constitutive equation. By numerical computations, the influences of electric oscillating Reynolds number Re and Deborah number De on velocity amplitude are presented. For small Re, results show that the larger velocity amplitude is confined to the region near the charged wall when De is small. With the increase of the Deborah number De, the velocity far away the charged wall becomes larger for large Deborah number De. However, for larger Re, the oscillating characteristic of the velocity amplitude occurs and becomes significant with the increase of De, especially for larger Deborah number.     

Bubble oscillation and inertial cavitation in viscoelastic fluids

Non-linear acoustic oscillations of gas bubbles immersed in viscoelastic fluids are theoretically studied. The problem is formulated by considering a constitutive equation of differential type with an interpolated time derivative. With the aid of this rheological model, fluid elasticity, shear thinning viscosity and extensional viscosity effects may be taken into account. Bubble radius evolution in time is analyzed and it is found that the amplitude of the bubble oscillations grows drastically as the Deborah number (the ratio between the relaxation time of the fluid and the characteristic time of the flow) increases, so that, even for moderate values of the external pressure amplitude, the behavior may become chaotic. The quantitative influence of the rheological fluid properties on the pressure thresholds for inertial cavitation is investigated. Pressure thresholds values in terms of the Deborah number for systems of interest in ultrasonic biomedical applications, are provided. It is found that these critical pressure amplitudes are clearly reduced as the Deborah number is increased.     

Grid and subgrid-scale interactions in viscoelastic turbulent flow and implications for modelling

Using direct numerical simulations of turbulent plane channel flow of homogeneous polymer solutions, described by the Finitely Extensible Nonlinear Elastic-Peterlin (FENE-P) rheological constitutive model, a-priori analyses of the filtered momentum and FENE-P constitutive equations are performed. The influence of the polymer additives on the subgrid-scale (SGS) energy is evaluated by comparing the Newtonian and the viscoelastic flows, and a severe suppression of SGS stresses and energy is observed in the viscoelastic flow. All the terms of the transport equation of the SGS kinetic energy for FENE-P fluids are analysed, and an approximated version of this equation for use in future large eddy simulation closures is suggested. The terms responsible for kinetic energy transfer between grid-scale (GS) and SGS energy (split into forward/backward energy transfer) are evaluated in the presence of polymers. It is observed that the probability and intensity of forward scatter events tend to decrease in the presence of polymers.     

高磁雷诺数下剪切流对双撕裂模中二级磁岛的影响

在高磁雷诺数下,当双撕裂模发展进入快速磁场重联阶段时,会发生二级磁岛不稳定性,加剧磁场能量的释放。本文基于扰动形式的守恒磁流体方程组发展高精度的数值模拟程序,在平板位形下研究反对称位形剪切流对双撕裂模中二级磁岛的影响。结果表明：随着剪切流强度和剪切梯度的增加,二级磁岛的数目以及电流片横纵比变小。此外,较强的极向剪切流能够抑制二级磁岛不稳定性的发生。     

On the effect of viscosity and thermal conductivity on sound propagation in ducts: A re-visit to the classical theory with extensions for higher order modes and presence of mean flow

The dispersion equation for the axisymmetric modes of viscothermal acoustic wave propagation in uniform hard-walled circular ducts containing a quiescent perfect gas is classical. This has been extended to cover the non-axisymmetric modes and real fluids in contemporary studies. The fundamental axisymmetric mode has been the subject of a large number of studies proposing approximate solutions and the characteristics of the propagation constants for narrow and wide ducts with or without mean flow is well understood. In contrast, there are only few publications on the higher order modes and the current knowledge about their propagation characteristics is rather poor. On the other hand, there is a void of papers in the literature on the effect of the mean flow on the quiescent modes of propagation. The present paper aims to contribute to the filling of these gaps to some extent. The classical theory is re-considered with a view to cover all modes of acoustic propagation in circular ducts carrying a real fluid moving axially with a uniform subsonic velocity. The analysis reveals a new branch of propagation constants for the axisymmetric modes, which appears to have escaped attention hitherto. The solution of the governing wave equation is expressed in a modal transfer matrix form in frequency domain and numerical results are presented to show the effects over wide ranges of frequency, viscosity and mean flow parameters on the propagation constants. The theoretical formulation allows for the duct walls to have finite impedance, but no numerical results are presented for lined ducts or ducts carrying a sheared mean flow.     

高精度有限差分在湍流直接数值模拟中的应用

本文采用高精度有限差分法对矩形管道内充分发展湍流换热进行了直接数值模拟,湍流雷诺数Rer和普朗特数分别为400(Rem=6200)和0.71,压力泊松方程分别采用两种不同的离散格式(二阶和四阶中心差分离散).结果表明,同二阶中心差分格式相比,高精度有限差分可以在较少的网格下得到较好的结果;压力泊松方程采用四阶中心差分或二阶中心差分对计算结果的影响甚小,但采用二阶中心差分时,可以节省大量的计算时间.     

Comparison of measured broadband noise attenuation spectra from circular flow ducts and from lined engine intakes with predictions

Sound propagation in lined circular ducts is investigated in the presence of uniform and sheared flow. The modal solutions are obtained by solving an eigenvalue equation which, in the case of sheared flow, is derived by using finite differences and by matching the pressure and the radial component of the particle velocity at the interface of the regions of uniform and sheared flow. For the uniform flow region, standard Bessel function solutions are used. The attenuation of acoustic energy at a given frequency and for a given liner length is computed on the assumption that at the inlet to the lined duct, the acoustic energy is equally distributed among the propagating modes. The total number of propagating modes is determined from the hard wall “cut off” condition. The failure to find some of the modal solutions on the attenuation computed in this way is discussed. It is shown that the reliability of this method of computing liner attenuation depends on the ability to successfully compute most of the modal solutions over a large range of frequencies, flow conditions and duct wall impedance values. A numerical technique is developed which uses a fraction of the total number of solutions to compute the total attenuations without appreciable loss of accuracy. Measured attenuation spectra from a flow duct facility and from lined intake ducts of the RB.211 engine are compared with predictions. In general very good agreement between predictions and measurements is obtained.     

Rigorous solutions for sound radiation from circular ducts with hyperbolic horns or infinite plane baffle

A rigorous treatment is presented of sound radiation from circular ducts with either a hyperbolic horn or an infinite plane baffle. In the analysis hyperboloidal wave functions are used, which are defined here, for the first time, as a class of eigensolutions of the wave equation for oblate spheroidal co-ordinates. The numerical results include the complex conversion (or reflection) coefficients and the radiation directivity for various incident wave modes, spinning modes as well as axisymmetric modes. The solutions are valid for the whole frequency range including frequencies above and below the cut-off frequencies of the duct modes involved.     

Multi-dimensional vortex quasi-simple waves

Multi-dimensional vortex modes of a quasi-simple wave solution is presented. These are constructed on the basis of vortex modes for ideal simple waves. A version of 2D Burgers equation is derived which is the same as that obtained for sound quasi-simple waves if neglecting the last term of the latter. Some solutions are explained in physical detail which have a localized traveling behavior. A numerical simulation is shown to support the obtained analytical solutions.