Automatic grid refinement criterion for lattice Boltzmann method

In all kinds of engineering problems, and in particular in methods for computational fluid dynamics based on regular grids, local grid refinement is of crucial importance. To save on computational expense, many applications require to resolve a wide range of scales present in a numerical simulation by locally adding more mesh points. In general, the need for a higher (or a lower) resolution is not known a priori, and it is therefore difficult to locate areas for which local grid refinement is required. In this paper, we propose a novel algorithm for the lattice Boltzmann method, based on physical concepts, to automatically construct a pattern of local refinement. We apply the idea to the two-dimensional lid-driven cavity and show that the automatically refined grid can lead to results of equal quality with less grid points, thus sparing computational resources and time. The proposed automatic grid refinement strategy has been implemented in the parallel open-source library Palabos

Human bony labyrinth is an indicator of population history and dispersal from Africa.

The dispersal of modern humans from Africa is now well documented with genetic data that track population history, as well as gene flow between populations. Phenetic skeletal data, such as cranial and pelvic morphologies, also exhibit a dispersal-from-Africa signal, which, however, tends to be blurred by the effects of local adaptation and in vivo phenotypic plasticity, and that is often deteriorated by postmortem damage to skeletal remains. These complexities raise the question of which skeletal structures most effectively track neutral population history. The cavity system of the inner ear (the so-called bony labyrinth) is a good candidate structure for such analyses. It is already fully formed by birth, which minimizes postnatal phenotypic plasticity, and it is generally well preserved in archaeological samples. Here we use morphometric data of the bony labyrinth to show that it is a surprisingly good marker of the global dispersal of modern humans from Africa. Labyrinthine morphology tracks genetic distances and geography in accordance with an isolation-by-distance model with dispersal from Africa. Our data further indicate that the neutral-like pattern of variation is compatible with stabilizing selection on labyrinth morphology. Given the increasingly important role of the petrous bone for ancient DNA recovery from archaeological specimens, we encourage researchers to acquire 3D morphological data of the inner ear structures before any invasive sampling. Such data will constitute an important archive of phenotypic variation in present and past populations, and will permit individual-based genotype-phenotype comparisons

Decline of genetic diversity in ancient domestic stallions in Europe.

Present-day domestic horses are immensely diverse in their maternally inherited mitochondrial DNA, yet they show very little variation on their paternally inherited Y chromosome. Although it has recently been shown that Y chromosomal diversity in domestic horses was higher at least until the Iron Age, when and why this diversity disappeared remain controversial questions. We genotyped 16 recently discovered Y chromosomal single-nucleotide polymorphisms in 96 ancient Eurasian stallions spanning the early domestication stages (Copper and Bronze Age) to the Middle Ages. Using this Y chromosomal time series, which covers nearly the entire history of horse domestication, we reveal how Y chromosomal diversity changed over time. Our results also show that the lack of multiple stallion lineages in the extant domestic population is caused by neither a founder effect nor random demographic effects but instead is the result of artificial selection-initially during the Iron Age by nomadic people from the Eurasian steppes and later during the Roman period. Moreover, the modern domestic haplotype probably derived from another, already advantageous, haplotype, most likely after the beginning of the domestication. In line with recent findings indicating that the Przewalski and domestic horse lineages remained connected by gene flow after they diverged about 45,000 years ago, we present evidence for Y chromosomal introgression of Przewalski horses into the gene pool of European domestic horses at least until medieval times

Predicting the Propagation of Acoustic Waves using Deep Convolutional Neural Networks

A novel approach for numerically propagating acoustic waves in two-dimensional quiescent media has been developed through a fully convolutional multi-scale neural network. This data-driven method managed to produce accurate results for long simulation times with a database of Lattice Boltzmann temporal simulations of propagating Gaussian Pulses, even in the case of initial conditions unseen during training time, such as the plane wave configuration or the two initial Gaussian pulses of opposed amplitudes. Two different choices of optimization objectives are compared, resulting in an improved prediction accuracy when adding the spatial gradient difference error to the traditional mean squared error loss function. Further accuracy gains are observed when performing an a posteriori correction on the neural network prediction based on the conservation of acoustic energy, indicating the benefit of including physical information in data-driven methods

Mitochondrial Genomes Reveal an Explosive Radiation of Extinct and Extant Bears near the Miocene-Pliocene Boundary

