group8
Differences
This shows you the differences between two versions of the page.
| Both sides previous revision Previous revision Next revision | Previous revision | ||
|
group8 [2014/06/05 23:01] bartl update June 5 |
group8 [2014/06/06 15:47] talwar |
||
|---|---|---|---|
| Line 12: | Line 12: | ||
| - parametric study: What is the kill radius at given explosion energy and solar wind strength? | - parametric study: What is the kill radius at given explosion energy and solar wind strength? | ||
| - production of the "Battle of the Winds" movie blockbuster | - production of the "Battle of the Winds" movie blockbuster | ||
| + | |||
| + | ==== Code ==== | ||
| + | A commented review of our code modifications can be found under [[group8_code|Inside the Code]]. | ||
| ==== Logbook ==== | ==== Logbook ==== | ||
| Line 35: | Line 38: | ||
| Date : 06/04/2014 | Date : 06/04/2014 | ||
| - | - Changed the 2D inputs to more natural units ([L]=1 AU, [T]=1 yr, [M]=10^12 g) because VH-1 would crash when using the cgs system for our problem. | + | - Changed the 2D inputs to more natural units ([L]=1 AU, [T]=1 yr, [M]=10^12 g=1 Mt) because VH-1 would crash when using the cgs system for our problem. |
| - Fixed an error in our input densities. | - Fixed an error in our input densities. | ||
| - Fixed the bug of density accumulating on the symmetry axis (we use cylindrical coordinates) by setting the density of the two innermost gridpoints in each column to the value of the third. This is a quick and dirty fix and unphysical! However, we think the problem is caused by some boundary condition/ghost zone issue that we were not able to resolve after spending several hours on it. We thank Richard for suggesting this fix! | - Fixed the bug of density accumulating on the symmetry axis (we use cylindrical coordinates) by setting the density of the two innermost gridpoints in each column to the value of the third. This is a quick and dirty fix and unphysical! However, we think the problem is caused by some boundary condition/ghost zone issue that we were not able to resolve after spending several hours on it. We thank Richard for suggesting this fix! | ||
| Line 42: | Line 45: | ||
| Date : 06/05/2014 | Date : 06/05/2014 | ||
| - | - Finalized our input parameters | + | - Finalized our input parameters based on [[http://dx.doi.org/10.1086/523622|Brian's paper]] (density of injected material set to approximately produce 4.5e4 Mt/AU^3 at 1AU): <code fortran> |
| - | - Generated high-res simulations (10x4 AU, 750x250 gridpoints) for the presentation | + | ! ambient material |
| + | pamb = 3.27e4 !pressure | ||
| + | ramb = 5.60e3 !density | ||
| + | |||
| + | ! solar wind properties | ||
| + | zinj = 4 ! z-coordinate of injection disk | ||
| + | radinj = 0.25 ! 1e18 ! injection disk radius | ||
| + | widthinj = 1.6 ! width of injection disk (radinj/widthinj <= r <= radinj) | ||
| + | vinj = 88.75 ! velocity of injected material | ||
| + | pinj = pamb ! pressure of injected material | ||
| + | rinj = 20*4.54e4 ! density of injected material | ||
| + | |||
| + | ! SNR shock | ||
| + | ncycleshock = 4500 !number of cycles to wait for solar wind to settle before shock | ||
| + | vshock = 128.6 !velocity of shock material | ||
| + | pshock = 1.23e8 !pressure of shock material | ||
| + | rshock = 2.24e4 !density of shock material | ||
| + | |||
| + | </code> | ||
| + | |||
| + | - 1D-Simulation use some settings of Cartesian coordinate, zero-gradient boundaries, planar and constant solar & supernova winds | ||
| + | |||
| + | Velocity Profile | ||
| + | |||
| + | {{::output_5yzrfo.gif?direct&200|high-res simulations}} | ||
| + | |||
| + | Density Profile | ||
| + | |||
| + | {{::output_14fvpq.gif?direct&200|high-res simulations}} | ||
| + | |||
| + | - Generated high-res simulations (10x4 AU, 750x250 gridpoints) for the presentation (click on the image to show gif (33MB): | ||
| + | {{::2dwind1.gif?direct&200|high-res simulations}} | ||
| - Did a parameter study on a reduced grid (6x2 AU, 300x100 gridpoints) to determine the dependence of the penetration radius on the densities of the SNR and solar wind. | - Did a parameter study on a reduced grid (6x2 AU, 300x100 gridpoints) to determine the dependence of the penetration radius on the densities of the SNR and solar wind. | ||
| - Figured out how to create animated gifs / videos of the NetCDF files using [[https://wci.llnl.gov/codes/visit/|VisIt]] | - Figured out how to create animated gifs / videos of the NetCDF files using [[https://wci.llnl.gov/codes/visit/|VisIt]] | ||
| - | - Created our presentation | + | - Created our presentation : Here's the link : https://drive.google.com/file/d/0B5VQV2FlDo8wQllHdG1NakVZV3c/edit?usp=sharing |
group8.txt · Last modified: 2014/06/06 15:47 by talwar
Page Tools
Except where otherwise noted, content on this wiki is licensed under the following license: CC Attribution-Share Alike 4.0 International
