Fast Numerical simulations of 2D turbulence using a dynamic model for Subgrid Motions

We present numerical simulation of 2D turbulent flow using a new model for the subgrid scales which are computed using a dynamic equation linking the subgrid scales with the resolved velocity. This equation is not postulated, but derived from the constitutive equations under the assumption that the non-linear interactions of subgrid scales between themselves are equivalent to a turbulent viscosity.The performances of our model are compared with Direct Numerical Simulations of decaying and forced turbulence. For a same resolution, numerical simulations using our model allow for a significant reduction of the computational time (of the order of 100 in the case we consider), and allow the achievement of significantly larger Reynolds number than the direct method.Comment: 35 pages, 9 figure

The energy budget in Rayleigh-Benard convection

It is shown using three series of Rayleigh number simulations of varying aspect ratio AR and Prandtl number Pr that the normalized dissipation at the wall, while significantly greater than 1, approaches a constant dependent upon AR and Pr. It is also found that the peak velocity, not the mean square velocity, obeys the experimental scaling of Ra^{0.5}. The scaling of the mean square velocity is closer to Ra^{0.46}, which is shown to be consistent with experimental measurements and the numerical results for the scaling of Nu and the temperature if there are strong correlations between the velocity and temperature.Comment: 5 pages, 3 figures, new version 13 Mar, 200

Forced Stratified Turbulence: Successive Transitions with Reynolds Number

Numerical simulations are made for forced turbulence at a sequence of increasing values of Reynolds number, R, keeping fixed a strongly stable, volume-mean density stratification. At smaller values of R, the turbulent velocity is mainly horizontal, and the momentum balance is approximately cyclostrophic and hydrostatic. This is a regime dominated by so-called pancake vortices, with only a weak excitation of internal gravity waves and large values of the local Richardson number, Ri, everywhere. At higher values of R there are successive transitions to (a) overturning motions with local reversals in the density stratification and small or negative values of Ri; (b) growth of a horizontally uniform vertical shear flow component; and (c) growth of a large-scale vertical flow component. Throughout these transitions, pancake vortices continue to dominate the large-scale part of the turbulence, and the gravity wave component remains weak except at small scales.Comment: 8 pages, 5 figures (submitted to Phys. Rev. E

Calibrating the Mixing Length Parameter for a Red Giant Envelope

Two-dimensional hydrodynamical simulations were made to calibrate the mixing length parameter for modeling red giant's convective envelope. As was briefly reported in Asida & Tuchman (97), a comparison of simulations starting with models integrated with different values of the mixing length parameter, has been made. In this paper more results are presented, including tests of the spatial resolution and Large Eddy Simulation terms used by the numerical code. The consistent value of the mixing length parameter was found to be 1.4, for a red giant of mass 1.2 solar-mass, core mass of 0.96 solar-mass, luminosity of 200 solar-luminosities, and metallicity Z=0.001.Comment: 18 pages, 1 table, 13 figures. Accepted for publication in Ap.

Stability and turbulent transport in Taylor–Couette flow from analysis of experimental data

International audienceThis paper provides a prescription for the turbulent viscosity in rotating shear flows for use e.g. in geophysical and astrophysical contexts. This prescription is the result of the detailed analysis of the experimental data obtained in several studies of the transition to turbulence and turbulent transport in Taylor-Couette flow. We first introduce a new set of control parameters, based on dynamical rather than geometrical considerations, so that the analysis applies more naturally to rotating shear flows in general and not only to Taylor-Couette flow. We then investigate the transition thresholds in the supercritical and the subcritical regime in order to extract their general dependencies on the control parameters. The inspection of the mean profiles provides us with some general hints on the mean to laminar shear ratio. Then the examination of the torque data allows us to propose a decomposition of the torque dependence on the control parameters in two terms, one completely given by measurements in the case where the outer cylinder is at rest, the other one being a universal function provided here from experimental fits. As a result, we obtain a general expression for the turbulent viscosity and compare it to existing prescription in the literature. Finally, throughout all the paper we discuss the influence of additional effects such as stratification or magnetic fields

Global effects of local sound-speed perturbations in the Sun: A theoretical study

We study the effect of localized sound-speed perturbations on global mode frequencies by applying techniques of global helioseismology on numerical simulations of the solar acoustic wave field. Extending the method of realization noise subtraction (e.g. Hanasoge et al. 2007) to global modes and exploiting the luxury of full spherical coverage, we are able to achieve very highly resolved frequency differences that are used to study sensitivities and the signatures of the thermal asphericities. We find that (1) global modes are almost twice as sensitive to sound-speed perturbations at the bottom of the convection zone as in comparison to anomalies well in the radiative interior ($r\lesssim0.55 R_\odot$), (2) the $m$-degeneracy is lifted ever so slightly, as seen in the $a$ coefficients, and (3) modes that propagate in the vicinity of the perturbations show small amplitude shifts ($\sim 0.5$).Comment: Submitted to Solar Physic

Scientific drilling of Lake Chalco, Basin of Mexico (MexiDrill)

The primary scientific objective of MexiDrill, the Basin of Mexico Drilling Program, is development of a continuous, high-resolution ∼400 kyr lacustrine record of tropical North American environmental change. The field location, in the densely populated, water-stressed Mexico City region gives this record particular societal relevance. A detailed paleoclimate reconstruction from central Mexico will enhance our understanding of long-term natural climate variability in the North American tropics and its relationship with changes at higher latitudes. The site lies at the northern margin of the Intertropical Convergence Zone (ITCZ), where modern precipitation amounts are influenced by sea surface temperatures in the Pacific and Atlantic basins. During the Last Glacial Maximum (LGM), more winter precipitation at the site is hypothesized to have been a consequence of a southward displacement of the mid-latitude westerlies. It thus represents a key spatial node for understanding large-scale hydrological variability of tropical and subtropical North America and is at an altitude (2240 m a.s.l.), typical of much of western North America. In addition, its sediments contain a rich record of pre-Holocene volcanic history; knowledge of the magnitude and frequency relationships of the area's explosive volcanic eruptions will improve capacity for risk assessment of future activity. Explosive eruption deposits will also be used to provide the backbone of a robust chronology necessary for full exploitation of the paleoclimate record. Here we report initial results from, and outreach activities of, the 2016 coring campaign

Use of 16S rRNA Gene Based Clone Libraries to Assess Microbial Communities Potentially Involved in Anaerobic Methane Oxidation in a Mediterranean Cold Seep

This study provides data on the diversities of bacterial and archaeal communities in an active methane seep at the Kazan mud volcano in the deep Eastern Mediterranean sea. Layers of varying depths in the Kazan sediments were investigated in terms of (1) chemical parameters and (2) DNA-based microbial population structures. The latter was accomplished by analyzing the sequences of directly amplified 16S rRNA genes, resulting in the phylogenetic analysis of the prokaryotic communities. Sequences of organisms potentially associated with processes such as anaerobic methane oxidation and sulfate reduction were thus identified. Overall, the sediment layers revealed the presence of sequences of quite diverse bacterial and archaeal communities, which varied considerably with depth. Dominant types revealed in these communities are known as key organisms involved in the following processes: (1) anaerobic methane oxidation and sulfate reduction, (2) sulfide oxidation, and (3) a range of (aerobic) heterotrophic processes. In the communities in the lowest sediment layer sampled (22–34 cm), sulfate-reducing bacteria and archaea of the ANME-2 cluster (likely involved in anaerobic methane oxidation) were prevalent, whereas heterotrophic organisms abounded in the top sediment layer (0–6 cm). Communities in the middle layer (6–22 cm) contained organisms that could be linked to either of the aforementioned processes. We discuss how these phylogeny (sequence)-based findings can support the ongoing molecular work aimed at unraveling both the functioning and the functional diversities of the communities under study

Relationship between optical coherence tomography sector peripapillary angioflow-density and Octopus visual field cluster mean defect values

PURPOSE: To compare the relationship of Octopus perimeter cluster mean-defect (cluster MD) values with the spatially corresponding optical coherence tomography (OCT) sector peripapillary angioflow vessel-density (PAFD) and sector retinal nerve fiber layer thickness (RNFLT) values. METHODS: High quality PAFD and RNFLT images acquired on the same day with the Angiovue/RTVue-XR Avanti OCT (Optovue Inc., Fremont, USA) on 1 eye of 27 stable early-to-moderate glaucoma, 22 medically controlled ocular hypertensive and 13 healthy participants were analyzed. Octopus G2 normal visual field test was made within 3 months from the imaging. RESULTS: Total peripapillary PAFD and RNFLT showed similar strong positive correlation with global mean sensitivity (r-values: 0.6710 and 0.6088, P<0.0001), and similar (P = 0.9614) strong negative correlation (r-values: -0.4462 and -0.4412, P</=0.004) with global MD. Both inferotemporal and superotemporal sector PAFD were significantly (</=0.039) lower in glaucoma than in the other groups. No significant difference between the corresponding inferotemporal and superotemporal parameters was seen. The coefficient of determination (R2) calculated for the relationship between inferotemporal sector PAFD and superotemporal cluster MD (0.5141, P<0.0001) was significantly greater than that between inferotemporal sector RNFLT and superotemporal cluster MD (0.2546, P = 0.0001). The R2 values calculated for the relationships between superotemporal sector PAFD and RNFLT, and inferotemporal cluster MD were similar (0.3747 and 0.4037, respectively, P<0.0001). CONCLUSION: In the current population the relationship between inferotemporal sector PAFD and superotemporal cluster MD was strong. It was stronger than that between inferotemporal sector RNFLT and superotemporal cluster MD. Further investigations are necessary to clarify if our results are valid for other populations and can be usefully applied for glaucoma research

On the Long Lasting “C‐Type” Structures in the Sodium Lidargram: The Lifetime of Kelvin‐Helmholtz Billows in the Mesosphere and Lower Thermosphere Region

In order to understand the characteristics of long‐lasting “C‐type” structure in the Sodium (Na) lidargram, six cases from different observational locations have been analyzed. The Na lidargram, collected from low‐, middle‐, and high‐latitude sites, show long lifetime of the C‐type structures which is believed to be the manifestation of Kelvin‐Helmholtz (KH) billows in the Mesosphere and Lower Thermosphere (MLT) region. In order to explore the characteristics of the long‐lasting C‐type structures, the altitude profile of square of Brunt‐Väisälä frequency in the MLT region has been derived using the temperature profile collected from the Na lidar instruments and the SABER instrument onboard TIMED satellite. It is found to be positive in the C‐type structure region for all the six cases which indicates that the regions are convectively stable. Simultaneous wind measurements, which allowed us to calculate the Richardson numbers and Reynolds numbers for three cases, suggest that the regions where the C‐type structure appeared were dynamically stable and nonturbulent. This paper brings out a hypothesis wherein the low temperature can increase the magnitude of the Prandtl number and convectively stable atmospheric region can cause the magnitude of Reynolds number to decrease. As a consequence, the remnant of previously generated KH billows in nearly “frozen‐in” condition can be advected through this conducive region to a different location by the background wind where they can sustain for a long time without much deformation. These long‐lived KH billows in the MLT region will eventually manifest the long‐lasting C‐type structures in the Na lidargram