跨流域高超声速绕流环境Boltzmann模型方程统一算法研究

如何准确可靠地模拟从外层空间高稀薄流到近地面连续流的航天器高超声速绕流环境与复杂流动变化机理是流体物理的前沿基础科学问题. 基于对Boltzmann方程碰撞积分的物理分析与可计算建模, 确立了可描述自由分子流到连续流区各流域不同马赫数复杂流动输运现象统一的Boltzmann模型速度分布函数方程, 发展了适于高、低不同马赫数绕流问题的离散速度坐标法和直接求解分子速度分布函数演化更新的气体动理论数值格式, 建立了模拟复杂飞行器跨流域高超声速飞行热环境绕流问题的气体动理论统一算法. 对稀薄流到连续流不同Knudsen数0.002 ≤Kn_∞ ≤1.618、不同马赫数下可重复使用卫星体再入过程(110–70 km)中高超声速绕流问题进行算法验证分析, 计算结果与典型文献的Monte Carlo直接模拟值及相关理论分析符合得较好. 研究揭示了飞行器跨流域不同高度高超声速复杂流动机理、绕流现象与气动力/热变化规律, 提出了一个通过数值求解介观Boltzmann模型方程, 可靠模拟高稀薄自由分子流到连续流跨流域高超声速气动力/热绕流特性统一算法.     

Effects of thermal radiation on Casson fluid flow and heat transfer over an unsteady stretching surface subjected to suction/blowing

The unsteady flow of a Casson fluid and heat transfer over a stretching surface in presence of suction/blowing are investigated.The transformed equations are solved numerically by using the shooting method.The exact solution corresponding to the momentum equation for the steady case is obtained.Fluid velocity initially decreases with the increase of unsteadiness parameter.Due to an increasing Casson parameter the velocity field is suppressed.Thermal radiation enhances the effective thermal diffusivity and the temperature rises.     

Aeroacoustics of duct junction flows merging at different angles

This paper reports an exploratory study of the aeroacoustics of a merging flow at a duct junction with the same width in all branches and different merging angles. The focus is put on the acoustic generation due to the flow unsteadiness. The study is carried out by the direct aeroacoustic simulation (DAS) approach, which solves the unsteady compressible Navier–Stokes equations and the perfect gas equation of state simultaneously using the conservation element and solution element (CE/SE) method. The Mach number based on the maximum inlet velocity of side branch is 0.1 and the Reynolds number of the flow based on duct width and this velocity is 2.3×10⁵. The numerical simulations are performed in two dimensions and the aeroacoustics at different merging angles (30°, 45°, 60° and 90°) are studied. Both the levels of unsteady interactions of merging flow structures and the efficiency of the acoustic generation are observed to increase with the merging angles, where the increase in acoustic efficiency can be up to three orders of magnitude. The major acoustic source is found to be the fluctuating wall pressure induced by the flow unsteadiness in the downstream branch. A scaling law between the wall fluctuating force and the acoustic efficiency is also derived.     

Influence of radiative heat flux in nonlinear convection of quartic order in Casson fluid past a vertical permeable plate with variable suction and Hall current

The current study centralizes on unsteady free convection slip flow of Casson fluid past a vertical permeable plate with Hall current, radiative heat flux, and variable suction. The nonlinear convection is subjected to quartic order. Perturbation method is used to convert the non-linear coupled partial differential equation of the momentum and energy to a system of ordinary differential equations. The dimensionless governing equations are solved analytically for velocity and temperature profiles. The graphs are plotted for sundry parameters for variations in the distinct flow fields w.r.t distance from the plate. Variation in the skin friction for the axial and transverse cases are presented in the form of graphs for various parameters. It is observed that with the increase in the order of non-linear convection and value of radiation parameter, the velocity field increases in Casson fluid. The increase in heat absorption parameter and Prandtl number decreases the temperature profile and increase in radiative heat flux parameter increases the temperature profile.     

二维流场、热结构松耦合模拟研究

本文针对轴对称型飞行器前缘部位,采用松耦合方法,开展了二维高超声速流固界面的热耦合计算。满足非稳态N-S方程的外部流场和内部非稳态热传导均采用商用软件FLUENT进行计算。通过对相同时刻,不同松耦合推进时间步长引起的壁面热流密度、壁面温度、表面传热系数的差异对比,得出了松耦合推进时间对以上参数的影响规律,这对有效实施松耦合方法进行高超声速飞行器模拟具有积极意义。     

Influences of heat flux on extinction characteristics of steady/unsteady premixed stagnation flames

This paper investigates the extinction characteristics of premixed stagnation flames (PSFs) with controlled heat losses and flow disturbances. The low-frequency air flow pulsations that imitate the operational transients in practical combustors were specially introduced. The tunable diode-laser absorption spectroscopy (TDLAS) measurement was applied to obtain the temperature profile and wall heat flux. It is found that, for steady flame with a fixed equivalence ratio, the extinction stretch rate dramatically increases as the wall heat flux decreases. The extinction criterion is summarized as a global Karlovitz number of 0.57 by establishing a relationship between the global and local stretch rates. Numerical simulations reveal that the local extinction Karlovitz number of steady PSFs is approximately 1.0 regardless of the conditions such as heat flux and equivalence ratio. Further experiments present that the air pulsations with a repetition of ～5 Hz significantly deteriorate the flame stability. Particularly, for unsteady perturbed flames, the extinction stretch rate exhibits a nonlinear trend, yielding two regimes with discrepant sensitivities to wall heat flux. The unsteady simulation then highlights a local stretch rate overshoot in the presence of pulsation. It is caused by the time delay between the inlet velocity and flame front movement that eventually leads to poor flame stability. Moreover, in the high heat-flux regime, a smaller local stretch rate overshoot results in the weak dependence of extinction limits on heat fluxes.     

Large eddy simulation of compressible turbulent channel flow with spanwise wall oscillation

The influences of the modification of turbulent coherent structures on temperature field and heat transfer in turbulent channel flow are studied using large eddy simulation (LES) of compressible turbulent channel flows with spanwise wall oscillation (SWO). The reliability of the LES on such problems is proved by the comparisons of the drag reduction data with those of other researches. The high consistency of coherent velocity structures and temperature structures is found based on the analyses of the turbulent flow field. When the coherent velocity structures are suppressed, the transportations of momentum and heat are reduced simultaneously, demonstrating the same trend. This shows that the turbulent coherent structures have the same effects on the transportations of momentum and heat. The averaged wall heat flux can be reduced with appropriate oscillating parameters. Supported by the Key Subjects of National Natural Science Foundation of China (Grant No. 10732090), the National Natural Science Foundation of China (Grant No. 50476004), and the 111 Project (Grant No. B08009)     

METHOD AND SOME RESULTS ON EXPERIMENTAL STUDY OF UNSTEADY TURBULENT HEAT TRANSFER

Abstract

An experimental setup is described for studying unsteady heat transfer where both the heat generation rate in the channel wall and the fluid flow rate vary in time. For measuring the surface temperature and the wall heat flux, the temperature field in the thermal sublayer is used. The results of the first experiments carried out for a sudden change in channel wall heat generation are discussed. Unsteadiness rates much greater than those in previous work known to the authors were achieved     

Variable fluid properties and variable heat flux effects on the flow and heat transfer in a non-Newtonian Maxwell fluid over an unsteady stretching sheet with slip velocity

The effects of variable fluid properties and variable heat flux on the flow and heat transfer of a non-Newtonian Maxwell fluid over an unsteady stretching sheet in the presence of slip velocity have been studied. The governing differential equations are transformed into a set of coupled non-linear ordinary differential equations and then solved with a numerical technique using appropriate boundary conditions for various physical parameters. The numerical solution for the governing non-linear boundary value problem is based on applying the fourth-order Runge-Kutta method coupled with the shooting technique over the entire range of physical parameters. The effects of various parameters like the viscosity parameter, thermal conductivity parameter, unsteadiness parameter, slip velocity parameter, the Deborah number, and the Prandtl number on the flow and temperature profiles as well as on the local skin-friction coefficient and the local Nusselt number are presented and discussed. Comparison of numerical results is made with the earlier published results under limiting cases.     

高超声速飞行器前缘流固耦合计算方法研究

对高超声速流场和结构温度场进行了耦合计算分析, 同时基于准静态假设对结构应力进行了分析. 流场部分采用基于非定常Navier-Stokes (N-S)方程的有限体积法, 湍流模型采用SST k-ω 模型, 固体部分采用基于非稳态热传导方程的有限元法, 同时基于准静态假设对固体结构的应力应变进行了分析. 在流固交界面处, 高速流体从固体结构得到温度边界条件, 固体结构从高速流体得到热流边界条件, 从而实现了流场和固体温度场的紧耦合计算.通过与超声速无限长圆管绕流试验结果进行对比, 验证了该方法的可靠性. 同时对二维圆管结构在气动加热过程中的温度、应力等的变化进行了比较详细的分析. 研究结果表明: 随着气动加热时间的推进, 由于圆管结构的高温区在不断扩大, 导致了结构的热变形在不断地增大; 圆管最小变形区出现在θ为60°处; 同时研究发现在计算时间内圆管热变形对外部流场的影响可以忽略不计.     

