Changes between Version 46 and Version 47 of ModulesOnAstroBear


Ignore:
Timestamp:
04/20/17 12:52:56 (8 years ago)
Author:
Jonathan
Comment:

Legend:

Unmodified
Added
Removed
Modified
  • ModulesOnAstroBear

    v46 v47  
    101101Anything that is defined in the {{{InfoDef}}} type is now available.  For example, {{{q}}} need not be defined...just reference it by {{{Info%q}}}.
    102102
    103 There are two major data arrays in {{{InfoDef}}}: the {{{q}}} array and the {{{aux}}} array.  {{{q}}} holds the volume averaged data and is used by all AstroBEAR simulations while {{{aux}}} holds face averaged data and is used only for MHD when nDim > 1.  Note that volume averaged data and cell-centered data are often used interchangeably, but there is an important distinction.  Take for instance a simple function f(x) defined on the interval [0:h].  The average of f(x) over the interval is
    104 
    105  f(0)+d,,1,,f(0)h/2+d,,2,,f(0)h^2^/6+...
    106 
    107 where the cell centered value is
    108 
    109  f(0)+d,,1,,f(0)h/2+d,,2,,f(0)h^2^/8+...
     103There are two major data arrays in {{{InfoDef}}}: the {{{q}}} array and the {{{aux}}} array.  {{{q}}} holds the volume averaged data and is used by all AstroBEAR simulations while {{{aux}}} holds face averaged data and is used only for MHD when nDim > 1.  Note that volume averaged data and cell-centered data are often used interchangeably, but there is an important distinction.  Take for instance a simple function $f(x)$ defined on the interval $[-h/2:h/2]$.  Applying a Taylor expansion to $f$ and finding the average of $f(x)$ gives
     104
     105$\bar{f}=f(0)+\frac{1}{24} f''(0) h^2 + ...$
     106
     107where the cell centered value is $f(0)$
    110108
    111109so the cell centered value is second order accurate for the volume average and usually is a quick way to estimate the volume average.  However if the function has large 2nd derivatives (or higher) this can lead to large errors in the volume average.  This is often apparent when modeling discontinuities along curved boundaries.  There are two ways to handle this problem: