Stability analysis of a Maxwell fluid in a porous medium heated from below

Based on a modified-Darcy–Brinkman–Maxwell model, stability analysis of a horizontal layer of Maxwell fluid in a porous medium heated from below is performed. By solving the eigenvalue problems, the critical Rayleigh number, wave number and frequency for overstability are determined. It is found that the critical Rayleigh number for overstability decreases as the relaxation time increases and the elasticity of a Maxwell fluid has a destabilizing effect on the fluid layer in porous media. On the other hand, the critical Rayleigh number for overstability increases by increasing the porous parameter which acts to stabilize the system. In limiting cases, some previous results for viscoelastic fluids in nonporous media are recovered from our results.     

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     

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浓度的增大,拟临界温度升高,换热系数峰值点也随之向温度升高的方向移动。混合工质的换热系数随质量流速的增大而增大。在拟临界点前,增大热流密度及降低压力对管内传热有利,而在拟临界点之后,换热系数随热流密度的升高以及压力的降低而降低。     

Comparison of multicomponent fuel droplet vaporization experiments in forced convection with the Sirignano model

The vaporization of multicomponent fuel droplets was studied experimentally in a heated flow and the results were compared to the model proposed by Abramzon and Sirignano. The droplet was suspended on a permanent holder which was set up in a thermal wind-tunnel. This wind-tunnel was fitted with a video recording system and an infra-red camera. The period during which the droplet was suspended on the holder before the opening of the hot air flow damper was recorded. This first sequence corresponds to the droplet vaporization in natural convection, whose initial experiment conditions, especially diameter, temperature, composition of the droplet, are well known. Then the damper was turn on, and the sequence of forced convection begun. The initial diameter of the droplet was recorded by the video system. The other initial conditions of this second sequence cannot be determined experimentally. The distribution of temperature in the droplet and the surface temperature, the mass fraction distribution in the droplet and the surface mass fraction were unknown. These unknown parameters were determined by coupling our experiment with a model using “the film concept” in natural convection. Experimental results were compared with the calculations and found satisfactory, in natural convection as well as in forced convection initiated by this method. The method was tested in the case of a fuel mixture droplets (heptane–decane) for different initial concentrations and variable durations of the sequence in natural convection.