Changes between Version 37 and Version 38 of PhysicsDataExplained


Ignore:
Timestamp:
10/11/12 19:22:29 (12 years ago)
Author:
ehansen
Comment:

Legend:

Unmodified
Added
Removed
Modified
  • PhysicsDataExplained

    v37 v38  
    101101
    102102[[BR]]
    103 == Cylindrical Symmetry (2.5D) ==
    104 This section gives some brief details on how you can take advantage of cylindrical coordinates if your problem has such symmetry.
     103== Cylindrical Coordinates with Axial Symmetry (2.5D) ==
     104This section gives some brief details on how you can take advantage of cylindrical coordinates with axial symmetry if your problem has such symmetry.
    105105
    106 * '''''{{{iCylindrical:}}}'''''  Controls whether astrobear is to use cylindrical coordinates.  The default is 0 which leaves everything in Cartesian, 1 is for no angular momentum, and 2 is for with angular momentum.
     106* '''''{{{iCylindrical:}}}'''''  Controls whether astrobear is to use cylindrical coordinates with axial symmetry.  The default is 0 which leaves everything in Cartesian, 1 is for no angular momentum, and 2 is for with angular momentum.
    107107
    108 Basically when iCylindrical is not 0, astrobear changes the differential operators and how each axis is treated.  In other words, astrobear switches to cylindrical coordinates.  The x-axis becomes the radial axis, the y-axis becomes the cylindrical z-axis, and the z-axis is nothing since this is only implemented for ndim = 2.
     108Basically when iCylindrical is not 0, astrobear changes the differential operators and how each axis is treated.  In other words, astrobear switches to a cylindrical coordinate system with axial symmetry (no phi dependence).  The x-axis becomes the radial axis, the y-axis becomes the cylindrical z-axis, and the z-axis is nothing since this is only makes sense when ndim = 2.
    109109
    110 Furthermore, if your problem has rotation then iCylindrical should be set to 2.  This will change the ivz component to the angular momentum component.  Otherwise, if iCylindrical=1 then there is only ivx (radial) and ivy (cylindrical z).
     110Furthermore, if your problem has rotation then iCylindrical should be set to 2.  This will change the ivz component to the angular momentum component (phi).  This does not mean that there is phi dependence, only that there is motion in the phi direction.  Otherwise, if iCylindrical=1 then there is only ivx (radial) and ivy (cylindrical z).
    111111
    112112One other important thing to note...if iCylindrical is not 0, then mydim (a global parameter in physics declarations) will be set to 3, otherwise it is set to ndim.