Commit 21de2cff authored by NUNEZ Arturo's avatar NUNEZ Arturo

Automatic commit vendredi 3 août 2018, 16:30:01 (UTC+0200)

parent b7b98cd0
......@@ -104,7 +104,7 @@
},
{
"cell_type": "code",
"execution_count": 30,
"execution_count": 67,
"metadata": {
"collapsed": false
},
......@@ -113,15 +113,15 @@
"name": "stdout",
"output_type": "stream",
"text": [
"With Delta_x = 152.587890625 as in Roskar 2013 we reproduce the value for p_j = n_h = 2.59866886906\n",
"and Teq = 3101.66245078\n"
"With Delta_x = 34 as in Roskar 2013 we reproduce the value for p_j = n_h = 19.23706849\n",
"and Teq = 1139.98873219\n"
]
}
],
"source": [
"# reproducing Roskar 2013\n",
"levelmax = 17\n",
"Delta_x = (20e6 / 2.**levelmax)\n",
"levelmax = 20\n",
"Delta_x = 34#36e6 / 2.**levelmax\n",
"\n",
"n_h, aux=5., 20.\n",
"i=0\n",
......@@ -131,15 +131,16 @@
" (G* m_mg*np.sqrt(aux)*(4*Delta_x*pctocm/100)**2)\n",
" n_h /= (m_p*1e6)#\n",
"\n",
"\n",
"Teq = 5000/np.sqrt(n_h)\n",
" \n",
"\n",
"print \"With Delta_x = {0} as in Roskar 2013 we reproduce the value for p_j = n_h = {1}\".format(Delta_x,n_h)\n",
"print \"and Teq = {0}\".format(Teq)"
]
},
{
"cell_type": "code",
"execution_count": 29,
"execution_count": 52,
"metadata": {
"collapsed": false
},
......@@ -148,8 +149,27 @@
"name": "stdout",
"output_type": "stream",
"text": [
"With Delta_x = 34.0595245361 as in Roskar 2013 we reproduce the value for p_j = n_h = 19.1922557107\n",
"and Teq = 1141.31885925\n"
"39.0\n",
"0.8670303260137321\n",
"1.7423551526078138\n",
"0.3961373743997542\n",
"0.28685888223305883\n",
"0.11847519410838325\n",
"0.06508122993995183\n",
"0.031033785686198457\n",
"0.01588765486721152\n",
"0.0078504068917792\n",
"0.003948466516169543\n",
"0.001968406033336382\n",
"0.0009856583875317458\n",
"0.0004924651718333048\n",
"0.00024632356898545565\n",
"0.00012313903592466217\n",
"6.157520475427508e-05\n",
"3.0786180635163585e-05\n",
"1.539344574734524e-05\n",
"With Delta_x = 34.0595245361 as in Roskar 2013 we reproduce the value for p_j = n_h = 19.1926988648\n",
"and Teq = 1141.30568281\n"
]
}
],
......@@ -158,10 +178,12 @@
"levelmax = 20\n",
"Delta_x = (35.714e6 / 2.**levelmax)\n",
"\n",
"n_h, aux=5., 20.\n",
"n_h, aux=5., 200.\n",
"i=0\n",
"while (np.abs(n_h-aux)/n_h ) > 1e-4:\n",
"while (np.abs(n_h-aux)/n_h ) > 1e-5:\n",
" print (np.abs(n_h-aux)/n_h )\n",
" aux = np.copy(n_h)\n",
" \n",
" n_h = (2*np.pi*k_b*1e4*np.sqrt(0.3))/\\\n",
" (G* m_mg*np.sqrt(aux)*(4*Delta_x*pctocm/100)**2)\n",
" n_h /= (m_p*1e6)#\n",
......@@ -322,12 +344,25 @@
},
{
"cell_type": "code",
"execution_count": null,
"execution_count": 57,
"metadata": {
"collapsed": true
"collapsed": false
},
"outputs": [],
"source": []
"outputs": [
{
"data": {
"text/plain": [
"5477.225575051661"
]
},
"execution_count": 57,
"metadata": {},
"output_type": "execute_result"
}
],
"source": [
"1e4*np.sqrt(.3)"
]
},
{
"cell_type": "code",
......
Markdown is supported
0% or
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment