The SKIRT project
advanced radiative transfer for astrophysics
Protected Member Functions | Private Types | Private Attributes | Friends | List of all members
AdaptiveMeshMedium Class Reference

#include <AdaptiveMeshMedium.hpp>

Inheritance diagram for AdaptiveMeshMedium:
Inheritance graph
[legend]

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 ()
 
virtual SnapshotcreateAndOpenSnapshot ()=0
 
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
 
virtual void setupSelfAfter ()
 
virtual void setupSelfBefore ()
 
- 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 = MeshMedium
 
using ItemType = AdaptiveMeshMedium
 

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
 
virtual Vec bulkVelocity (Position bfr) const =0
 
virtual int dimension () const =0
 
virtual Position generatePosition () const =0
 
virtual bool hasMagneticField () const =0
 
virtual bool hasMetallicity () const =0
 
virtual bool hasParameters () const =0
 
virtual bool hasTemperature () const =0
 
virtual bool hasVariableMix () const =0
 
virtual bool hasVelocity () const =0
 
virtual Vec magneticField (Position bfr) const =0
 
virtual double mass () const =0
 
virtual double massDensity (Position bfr) const =0
 
virtual double metallicity (Position bfr) const =0
 
virtual const MaterialMixmix () const =0
 
virtual const MaterialMixmix (Position bfr) const =0
 
virtual double number () const =0
 
virtual double numberDensity (Position bfr) const =0
 
virtual double opticalDepthX (double lambda) const =0
 
virtual double opticalDepthY (double lambda) const =0
 
virtual double opticalDepthZ (double lambda) const =0
 
virtual void parameters (Position bfr, Array &params) const =0
 
virtual double temperature (Position bfr) const =0
 
- 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 ()
 
virtual int numSites () const =0
 
virtual Position sitePosition (int index) const =0
 
- Public Member Functions inherited from AdaptiveMeshInterface
virtual ~AdaptiveMeshInterface ()
 
virtual AdaptiveMeshSnapshotadaptiveMesh () const =0
 

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: