The SKIRT project
advanced radiative transfer for astrophysics
AdaptiveMeshMedium Class Reference

#include <AdaptiveMeshMedium.hpp>

Inheritance diagram for AdaptiveMeshMedium:

Protected Member Functions

 AdaptiveMeshMedium ()
AdaptiveMeshSnapshotadaptiveMesh () const override
SnapshotcreateAndOpenSnapshot () override
Protected Member Functions inherited from MeshMedium
 MeshMedium ()
const Boxdomain () const
void setupSelfBefore () override
Protected Member Functions inherited from ImportedMedium
 ImportedMedium ()
 ~ImportedMedium ()
void setupSelfAfter () override
Protected Member Functions inherited from Medium
 Medium ()
Protected Member Functions inherited from SimulationItem
 SimulationItem ()
virtual bool offersInterface (const std::type_info &interfaceTypeInfo) const
Protected Member Functions inherited from Item
 Item ()
Protected Member Functions inherited from SiteListInterface
 SiteListInterface ()
Protected Member Functions inherited from AdaptiveMeshInterface
 AdaptiveMeshInterface ()

Private Types

using BaseType
using ItemType

Private Attributes

AdaptiveMeshSnapshot_adaptiveMeshSnapshot

Friends

class ItemRegistry

Additional Inherited Members

Public Types inherited from MeshMedium
enum class  MassType : int {
  MassDensity , Mass , MassDensityAndMass , NumberDensity ,
  Number , NumberDensityAndNumber
}
Public Member Functions inherited from MeshMedium
MassType massType () const
double maxX () const
double maxY () const
double maxZ () const
double minX () const
double minY () const
double minZ () const
Public Member Functions inherited from ImportedMedium
Vec bulkVelocity (Position bfr) const override
int dimension () const override
string filename () const
Position generatePosition () const override
bool hasMagneticField () const override
bool hasMetallicity () const override
bool hasParameters () const override
bool hasTemperature () const override
bool hasVariableMix () const override
bool hasVelocity () const override
bool importMagneticField () const
bool importMetallicity () const
bool importTemperature () const
bool importVariableMixParams () const
bool importVelocity () const
Vec magneticField (Position bfr) const override
double mass () const override
double massDensity (Position bfr) const override
double massFraction () const
MaterialMixmaterialMix () const
MaterialMixFamilymaterialMixFamily () const
double maxTemperature () const
double metallicity (Position bfr) const override
const MaterialMixmix () const override
const MaterialMixmix (Position bfr) const override
double number () const override
double numberDensity (Position bfr) const override
int numSites () const override
double opticalDepthX (double lambda) const override
double opticalDepthY (double lambda) const override
double opticalDepthZ (double lambda) const override
void parameters (Position bfr, Array &params) const override
Position sitePosition (int index) const override
const Snapshotsnapshot () const
double temperature (Position bfr) const override
string useColumns () const
Public Member Functions inherited from SimulationItem
template<class T>
T * find (bool setup=true) const
template<class T>
T * interface (int levels=-999999, bool setup=true) const
virtual string itemName () const
void setup ()
string typeAndName () const
Public Member Functions inherited from Item
 Item (const Item &)=delete
virtual ~Item ()
void addChild (Item *child)
const vector< Item * > & children () const
virtual void clearItemListProperty (const PropertyDef *property)
void destroyChild (Item *child)
virtual bool getBoolProperty (const PropertyDef *property) const
virtual vector< double > getDoubleListProperty (const PropertyDef *property) const
virtual double getDoubleProperty (const PropertyDef *property) const
virtual string getEnumProperty (const PropertyDef *property) const
virtual int getIntProperty (const PropertyDef *property) const
virtual vector< Item * > getItemListProperty (const PropertyDef *property) const
virtual ItemgetItemProperty (const PropertyDef *property) const
virtual string getStringProperty (const PropertyDef *property) const
int getUtilityProperty (string name) const
virtual void insertIntoItemListProperty (const PropertyDef *property, int index, Item *item)
Itemoperator= (const Item &)=delete
Itemparent () const
virtual void removeFromItemListProperty (const PropertyDef *property, int index)
virtual void setBoolProperty (const PropertyDef *property, bool value)
virtual void setDoubleListProperty (const PropertyDef *property, vector< double > value)
virtual void setDoubleProperty (const PropertyDef *property, double value)
virtual void setEnumProperty (const PropertyDef *property, string value)
virtual void setIntProperty (const PropertyDef *property, int value)
virtual void setItemProperty (const PropertyDef *property, Item *item)
virtual void setStringProperty (const PropertyDef *property, string value)
void setUtilityProperty (string name, int value)
virtual string type () const
Public Member Functions inherited from SiteListInterface
virtual ~SiteListInterface ()
Public Member Functions inherited from AdaptiveMeshInterface
virtual ~AdaptiveMeshInterface ()

