
Overview
The general nonreflecting boundary conditions in ANSYS FLUENT are based on characteristic wave relations derived from the Euler equations, and applied only on pressureoutlet boundary conditions. To obtain the primitive flow quantities ( ) on the pressureoutlet, reformulated Euler equations are solved on the boundary of the domain in an algorithm similar to the flow equations applied to the interior of the domain.
Unlike the turbospecific NRBC method presented in the previous section, the general NRBC method is not restricted by geometric constraints or mesh type. However, good cell skewness near the boundaries where the NRBC will be applied is advisable for a more stable, converged solution. In addition, the general NRBCs can be applied to steady or transient flows as long as the compressible idealgas law is used.
Restrictions and Limitations
Note the following restrictions and limitations on the general NRBCs:

The general NRBC should not be used with the wet steam or real gas models.


If you switch from Turbospecific NRBC to general NRBC or vice versa, then make sure you switch off one NRBC model before turning on the next. You cannot use both NRBC models at the same time.

Theory
General NRBCs are derived by first recasting the Euler equations in an orthogonal coordinate system ( ) such that one of the coordinates, , is normal to the boundary Figure 7.4.5. The characteristic analysis [ 88, 89] is then used to modify terms corresponding to waves propagating in the normal direction. When doing so, a system of equations can be written to describe the wave propagation as follows:
Where , and and , and are the velocity components in the coordinate system ( , , ). The equations above are solved on pressureoutlet boundaries, along with the interior governing flow equations, using similar time stepping algorithms to obtain the values of the primitive flow variables ( ).

Note that a transformation between the local orthogonal coordinate system (
,
,
) and the global Cartesian system (X, Y, Z) must be defined on each face on the boundary to obtain the velocity components (
,
,
) in a global Cartesian system.

The terms in the transformed Euler equations contain the outgoing and incoming characteristic wave amplitudes, , and are defined as follows:
From characteristic analyses, the wave amplitudes, , are given by:
The outgoing and incoming characteristic waves are associated with the characteristic velocities of the system (i.e eigenvalues), , as seen in Figure 7.4.6. These eigenvalues are given by:
For subsonic flow leaving a pressureoutlet boundary, four waves leave the domain (associated with positive eigenvalues , , , and ) and one enters the domain (associated with negative eigenvalue ).
To solve Equations 7.444 on a pressureoutlet boundary, the values of , , and must be first determined from Equations 7.446 by using extrapolated values of , , , and from inside the domain. Then, for the lone incoming wave, the Linear Relaxation Method (LRM) of Poinsot [ 61, 62] is used to determine the value of the wave amplitude. The LRM method sets the value of the incoming wave amplitude to be proportional to the differences between the local pressure on a boundary face and the imposed exit pressure. Therefore, is given by
where is the imposed pressure at the exit boundary, is the relaxation factor, and is the local pressure value at the boundary.
In general, the desirable average pressure on a nonreflecting boundary can be either relaxed toward a pressure value at infinity or enforced to be equivalent to some desired pressure at the exit of the boundary.
If you want the average pressure at the boundary to relax toward at infinity (i.e. ), the suggested factor is given by:
where is the acoustic speed , is the domain size, is the maximum Mach number in the domain, and is the underrelaxation factor (default value is 0.15) . On the other hand, if the desired average pressure at the boundary is to approach a specific imposed value at the boundary, then the factor is given by:
where the default value for is 5.0
Using General NonReflecting Boundary Conditions
The general NRBC is available for use in the Pressure Outlet dialog box when either the densitybased explicit or the densitybased implicit solvers are activated to solve for compressible flows using the idealgas law.
To activate the general NRBC
define boundaryconditions nonreflectingbc generalnrbc set
In the set/ submenu, you can set the sigma value. The default value for sigma is 0.15.
define boundaryconditions nonreflectingbc generalnrbc set
In the set/ submenu, you can set the sigma2 value. The default value for sigma2 is 5.0 .

There is no guarantee that the
sigma2 value of 5.0 will force the average boundary pressure to match the specified exit pressure in all flow situations. In the case where the desired average boundary pressure has not been achieved, the user can intervene to adjust the
sigma2 value so that the desired average pressure on the boundary is approached.

Usually, the solver can operate at higher CFL values without the NRBCs being turned on. Therefore, for steadystate solutions the best practice is to first achieve a good stable solution (not necessarily converged) before activating the nonreflecting boundary condition. In many flow situations, the CFL value must be reduced from the normal operation to keep the solution stable. This is particularly true with the densitybased implicit solver since the boundary update is done in an explicit manner. A typical CFL value in the densitybased implicit solver, with the NRBC activated, is 2.0.