Update Slope-Limiters-and-Shock-Resolution.md

_data/docs_v7.yml

@@ -30,6 +30,7 @@
   docs_v7:
   - Theory
   - Streamwise-Periodicity
+  - Thermochemical-Nonequilibrium
 
 - title: Users Guide
   docs_v7:

_docs_v7/Custom-Output.md

@@ -218,6 +218,7 @@ The available types are:
 - `Function`: Introduces a new scalar output that is a function of other scalar outputs, it cannot reference fields (e.g. velocity).
 - `AreaAvg` and `AreaInt`: Computes an area average or integral of a field (the expression) over the list of markers.
 - `MassFlowAvg` and `MassFlowInt`: Computes a mass flow average or integral.
+- `Probe`: Evaluates the expression using the values of the mesh point closest to the coordinates specified inside "[]", [x, y] or [x, y, z] (2 or 3D).
 
 **Note:** Each custom output can only use one type, e.g. it is not possible to write `p_drop : AreaAvg{PRESSURE}[inlet] - AreaAvg{PRESSURE}[outlet]`. This would need to be separated into two `AreaAvg` outputs and one `Function` to compute their difference.
 
@@ -226,7 +227,8 @@ The available types are:
 CUSTOM_OUTPUTS= 'velocity : Macro{sqrt(pow(VELOCITY_X, 2) + pow(VELOCITY_Y, 2) + pow(VELOCITY_Z, 2))};\
                  avg_vel : AreaAvg{$velocity}[z_minus, z_plus];\
                  var_vel : AreaAvg{pow($velocity - avg_vel, 2)}[z_minus, z_plus];\
-                 dev_vel : Function{sqrt(var_vel) / avg_vel}'
+                 dev_vel : Function{sqrt(var_vel) / avg_vel};\
+                 probe1 : Probe{$velocity}[0.005, 0.005, 0.05]'
 ```
 
 To obtain the list of solver variables that can be used, write an invalid expression (e.g. 'x : AreaAvg{INVALID}[]') and run SU2.

_docs_v7/Physical-Definition.md

@@ -19,6 +19,7 @@ SU2 offers different ways of setting and computing this definition. This documen
 - [Free-Stream Definition (Thermochemical Nonequilibrium)](#free-stream-definition-thermochemical-nonequilibrium)
   - [Free-Stream Temperatures](#free-stream-temperatures)
   - [Chemical Composition and Mass Fractions](#chemical-composition-and-mass-fractions)
+- [Transport Coefficient Model (Thermochemical Nonequilibrium)](#transport-coefficient-model-thermochemical-nonequilibrium)
 - [Flow Condition (Incompressible)](#flow-condition-incompressible)
   - [Thermodynamic and Gauge Pressure](#thermodynamic-and-gauge-pressure)
   - [Initial State and Non-Dimensionalization](#initial-state-and-non-dimensionalization)
@@ -128,6 +129,18 @@ A chemistry model, consisting of a set of flow species, thermochemical propertie
 
 Free-stream mass fractions must also be specified in list using the option `GAS_COMPOSITION`. The mass fractions are specified as decimal values in the order of the species in the gas model. For example, an AIR-5 mixture of 77% oxygen and 23% nitrogen would be expressed as (0.77, 0.23, 0.00, 0.00, 0.00).
 
+## Transport Coefficient Model (Thermochemical Nonequilibrium) ##
+
+| Solver | Version |
+| --- | --- |
+| `NEMO_NAVIER_STOKES` | 7.0.0 |
+
+A transport coefficient model must be specified for viscous simulations with the NEMO solver, using the `TRANSPORT_COEFF_MODEL` config option. Available options when using the SU2TCLib thermochemical library are the Wilkes-Blottner-Eucken, Gupta-Yos, and Sutherland viscosity models, specified by `WILKE`, `GUPTA-YOS`, and `SUTHERLAND`, respectively. The default option for transport coefficient model is Wilkes-Blottner-Eucken. 
+
+It should be noted the Sutherland model is only used to evaluate viscosity, and the Wilkes-Blottner-Eucken model is used to evaluate diffusion coefficient and thermal conductivity.
+
+---
+
 ## Flow Condition (Incompressible) ##
 
 | Solver | Version |

_docs_v7/Slope-Limiters-and-Shock-Resolution.md

@@ -96,6 +96,9 @@ The `SLOPE_LIMITER_` options above may each be changed to use different limiters
 | `BARTH_JESPERSEN`       | Barth-Jespersen                                 | This limiter is a smooth version of the commonly seen Barth-Jespersen limiter seen in the literature  |
 | `VENKATAKRISHNAN`       | Venkatakrishnan                                 |  |
 | `VENKATAKRISHNAN_WANG`  | Venkatakrishnan-Wang                            |  |
+| `NISHIKAWA_R3`          | Nishikawa-R3                                    |  |
+| `NISHIKAWA_R4`          | Nishikawa-R4                                    |  |
+| `NISHIKAWA_R5`          | Nishikawa-R5                                    |  |
 | `SHARP_EDGES`           | Venkatakrishnan with sharp-edge modification    | This limiter should not be used for flow solvers |
 | `WALL_DISTANCE`         | Venkatakrishnan with wall distance modification | This limiter should not be used for flow solvers |
 | `VAN_ALBADA_EDGE`       | Van Albada (edge formulation)                   | This limiter is only implemented for flow solvers and does not output limiter values when using the VOLUME_OUTPUT option |
@@ -106,7 +109,7 @@ The default limiter is `VENKATAKRISHNAN`.
 
 The `VENKAT_LIMITER_COEFF` parameter is generally a small constant, defaulting to $$0.05$$, but its specific definition depends on the limiter being used.
 
-For the `VENKATAKRISHNAN`, `SHARP_EDGES`, and `WALL_DISTANCE` limiters, the `VENKAT_LIMITER_COEFF` parameter refers to $$K$$ in $$\epsilon^2=\left(K\bar{\Delta} \right)^3$$, where $$\bar{\Delta}$$ is an average grid size (this is hardcoded as 1m and thus all tuning is via $$K$$).
+For the `VENKATAKRISHNAN`, `SHARP_EDGES`, and `WALL_DISTANCE` limiters, the `VENKAT_LIMITER_COEFF` parameter refers to $$K$$ in $$\epsilon^2=\left(K\bar{\Delta} \right)^3$$, where $$\bar{\Delta}$$ is an average grid size (this is hardcoded as 1m and thus all tuning is via $$K$$). For NISHIKAWA_Rp limiters, $$\epsilon^p=\left(K\bar{\Delta} \right)^{p+1}$$ (p = 3, 4 or 5).
 The $$K$$ parameter defines a threshold, below which oscillations are not damped by the limiter, as described by [Venkatakrishnan](https://doi.org/10.1006/jcph.1995.1084).
 Thus, a large value will approach the case of using no limiter with undamped oscillations, while too small of a value will slow the convergence and add extra diffusion.
 The SU2 implementation of the `BARTH_JESPERSEN` limiter actually uses `VENKATAKRISHNAN` with $$K=0$$.
@@ -115,10 +118,10 @@ The SU2 implementation of the `BARTH_JESPERSEN` limiter actually uses `VENKATAKR
 When using the `VENKATAKRISHNAN_WANG` limiter, `VENKAT_LIMITER_COEFF` is instead $$\varepsilon '$$ in $$\varepsilon = \varepsilon ' (q_{max} - q_{min})$$, where $$q_{min}$$ and $$q_{max}$$ are the respective *global* minimum and maximum of the field variable being limited.
 This global operation incurs extra time costs due to communication between MPI ranks.
 The original work by [Wang](https://doi.org/10.2514/6.1996-2091) suggests using `VENKAT_LIMITER_COEFF` in the range of $$[0.01, 0.20]$$, where again larger values approach the case of using no limiter.
-**Note:** unlike the aforementioned `VENKATAKRISHNAN` limiter, the `VENKATAKRISHNAN_WANG` limiter does not depend directly on the mesh size and can thus be used without non-dimensionalization. If the `VENKATAKRISHNAN` limiter is used outside of non-dimensional mode, the fields with larger values (pressure and temperature) will generally be limited more aggressively than velocity.
+**Note:** unlike the aforementioned `VENKATAKRISHNAN` limiter and NISHIKAWA_Rp limiter, the `VENKATAKRISHNAN_WANG` limiter does not depend directly on the mesh size and can thus be used without non-dimensionalization. If the `VENKATAKRISHNAN` limiter is used outside of non-dimensional mode, the fields with larger values (pressure and temperature) will generally be limited more aggressively than velocity.
 
 
-The `NONE`, `BARTH_JESPERSEN`, `VENKATAKRISHNAN`, and `VENKATAKRISHNAN_WANG` limiter options all have no **geometric modifier**.
+The `NONE`, `BARTH_JESPERSEN`, `VENKATAKRISHNAN`, `VENKATAKRISHNAN_WANG`, and NISHIKAWA_Rp limiter options all have no **geometric modifier**.
 A geometric modifier increases limiting near walls or sharp edges. This is done by multiplying the limiter value by a **geometric factor**. 
 
 For both the `SHARP_EDGES` and `WALL_DISTANCE` limiters, the influence of the geometric modifier is controlled with `ADJ_SHARP_LIMITER_COEFF` which defaults to 3.0.

_docs_v7/Style-Guide.md

@@ -25,6 +25,8 @@ Patch releases cannot break backward compatibility.
 SU2 is written for C++11, the formatting rules are defined in a `clang-format` file located in the root of the repository.
 **New files must follow the formatting rules exactly.**
 
+SU2 uses pre-commit to enforce a consistent formatting. To use, [install pre-commit](https://pre-commit.com/#install) and run `pre-commit install` at the root of the project. You can now force the formatting on all files with `pre-commit run -a`. This will also run all pre-commit hooks before each commit, preventing dirty commits in the repository.
+
 ### Files, functions, and variables
 
 Basic recommendations for creating files, functions, and variables:

_docs_v7/Thermochemical-Nonequilibrium.md

@@ -0,0 +1,400 @@
+---
+title: Thermochemical Nonequilibrium
+permalink: /docs_v7/Thermochemical-Nonequilibrium/
+---
+
+This page contains a summary of the physical models implemented in the NEMO solvers in SU2 designed ot simulate hypersonic flows in thermochemical nonequilibrium. This includes detials on thermodynamic and chemistry models, as well as transport properties and boundary conditions. 
+
+---
+
+- [Thermodynamic Model](#thermodynamic-model)
+- [Finite Rate Chemistry](#finite-rate-chemistry)
+- [Vibrational Relaxation](#vibrational-relaxation)
+- [Viscous Phenomena and Transport Coefficients](#viscous-phenomena-and-transport-coefficients)
+    - [Wilkes-Blottner-Eucken](#wilkes-blottner-eucken)
+    - [Gupta-Yos](#gupta-yos)
+    - [Sutherland Viscosity Model](#sutherland-viscosity-model)
+- [Slip Flow](#slip-flow)
+- [Gas-surface Interaction](#gas-surface-interaction)
+
+---
+
+# Thermodynamic Model #
+
+| Solver | Version | 
+| --- | --- |
+| `NEMO_EULER`, `NEMO_NAVIER_STOKES` | 7.0.0 |
+
+A rigid-rotor harmonic oscillator (RRHO) two-temperature model is used to model the thermodynamic state of continuum hypersonic flows. Through the independence of the energy levels, the total energy and vibrational--electronic energy per unit volume can be expressed as
+$$   \rho e = \sum_s \rho_s \left(e_s^{tr} + e_s^{rot} + e_s^{vib} + e_s^{el} + e^{\circ}_s + \frac{1}{2} \bar{v}^{\top} \bar{v}\right),
+$$
+and
+$$
+        \rho e^{ve} = \sum_s \rho_{s} \left(e_s^{vib} + e_s^{el}\right).
+$$
+
+Considering a general gas mixture consisting of polyatomic, monatomic, and free electron species, expressions for the energy stored in the translational, rotational, vibrational, and electronic modes are given as
+$$
+    e^{tr}_s =\begin{cases}
+    \frac{3}{2} \frac{R}{M_s} T & \text{for monatomic and polyatomic species,}\\
+    0 & \text{for electrons,}
+    \end{cases}
+$$
+$$
+    e^{rot}_s =\begin{cases}
+    \frac{\xi  }{2} \frac{R}{M_s} T & \text{for polyatomic species,}\\
+    0 & \text{for monatomic species and electrons,}
+    \end{cases}
+$$
+\
+where $$\xi$$ is an integer specifying the number of axes of rotation,
+$$
+    e^{vib}_s =\begin{cases}
+    \frac{R}{M_s} \frac{\theta^{vib}_s}{exp\left( \theta^{vib}_s / T^{ve}\right) - 1} & \text{for polyatomic species,}\\
+    0 & \text{for monatomic species and electrons,}
+    \end{cases}
+$$
+\
+where $$\theta^{vib}_s$$ is the characteristic vibrational temperature of the species, and 
+
+
+$$
+    e^{el}_s =\begin{cases}
+    \frac{R}{M_s}\frac{\sum_{i=1}^{\infty} g_{i,s}{\theta^{el}_{i,s} \exp(-\theta^{el}_{i,s}/T_{ve})}}{\sum_{i=0}^{\infty} g_{i,s} exp(-\theta^{el}_{i,s}/T_{ve})} & \text{for polyatomic and monatomic species,}\\
+    \frac{3}{2} \frac{R}{M_s} T^{ve} & \text{for electrons,}
+    \end{cases}
+$$
+
+where $$\theta^{el}_s$$ is the characteristic electronic temperature of the species and $$g_i$$ is the degeneracy of the $$i^{th}$$ state.
+
+---
+
+# Finite Rate Chemistry #
+
+| Solver | Version | 
+| --- | --- |
+| `NEMO_EULER`, `NEMO_NAVIER_STOKES` | 7.0.0 |
+
+The source terms in the species conservation equations are the  volumetric mass production rates which are governed by the forward and backward reaction rates, $$R^f$$ and $$R^b$$, for a given reaction $$r$$, and can be expressed as
+$$
+    \dot{w}_s = M_s \sum_r (\beta_{s,r} - \alpha_{s,r})(R_{r}^{f} - R_{r}^{b}). 
+$$
+
+From kinetic theory, the forward and backward reaction rates are dependent on the molar concentrations of the reactants and products, as well as the forward and backward reaction rate coefficients, $$k^f$$ and $$k^b$$, respectively, and can be expressed as
+$$
+    R_{r}^f = k_{r}^f \prod_s (\frac{\rho_s}{M_s})^{\alpha_{s,r}},
+$$
+and
+$$
+    R_{r}^b = k_{r}^b \prod_s (\frac{\rho_s}{M_s})^{\beta_{s,r}}.
+$$
+
+For an Arrhenius reaction, the forward reaction rate coefficient can be computed as
+$$
+    k_{r}^f = C_r(T_r)^{\eta_r} exp\left(- \frac{\epsilon_r}{k_B T_r}\right),
+$$
+where $$C_r$$ is the pre-factor, $$T_r$$ is the rate-controlling temperature for the reaction, $$\eta_r$$ is an empirical exponent, and $$\epsilon_r$$ is the activation energy per molecule.
+
+The rate-controlling temperature of the reaction is calculated as a geometric average of the translation-rotational and vibrational-electronic temperatures,
+$$
+    T_r = (T)^{a_r}(T^{ve})^{b_r}.
+$$
+
+The value of he equilibrium constant $$K_{eq}$$ is expressed as 
+
+$$
+    K_{eq} = \exp( A_0 \left(\frac{T^c}{10,000}\right) + A_1 + A_2 \log \left( \frac{10,000}{T^c} \right) + A_3 \left( \frac{10,000}{T^c} \right) + A_4 \left( \frac{10,000}{T^c} \right)^2  ),
+$$
+
+where $$T^c$$ is a controlling temperature and $$A_0 - A_4$$ are constants dependent on the reaction. These reaction constants, the rate constrolling temperature and Arrhenius parameters are stored within the fluid model class in SU2 NEMO.
+
+---
+
+# Vibrational Relaxation #
+
+| Solver | Version | 
+| --- | --- |
+| `NEMO_EULER`, `NEMO_NAVIER_STOKES` | 7.0.0 |
+
+Vibrational relaxation is computed using a standard Landau-Teller relaxation time with a Park high-temperature correction
+$$
+     \dot{\Theta}^{tr:ve} = \sum _s \rho_s \frac{de^{ve}_{s}}{dt} = \sum _s \rho_s \frac{e^{ve*}_{s} - e^{ve}_{s}}{\tau_s},
+$$
+where $$\tau_s$$ is computed using a combination of the Landau-Teller relaxation time, $$\langle \tau_s \rangle_{L-T}$$, and a limiting relaxation time from Park, $$\tau_{ps}$$ using
+$$
+    \tau_s = \langle \tau_s \rangle_{L-T} + \tau_{ps},
+$$
+and
+$$
+    \langle \tau_s \rangle_{L-T} = \frac{\sum_r X_r}{\sum_r X_r/\tau_{sr}}.
+$$
+The interspecies relaxation times are taken from experimental data from Millikan and White, expressed as
+$$
+    \tau_{sr} = \frac{1}{P}exp\left[A_sr\left(T^{-1/3} - 0.015\mu_{sr}^{1/4}\right) - 18.42\right].
+$$
+A limiting relaxation time, $$\tau_{ps}$$, is used to correct for under-prediction of the Millikan--White model at high temperatures. $$\tau_{ps}$$ is defined as
+$$
+    \tau_{ps} = \frac{1}{\sigma_s c_s n},
+$$
+
+where $$\sigma_s$$ is the effective collision~cross-section.
+
+---
+
+# Viscous Phenomena and Transport Coefficients #
+
+| Solver | Version | 
+| --- | --- |
+| `NEMO_NAVIER_STOKES` | 7.0.0 |
+
+
+Mass, momentum, and  energy transport in fluids are all governed by molecular collisions, and  expressions for these transport properties can be derived from the kinetic theory. The  mass diffusion fluxes, $$\mathbf{J}_s$$, are computed using Fick's Law of Diffusion:
+$$
+    \mathbf{J}_s = - \rho D_s \nabla(Y_s) + Y_s \sum_k \rho D_k \nabla(Y_k)
+$$
+
+where $$c_s$$ is the species mass fraction and $$D_s$$ is the species multi-component diffusion coefficient. The  values of $$D_s$$ are computed as a weighted sum of binary diffusion coefficients between all species in the mixture. These are obtained by solving the Stefan--Maxwell equations under the Ramshaw approximations. The  viscous stress tensor is written as
+$$
+   \boldsymbol{\sigma} = \mu \left( \nabla \mathbf{u} + \nabla {\mathbf{u}}^\mathsf{T} - \frac{2}{3} \mathbf{I} (\nabla \cdot \mathbf{u}) \right),
+$$
+where $$\mu$$ is the mixture viscosity coefficient. The  conduction heat flux for each thermal energy mode, $$\mathbf{q}^{k}$$, is modeled by Fourier’s Law of heat conduction:
+$$
+\mathbf{q}^{k} = \kappa^{k}  \nabla(T^k),
+$$
+
+where $$\kappa^{k}$$ is the thermal conductivity associated with energy mode $$k$$.
+
+$$D_s$$, $$\mu$$, and $$\kappa$$ can be evaluated using the models discussed below by selecting the appropriate options in the configuration file.
+
+
+## Wilkes-Blottner-Eucken ##
+
+The mixture dynamic viscosity and thermal conductivity are computed using Wilke's semi-empirical mixing rule as 
+
+$$
+\mu = \sum_s \frac{X_s \mu_s}{\phi_s},
+$$
+
+and 
+
+$$
+\kappa = \sum_s \frac{X_s \kappa_s}{\phi_s},
+$$
+
+where $$X_s$$ is the mole fraction of species $$s$$. The species dynamic viscosity is computed using Blottner's three paramter curve fit for high temperature air,
+
+$$
+\mu_s = 0.1 \exp [(A_s\log(T) + B_s)\log(T) + C_s].
+$$
+
+The species thermal conductivities are computed according to Eucken's formula as
+
+$$
+\kappa^{tr}_s = \mu_s \left( \frac{5}{2} C_{v_s}^{trans} + C_{v_s}^{rot} \right),
+$$
+
+$$
+\kappa^{ve}_s = \mu_s C^{ve}_{v_s}.
+$$
+
+And the term $$\phi_s$$ is given by
+$$
+\phi_s = \sum_r X_r \left[ 1 + \sqrt{\frac{\mu_r}{\mu_s}}\left( \frac{M_r}{M_s} \right)^{1/4} \right]^{2} \left[ \sqrt{8 \left(1 + \frac{M_s}{M_r} \right)} \right]^{-1}.
+$$
+
+The effective species diffusion coefficeint is copmuted as a weighted sum of the species binary diffusion coefficients
+
+$$
+\frac{(1 - X_i)}{D_i} = \sum_{i \neq j} \frac{X_j}{D_{ij}},
+$$
+
+where the binary diffusion coefficients are computed as
+
+$$
+\rho D_{ij} = 1.1613 \times 10^{-25} \frac{M \sqrt{T \left( \frac{1}{M_i} + \frac{1}{M_j} \right) }}{\Omega_{ij}^{(1,1)}}.
+$$
+
+Collision integrals are computed using a four parameter curve fit for neutral-neutral, neutral-ion, and electron-ion collisions
+
+$$
+\pi \Omega_{ij}^{(n,n)} = D T^{A(\log(T))^2 + B \log(T) + C},
+$$
+
+where A-D are constants. Ion-ion, electron-ion, and electron-electron collisions modeled using a shielded Coulomb potential as
+
+$$
+\pi \Omega_{ij}^{(n,n)} = 5.0 \times 10^{15} \pi (\lambda_D / T)^2 \log \{D_n T^{*} \left[ 1 - C_n \exp\left( -c_n T^{*} \right) \right] + 1 \}
+$$
+
+where
+
+$$
+T^{*} = \frac{\lambda_D}{e^2_{CGS} / (k_{B,CGS} T) }
+$$
+
+and the Debye length $$\lambda_D$$ is defined as
+
+$$
+\lambda_D = \sqrt{\frac{k_{B,CGS} T}{4 \pi n_{e,CGS} e^2_{CGS}}}.
+$$
+
+The Wilkes-Blottner-Eucken model is generally efective up to temperatures of 10,000 K. Above these temperatures it is recommended to use the Gupta-Yos model.
+
+## Gupta-Yos ##
+
+Aother model develped by Gupta focuses on the transport properties of weakly ionized flows, and is generally more  accurate than the Wilkes-Blottner-Eucken model at temperatures above 10,000 K.
+
+The forumalae for the transport coefficients are dependent on the collision terms
+
+$$
+\Delta_{s,r}^{(1)}(T) = \frac{8}{3} \left[ \frac{2M_s M_r}{\pi R T (M_s + M_r)}  \right]^{1/2} \pi {\Omega_{s,r}^{(1,1)}}
+$$
+
+and 
+$$
+\Delta_{s,r}^{(2)}(T) = \frac{16}{5} \left[ \frac{2M_s M_r}{\pi R T (M_s + M_r)}  \right]^{1/2} \pi {\Omega_{s,r}^{(2,2)}},
+$$
+
+where the collision cross-sections are computed as described in the Wilkes-Blottner-Eucken section.
+
+The mixutre viscoisty is computed as
+
+$$
+\mu = \sum_{s \neq e} \frac{m_s \gamma_s}{\sum_{r \neq e} \gamma_r \Delta_{s,r}^{(2)}(T_{tr}) + \gamma_r \Delta_{e,r}^{(2)}(T_{ve})} + \frac{m_e \gamma_e}{\sum_r \gamma_r \Delta_{e,r}^{(2)}(T_{ve}) }
+$$
+
+where 
+
+$$
+\gamma_s = \frac{\rho_s}{\rho M_s}.
+$$
+
+Thermal conductivity is computed in terms of different energy modes. The contribution due to translation  modes is expressed as
+
+$$
+\kappa_t = \frac{15}{4} k_{B} \sum_{s \neq e} 
+\frac{\gamma_s}{\sum_{r \neq e} a_{s,r} \gamma_r \Delta_{s,r}^{(2)}(T_{tr}) + 3.54 \gamma_e \Delta_{s,e}^{(2)}(T_{ve})},
+$$
+
+where
+
+$$
+a_{s,r} = 1 + \frac{\left[1 - (m_s/m_r) \right] \left[ 0.45 - 2.54(m_s/m_r)  \right]  }{\left[1 + (m_s/m_r) \right]^2}
+$$
+
+and where
+
+$$
+m_s = \frac{M_s}{N_{av}}
+$$
+
+with $$N_{av}$$ being Avogadro's Number. The thermal conductivity for the rotational modes is expressed as
+
+$$
+\kappa_r = k_{B} \sum_{s \neq e} 
+\frac{\gamma_s}{\sum_{r \neq e} \gamma_r \Delta_{s,r}^{(1)}(T_{tr}) + \gamma_e \Delta_{s,e}^{(1)}(T_{ve})}.
+$$
+
+The mixture translational/rotational thermal conductivity can then be expressed as
+
+$$
+\kappa_{tr} = \kappa_t + \kappa_r.
+$$
+
+The vibrational/electronic mode thermal conductivity is 
+
+$$
+\kappa_{ve} = k_{B} \frac{C_{ve}}{R} \sum_{s \in molecules} \frac{\gamma_s}
+{\sum_{r \neq e} \gamma_r \Delta_{s,r}^{(1)}(T_{tr}) + \gamma_e \Delta_{s,r}^{(1)}(T_{ve}) }
+$$
+
+and the thermal conductivity for electrons is given by
+
+$$
+\kappa_e = \frac{15}{4} k_{B} \frac{\gamma_e}{\sum_r 1.45 \gamma_r \Delta_{e,r}^{(2)}(T_{ve})}.
+$$
+
+Finally, the binary diffusion coefficient for heavy particles is given by
+
+$$
+D_{s,r} = \frac{k_{B} T_{tr}}{p \Delta_{s,r}^{(1)}(T_{tr})},
+$$
+
+and for electrons,
+
+$$
+D_{e,r} = \frac{k_{B} T_{ve}}{p \Delta_{e,r}^{(1)}(T_{ve})}.
+$$
+
+## Sutherland Viscosity Model ##
+
+In addition to the two models discussed above, there is the option to use a Sutherland model to calculate the flow viscosity. The Sutherland model is not applicable at high temperatures.
+
+In this case, the viscosity is computed as
+
+$$
+\mu = \mu_{0} \left( \frac{T}{T_{0}} \right)^{3/2} \frac{T_0 + S_{\mu}}{T + S_{\mu}},
+$$
+
+where $$T_0$$ is a reference temperature (273.15 K), $$\mu_0$$ is a reference viscosity, and $$S_{\mu}$$ is the Sutherland constant.
+
+If the Sutherland model is selected with a NEMO solver, species diffusion coefficients and thermal conductivity are computed using the models described in the Wilkes-Blottner-Eucken section.
+
+---
+
+# Slip Flow #
+
+| Solver | Version | 
+| --- | --- |
+| `NEMO_NAVIER_STOKES` | 7.0.0 |
+
+SU2-NEMO uses the Maxwell velocity and Smoluchowski temperature jump equations to compute the velocity and temperature of the gas in contact with the surface. The equations are given as
+$$
+v_s =  \frac{2 - \sigma}{\sigma} \lambda \frac{\partial v}{\partial n } + \\
+\frac{3}{4} \frac{\mu}{\rho T} \frac{\partial T}{\partial x},
+$$
+
+and
+$$
+T - T_w =  \frac{2 - \alpha}{\alpha} \lambda \frac{2\gamma}{(\gamma + 1 )Pr} \frac{\partial T}{\partial n}, 
+$$
+
+respectively, where $$\mu$$ is the flow viscosity, $$\rho$$ is the mixture density, $$Pr$$ is the Prandtl number, $$\gamma$$ is the specific heat ratio, $$T$$ is the temperature of the gas, $$T_w$$ is the temperature of the surface, and $$\lambda$$ is the mean free path, calculated as
+$$
+    \lambda = \frac{\mu}{\rho} \frac{\pi}{\sqrt{2RT}}.
+$$
+
+The coefficients $$\sigma$$ and $$\alpha$$ are referred to as the Tangential Momentum Accommodation Coefficient (TMAC) and the Thermal Accommodation Coefficient (TAC), respectively. The values of the accommodation coefficients depend on the physical characteristics of the surface, and are usually determined empirically.
+
+
+---
+
+# Gas-surface Interaction #
+
+| Solver | Version | 
+| --- | --- |
+| `NEMO_NAVIER_STOKES` | 7.0.0 |
+
+Mechanisms of gas-surface interaction are implemented as specific boundary conditions within the SU2-NEMO computational suite. The net result of recombination reactions occurring on the surface is a production of chemical species due to catalytic reactions, $$\dot{\omega}_s^{cat}$$, that must be balanced by the normal diffusive and convective flux at the wall. For steady flow and a no-slip boundary, this can be expressed as
+
+$$
+    \mathbf{J}_s  \cdot \mathbf{n}  = \dot{\omega}_s^{cat}.
+$$
+
+In SU2-NEMO, the chemical production of species due to catalytic processes is included in the computation of the viscous component of the residual, as an additional diffusive flux equivalent to the chemical source term computed due to catalytic reactions. Gradients of species density are then computed directly as part of the SU2-NEMO computational routine, which are used to compute gradients of species mass fraction at wall vertices. 
+
+Options in SU2-NEMO include a super-catalytic wall in which species concentrations are set to specify full recombination to a specified equilibrium concentration (typically the free-stream conditions)
+
+$$
+Y_{w,s} = Y_{eq,s},
+$$
+
+as well as a partiall catalytic wall using a specified reaction efficiency model 
+
+$$
+    \dot{\omega}_s^{cat} = \gamma_s Y_s \rho_w \sqrt{\frac{R_s T_w}{2\pi}},
+$$
+
+where $$\gamma_{s}$$ is the species catalytic efficiency, and represents the proportion of incident mass flux of monatomic species $$s$$ which recombines into its heteronuclear diatomic molecule at the wall.
+
+---

_tutorials/incompressible_flow/Inc_Turbulent_Bend/Inc_Turbulent_Bend.md

@@ -27,7 +27,7 @@ In this tutorial we will simulate the turbulent flow in a 90 degree pipe bend. W
 
 ## Resources
 
-The resources for this tutorial can be found in the [incompressible_flow/Inc_Turbulent_Bend_Wallfunctions](https://github.com/su2code/Tutorials/tree/master/incompressible_flow/Inc_Turbulent_Bend_Wallfunctions) directory in the [tutorial repository](https://github.com/su2code/Tutorials). You will need the configuration file ([sudo.cfg](https://github.com/su2code/Tutorials/tree/master/incompressible_flow/Inc_Turbulent_Bend_Wallfunctions/sudo.cfg)) and the mesh file ([sudo.su2](https://github.com/su2code/Tutorials/tree/master/incompressible_flow/Inc_Turbulent_Bend_Wallfunctions/sudo.su2)). Additionally, the Gmsh geometry is also provided so you can recreate the mesh yourself: [sudo_coarse.geo](https://github.com/su2code/Tutorials/tree/master/incompressible_flow/Inc_Turbulent_Bend_Wallfunctions/sudo_coarse.geo).
+The resources for this tutorial can be found in the [incompressible_flow/Inc_Turbulent_Bend_Wallfunctions](https://github.com/su2code/Tutorials/tree/master/incompressible_flow/Inc_Turbulent_Bend_Wallfunctions) directory in the [tutorial repository](https://github.com/su2code/Tutorials). You will need the configuration file ([sudo.cfg](https://github.com/su2code/Tutorials/tree/master/incompressible_flow/Inc_Turbulent_Bend_Wallfunctions/sudo.cfg)) and the mesh file ([sudo_coarse.su2](https://github.com/su2code/Tutorials/tree/master/incompressible_flow/Inc_Turbulent_Bend_Wallfunctions/sudo_coarse.su2)). Additionally, the Gmsh geometry is also provided so you can recreate the mesh yourself: [sudo_coarse.geo](https://github.com/su2code/Tutorials/tree/master/incompressible_flow/Inc_Turbulent_Bend_Wallfunctions/sudo_coarse.geo).
 
 
 ## Tutorial