Différences

Ci-dessous, les différences entre deux révisions de la page.

--- sunfluidh:radiative_heat_transfer_dom_setup_namelist [2016/12/13 12:11] – [ka_min , ka_max] cadet
+++ sunfluidh:radiative_heat_transfer_dom_setup_namelist [2016/12/14 18:42] (Version actuelle) – [Namelist Radiative_Heat_Transfer_DOM] cadet
@@ Ligne 1: / Ligne 1: @@
 ===== Namelist "Radiative_Heat_Transfer_DOM" =====
-This data set is used to define the radiative problem. Otherwise, it can be omitted.
-This module considers the Radiative Transfer Equation (RTE) for an emitting-absorbing non-scattering medium enclosed by difuse boundaries.
-To take into account the gas behavior, it considers both gray-gas assumption as well as real gas behavior through the Spectral-Line-Weighted-Sum-of-Gray-Gases (SLW) model. The final RTE-SLW problem is then discretize with the Discrete Ordinates Method (DOM).
-**The DOM** discretize the $4\pi$ steradians angular integration in a set of $M$ discrete directions represented by their direct cosines and corresponding weights $\vec{q_m} = (\vec{s_m},\omega_m) = (\mu_m,\eta_m,\xi_m,\omega_m)$ for all $m \in [1,M]$.
+__** Not for the release SUNFLUIDH_EDU**__ .\\
+This data set is used to define the radiative problem. Otherwise, it can be omitted.\\
+This module considers the Radiative Transfer Equation (RTE) for an emitting-absorbing non-scattering medium enclosed by diffuse boundaries.
+To take into account the gas behavior, it considers both gray-gas assumption as well as real gas behavior through the Spectral-Line-Weighted-Sum-of-Gray-Gases (SLW) model. \\
+The final RTE-SLW problem is then discretize with the Discrete Ordinates Method (DOM).
+**The DOM** discretize the $4\pi$ steradians integration in a set of $M$ discrete directions represented by their direct cosines and corresponding weights $\vec{q_m} = (\vec{s_m},\omega_m) = (\mu_m,\eta_m,\xi_m,\omega_m)$ for all $m \in [1,M]$.\\
 **The SLW model** will change the spectral integration in a weighted sum of $N_g$ gray-gases represented by their absorption coefficient and corresponding weights $(\kappa_j,a_j)$ for all $j \in [1,N_g]$.
 Thus, the resulting **RTE-SLW-DOM** problem for emitting-absorbing non-scattering medium stands as below :
-$\vec{s}_m \cdot \nabla I_j^m (x_i,\vec{s}_m) = \kappa_j \left[ a_j I_b({T}(x_i)) - I_j^m(x_i,\vec{s}_m) \right]$     $; \quad \forall (m,j) \in [M,N_g]$
+\begin{equation}
+\vec{s}_m \cdot \nabla I_j^m (x_i,\vec{s}_m) = \kappa_j \left[ a_j I_b({T}(x_i)) - I_j^m(x_i,\vec{s}_m) \right]; \quad \forall (m,j) \in [M,N_g]
+\end{equation}
-where $I_j^m$ is the radiative intensity for the virtual gray-gas $j$ in direction $m$ and $I_b$ is the blackbody radiative intensity and $(\kappa_j,a_j)$.
+where $I_j^m$ is the radiative intensity for the virtual gray-gas $j$ in direction $m$ and $I_b$ is the blackbody radiative intensity.
-The **dimensional** radiative source term $S_r$ and the boundary net radiative heat flux $q_r^{net}$ are defined as :
+The **dimensional** radiative source term $S_r$ and boundary net radiative heat flux $q_r^{net}$ are defined as :
-$S_r(x_i,{T}) = - \sum_{j=0}^{N_g} \kappa_j \left[ \sum_{m=1}^{M} \omega_m I_j^m (x_i,\vec{s}_m)  - 4 a_j \sigma_B ({T}(x_i))^4 \right] $
-${q}_r^{net}(x_i^{wall}) = \varepsilon_{wall} \left[ \sigma_B ({T}(x_i^{wall}))^4 - \sum_{j=0}^{N_g} \sum_{m:\vec{s}_m \cdot \vec{n} > 0} \omega_m |\vec{s}_m\cdot \vec{n}| I_j^m (x_i^{wall},\vec{s}_m) \right]$
+\begin{eqnarray}
+S_r(x_i,{T}) & = & - \sum_{j=0}^{N_g} \kappa_j \left[ \sum_{m=1}^{M} \omega_m I_j^m (x_i,\vec{s}_m)  - 4 a_j \sigma_B ({T}(x_i))^4 \right] \\
+{q}_r^{net}(x_i^{wall}) & = & \varepsilon_{wall} \left[ \sigma_B ({T}(x_i^{wall}))^4 - \sum_{j=0}^{N_g} \sum_{m:\vec{s}_m \cdot \vec{n} > 0} \omega_m |\vec{s}_m\cdot \vec{n}| I_j^m (x_i^{wall},\vec{s}_m) \right]
+\end{eqnarray}
 where $\sigma_b$ is the Stefan-Boltzmann constant, $\varepsilon$ is the boundary emissivity and $\vec{n}$ is the normal to the wall pointing out of the domain.
