Finite element modelling of three-dimensional convection problems in fluid-saturated porous media heated from below

We use the finite element method to model three-dimensional convective pore-fluid flow in fluid-saturated porous media when they are heated from below. In particular, we employ the particle-tracking technique to mimic the trajectories of particles in three-dimensional fluid flow problems. The related numerical results demonstrated that: (1) The progressive asymptotic approach procedure, which was previously developed for the finite element modelling of two-dimensional convective pore-fluid flow problems, is equally applicable to the finite element modelling of three-dimensional convective pore-fluid flow in fluid-saturated porous media heated from below. (2) The perturbation of gravity at different planes has a significant effect on the pattern of three-dimensional convective pore-fluid flow and therefore, may influence the pattern of orebody formation and mineralization in three-dimensional hydrothermal systems. Copyright (C) 2001 John Wiley & Sons, Ltd.

SURFACE SOLITARY WAVES IN A DOUBLE DIFFUSIVE SYSTEM

We show that a surface solitary wave governed by the Korteweg-de Vries equation can develop in a fluid acted upon by fluxes of heat and of a second diffusive element. This solitary wave appears as a manifestation of a hydrodynamical instability which sets in only when a certain relation involving the parameters of the system is satisfied.

An experimental study of the effects of double-diffusive convention processes during the solification of binary alloys

Experiments were conducted to show the effects of thermal and geometric boundary conditions on the liquid pool of a binary alloy system which is undergoing phase change, solidification. Transparent analogue solutions were selected for study and experimental apparatus were designed and built. Thermal distribution and concentration data were collected and analysed for the melt pool of various selected geometries and boundary conditions of the systems under study. The data indicate-that characteristic flows develop for both Hypereutectic and Hypoeutectic concentration levels and that the development of macrosegregation and microsegregation defects in continuous casting materials can be minimised by the adjustment of the process variables.

Effects of geological inhomogeneity on high Rayleigh number steady state heat and mass transfer in fluid-saturated porous media heated from below

A parametric study is carried out to investigate how geological inhomogeneity affects the pore-fluid convective flow field, the temperature distribution, and the mass concentration distribution in a fluid-saturated porous medium. The related numerical results have demonstrated that (1) the effects of both medium permeability inhomogeneity and medium thermal conductivity inhomogeneity are significant on the pore-fluid convective flow and the species concentration distribution in the porous medium; (2) the effect of medium thermal conductivity inhomogeneity is dramatic on the temperature distribution in the porous medium, but the effect of medium permeability inhomogeneity on the temperature distribution may be considerable, depending on the Rayleigh number involved in the analysis; (3) if the coupling effect between pore-fluid flow and mass transport is weak, the effect of the Lewis number is negligible on the pore-fluid convective flow and temperature distribution, hut it is significant on the species concentration distribution in the medium.

Microstructures layering, mechanisms in the middle waters of the Persian Gulf

Layered structures, known as micro structures in marine environments are common features of which their formation mechanisms are first reviewed. Some methods of measuring such features based on the measurements and theories are presented for the Persian Gulf. This includes determination of layers with temperature inversion (TI) associated with double diffusive convection (DDC). The relevant associated parameters are estimated from ROPME CTD data for late winter and early summer of 1992. Only in certain parts temperature inversion and DDC are observed which seem to produce layered structures. Observations show that the places with TI and DDC are mainly confined to the frontal regions where the water entering the Persian Gulf and water exiting it meet, nearly along the axis of the Gulf. TI and DDC is mainly observer in the northern bound of the front. Typical density ratio for regions with TI and DDC is 0.7 to 0.2 and the mean depth is at about 37 ± 3 m for the Persian Gulf. TI and DDC are also found in the outflow from the Persian Gulf to the Oman Gulf which is found to be at a depth of about 250 m. Horizontal addiction and reduction of solar heating seem to be the main reasons in producing layers with TI and DDC. It is also found that the regime of DDC in the Persian Gulf is more diffusive and the flow associated with intrusion layers with TI is non-isopycnal (more unstable). However for the Oman sea both diffusive and finger regime are observed and the flow is inferred to be isopycnal (more stable statically). Typical heat and salt fluxes due to DDC are found to be 6 W/m2 and 0.36 W/m2 respectively. Effective salinity diffusivity, Ks and heat diffusivity, Kr have been estimated for the places with DDC in the Persian Gulf and Oman Gulf (Ks=1.1 *10-7 m2/s, KT= 1.88*10-6 m2/s). Their values are within the values obtained by others. The buoyancy frequency for the Persian Gulf with typical mean value of 0.05s-1 is much higher than these of the free Oceans. Such large values of N (typically 0.05 s-1) indicate that processes such as tide can produce strong internal waves which may be another factor in producing layered structures. This requires separate study.

