A numerical study of turbulent Rayleigh-Bénard convection

Il flusso di Rayleigh-Bénard, costituito da un fluido racchiuso fra due pareti a diversa temperatura, rappresenta il paradigma della convezione termica. In natura e nelle applicazioni industriali, il moto convettivo avviene principalmente in regime turbolento, rivelando un fenomeno estremamente complesso. L'obiettivo principale di questo elaborato di tesi consiste nell'isolare e descrivere gli aspetti salienti di un flusso turbolento di Rayleigh-Bénard. L'analisi è applicata a dati ottenuti da tre simulazioni numeriche dirette effettuate allo stesso numero di Rayleigh (Ra=10^5) e a numeri di Prandtl differenti (Pr=0.7,2,7). Sulla base di alcune statistiche a singolo punto, vengono definite nel flusso tre regioni caratteritiche: il bulk al centro della cella, lo strato limite termico e quello viscoso in prossimità delle pareti. Grazie all'analisi dei campi istantanei e delle correlazioni spaziali a due punti, sono state poi individuate due strutture fondamentali della convezione turbolenta: le piume termiche e la circolazione a grande scala. L'equazione generalizzata di Kolmogorov, introdotta nell'ultima parte della trattazione, permette di approcciare il problema nella sua complessità, visualizzando come l'energia cinetica viene immessa, si distribuisce e viene dissipata sia nello spazio fisico, sia in quello delle scale turbolente. L'immagine che emerge dall'analisi complessiva è quella di un flusso del tutto simile a una macchina termica. L'energia cinetica viene prodotta nel bulk, considerato il motore del flusso, e da qui fluisce verso le pareti, dove viene infine dissipata.

Range of validity of weakly non-linear theory in Rayleigh-Bénard convection

In this paper we examine the equilibrium states of periodic finite amplitude flow in a horizontal channel with differential heating between the two rigid boundaries. The solutions to the Navier-Stokes equations are obtained by means of a perturbation method for evaluating the Landau coefficients and through a Newton-Raphson iterative method that results from the Fourier expansion of the solutions that bifurcate above the linear stability threshold of infini- tesimal disturbances. The results obtained from these two different methods of evaluating the convective flow are compared in the neighbourhood of the critical Rayleigh number. We find that for small Prandtl numbers the discrepancy of the two methods is noticeable.

Dynamic states in rotating Rayleigh-Bénard convection systems

A new geometry-independent state - a traveling-wave wall state - is proposed as the mechanism whereby which the experimentally observed wall-localized states in rotating Rayleigh-Bénard convection systems preempt the bulk state at large rotation rates. Its properties are calculated for the illustrative case of free-slip top and bottom boundary conditions. At small rotation rates, this new wall state is found to disappear. A detailed study of the dynamics of the wall state and the bulk state in the transition region where this disappearance occurs is conducted using a Swift-Hohenberg model system. The Swift-Hohenberg model, with appropriate reflection-symmetry- breaking boundary conditions, is also shown to exhibit traveling-wave wall states, further demonstrating that traveling-wave wall states are a generic feature of nonequilibrium pattern-forming systems. A numerical code for the Swift-Hohenberg model in an annular geometry was written and used to investigate the dynamics of rotating Rayleigh-Bénard convection systems.

Large Rayleigh number thermal convection: heat flux predictions and strongly nonlinear solutions

We investigate the structure of strongly nonlinear Rayleigh–Bénard convection cells in the asymptotic limit of large Rayleigh number and fixed, moderate Prandtl number. Unlike the flows analyzed in prior theoretical studies of infinite Prandtl number convection, our cellular solutions exhibit dynamically inviscid constant-vorticity cores. By solving an integral equation for the cell-edge temperature distribution, we are able to predict, as a function of cell aspect ratio, the value of the core vorticity, details of the flow within the thin boundary layers and rising/falling plumes adjacent to the edges of the convection cell, and, in particular, the bulk heat flux through the layer. The results of our asymptotic analysis are corroborated using full pseudospectral numerical simulations and confirm that the heat flux is maximized for convection cells that are roughly square in cross section.

Dynamos driven by weak thermal convection and heterogeneous outer boundary heat flux

We use numerical dynamo models with heterogeneous core-mantle boundary (CMB) heat flux to show that lower mantle lateral thermal variability may help support a dynamo under weak thermal convection. In our reference models with homogeneous CMB heat flux, convection is either marginally supercritical or absent, always below the threshold for dynamo onset. We find that lateral CMB heat flux variations organize the flow in the core into patterns that favour the growth of an early magnetic field. Heat flux patterns symmetric about the equator produce non-reversing magnetic fields, whereas anti-symmetric patterns produce polarity reversals. Our results may explain the existence of the geodynamo prior to inner core nucleation under a tight energy budget. Furthermore, in order to sustain a strong geomagnetic field, the lower mantle thermal distribution was likely dominantly symmetric about the equator. (C) 2015 Elsevier B.V. All rights reserved.

Thermal Convection in Ice Sheets: We Look But Do Not See

Thermal convection in the Antarctic and Greenland ice sheets has been dismissed on the grounds that radio-echo stratigraphy is undisturbed for long distances. However, the undisturbed stratigraphy lies, for the most part, above the density inversion in polar ice sheets and therefore does not disprove convection. An echo-free zone is widespread below the density inversion, yet nobody has cited this as a strong indication that convection is indeed present at d�pth. A generalized Rayleigh criterion for thermal convection in e1astic-viscoplastic polycrystalline solids heated from below is developed and applied to ice-sheet convection. An infinite Rayleigh number at the onset of primary creep decreases with time and becomes constant when secondary creep dominates, suggesting that any thermal buoyancy stress can initiate convection but convection cannot be sustained below a buoyancy stress of about 3 kPa. An analysis of the temperature profile down the Byrd Station core hole suggests that about 1000 m of ice below the density inversion will sustain convection. Creep along the Byrd Station strain network, radar sounding in East Antarctica, and seismic sounding in West Antarctica are examined for evidence of convective creep superimposed on advective creep. It is concluded that the evidence for convection is there, if we look for it with the intention offinding it.

Tropical deep convection, convective available potential energy and sea surface temperature

The convective available potential energy (CAFE) based on monthly mean sounding has been shown to be relevant to deep convection in the tropics. The variation of CAFE with SST has been found to be similar to the variation of the frequency of deep convection at one station each in the tropical Atlantic and W. Pacific oceans. This suggests a strong link between the frequency of tropical convection and CAFE. It has been shown that CAFE so derived can be interpreted as the work potential of the atmosphere above the boundary layer with ascent in the convective region and subsidence in the surrounding cloud-free region.

Experimental investigation on the chaotic phenomena in the wake of a natural thermal convection flow

Chaotic phenomena in the wake of thermal convection flow fields above a heating flat plate were investigated experimentally. A newly developed electron beam fluorescence technique (EBF) was used to simultaneously measure density fluctuation at 7 points in a cross section above the plate. Correlation dimensions, intermittence coefficients, Fourier spectrum have been obtained for different Grashof numbers. Spatial distribution of correlation dimensions are presented. The experimental result shows that there is a certain relationship between the density fluctuation and the Gr number. And time-spacial characteristic of chaos evolution is also given.

Rayleigh-Bénard热对流的动理论分析

利用基于分子动理论的直接模拟Monte Carlo(DSMC)方法,研究了 Rayleigh-Bènard问题.计算中,上下平板表面温度之比固定为0.1.Kn=0.01时,随着Ra数的增大,大约在1700附近,流动从热传导状态转变为热对流状态,DSMC计算得到的下平板热流与Ra数的关系与经典实验和理论结果相符.Kn=0.05时,流动保持稳定的热传导状态,Ra数的增大并不能引发热对流现象.

Estimation thermal expansion coefficient from lattice energy for inorganic crystals

By using the study of the lattice energy and the structural parameters of binary inorganic crystals, a new parameter reflecting the thermal expansion property has been found, the relation between the linear expansion coefficient and new parameter has been established. A semiempirical method for evaluation of linear expansion coefficient from the lattice energy is presented, and developed to the complex crystals. The estimated values of the linear expansion coefficients of both simple and complex crystals are in good agreement with the experimental values.

Two thermal methods to measure the energy fluence of a brief exposure of diagnostic x rays

This paper describes two simple thermal methods for measuring the energy fluence in J/cm 2 from a diagnostic x-ray exposure. Both detectors absorb essentially 100% of the radiation and give a signal that is directly proportional to the energy fluence of the x-ray beam. One detector measures the thermal effect when a pulse of x rays is totally absorbed in the pyroelectric detector of lead-zirconium-titanate (PZT). The other detector measures the expansion of a gas surrounding a lead disk detector in a photoacoustic chamber. The increased pressure of the gas is transmitted through a 1-mm duct to a sensitive microphone. Both detectors have previously been used to measure the energy fluence rate of continuous x-ray beams in the same energy region using a chopped beam and a lock-in amplifier. Measurement of the energy fluence of a pulse of radiation eliminates the need for the beam chopper and lock-in amplifier and results in a simple, rugged, and inexpensive dosimeter. Either method can be combined with the area of the beam to give an estimate of the imparted energy to the patient from a diagnostic x-ray exposure.

Numerical coupling of thermal-electric network models and energy-transport equations including optoelectronic semiconductor devices

In this work the numerical coupling of thermal and electric network models with model equations for optoelectronic semiconductor devices is presented. Modified nodal analysis (MNA) is applied to model electric networks. Thermal effects are modeled by an accompanying thermal network. Semiconductor devices are modeled by the energy-transport model, that allows for thermal effects. The energy-transport model is expandend to a model for optoelectronic semiconductor devices. The temperature of the crystal lattice of the semiconductor devices is modeled by the heat flow eqaution. The corresponding heat source term is derived under thermodynamical and phenomenological considerations of energy fluxes. The energy-transport model is coupled directly into the network equations and the heat flow equation for the lattice temperature is coupled directly into the accompanying thermal network. The coupled thermal-electric network-device model results in a system of partial differential-algebraic equations (PDAE). Numerical examples are presented for the coupling of network- and one-dimensional semiconductor equations. Hybridized mixed finite elements are applied for the space discretization of the semiconductor equations. Backward difference formluas are applied for time discretization. Thus, positivity of charge carrier densities and continuity of the current density is guaranteed even for the coupled model.