[Fluent Inc. Logo] return to home search
next up previous contents index

6.13.1 Inputs at Wall Boundaries


You will enter the following information for a wall boundary:

Defining Thermal Boundary Conditions at Walls

When you are solving the energy equation, you need to define thermal boundary conditions at wall boundaries. Five types of thermal conditions are available:

If the wall zone is a ``two-sided wall'' (a wall that forms the interface between two regions, such as the fluid/solid interface for a conjugate heat transfer problem) a subset of these thermal conditions will be available, but you will also be able to choose whether or not the two sides of the wall are ``coupled''. See below for details.

The inputs for each type of thermal condition are described below. If the wall has a non-zero thickness, you should also set parameters for calculating thin-wall thermal resistance and heat generation in the wall, as described below.

You can model conduction within boundary walls and internal (i.e., two-sided) walls of your model. This type of conduction, called shell conduction, allows you to more conveniently model heat conduction on walls where the wall thickness is small with respect to the overall geometry (e.g., finned heat exchangers or sheet metal in automobile underhoods). Meshing these walls with solid cells would lead to high-aspect-ratio meshes and a significant increase in the total number of cells. See below for details about shell conduction.

Thermal conditions are entered in the Thermal section of the Wall panel (Figure  6.13.1), which is opened from the Boundary Conditions panel (as described in Section  6.1.4). To view the thermal conditions, click the Thermal tab.

Figure 6.13.1: The Wall Panel ( Thermal Section)
\begin{figure} \psfig{file=figures/BCWall.ps} \end{figure}

Heat Flux Boundary Conditions

For a fixed heat flux condition, choose the Heat Flux option under Thermal Conditions. You will then need to set the appropriate value for the heat flux at the wall surface in the Heat Flux field. You can define an adiabatic wall by setting a zero heat flux condition. This is the default condition for all walls.

Temperature Boundary Conditions

To select the fixed temperature condition, choose the Temperature option under Thermal Conditions in the Wall panel. You will need to specify the temperature at the wall surface ( Temperature). The heat transfer to the wall is computed using Equation  6.13-8 or Equation  6.13-9.

Convective Heat Transfer Boundary Conditions

For a convective heat transfer wall boundary, select Convection under Thermal Conditions. Your inputs of Heat Transfer Coefficient and Free Stream Temperature will allow FLUENT to compute the heat transfer to the wall using Equation  6.13-12.

External Radiation Boundary Conditions

If radiation heat transfer from the exterior of your model is of interest, you can enable the Radiation option in the Wall panel and set the External Emissivity and External Radiation Temperature.

Combined Convection and External Radiation Boundary Conditions

You can choose a thermal condition that combines the convection and radiation boundary conditions by selecting the Mixed option. With this thermal condition, you will need to set the Heat Transfer Coefficient, Free Stream Temperature, External Emissivity, and External Radiation Temperature.

Thin-Wall Thermal Resistance Parameters

By default, a wall will have a thickness of zero. You can, however, in conjunction with any of the thermal conditions, model a thin layer of material between two zones. For example, you can model the effect of a piece of sheet metal between two fluid zones, a coating on a solid zone, or contact resistance between two solid regions. FLUENT will solve a 1D conduction equation to compute the thermal resistance offered by the wall and the heat generation in the wall.

To include these effects in the heat transfer calculation you will need to specify the type of material, the thickness of the wall, and the heat generation rate in the wall. Select the material type in the Material Name drop-down list, and specify the thickness in the Wall Thickness field. If you want to check or modify the properties of the selected material, you can click Edit... to open the Material panel; this panel contains just the properties of the selected material, not the full contents of the standard Materials panel.

The thermal resistance of the wall is $\Delta x/k$, where $k$ is the conductivity of the wall material and $\Delta x$ is the wall thickness. The thermal boundary condition you set will be specified on the outside of the thin wall, as shown in Figure  6.13.2, where $T_b$ is the fixed Temperature specified at the wall.

Figure 6.13.2: Thermal Conditions are Specified on the Outside of a Thin Wall
\begin{figure} \psfig{file=figures/bc-wall-outside-temp.ps,height=2.75in} \end{figure}

!! Note that for thin walls, you can specify only a constant thermal conductivity. If you want to use a non-constant thermal conductivity for a wall with non-zero thickness, you should use the shell conduction model (see below for details).

Specify the heat generation rate inside the wall in the Heat Generation Rate field. This option is useful if, for example, you are modeling printed circuit boards where you know the electrical power dissipated in the circuits.

Thermal Conditions for Two-Sided Walls

If the wall zone has a fluid or solid region on each side, it is called a ``two-sided wall''. When you read a grid with this type of wall zone into FLUENT, a ``shadow'' zone will automatically be created so that each side of the wall is a distinct wall zone. In the Wall panel, the shadow zone's name will be shown in the Shadow Face Zone field. You can choose to specify different thermal conditions on each zone, or to couple the two zones:

Shell Conduction in Walls

To enable shell conduction for a wall, turn on the Shell Conduction option in the Wall boundary condition panel. When this option is enabled, FLUENT will compute heat conduction within the wall, in addition to conduction across the wall (which is always computed when the energy equation is solved). The Shell Conduction option will appear in the Wall panel for all walls when solution of the energy equation is active. For a boundary wall, the thermal conditions are applied as described above for thin walls.

!! You must specify a non-zero Wall Thickness in the Wall panel, because the shell conduction model is relevant only for walls with non-zero thickness.

!! Note that the shell conduction model has several limitations:

Defining Wall Motion

Wall boundaries can be either stationary or moving. The stationary boundary condition specifies a fixed wall, whereas the moving boundary condition can be used to specify the translational or rotational velocity of the wall, or the velocity components.

Wall motion conditions are entered in the Momentum section of the Wall panel (Figure  6.13.4), which is opened from the Boundary Conditions panel (as described in Section  6.1.4). To view the wall motion conditions, click the Momentum tab.

Figure 6.13.4: The Wall Panel for a Moving Wall
\begin{figure} \psfig{file=figures/BCWall-move.ps} \end{figure}

Defining a Stationary Wall

For a stationary wall, choose the Stationary Wall option under Wall Motion.

Defining Velocity Conditions for Moving Walls

If you wish to include tangential motion of the wall in your calculation, you need to define the translational or rotational velocity, or the velocity components. Select the Moving Wall option under Wall Motion. The Wall panel will expand, as shown in Figure  6.13.4, to show the wall velocity conditions.

Note that you cannot use the moving wall condition to model problems where the wall has a motion normal to itself. FLUENT will neglect any normal component of wall motion that you specify using the methods below.

Specifying Relative or Absolute Velocity

If the cell zone adjacent to the wall is moving (e.g., if you are using a moving reference frame or a sliding mesh), you can choose to specify velocities relative to the zone motion by enabling the Relative to Adjacent Cell Zone option. If you choose to specify relative velocities, a velocity of zero means that the wall is stationary in the relative frame, and therefore moving at the speed of the adjacent cell zone in the absolute frame. If you choose to specify absolute velocities (by enabling the Absolute option), a velocity of zero means that the wall is stationary in the absolute frame, and therefore moving at the speed of the adjacent cell zone--but in the opposite direction--in the relative reference frame.

!! If you are using one or more moving reference frames, sliding meshes, or mixing planes, and you want the wall to be fixed in the moving frame, it is recommended that you specify relative velocities (the default) rather than absolute velocities. Then, if you modify the speed of the adjacent cell zone, you will not need to make any changes to the wall velocities, as you would if you specified absolute velocities.

Note that if the adjacent cell zone is not moving, the absolute and relative options are equivalent.

Translational Wall Motion

For problems that include linear translational motion of the wall boundary (e.g., a rectangular duct with a moving belt as one wall) you can enable the Translational option and specify the wall's Speed and Direction ( X,Y,Z vector). By default, wall motion is ``disabled'' by the specification of Translational velocity with a Speed of zero.

If you need to define non-linear translational motion, you will need to use the Components option, described below.

Rotational Wall Motion