Small scale turbulence and mixing in semi-enclosed seas (Persian Gulf)

Turbulence and internal waves are probably important in generating layered structures in frontal region of marine environments (e.g. near river plumes outflow into the sea). Here we investigate the role of normal modes of internal waves in generation of layered structure in a part of Persian Gulf where river plume inters and in some laboratory experiments. The model prediction and observations show that layers so formed have a thickness of about 2m based on salinity variations with depth, but layers (about 5m) based on horizontal velocity profiles. Laboratory experiments with a plume outflow in a Filling Box profile also generate normal mode layered structure with 21H=0.5 (where A is layer thickness and H is the plume depth). In these experiments as Re of the flow is smaller than the Re of field flow. The normal modes are substantially dissipated with depth. Typical values of dissipation factor is about 0(100). This factor for field observation is 0(10) which is still substantial. Qualitative comparison between layered structure in field and laboratory is good. It should be emphasized that field observation is for semi-enclosed seas but the laboratory experiments are for enclosed region. Hence some of the discrepancies in the results of two cases are inevitable. Layered structures in marine environments are also produced by double diffusive convection. In this region this should be studied separately.

Concentrations of organic carbon, nitrogen and phosphorus on vertical profiles in waters of the tropical North Atlantic in March 2002

In the literature, an inconsistency exists between estimates of biotically-effected carbon export inferred from large-scale geochemical studies (Jenkins 1982; 47 gC m-2 a-1) and local measurements of turbulent nutrient supply (Lewis et al. 1986; 4 gC m-2 a-1) in the eastern subtropical North Atlantic. Nutrient supply to the upper ocean by turbulent mixing is reexamined using local standard oceanographic measurements and high-resolution vertical profiles of nutrients averaged over a large region directly comparable to that investigated by Jenkins (1982). Turbulent fluxes induced by internal waves and salt fingering, respectively, are separated according to Gregg (1989) and Zhang et al. (1998). Nutrient transport into the nutrient-consuming surface layer by salt fingering is more than fivefold higher than transport due to internal-wave induced turbulence. Still, this cannot resolve the above- mentioned apparent inconsistency, even if additional physical transport mechanisms such as eddy pumping, advection and horizontal diffusion are accounted for. Estimated nitrate fluxes due to vertical turbulent diffusion are 0.05-0.15 mol m-2 a-1, corresponding to 4-11 gC m-2 a-1. Observed NO3/PO4 turbulent flux ratios of up to 23 are interpreted as the imprint of N2 fixation.

Numerical modelling of double diffusion driven reactive flow transport in deformable fluid-saturated porous media with particular consideration of temperature-dependent chemical reaction rates

Numerical methods ave used to solve double diffusion driven reactive flow transport problems in deformable fluid-saturated porous media. in particular, thp temperature dependent reaction rate in the non-equilibrium chemical reactions is considered. A general numerical solution method, which is a combination of the finite difference method in FLAG and the finite element method in FIDAP, to solve the fully coupled problem involving material deformation, pore-fluid flow, heat transfer and species transport/chemical reactions in deformable fluid-saturated porous media has been developed The coupled problem is divided into two subproblems which are solved interactively until the convergence requirement is met. Owing to the approximate nature of the numerical method, if is essential to justify the numerical solutions through some kind of theoretical analysis. This has been highlighted in this paper The related numerical results, which are justified by the theoretical analysis, have demonstrated that the proposed solution method is useful for and applicable to a wide range of fully coupled problems in the field of science and engineering.

