Changes between Version 2 and Version 3 of AstroBearProjects/WindDiskInteraction


Ignore:
Timestamp:
09/08/11 12:11:49 (13 years ago)
Author:
Jonathan
Comment:

Legend:

Unmodified
Added
Removed
Modified
  • AstroBearProjects/WindDiskInteraction

    v2 v3  
     1[[PageOutline]]
     2
    13 = Wind Disk Interacion =
    24
    35= Conservation of Angular Momentum =
     6
     7== 2D Results ==
    48  While the code conserves momentum, it does not explicitly conserve angular momentum.  To test this we setup a 2D disk around a point gravity source with a density contrast of 100.  Pictured below is the angular momentum.  The left panel has a truncated color scale to show the effects of periodic boundaries as well as the transport of angular momentum from the disk into the ambient.  Also shown are contours of angular velocity corresponding to .5, 1, and 2 rotations over the course of the simulation.  The right panel is a zoom of the disk itself without a truncated color scale.  Here is the [attachment:AngMom.gif movie]
    59
    610[[Image(AngMom0100.png, width=800)]]
    711
     12
     13
     14
     15
     16
     17
     18
     19
     20
     21
     22
     23
     24
     25
     26
     27
     28
     29
     30
     31
     32
     33
     34
     35
     36
    837We also did an explicit query of the total angular momentum as a function of time.  Over the course of the simulation the total fractional change in angular momentum is only .00081 or less then .1% It is likely that the loss of angular momentum is due to numerical dissipation at the center of the disk, although with levels so small it is difficult to rule out boundary effects.  Although the symmetry of the problem should keep the net angmom transport from the boundaries to zero.
    938
    10 [[Image(AngMom.png, width=800)]]
     39[[Image(AngMom.png, width=400)]]
     40
     41
     42
     43
     44
     45
     46
     47
     48
     49
     50
     51
     52
     53
     54
     55
     56
     57
     58
     59
     60
     61
     62
     63
     64
     65
     66
     67
     68
     69
     70
     71
     72
     73
     74
     75== 3D Results ==
     76
     77We then extended the simple disk into 3D neglecting pressure support in the z direction.  First we used no gravitational softening.  The large diffusion of angular momentum at the origin produces heat that leads to jets.
     78
     79||  [attachment:3DDensVel.gif movie]  ||  [attachment:3DContours.gif movie]  ||
     80||  [[Image(3DDensVel0100.png, width=400)]]  ||  [[Image(3DContours0100.png, width=400)]]  ||
     81
     82
     83
     84
     85
     86
     87
     88
     89
     90
     91
     92
     93
     94
     95
     96
     97
     98
     99
     100
     101
     102
     103
     104Next we used gravitational softening at a radius of 20 computational units or 8 cells.  This suppressed any jet
     105
     106||  [attachment:3DDensVelSL20.gif movie]  ||  [attachment:3DContourSL20.gif movie]  ||
     107||  [[Image(3DDensVelSL200100.png, width=400)]]  ||  [[Image(3DContourSL200100.png, width=400)]]  ||
     108
     109
     110
     111
     112We then increased the resolution by 2 and decreased the softening radius to 5 computational units or 4 cells.  No jet forms although there is a hot halo.
     113
     114||  [attachment:3DDensVelSL10Res64.gif movie]  ||  [attachment:3DContourSL10Res64.gif movie]  ||
     115||  [[Image(3DDensVelSL10Res640100.png, width=400)]]  ||  [[Image(3DContourSL10Res640100.png, width=400)]]  ||
     116
     117
     118We then compared the 3D results to a 2.5D run with the same setup. 
     119
     120||  [attachment:2D3DSL5Res64.gif movie] (left panel is 3D, right panel is 2.5D)  ||  [attachment:2D3DSL5Res64LineOuts.2.gif movie] Red is density and blue is Angular Momentum.  Solid is 2.5D and dashed is 3D  ||
     121||  [[Image(2D3DSL5Res640100.png, width=400)]]  ||  [[Image(2D3DSL5Res64LineOuts0100.2.png, width=400)]]  ||
     122
     123
     124And finally compared the total z component of the angular momentum of the 2.5D run (red) with the 3D runs with no gravitational softening (blue) and with 4 cells of gravitational softening (green).  They are all within a few percent and the boundaries may be playing a small role as angular momentum can be lost/gained from the reflecting boundaries.  However the increased loss of angular momentum for the 3D run with no softening (blue) is consistent with the presence of a jet and the diffusion of angular momentum at the origin.
     125
     126[[Image(AngMomComparison.png)]]
     127
     128
     129