Towards higher order lattice Boltzmann schemes

In this contribution we extend the Taylor expansion method proposed previously by one of us and establish equivalent partial differential equations of DDH lattice Boltzmann scheme at an arbitrary order of accuracy. We derive formally the associated dynamical equations for classical thermal and linear fluid models in one to three space dimensions. We use this approach to adjust relaxation parameters in order to enforce fourth order accuracy for thermal model and diffusive relaxation modes of the Stokes problem. We apply the resulting scheme for numerical computation of associated eigenmodes and compare our results with analytical references

Optimal low-dispersion low-dissipation LBM schemes for computational aeroacoustics

Lattice Boltmzmann Methods (LBM) have been proved to be very effective methods for computational aeroacoustics (CAA), which have been used to capture the dynamics of weak acoustic fluctuations. In this paper, we propose a strategy to reduce the dispersive and disspative errors of the two-dimensional (2D) multi-relaxation-time lattice Boltzmann method (MRT-LBM). By presenting an effective algorithm, we obtain a uniform form of the linearized Navier-Stokes equations corresponding to the MRT-LBM in wave-number space. Using the matrix perturbation theory and the equivalent modified equation approach for finite difference methods, we propose a class of minimization problems to optimize the free-parameters in the MRT-LBM. We obtain this way a dispersion-relation-preserving LBM (DRP-LBM) to circumvent the minimized dispersion error of the MRT-LBM. The dissipation relation precision is also improved.And the stability of the MRT-LBM with the small bulk viscosity is guaranteed. Von Neuman analysis of the linearized MRT-LBM is performed to validate the optimized dispersion/dissipation relations considering monochromatic wave solutions. Meanwhile, dispersion and dissipation errors of the optimized MRT-LBM are quantitatively compared with the original MRT-LBM . Finally, some numerical simulations are carried out to assess the new optimized MRT-LBM schemes.Comment: 33 page

Reconstructing past terrace fields in the Pyrenees: Insights into land management and settlement from the Bronze Age to the Early Modern era at Vilalta (1650 masl, Cerdagne, France)

© Trustees of Boston University 2015. The building of a solar power station at Thémis, at 1650 masl on the south-facing slope of the Carlit massif in the eastern Pyrenees, led to an archaeological evaluation from April-June 2009. This evaluation covered a surface of 10 ha that included a medieval village as well as the surrounding agricultural land in terraces. Non-destructive archaeological methods were used for the village. A detailed study of the 6 ha of terraces began with a fieldwalking survey, mapping every visible feature, followed by systematic trial trenches. Fifty-five trenches, 11 in the village and 44 in the fields, were opened. The stratigraphies were then compared with a series of 22 radiocarbon dates and eight relative dates provided by ceramic typologies. This combination of surface and buried evidence supported our preliminary hypothesis about the dynamics of the slope. The results suggest the existence of agrarian features beginning in the Bronze Age and reveal that the field patterns were frequently transformed, both in the Medieval and Early Modern periods. The transformations in the terrace fields after the village was abandoned are as interesting as those during occupation because, contrary to the idea of a fixed, unchanging landscape after the end of the Middle Ages, they challenge the idea that mountain zones are marginal spaces by nature, or were marginalized later.Peer Reviewe

On the Three-dimensional Central Moment Lattice Boltzmann Method

A three-dimensional (3D) lattice Boltzmann method based on central moments is derived. Two main elements are the local attractors in the collision term and the source terms representing the effect of external and/or self-consistent internal forces. For suitable choices of the orthogonal moment basis for the three-dimensional, twenty seven velocity (D3Q27), and, its subset, fifteen velocity (D3Q15) lattice models, attractors are expressed in terms of factorization of lower order moments as suggested in an earlier work; the corresponding source terms are specified to correctly influence lower order hydrodynamic fields, while avoiding aliasing effects for higher order moments. These are achieved by successively matching the corresponding continuous and discrete central moments at various orders, with the final expressions written in terms of raw moments via a transformation based on the binomial theorem. Furthermore, to alleviate the discrete effects with the source terms, they are treated to be temporally semi-implicit and second-order, with the implicitness subsequently removed by means of a transformation. As a result, the approach is frame-invariant by construction and its emergent dynamics describing fully 3D fluid motion in the presence of force fields is Galilean invariant. Numerical experiments for a set of benchmark problems demonstrate its accuracy.Comment: 55 pages, 8 figure

Cytoplasmic PML promotes TGF-β-associated epithelial–mesenchymal transition and invasion in prostate cancer

Epithelial–mesenchymal transition (EMT) is a key event that is involved in the invasion and dissemination of cancer cells. Although typically considered as having tumour-suppressive properties, transforming growth factor (TGF)-β signalling is altered during cancer and has been associated with the invasion of cancer cells and metastasis. In this study, we report a previously unknown role for the cytoplasmic promyelocytic leukaemia (cPML) tumour suppressor in TGF-β signalling-induced regulation of prostate cancer-associated EMT and invasion. We demonstrate that cPML promotes a mesenchymal phenotype and increases the invasiveness of prostate cancer cells. This event is associated with activation of TGF-β canonical signalling pathway through the induction of Sma and Mad related family 2 and 3 (SMAD2 and SMAD3) phosphorylation. Furthermore, the cytoplasmic localization of promyelocytic leukaemia (PML) is mediated by its nuclear export in a chromosomal maintenance 1 (CRM1)-dependent manner. This was clinically tested in prostate cancer tissue and shown that cytoplasmic PML and CRM1 co-expression correlates with reduced disease-specific survival. In summary, we provide evidence of dysfunctional TGF-β signalling occurring at an early stage in prostate cancer. We show that this disease pathway is mediated by cPML and CRM1 and results in a more aggressive cancer cell phenotype. We propose that the targeting of this pathway could be therapeutically exploited for clinical benefit

Dynamics of fluctuations in a fluid below the onset of Rayleigh-B\'enard convection

We present experimental data and their theoretical interpretation for the decay rates of temperature fluctuations in a thin layer of a fluid heated from below and confined between parallel horizontal plates. The measurements were made with the mean temperature of the layer corresponding to the critical isochore of sulfur hexafluoride above but near the critical point where fluctuations are exceptionally strong. They cover a wide range of temperature gradients below the onset of Rayleigh-B\'enard convection, and span wave numbers on both sides of the critical value for this onset. The decay rates were determined from experimental shadowgraph images of the fluctuations at several camera exposure times. We present a theoretical expression for an exposure-time-dependent structure factor which is needed for the data analysis. As the onset of convection is approached, the data reveal the critical slowing-down associated with the bifurcation. Theoretical predictions for the decay rates as a function of the wave number and temperature gradient are presented and compared with the experimental data. Quantitative agreement is obtained if allowance is made for some uncertainty in the small spacing between the plates, and when an empirical estimate is employed for the influence of symmetric deviations from the Oberbeck-Boussinesq approximation which are to be expected in a fluid with its density at the mean temperature located on the critical isochore.Comment: 13 pages, 10 figures, 52 reference

Prehistory of Transit Searches

