磁约束核聚变关键能量转换部件的磁流体力学研究

简述了磁约束热核聚变反应堆关键能量转换部件包层相关的流体力学及传热学研究背景, 重点介绍了强磁场作用下的金属流体流动与传热问题的实验及数值研究进展.     

基于蒙特卡罗随机有限元方法的随机多孔介质内流体自然对流不确定性研究

为分析孔隙率不确定性对多孔介质方腔内自然对流换热的影响,发展了一种基于KL（Karhunen-Loeve展开）-蒙特卡罗随机有限元算法的随机多孔介质内自然对流不确定性分析数理模型及有限元数值模拟程序框架。通过K-L展开及基于拉丁抽样法生成多孔介质孔隙率随机实现,并耦合多孔介质自然对流有限元程序,进行随机多孔介质内自然对流传热数值模拟,得出了多孔介质内流场与温度场平均值与标准偏差,并分析了孔隙率不确定性条件下Da数对Nu数的影响。结果表明,孔隙率不确定性对多孔介质方腔内自然对流有重要影响。随机多孔介质内流场及温度场与确定性条件下的流场及温度场存在一定偏差,Nu数标准偏差随着Da的增大先增大后减小。     

NUMERICAL INVESTIGATION ON FLOW STABILITY OF RAYLEIGH-B?NARD CONVECTION OF COLD WATER IN A RECTANGULAR CAVITY COOLED AND HEATED SYMMETRICALLY RELATIVE TO THE TEMPERATURE OF DENSITY MAXIMUM

为了解具有密度极值流体瑞利-贝纳德对流特有现象和规律,利用有限容积法对长方体腔内关于密度极值温度对称加热-冷却时冷水瑞利-贝纳德对流的分岔特性进行了三维数值模拟,得到了不同条件下的对流结构型态及其分岔序列,分析了密度极值特性、瑞利数、热边界条件以及宽深比对瑞利-贝纳德对流的影响. 结果表明:具有密度极值冷水瑞利-贝纳德对流系统较常规流体更加稳定,且流动型态及其分岔序列更加复杂;相同瑞利数下多种流型可以稳定共存,各流型在相互转变中存在滞后现象;随着宽深比的增加,流动更易失稳,对流传热能力增强;系统在导热侧壁时比绝热侧壁更加稳定,对流传热能力有所减弱;基于计算结果,采用线性回归方法,得到了热壁传热关联式.     

Magnetic Heat Transfer Enhancements on Fin-Tube Heat Exchangers

通过DNS方法解耦合的三维非稳态流动和固流体能量方程组,本文研究了两平行磁质平板和圆管所组成的肋片式圆管换热器单元与震荡流体间的传热过程.对不同的磁场频率和振幅的三维动态流热场的模拟结果表明增强磁场频率和振幅能很有效地增加周期平均传热强度达到强化传热的目的.     

半浮区热毛细对流的不稳定性与转捩

近二十年来,微重力流体开展了半浮区液桥热毛细对流的不稳定性与转捩的研究．文中给出了热毛细振荡对流发生的临界参数,分析了液桥几何位形（尺度比,体积比）、物理参数及传热参数对临界Ｍａｘａｎｇｏｎｉ的影响．报导了有关的地面模拟实验,微重力实验以及本问题的线性稳定性分析、能量分析和数值模拟结果,并介绍了定常轴对称热毛细对流通过非定常振荡热毛细对流到湍流的转捩过程和三种热毛细振荡对流的产生机理．     

半浮区热毛细对流的不稳定性与转捩

近二十年来,微重力流体开展了半浮区液桥热毛细对流的不稳定性与转捩的研究．文中给出了热毛细振荡对流发生的临界参数,分析了液桥几何位形（尺度比,体积比）、物理参数及传热参数对临界Ｍａｘａｎｇｏｎｉ的影响．报导了有关的地面模拟实验,微重力实验以及本问题的线性稳定性分析、能量分析和数值模拟结果,并介绍了定常轴对称热毛细对流通过非定常振荡热毛细对流到湍流的转捩过程和三种热毛细振荡对流的产生机理．      

多孔介质内黏弹性流体的热对流稳定性研究

基于修正的Darcy模型, 介绍了多孔介质内黏弹性流体热对流稳定性研究的现状和主要进展. 通过线性稳定性理论, 分析计算多孔介质几何形状(水平多孔介质层、多孔圆柱以及多孔方腔)、热边界条件(底部等温加热、底部等热流加热、底部对流换热以及顶部自由开口边界)、黏弹性流体的流动模型(Darcy-Jeffrey, Darcy-Brinkman-Oldroyd以及Darcy-Brinkman -Maxwell模型)、局部热非平衡效应以及旋转效应对黏弹性流体热对流失稳的临界Rayleigh数的影响. 利用弱非线性分析方法, 揭示失稳临界点附近热对流流动的分叉情况, 以及失稳临界点附近黏弹性流体换热Nusselt数的解析表达式. 采用数值模拟方法, 研究高Rayleigh数下黏弹性流体换热Nusselt数和流场的演化规律,分析各参数对黏弹性流体热对流失稳和对流换热速率的影响.主要结果:(1)流体的黏弹性能够促进振荡对流的发生;(2)旋转效应、流体与多孔介质间的传热能够抑制黏弹性流体的热对流失稳;(3)在临界Rayleigh数附近,静态对流分叉解是超临界稳定的, 而振荡对流分叉可能是超临界或者亚临界的,主要取决于流体的黏弹性参数、Prandtl数以及Darcy数;(4)随着Rayleigh数的增加,热对流的流场从单个涡胞逐渐演化为多个不规则单元涡胞, 最后发展为混沌状态.      

浮力对混合对流流动及换热特性的影响

用热线和冷线相结合的技术测量垂直圆管内逆混合对流流体的平均速度、 温度以及它们的脉动. 较详细地研究了浮力对逆混合对流的流动特性和传热特性的影响. 评 估了实验中采用的冷线测量温度补偿速度探头温度敏感的影响. 逆混合对流的传热结果用无 量纲参数Ω (Ω= Grd / Red2 )来表示，其中，基于管道直 径的雷诺数Red变化范围为900~18000, 浮力参数Ω变化范围为 0.004899~0.5047. 研究结果表明，浮力对逆混合对流的换热有强化作用. 随着葛拉晓夫数Grd的增加，温度脉动，流向雷诺正应力和流向温度通量增 大，并且在靠近壁面的流体区域尤其明显. 热线与冷线相结合的技术适合于研究非绝热的流 动测量，可以用于研究浮力对流动和换热特性的影响.     

多孔介质内黏弹性流体的热对流稳定性研究

基于修正的Darcy模型,介绍了多孔介质内黏弹性流体热对流稳定性研究的现状和主要进展.通过线性稳定性理论,分析计算多孔介质几何形状(水平多孔介质层、多孔圆柱以及多孔方腔)、热边界条件(底部等温加热、底部等热流加热、底部对流换热以及顶部自由开口边界)、黏弹性流体的流动模型(Darcy-Jeffrey, DarcyBrinkman-Oldroyd以及Darcy-Brinkman-Maxwell模型)、局部热非平衡效应以及旋转效应对黏弹性流体热对流失稳的临界Rayleigh数的影响.利用弱非线性分析方法,揭示失稳临界点附近热对流流动的分叉情况,以及失稳临界点附近黏弹性流体换热Nusselt数的解析表达式.采用数值模拟方法,研究高Rayleigh数下黏弹性流体换热Nusselt数和流场的演化规律,分析各参数对黏弹性流体热对流失稳和对流换热速率的影响.主要结果:(1)流体的黏弹性能够促进振荡对流的发生;(2)旋转效应、流体与多孔介质间的传热能够抑制黏弹性流体的热对流失稳;(3)在临界Rayleigh数附近,静态对流分叉解是超临界稳定的,而振荡对流分叉可能是超临界或者亚临界的,主要取决于流体的黏弹性参数、Prandtl数以及Darcy数;(4)随着Rayleigh数的增加,热对流的流场从单个涡胞逐渐演化为多个不规则单元涡胞,最后发展为混沌状态.     

基于蒙特卡罗随机有限元方法的随机多孔介质内流体自然对流不确定性研究

为分析孔隙率不确定性对多孔介质方腔内自然对流换热的影响,发展了一种基于KL(Karhunen-Loeve展开)-蒙特卡罗随机有限元算法的随机多孔介质内自然对流不确定性分析数理模型及有限元数值模拟程序框架。通过K-L展开及基于拉丁抽样法生成多孔介质孔隙率随机实现,并耦合多孔介质自然对流有限元程序,进行随机多孔介质内自然对流传热数值模拟,得出了多孔介质内流场与温度场平均值与标准偏差,并分析了孔隙率不确定性条件下Da数对Nu数的影响。结果表明,孔隙率不确定性对多孔介质方腔内自然对流有重要影响。随机多孔介质内流场及温度场与确定性条件下的流场及温度场存在一定偏差,Nu数标准偏差随着Da的增大先增大后减小。     

基于概率配点法的岩土材料参数随机场及其响应分析

概率配点法是进行不确定性问题分析的一种有效方法。通过对输入参数场进行Karhunen-Loeve展开,将其表示为一系列独立随机变量在不同权重下的线性组合,再以与之相同的随机变量组合形成混沌多项式展开对输出随机场进行分解,通过某种算法选取随机变量的值,将其作为插值点的组合(配点),在这些配点上,概率方程演化为一个确定性问题方程。由此,可以直接利用现有软件或者确定性问题计算程序进行求解,生成混沌多项式的系数矩阵后,即可得到该随机问题的各阶统计矩,从而实现参数随机场的不确定性分析。本文将该方法引进岩土工程材料参数随机场,将体积模量视为输入随机场,位移视为输出场,分别进行了弹性及塑性变形计算。结果表明该方法极大地降低了随机问题的求解难度,与MC法(Mento Carlo)相比,减少了运算消耗,提高了计算效率,具有明显的优越性。     

Efficient stochastic FEM for flow in heterogeneous porous media. Part 1: random Gaussian conductivity coefficients

This paper is concerned with the development of efficient iterative methods for solving the linear system of equations arising from stochastic FEMs for single‐phase fluid flow in porous media. It is assumed that the conductivity coefficient varies randomly in space according to some given correlation function and is approximated using a truncated Karhunen–Loève expansion. Distinct discretizations of the deterministic and stochastic spaces are required for implementations of the stochastic FEM. In this paper, the deterministic space is discretized using classical finite elements and the stochastic space using a polynomial chaos expansion. The highly structured linear systems which result from this discretization mean that Krylov subspace iterative solvers are extremely effective. The performance of a range of preconditioned iterative methods is investigated and evaluated in terms of robustness with respect to mesh size and variability of the conductivity coefficient. An efficient symmetric block Gauss–Seidel preconditioner is proposed for problems in which the conductivity coefficient has a large standard deviation.The companion paper, herein, referred to as Part 2, considers the situation in which Gaussian random fields are transformed into lognormal ones by projecting the truncated Karhunen–Loève expansion onto a polynomial chaos basis. This results in a stochastic nonlinear problem because the random fields are represented using polynomial chaos containing terms that are generally nonlinear in the random variables. Copyright © 2013 John Wiley & Sons, Ltd.     

A computational strategy for the random response of assemblies of structures

This paper presents a dedicated approach to the calculation of the random response of assemblies with uncertain interface characteristics. The random response is constructed using a polynomial chaos expansion (PCE). A decomposition of the assemblies into substructures and interfaces is defined and associated with a dedicated computational strategy which leads to a local/global algorithm enabling the treatments of the substructure and of the interface problems to be uncoupled. Since the only uncertain parameters are those which appear in the interface equations, this approach results in a drastic reduction of the computational costs. This paper first presents the classical stochastic finite element strategy for this kind of problem, then details the proposed dedicated approach. The applications concern structures assembled with uncertain elastic bonded joints. The proposed approach is compared to the Monte Carlo method and to the stochastic finite element method.     

On the propagation of statistical model parameter uncertainty in CFD calculations

This work considers a new class of finite-volume approximations for scalar and system nonlinear conservation laws with multiple sources of stochastic model parameter uncertainty. The deterministic propagation of model parameter uncertainty is achieved through the introduction of additional stochastic coordinates. Particular attention is given to the construction of specialized piecewise polynomial approximation spaces well suited to the high-order accurate approximation of solution discontinuities in both physical and stochastic dimensions without exhibiting Gibbs-like oscillations characteristic of polynomial approximation. The proposed discretization easily retrofits existing finite-volume CFD codes in use today. Numerical results are presented for inviscid Burgers equation with uncertain initial data as well as the compressible Reynolds-averaged Navier–Stokes equations with uncertain boundary data and turbulence model parameters.     

Bifurcation and Chaos Analysis of Stochastic Duffing System Under Harmonic Excitations

The Chebyshev polynomial approximation is applied to the dynamic response problem of a stochastic Duffing system with bounded random parameters, subject to harmonic excitations. The stochastic Duffing system is first reduced into an equivalent deterministic non-linear one for substitution. Then basic non-linear phenomena, such as stochastic saddle-node bifurcation, stochastic symmetry-breaking bifurcation, stochastic period-doubling bifurcation, coexistence of different kinds of steady-state stochastic responses, and stochastic chaos, are studied by numerical simulations. The main feature of stochastic chaos is explored. The suggested method provides a new approach to stochastic dynamic response problems of some dissipative stochastic systems with polynomial non-linearity.     

Efficient stochastic finite element methods for flow in heterogeneous porous media. Part 2: Random lognormal permeability

Efficient and robust iterative methods are developed for solving the linear systems of equations arising from stochastic finite element methods for single phase fluid flow in porous media. Permeability is assumed to vary randomly in space according to some given correlation function. In the companion paper, herein referred to as Part 1, permeability was approximated using a truncated Karhunen-Loève expansion (KLE). The stochastic variability of permeability is modeled using lognormal random fields and the truncated KLE is projected onto a polynomial chaos basis. This results in a stochastic nonlinear problem since the random fields are represented using polynomial chaos containing terms that are generally nonlinear in the random variables. Symmetric block Gauss-Seidel used as a preconditioner for CG is shown to be efficient and robust for stochastic finite element method.     

混合不确定度量化方法及其在计算流体动力学迎风格式中的应用

针对流体力学数值求解间断问题时,初始状态含有偶然和认知混合型的不确定性,将认知不确定度作为外层,偶然不确定度作为内层,分别使用非嵌入多项式混沌方法(non-intrusive polynomial chaos, NIPC)和概率盒(P-box)理论处理偶然不确定度和认知不确定度,发展了流体力学数值求解过程中,初始状态含有混合不确定度传播量化的一种方法。以迎风格式和黎曼解法器求解Sod问题为例,评估了左状态密度(偶然不确定度)和理想气体多方指数(认知不确定度)对模型输出结果的影响,有效验证了该方法的可行性。     

不确定结构时域响应分析的多项式维数分解法

针对具有不确定参数结构,提出时域不确定性传播和量化的多项式维数分解法,确定了结构响应统计量的演变过程.首先,采用参数概率模型来描述结构参数的不确定性,建立结构动力学方程,将结构响应表达为不确定参数的函数;进一步,将所关心的结构响应采用成员函数进行维数分解,并利用正交多项式基底对成员函数进行Fourier展开;最后,应用降维积分方法进行展开系数的求解,给出了响应均值和标准差的计算表达式.在数值算例中,将本文方法与蒙特卡洛方法进行对比,结果表明所建立方法具有较高的求解效率和计算精度.     

On the uncertainty quantification of the unsteady aerodynamics of 2D free falling plates

The aim of this paper is to conduct a statistical analysis of the effects of the fillet radii on the dynamics of the falling plate using the nonintrusive spectral projection (NISP) method. The free fall of two-dimensional cards immersed in a fluid was studied using a deterministic and stochastic numerical approach. The motion is characterized by the fluid-body interaction described by coupling the Navier–Stokes and rigid body dynamic equations. The model’s predictions have been validated using both experimental and numerical data available in the literature. In the stochastic simulations, the fillet radius of the plate was considered a random variable characterized by a uniform probability density function introducing, in this way, some uncertainties in the plate’s trajectory. To take into account the uncertainties, we employed the NISP method based on polynomial chaos expansion. The analysis was focused on finding the ensemble mean trajectory and error bar for a confidence interval of 95 % for both tumbling and fluttering regimes.