Numerical simulation of double-diffusion driven convective flow and rock alteration in three-dimensional fluid-saturated geological fault zones

Numerical methods are used to simulate the double-diffusion driven convective pore-fluid flow and rock alteration in three-dimensional fluid-saturated geological fault zones. The double diffusion is caused by a combination of both the positive upward temperature gradient and the positive downward salinity concentration gradient within a three-dimensional fluid-saturated geological fault zone, which is assumed to be more permeable than its surrounding rocks. In order to ensure the physical meaningfulness of the obtained numerical solutions, the numerical method used in this study is validated by a benchmark problem, for which the analytical solution to the critical Rayleigh number of the system is available. The theoretical value of the critical Rayleigh number of a three-dimensional fluid-saturated geological fault zone system can be used to judge whether or not the double-diffusion driven convective pore-fluid flow can take place within the system. After the possibility of triggering the double-diffusion driven convective pore-fluid flow is theoretically validated for the numerical model of a three-dimensional fluid-saturated geological fault zone system, the corresponding numerical solutions for the convective flow and temperature are directly coupled with a geochemical system. Through the numerical simulation of the coupled system between the convective fluid flow, heat transfer, mass transport and chemical reactions, we have investigated the effect of the double-diffusion driven convective pore-fluid flow on the rock alteration, which is the direct consequence of mineral redistribution due to its dissolution, transportation and precipitation, within the three-dimensional fluid-saturated geological fault zone system. (c) 2005 Elsevier B.V. All rights reserved.

Numerical study on mixed convection in a horizontal flow past a square porous cavity using UNIFAES scheme

Mixed convection on the flow past a heated length and past a porous cavity located in a horizontal wall bounding a saturated porous medium is numerically simulated. The cavity is heated from below. The steady-state regime is studied for several intensities of the buoyancy effects due to temperature variations. The influences of Péclet and Rayleigh numbers on the flow pattern and the temperature distributions are examined. Local and global Nusselt numbers are reported for the heated surface. The convective-diffusive fluxes at the volume boundaries are represented using the UNIFAES, Unified Finite Approach Exponential-type Scheme, with the Power-Law approximation to reduce the computing time. The conditions established by Rivas for the quadratic order of accuracy of the central differencing to be maintained in irregular grids are shown to be extensible to other quadratic schemes, including UNIFAES, so that accuracy estimates could be obtained.

Using a Double-pass solar drier for drying of bamboo shoots

Three different drying methods, a forced convection double-pass solar drier (DPSD), typical cabinet type natural convection solar drier (CD) and traditional open-sun drying (OSD) were used for draying of bamboo shoots in central Vietnam. During drying the operational parameters such as drying temperature, relative humidity, air velocity, insolation and water evaporation have been recorded hourly. The mean drying temperatures and relative humidity in the drying chamber were 55.2°C, 23.7%; 47.5°C, 37,6%; 36.2°C, 47.8% in DPSD, CD and OSD, respectively. The mean global radiation during all experimental runs was 670 Wm^−2. The result also shows that fastest drying process was occurred in DPSD where the falling-rate period was achieved after 7 hours, in change to OSD where it took 16 hours. The overall drying efficiency was 23.11%, 15.83% and 9.73% in case of DPSD, CD and OSD, respectively. Although the construction cost of DPSD was significantly higher than in CD, the drying costs per one kilogram of bamboo shoots were by 42.8% lower in case of DPSD as compared to CD. Double-pass solar drier was found to be technically and economically suitable for drying of bamboo shoots under the specific conditions in central Vietnam and in all cases, the use of this drier led to considerable reduction in drying time in comparison to traditional open-sun drying.

A case study of boundary layer ventilation by convection and coastal processes