Background: Despite being one of the most studied families within the Carnivora, the phylogenetic relationships among the members of the bear family (Ursidae) have long remained unclear. Widely divergent topologies have been suggested based on various data sets and methods. Results: We present a fully resolved phylogeny for ursids based on ten complete mitochondrial genome sequences from all eight living and two recently extinct bear species, the European cave bear (Ursus spelaeus) and the American giant short-faced bear (Arctodus simus). The mitogenomic data yield a well-resolved topology for ursids, with the sloth bear at the basal position within the genus Ursus. The sun bear is the sister taxon to both the American and Asian black bears, and this clade is the sister clade of cave bear, brown bear and polar bear confirming a recent study on bear mitochondrial genomes. Conclusion: Sequences from extinct bears represent the third and fourth Pleistocene species for which complete mitochondrial genomes have been sequenced. Moreover, the cave bear specimen demonstrates that mitogenomic studies can be applied to Pleistocene fossils that have not been preserved in permafrost, and therefore have a broad application within ancient DNA research. Molecular dating of the mtDNA divergence times suggests a rapid radiation of bears in both the Old and New Worlds around 5 million years ago, at the Miocene-Pliocene boundary. This coincides with major global changes, such as the Messinian crisis and the first opening of the Bering Strait, and suggests a global influence of such events on species radiations

Leveraging the performance of LBM-HPC for large sizes on GPUs using ghost cells

Today, we are living a growing demand of larger and more efficient computational resources from the scientific community. On the other hand, the appearance of GPUs for general purpose computing supposed an important advance for covering such demand. These devices offer an impressive computational capacity at low cost and an efficient power consumption. However, the memory available in these devices is (sometimes) not enough, and so it is necessary computationally expensive memory transfers from (to) CPU to (from) GPU, causing a dramatic fall in performance. Recently, the Lattice-Boltzmann Method has positioned as an efficient methodology for fluid simulations. Although this method presents some interesting features particularly amenable to be efficiently exploited on parallel computers, it requires a considerable memory capacity, which can suppose an important drawback, in particular, on GPUs. In the present paper, it is proposed a new GPU-based implementation, which minimizes such requirements with respect to other state-of-the-art implementations. It allows us to execute almost 2$x$ bigger problems without additional memory transfers, achieving faster executions when dealing with large problems

Genomic insights into the origin and diversification of late maritime hunter-gatherers from the Chilean Patagonia.

Patagonia was the last region of the Americas reached by humans who entered the continent from Siberia ∼15,000-20,000 y ago. Despite recent genomic approaches to reconstruct the continental evolutionary history, regional characterization of ancient and modern genomes remains understudied. Exploring the genomic diversity within Patagonia is not just a valuable strategy to gain a better understanding of the history and diversification of human populations in the southernmost tip of the Americas, but it would also improve the representation of Native American diversity in global databases of human variation. Here, we present genome data from four modern populations from Central Southern Chile and Patagonia ( <i>n</i> = 61) and four ancient maritime individuals from Patagonia (∼1,000 y old). Both the modern and ancient individuals studied in this work have a greater genetic affinity with other modern Native Americans than to any non-American population, showing within South America a clear structure between major geographical regions. Native Patagonian Kawéskar and Yámana showed the highest genetic affinity with the ancient individuals, indicating genetic continuity in the region during the past 1,000 y before present, together with an important agreement between the ethnic affiliation and historical distribution of both groups. Lastly, the ancient maritime individuals were genetically equidistant to a ∼200-y-old terrestrial hunter-gatherer from Tierra del Fuego, which supports a model with an initial separation of a common ancestral group to both maritime populations from a terrestrial population, with a later diversification of the maritime groups

Wall model for large-eddy simulation based on the lattice Boltzmann method

International audienceAn original wall-model for lattice Boltzmann based large-eddy simulation is presented. It relies on either an analytical profile for the velocity profile within the first off-wall cell or the solution of Turbulent Boundary Layer Equations. The latter are solved thanks to a meshless approach to obtain a user-friendly method. The proposed method is observed to be both robust and accurate, yielding very satisfactory results in plane channel flow up to Reτ=20 000Reτ=20 000 on uniform grids with Δx+=Δy+=Δz+=1000Δx+=Δy+=Δz+=1000. The computed maximum error on skin friction is about 4% in most investigated cases on the coarsest grids (10 grid points across the half-channel height). It lies within 1–2% on finer grids with 20–60 uniformly distributed grid points in the channel half-height

Simulation of a generalized Newtonian fluid by the Lattice Boltzmann Method

articl

A lattice Boltzmann method for nonlinear disturbances around an arbitrary base flow

International audienceIn this paper we address the problem of the time evolution of a perturbation around a steady base flow with the use of the lattice Boltzmann method (LBM). This approach, named base flow lattice Boltzmann method, is of great interest in particular for aeroacoustic fields where the acoustic perturbation, on the one hand, is almost exclusively influenced by the large scale average structures of the underlying flow, and on the other hand, has a low effect on the large structures. The method is implemented for weakly compressible flows and the results of the base flow lattice Boltzmann are compared with the standard single relaxation time LBM. The boundary conditions for the base flow lattice Boltzmann method are discussed, as well as the implementation of outflow conditions for acoustic waves