The SKIRT project
advanced radiative transfer for astrophysics
ModifiedLogNormalGrainSizeDistribution Class Reference

#include <ModifiedLogNormalGrainSizeDistribution.hpp>

Inheritance diagram for ModifiedLogNormalGrainSizeDistribution:

Public Member Functions

double dnda (double a) const override
double firstMixingParameter () const
double secondMixingParameter () const
Public Member Functions inherited from LogNormalGrainSizeDistribution
 LogNormalGrainSizeDistribution (SimulationItem *parent, double minSize, double maxSize, double centroid, double width)
double centroid () const
double width () const
Public Member Functions inherited from RangeGrainSizeDistribution
double amax () const override
double amin () const override
double maxSize () const
double minSize () 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

Protected Member Functions

 ModifiedLogNormalGrainSizeDistribution ()
Protected Member Functions inherited from LogNormalGrainSizeDistribution
 LogNormalGrainSizeDistribution ()
Protected Member Functions inherited from RangeGrainSizeDistribution
 RangeGrainSizeDistribution ()
 RangeGrainSizeDistribution (double minSize, double maxSize)
void setupSelfBefore () override
Protected Member Functions inherited from GrainSizeDistribution
 GrainSizeDistribution ()
Protected Member Functions inherited from SimulationItem
 SimulationItem ()
virtual bool offersInterface (const std::type_info &interfaceTypeInfo) const
virtual void setupSelfAfter ()
Protected Member Functions inherited from Item
 Item ()

Private Types

using BaseType
using ItemType

Private Attributes

double _firstMixingParameter
double _secondMixingParameter
const double & _y0
const double & _y1

Friends

class ItemRegistry

Detailed Description

ModifiedLogNormalGrainSizeDistribution is a GrainSizeDistribution subclass that represents a modified log-normal dust grain size distribution of the form

\[\frac{\text{d}n_\text{D}}{\text{d}a} \propto \frac{1}{a} \,\exp\left[ - \frac{(\ln(a/a_0))^2}{2\sigma^2} \right] \, M(a) \qquad \text{for}\quad a_\text{min} \leq a \leq a_\text{max} \]

with a mixing term

\[ M(a) = y_0 + (y_1-y_0)\frac{\ln(a/a_\text{min})}{\ln(a_\text{max}/a_\text{min})}. \]

The size range of the distribution can be configured in the RangeGrainSizeDistribution base class. The centroid \(a_0\) and the width \(\sigma\) can be configured in the LogNormalGrainSizeDistribution base class. The remaining two parameters \(y_0\) and \(y_1\) can be configured as attributes in this class. The function is scaled arbitrarily.

The functional form for the grain size distribution implemented by this class is inspired by the DustEM code, which is described in Compiègne et al. 2011 (AA, 525, A103) and can be downloaded from http://www.ias.u-psud.fr/DUSTEM/.

Constructor & Destructor Documentation

◆ ModifiedLogNormalGrainSizeDistribution()

ModifiedLogNormalGrainSizeDistribution::ModifiedLogNormalGrainSizeDistribution ( )
inlineprotected

Default constructor for concrete Item subclass ModifiedLogNormalGrainSizeDistribution: "a modified log-normal dust grain size distribution".

Member Function Documentation

◆ dnda()

double ModifiedLogNormalGrainSizeDistribution::dnda ( double a) const
overridevirtual

This function returns the value of \(\frac{\text{d}n_\text{D}}{\text{d}a}\) as described in the header for this class (with an arbitrary proportionality factor of one).

Reimplemented from LogNormalGrainSizeDistribution.

◆ firstMixingParameter()

double ModifiedLogNormalGrainSizeDistribution::firstMixingParameter ( ) const
inline

This function returns the value of the discoverable double property firstMixingParameter: "the first mixing parameter y0".

The minimum value for this property is "[0".

The maximum value for this property is "1]".

The default value for this property is given by the conditional value expression "1".

◆ secondMixingParameter()

double ModifiedLogNormalGrainSizeDistribution::secondMixingParameter ( ) const
inline

This function returns the value of the discoverable double property secondMixingParameter: "the second mixing parameter y1".

The minimum value for this property is "[0".

The maximum value for this property is "1]".

The default value for this property is given by the conditional value expression "1".


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