Complex Transition of Double-Diffusive Convection in a Rectangular Enclosure with Height-to-Length Ratio Equal to 4: Part I

This is the first part of direct numerical simulation (DNS) of double-diffusive convection in a slim rectangular enclosure with horizontal temperature and concentration gradients. We consider the case with the thermal Rayleigh number of 10^5, the Pradtle number of 1, the Lewis number of 2, the buoyancy ratio of composition to temperature being in the range of [0,1], and height-to-width aspect ration of 4. A new 7th order upwind compact scheme was developed for approximation of convective terms, and a three-stage third-order Runge-Kutta method was employed for time advancement. Our DNS suggests that with the buoyancy ratio increasing form 0 to 1, the flow of transition is a complex series changing fromthe steady to periodic, chaotic, periodic, quasi-periodic, and finally back to periodic. There are two types of periodic flow, one is simple periodic flow with single fundamental frequency (FF), and another is complex periodic flow with multiple FFs. This process is illustrated by using time-velocity histories, Fourier frequency spectrum analysis and the phase-space rajectories.

Kramers Escape Rate in Nonlinear Diffusive Media

In this paper, we study nonlinear Kramers problem by investigating overdamped systems ruled by the one-dimensional nonlinear Fokker-Planck equation. We obtain an analytic expression for the Kramers escape rate under quasistationary conditions by employing

Diffusive Wave Solutions for Open Channel Flows with Uniform and Concentrated Lateral Inflow

The diffusive wave equation with inhomogeneous terms representing hydraulics with uniform or concentrated lateral inflow intoa river is theoretically investigated in the current paper. All the solutions have been systematically expressed in a unified form interms of response function or so called K-function. The integration of K-function obtained by using Laplace transform becomesS-function, which is examined in detail to improve the understanding of flood routing characters. The backwater effects usuallyresulting in the discharge reductions and water surface elevations upstream due to both the downstream boundary and lateral infloware analyzed. With a pulse discharge in upstream boundary inflow, downstream boundary outflow and lateral inflow respectively,hydrographs of a channel are routed by using the S-functions. Moreover, the comparisons of hydrographs in infinite, semi-infiniteand finite channels are pursued to exhibit the different backwater effects due to a concentrated lateral inflow for various channeltypes.

Experimental Study on Convective Boiling Heat Transfer in Narrow-Gap Annulus Tubes

Since convective boiling or highly subcooled single-phase forced convection in micro-channels is an effective cooling mechanism with a wide range of applications, more experimental and theoretical studies are required to explain and verify the forced convection heat transfer phenomenon in narrow channels. In this experimental study, we model the convective boiling behavior of water with low latent heat substance Freon 113 (R-113), with the purpose of saving power consumption and visualizing experiments. Both heat transfer and pressure drop characteristics were measured in subcooled and saturated concentric narrow gap forced convection boiling. Data were obtained to qualitatively identify the effects of gap size, pressure, flow rate and wall superheat on boiling regimes and the transition between various regimes. Some significant differences from unconfined forced convection boiling were found,and also, the flow patterns in narrow vertical annulus tubes have been studied quantitatively.

Convective coupled map for simulating spatiotemporal chaos in flows

A coupled map lattices with convective nonlinearity or, for short, Convective Coupled Map (CCM) is proposed in this paper to simulate spatiotemporal chaos in fluid hows. It is found that the parameter region of spatiotemporal chaos can be determined by the maximal Liapunov exponent of its complexity time series. This simple model implies a similar physical mechanism for turbulence such that the route to spatiotemporal chaos in fluid hows can be envisaged.

Air flow over a mountain and the convective boundary-layer