Detailed Description

An AdaptiveMeshMedium instance represents a transfer medium with a spatial density distribution described by an Adaptive Mesh Refinement (AMR) grid partitioning a cuboidal domain. The data is usually extracted from a cosmological simulation snapshot, and it must be provided in a column text file formatted as described below.

Refer to the description of the AdaptiveMeshSnapshot class for information on the structure of an adaptive mesh and on how to represent it in text column file format. Refer to the description of the TextInFile class for information on overall formatting and on how to include header lines specifying the units for each column in the input file. In case the input file has no unit specifications, the default units mentioned below are used instead. The number of columns expected in the input file depends on the options configured by the user for this AdaptiveMeshMedium instance:

\[ \{\, \rho\,(\text{M}_\odot\,\text{pc}^{-3}) \;\;|\;\; M\,(\text{M}_\odot) \;\;|\;\; n\,(\text{cm}^{-3}) \;\;|\;\; N\,(1) \,\} \quad [Z\,(1)] \quad [T\,(\mathrm{K})] \quad [ v_x\,(\mathrm{km/s}) \quad v_y\,(\mathrm{km/s}) \quad v_z\,(\mathrm{km/s}) ] \quad [ B_x\,(\mu\mathrm{G}) \quad B_y\,(\mu\mathrm{G}) \quad B_z\,(\mu\mathrm{G}) ] \quad [ \dots\text{mix family params}\dots ] \]

Depending on the value of the massType option, the first column lists the average mass density \(\rho\), the integrated mass \(M\), the average number density \(n\), or the integrated number density \(N\) for the cell. This quantity is multiplied by the value of the massFraction option.

If the importMetallicity option is enabled, the next column specifies a "metallicity" fraction, which is multiplied with the mass/density column to obtain the actual mass/density of the cell. If the importTemperature option is enabled, the next column specifies a temperature. If this temperature is higher than the value of the maxTemperature option, the mass and density for the cell are set to zero, regardless of the mass or density specified in the first column. If the importTemperature option is disabled, or the maximum temperature value is set to zero, such a cutoff is not applied.

If both the importMetallicity and importTemperature options are enabled, this leads to the following expression for the density of an imported cell (or a simular formula for the other mass quantity types):

\[ \rho_\mathrm{imported} = \begin{cases} f_\mathrm{mass}\,Z\,\rho & \mathrm{if}\; T<T_\mathrm{max} \;\mathrm{or}\; T_\mathrm{max}=0 \\ 0 & \mathrm{otherwise} \end{cases} \]

If the importVelocity option is enabled, the subsequent three columns specify the \(v_x\), \(v_y\), \(v_z\) velocity components of the cell (considered as the average bulk velocity for the mass in the cell).

If the importMagneticField option is enabled, the subsequent three columns specify the \(B_x\), \(B_y\), \(B_z\) magnetic field vector components for the cell.

Finally, if the importVariableMixParams option is enabled, the remaining columns specify the parameters used by the configured material mix family to select a particular material mix for the cell.

When an item of this type is used, the names provided by the conditional value expression "AdaptiveMeshInterface" are inserted into the name sets used for evaluating Boolean expressions.

Constructor & Destructor Documentation

◆ AdaptiveMeshMedium()

AdaptiveMeshMedium::AdaptiveMeshMedium ( )
inlineprotected

Default constructor for concrete Item subclass AdaptiveMeshMedium: "a transfer medium imported from data represented on an adaptive mesh (AMR grid)".

Member Function Documentation

◆ adaptiveMesh()

AdaptiveMeshSnapshot * AdaptiveMeshMedium::adaptiveMesh ( ) const
overrideprotectedvirtual

This function implements the AdaptiveMeshInterface interface. It returns a pointer to the adaptive mesh snapshot maintained by this geometry.

Implements AdaptiveMeshInterface.

◆ createAndOpenSnapshot()

Snapshot * AdaptiveMeshMedium::createAndOpenSnapshot ( )
overrideprotectedvirtual

This function constructs a new AdaptiveMeshSnapshot object, calls its open() function, passes it the domain extent configured by the user, configures it to import a mass or a density column, and finally returns a pointer to the object. Ownership of the Snapshot object is transferred to the caller.

Implements ImportedMedium.


The documentation for this class was generated from the following file: