You signed in with another tab or window. Reload to refresh your session.You signed out in another tab or window. Reload to refresh your session.You switched accounts on another tab or window. Reload to refresh your session.Dismiss alert
all the upwind intensities are determined from solid or gas phase boundaries. In theory, iterations are needed if the reflections
536
536
or scattering are important, or if the scenario is optically very thick. Currently, no iterations are made.
537
537
538
+
To mitigate ray effects in the wall radiative heat flux, an option is implemented to randomly rotate the local meridian ($x$), azimuthal ($y$), and axial ($z$) coordinate system at each selected time step with respect to which the angular coordinates, $\theta$ and $\phi$ are evaluated. When random rotation is enabled, the discrete angular integration coefficients (e.g., $D_{xu}^l$) are no longer exactly symmetric with respect to the mesh faces. Consequently, the sum of $D_{xu}^l$ (and the corresponding coefficients for the other coordinate directions) over a face may deviate slightly from the theoretical value of $\pi$.
539
+
538
540
\subsection{Spatial Discretization in Cut Cells}
539
541
540
542
In the case of complex geometry where irregular cut-cells may exist (see Sec.~\ref{sec:unstructured_geometry}) an unstructured version of the RTE in Eq.~\ref{RTEdiscrete3} is adopted
Copy file name to clipboardExpand all lines: Manuals/FDS_User_Guide/FDS_User_Guide.tex
+28-10Lines changed: 28 additions & 10 deletions
Original file line number
Diff line number
Diff line change
@@ -2434,23 +2434,32 @@ \subsection{Walls with Different Materials Front and Back}
2434
2434
\subsection{Specified Internal Heat Source}
2435
2435
2436
2436
\label{info:INTERNAL_HEAT_SOURCE}
2437
+
\label{internal_heating}
2437
2438
2438
2439
The condensed phase heat conduction equation has a source term that describes the internal sources and sinks of energy.
2439
2440
There are three types of sources that contribute to this term:
2440
2441
heats of reaction for the pyrolysis (see Sec.~\ref{info:solid_pyrolysis}),
2441
2442
internal absorption and emission of radiation (see Sec.~\ref{info:liquid_fuels}), and the source specified by the user.
2442
2443
An example of the case where specified heat source could be needed is the heating of electrical cables due to internal current.
2443
2444
2444
-
The internal source term for each layer of the surface is specified using \ct{INTERNAL_HEAT_SOURCE} on the \ct{SURF} line. Its units are \unit{kW/m^3} and the default value is zero. An optional time ramp can be specified for each layer's heat source using \ct{RAMP_IHS}. In the example below, the cylindrical surface describing a cable consists of an outer plastic layer and inner core of metal. The metal core is heated with a power of 300~\unit{kW/m^3}.
2445
+
The internal source term for each layer of the surface is specified using \ct{INTERNAL_HEAT_SOURCE} on the \ct{SURF} line. Its units are \unit{kW/m^3} and the default value is zero. An optional time ramp can be specified for each layer's heat source using \ct{RAMP_IHS}. In the example below, the cylindrical surface describing a 10~cm long cable segment consists of an outer plastic layer and inner core of metal. The metal core is heated with a power of 300~\unit{kW/m^3}.
2445
2446
\begin{lstlisting}
2446
-
&SURF ID = 'Cable'
2447
-
THICKNESS = 0.002,0.008
2448
-
MATL_ID(1,1) = 'PLASTIC'
2449
-
MATL_ID(2,1) = 'METAL'
2450
-
GEOMETRY = 'CYLINDRICAL'
2451
-
LENGTH = 0.1
2452
-
INTERNAL_HEAT_SOURCE = 0.,300. /
2447
+
&SURF ID = 'Cable'
2448
+
THICKNESS = 0.002,0.008
2449
+
MATL_ID(1,1) = 'PLASTIC'
2450
+
MATL_ID(2,1) = 'METAL'
2451
+
GEOMETRY = 'CYLINDRICAL'
2452
+
LENGTH = 0.1
2453
+
INTERNAL_HEAT_SOURCE = 0.,300. /
2453
2454
\end{lstlisting}
2455
+
Figure~\ref{fig:internal_heating} displays the heat generated by 10 of these cable segments. The exact value is 300~\unit{kW/m^3} multiplied by the volume of the metal within the cable segment, $2 \times 10^{-5}$~\unit{m^3}, multiplied by the number of segments, 10, which equals approximately 0.06~kW.
\caption[Results of the \ct{internal_heating} test case]{Heating rate of a set of 10 cable segments. The dashed line represents the heat generated by the cable and the dotted line the heat flowing out of the computational domain.}
2460
+
\label{fig:internal_heating}
2461
+
\end{figure}
2462
+
2454
2463
2455
2464
\subsection{Non-Planar Walls and Targets}
2456
2465
@@ -2533,7 +2542,8 @@ \subsection{Limitations}
2533
2542
\item Avoid contact between 3-D and 1-D solids. If two sides of a 3-D solid touch 1-D solids, there will be no lateral heat conduction computed in that particular direction.
2534
2543
\item If your 3-D obstruction extends beyond meshes that abut, FDS uses a special algorithm to identify all the meshes where this obstruction lives, and also those obstructions connected to it. However, if this algorithm fails to detect all the meshes and you receive an error stating that there is a problem, add the parameter \ct{NEIGHBOR_SEPARATION_DISTANCE} to the \ct{MISC} line. Any mesh within this distance of another mesh will share geometry information for use in the 3-D heat conduction calculation. If you set this parameter to a value larger than the width of the computational domain, then all meshes will establish communication channels for exchanging boundary information.
2535
2544
\item By default, the interior nodes are clustered near the surface and stretched out deeper within the solid. If you want to maintain uniform spacing, set \ct{CELL_SIZE} on the \ct{SURF} or \ct{OBST} line to indicate the desired interior node spacing. The \ct{CELL_SIZE} is typically chosen to be comparable to the gas phase cells. If the obstruction is thin; that is, less than one gas phase cell thick, the specified \ct{CELL_SIZE} will only apply to the heat conduction in the transverse, not normal, direction. The normal direction gridding will be controlled by the parameters \ct{STRETCH_FACTOR} and \ct{CELL_SIZE_FACTOR}. This may be useful in cases where the specified \ct{CELL_SIZE} is too coarse to resolve variations in surface definition along the normal dimension (see Section~\ref{checkerboard} for an example).
2536
-
\item \ct{HT3D} cannot be applied to an \ct{OBST} that is to \ct{BURN_AWAY}. In addition, if the solid undergoes significant shrinking or swelling, do not use \ct{HT3D}. ``Significant'' means that the number of internal cells changes, in which case the 3-D nodal structure breaks down. To determine if this happens, use the \ct{PROF}ile output feature (Sec.~\ref{info:PROF}) to visualize profiles of internal temperature or other solid phase quantities. These output files contain the internal node coordinates as a function of time.
2545
+
\item \ct{HT3D=T} cannot be applied to an \ct{OBST} that is to \ct{BURN_AWAY}. In addition, if the solid undergoes significant shrinking or swelling, do not use \ct{HT3D}. ``Significant'' means that the number of internal cells changes, in which case the 3-D nodal structure breaks down. To determine if this happens, use the \ct{PROF}ile output feature (Sec.~\ref{info:PROF}) to visualize profiles of internal temperature or other solid phase quantities. These output files contain the internal node coordinates as a function of time.
2546
+
\item \ct{HT3D=T} or \ct{VARIABLE_THICKNESS=T} cannot be applied to an \ct{OBST} that is to be created or removed during the simulation, unless the \ct{OBST} does not touch any other \ct{OBST} or the boundary of the domain. This rule does not apply to two \ct{OBST}s that are controlled by the exact same device or controller; that is, the two \ct{OBST}s are created or removed at the same time.
2537
2547
\item If a \ct{SURF} line specifies either \ct{HT3D=T} or \ct{VARIABLE_THICKNESS=T} and also specifies an \ct{HRRPUA}, \ct{MLPUA}, or \ct{MASS_FLUX}, you must specify a \ct{MATL_ID} on the \ct{SURF} line and the appropriate \ct{OBST} lines. The reason for doubly specifying the \ct{MATL_ID} is so that the \ct{SURF} line can be set up properly. Note that the specification of \ct{HRRPUA} or similar on the surface of a 3-D or variably thick solid means that no obstruction making up the solid can have specified internal reactions, i.e. pyrolysis.
2538
2548
\item Issues regarding pyrolysis and material transport in 3-D solids is discussed in Section~\ref{info:LAYER_DIVIDE}.
2539
2549
\end{enumerate}
@@ -3814,6 +3824,7 @@ \subsection{HVAC Fan Parameters}
3814
3824
\end{lstlisting}
3815
3825
where:
3816
3826
\begin{itemize}
3827
+
\item \ct{CURVE_TEMP} is the temperature the fan curve is based upon. FDS assumes a default value of \SI{20}{\degreeCelsius}.
3817
3828
\item \ct{LOSS} is the loss coefficient for flow through the fan when it is not operational. FDS assumes a default value of 1.
3818
3829
\item \ct{MAX_FLOW} is the maximum volumetric flow of the fan in \unit{m^3/s}. This input activates a quadratic fan model. Value must be > 0.
3819
3830
\item \ct{MAX_PRESSURE} is the stall pressure of the fan in units of Pa. This input activates a quadratic fan model. Value must be > 0.
@@ -3860,6 +3871,8 @@ \subsubsection{Fan Curves}
3860
3871
\end{center}
3861
3872
\end{figure}
3862
3873
3874
+
FDS will automatically adjust the pressures of the fan curve to reflect the fan affinity law for density. This law states that the fan curve will shift in proportion to the ratio of the current duct density to the as tested duct density. The as tested duct density is computed assuming air as the gas, standard pressure, and \ct{CURVE_TEMP}.
3875
+
3863
3876
\subsubsection{Jet Fans}
3864
3877
3865
3878
Fans do not have to be mounted on a solid wall, like a supply or an exhaust fan. If you just want to blow gases in a particular direction, create an obstruction \ct{OBST}, at least one cell thick, and apply to it \ct{VENT} lines that are associated with a simple HVAC system. This allows hot, smokey gases to pass through the obstruction, much like a free-standing fan. See the example case \ct{jet_fan.fds} which places a louvered fan (blowing diagonally down) near a fire (see Fig.~\ref{fig:Jet_Fan}).
@@ -5853,6 +5866,8 @@ \section{Spatial and Temporal Resolution of the Radiation Transport Solver}
5853
5866
5854
5867
The frequency of calls to the radiation solver can be changed from every 3 time steps with an integer called \ct{TIME_STEP_INCREMENT}. The increment over which the angles are updated can be changed from 5 with the integer called \ct{ANGLE_INCREMENT}. If \ct{TIME_STEP_INCREMENT} and \ct{ANGLE_INCREMENT} are both set to 1, the radiation field is completely updated in a single time step, but the cost of the calculation increases significantly. By default, the radiation transport equation is fully updated every 15 time steps. Given the relatively small time steps typically dictated by the CFL constraint, 15 time steps is usually still a relatively short interval of time. Increasing the temporal resolution of the radiation solver rarely adds to the overall accuracy of the calculation. Spatial resolution is far more important.
5855
5868
5869
+
Additionally, to mitigate ray effects in the wall radiative heat flux, the option \ct{RANDOMIZE_RADIATION_DIRECTIONS} randomly rotates the local meridian ($x$), azimuthal ($y$), and axial ($z$) coordinate system used to define the angular coordinates $\theta$ and $\phi$. If \ct{RANDOMIZE_RADIATION_DIRECTIONS} is true, this rotation is performed after each complete sweep of the angular directions.
5870
+
5856
5871
If you are using multiple meshes, the radiation solver cannot transfer energy from mesh to mesh within a single time step. If you notice an obvious delay in the propagation of radiative intensity from one mesh to another, you can increase the number of times the radiative intensity is updated within a single time step using \ct{RADIATION_ITERATIONS}, which is 1 by default. You rarely need to do this, unless it is obvious from the animation of various slice files.
5857
5872
5858
5873
The radiation solver is called before the start of the calculation to establish the radiation field in the event that you specify something to have a non-ambient temperature initially. By default, the radiation and wall boundary routines are iterated three times to establish thermal equilibrium. To change the number of iterations, set \ct{INITIAL_RADIATION_ITERATIONS} on the \ct{RADI} line.
@@ -12402,6 +12417,7 @@ \section{\texorpdfstring{{\tt HVAC}}{HVAC} (HVAC System Definition)}
12402
12417
\ct{COOLANT_SPECIFIC_HEAT} & Real & Section~\ref{info:HVACaircoil} & \unit{kJ/(kg.K)} & \\ \hline
0 commit comments