Merge pull request #8582 from rmcdermo/master

FDS Verification Guide: clarifications to mass balance reac section

Manuals/FDS_Verification_Guide/FDS_Verification_Guide.tex

@@ -2127,7 +2127,7 @@ \subsection{Mass Flux through Domain Boundaries (\texorpdfstring{\textct{mass\_f
 \subsection{Mass Balance for Species in a Reacting Flow (\texorpdfstring{\textct{mass\_balance\_reac}}{mass\_balance\_reac})}
 \label{mass_balance_reac}
 
-A 0.25 m$^2$ propane burner ramps up to 500 kW linearly in 10 s and holds this heat release rate for another 10 s.  Using a mass file, appropriate wall flux devices, surface integrals, and reaction source term devices, we perform a detailed mass balance on all the primitive species in the reaction.  The results are presented in Fig.~\ref{fig:mass_balance_reac}.  The mass balance can be understood through the mnemonic, ``Accumulation equals In minus Out plus Generation''.  For each species, we plot the ``accumulation'' term (first-order numerical differentiation of the columns in the {\ct MASS\_FILE}, the ``in-out'' term (surface intgral over the domain boundaries), ``generation'' term (reaction source term), and the sum of these as the ``balance''.
+A 0.25 m$^2$ propane burner ramps up to 500 kW linearly in 10~s and holds this heat release rate for another 10~s.  Using a mass file, appropriate wall flux devices with surface integrals, and reaction source term devices with volume integrals, we perform a detailed mass balance on all the primitive species in the reaction.  The results are presented in Fig.~\ref{fig:mass_balance_reac}.  The mass balance can be understood through the mnemonic, ``Accumulation equals In minus Out plus Generation''.  For each species, we plot the ``accumulation'' term (first-order numerical differentiation of the columns in the {\ct MASS\_FILE}, which are volume integrated over the domain), the ``in-out'' term (surface intgral of the mass fluxes over the domain boundaries), the ``generation'' term (volume integral of the reaction source term over the domain), and the sum of these as the ``balance'' (the balance should be zero).
 
 \begin{figure}[ht]
 \begin{tabular*}{\textwidth}{lr}

Utilities/Matlab/scripts/mass_balance_reac.m

@@ -41,11 +41,11 @@
 
 bal = -dmdt + mdot_out + gen;
 
-plot(t,zeros(1,length(t)),'k-'); hold on
-H(1)=plot(t,dmdt,'g-');
+%plot(t,zeros(1,length(t)),'k-'); hold on
+H(1)=plot(t,dmdt,'g-'); hold on
 H(2)=plot(t,mdot_out,'b-');
 H(3)=plot(t,gen,'r-');
-H(4)=plot(t,bal,'c--');
+H(4)=plot(t,bal,'k-');
 
 ylabel('mass flow (kg/s)', 'FontSize',Label_Font_Size)
 xlabel('time (s)', 'FontSize',Label_Font_Size)