It is often assumed that ventilation of the atmospheric boundary layer is weak in the absence of fronts, but is this always true? In this paper we investigate the processes responsible for ventilation of the atmospheric boundary layer during a nonfrontal day that occurred on 9 May 2005 using the UK Met Office Unified Model. Pollution sources are represented by the constant emission of a passive tracer everywhere over land. The ventilation processes observed include shallow convection, turbulent mixing followed by large-scale ascent, a sea breeze circulation and coastal outflow. Vertical distributions of tracer are validated qualitatively with AMPEP (Aircraft Measurement of chemical Processing Export fluxes of Pollutants over the UK) CO aircraft measurements and are shown to agree impressively well. Budget calculations of tracers are performed in order to determine the relative importance of these ventilation processes. Coastal outflow and the sea breeze circulation were found to ventilate 26% of the boundary layer tracer by sunset of which 2% was above 2 km. A combination of coastal outflow, the sea breeze circulation, turbulent mixing and large-scale ascent ventilated 46% of the boundary layer tracer, of which 10% was above 2 km. Finally, coastal outflow, the sea breeze circulation, turbulent mixing, large-scale ascent and shallow convection together ventilated 52% of the tracer into the free troposphere, of which 26% was above 2 km. Hence this study shows that significant ventilation of the boundary layer can occur in the absence of fronts (and thus during high-pressure events). Turbulent mixing and convection processes can double the amount of pollution ventilated from the boundary layer.

Experimental and Numerical Analysis of Gas Forced Convection through Microtubes and Micro Heat Exchangers

The last decade has witnessed very fast development in microfabrication technologies. The increasing industrial applications of microfluidic systems call for more intensive and systematic knowledge on this newly emerging field. Especially for gaseous flow and heat transfer at microscale, the applicability of conventional theories developed at macro scale is not yet completely validated; this is mainly due to scarce experimental data available in literature for gas flows. The objective of this thesis is to investigate these unclear elements by analyzing forced convection for gaseous flows through microtubes and micro heat exchangers. Experimental tests have been performed with microtubes having various inner diameters, namely 750 m, 510 m and 170 m, over a wide range of Reynolds number covering the laminar region, the transitional zone and also the onset region of the turbulent regime. The results show that conventional theory is able to predict the flow friction factor when flow compressibility does not appear and the effect of fluid temperature-dependent properties is insignificant. A double-layered microchannel heat exchanger has been designed in order to study experimentally the efficiency of a gas-to-gas micro heat exchanger. This microdevice contains 133 parallel microchannels machined into polished PEEK plates for both the hot side and the cold side. The microchannels are 200 µm high, 200 µm wide and 39.8 mm long. The design of the micro device has been made in order to be able to test different materials as partition foil with flexible thickness. Experimental tests have been carried out for five different partition foils, with various mass flow rates and flow configurations. The experimental results indicate that the thermal performance of the countercurrent and cross flow micro heat exchanger can be strongly influenced by axial conduction in the partition foil separating the hot gas flow and cold gas flow.

Charge effects on hindrance factors for diffusion and convection of solute in pores I

The transport of a spherical solute through a long circular cylindrical pore filled with an electrolyte solution is studied numerically, in the presence of constant surface charge on the solute and the pore wall. Fluid dynamic analyses were carried out to calculate the flow field around the solute in the pore to evaluate the drag coefficients exerted on the solute. Electrical potentials around the solute in the electrolyte solution were computed based on a mean-field theory to provide the interaction energy between the charged solute and the pore wall. Combining the results of the fluid dynamic and electrostatic analyses, we estimated the rate of the diffusive and convective transport of the solute across the pore. Although the present estimates of the drag coefficients on the solute suggest more than 10% difference from existing studies, depending on the radius ratio of the solute relative to the pore and the radial position of the solute center in the pore, this difference leads to a minor effect on the hindrance factors. It was found that even at rather large ion concentrations, the repulsive electrostatic interaction between the charged solute and the pore wall of like charge could significantly reduce the transport rate of the solute.