Numerical Analysis of Unsteady Magneto-Biphase Williamson Fluid Flow with Time Dependent Magnetic Field

Numerical investigation of the dusty Williamson fluid with the dependency of time has been done in current disquisition. The flow of multiphase liquid/particle suspension saturating the medium is caused by stretching of porous surface. The influence of magnetic field and heat generation/absorption is observed. It is assumed that particle has a spherical shape and distributed uniformly in fluid matrix. The unsteady two-dimensional problems are modeled for both fluid and particle phase using conservation of mass, momentum and heat transfer. The finalized model generates the non-dimensioned parameters, namely Weissenberg number, unsteadiness parameter, magnetic parameter,heat generation/absorption parameter, Prandtl number, fluid particle interaction parameter, and mass concentration parameters. The numerical solution is obtained. Locality of skin friction and Nusselt number is deliberately focused to help of tables and graphs. While inferencing the current article it is clearly observed that increment of Williamson parameter, unsteadiness parameter, magnetic parameter, volume fraction parameter, and mass concentration parameter reduces the velocity profile of fluid and solid particles as well. And increment of Prandtl number, unsteadiness parameter,volume fraction parameter, and mass concentration parameter reduces the temperature profile of fluid and solid particles as well.     

微通道内给定壁面热流的气体换热研究

本文利用实施给定热流边界条件的DSMC方法,对短通道内给定壁面热流边界条件下的气体换热情况进行了模拟.结果表明,壁面热流密度增大导致通道内压力分布非线性程度增加.随着热流密度的增大,截面速度分布趋于平缓,滑移速度增大.给定热流密度的通道壁面温度与气流截面平均温度的差值沿程增大,温度梯度沿程下降,气体稀薄性增大时,通道换热减弱.     

Advanced Study of Unsteady Heat and Chemical Reaction with Ramped Wall and Slip Effect on a Viscous Fluid

This paper investigate the effect of slip boundary condition, thermal radiation, heat source, Dufour number,chemical reaction and viscous dissipation on heat and mass transfer of unsteady free convective MHD flow of a viscous fluid past through a vertical plate embedded in a porous media. Numerical results are obtained for solving the nonlinear governing momentum, energy and concentration equations with slip boundary condition, ramped wall temperature and ramped wall concentration on the surface of the vertical plate. The influence of emerging parameters on velocity,temperature and concentration fields are shown graphically.     

Investigation of the benefits of unsteady blowing actuation on a 2D wind turbine blade

This paper investigates the benefit of unsteady blowing actuation over a two-dimensional (2D) airfoil specially designed for wind turbine applications. The experiments were carried out in Syracuse University’s anechoic wind tunnel, both with and without large-scale unsteadiness in the free stream generated by a 2D cylinder upstream of the airfoil. By analyzing both surface pressure through wavelet analysis and Particle Image Velocimetry (PIV) velocity field measurements, we found a drastic change in the flow physics and the aerodynamic loading on the airfoil between steady and unsteady free-stream conditions. When there was no large-scale unsteadiness introduced in the flow, under open-loop flow control conditions with unsteady blowing, the leading-edge separation was delayed and the maximum lift coefficient was increased. For the cases where large-scale unsteadiness was introduced into the flow, the experiments showed that both open-loop and closed-loop control cases were capable of reducing load fluctuations by a measurable amount. However, only the closed-loop control case that utilized dynamic surface pressure information from the airfoil suction side near the leading edge was capable of consistently mitigating the fluctuating load.     

Handling of temperature dependence of viscosity in problems of incompressible medium flow around a cylinder

