Multigroup radiation hydrodynamics with flux-limited diffusion and adaptive mesh refinement

International audienceContext. Radiative transfer plays a crucial role in the star formation process. Because of the high computational cost, radiation-hydrodynamics simulations performed up to now have mainly been carried out in the grey approximation. In recent years, multifrequency radiation-hydrodynamics models have started to be developed in an attempt to better account for the large variations in opacities as a function of frequency.Aims. We wish to develop an efficient multigroup algorithm for the adaptive mesh refinement code RAMSES which is suited to heavy proto-stellar collapse calculations.Methods. Because of the prohibitive timestep constraints of an explicit radiative transfer method, we constructed a time-implicit solver based on a stabilized bi-conjugate gradient algorithm, and implemented it in RAMSES under the flux-limited diffusion approximation.Results. We present a series of tests that demonstrate the high performance of our scheme in dealing with frequency-dependent radiation-hydrodynamic flows. We also present a preliminary simulation of a 3D proto-stellar collapse using 20 frequency groups. Differences between grey and multigroup results are briefly discussed, and the large amount of information this new method brings us is also illustrated.Conclusions. We have implemented a multigroup flux-limited diffusion algorithm in the RAMSES code. The method performed well against standard radiation-hydrodynamics tests, and was also shown to be ripe for exploitation in the computational star formation context

Modeling electricity loads in California: a continuous-time approach

In this paper we address the issue of modeling electricity loads and prices with diffusion processes. More specifically, we study models which belong to the class of generalized Ornstein-Uhlenbeck processes. After comparing properties of simulated paths with those of deseasonalized data from the California power market and performing out-of-sample forecasts we conclude that, despite certain advantages, the analyzed continuous-time processes are not adequate models of electricity load and price dynamics.Comment: To be published in Physica A (2001): Proceedings of the NATO ARW on Application of Physics in Economic Modelling, Prague, Feb. 8-10, 200

Scaling Behaviour and Complexity of the Portevin-Le Chatelier Effect

The plastic deformation of dilute alloys is often accompanied by plastic instabilities due to dynamic strain aging and dislocation interaction. The repeated breakaway of dislocations from and their recapture by solute atoms leads to stress serrations and localized strain in the strain controlled tensile tests, known as the Portevin-Le Chatelier (PLC) effect. In this present work, we analyse the stress time series data of the observed PLC effect in the constant strain rate tensile tests on Al-2.5%Mg alloy for a wide range of strain rates at room temperature. The scaling behaviour of the PLC effect was studied using two complementary scaling analysis methods: the finite variance scaling method and the diffusion entropy analysis. From these analyses we could establish that in the entire span of strain rates, PLC effect showed Levy walk property. Moreover, the multiscale entropy analysis is carried out on the stress time series data observed during the PLC effect to quantify the complexity of the distinct spatiotemporal dynamical regimes. It is shown that for the static type C band, the entropy is very low for all the scales compared to the hopping type B and the propagating type A bands. The results are interpreted considering the time and length scales relevant to the effect.Comment: 35 pages, 6 figure

Protostellar birth with ambipolar and ohmic diffusion

This is the final version of the article. Available from EDP Sciences via the DOI in this record.The transport of angular momentum is capital during the formation of low-mass stars; too little removal and rotation ensures stellar densities are never reached, too much and the absence of rotation means no protoplanetary disks can form. Magnetic diffusion is seen as a pathway to resolving this long-standing problem. We investigate the impact of including resistive MHD in simulations of the gravitational collapse of a 1 solar mass gas sphere, from molecular cloud densities to the formation of the protostellar seed; the second Larson core. We used the AMR code RAMSES to perform two 3D simulations of collapsing magnetised gas spheres, including self-gravity, radiative transfer, and a non-ideal gas equation of state to describe H2 dissociation which leads to the second collapse. The first run was carried out under the ideal MHD approximation, while ambipolar and ohmic diffusion was incorporated in the second calculation. In the ideal MHD simulation, the magnetic field dominates the energy budget everywhere inside and around the first core, fueling interchange instabilities and driving a low-velocity outflow. High magnetic braking removes essentially all angular momentum from the second core. On the other hand, ambipolar and ohmic diffusion create a barrier which prevents amplification of the magnetic field beyond 0.1 G in the first Larson core which is now fully thermally supported. A significant amount of rotation is preserved and a small Keplerian-like disk forms around the second core. When studying the radiative efficiency of the first and second core accretion shocks, we found that it can vary by several orders of magnitude over the 3D surface of the cores. Magnetic diffusion is a pre-requisite to star-formation; it enables the formation of protoplanetary disks in which planets will eventually form, and also plays a determinant role in the formation of the protostar itself.We are indebted to the anonymous referee for his/her insightful comments that have vastly improved the solidity of our study, with no stones left unturned. We also thank Troels Haugbølle for very useful discussions during the writing of this paper. NV gratefully acknowledges support from the European Commission through the Horizon 2020 Marie Skłodowska-Curie Actions Individual Fellowship 2014 programme (Grant Agreement no. 659706). The research leading to these results has also received funding from the European Research Council under the European Community’s Seventh Framework Programme (FP7/2007-2013 Grant Agreement no. 247060). We acknowledge financial support from "Programme National de Physique Stellaire" (PNPS) of CNRS/INSU, CEA and CNES, France. This work was granted access to the HPC resources of CINES (Occigen) under the allocation 2016-047247 made by GENCI. We also made use of the astrophysics HPC facility at the University of Copenhagen, which is supported by a research grant (VKR023406) from Villum Fonden. In addition, we thank the Service d’Astrophysique, IRFU, CEA Saclay, and the Laboratoire Astrophysique Instrumentation Modélisation, France, for granting us access to the supercomputer IRFUCOAST where the groundwork with many test calculations were performed. All the figures were created using the OSIRIS8 visualization package for RAMSES, except Fig. 4 which was rendered with the PARAVIEW9 software

La prise en compte de la durabilité dans les stratégies marketing des PME : apports et limites

This paper examines how SME managers in the agro-food sector include sustainable development into their marketing strategy. We focus on the link between sustainable development and performance, the impact of environmental factors on their practices and the limits that they perceive when trying to include sustainability into their marketing strategy. More than a marketing strength, sustainability appears as important for the people who work in the SME. However, there are some limits. ...French Abstract : Nous examinons dans cet article comment des directeurs de PME agroalimentaires à fort ancrage territorial s’approprient la notion de développement durable et l’expriment dans leurs stratégies marketing. Nous examinons successivement leur vision du lien entre développement durable et performance, l’impact des caractéristiques environnementales sur leurs pratiques et les limites qu’ils perçoivent à cette volonté d’intégrer la durabilité dans leur stratégie au moment de sa mise en œuvre. Au-delà de l’argument commercial, la durabilité apparaît comme un levier intéressant du projet d’entreprise et de la mobilisation des salariés. Néanmoins, la démarche comporte aussi des limites.SUSTAINABILITY; MARKETING STRATEGY; FOOD SECTOR; SMALL FIRMS

Simulations of protostellar collapse using multigroup radiation hydrodynamics. II. The second collapse

15 pages, 11 figures, accepted for publication in A&AStar formation begins with the gravitational collapse of a dense core inside a molecular cloud. As the collapse progresses, the centre of the core begins to heat up as it becomes optically thick. The temperature and density in the centre eventually reach high enough values where fusion reactions can ignite; the protostar is born. This sequence of events entail many physical processes, of which radiative transfer is of paramount importance. Many simulations of protostellar collapse make use of a grey treatment of radiative transfer coupled to the hydrodynamics. However, interstellar gas and dust opacities present large variations as a function of frequency. In this paper, we follow-up on a previous paper on the collapse and formation of Larson's first core using multigroup radiation hydrodynamics (Paper I) by extending the calculations to the second phase of the collapse and the formation of Larson's second core. We have made the use of a non-ideal gas equation of state as well as an extensive set of spectral opacities in a spherically symmetric fully implicit Godunov code to model all the phases of the collapse of a 0.1, 1 and 10 solar mass cloud cores. We find that, for a same central density, there are only small differences between the grey and multigroup simulations. The first core accretion shock remains supercritical while the shock at the second core border is found to be strongly subcritical with all the accreted energy being transfered to the core. The size of the first core was found to vary somewhat in the different simulations (more unstable clouds form smaller first cores) while the size, mass and temperature of the second cores are independent of initial cloud mass, size and temperature. Our simulations support the idea of a standard (universal) initial second core size of 0.003 AU and mass 0.0014 solar masses

RANS simulations of wind flow at the Bolund experiment

As part of their development, the predictions of numerical wind flow models must be compared with measurements in order to estimate the uncertainty related to their use. Of course, the most rigorous such comparison is under blind conditions. The following paper includes a detailed description of three different wind flow models, all based on a Reynolds-averaged Navier-Stokes approach and two-equation k-ε closure, that were tested as part of the Bolund blind comparison (itself based on the Bolund experiment which measured the wind around a small coastal island). The models are evaluated in terms of predicted normalized wind speed and turbulent kinetic energy at 2 m and 5 m above ground level for a westerly wind direction. Results show that all models predict the mean velocity reasonably well; however accurate prediction of the turbulent kinetic energy remains achallenge

Chemical solver to compute molecule and grain abundances and non-ideal MHD resistivities in prestellar core collapse calculations

This is the author accepted manuscript. The final version is available from the publisher via the DOI in this record.We develop a detailed chemical network relevant to the conditions characteristic of prestellar core collapse. We solve the system of time-dependent differential equations to calculate the equilibrium abundances of molecules and dust grains, with a size distribution given by size-bins for these latter. These abundances are used to compute the different non-ideal magneto-hydrodynamics resistivities (ambipolar, Ohmic and Hall), needed to carry out simulations of protostellar collapse. For the first time in this context, we take into account the evaporation of the grains, the thermal ionisation of Potassium, Sodium and Hydrogen at high temperature, and the thermionic emission of grains in the chemical network, and we explore the impact of various cosmic ray ionisation rates. All these processes significantly affect the non-ideal magneto-hydrodynamics resistivities, which will modify the dynamics of the collapse. Ambipolar diffusion and Hall effect dominate at low densities, up to n_H = 10^12 cm^-3, after which Ohmic diffusion takes over. We find that the time-scale needed to reach chemical equilibrium is always shorter than the typical dynamical (free fall) one. This allows us to build a large, multi-dimensional multi-species equilibrium abundance table over a large temperature, density and ionisation rate ranges. This table, which we make accessible to the community, is used during first and second prestellar core collapse calculations to compute the non-ideal magneto-hydrodynamics resistivities, yielding a consistent dynamical-chemical description of this process