In this paper, the analytical model coupling the convective boundary layer (CBL) with the free atmosphere developed by Qi and Fu (1992) is improved. And by this improved model, the interaction between airflow over a mountain and the CBL is further discussed. The conclusions demonstrate: (1) The perturbation potential temperatures in the free atmosphere can counteract the effect of orographic thermal forcing through entraining and mixing in the CBL. If u(M)BAR > u(F)BAR, the feedback of the perturbation potential temperatures in the free atmosphere is more important than orographic thermal forcing, which promotes the effect of interfacial waves. If u(M)BAR < u(F)BAR, orographic thermal forcing is more important, which makes the interfacial height and the topographic height identical in phase, and the horizontal speeds are a maximum at the top of the mountain. (2) The internal gravity waves propagating vertically in the free atmosphere cause a strong downslope wind to become established above the lee slope in the CBL and result in the hydraulic jump at the top of the CBL. (3) With the CBL deepening, the interfacial gravity waves induced by the potential temperature jump at the top of the CBL cause the airflow in the CBL to be subcritical.

A perturbational h4 exponential finite-difference scheme for the convective diffusion equation

A perturbational h4 compact exponential finite difference scheme with diagonally dominant coefficient matrix and upwind effect is developed for the convective diffusion equation. Perturbations of second order are exerted on the convective coefficients and source term of an h2 exponential finite difference scheme proposed in this paper based on a transformation to eliminate the upwind effect of the convective diffusion equation. Four numerical examples including one- to three-dimensional model equations of fluid flow and a problem of natural convective heat transfer are given to illustrate the excellent behavior of the present exponential schemes, the h4 accuracy of the perturbational scheme is verified using double precision arithmetic.

H4 exponential finite-difference scheme for convective diffusion-equations

Perturbations are applied to the convective coefficients and source term of a convection-diffusion equation so that second-order corrections may be applied to a second-order exponential scheme. The basic Structure of the equations in the resulting fourth-order scheme is identical to that for the second order. Furthermore, the calculations are quite simple as the second-order corrections may be obtained in a single pass using a second-order scheme. For one to three dimensions, the fourth-order exponential scheme is unconditionally stable. As examples, the method is applied to Burgers' and other fluid mechanics problems. Compared with schemes normally used, the accuracies are found to be good and the method is applicable to regions with large gradients.

Turbulent passive scalar field of a small prandtl number

The statistical-mechanics theory of the passive scalar field convected by turbulence, developed in an earlier paper [Phys. Fluids 28, 1299 (1985)], is extended to the case of a small molecular Prandtl number. The set of governing integral equations is solved by the equation-error method. The resultant scalar-variance spectrum for the inertial range is F(k)~x−5/3/[1+1.21x1.67(1+0.353x2.32)], where x is the wavenumber scaled by Corrsin's dissipation wavenumber. This result reduces to the − (5)/(3) law in the inertial-convective range. It also approximately reduces to the − (17)/(3) law in the inertial-diffusive range, but the proportionality constant differs from Batchelor's by a factor of 3.6.

Rayleigh-Marangoni-Bénard Convective Instability

The Rayleigh-Marangoni-Benard convective instability (R-M-B instability) and flow patterns in the two-layer system of silicon oil 10cSt and Fluorinert FC70 liquids are studied theoretically and experimentally. Both linear instability analysis and 2D numerical simulation (A=L/H=10) were performed to study the influence of thermocapillary force on the convective instability of the two-layer system. Time-dependent oscillations arising at the onset of convection were investigated in a larger various range of two-layer depth ratios (Hr=H1/H2) from 0.2 to 5.0 for different total depth less than 12mm. Our results are different from the previous study on the Rayleig-B閚ard instability and show the strong effects of thermocapillary force at the interface on the time-dependent oscillations at the onset of instability convection. Primary experimental results of the critical instability parameters and the convective structure in the R-M-B convection have been obtained by using the digital particle image velocimetry (DPIV) system, and a good agreement in comparison with the results of numerical simulation was obtained.

Numerical-Value Perturbation Method with Application of Solving Convective-Diffusion Integral Equation