The viscous incompressible medium (water, air) flow past a circular cylinder is considered with regard for the temperature T dependent viscosity v. The influence of different boundary conditions for temperature on flow structure, the drag coefficient and its components due to the pressure and viscosity is investigated in the problem of the flow past a cylinder at rest for the (diameter-based) Reynolds number Re_D = 40. A relation between the viscosity gradient along a normal to the body surface and the integral vorticity flux from the body surface into the boundary layer is discussed. Unlike the constant viscosity case the vorticity flux may be different from zero, which must lead because of the integral conservation law for the vorticity to an alteration of the far-field boundary conditions for the velocity. In the same connection, the problem is analysed on the heat spot entry into the computational region under consideration for the flow past a circular cylinder. The examples of the symmetrization of separated flow past a cylinder performing rotation oscillations in a uniform free stream (the Taneda problem) are considered. A comparison with flow computations for low Mach numbers M « 1 for the flow of a medium past a cylinder at rest is carried out. At the computation of the equation for heat transfer under the assumption of incompressibility of such media as air, it is proposed to retain the pressure derivative, which is typical of gases. In this case, a better agreement with the computations of compressible flows (for M « 1) is achieved, for example, at the determination of the sizes of a symmetric zone of flow separation past a circular cylinder. An unsteady flow in the neighborhood of the point of joining the zero streamline bounding a closed region of separated flow (the cavity) in a wake of the cylinder at rest is obtained by a numerical simulation at the Reynolds number equal to 40.     

用格子Boltzmann模型模拟非等温流场

根据微观和宏观之间的质量、动量、能量守恒准则和在原格子Boltzmann模型基础上,建立了几个新的格子Boltzmann模型,使得在外力场中的格子Boltzmann模型得到进一步完善.通过还原宏观流体力学方程,捕捉到了浮力强迫系数与Grashof数之间的关系.所得动量方程和Navier Stokes方程相比,在黏性输运项上有明显的改进,说明黏性应力不但与流体的速度梯度和流体的压缩性有关,而且还与非定常的内能梯度和动量通量有关.该模型对非等温流场的数值结果证明了其具有很好的数值稳定性和适用性. 关键词： Boltzmann模型 平衡分布函数 流体力学方程     

圆管层流脉冲流动对流换热数值分析

对等热流和等壁温边界条件下圆管内层流脉冲流动对流换热问题进行了数值模拟。在等热流边界条件下的数值计算结果与理论解吻合很好。计算结果表明:在等热流和等壁温边界下脉冲流动可引起速度、温度以及努塞尔数随时间波动,振幅越大,脉冲频率越小,波动越大。但它们的时均值均等于在相同雷诺数下稳态流动的值,脉冲流动不能强化换热。     

Serrin-Type Blowup Criterion for Viscous, Compressible, and Heat Conducting Navier-Stokes and Magnetohydrodynamic Flows

This paper establishes a blowup criterion for the three-dimensional viscous, compressible, and heat conducting magnetohydrodynamic (MHD) flows. It is essentially shown that for the Cauchy problem and the initial-boundary-value one of the three-dimensional compressible MHD flows with initial density allowed to vanish, the strong or smooth solution exists globally if the density is bounded from above and the velocity satisfies Serrin’s condition. Therefore, if the Serrin norm of the velocity remains bounded, it is not possible for other kinds of singularities (such as vacuum states vanishing or vacuum appearing in the non-vacuum region or even milder singularities) to form before the density becomes unbounded. This criterion is analogous to the well-known Serrin’s blowup criterion for the three-dimensional incompressible Navier-Stokes equations, in particular, it is independent of the temperature and magnetic field and is just the same as that of the barotropic compressible Navier-Stokes equations. As a direct application, it is shown that the same result also holds for the strong or smooth solutions to the three-dimensional full compressible Navier-Stokes system describing the motion of a viscous, compressible, and heat conducting fluid.     

General Solution for Unsteady Natural Convection Flow with Heat and Mass in the Presence of Wall Slip and Ramped Wall Temperature

This work is focused on the effect of heat and mass transfer with unsteady natural convection flow of viscous fluid along with ramped wall temperature under the assumption of the slip wall condition at the boundary. Analytical solutions are obtained by using Laplace transformation to the non-dimensional set of governing equations containing velocity, temperature and concentration. Moreover, the expression for skin-friction is derived by differentiating the analytical solutions of fluid velocity. Numerical tables for Skin-friction, Sherwood number and Nusselt-number are examined. For the physical aspects of the flow, we use various values of involved physical parameters such as Prandtl number (Pr), slip parameter ($\eta$), Schmidt number (Sc), buoyancy ratio parameter ($N$), Sherwood number (Sh), and time $(t)$. Additionally, the general solutions are plotted graphically and a comprehensive theoretical section of numerical discussions is included.