astrobear - Blog
https://bluehound2.circ.rochester.edu/astrobear/blog
About blog postsen-USTrac 1.4.1Meeting Update-08/05/13 - ErinielambridMon, 05 Aug 2013 18:02:49 GMT
https://bluehound2.circ.rochester.edu/astrobear/blog/elambrid08052013
https://bluehound2.circ.rochester.edu/astrobear/blog/elambrid08052013<p>
-Running a large polytrope run atm (Gmx =128<sup>3, gxbounds=-2-22222, t_final=.04d0 with timescale ~ 7.6 *10</sup>6
</p>
<p>
-lineout at 2.5 sound crossing time
</p>
<p>
<a href="https://bluehound2.circ.rochester.edu/astrobear/attachment/blog/elambrid08052013/2.5_sound_crossing_time_lineout.gif" style="padding:0; border:none"><img crossorigin="anonymous" src="https://bluehound2.circ.rochester.edu/astrobear/raw-attachment/blog/elambrid08052013/2.5_sound_crossing_time_lineout.gif" width="50%" /></a>
</p>
<p>
-pseudo
<a href="https://bluehound2.circ.rochester.edu/astrobear/attachment/blog/elambrid08052013/2.5_sound_crossing_time_pseudo.gif" style="padding:0; border:none"><img crossorigin="anonymous" src="https://bluehound2.circ.rochester.edu/astrobear/raw-attachment/blog/elambrid08052013/2.5_sound_crossing_time_pseudo.gif" width="50%" /></a>
</p>
<p>
No explosion persay, but definitley some odd boundary effects. This is the longest run I have done at this resolution.
</p>
<p>
-check out <a class="wiki" href="https://bluehound2.circ.rochester.edu/astrobear/wiki/u/Polytropes">u/Polytropes</a> - for it is updated
</p>
Meeting update - Erini 7/29/13elambridMon, 29 Jul 2013 19:35:58 GMT
https://bluehound2.circ.rochester.edu/astrobear/blog/elambrid07292013
https://bluehound2.circ.rochester.edu/astrobear/blog/elambrid07292013<p>
<a href="https://bluehound2.circ.rochester.edu/astrobear/attachment/blog/elambrid07292013/n3_646464_-11_.005_extrap.gif" style="padding:0; border:none"><img crossorigin="anonymous" src="https://bluehound2.circ.rochester.edu/astrobear/raw-attachment/blog/elambrid07292013/n3_646464_-11_.005_extrap.gif" width="99%" /></a>
</p>
Meeting Update 7/23/13 -ErinielambridTue, 23 Jul 2013 19:14:05 GMT
https://bluehound2.circ.rochester.edu/astrobear/blog/elambrid07232013
https://bluehound2.circ.rochester.edu/astrobear/blog/elambrid07232013<ul><li>polytrope wiki current version <a class="wiki" href="https://bluehound2.circ.rochester.edu/astrobear/wiki/u/Polytropes">u/Polytropes</a>
</li></ul><ul><li>I have been working on my ticket <a href="https://bluehound2.circ.rochester.edu/astrobear/query?id=" title="Tickets ">#303</a>
</li></ul><ul><li>I had some confusion in terms of various global.data parameters that I have cleaned through, thus the following two gifs represent my issues
</li></ul><p>
n=3 res:80^3 comp cell size: 2's - notice the asymmetry
<a href="https://bluehound2.circ.rochester.edu/astrobear/attachment/blog/elambrid07232013/n3_808080_-22_.001.gif" style="padding:0; border:none"><img crossorigin="anonymous" src="https://bluehound2.circ.rochester.edu/astrobear/raw-attachment/blog/elambrid07232013/n3_808080_-22_.001.gif" width="80%" /></a>
n=1.5 res:80^3 comp cell size: 2's
<a href="https://bluehound2.circ.rochester.edu/astrobear/attachment/blog/elambrid07232013/n1.5_808080_-22_.001.gif" style="padding:0; border:none"><img crossorigin="anonymous" src="https://bluehound2.circ.rochester.edu/astrobear/raw-attachment/blog/elambrid07232013/n1.5_808080_-22_.001.gif" width="80%" /></a>
</p>
<p>
-beginning to think that some of the issues might actually be with the numerics
<a class="ext-link" href="http://adsabs.harvard.edu/full/2000PASJ...52..217H"><span class="icon"></span>http://adsabs.harvard.edu/full/2000PASJ...52..217H</a>
</p>
<p>
Horedt lists that the density ratio is a dependent criteria in the stability equation for embedded polytropes - perhaps the density ratio being used is too high, e-mailed prof blackman and perhaps he can provide some insight on the usefulness of the paper
</p>
Erini - Meeting Update - July 1st 2013elambridMon, 01 Jul 2013 19:20:27 GMT
https://bluehound2.circ.rochester.edu/astrobear/blog/elambrid07012013
https://bluehound2.circ.rochester.edu/astrobear/blog/elambrid07012013<ul><li>updated a couple of wiki pages <a class="wiki" href="https://bluehound2.circ.rochester.edu/astrobear/wiki/LE_Module">LE_Module</a>
</li></ul><ul><li>polytrope is on grid<br />
</li><li>Parameters: <br />
grid resolution: 64 cubed <br />
domain: -2d0 -2d0 2d0 2d0 2d0 <br />
final times used: (.001, .01)<br />
frames used: (10,20)<br />
timescale: 7655157.60438493<br />
</li></ul><p>
→ so for the .01 final time, we have 76551.57…yikes too small <br />
→ so for the .01 t final time and 20 frames case we get the following
<a href="https://bluehound2.circ.rochester.edu/astrobear/attachment/blog/elambrid07012013/rho_070113.gif" style="padding:0; border:none"><img alt="density" crossorigin="anonymous" src="https://bluehound2.circ.rochester.edu/astrobear/raw-attachment/blog/elambrid07012013/rho_070113.gif" title="density" width="50" /></a>
<a href="https://bluehound2.circ.rochester.edu/astrobear/attachment/blog/elambrid07012013/pressure_070113.gif" style="padding:0; border:none"><img alt="pressure" crossorigin="anonymous" src="https://bluehound2.circ.rochester.edu/astrobear/raw-attachment/blog/elambrid07012013/pressure_070113.gif" title="pressure" width="50" /></a>
</p>
Erini - Meeting Update- June 24th 2013elambridMon, 24 Jun 2013 19:18:05 GMT
https://bluehound2.circ.rochester.edu/astrobear/blog/elambrid06242013
https://bluehound2.circ.rochester.edu/astrobear/blog/elambrid06242013<ul><li>having scaling errors with polytrope profile
</li><li>I initially got it to produce a first cut off polytropic profile (domain size too small) that agreed with literature for the well documented case of n=3
</li><li>upon trying to fix the scaling errors - I am now getting NANS
</li></ul><ul><li>still chugging away, but ever closer each day
</li></ul>Meeting Update June-17th ErinielambridMon, 17 Jun 2013 19:28:23 GMT
https://bluehound2.circ.rochester.edu/astrobear/blog/elambrid06172013
https://bluehound2.circ.rochester.edu/astrobear/blog/elambrid06172013<ul><li>Got the code to compile
</li></ul><ul><li>It's freezing and doesn't run - trying to figure out where it's crashing!
</li></ul><ul><li>So here is what I'm trying to change in terms of initializing:
</li></ul><pre class="wiki">putting the profile object on the ambient object
</pre><pre class="wiki">having the ambient density and pressure be the boundary density and pressure for the polytrope... ideas?
</pre>Meeting Update June 10th - ErinielambridMon, 10 Jun 2013 18:58:19 GMT
https://bluehound2.circ.rochester.edu/astrobear/blog/elambrid06102013
https://bluehound2.circ.rochester.edu/astrobear/blog/elambrid06102013<p>
Almost done with polytrope stuff<br />
</p>
<p>
I reformatted everything so that input params were central density, Mass and polytropic index as opposed to central density, star radius and polytropic index…but i have now run into an issue
<br />
</p>
<p>
Trying to figure out a way to approximate cell length - I was previously difining it as: cell length = R_outer/npoints <br />
</p>
<p>
Code finds R_outer AFTER the integration, any ideas as to how to find a good bound for cell length without radius being an input?
</p>
Erini's Update -4/23/13elambridTue, 23 Apr 2013 17:15:37 GMT
https://bluehound2.circ.rochester.edu/astrobear/blog/elambrid04232013
https://bluehound2.circ.rochester.edu/astrobear/blog/elambrid04232013<ul><li>still working on polytrope module
</li><li>almost done writing up generalized polytropes for wiki
</li></ul>Erini's Update 4/16/2013elambridTue, 16 Apr 2013 19:39:20 GMT
https://bluehound2.circ.rochester.edu/astrobear/blog/elambrid04162013
https://bluehound2.circ.rochester.edu/astrobear/blog/elambrid04162013<ul><li>finished LE_module
</li><li>made a project page for it
</li></ul>Erini's Update 4/1/13elambridMon, 01 Apr 2013 05:18:53 GMT
https://bluehound2.circ.rochester.edu/astrobear/blog/elambrid04012013
https://bluehound2.circ.rochester.edu/astrobear/blog/elambrid04012013<ul><li>i have made a test polytrope code, it is located in my test directory if anyone is interested
</li></ul><p>
-i have tried to upload a couple of images, but embarrassingly enough I cannot get it to work on the wiki
</p>
<ul><li>my dimensionless total radius matches to at least 4 decimal places to Chandrasekhar's table and current tables
</li></ul>Erini's Update 3/25-2012elambridMon, 25 Mar 2013 05:53:30 GMT
https://bluehound2.circ.rochester.edu/astrobear/blog/elambrid03252013
https://bluehound2.circ.rochester.edu/astrobear/blog/elambrid03252013<p>
-met with Erica, and Professor Blackman, and Jonathon
</p>
<ul><li> coding numerically the lane-emden equation as opposed to the approximate analytic paper I initially started with
</li><li> Jonathon led me through using various read-in profiles, which I will ultimately need to do with lane emden module
<ul><li>option of using stellar evolution code data from public tables
</li></ul></li></ul><p>
</p>
<p>
</p>
Erini's Update -3/18/13elambridMon, 18 Mar 2013 05:08:40 GMT
https://bluehound2.circ.rochester.edu/astrobear/blog/elambrid03182013
https://bluehound2.circ.rochester.edu/astrobear/blog/elambrid03182013<p>
Began coding the analytic approximate solution for the general case of the Lane Emden equation, I am using the solution derived here:<a class="ext-link" href="http://articles.adsabs.harvard.edu/cgi-bin/nph-iarticle_query?db_key=AST&bibcode=1996MNRAS.281.1197L&letter=0&classic=YES&defaultprint=YES&whole_paper=YES&page=1204&epage=1204&send=Send+PDF&filetype=.pdf"><span class="icon"></span>F.K.Liu's Polytropic Gas Spheres: An Analytic Solution to the Lane Emden Equation</a>(this is the same paper Erica used for her isothermal case, n= \(\infty\)).
</p>