Numerical dissipation vs. subgrid-scale modelling for large eddy simulation

This study presents an alternative way to perform large eddy simulation based on a targeted numerical dissipation introduced by the discretization of the viscous term. It is shown that this regularisation technique is equivalent to the use of spectral vanishing viscosity. The flexibility of the method ensures high-order accuracy while controlling the level and spectral features of this purely numerical viscosity. A Pao-like spectral closure based on physical arguments is used to scale this numerical viscosity a priori. It is shown that this way of approaching large eddy simulation is more efficient and accurate than the use of the very popular Smagorinsky model in standard as well as in dynamic version. The main strength of being able to correctly calibrate numerical dissipation is the possibility to regularise the solution at the mesh scale. Thanks to this property, it is shown that the solution can be seen as numerically converged. Conversely, the two versions of the Smagorinsky model are found unable to ensure regularisation while showing a strong sensitivity to numerical errors. The originality of the present approach is that it can be viewed as implicit large eddy simulation, in the sense that the numerical error is the source of artificial dissipation, but also as explicit subgrid-scale modelling, because of the equivalence with spectral viscosity prescribed on a physical basis

Physical scaling of numerical dissipation for LES

In this work, we are interested in an alternative way to perform LES using a numerical substitute of a subgrid-scale model with a calibration based on physical inputs

Xcompact3D: An open-source framework for solving turbulence problems on a Cartesian mesh

Xcompact3D is a Fortran 90–95 open-source framework designed for fast and accurate simulations of turbulent flows, targeting CPU-based supercomputers. It is an evolution of the flow solver Incompact3D which was initially designed in France in the mid-90’s for serial processors to solve the incompressible Navier–Stokes equations. Incompact3D was then ported to parallel High Performance Computing (HPC) systems in the early 2010’s. Very recently the capabilities of Incompact3D have been extended so that it can now tackle more flow regimes (from incompressible flows to compressible flows at low Mach numbers), resulting in the design of a new user-friendly framework called Xcompact3D. The present manuscript presents an overview of Xcompact3D with a particular focus on its functionalities, its ready-to-run simulations and a few case studies to demonstrate its impact

One-Point Probability Distribution Functions of Supersonic Turbulent Flows in Self-Gravitating Media

