1
|
Janardhana Reddy G, Basha H, Venkata Narayanan N. A numerical investigation of transient natural convective heat transfer to isobutane in the supercritical region. J Mol Liq 2018. [DOI: 10.1016/j.molliq.2017.11.074] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/18/2022]
|
2
|
Tilgner A. Convection in an ideal gas at high Rayleigh numbers. PHYSICAL REVIEW. E, STATISTICAL, NONLINEAR, AND SOFT MATTER PHYSICS 2011; 84:026323. [PMID: 21929106 DOI: 10.1103/physreve.84.026323] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 03/13/2011] [Revised: 06/23/2011] [Indexed: 05/31/2023]
Abstract
Numerical simulations of convection in a layer filled with ideal gas are presented. The control parameters are chosen such that there is a significant variation of density of the gas in going from the bottom to the top of the layer. The relations between the Rayleigh, Peclet, and Nusselt numbers depend on the density stratification. It is proposed to use a data reduction which accounts for the variable density by introducing into the scaling laws an effective density. The relevant density is the geometric mean of the maximum and minimum densities in the layer. A good fit to the data is then obtained with power laws with the same exponent as for fluids in the Boussinesq limit. Two relations connect the top and bottom boundary layers: The kinetic energy densities computed from free fall velocities are equal at the top and bottom, and the products of free fall velocities and maximum horizontal velocities are equal for both boundaries.
Collapse
Affiliation(s)
- A Tilgner
- Institute of Geophysics, University of Göttingen, Friedrich-Hund-Platz 1, D-37077 Göttingen, Germany
| |
Collapse
|
3
|
Shen B, Zhang P. Thermoacoustic waves along the critical isochore. PHYSICAL REVIEW. E, STATISTICAL, NONLINEAR, AND SOFT MATTER PHYSICS 2011; 83:011115. [PMID: 21405669 DOI: 10.1103/physreve.83.011115] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 08/26/2010] [Indexed: 05/30/2023]
Abstract
Near the liquid-gas critical point, thermal disturbances can generate sounds. We study the acoustic emission over four decades of reduced temperatures [defined as ɛ=(T-T(c))/T(c), with T(c) the critical temperature] along the critical isochore, under linear and nonlinear temperature perturbations, respectively. We identify various thermoacoustic behaviors by numerically solving the governing equations. It is shown that a homogeneous thermoacoustic-wave pattern dominates in the linear case, largely independent of ɛ; whereas under the nonlinear perturbation, variation in ɛ could lead to severe wavefront deformation. The strong nonlinear effect is found to be of a transient nature because, in due time, both cases tend to converge in terms of the energy yield of the adiabatic process.
Collapse
Affiliation(s)
- B Shen
- Institute of Refrigeration and Cryogenics, Shanghai Jiao Tong University, Shanghai, China
| | | |
Collapse
|
4
|
Házi G, Márkus A. Modeling heat transfer in supercritical fluid using the lattice Boltzmann method. PHYSICAL REVIEW. E, STATISTICAL, NONLINEAR, AND SOFT MATTER PHYSICS 2008; 77:026305. [PMID: 18352120 DOI: 10.1103/physreve.77.026305] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 10/03/2006] [Revised: 11/26/2007] [Indexed: 05/26/2023]
Abstract
A lattice Boltzmann model has been developed to simulate heat transfer in supercritical fluids. A supercritical viscous fluid layer between two plates heated from the bottom has been studied. It is demonstrated that the model can be used to study heat transfer near the critical point where the so-called piston effect speeds up the transfer of heat and results in homogeneous heating in the bulk of the layer. We have also studied the onset of convection in a Rayleigh-Bénard configuration. It is shown that our model can well predict qualitatively the onset of convection near the critical point, where there is a crossover between the Rayleigh and Schwarzschild criteria.
Collapse
Affiliation(s)
- Gábor Házi
- Simulator Development Department, KFKI Atomic Energy Research Institute, H-1525 Budapest, Hungary.
| | | |
Collapse
|
5
|
Onuki A. Thermoacoustic effects in supercritical fluids near the critical point: Resonance, piston effect, and acoustic emission and reflection. PHYSICAL REVIEW. E, STATISTICAL, NONLINEAR, AND SOFT MATTER PHYSICS 2007; 76:061126. [PMID: 18233833 DOI: 10.1103/physreve.76.061126] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 08/31/2007] [Indexed: 05/25/2023]
Abstract
We present a general theory of thermoacoustic phenomena in one phase states of one-component fluids. Singular behavior is predicted in supercritical fluids near the critical point. In a one-dimensional geometry we start with linearized hydrodynamic equations taking into account the effects of heat conduction in the boundary walls and the bulk viscosity. We introduce a coefficient Z(omega) characterizing reflection of sound with frequency omega at the boundary in a rigid cell. As applications, we examine acoustic eigenmodes, response to time-dependent perturbations, and sound emission and reflection. Resonance and rapid adiabatic changes are noteworthy. In these processes, the role of the thermal diffusion layers is enhanced near the critical point because of the strong critical divergence of the thermal expansion.
Collapse
Affiliation(s)
- Akira Onuki
- Department of Physics Kyoto University, Kyoto 606-8502, Japan
| |
Collapse
|
6
|
Onuki A. Dynamic van der Waals theory. PHYSICAL REVIEW. E, STATISTICAL, NONLINEAR, AND SOFT MATTER PHYSICS 2007; 75:036304. [PMID: 17500788 DOI: 10.1103/physreve.75.036304] [Citation(s) in RCA: 45] [Impact Index Per Article: 2.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 10/26/2006] [Indexed: 05/15/2023]
Abstract
We present a dynamic van der Waals theory starting with entropy and energy functional with gradient contributions. The resultant hydrodynamic equations contain the stress arising from the density gradient. It provides a general scheme of two-phase hydrodynamics involving the gas-liquid transition in nonuniform temperature. Some complex hydrodynamic processes with evaporation and condensation are examined numerically. They are (i) adiabatically induced spinodal decomposition, (ii) piston effect with a bubble in liquid, (iii) temperature and velocity profiles around a droplet in heat flow, (iv) efficient latent heat transport at small liquid densities (the mechanism of heat pipes), (v) boiling in gravity with continuous bubble formation and rising, and (vi) spreading and evaporation of liquid on a heated boundary wall.
Collapse
Affiliation(s)
- Akira Onuki
- Department of Physics, Kyoto University, Kyoto 606-8502, Japan
| |
Collapse
|
7
|
Accary G, Meyer H. Perturbation-controlled numerical simulations of the convection onset in a supercritical fluid layer. PHYSICAL REVIEW. E, STATISTICAL, NONLINEAR, AND SOFT MATTER PHYSICS 2006; 74:046308. [PMID: 17155174 DOI: 10.1103/physreve.74.046308] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 04/01/2006] [Indexed: 05/12/2023]
Abstract
This paper reports 2D and 3D direct numerical simulations of the Rayleigh-Bénard convection in a fluid layer, 3He, close to its gas-liquid critical point. The main quantity of interest is the time when the convective motion becomes appreciable after the heat current has been started. A space periodic, time-independent temperature perturbation is applied to the top plate aiming to represent the various sources of noise in the experiments. A single amplitude of this additional perturbation reproduces the noise level of the experiments for all the used combinations of the critical point proximity and of the heating power. Over four decades of the Rayleigh number, this simple operation removes the systematic discrepancy reported in an earlier paper and brings the simulations into good agreement with the measurements. Scaling of the exponential convection growth rate is presented and discussed.
Collapse
Affiliation(s)
- G Accary
- MSNM-GP, UMR 6181 CNRS, IMT - La Jetée - Technopôle de Château-Gombert, 38 rue Frédéric Joliot Curie, 13451 Marseille, France.
| | | |
Collapse
|
8
|
Meyer H. Onset of the convection in a supercritical fluid. PHYSICAL REVIEW. E, STATISTICAL, NONLINEAR, AND SOFT MATTER PHYSICS 2006; 73:016311. [PMID: 16486281 DOI: 10.1103/physreve.73.016311] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 08/30/2005] [Indexed: 05/06/2023]
Abstract
A model is proposed that leads to the scaled relation tp/tau D=Ftp(Ra-Rac) for the development of convection in a pure fluid in a Rayleigh-Bénard cell after the start of the heat current at t=0. Here tp is the time of the first maximum of the temperature drop DeltaT(t) across the fluid layer, the signature of rapidly growing convection, tau D is the diffusion relaxation time, and Rac is the critical Rayleigh number. Such a relation was first obtained empirically from experimental data. Because of the unknown perturbations in the cell that lead to convection development beyond the point of the fluid instability, the model determines tp/tau D within a multiplicative factor Psi square root Rac(HBL), the only fit parameter product. Here Rac(HBL), of the order 10(3), is the critical Rayleigh number of the hot boundary layer and Psi is a fit parameter. There is then good agreement over more than four decades of Ra-Rac between the model and the experiments on supercritical 3He at various heat currents and temperatures. The value of the parameter Psi, which phenomenologically represents the effectiveness of the perturbations, is discussed in connection with predictions by El Khouri and Carlès of the fluid instability onset time.
Collapse
Affiliation(s)
- H Meyer
- Department of Physics, Duke University, Durham, North Carolina 27708-0305, USA
| |
Collapse
|
9
|
Onuki A, Kanatani K. Droplet motion with phase change in a temperature gradient. PHYSICAL REVIEW. E, STATISTICAL, NONLINEAR, AND SOFT MATTER PHYSICS 2005; 72:066304. [PMID: 16486056 DOI: 10.1103/physreve.72.066304] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 06/20/2005] [Indexed: 05/06/2023]
Abstract
We examine the droplet motion in one-component fluids in a small temperature gradient by solving linearized hydrodynamic equations supplemented with appropriate surface boundary conditions. We show that the velocity field and the temperature around the droplet are strongly influenced by first-order phase transition taking place at the interface. Latent heat released or absorbed at the interface drastically changes the hydrodynamic flow around the droplet. As a result, the temperature becomes almost homogeneous inside the droplet and the Marangoni effect arising from the surface tension gradient is much suppressed. The droplet velocity is also much decelerated.
Collapse
Affiliation(s)
- Akira Onuki
- Department of Physics, Kyoto University, Kyoto 606-8502, Japan
| | | |
Collapse
|
10
|
Hwa RC, Yang CB, Bershadskii S, Niemela JJ, Sreenivasan KR. Critical fluctuation of wind reversals in convective turbulence. PHYSICAL REVIEW. E, STATISTICAL, NONLINEAR, AND SOFT MATTER PHYSICS 2005; 72:066308. [PMID: 16486060 DOI: 10.1103/physreve.72.066308] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 10/21/2005] [Indexed: 05/06/2023]
Abstract
The irregular reversals of wind direction in convective turbulence are found to have fluctuating intervals that can be related, under certain circumstances, to critical behavior. In particular, by focusing on its temporal evolution, the net magnetization of a two-dimensional Ising lattice of finite size is observed to fluctuate in the same way. Detrended fluctuation analysis of the wind reversal time series results in a scaling behavior that agrees remarkably well with that of the Ising problem. The specific properties found here, as well as the lack of an external tuning parameter, also suggest that the wind reversal phenomenon exhibits signs of self-organized criticality.
Collapse
Affiliation(s)
- Rudolph C Hwa
- Institute of Theoretical Science and Department of Physics, University of Oregon, Eugene, Oregon 97403-5203, USA
| | | | | | | | | |
Collapse
|
11
|
Accary G, Raspo I, Bontoux P, Zappoli B. Reverse transition to hydrodynamic stability through the Schwarzschild line in a supercritical fluid layer. PHYSICAL REVIEW. E, STATISTICAL, NONLINEAR, AND SOFT MATTER PHYSICS 2005; 72:035301. [PMID: 16241504 DOI: 10.1103/physreve.72.035301] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 10/18/2004] [Indexed: 05/05/2023]
Abstract
We consider a fluid close to its gas-liquid critical point in a bottom-heated cavity. Due to strong density stratification, both the Schwarzschild and the Rayleigh stability criteria are relevant at the same space scale. Taking advantage of the competition between these two limits of the convection-onset criterion, we numerically exhibit striking non-Boussinesq effects: the reverse transition to stability through the Schwarzschild line of a heat diffusing layer subjected to convection, and the convection onset inside a few-millimeters-thick layer according to the Schwarzschild criterion.
Collapse
Affiliation(s)
- Gilbert Accary
- MSNM-GP UMR 6181 CNRS, Les Universités Aix-Marseille, IMT-La Jetée, Technopôle de Château Gombert, 38 rue Frédéric Joliot Curie, 13451 Marseille Cedex 20, France
| | | | | | | |
Collapse
|
12
|
Araujo FF, Grossmann S, Lohse D. Wind reversals in turbulent Rayleigh-Bénard convection. PHYSICAL REVIEW LETTERS 2005; 95:084502. [PMID: 16196862 DOI: 10.1103/physrevlett.95.084502] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 07/13/2004] [Indexed: 05/04/2023]
Abstract
The phenomenon of irregular cessation and subsequent reversal of the large-scale circulation in turbulent Rayleigh-Bénard convection is theoretically analyzed. The force and thermal balance on a single plume detached from the thermal boundary layer yields a set of coupled nonlinear equations, whose dynamics is related to the Lorenz equations. For Prandtl and Rayleigh numbers in the range 10(-2) < or = Pr < or = 10(3) and 10(7) < or = Ra < or = 10(12), the model has the following features: (i) chaotic reversals may be exhibited at Ra > or = 10(7); (ii) the Reynolds number based on the root mean square velocity scales as Re(rms) approximately Ra([0.41...0.47]) (depending on Pr), and as Re(rms) approximately Pr(-[0.66...0.76]) (depending on Ra); and (iii) the mean reversal frequency follows an effective scaling law omega/(nu L(-2)) approximately Pr(-(0.64 +/- 0.01))Ra(0.44 +/- 0.01). The phase diagram of the model is sketched, and the observed transitions are discussed.
Collapse
Affiliation(s)
- Francisco Fontenele Araujo
- Department of Applied Physics and J. M. Burgers Centre for Fluid Dynamics, University of Twente, 7500 AE Enschede, The Netherlands
| | | | | |
Collapse
|
13
|
Furukawa A, Meyer H, Onuki A. Numerical simulation studies of the convective instability onset in a supercritical fluid. PHYSICAL REVIEW. E, STATISTICAL, NONLINEAR, AND SOFT MATTER PHYSICS 2005; 71:067301. [PMID: 16089918 DOI: 10.1103/physreve.71.067301] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 11/18/2004] [Indexed: 05/03/2023]
Abstract
Numerical simulation studies in 2D with the addition of noise are reported for the convection of a supercritical fluid, 3He , in a Rayleigh-Bénard cell. The noise addition is to accelerate the instability growth after starting the heat flow across the fluid, so as to bring simulations into better agreement with published experimental observations. Homogeneous temperature noise and spatial longitudinal periodic temperature variations in either top or bottom plates were programmed into the simulations. The second method was the most effective in speeding up the instability onset. For a small amplitude of the longitudinal perturbations, a semiquantitative agreement with the observations was obtained. The results are discussed in relation to predictions by El Khouri and Carlès.
Collapse
Affiliation(s)
- A Furukawa
- Department of Physics, Kyoto University, Kyoto 606-8502, Japan
| | | | | |
Collapse
|
14
|
Maekawa T, Ishii K, Shiroishi Y, Azuma H. Onset of buoyancy convection in a horizontal layer of a supercritical fluid heated from below. ACTA ACUST UNITED AC 2004. [DOI: 10.1088/0305-4470/37/32/009] [Citation(s) in RCA: 17] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/11/2022]
|
15
|
Furukawa A, Meyer H, Onuki A, Kogan AB. Convection in a very compressible fluid: comparison of simulations with experiments. PHYSICAL REVIEW. E, STATISTICAL, NONLINEAR, AND SOFT MATTER PHYSICS 2003; 68:056309. [PMID: 14682886 DOI: 10.1103/physreve.68.056309] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 05/21/2003] [Indexed: 05/24/2023]
Abstract
The time profile DeltaT(t) of the temperature difference, measured across a very compressible fluid layer of supercritical 3He after the start of a heat flow, shows a damped oscillatory behavior before steady-state convection is reached. The results for DeltaT(t) obtained from numerical simulations and from laboratory experiments are compared over a temperature range where the compressibility varies by a factor of approximately 40. First the steady-state convective heat current j(conv) as a function of the Rayleigh number Ra is presented, and the agreement is found to be good. Second, the shape of the time profile and two characteristic times in the transient part of DeltaT(t) from simulations and experiments are compared, namely (1) t(osc), the oscillatory period, and (2) t(p), the time of the first peak after starting the heat flow. These times, scaled by the diffusive time tau(D) versus Ra, are presented. The agreement is good for t(osc)/tau(D), where the results collapse on a single curve showing a power-law behavior. The simulation hence confirms the universal scaling behavior found experimentally. However for t(p)/tau(D), where the experimental data also collapse on a single curve, the simulation results show systematic departures from such a behavior. A possible reason for some of the disagreements, both in the time profile and in t(p), is discussed.
Collapse
Affiliation(s)
- A Furukawa
- Department of Physics, Kyoto University, Kyoto 606-8502, Japan
| | | | | | | |
Collapse
|
16
|
Soboleva EB. Adiabatic heating and convection caused by a fixed-heat-flux source in a near-critical fluid. PHYSICAL REVIEW. E, STATISTICAL, NONLINEAR, AND SOFT MATTER PHYSICS 2003; 68:042201. [PMID: 14682982 DOI: 10.1103/physreve.68.042201] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 08/01/2002] [Revised: 07/11/2003] [Indexed: 05/24/2023]
Abstract
Dynamics and heat transfer in a near-critical fluid in a square cavity with a finite heat source located at the bottom are studied numerically. A thermally insulated enclosure and a fixed-heat-flux source are considered. The two-dimensional simulation is based on the full Navier-Stokes equations with two-scale splitting of the pressure and the van der Waals equation of state. It is shown that the piston effect is independent of convection. Near the critical point, this effect becomes independent of criticality and convective motions are damped.
Collapse
Affiliation(s)
- E B Soboleva
- Institute for Problems in Mechanics, Russian Academy of Sciences, Prospect Vernadskogo, 101, b.1, 119526 Moscow, Russia.
| |
Collapse
|
17
|
Nikolayev VS, Dejoan A, Garrabos Y, Beysens D. Fast heat transfer calculations in supercritical fluids versus hydrodynamic approach. PHYSICAL REVIEW. E, STATISTICAL, NONLINEAR, AND SOFT MATTER PHYSICS 2003; 67:061202. [PMID: 16241213 DOI: 10.1103/physreve.67.061202] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 10/29/2002] [Indexed: 05/04/2023]
Abstract
This study investigates the heat transfer in a simple pure fluid whose temperature is slightly above its critical temperature. We propose an efficient numerical method to predict the heat transfer in such fluids when the gravity can be neglected. The method, based on a simplified thermodynamic approach, is compared with direct numerical simulations of the Navier-Stokes and energy equations performed for CO2 and SF6. A realistic equation of state is used to describe both fluids. The proposed method agrees with the full hydrodynamic solution and provides a huge gain in computation time. The connection between the purely thermodynamic and hydrodynamic descriptions is also discussed.
Collapse
Affiliation(s)
- V S Nikolayev
- DSM-DRFMC-Service des Basses Températures, CEA Grenoble, France.
| | | | | | | |
Collapse
|
18
|
Meyer H, Kogan AB. Onset of convection in a very compressible fluid: the transient toward steady state. PHYSICAL REVIEW. E, STATISTICAL, NONLINEAR, AND SOFT MATTER PHYSICS 2002; 66:056310. [PMID: 12513599 DOI: 10.1103/physreve.66.056310] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 02/25/2002] [Revised: 07/11/2002] [Indexed: 05/24/2023]
Abstract
We analyze the time profile deltaT(t) of the temperature difference, measured across a very compressible supercritical 3He fluid layer in its convective state. The experiments were done along the critical isochore in a Rayleigh-Bénard cell after starting the vertical constant heat flow q. For q sufficiently well above that needed for the convection onset, the transient deltaT(t) for a given epsilon identical with (T-T(c))/T(c), with T(c)=3.318 K, shows a damped oscillatory profile with period t(osc) modulating a smooth base profile. The smooth profile forms the exponential tail of the transient which tends to the steady-state deltaT( infinity ) with a time constant tau(tail). The scaled times t(osc)/t(D) and tau(tail)/t(D) from all the data could be collapsed onto two curves as a function of the Rayleigh number over approximately 3.5 decades. Here t(D) is the characteristic thermal diffusion time. Furthermore, comparisons are made between measurements of a third characteristic time t(m) between the first peak and the first minimum in the deltaT(t) profile and its estimation by Onuki et al. Also, comparisons are made between the observed oscillations and the two-dimensional simulations by Onuki et al. and by Amiroudine and Zappoli. For epsilon<9 x 10(-3), the experiments show a crossover to a different transient regime. This regime, which we briefly describe, is not understood at present.
Collapse
Affiliation(s)
- Horst Meyer
- Department of Physics, Duke University, Durham, North Carolina 27708-0305, USA
| | | |
Collapse
|