cwpearson/astaroth
1
0
Fork 0
Code Issues Pull Requests Actions Packages Projects Releases Wiki Activity

Had some rough draft of accretion criterion.

Browse Source
This commit is contained in:
JackHsu
2019-08-13 18:31:57 +08:00
parent 927a5e74f2
commit 0587d08b01
1 changed files with 4 additions and 1 deletions
Download Patch File Download Diff File Expand all files Collapse all files

5
acc/mhd_solver/stencil_process.sps
View File

@@ -118,7 +118,10 @@ accretion_profile(int3 globalVertexIdx, in Scalar lnrho){
// multiplying the truelove density by a wave function to avoid step-function like accretion profile. // multiplying the truelove density by a wave function to avoid step-function like accretion profile.
const Scalar weight = exp(-(accretion_distance/profile_range)); const Scalar weight = exp(-(accretion_distance/profile_range));
const Scalar rate = truelove_density(lnrho); // const Scalar rate = truelove_density(lnrho);
const Scalar B = Scalar(3.0);
const Scalar k = Scalar(2.0);
const Scalar rate = B * (pow(sink_mass, k) / (dsx * dsy * dsz));
accretion_density = weight * rate; accretion_density = weight * rate;
return accretion_density; return accretion_density;