-<note important>This radiative solver implementation considers only **3D cartesian** problems under **MPI cartesian domain decomposition** approach.</note>
+<note warning>This radiative solver implementation considers only **cartesian** problems and **does not support immersed bodies**.</note>
 ===== Full data set of the namelist =====
@@ Ligne 38: / Ligne 43: @@
 -----
-===== Definition of the data set for the RTE problem =====
+===== Definition of the data set for the DOM-RTE problem =====
 -----
 ==== activateRadiation ====
@@ Ligne 64: / Ligne 69: @@
 ==== WallRadcoeff ====
    * Type : Real value
-   * Prescaler on the net radiative heat flux at walls.
+   * Prescaler on the net radiative heat flux $q_r^{net}$ at walls.
    * **For debugging only**.
    * Default value = 1.0
 ==== VolRadcoeff ====
    * Type : Real value
-   * Prescaler on the radiative source term.
+   * Prescaler on the radiative source term $S_r$.
    * **For debugging only**.
    * Default value = 1.0
@@ Ligne 88: / Ligne 93: @@
 ==== Tref ====
    * Type : Real value.
-   * Reference temperature $T_{ref}$.
+   * Reference temperature $T_{ref}$ in [$K$].
    * Default value = Fluid_Properties%Reference_Temperature
 ==== Href ====
    * Type : Real value.
-   * Reference Length $H_{ref}$.
+   * Reference Length $H_{ref}$ in [$m$].
    * Default value = Nondimensionalization%Reference_Length
 ==== WallEmissivity ====
@@ Ligne 98: / Ligne 103: @@
    * Boundaries emissivities $\varepsilon$ sorted as (x-,x+,y-,y+,z-,z+).
    * Default values = 0.0 0.0 0.0 0.0 0.0 0.0
------
 ===== Definition of the data set for the SLW model =====
 -----
@@ Ligne 116: / Ligne 121: @@
      * NbGas $\ge$ 2  : SLW model is employed
    * Default value = 1
+<note> Setting **NbGas $=$ 1 and ka_min = 0** is equivalent to **wall-to-wall radiation** du to the presence of transparent medium </note>
 ==== ka_min , ka_max ====
    * Type : Real values
-   * These options set the lower and upper bounds of absorbing coefficient for the SLW model.
+   * These options set the lower and upper bounds of dimensional absorbing coefficient [$m^{-1}$] for the SLW model.
      * if NbGas $=$ 1 : $\kappa$ = //ka_min//, **//ka_max// is useless**
      * if NbGas $\ge$ 2 : //ka_min// < $\kappa_j$ < //ka_max//  for all $j \in [1,N_g]$
@@ Ligne 125: / Ligne 132: @@
    * Type : Character string with a maximum size of 3
    * Name of the absorbing species when NbGas $\ge$ 2 (SLW model).
+     * **if "NbGas" $=$ 1 : useless**
    * Available values :
      * "H2O" : $air-H_2O$ mixture
@@ Ligne 132: / Ligne 140: @@
    * Type : Real value
    * This option set the reference molar fraction $x_{ref}$ of the absorbing species for the SLW model.
-     * if "NbGas" $=$ 1 : useless
+     * **if "NbGas" $=$ 1 : useless**
    * Default value = 0.07
 ==== xaUniform ====
    * Type : Real value
    * As long as the SLW model is not coupled with species equations, this option set a uniform molar fraction $x_{a}$ of the absorbing species in the overall domain.
-     * if NbGas $=$ 1 : useless
+     * **if NbGas $=$ 1 : useless**
    * Default value = 0.07
 ==== Pref ====
    * Type : Real value.
-   * Reference pressure $P_{ref}$.
+   * Reference pressure $P_{ref}$ in [$Pa$].
+     * **if "NbGas" $=$ 1 : useless**
    * Default value = obtained from Fluid_Properties quantities