- Mail actions: [ respond to this message ] [ mail a new topic ]
- Contemporary messages sorted: [ by date ] [ by thread ] [ by subject ] [ by author ] [ by messages with attachments ]

From: Yves Revaz <yves.revaz_at_obspm.fr>

Date: Thu, 01 Feb 2007 18:17:28 +0100

Dear Gadget list,

I'm trying to understand some points in the comoving

implementation of the hydro part of Gadget2.

In "timestep.c" when ComovingIntegrationOn is true,

dt_entr = (tend - tstart) * All.Timebase_interval;

which in fact corresponds to dt_entr = da/a (where a is the scale factor).

Then, the evolution of the entropy is (as for non-comoving simulations)

SphP[i].Entropy += SphP[i].DtEntropy * dt_entr;

SphP[i].DtEntropy is computed in "hydra.c", where we find in

/* do final operations on results */ :

SphP[i].DtEntropy *= GAMMA_MINUS1 / (hubble_a2 * pow(SphP[i].Density,

GAMMA_MINUS1));

which means, that :

SphP[i].DtEntropy)comoving = DEntropy/da

= DEntropy/dt * 1/(H(a) a^2)

= SphP[i].DtEntropy)non-comoving * 1/(H(a) a^2)

Here, I expected to have :

DtEntropy/da = DtEntropy/dt * 1/(H(a))

in order to have then :

SphP[i].DtEntropy * dt_entr = DEntropy/dt * 1/(H(a)) * da/a

= DEntropy

because dt*H(a)*a = da

Probably I missed a point here.

It is also difficult to understand the fac_mu correction in "hydra.c":

1) 0.0001 * soundspeed_i / SphP[i].Hsml / fac_mu

with fac_mu = a^(3/2(GAMMA - 1))/a

Here I expected to find

fac_mu = a^(1/2(GAMMA - 1))/a

in order to write soundspeed_i / SphP[i].Hsml in non-comoving units.

2) mu_ij = fac_mu * vdotr2 / r

is also unclear to me

Another question conserns the

dt_hydrokick factor which corresponds to the integral

Int_a1^a2 da/( H(a)*a*a^3(gamma-1) )

If, following Quinn et al. 97, the origin of dt_gravkick is clear,

I fail to understand the origin of dt_hydrokick.

Thanks is advance.

Date: Thu, 01 Feb 2007 18:17:28 +0100

Dear Gadget list,

I'm trying to understand some points in the comoving

implementation of the hydro part of Gadget2.

In "timestep.c" when ComovingIntegrationOn is true,

dt_entr = (tend - tstart) * All.Timebase_interval;

which in fact corresponds to dt_entr = da/a (where a is the scale factor).

Then, the evolution of the entropy is (as for non-comoving simulations)

SphP[i].Entropy += SphP[i].DtEntropy * dt_entr;

SphP[i].DtEntropy is computed in "hydra.c", where we find in

/* do final operations on results */ :

SphP[i].DtEntropy *= GAMMA_MINUS1 / (hubble_a2 * pow(SphP[i].Density,

GAMMA_MINUS1));

which means, that :

SphP[i].DtEntropy)comoving = DEntropy/da

= DEntropy/dt * 1/(H(a) a^2)

= SphP[i].DtEntropy)non-comoving * 1/(H(a) a^2)

Here, I expected to have :

DtEntropy/da = DtEntropy/dt * 1/(H(a))

in order to have then :

SphP[i].DtEntropy * dt_entr = DEntropy/dt * 1/(H(a)) * da/a

= DEntropy

because dt*H(a)*a = da

Probably I missed a point here.

It is also difficult to understand the fac_mu correction in "hydra.c":

1) 0.0001 * soundspeed_i / SphP[i].Hsml / fac_mu

with fac_mu = a^(3/2(GAMMA - 1))/a

Here I expected to find

fac_mu = a^(1/2(GAMMA - 1))/a

in order to write soundspeed_i / SphP[i].Hsml in non-comoving units.

2) mu_ij = fac_mu * vdotr2 / r

is also unclear to me

Another question conserns the

dt_hydrokick factor which corresponds to the integral

Int_a1^a2 da/( H(a)*a*a^3(gamma-1) )

If, following Quinn et al. 97, the origin of dt_gravkick is clear,

I fail to understand the origin of dt_hydrokick.

Thanks is advance.

-- (o o) --------------------------------------------oOO--(_)--OOo------- Yves Revaz Lerma Batiment A Tel : ++ 33 (0) 1 40 51 20 79 Observatoire de Paris Fax : ++ 33 (0) 1 40 51 20 02 77 av Denfert-Rochereau e-mail : yves.revaz_at_obspm.fr F-75014 Paris Web : http://obswww.unige.ch/~revaz/ FRANCE ----------------------------------------------------------------Received on 2007-02-01 18:17:29

*
This archive was generated by hypermail 2.3.0
: 2023-01-10 10:01:30 CET
*