Gliese 581g as a scaled-up version of Earth: atmospheric circulation simulations

We use three-dimensional simulations to study the atmospheric circulation on the first Earth-sized exoplanet discovered in the habitable zone of an M star. We treat Gliese 581g as a scaled-up version of Earth by considering increased values for the exoplanetary radius and surface gravity, while retaining terrestrial values for parameters which are unconstrained by current observations. We examine the long-term, global temperature and wind maps near the surface of the exoplanet --- the climate. The specific locations for habitability on Gliese 581g depend on whether the exoplanet is tidally-locked and how fast radiative cooling occurs on a global scale. Independent of whether the existence of Gliese 581g is confirmed, our study highlights the use of general circulation models to quantify the atmospheric circulation on potentially habitable, Earth-sized exoplanets, which will be the prime targets of exoplanet discovery and characterization campaigns in the next decade.Comment: Accepted by MNRAS. 15 pages, 13 figures. Sample movies of simulations are available at http://www.phys.ethz.ch/~kheng/fms

Large-Eddy simulation of pulsatile blood flow

Large-Eddy simulation (LES) is performed to study pulsatile blood flow through a 3D model of arterial stenosis. The model is chosen as a simple channel with a biological type stenosis formed on the top wall. A sinusoidal non-additive type pulsation is assumed at the inlet of the model to generate time dependent oscillating flow in the channel and the Reynolds number of 1200, based on the channel height and the bulk velocity, is chosen in the simulations. We investigate in detail the transition-to-turbulent phenomena of the non-additive pulsatile blood flow downstream of the stenosis. Results show that the high level of flow recirculation associated with complex patterns of transient blood flow have a significant contribution to the generation of the turbulent fluctuations found in the post-stenosis region. The importance of using LES in modelling pulsatile blood flow is also assessed in the paper through the prediction of its sub-grid scale contributions. In addition, some important results of the flow physics are achieved from the simulations, these are presented in the paper in terms of blood flow velocity, pressure distribution, vortices, shear stress, turbulent fluctuations and energy spectra, along with their importance to the relevant medical pathophysiology

Workshop on Drought Forecasting for Northeast Brazil

Precipitation forecasting parameters for northeast Brazil were developed. Hydrological, sociological, and economic aspects were examined. A drought forecasting model is presented

Large eddy simulation of two-dimensional isotropic turbulence

Large eddy simulation (LES) of forced, homogeneous, isotropic, two-dimensional (2D) turbulence in the energy transfer subrange is the subject of this paper. A difficulty specific to this LES and its subgrid scale (SGS) representation is in that the energy source resides in high wave number modes excluded in simulations. Therefore, the SGS scheme in this case should assume the function of the energy source. In addition, the controversial requirements to ensure direct enstrophy transfer and inverse energy transfer make the conventional scheme of positive and dissipative eddy viscosity inapplicable to 2D turbulence. It is shown that these requirements can be reconciled by utilizing a two-parametric viscosity introduced by Kraichnan (1976) that accounts for the energy and enstrophy exchange between the resolved and subgrid scale modes in a way consistent with the dynamics of 2D turbulence; it is negative on large scales, positive on small scales and complies with the basic conservation laws for energy and enstrophy. Different implementations of the two-parametric viscosity for LES of 2D turbulence were considered. It was found that if kept constant, this viscosity results in unstable numerical scheme. Therefore, another scheme was advanced in which the two-parametric viscosity depends on the flow field. In addition, to extend simulations beyond the limits imposed by the finiteness of computational domain, a large scale drag was introduced. The resulting LES exhibited remarkable and fast convergence to the solution obtained in the preceding direct numerical simulations (DNS) by Chekhlov et al. (1994) while the flow parameters were in good agreement with their DNS counterparts. Also, good agreement with the Kolmogorov theory was found. This LES could be continued virtually indefinitely. Then, a simplifiedComment: 34 pages plain tex + 18 postscript figures separately, uses auxilary djnlx.tex fil

LES of non-Newtonian physiological blood flow in a model of arterial stenosis

Large Eddy Simulation (LES) is performed to study the physiological pulsatile transition-to-turbulent non-Newtonian blood flow through a 3D model of arterial stenosis by using five different blood viscosity models: (i) Power-law, (ii) Carreau, (iii) Quemada, (iv) Cross and (v) modified-Casson. The computational domain has been chosen is a simple channel with a biological type stenosis formed eccentrically on the top wall. The physiological pulsation is generated at the inlet of the model using the first four harmonic series of the physiological pressure pulse (Loudon and Tordesillas [1]). The effects of the various viscosity models are investigated in terms of the global maximum shear rate, post-stenotic re-circulation zone, mean shear stress, mean pressure, and turbulent kinetic energy. We find that the non-Newtonian viscosity models enlarge the length of the post-stenotic re-circulation region by moving the reattachment point of the shear layer separating from the upper wall further downstream. But the turbulent kinetic energy at the immediate post-lip of the stenosis drops due to the effects of the non-Newtonian viscosity. The importance of using LES in modelling the non-Newtonian physiological pulsatile blood flow is also assessed for the different viscosity models in terms of the results of the dynamic subgrid-scale (SGS) stress Smagorinsky model constant, C<sub>s</sub>, and the corresponding SGS normalised viscosity

Atmospheric circulation of tidally locked exoplanets: a suite of benchmark tests for dynamical solvers

The complexity of atmospheric modelling and its inherent non-linearity, together with the limited amount of data of exoplanets available, motivate model intercomparisons and benchmark tests. In the geophysical community, the Held-Suarez test is a standard benchmark for comparing dynamical core simulations of the Earth's atmosphere with different solvers, based on statistically-averaged flow quantities. In the present study, we perform analogues of the Held-Suarez test for tidally-locked exoplanets with the GFDL-Princeton Flexible Modeling System (FMS) by subjecting both the spectral and finite difference dynamical cores to a suite of tests, including the standard benchmark for Earth, a hypothetical tidally-locked Earth, a "shallow" hot Jupiter model and a "deep" model of HD 209458b. We find qualitative and quantitative agreement between the solvers for the Earth, tidally-locked Earth and shallow hot Jupiter benchmarks, but the agreement is less than satisfactory for the deep model of HD 209458b. Further investigation reveals that closer agreement may be attained by arbitrarily adjusting the values of the horizontal dissipation parameters in the two solvers, but it remains the case that the magnitude of the horizontal dissipation is not easily specified from first principles. Irrespective of radiative transfer or chemical composition considerations, our study points to limitations in our ability to accurately model hot Jupiter atmospheres with meteorological solvers at the level of ten percent for the temperature field and several tens of percent for the velocity field. Direct wind measurements should thus be particularly constraining for the models. Our suite of benchmark tests also provides a reference point for researchers wishing to adapt their codes to study the atmospheric circulation regimes of tidally-locked Earths/Neptunes/Jupiters.Comment: Accepted by MNRAS, 23 pages, 17 figures, 2 tables. No changes from previous version, except MNRAS wants no hyphen in the title. Sample movies of simulations are available at http://www.phys.ethz.ch/~kheng/fms

Large-eddy simulation of the lid-driven cubic cavity flow by the spectral element method

This paper presents the large-eddy simulation of the lid-driven cubic cavity flow by the spectral element method (SEM) using the dynamic model. Two spectral filtering techniques suitable for these simulations have been implemented. Numerical results for Reynolds number $\text{Re}=12'000$ are showing very good agreement with other experimental and DNS results found in the literature

Weather, climate and the nature of predictability

The prediction and simulation of future weather and climate is a key ingredient in good weather risk management. This chapter briefly reviews the nature and underlying sources of predictability on timescales from hours-ahead to centuries-ahead. The traditional distinction between ‘weather’ and ‘climate’ predictions is described, and the role of recent scientific developments in driving a convergence of these two classic problems is highlighted. The chapter concludes by outlining and comparing the two main strategies used for creating weather and climate predictions, and discussing the challenges of using predictions in quantitative applications

Large Eddy Simulation of Turbulent Channel Flows by the Rational LES Model

The rational large eddy simulation (RLES) model is applied to turbulent channel flows. This approximate deconvolution model is based on a rational (subdiagonal Pade') approximation of the Fourier transform of the Gaussian filter and is proposed as an alternative to the gradient (also known as the nonlinear or tensor-diffusivity) model. We used a spectral element code to perform large eddy simulations of incompressible channel flows at Reynolds numbers based on the friction velocity and the channel half-width Re{sub tau} = 180 and Re{sub tau} = 395. We compared the RLES model with the gradient model. The RLES results showed a clear improvement over those corresponding to the gradient model, comparing well with the fine direct numerical simulation. For comparison, we also present results corresponding to a classical subgrid-scale eddy-viscosity model such as the standard Smagorinsky model.Comment: 31 pages including 15 figure

Numerical studies towards practical large-eddy simulation

Large-eddy simulation developments and validations are presented for an improved simulation of turbulent internal flows. Numerical methods are proposed according to two competing criteria: numerical qualities (precision and spectral characteristics), and adaptability to complex configurations. First, methods are tested on academic test-cases, in order to abridge with fundamental studies. Consistent results are obtained using adaptable finite volume method, with higher order advection fluxes, implicit grid filtering and "low-cost" shear-improved Smagorinsky model. This analysis particularly focuses on mean flow, fluctuations, two-point correlations and spectra. Moreover, it is shown that exponential averaging is a promising tool for LES implementation in complex geometry with deterministic unsteadiness. Finally, adaptability of the method is demonstrated by application to a configuration representative of blade-tip clearance flow in a turbomachine