Nowadays the more powerful method to detect extrasolar planets is the transit method. We review the planet transits which were anticipated, searched, and the first ones which were observed all through history. Indeed transits of planets in front of their star were first investigated and studied in the solar system. The first observations of sunspots were sometimes mistaken for transits of unknown planets. The first scientific observation and study of a transit in the solar system was the observation of Mercury transit by Pierre Gassendi in 1631. Because observations of Venus transits could give a way to determine the distance Sun-Earth, transits of Venus were overwhelmingly observed. Some objects which actually do not exist were searched by their hypothetical transits on the Sun, as some examples a Venus satellite and an infra-mercurial planet. We evoke the possibly first use of the hypothesis of an exoplanet transit to explain some periodic variations of the luminosity of a star, namely the star Algol, during the eighteen century. Then we review the predictions of detection of exoplanets by their transits, those predictions being sometimes ancient, and made by astronomers as well as popular science writers. However, these very interesting predictions were never published in peer-reviewed journals specialized in astronomical discoveries and results. A possible transit of the planet beta Pic b was observed in 1981. Shall we see another transit expected for the same planet during 2018? Today, some studies of transits which are connected to hypothetical extraterrestrial civilisations are published in astronomical refereed journals. Some studies which would be classified not long ago as science fiction are now considered as scientific ones.Comment: Submiited to Handbook of Exoplanets (Springer

Back-arc strain in subduction zones: Statistical observations versus numerical modeling

International audience1] Recent statistical analysis by Lallemand et al. (2008) of subduction zone parameters revealed that the back-arc deformation mode depends on the combination between the subducting (nu(sub)) and upper (nu(up)) plate velocities. No significant strain is recorded in the arc area if plate kinematics verifies nu(up) = 0.5 vsub - 2.3 (cm/a) in the HS3 reference frame. Arc spreading ( shortening) occurs if nu(up) is greater ( lower) than the preceding relationship. We test this statistical law with numerical models of subduction, by applying constant plate velocities far away from the subduction zone. The subducting lithosphere is free to deform at all depths. We quantify the force applied on the two converging plates to sustain constant surface velocities. The simulated rheology combined viscous (non-Newtonian) and brittle behaviors, and depends on water content. The influence of subduction rate vs is first studied for a fixed upper plate. After 950 km of convergence ( steady state slab pull), the transition from extensional to compressive stresses in the upper plate occurs for vs similar to 1.4 cm/a. The effect of upper plate velocity is then tested at constant subduction rate. Upper plate retreat ( advance) with respect to the trench increases extension ( compression) in the arc lithosphere and increases ( decreases) the subducting plate dip. Our modeling confirms the statistical kinematic relationship between vsub and nu(up) that describes the transition from extensional to compressive stresses in the arc lithosphere, even if the modeled law is shifted toward higher rates of upper plate retreat, using our set of physical parameters ( e. g., 100 km thick subducting oceanic plate) and short- term simulations. Our results make valid the choice of the HS3 reference frame for assessing plate velocity influence on arc tectonic regime. The subduction model suggests that friction along the interplate contact and the mantle Stokes reaction could be the two main forces competing against slab pull for upper mantle subductions. Besides, our simulations show that the arc deformation mode is strongly time dependent

Lattice Boltzmann simulations of soft matter systems

This article concerns numerical simulations of the dynamics of particles immersed in a continuum solvent. As prototypical systems, we consider colloidal dispersions of spherical particles and solutions of uncharged polymers. After a brief explanation of the concept of hydrodynamic interactions, we give a general overview over the various simulation methods that have been developed to cope with the resulting computational problems. We then focus on the approach we have developed, which couples a system of particles to a lattice Boltzmann model representing the solvent degrees of freedom. The standard D3Q19 lattice Boltzmann model is derived and explained in depth, followed by a detailed discussion of complementary methods for the coupling of solvent and solute. Colloidal dispersions are best described in terms of extended particles with appropriate boundary conditions at the surfaces, while particles with internal degrees of freedom are easier to simulate as an arrangement of mass points with frictional coupling to the solvent. In both cases, particular care has been taken to simulate thermal fluctuations in a consistent way. The usefulness of this methodology is illustrated by studies from our own research, where the dynamics of colloidal and polymeric systems has been investigated in both equilibrium and nonequilibrium situations.Comment: Review article, submitted to Advances in Polymer Science. 16 figures, 76 page

Ebf factors and MyoD cooperate to regulate muscle relaxation via Atp2a1

Jin, Saihong et al.Myogenic regulatory factors such as MyoD and Myf5 lie at the core of vertebrate muscle differentiation. However, E-boxes, the cognate binding sites for these transcription factors, are not restricted to the promoters/enhancers of muscle cell-specific genes. Thus, the specificity in myogenic transcription is poorly defined. Here we describe the transcription factor Ebf3 as a new determinant of muscle cell-specific transcription. In the absence of Ebf3 the lung does not unfold at birth, resulting in respiratory failure and perinatal death. This is due to a hypercontractile diaphragm with impaired Ca2+ efflux-related muscle functions. Expression of the Ca2+ pump Serca1 (Atp2a1) is downregulated in the absence of Ebf3, and its transgenic expression rescues this phenotype. Ebf3 binds directly to the promoter of Atp2a1 and synergises with MyoD in the induction of Atp2a1. In skeletal muscle, the homologous family member Ebf1 is strongly expressed and together with MyoD induces Atp2a1. Thus, Ebf3 is a new regulator of terminal muscle differentiation in the diaphragm, and Ebf factors cooperate with MyoD in the induction of muscle-specific genes. © 2014 Macmillan Publishers Limited.This work was supported by grants from the German Research Foundation (DFG, TRR54; FOR1586; FOR2033) and by a stipend of the Max Planck SocietyPeer Reviewe