For problems that include rotational wall motion you can enable the Rotational option and define the rotational Speed about a specified axis. To define the axis, set the Rotation-Axis Direction and Rotation-Axis Origin. This axis is independent of the axis of rotation used by the adjacent cell zone, and independent of any other wall rotation axis. For 3D problems, the axis of rotation is the vector passing through the specified Rotation-Axis Origin and parallel to the vector from (0,0,0) to the ( X,Y,Z) point specified under Rotation-Axis Direction. For 2D problems, you will specify only the Rotation-Axis Origin; the axis of rotation is the $z$-direction vector passing through the specified point. For 2D axisymmetric problems, you will not define the axis: the rotation will always be about the $x$ axis, with the origin at (0,0).

Note that the modeling of tangential rotational motion will be correct only if the wall bounds a surface of revolution about the prescribed axis of rotation (e.g., a circle or cylinder). Note also that rotational motion can be specified for a wall in a stationary reference frame.

Wall Motion Based on Velocity Components

For problems that include linear or non-linear translational motion of the wall boundary you can enable the Components option and specify the X-Velocity, Y-Velocity, and Z-Velocity of the wall. You can define non-linear translational motion using a boundary profile or a user-defined function for the X-Velocity, Y-Velocity, and/or Z-Velocity of the wall.

Wall Motion for Two-Sided Walls

As discussed earlier in this section, when you read a grid with a two-sided wall zone (which forms the interface between fluid/solid regions) into FLUENT, a ``shadow'' zone will automatically be created so that each side of the wall is a distinct wall zone. For two-sided walls, it is possible to specify different motions for the wall and shadow zones, whether or not they are coupled. Note, however, that you cannot specify motion for a wall (or shadow) that is adjacent to a solid zone.

Defining Shear Conditions at Walls

Three types of shear conditions are available:

Note that all moving walls are no-slip; the other shear conditions are relevant only for stationary walls.

The no-slip condition is the default, and it indicates that the fluid sticks to the wall and moves with the same velocity as the wall, if it is moving. The specified shear and Marangoni stress boundary conditions are useful in modeling situations in which the shear stress (rather than the motion of the fluid) is known. Examples of such situations are applied shear stress, slip wall (zero shear stress), and free surface conditions (zero shear stress or shear stress dependent on surface tension gradient). The specified shear boundary condition allows you to specify the $x$, $y$, and $z$ components of the shear stress as constant values or boundary profiles. The Marangoni stress boundary condition allows you to specify the gradient of the surface tension with respect to the temperature at this surface. The shear stress is calculated based on the surface gradient of the temperature and the specified surface tension gradient. The Marangoni stress option is available only for calculations in which the energy equation is being solved.

Shear conditions are entered in the Momentum section of the Wall panel, which is opened from the Boundary Conditions panel (as described in Section  6.1.4). To view the shear conditions, click the Momentum tab.

Modeling No-Slip Walls

You can model a no-slip wall by selecting the No Slip option under Shear Condition. This is the default for all walls in viscous flows.

Specified Shear

In addition to the no-slip wall that is the default for viscous flows, you can model a slip wall by specifying zero or non-zero shear. To specify the shear, select the Specified Shear option under Shear Condition (see Figure  6.13.5). You can then enter $x$, $y$, and $z$ components of shear under Shear Stress. Wall functions for turbulence are not used with the Specified Shear option.

Figure 6.13.5: The Wall Panel for Specified Shear
\begin{figure} \psfig{file=figures/BCWall-shear.ps} \end{figure}

Marangoni Stress

FLUENT can also model shear stresses caused by the variation of surface tension due to temperature. The shear stress applied at the wall is given by

\begin{displaymath} \tau = \frac{d\sigma}{dT} \nabla_s T \end{displaymath} (6.13-1)

where $d \sigma/dT$ is the surface tension gradient with respect to temperature, and $\nabla_s T$ is the surface gradient. This shear stress is then applied to the momentum equation.

To model Marangoni stress for the wall, select the Marangoni Stress option under Shear Condition (see Figure  6.13.6). This option is available only for calculations in which the energy equation is being solved. You can then enter the surface tension gradient ( $d \sigma/dT$ in Equation  6.13-1) in the Surface Tension Gradient field. Wall functions for turbulence are not used with the Marangoni Stress option.

Figure 6.13.6: The Wall Panel for Marangoni Stress
\begin{figure} \psfig{file=figures/BCWall-marangoni.ps} \end{figure}

Modeling Wall Roughness Effects in Turbulent Wall-Bounded Flows

Fluid flows over rough surfaces are encountered in diverse situations. Examples are, among many others, flows over the surfaces of airplanes, ships, turbomachinery, heat exchangers, and piping systems, and atmospheric boundary layers over terrain of varying roughness. Wall roughness affects drag (resistance) and heat and mass transfer on the walls.

If you are modeling a turbulent wall-bounded flow in which the wall roughness effects are considered to be significant, you can include the wall roughness effects through the law-of-the-wall modified for roughness.

Law-of-the-Wall Modified for Roughness

Experiments in roughened pipes and channels indicate that the mean velocity distribution near rough walls, when plotted in the usual semi-logarithmic scale, has the same slope ( $1/\kappa$) but a different intercept (additive constant $B$ in the log-law). Thus, the law-of-the-wall for mean velocity modified for roughness has the form

\begin{displaymath} \frac{u_p u^{*}}{\tau_w/\rho} = \frac{1}{\kappa} \ln (E \frac{\rho u^{*} y_p}{\mu}) - \Delta B \end{displaymath} (6.13-2)

where $u^{*} = C_\mu^{1/4}k^{1/2}$ and

\begin{displaymath} \Delta B = \frac{1}{\kappa} \ln f_r \end{displaymath} (6.13-3)

where $f_r$ is a roughness function that quantifies the shift of the intercept due to roughness effects.

$\Delta B$ depends, in general, on the type (uniform sand, rivets, threads, ribs, mesh-wire, etc.) and size of the roughness. There is no universal roughness function valid for all types of roughness. For a sand-grain roughness and similar types of uniform roughness elements, however, $\Delta B$ has been found to be well-correlated with the nondimensional roughness height, $K_s ^{+} = \rho K_s u^{*}/\mu$, where $K_s$ is the physical roughness height and $u^{*} = C_\mu^{1/4}k^{1/2}$. Analyses of experimental data show that the roughness function is not a single function of $K_s^{+}$, but takes different forms depending on the $K_s^{+}$ value. It has been observed that there are three distinct regimes:

According to the data, roughness effects are negligible in the hydrodynamically smooth regime, but become increasingly important in the transitional regime, and take full effect in the fully rough regime.

In FLUENT, the whole roughness regime is subdivided into the three regimes, and the formulas proposed by Cebeci and Bradshaw based on Nikuradse's data [ 34] are adopted to compute $\Delta B$ for each regime.

For the hydrodynamically smooth regime ( $K_s^{+} \leq 2.25$):

\begin{displaymath} \Delta B = 0 \end{displaymath} (6.13-4)

For the transitional regime ( $2.25 < K_s^{+} \leq 90$):

\begin{displaymath} \Delta B = \frac{1}{\kappa} \ln \left[\frac{K_s^+ -2.25}{87.... ...}\right] \times \sin\left\{0.4258(\ln K_s^{+} - 0.811)\right\} \end{displaymath} (6.13-5)

where $C_{s}$ is a roughness constant, and depends on the type of the roughness.

In the fully rough regime ( $K_s^{+} > 90$):

\begin{displaymath} \Delta B = \frac{1}{\kappa} \ln (1 + C_{s} K_s^{+}) \end{displaymath} (6.13-6)

In the solver, given the roughness parameters, $\Delta B(K_s^{+})$ is evaluated using the corresponding formula (Equation  6.13-4, 6.13-5, or 6.13-6). The modified law-of-the-wall in Equation  6.13-2 is then used to evaluate the shear stress at the wall and other wall functions for the mean temperature and turbulent quantities.

Setting the Roughness Parameters

The roughness parameters are in the Momentum section of the Wall panel (see Figure  6.13.6), which is opened from the Boundary Conditions panel (as described in Section  6.1.4). To view the wall roughness parameters, click the Momentum tab.

To model the wall roughness effects, you must specify two roughness parameters: the Roughness Height, $K_s$, and the Roughness Constant, $C_{s}$. The default roughness height ( $K_s$) is zero, which corresponds to smooth walls. For the roughness to take effect, you must specify a non-zero value for $K_s$. For a uniform sand-grain roughness, the height of the sand-grain can simply be taken for $K_s$. For a non-uniform sand-grain, however, the mean diameter ( $D_{50}$) would be a more meaningful roughness height. For other types of roughness, an ``equivalent'' sand-grain roughness height should be used for $K_s$.

Choosing a proper roughness constant ( $C_{s}$) is dictated mainly by the type of the given roughness. The default roughness constant ( $C_{s} = 0.5$) was determined so that, when used with $k$- $\epsilon$ turbulence models, it reproduces Nikuradse's resistance data for pipes roughened with tightly-packed, uniform sand-grain roughness. You may need to adjust the roughness constant when the roughness you want to model departs much from uniform sand-grain. For instance, there is some experimental evidence that, for non-uniform sand-grains, ribs, and wire-mesh roughness, a higher value ( $C_{s} = 0.5 \sim 1.0$) is more appropriate. Unfortunately, a clear guideline for choosing $C_{s}$ for arbitrary types of roughness is not available.

Note that it is not physically meaningful to have a mesh size such that the wall-adjacent cell is smaller than the roughness height. For best results, make sure that the distance from the wall to the centroid of the wall-adjacent cell is greater than $K_s$.

Defining Species Boundary Conditions for Walls

By default, a zero-gradient condition for all species is assumed at walls (except for species that participate in surface reactions), but it is also possible to specify species mass fractions at walls. That is, Dirichlet boundary conditions such as those that are specified at inlets can be used at walls as well.

If you wish to retain the default zero-gradient condition for a species, no inputs are required. If you want to specify the mass fraction for a species at the wall, the steps are as follows:

Click the Species tab in the Wall panel to view the species boundary conditions for the wall (see Figure  6.13.7).

Under Species Boundary Condition, select Specified Mass Fraction (rather than Zero Diffusive Flux) in the drop-down list to the right of the species name. The panel will expand to include space for Species Mass Fractions.

Figure 6.13.7: The Wall Panel for Species Boundary Condition Input
\begin{figure} \psfig{file=figures/spmf-wall.ps} \end{figure}

Under Species Mass Fractions, specify the mass fraction for the species.

The boundary condition type for each species is specified separately, so you can choose to use different methods for different species.

If you are modeling species transport with reactions, you can, alternatively, enable a reaction mechanism at a wall by turning on the Reaction option and selecting an available mechanism from the Reaction Mechanisms drop-down list. See Section  13.1.4 more information about defining reaction mechanisms.

Defining Reaction Boundary Conditions for Walls

If you have enabled the modeling of wall surface reactions in the Species Model panel, you can indicate whether or not surface reactions should be activated for the wall. In the Species section of the Wall panel (Figure  6.13.7), turn the Surface Reactions option on or off.

Note that a zero-gradient condition is assumed at the wall for species that do not participate in any surface reactions.

Defining Radiation Boundary Conditions for Walls

If you are using the P-1 radiation model, the DTRM, the DO model, or the surface-to-surface model, you will need to set the emissivity of the wall ( Internal Emissivity) in the Radiation section of the Wall panel. If you are using the Rosseland model you do not need to set the emissivity, because FLUENT assumes the emissivity is 1. If you are using the DO model you will also need to define the wall as diffuse, specular, or semi-transparent. See Section  11.3.16 for details.

Defining Discrete Phase Boundary Conditions for Walls

If you are modeling a discrete phase of particles, you can set the fate of particle trajectories at the wall in the DPM section of the Wall panel. See Section  21.10 for details.

Defining Multiphase Boundary Conditions for Walls

If you are using the VOF model and you are modeling wall adhesion, you can specify the contact angle for each pair of phases at the wall in the Momentum section of the Wall panel. See Section  22.6.16 for details.

Defining Boundary Conditions for User-Defined Scalars

If you are using any user-defined scalars in your simulation, you can specify boundary conditions for them in the UDS section of the Wall panel. See the separate UDF Manual for details about user-defined scalars.

next up previous contents index Previous: 6.13 Wall Boundary Conditions
Up: 6.13 Wall Boundary Conditions
Next: 6.13.2 Default Settings at
© Fluent Inc. 2003-01-25