32.3.2 Boundary conditions in Abaqus/CFD

Products: Abaqus/CFD  Abaqus/CAE  

References

Overview

Boundary conditions:

Computational fluid dynamics problems typically require the prescription of multiple variables such as pressure, temperature, and velocity for boundary conditions. In practice, boundary conditions tend to appear together to collectively define a physical behavior; e.g., no-slip/no-penetration conditions at a wall. In contrast, Neumann conditions (e.g., prescribed heat flux) are specified as loads (see Specifying surface-based distributed heat fluxes” in “Thermal loads, Section 32.4.4). In the absence of a prescribed boundary condition or load, the default behavior for Abaqus/CFD is to enforce a homogeneous (zero) Neumann condition. For example, if the temperature is not specified at a wall, the default behavior is to automatically specify a perfectly insulated boundary; i.e., zero normal heat flux. Similarly, if the velocity is not prescribed, the normal derivative of the velocity is set to zero.

In Abaqus/CFD combinations of boundary conditions that represent an inflow, outflow, or wall behavior are grouped collectively for ease of use. These types of boundary conditions are prescribed using Abaqus/CAE (for more information, see Using the boundary condition editors, Section 16.10 of the Abaqus/CAE User's Manual).

Active degrees of freedom

In Abaqus/CFD the active fields (degrees of freedom) are determined by the analysis procedure and the options specified, such as turbulence models and auxiliary transport equations. You specify a boundary condition type to identify the degree of freedom for a fluid boundary condition. Element-based and node-based degrees of freedom and the analysis procedure and additional options required for activation, if any, are listed in Table 32.3.2–1 and Table 32.3.2–2, respectively.

Table 32.3.2–1 Element-based degrees of freedom and activation options for fluid boundary conditions.

Boundary condition type DescriptionIncompressible flow
TEMPFluid temperatureEnergy equation
TEMPnFluid temperature on face nEnergy equation
TURBEPSTurbulent energy dissipation rate ()RNG - model
TURBEPSnTurbulent energy dissipation rate () on face nRNG - model
TURBKETurbulent kinetic energy ()RNG - model
TURBKEnTurbulent kinetic energy () on face nRNG - model
TURBNUTurbulent kinematic eddy viscositySpalart-Allmaras model
TURBNUnTurbulent kinematic eddy viscosity on face nSpalart-Allmaras model
VELXx-velocity
VELXnx-velocity on face n
VELYy-velocity
VELYny-velocity on face n
VELZz-velocity
VELZnz-velocity on face n

Table 32.3.2–2 Node-based degrees of freedom and activation options for fluid boundary conditions.

Boundary condition type DescriptionIncompressible flow
PFluid pressure
DISTWall-distance normal function

Prescribing inflow and outflow boundary conditions

You can specify boundary conditions to describe the flow behavior where fluid enters the analysis domain and where the fluid leaves the analysis domain.

Input File Usage:          Use the following option to define the inflow and outflow boundary conditions at surfaces:
*FLUID BOUNDARY, TYPE=SURFACE
surface name, boundary condition type label, magnitude

where boundary condition type label is VELX, VELY, VELZ, TEMP, TURBKE, TURBEPS, TURBNU, or P.

Use the following option to define distributed inflow and outflow boundary conditions at element faces:

*FLUID BOUNDARY, TYPE=ELEMENT
element set label, boundary condition type label, magnitude

where boundary condition type label is VELXn, VELYn, VELZn, TEMPn, TURBKEn, TURBEPSn, or TURBNUn.

Use the following option to define distributed inflow and outflow boundary conditions at nodes:

*FLUID BOUNDARY, TYPE=NODE
node set label, P, magnitude

Abaqus/CAE Usage:   Use the following option to define the inflow and outflow boundary conditions at surfaces:

Load module: Create Boundary Condition: Step: flow_step: Category: Fluid: Fluid inlet/outlet: select inlet regions or outlet regions; and specify momentum (pressure or velocity), thermal energy (temperature), and turbulence conditions at the inlet or outlet

Defining distributed inflow and outflow boundary conditions at element faces is supported in Abaqus/CAE only for velocity boundary conditions.

Use the following option:

Load module: Create Boundary Condition: Step: flow_step: Category: Fluid: Fluid inlet/outlet: select inlet regions or outlet regions; Momentum: toggle on Specify, and choose Velocity; Distribution: select an analytical field

Defining distributed inflow and outflow boundary conditions at nodes is not supported in Abaqus/CAE.


Inflow boundary conditions

An inflow boundary condition is used to describe the flow behavior at a surface where fluid enters the analysis domain. For incompressible flows, inflow conditions can be prescribed for velocity or pressure, temperature, and turbulence variables. If boundary conditions are not specified explicitly for a variable, a homogeneous Neumann condition is assumed automatically. This corresponds to permitting the variable (e.g., temperature) to vary at the inflow and the incoming fluid to correspond to that local variable. Similarly, if pressure is not specified, its normal derivative at the inflow surface is automatically set to zero. The velocity components can be prescribed independently.

Outflow boundary conditions

An outflow boundary corresponds to a surface where the fluid flow leaves the analysis domain. In Abaqus/CFD outflow conditions are most frequently associated with a specified pressure. However, all other flow variables can be prescribed at an outflow boundary as well. Similar to an inflow boundary, when a variable is not specified, its normal derivative is assumed to be zero. As such, convective outflows carry their quantities out of the domain at a fixed level, resulting in essentially nonreflecting boundaries.

Prescribing wall boundary conditions

Wall boundary conditions are typically associated with the no-slip/no-penetration behavior at a solid surface. However, the behavior at a solid wall may also require the prescription of temperature and, optionally, turbulence variables depending on the flow conditions. In situations where a wall heat flux is required, a heat flux loading must be prescribed in addition to the wall boundary conditions.

Depending on the physical properties of the wall, the wall boundary conditions can be modified to achieve a variety of physical behaviors that include slip, no-slip, infiltration, symmetry, etc.

Input File Usage:          Use the following option to define wall boundary conditions at surfaces:
*FLUID BOUNDARY, TYPE=SURFACE
surface name, boundary condition type label, magnitude

where boundary condition type label is VELX, VELY, VELZ, TEMP, TURBKE, TURBEPS, TURBNU, P or DIST.

Use the following option to define distributed wall boundary conditions at element faces:

*FLUID BOUNDARY, TYPE=ELEMENT
element set label, boundary condition type label, magnitude

where boundary condition type label is VELXn, VELYn, VELZn, TEMPn, TURBKEn, TURBEPSn, or TURBNUn.

Use the following option to define distributed wall boundary conditions at nodes:

*FLUID BOUNDARY, TYPE=NODE
node set label, P, magnitude

For example, use the following settings for a no-slip/no-penetration wall that is not moving and with the Spalart-Allmaras turbulence model active (wall-normal distance boundary condition and turbulent eddy viscosity set to zero at the wall):

*FLUID BOUNDARY, TYPE=SURFACE 
surface name, DIST, 0
surface name, VELX, 0
surface name, VELY, 0
surface name, VELZ, 0
surface name, TURBNU, 0

Abaqus/CAE Usage:   Use the following option to define wall boundary conditions at surfaces:

Load module: Create Boundary Condition: Step: flow_step: Category: Fluid: Fluid wall condition: select regions; select Condition: No slip, Shear, or Infiltration; and specify velocity, thermal energy (temperature), and turbulence conditions at the wall

Defining distributed wall boundary conditions at elements is supported in Abaqus/CAE only for velocity boundary conditions at a slip wall or infiltration wall.

Use the following option to define distributed wall boundary conditions at elements:

Load module: Create Boundary Condition: Step: flow_step: Category: Fluid: Fluid wall condition: select regions; Velocity: Distribution: select an analytical field

Defining distributed wall boundary conditions at nodes is not supported in Abaqus/CAE.


No-slip/no-penetration wall

A no-slip (and no-penetration) wall is a surface where the fluid adheres to the wall without penetrating it. No-slip/no-penetration conditions are prescribed by setting all velocity components equal to the wall velocity (zero if the wall is not moving). If a turbulence model is specified, the wall-normal distance boundary condition must be set to zero at the wall. The boundary conditions for the different turbulence variables depend on the model selected. For the Spalart-Allmaras model, the turbulent eddy viscosity, , is set to zero at the wall. For the RNG k model, the wall boundary conditions are automatically implemented by the solver using the wall-function approach; no user settings for k or are required because they are prescribed automatically.

Slip wall

A slip wall is a surface where the fluid does not adhere to the wall but cannot penetrate it. This wall condition is modeled by specifying the wall-normal fluid velocity equal to the wall velocity (zero if the wall is not moving). This situation also represents a symmetry condition for fluid flow since the in-plane velocities can vary, but the out-of-plane velocity is zero. In cases where a moving boundary is being considered, an associated set of mesh displacement boundary conditions must be prescribed in conjunction with the surface fluid velocity to achieve the proper behavior. If a turbulence model is specified, the wall-normal distance boundary condition must be set to zero at the wall.

Infiltration wall

Infiltration at a surface permits the fluid to penetrate the surface while maintaining the no-slip condition. This wall condition is modeled by specifying the wall-normal velocity equal to the velocity representing the infiltration velocity, while the wall-tangent fluid velocity is equal to the wall velocity (zero if the wall is not moving). In the special case when a turbulence model is implemented, the wall-normal distance boundary condition must be set to zero at the wall. If the Spalart-Allmaras turbulence model is enabled, you can specify the value of the Spalart-Allmaras turbulent eddy viscosity, , that is allowed at the wall due to infiltration. If the RNG k model is implemented, you can prescribe values at the wall for the turbulent kinetic energy, k, and the dissipation rate, .

Prescribed temperature

Temperatures can be prescribed at a wall. By default, if no temperature is prescribed at a wall, a perfectly insulated boundary is specified automatically. For multiphysics applications such as conjugate heat transfer, a variable temperature condition is imposed automatically using a co-simulation region (for more information, see Preparing an Abaqus analysis for co-simulation, Section 16.1.2).

Prescribed displacement

Abaqus/CFD provides the capability to perform both deforming-mesh and fluid-structure interaction (FSI) simulations using an arbitrary Lagrangian Eulerian (ALE) methodology for the fluid flow. For FSI and deforming-mesh problems, typically some portion of the fluid domain is deformed consistent with a boundary motion. To manage the mesh motion, you must prescribe displacement boundary conditions on the mesh. For FSI problems, displacement boundary conditions are not permitted at the co-simulation region because these conditions are prescribed automatically.

Input File Usage:          
*BOUNDARY
node or node set,  first degree of freedom, last degree of freedom, magnitude

where first degree of freedom is 1 for the x-displacement, 2 for the y-displacement, or 3 for the z-displacement.

Abaqus/CAE Usage:   

Load module: Create Boundary Condition: Step: flow_step: Category: Mechanical: Displacement/Rotation: select regions and toggle on the degree or degrees of freedom


Defining boundary conditions that vary with time

The prescribed magnitude of the boundary conditions can vary with time during a step according to an amplitude definition (Amplitude curves, Section 32.1.2).

Input File Usage:          Use both of the following options to define the prescribed displacement at a moving boundary:
*AMPLITUDE, NAME=name
*BOUNDARY, AMPLITUDE=name

Use both of the following options to define inflow and outflow boundary conditions and wall boundary conditions that vary with time:

*AMPLITUDE, NAME=name
*FLUID BOUNDARY, AMPLITUDE=name

Abaqus/CAE Usage:   

Load or Interaction module: Create Amplitude: Name: amplitude_name
Load module: Create Boundary Condition: Step: flow_step: boundary condition; Amplitude: amplitude_name


Boundary condition propagation

By default, all boundary conditions defined in the previous general analysis step remain unchanged in the subsequent general step. You define the boundary conditions in effect for a given step relative to the preexisting boundary conditions. At each new step the existing boundary conditions can be modified and additional boundary conditions can be specified. Alternatively, you can release all previously applied boundary conditions in a step and specify new ones. In this case any boundary conditions that are to be retained must be respecified.

Modifying boundary conditions

When you modify an existing boundary condition, the node or node set must be specified in exactly the same way as previously. For example, if a boundary condition is specified for a node set in one step and for an individual node contained in the set in another step, Abaqus issues an error. You must remove the boundary condition and respecify it to change the way the node or node set is specified.

Input File Usage:          Use one of the following options to modify an existing boundary condition or to specify an additional boundary condition:
*BOUNDARY
*BOUNDARY, OP=MOD
*FLUID BOUNDARY
*FLUID BOUNDARY, OP=MOD

Abaqus/CAE Usage:   

Load module: Create Boundary Condition or Boundary Condition Manager: Edit


Removing boundary conditions

If you choose to remove any boundary condition in a step, no boundary conditions will be propagated from the previous general step. Therefore, all boundary conditions that are in effect during this step must be respecified.

Setting a boundary condition to zero is not the same as removing it.

Input File Usage:          Use one of the following options to release all previously applied boundary conditions and to specify new boundary conditions:
*BOUNDARY, OP=NEW

If the OP=NEW parameter is used on any *BOUNDARY option within a step, it must be used on all *BOUNDARY options in the step.

*FLUID BOUNDARY, OP=NEW

If the OP=NEW parameter is used on any *FLUID BOUNDARY option within a step, it must be used on all *FLUID BOUNDARY options in the step.


Abaqus/CAE Usage:   Use the following option to remove a boundary condition within a step:

Load module: Boundary Condition Manager: Deactivate

Abaqus/CAE automatically respecifies any boundary conditions that should remain in effect during this step.