å®å®ã«æ¨ã¦ããå¢æ®ããæ é¥ é ã¯ã©ãããæã«ãªãã®ãããã©ãã¯ãã¼ã«ããä¸æ§åæããè¶ æ°æããç½è²ç®æãï¼ããã®ã¤ã¥ã http://d.hatena.ne.jp/ita/20081223/p1
æ°ä½ã¯å§ç¸®ããã¨ãã®ä»äºãç±ã«å¤æãããããåºä½ã®å ´åã¯åçºããååå士ãç¡çããæ¼ãã¤ã¶ãã®ã«ä½¿ãããããªã®ã§æ¸©åº¦ã¯ä¸å®ã¨ãã¦ç¡è¦ã
æ é¥
é ã®å¯åº¦ã¨å§åã®é¢ä¿ã¯ã常å§ã§ã¯æ°´ã®ãã¼ã¿ã«åããããã£ã¨å§ç¸®ããã¨ç½è²ç®æã®ããã«é»åã®ç¸®éå§ãåºã¦ããããã«é©å½ã«æ±ºããããããªæã
p=ACONST *(rho/RHO0-1) + FERMIC*pow(rho/RHO0-1.0, 1.66666)
ããã§åç¹ã§å§åå¾é
ï¼å¯åº¦ï¼éåå é度ã¨ãªããããªå¹³è¡¡ç¶æ
ãè¨ç®ããã
ã¾ãæ é¥
é ã¯é£ã¹ãã¨å¢æ®ããªããªãã®ã§ãé¥
é ä¸åã«1kgã®éããä¹ãããããã®å§å(ç´4000 Pa)ããããã¨æ½°ãã¦å¢æ®ãåæ¢ããã¨ä»®å®ãããã¾ãå°ä¸æ·±ãã«ããé¥
é ã¯å¨ãã«ç©ºéããªãã¦å¢æ®ã§ããªãã®ã§ãæ ¹æ ã¯ãªããã©å°è¡¨ããï¼ï¼ï½ï½ã¾ã§ã®é¥
é ã®ã¿ï¼åã«ä¸åº¦å¢æ®ããã¨ä»®å®ããã
ãããã®ä»®å®ã®å
ã«é¥
é æã®è³ªéãåå¾ã表é¢ã§ã®éåãæéã¨å
±ã«ã©ãå¤ããããè¨ç®ãå³ã¯ãã¹ã¦å°çã®å¤ãï¼ã¨ãã¦ããã
éä¸ããç½è²ç®æã®ç¹å¾´ã§ããã質éãå¢ããã»ã©åå¾ãå°ãããªãã¨ããç¾è±¡ãè¦ããããããããã£ã¨é²ãã¨è¶ æ°æççºãèµ·ããããã質éã太é½ã®1.4åç¨åº¦ã®ãã£ã³ãã©ã»ã«ã¼ã«éçã¨å¼ã°ããã¨ãããè¶ ããã¨ç¸®éã¨ãã«ã®ã¼ãé«ããªãããã¦é»åã¨é½åãã²ã£ã¤ãã¦ä¸æ§åã«ãªã£ãã»ãããã·ãã¨ããç¶æ ã«ãªãä¸æ§åæã«ãªã£ããããããããããã«éããåã«å°è¡¨ã§ã®éåãå¼·ããªã£ã¦è¡¨é¢ã®æ é¥ é ãå ¨ã¦æ½°ãã¦æã®é²åãåæ¢ãããã ãããæ°ä¸å¹´ã®æªæ¥ã
æè¡ç詳細
å¨æããã³ã·ã£ã«ä¸ã®ã¨ãã«ã®ã¼åºæç¶æ ã¯ããããã®å®çããexp(ikx)ï¼ãªã«ããåºæå¤ã¯ï¼ãªã«ããã®åºæå¤ï¼ï¼k^2ã§æ³¢æ°ã®æ大å¤ã¯Ï/å¨æãé»åãä¸ããè©°ããã¨åè¨ã§1/ååéé^2 ã®ãªã¼ãã¼ã§ã¯ã¼ãã³ãä¸åããå¯åº¦ã§è¡¨ãã¨å¯åº¦ã®2/3ä¹ã
平衡ç¶æ ã¯ã¾ãä¸å¿ã®å¯åº¦rhoãè¨å®ãã
r=0, m=0 while(1) { r += dr m += 4*pi*r*r*dr*rho g = G*m/(r*r) // exact in 3D //dP/dx = dP/drho * drho/dx and rho*g= -dP/dx // -> drho/dx = -rho*g / (dP/drho) rho -= g * rho / dPdrho(rho) * dr if (P(rho)<0.0) break }
ã®ããã«å¤å´ã«åãã£ã¦è§£ãã¦ããå§åãï¼ã®å¤ç¸ã§çµäºã次ã«rhoã®åæå¤ãå°ãå¢ããã¦ã¾ã解ããå ¨è³ªéã®ååã¨ã®å·®åã質éå¢å é度ã§å²ã£ã¦æéã«å ç®ãããã¨ãç¹°ãè¿ãã