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

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

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

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

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

Modular Groups, Visibility Diagram and Quantum Hall Effect

We consider the action of the modular group $\Gamma (2)$ on the set of positive rational fractions. From this, we derive a model for a classification of fractional (as well as integer) Hall states which can be visualized on two ``visibility" diagrams, the first one being associated with even denominator fractions whereas the second one is linked to odd denominator fractions. We use this model to predict, among some interesting physical quantities, the relative ratios of the width of the different transversal resistivity plateaus. A numerical simulation of the tranversal resistivity plot based on this last prediction fits well with the present experimental data.Comment: 17 pages, plain TeX, 4 eps figures included (macro epsf.tex), 1 figure available from reques

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

Analysis of model rotor blade pressures during parallel interaction with twin vortices

This paper presents and provides analysis of unsteady surface pressures measured on a model rotor blade as the blade experienced near parallel blade vortex interaction with a twin vortex system. To provide a basis for analysis, the vortex system was characterized by hot-wire measurements made in the interaction plane but in the absence of the rotor. The unsteady pressure response resulting from a single vortex interaction is then presented to provide a frame of reference for the twin vortex results. A series of twin vortex interaction cases are then presented and analyzed. It is shown that the unsteady blade pressures and forces are very sensitive to the inclination angle and separation distance of the vortex pair. When the vortex cores lie almost parallel to the blade chord, the interaction is characterized by a two-stage response associated with the sequential passage of the two cores. Conversely, when the cores lie on a plane that is almost perpendicular to the blade chord, the response is similar to that of a single vortex interaction. In all cases, the normal force response is consistent with the distribution of vertical velocity in the flow field of the vortex system. The pitching moment response, on the other hand, depends on the localized suction associated with the vortex cores as they traverse the blade chord