The convective--diffusion equation is of primary importance in such fields as fluid dynamics and heat transfer hi the numerical methods solving the convective-diffusion equation, the finite volume method can use conveniently diversified grids (structured and unstructured grids) and is suitable for very complex geometry The disadvantage of FV methods compared to the finite difference method is that FV-methods of order higher than second are more difficult to develop in three-dimensional cases. The second-order central scheme (2cs) offers a good compromise among accuracy, simplicity and efficiency, however, it will produce oscillatory solutions when the grid Reynolds numbers are large and then very fine grids are required to obtain accurate solution. The simplest first-order upwind (IUW) scheme satisfies the convective boundedness criteria, however. Its numerical diffusion is large. The power-law scheme, QMCK and second-order upwind (2UW) schemes are also often used in some commercial codes. Their numerical accurate are roughly consistent with that of ZCS. Therefore, it is meaningful to offer higher-accurate three point FV scheme. In this paper, the numerical-value perturbational method suggested by Zhi Gao is used to develop an upwind and mixed FV scheme using any higher-order interpolation and second-order integration approximations, which is called perturbational finite volume (PFV) scheme. The PFV scheme uses the least nodes similar to the standard three-point schemes, namely, the number of the nodes needed equals to unity plus the face-number of the control volume. For instanc6, in the two-dimensional (2-D) case, only four nodes for the triangle grids and five nodes for the Cartesian grids are utilized, respectively. The PFV scheme is applied on a number of 1-D problems, 2~Dand 3-D flow model equations. Comparing with other standard three-point schemes, The PFV scheme has much smaller numerical diffusion than the first-order upwind (IUW) scheme, its numerical accuracy are also higher than the second-order central scheme (2CS), the power-law scheme (PLS), the QUICK scheme and the second-order upwind(ZUW) scheme.

Diffusive waves in a channel with concentrated lateral inflow

In the current paper an analytical solution for diffusive wave equation with the concentrate-distributed lateral inflow is yielded. Finite-difference numerical method is also employed to validate this model. The backwater effects drawn from lateral inflow on the mainstream are examined finally.

The Convective Instabilities in a Liquid-Vapor System with a Non-equilibrium Evaporation Interface

Rayleigh-Marangoni-B,nard instability in a system consisting of a horizontal liquid layer and its own vapor has been investigated. The two layers are separated by a deformable evaporation interface. A linear stability analysis is carried out to study the convective instability during evaporation. In previous works, the interface is assumed to be under equilibrium state. In contrast with previous works, we give up the equilibrium assumption and use Hertz-Knudsen's relation to describe the phase change under non-equilibrium state. The influence of Marangoni effect, gravitational effect, degree of non-equilibrium and the dynamics of the vapor on the instability are discussed.

Diffusive coevolution and mutualism maintenance mechanisms in a fig-fig wasp system

In reciprocal mutualism systems, the exploitation events by exploiters might disrupt the reciprocal mutualism, wherein one exploiter species might even exclude other coexisting exploiter species over an evolutionary time frame. What remains unclear is how such a community is maintained. Niche partitioning, or spatial heterogeneity among the mutualists and exploiters, is generally believed to enable stability within a mutualistic system. However, our examination of a reciprocal mutualism between a fig species (Ficus racemosa) and its pollinator wasp (Ceratosolen fusciceps) shows that spatial niche partitioning does not sufficiently prevent exploiters from overexploiting the common resource (i.e., the female flowers), because of the considerable niche overlap between the mutualists and exploiters. In response to an exploiter, our experiment shows that the fig can (1) abort syconia-containing flowers that have been galled by the exploiter, Apocryptophagus testacea, which oviposits before the pollinators do; and (2) retain syconia-containing flowers galled by Apocryptophagus mayri, which oviposit later than pollinators. However, as a result of (2), there is decreased development of adult non-pollinators or pollinator species in syconia that have not been sufficiently pollinated, but not aborted. Such discriminative abortion of figs or reduction in offspring development of exploiters while rewarding cooperative individuals with higher offspring development by the fig will increase the fitness of cooperative pollinating wasps, but decrease the fitness of exploiters. The fig fig wasp interactions are diffusively coevolved, a case in which fig wasps diversify their genotype, phenotype, or behavior as a result of competition between wasps, while figs diverge their strategies to facilitate the evolution of cooperative fig waps or lessen the detrimental behavior by associated fig wasps. In habitats or syconia that suffer overexploitation, discriminative abortion of figs or reduction in the offspring development of exploiters in syconia that are not or not sufficiently pollinated will decrease exploiter fitness and perhaps even drive the population of exploiters to local extinction, enabling the evolution and maintenance of cooperative pollinators through the movement between habitats or syconia (i.e., the metapopulations).