Turbulence is essential for understanding the structure and dynamics of molecular clouds and star-forming regions. There is a need for adequate tools to describe and characterize the properties of turbulent flows. One-point probability distribution functions (pdf's) of dynamical variables have been suggested as appropriate statistical measures and applied to several observed molecular clouds. However, the interpretation of these data requires comparison with numerical simulations. To address this issue, SPH simulations of driven and decaying, supersonic, turbulent flows with and without self-gravity are presented. In addition, random Gaussian velocity fields are analyzed to estimate the influence of variance effects. To characterize the flow properties, the pdf's of the density, of the line-of-sight velocity centroids, and of the line centroid increments are studied. This is supplemented by a discussion of the dispersion and the kurtosis of the increment pdf's, as well as the spatial distribution of velocity increments for small spatial lags. From the comparison between different models of interstellar turbulence, it follows that the inclusion of self-gravity leads to better agreement with the observed pdf's in molecular clouds. The increment pdf's for small spatial lags become exponential for all considered velocities. However, all the processes considered here lead to non-Gaussian signatures, differences are only gradual, and the analyzed pdf's are in addition projection dependent. It appears therefore very difficult to distinguish between different physical processes on the basis of pdf's only, which limits their applicability for adequately characterizing interstellar turbulence.Comment: 38 pages (incl. 17 figures), accepted for publication in ApJ, also available with full resolution figures at http://www.strw.leidenuniv.nl/~klessen/Preprint

The association between body mass index and brain morphology in children: a population-based study

Brain morphology is altered in both anorexia nervosa and obesity. However, it is yet unclear if the relationship between Body Mass Index-Standard Deviation Score (BMI-SDS) and brain morphology exists across the BMI-SDS spectrum, or is present only in the extremes. The study involved 3160 9-to-11 year-old children (50.3% female) who participate in Generation R, a population-based study. Structural MRI scans were obtained from all children and FreeSurfer was used to quantify both global and surface-based measures of gyrification and cortical thickness. Body length and weight were measured to calculate BMI. Dutch growth curves were used to calculate BMI-SDS. BMI-SDS was analyzed continuously and in two categories (median split). The relationship between BMI-SDS (range − 3.82 to 3.31) and gyrification showed an inverted-U shape curve in children with both lower and higher BMI-SDS values having lower gyrification in widespread areas of the brain. BMI-SDS had a positive linear association with cortical thickness in multiple brain regions. This study provides evidence for an association between BMI-SDS and brain morphology in a large sample of children from the general population and suggests that a normal BMI during childhood is important for brain development. Future studies could determine whether lifestyle modifications optimize BMI-SDS result in return to more typical patterns of brain morphology

Genetic variants for head size share genes and pathways with cancer

The size of the human head is highly heritable, but genetic drivers of its variation within the general population remain unmapped. We perform a genome-wide association study on head size (N = 80,890) and identify 67 genetic loci, of which 50 are novel. Neuroimaging studies show that 17 variants affect specific brain areas, but most have widespread effects. Gene set enrichment is observed for various cancers and the p53, Wnt, and ErbB signaling pathways. Genes harboring lead variants are enriched for macrocephaly syndrome genes (37-fold) and high-fidelity cancer genes (9-fold), which is not seen for human height variants. Head size variants are also near genes preferentially expressed in intermediate progenitor cells, neural cells linked to evolutionary brain expansion. Our results indicate that genes regulating early brain and cranial growth incline to neoplasia later in life, irrespective of height. This warrants investigation of clinical implications of the link between head size and cancer.</p

Genetic Variants For Head Size Share Genes and Pathways With Cancer

The size of the human head is highly heritable, but genetic drivers of its variation within the general population remain unmapped. We perform a genome-wide association study on head size (N = 80,890) and identify 67 genetic loci, of which 50 are novel. Neuroimaging studies show that 17 variants affect specific brain areas, but most have widespread effects. Gene set enrichment is observed for various cancers and the p53, Wnt, and ErbB signaling pathways. Genes harboring lead variants are enriched for macrocephaly syndrome genes (37-fold) and high-fidelity cancer genes (9-fold), which is not seen for human height variants. Head size variants are also near genes preferentially expressed in intermediate progenitor cells, neural cells linked to evolutionary brain expansion. Our results indicate that genes regulating early brain and cranial growth incline to neoplasia later in life, irrespective of height. This warrants investigation of clinical implications of the link between head size and cancer

Etude expérimentale et numérique de l'écoulement autour d'un cylindre circulaire à un nombre de Reynolds égal à 3900

International audienceThis work contributes to the study of flow over a circular cylinder at Reynolds number Re= 3900. Although this classical flow is widely documented in the literature, especially for this precise Reynolds number that leads to a subcritical flow regime, there is no consensus about the turbulence statistics immediately just behind the obstacle. Here, the flow is investigated both numerically with large eddy simulation and experimentally with hot-wire anemometry and particle image velocimetry. The numerical simulation is performed using high-order schemes and a specific immersed boundary method. The present study focuses on turbulence statistics and power spectra in the near wake up to ten diameters. Statistical estimation is shown to need large integration times increasing the computational cost and leading to an uncertainty of about 10% for most flow characteristics considered in this study. The present numerical and experimental results are found to be in good agreement with previous large eddy simulation data. Contrary to this, the present results show differences compared to the experimental data found in the literature, the differences being larger than the estimated uncertainty range. Therefore, previous numerical-experimental controversy for this flow seems to be reduced with the data presented in this article

Etude numérique et expérimentale de l'écoulement autour d'un cylindre circulaire à ReD=3900

PIV measurements of the wake behind a circular cylinder at ReD =3900 are confronted with data of Large Eddy Simulation (LES) based on a modelling of the cylinder by a Direct Forcing method. New results, numerous and with a good degree of cofidence, are provided and offer a very good numerical/experiment agreement in the field of the turbulent statistics. The disparity with some results found in the literature is discussed.Des mesures PIV du sillage derrière un cylindre circulaire à ReD =3900 sont confrontées à des données de Simulation aux Grandes Echelles (SGE) basées sur une modélisation du cylindre par la méthode de Forçage Direct. Ces nouveaux résultats, complets et avec un bon niveau de confiance, offrent une très bonne concordance numérique/expérience sur le plan des statistiques turbulentes. La disparité avec certains résultats de la littérature est discutée