竖壁传热对方腔内Benard对流的影响

本文采用具有QUICK差分格式的SIMPLE算法对边界竖壁传热的方腔内空气Benard对流进行了数值计算,根据计算结果探讨了竖壁传热对Benard对流的影响.计算表明,在所考虑的几何和物理条件下,所有竖壁绝热时,腔内流体形成平行与短轴的多个涡卷;竖壁存在传热时,腔内流体形成平行于长轴的两个涡卷,并且平行于长轴竖壁的传热热流方向相反时,涡卷的旋转方向也相反;垂直于长轴蛏壁的传热对近壁附近的流动有一定影响.     

Time-dependent natural convection in a square cavity: Application of a new finite volume method

A new finite volume (FV) approach with adaptive upwind convection is used to predict the two-dimensional unsteady flow in a square cavity. The fluid is air and natural convection is induced by differentially heated vertical walls. The formulation is made in terms of the vorticity and the integral velocity (induction) law. Biquadratic interpolation formulae are used to approximate the temperature and vorticity fields over the finite volumes, to which the conservation laws are applied in integral form. Image vorticity is used to enforce the zero-penetration condition at the cavity walls. Unsteady predictions are carried sufficiently forward in time to reach a steady state. Results are presented for a Prandtl number (Pr) of 0-71 and Rayleigh numbers equal to 10³, 10⁴ and 10⁵. Both 11 × 11 and 21 × 21 meshes are used. The steady state predictions are compared with published results obtained using a finite difference (FD) scheme for the same values of Pr and Ra and the same meshes, as well as a numerical bench-mark solution. For the most part the FV predictions are closer to the bench-mark solution than are the FD predictions.     

Aspect-ratio dependent oscillatory thermocapillary convection in the floating half zone

The dependency of the critical Marangoni number on the geometrical aspect ratio of the floating half zone is essential to predict the onset of oscillatory thermocapillary convection.The experimental studies in the microgravity conditions on floating half zones of several centimeters in diameter have predicted that the critical Marangoni number increases with the increasing aspect ratio,and the terrestrial experimental studies have predicted the contradictory conclusion for floating half zones of several mil...     

Consistency of iterative operator‐splitting methods: Theory and applications

In this article, we describe a different operator‐splitting method for decoupling complex equations with multidimensional and multiphysical processes for applications for porous media and phase‐transitions. We introduce different operator‐splitting methods with respect to their usability and applicability in computer codes. The error‐analysis for the iterative operator‐splitting methods is discussed. Numerical examples are presented. © 2009 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2010     

流体饱和多孔介质中热发展强迫对流的熵产分析

利用Darcy流动模型,同时考虑粘性耗散效应,分析研究了以两等温平板为边界的多孔介质中,热发展强迫对流的熵产.参数研究表明,组参数和Péclet数减小,而Brinkman数增大时,熵产增大.形象化的热线表示方法着重应用于Br＜0的情况,在这种情况下,在顺流的某些位置上,存在热传递方向的改变,即从原来的壁面向外传热变为向壁面传热.     

A characteristic nonoverlapping domain decomposition method for multidimensional convection‐diffusion equations

We develop a quasi‐two‐level, coarse‐mesh‐free characteristic nonoverlapping domain decomposition method for unsteady‐state convection‐diffusion partial differential equations in multidimensional spaces. The development of the domain decomposition method is carried out by utilizing an additive Schwarz domain decomposition preconditioner, by using an Eulerian‐Lagrangian method for convection‐diffusion equations and by delicately choosing appropriate interface conditions that fully respect and utilize the hyperbolic nature of the governing equations. Numerical experiments are presented to illustrate the method. © 2004 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2005     

Solution to transient Navier–Stokes equations by the coupling of differential quadrature time integration scheme with dual reciprocity boundary element method

The two‐dimensional time‐dependent Navier–Stokes equations in terms of the vorticity and the stream function are solved numerically by using the coupling of the dual reciprocity boundary element method (DRBEM) in space with the differential quadrature method (DQM) in time. In DRBEM application, the convective and the time derivative terms in the vorticity transport equation are considered as the nonhomogeneity in the equation and are approximated by radial basis functions. The solution to the Poisson equation, which links stream function and vorticity with an initial vorticity guess, produces velocity components in turn for the solution to vorticity transport equation. The DRBEM formulation of the vorticity transport equation results in an initial value problem represented by a system of first‐order ordinary differential equations in time. When the DQM discretizes this system in time direction, we obtain a system of linear algebraic equations, which gives the solution vector for vorticity at any required time level. The procedure outlined here is also applied to solve the problem of two‐dimensional natural convection in a cavity by utilizing an iteration among the stream function, the vorticity transport and the energy equations as well. The test problems include two‐dimensional flow in a cavity when a force is present, the lid‐driven cavity and the natural convection in a square cavity. The numerical results are visualized in terms of stream function, vorticity and temperature contours for several values of Reynolds (Re) and Rayleigh (Ra) numbers. Copyright © 2008 John Wiley & Sons, Ltd.     

A convective weakly viscoelastic rotating flow with pressure Neumann condition

The objective of this work is to investigate through the numeric simulation, the effects of the weakly viscoelastic flow within a rotating rectangular duct subject to a buoyancy force due to the heating of one of the walls of the duct. A direct velocity–pressure algorithm in primitive variables with a Neumann condition for the pressure is employed. The spatial discretization is made with finite central differences on a staggered grid. The pressure field is directly updated without any iteration. Numerical simulations were done for several Weissemberg numbers (We) and Grashof numbers (Gr) . The numerical results show that for high Weissemberg numbers (We>7.4 × 10^?5) and for ducts with aspect ratio 2:1 and 8:1, the secondary flow is restabilized with a stretched double vortex configuration. It is also observed that when the Grashof number is increased (Gr>17 × 10^?4) , the buoyancy force neutralizes the effects of the Coriolis force for ducts with aspect ratio 8:1. Copyright © 2008 John Wiley & Sons, Ltd.     

Topology optimization of microfluidic mixers

This paper demonstrates the application of the topology optimization method as a general and systematic approach for microfluidic mixer design. The mixing process is modeled as convection dominated transport in low Reynolds number incompressible flow. The mixer performance is maximized by altering the layout of flow/non‐flow regions subject to a constraint on the pressure drop between inlet and outlet. For a square cross‐sectioned pipe the mixing is increased by 70% compared with a straight pipe at the cost of a 2.5 fold increase in pressure drop. Another example where only the bottom profile of the channel is a design domain results in intricate herring bone patterns that confirm findings from the literature. Copyright © 2008 John Wiley & Sons, Ltd.     

超临界CO2-DME混合物竖直圆管内传热特性数值研究

本文采用RNG k-ε湍流模型对超临界CO₂/DME(二甲醚)二元混合工质在竖直圆管内的传热特性进行了数值模拟研究。管径4 mm,管长为1000 mm;CO₂/DME浓度配比分别为97/3、95/5、92/8、90/10、85/15、以及70/30;质量流速为125~200 kg·m^-2.s^-1;热流密度为15~30 kW.m^-2,入口温度295~308 K,入口压力8~15 MPa。不同浓度配比的混合工质在各自临界压力下应用时,随着DME浓度的增加,换热系数的峰值逐渐减低,但在温度大于310 K时混合工质的换热系数会高于纯CO₂。压力相同时,随着DME浓度的增大,拟临界温度升高,换热系数峰值点也随之向温度升高的方向移动。混合工质的换热系数随质量流速的增大而增大。在拟临界点前,增大热流密度及降低压力对管内传热有利,而在拟临界点之后,换热系数随热流密度的升高以及压力的降低而降低。