astromugs.modeling.Disk#

class astromugs.modeling.Disk(params: DiskParams, dust: object | None = None)[source]#

Static flared disk model for radiative transfer and chemistry simulations.

Provides methods to compute gas and dust physical quantities (density, temperature, scale height, surface density, extinction) on various coordinate grids.

Parameters:

  • params (DiskParams) – Dataclass holding all disk model parameters.

  • dust (object, optional) – Dust model object providing grain size distributions and opacities.

__init__(params: DiskParams, dust: object | None = None)[source]#

Initialize the Disk model.

Parameters:

  • params (DiskParams) – Dataclass holding all disk model parameters. Its fields are also copied as instance attributes for convenience.

  • dust (object, optional) – Dust model object providing grain size distributions and opacities.

Methods

__init__(params[, dust])

Initialize the Disk model.

av_z(lam, nd, r, z)

Compute the vertical visual extinction profile.

avz_ism(lam[, A, q_c])

Compute visual extinction using the ISM H-column-density-to-Av conversion.

density(x1, x2[, x3])

Compute the gas mass density assuming hydrostatic equilibrium.

density_d(x1, x2[, x3])

Compute the dust mass density for each grain-size bin.

numberdensity(x1, x2[, x3])

Compute the gas number density assuming hydrostatic equilibrium.

numberdensity_d(x1, x2[, x3])

Compute the dust number density for each grain-size bin.

numberdensity_d_single(x1, x2[, x3])

Compute the dust number density for a single representative grain.

omega2(r)

Compute the square of the Keplerian angular velocity.

q_ext(lam[, A, q_c])

Compute the extinction efficiency for each grain size and wavelength.

scaleheight(r)

Compute the gas scale height.

scaleheight_d(r)

Compute the dust scale height for each grain size.

surfacedensity(r)

Compute the gas surface density.

surfacedensity_d(r)

Compute the dust surface density per grain-size bin.

temp_altitude(r, z)

Compute the vertical temperature profile following Williams & Best (2014).

temp_atmos(r)

Compute the atmospheric temperature radial profile.

temp_mid(r)

Compute the midplane temperature radial profile.

verticaldensity_gauss(r)

Compute the gas number density assuming a Gaussian vertical profile.

viscous_accretion_heating(x1, x2[, x3])

Compute the viscous accretion heating rate.