-
New insight on Young Stellar Objects accretion shocks -- a claim for NLTE opacities
Authors:
Lionel de Sá,
Jean-Pierre Chièze,
Chantal Stehlé,
Ivan Hubeny,
Thierry Lanz,
Véronique Cayatte
Abstract:
Context. Accreted material onto CTTSs is expected to form a hot quasi-periodic plasma structure that radiates in X-rays. Simulations of this phenomenon only partly match with observations. They all rely on a static model for the chromosphere model and on the assumption that radiation and matter are decoupled. Aims. We explore the effects on the structure and on the dynamics of the accretion flow o…
▽ More
Context. Accreted material onto CTTSs is expected to form a hot quasi-periodic plasma structure that radiates in X-rays. Simulations of this phenomenon only partly match with observations. They all rely on a static model for the chromosphere model and on the assumption that radiation and matter are decoupled. Aims. We explore the effects on the structure and on the dynamics of the accretion flow of both a shock-heated chromosphere and of the coupling between radiation and hydrodynamics. Methods. We simulate accretion columns falling onto a stellar chromosphere using the 1D ALE code AstroLabE. This code solves the hydrodynamics equations along with the two first momenta equations for radiation transfer, with the help of a dedicated opacity table for the coupling between matter and radiation. We derive the total electron and ions densities from collisional-radiative NLTE ionisation equilibrium. Results. The chromospheric acoustic heating has an impact on the duration of the cycle and on the structure of the heated slab. In addition, the coupling between radiation and hydrodynamics leads to a heating of the accretion flow and the chromosphere, inducing a possible unburial of the whole column. These two last conclusions are in agreement with the computed monochromatic intensity. Both effects (acoustic heating and radiation coupling) have an influence on the amplitude and temporal variations of the net X-ray luminosity, which varies between 30 and 94% of the incoming mechanical energy flux, depending on the model considered.
△ Less
Submitted 10 August, 2019; v1 submitted 19 April, 2019;
originally announced April 2019.
-
Young stellar object jet models: From theory to synthetic observations
Authors:
O. Tesileanu,
T. Matsakos,
S. Massaglia,
E. Trussoni,
A. Mignone,
N. Vlahakis,
K. Tsinganos,
M. Stute,
V. Cayatte,
C. Sauty,
C. Stehle,
J. -P. Chieze
Abstract:
Astronomical observations, analytical solutions and numerical simulations have provided the building blocks to formulate the current theory of young stellar object jets. Although each approach has made great progress independently, it is only during the last decade that significant efforts are being made to bring the separate pieces together. Building on previous work that combined analytical solu…
▽ More
Astronomical observations, analytical solutions and numerical simulations have provided the building blocks to formulate the current theory of young stellar object jets. Although each approach has made great progress independently, it is only during the last decade that significant efforts are being made to bring the separate pieces together. Building on previous work that combined analytical solutions and numerical simulations, we apply a sophisticated cooling function to incorporate optically thin energy losses in the dynamics. On the one hand, this allows a self-consistent treatment of the jet evolution and on the other, it provides the necessary data to generate synthetic emission maps. Firstly, analytical disk and stellar outflow solutions are properly combined to initialize numerical two-component jet models inside the computational box. Secondly, magneto-hydrodynamical simulations are performed in 2.5D, following properly the ionization and recombination of a maximum of $29$ ions. Finally, the outputs are post-processed to produce artificial observational data. The first two-component jet simulations, based on analytical models, that include ionization and optically thin radiation losses demonstrate promising results for modeling specific young stellar object outflows. The generation of synthetic emission maps provides the link to observations, as well as the necessary feedback for the further improvement of the available models.
△ Less
Submitted 16 December, 2013;
originally announced December 2013.
-
Radiative accretion shocks along nonuniform stellar magnetic fields in classical T Tauri stars
Authors:
S. Orlando,
R. Bonito,
C. Argiroffi,
F. Reale,
G. Peres,
M. Miceli,
T. Matsakos,
C. Stehlé,
L. Ibgui,
L. de Sá,
J. P. Chièze,
T. Lanz
Abstract:
(abridged) AIMS. We investigate the dynamics and stability of post-shock plasma streaming along nonuniform stellar magnetic fields at the impact region of accretion columns. We study how the magnetic field configuration and strength determine the structure, geometry, and location of the shock-heated plasma. METHODS. We model the impact of an accretion stream onto the chromosphere of a CTTS by 2D a…
▽ More
(abridged) AIMS. We investigate the dynamics and stability of post-shock plasma streaming along nonuniform stellar magnetic fields at the impact region of accretion columns. We study how the magnetic field configuration and strength determine the structure, geometry, and location of the shock-heated plasma. METHODS. We model the impact of an accretion stream onto the chromosphere of a CTTS by 2D axisymmetric magnetohydrodynamic simulations. Our model takes into account the gravity, the radiative cooling, and the magnetic-field-oriented thermal conduction. RESULTS. The structure, stability, and location of the shocked plasma strongly depend on the configuration and strength of the magnetic field. For weak magnetic fields, a large component of B may develop perpendicular to the stream at the base of the accretion column, limiting the sinking of the shocked plasma into the chromosphere. An envelope of dense and cold chromospheric material may also develop around the shocked column. For strong magnetic fields, the field configuration determines the position of the shock and its stand-off height. If the field is strongly tapered close to the chromosphere, an oblique shock may form well above the stellar surface. In general, a nonuniform magnetic field makes the distribution of emission measure vs. temperature of the shocked plasma lower than in the case of uniform magnetic field. CONCLUSIONS. The initial strength and configuration of the magnetic field in the impact region of the stream are expected to influence the chromospheric absorption and, therefore, the observability of the shock-heated plasma in the X-ray band. The field strength and configuration influence also the energy balance of the shocked plasma, its emission measure at T > 1 MK being lower than expected for a uniform field. The above effects contribute in underestimating the mass accretion rates derived in the X-ray band.
△ Less
Submitted 19 September, 2013;
originally announced September 2013.
-
YSO accretion shocks: magnetic, chromospheric or stochastic flow effects can suppress fluctuations of X-ray emission
Authors:
T. Matsakos,
J. -P. Chièze,
C. Stehlé,
M. González,
L. Ibgui,
L. de Sá,
T. Lanz,
S. Orlando,
R. Bonito,
C. Argiroffi,
F. Reale,
G. Peres
Abstract:
Context. Theoretical arguments and numerical simulations of radiative shocks produced by the impact of the accreting gas onto young stars predict quasi-periodic oscillations in the emitted radiation. However, observational data do not show evidence of such periodicity. Aims. We investigate whether physically plausible perturbations in the accretion column or in the chromosphere could disrupt the s…
▽ More
Context. Theoretical arguments and numerical simulations of radiative shocks produced by the impact of the accreting gas onto young stars predict quasi-periodic oscillations in the emitted radiation. However, observational data do not show evidence of such periodicity. Aims. We investigate whether physically plausible perturbations in the accretion column or in the chromosphere could disrupt the shock structure influencing the observability of the oscillatory behavior. Methods. We performed local 2D magneto-hydrodynamical simulations of an accretion shock impacting a chromosphere, taking optically thin radiation losses and thermal conduction into account. We investigated the effects of several perturbation types, such as clumps in the accretion stream or chromospheric fluctuations, and also explored a wide range of plasma-βvalues. Results. In the case of a weak magnetic field, the post-shock region shows chaotic motion and mixing, smoothing out the perturbations and retaining a global periodic signature. On the other hand, a strong magnetic field confines the plasma in flux tubes, which leads to the formation of fibrils that oscillate independently. Realistic values for the amplitude, length, and time scales of the perturbation are capable of bringing the fibril oscillations out of phase, suppressing the periodicity of the emission. Conclusions. The strength of a locally uniform magnetic field in YSO accretion shocks determines the structure of the post-shock region, namely, whether it will be somewhat homogeneous or if it will split up to form a collection of fibrils. In the second case, the size and shape of the fibrils is found to depend strongly on the plasma-βvalue but not on the perturbation type. Therefore, the actual value of the protostellar magnetic field is expected to play a critical role in the time dependence of the observable emission.
△ Less
Submitted 20 July, 2013;
originally announced July 2013.
-
Deflagration to detonation transition by amplification of acoustic waves in type Ia supernovae
Authors:
Camille Charignon,
Jean-Pierre Chièze
Abstract:
We study a new mechanism for deflagration to detonation transition in thermonuclear supernovae (SNe Ia), based on the formation of shocks by amplification of sound waves in the steep density gradients of white dwarfs envelopes. Given a large enough jump in density a small pressure and velocity perturbation, produced by the turbulent deflagration, turns into a shock down of the gradient, where it w…
▽ More
We study a new mechanism for deflagration to detonation transition in thermonuclear supernovae (SNe Ia), based on the formation of shocks by amplification of sound waves in the steep density gradients of white dwarfs envelopes. Given a large enough jump in density a small pressure and velocity perturbation, produced by the turbulent deflagration, turns into a shock down of the gradient, where it will dissipate and heat up the media. With the right frequency and amplitude the heating can be enough to initiate a detonation, which can propagate backward and up the density gradient. We studied planar and spherical geometry. In the planar case we made a parametric study of the frequency and amplitude. We found it possible to obtain a detonation for perturbations down to Mach number M=0.003. In the spherical case, geometrical damping makes it harder to initiate a detonation, but considering a small He atmosphere (<0.01 Msol) makes it possible again to obtain a detonation down to small perturbation (M=0.002). In the context of thermonuclear supernovae, this could be a mean to turn a turbulent flame producing sound waves to a detonation.
△ Less
Submitted 8 January, 2013;
originally announced January 2013.
-
Physical and radiative properties of the first core accretion shock
Authors:
Benoît Commerçon,
Edouard Audit,
Gilles Chabrier,
Jean-Pierre Chièze
Abstract:
Radiative shocks play a dominant role in star formation. The accretion shocks on the first and second Larson's cores involve radiative processes and are thus characteristic of radiative shocks. In this study, we explore the formation of the first Larson's core and characterize the radiative and dynamical properties of the accretion shock, using both analytical and numerical approaches. We develop…
▽ More
Radiative shocks play a dominant role in star formation. The accretion shocks on the first and second Larson's cores involve radiative processes and are thus characteristic of radiative shocks. In this study, we explore the formation of the first Larson's core and characterize the radiative and dynamical properties of the accretion shock, using both analytical and numerical approaches. We develop both numerical RHD calculations and a semi-analytical model that characterize radiative shocks in various physical conditions, for radiating or barotropic fluids. Then, we perform 1D spherical collapse calculations of the first Larson's core, using a grey approximation for the opacity of the material. We consider three different models for radiative transfer, namely: the barotropic approximation, the FLD approximation and the more complete M1 model. We investigate the characteristic properties of the collapse and of the first core formation. Comparison between the numerical results and our semi-analytical model shows that this latter reproduces quite well the core properties obtained with the numerical calculations. The accretion shock on the first Larson core is found to be supercritical, implying that all the accretion shock energy on the core is radiated away. The FLD approximation is found to agree quite well with the results based on the M1 model, and is thus appropriate to study the star formation process. In contrast, the barotropic approximation does not correctly describe the thermal properties of the gas during the collapse. We show that a consistent treatment of radiation and hydrodynamics is mandatory to correctly handle the cooling of the gas during the core formation and thus to obtain the correct mechanical and thermal properties for this latter.
△ Less
Submitted 14 February, 2011;
originally announced February 2011.
-
The dynamical influence of cooling in the envelope of prestellar and protostellar cores
Authors:
P. Lesaffre,
A. Belloche,
J. -P. Chieze,
P. Andre
Abstract:
We compute numerical simulations of spherical collapse triggered by a slow increase in external pressure. We compare isothermal models to models including cooling with a simple but self-consistent treatment of the coupling between gas, grains and radiation field temperatures. The hydrostatic equilibrium appears to hold past the marginally stable state, until the collapse proceeds. The last hydro…
▽ More
We compute numerical simulations of spherical collapse triggered by a slow increase in external pressure. We compare isothermal models to models including cooling with a simple but self-consistent treatment of the coupling between gas, grains and radiation field temperatures. The hydrostatic equilibrium appears to hold past the marginally stable state, until the collapse proceeds. The last hydrostatic state before collapse has a lower central gas temperature in the centre due to the enhanced coupling between gas, grains and radiation field. This results in slightly lower pressure gradients in the bulk of the envelope which is hence slightly more extended than in the isothermal case. Due to the sensitivity of the collapse on these initial conditions, protostellar infall velocities in the envelope turn out to be much slower in the case with cooling.Our models also compute the radiative transfer and a rather large chemical network coupled to gas dynamics. However, we note that the steady-state chemisorption of CO is sufficient to provide an accurate cooling function of the gas. This justifies the use of post-processing techniques to account for the abundance of observed molecules. Existing observations of infall signatures put very stringent constraints on the kinematics and temperature profile of the class 0 protostar IRAM 04191+1522. We show that isothermal models fail to account for the innermost slow infall motions observed, even with the most hydrostatic initial conditions. In contrast, models with cooling reproduce the general shape of the temperature profile inferred from observations and are in much better agreement with the infall signatures in the inner 3000 AU.
△ Less
Submitted 5 August, 2005;
originally announced August 2005.
-
Temporal evolution of magnetic molecular shocks II. Analytics of the steady state and semi-analytical construction of intermediate ages
Authors:
P. Lesaffre,
J. -P. Chièze,
S. Cabrit,
G. Pineau des Forêts
Abstract:
In the first paper of this series (Paper I) we computed time dependent simulations of multifluid shocks with chemistry and a transverse magnetic field frozen in the ions, using an adaptive moving grid. In this paper, we present new analytical results on steady-state molecular shocks. Relationships between density and pressure in the neutral fluid are derived for the cold magnetic precursor, hot…
▽ More
In the first paper of this series (Paper I) we computed time dependent simulations of multifluid shocks with chemistry and a transverse magnetic field frozen in the ions, using an adaptive moving grid. In this paper, we present new analytical results on steady-state molecular shocks. Relationships between density and pressure in the neutral fluid are derived for the cold magnetic precursor, hot magnetic precursor, adiabatic shock front, and the following cooling layer. The compression ratio and temperature behind a fully dissociative adiabatic shock is also derived. To prove that these results may even hold for intermediate ages, we design a test to locally characterise the validity of the steady state equations in a time-dependent shock simulation. Applying this tool to the results of Paper I, we show that most of these shocks (all the stable ones) are indeed in a quasi-steady state at all times, i.e. : a given snapshot is composed of one or more truncated steady shock. Finally, we use this property to produce a construction method of any intermediate time of low velocity shocks (u < 20 km/s) with only a steady-state code. In particular, this method allows one to predict the occurrence of steady CJ-type shocks more accurately than previously proposed criteria.
△ Less
Submitted 28 July, 2004;
originally announced July 2004.
-
Temporal evolution of magnetic molecular shocks I. Moving grid simulations
Authors:
P. Lesaffre,
J. -P. Chièze,
S. Cabrit,
G. Pineau des Forêts
Abstract:
We present time-dependent 1D simulations of multifluid magnetic shocks with chemistry resolved down to the mean free path. They are obtained with an adaptive moving grid implemented with an implicit scheme. We examine a broad range of parameters relevant to conditions in dense molecular clouds, with preshock densities between 10^3 and 10^5 cm-3, velocities between 10 and 40 km/s, and three diffe…
▽ More
We present time-dependent 1D simulations of multifluid magnetic shocks with chemistry resolved down to the mean free path. They are obtained with an adaptive moving grid implemented with an implicit scheme. We examine a broad range of parameters relevant to conditions in dense molecular clouds, with preshock densities between 10^3 and 10^5 cm-3, velocities between 10 and 40 km/s, and three different scalings for the transverse magnetic field: B=0,0.1,1 μG \sqrt{n.cm3}. We first use this study to validate the results of Chièze, Pineau des Forêts & Flower (1998), in particular the long delays necessary to obtain steady C-type shocks, and we provide evolutionary time-scales for a much greater range of parameters. We also present the first time-dependent models of dissociative shocks with a magnetic precursor, including the first models of stationary CJ shocks in molecular conditions. We find that the maximum speed for steady C-type shocks is reached before the occurrence of a sonic point in the neutral fluid, unlike previously thought. As a result, the maximum speed for C-shocks is lower than previously believed. Finally, we find a large amplitude bouncing instability in J-type fronts near the H2 dissociation limit (u ~ 25-30 km/s), driven by H2 dissociation/reformation. At higher speeds, we find an oscillatory behaviour of short period and small amplitude linked to collisional ionisation of H. Both instabilities are suppressed after some time when a magnetic field is present. In a companion paper, we use the present simulations to validate a new semi-analytical construction method for young low-velocity magnetic shocks based on truncated steady-state models.
△ Less
Submitted 28 July, 2004;
originally announced July 2004.
-
A radiation-hydrodynamics scheme valid from the transport to the diffusion limit
Authors:
E. Audit,
P. Charrier,
J. -P. Chièze,
B. Dubroca
Abstract:
We present in this paper the numerical treatment of the coupling between hydrodynamics and radiative transfer. The fluid is modeled by classical conservation laws (mass, momentum and energy) and the radiation by the grey moment $M_1$ system. The scheme introduced is able to compute accurate numerical solution over a broad class of regimes from the transport to the diffusive limits. We propose an…
▽ More
We present in this paper the numerical treatment of the coupling between hydrodynamics and radiative transfer. The fluid is modeled by classical conservation laws (mass, momentum and energy) and the radiation by the grey moment $M_1$ system. The scheme introduced is able to compute accurate numerical solution over a broad class of regimes from the transport to the diffusive limits. We propose an asymptotic preserving modification of the HLLE scheme in order to treat correctly the diffusion limit. Several numerical results are presented, which show that this approach is robust and have the correct behavior in both the diffusive and free-streaming limits. In the last numerical example we test this approach on a complex physical case by considering the collapse of a gas cloud leading to a proto-stellar structure which, among other features, exhibits very steep opacity gradients.
△ Less
Submitted 17 June, 2002;
originally announced June 2002.
-
Non equilibrium thermodynamics and cosmological pancakes formation
Authors:
Romain Teyssier,
Jean-Pierre Chieze,
Jean-Michel Alimi
Abstract:
We investigate the influence of non equilibrium thermodynamics on cosmological structure formation. In this paper, we consider the collapse of planar perturbations usually called "Zel'dovich pancakes". We have developed for that purpose a new two fluids (gas and dark matter) hydrodynamical code, with three different thermodynamical species: electrons, ions and neutral particles (T_e\ne T_i \ne T…
▽ More
We investigate the influence of non equilibrium thermodynamics on cosmological structure formation. In this paper, we consider the collapse of planar perturbations usually called "Zel'dovich pancakes". We have developed for that purpose a new two fluids (gas and dark matter) hydrodynamical code, with three different thermodynamical species: electrons, ions and neutral particles (T_e\ne T_i \ne T_n). We describe in details the complex structure of accretion shock waves. We include several relevant processes for a low density, high temperature, collisional plasma such as non-equilibrium chemical reactions, cooling, shock heating, thermal energy equipartition between electrons, ions and neutral particles and electronic conduction. We find two different regions in the pancake structure: a thermal precursor ahead of the compression front and an equipartition wave after the compression front where electrons and ions temperatures differ significantly. This complex structure may have two interesting consequences: pre-heating of unshocked regions in the vicinity of massive X-ray clusters and ions and electrons temperatures differences in the outer regions of X-rays clusters.
△ Less
Submitted 12 August, 1998;
originally announced August 1998.
-
Multi-phase Hydrodynamics and X-ray Clusters Formation
Authors:
Romain Teyssier,
Jean-Pierre CHIEZE,
Jean-Michel Alimi
Abstract:
We investigate the role of radiative cooling within the core of large X-ray clusters using multi-phase hydrodynamics. We developed for that purpose a spherically symmetric hydrodynamical code, coupled to a "fluid model" that describes accurately the dark matter component. Cooling is included using a self-consistent multi-phase approach, leading to "cooled gas" mass deposition throughout the flow…
▽ More
We investigate the role of radiative cooling within the core of large X-ray clusters using multi-phase hydrodynamics. We developed for that purpose a spherically symmetric hydrodynamical code, coupled to a "fluid model" that describes accurately the dark matter component. Cooling is included using a self-consistent multi-phase approach, leading to "cooled gas" mass deposition throughout the flow. We simulate the collapse and the subsequent evolution of a Coma-like X-ray cluster, avoiding the well-known "cooling catastrophe". The total mass profile of our simulated cluster is very similar to the "universal" profile proposed by Navarro, Frenk & White (1995). More interestingly, we also obtain a quasi-isothermal temperature profile, which is a direct consequence of multi-phase cooling within such a potential well.
△ Less
Submitted 1 October, 1997;
originally announced October 1997.
-
Are Large X-ray Clusters at Thermal Equilibrium ?
Authors:
Jean-Pierre Chieze,
Romain Teyssier,
Jean-Michel Alimi
Abstract:
We simulate the formation of a large X-ray cluster using a fully 3D hydrodynamical code coupled to a Particle-Mesh scheme which models the dark matter component. We focus on a possible decoupling between electrons and ions temperatures. We then solve the energy transfer equations between electrons, ions and neutrals without assuming thermal equilibrium between the three gases (T_e <> T_i <> T_n)…
▽ More
We simulate the formation of a large X-ray cluster using a fully 3D hydrodynamical code coupled to a Particle-Mesh scheme which models the dark matter component. We focus on a possible decoupling between electrons and ions temperatures. We then solve the energy transfer equations between electrons, ions and neutrals without assuming thermal equilibrium between the three gases (T_e <> T_i <> T_n). We solve self-consistently the chemical equations for an hydrogen/helium primordial plasma without assuming ionization-recombination equilibrium. We find that the electron temperature differs from the true dynamical temperature by 20% at the Virial radius of our simulated cluster. This could lead marginally to an underestimate of the total mass in the outer regions of large X-ray clusters.
△ Less
Submitted 18 September, 1997;
originally announced September 1997.
-
Self Similar Spherical Collapse Revisited: a Comparison between Gas and Dark Matter Dynamics
Authors:
Romain Teyssier,
Jean-Pierre CHIEZE,
Jean-Michel Alimi
Abstract:
We reconsider the collapse of cosmic structures in an Einstein-de Sitter Universe, using the self similar initial conditions of Fillmore & Goldreich (1984). We first derive a new approximation to describe the dark matter dynamics in spherical geometry, that we refer to the "fluid approach". This method enables us to recover the self-similarity solutions of Fillmore & Goldreich for dark matter. W…
▽ More
We reconsider the collapse of cosmic structures in an Einstein-de Sitter Universe, using the self similar initial conditions of Fillmore & Goldreich (1984). We first derive a new approximation to describe the dark matter dynamics in spherical geometry, that we refer to the "fluid approach". This method enables us to recover the self-similarity solutions of Fillmore & Goldreich for dark matter. We derive also new self-similarity solutions for the gas. We thus compare directly gas and dark matter dynamics, focusing on the differences due to their different dimensionalities in velocity space. This work may have interesting consequences for gas and dark matter distributions in large galaxy clusters, allowing to explain why the total mass profile is always steeper than the X-ray gas profile. We discuss also the shape of the dark matter density profile found in N-body simulations in terms of a change of dimensionality in the dark matter velocity space. The stable clustering hypothesis has been finally considered in the light of this analytical approach.
△ Less
Submitted 3 April, 1997;
originally announced April 1997.
-
Gas and Dark Matter Spherical Dynamics
Authors:
Jean-Pierre CHIEZE,
Romain Teyssier,
Jean-Michel Alimi
Abstract:
We investigate the formation of spherical cosmological structures following both dark matter and gas components. We focus on the dynamical aspect of the collapse assuming an adiabatic, $γ= 5/3$, fully ionized primordial plasma. We use for that purpose a fully Lagrangian hydrodynamical code designed to describe highly compressible flows in spherical geometry. We investigate also a "fluid approach…
▽ More
We investigate the formation of spherical cosmological structures following both dark matter and gas components. We focus on the dynamical aspect of the collapse assuming an adiabatic, $γ= 5/3$, fully ionized primordial plasma. We use for that purpose a fully Lagrangian hydrodynamical code designed to describe highly compressible flows in spherical geometry. We investigate also a "fluid approach" to describe the mean physical quantities of the dark matter flow. We test its validity for a wide range of initial density contrast. We show that an homogeneous isentropic core forms in the gas distribution, surrounded by a self-similar hydrostatic halo, with much higher entropy generated by shock dissipation. We derive analytical expressions for the size, density and temperature of the core, as well as for the surrounding halo. We show that, unless very efficient heating processes occur in the intergalactic medium, we are unable to reproduce within adiabatic models the typical core sizes in X-ray clusters. We also show that, for dynamical reasons only, the gas distribution is naturally antibiased relative to the total mass distribution, without invoking any reheating processes. This could explain why the gas fraction increases with radius in very large X-ray clusters. As a preparation for the next study devoted to the thermodynamical aspect of the collapse, we investigate the initial entropy level required to solve the core problem in X-ray clusters.
△ Less
Submitted 3 April, 1997;
originally announced April 1997.
-
Adaptative Smooth Particle Hydrodynamics and Particle-Particle coupled codes: Energy and Entropy Conservation
Authors:
Arturo Serna,
Jean-Michel Alimi,
Jean-Pierre Chieze
Abstract:
We present and test a general-purpose code, called PPASPH, for evolving self-gravitating fluids in astrophysics, both with and without a collisionless component. In PPASPH, hydrodynamical properties are computed by using the SPH (Smoothed Particle Hydrodynamics) method while, unlike most previous implementations of SPH, gravitational forces are computed by a PP (Particle-Particle) approach. anot…
▽ More
We present and test a general-purpose code, called PPASPH, for evolving self-gravitating fluids in astrophysics, both with and without a collisionless component. In PPASPH, hydrodynamical properties are computed by using the SPH (Smoothed Particle Hydrodynamics) method while, unlike most previous implementations of SPH, gravitational forces are computed by a PP (Particle-Particle) approach. another important feature of this code is that hydrodynamics equations optionally include the correction terms appearing when $h(t,\br)$ is not constant. Our code has been implemented by using the data parallel programming model on CM5(CM). PPASPH has been applied to study the importance of adaptative smoothing correction terms on the entropy conservation. We confirm Hernquist's (1993) interpretation of the entropy violation observed in previous SPH simulations as a result of having neglected these terms. An improvement on the entropy conservation is not found by merely considering larger numbers of particles or different $N_S$ choices. The correct continuum description is only obtained if the $\bn h$ correction terms are included. Otherwise, the entropy conservation is always rather poor as compared to that found for the total energy.
△ Less
Submitted 24 November, 1995;
originally announced November 1995.