转载自：http://blog.sina.com.cn/s/blog_a0b4201d0102v7jt.html

用习惯了FLUENT的操作界面，再使用OpenFOAM就会觉得非常繁琐。遇到的第一个问题就是边界条件的设定。在FLUENT中，例如给定了一个进口边界（质量或者压力或者速度）和一个出口边界（一般是压力），并在进口处设定湍流强度就可以了，不必关心这三个量在求解的各个方程中的分配。而在OpenFOAM中，这个问题就要显得复杂的多，以k-epsl模型求解可压缩流为例为例，要求解质量守恒方程，动量守恒方程和能量守恒方程，湍动能方程以及耗散率方程，分别对应P方程，U方程，T方程，k方程和epsl方程，其中U方程是矢量方程，包含三个标量方程。在求解之前要为每个方程指定边界条件。现在的问题就是：只知道压力进口条件和压力出口回流条件，以及进出口都是室温条件，也就是说只能显式的给定P方程的进口条件和出口回流条件，T方程的进口和出口条件，而U方程的边界条件需要由已知的P方程和T方程的边界条件间接的得到。对于这个问题，OpenFOAM提供了各种各样的边界条件给定方式，多的足够让你眼花缭乱，所以非常有必要深入的扒一扒。先从OpenFOAM在tutorial中给出的例子入手。

以rhoSimplecFoam中的squareBend为例，给出了alphat,epsilon, k, mut, p, T, U 等7个量的初始条件，当然也包括了边界条件。其中alphat指的是turbulence thermal

diffusivity， 例子中给定的壁面边界是

compressible::alphatWallFunction; 进口边界是calculated, value uniform 0;

出口边界是 calculated, value uniform 0； epsilon指的是turbulence

kinetic energy dissipation

rate， k指的是turbulence

kinetic energy，mut指的是turbulence

viscosity，同时alphat, epsilon, k, mut

这个几个参数都是要给定壁面条件模型的。像其他的tutorial例子中给出的也都是速度进口边界条件，跟我的算例所需要的压力进口边界条件不太吻合。

好在OpenFoam给出了各种不同类型的边界条件类型供大家挑选使用，（当然是在不行就只能自己写了，这才是OF最强大的的地方嘛）。具体的不同内容可以查询页面： http://openfoam.org/docs/cpp/

最终我选择的边界条件是：

历经千辛万苦终于能够无报错运行啦

下面是openFoam给出的根据各种进出口条件分类的边界条件：

inlet boundary

conditions:

turbulentInletFvPatchField< Type >: This boundary condition generates a fluctuating inlet condition by adding a random component to a reference (mean) field.

surfaceNormalFixedValueFvPatchVectorField：This boundary condition provides a surface-normal vector boundary condition by its magnitude

outletMappedUniformInletFvPatchField< Type >：This boundary conditon averages the field over the "outlet" patch specified by name "outletPatchName" and applies this as the uniform value of the field over this patch

outletInletFvPatchField< Type >：This boundary condition provides a generic inflow condition, with specified outflow for the case of return flow

freestreamFvPatchField< Type >：This boundary condition provides a free-stream condition. It is a 'mixed' condition derived from the inletOutlet condition, whereby the mode of operation switches between fixed (free stream) value and zero gradient based on the sign of the flux

压力进口条件：

1: fixedPressureCompressibleDensityFvPatchScalarField: This boundary condition calculates a (liquid) compressible density as a function of pressure and fluid properties （这是属于密度进口条件）

2: uniformTotalPressureFvPatchScalarField: This boundary condition provides a time-varying form of the uniform total pressure boundary condition

3: totalPressureFvPatchScalarField: This boundary condition provides a total pressure condition. Four variants are possible:

4：freestreamPressureFvPatchScalarField：This boundary condition provides a free-stream condition for pressure. It is a zero-gradient condition that constrains the flux across the patch based on the free-stream velocity.

速度进口条件：

1: flowRateInletVelocityFvPatchVectorField: This boundary condition provides a velocity boundary condition, derived from the flux (volumetric or mass-based), whose direction is assumed to be normal to the patch.

2：atmBoundaryLayerInletVelocityFvPatchVectorField：This boundary condition specifies a velocity inlet profile appropriate for atmospheric boundary layers (ABL). The profile is derived from the friction velocity, flow direction and the direction of the parabolic co-ordinate z.

3: pressureNormalInletOutletVelocityFvPatchVectorField: This velocity inlet/outlet boundary condition is applied to patches where the pressure is specified. A zero-gradient condition is applied for outflow (as defined by the flux); for inflow, the velocity is obtained from the flux with a direction normal to the patch faces.

4：pressureInletVelocityFvPatchVectorField： This velocity inlet boundary condition is applied to patches where the pressure is specified. The inflow velocity is obtained from the flux with a direction normal to the patch faces

5： pressureInletUniformVelocityFvPatchVectorField： This velocity inlet boundary condition is applied to patches where the pressure is specified. The uniform inflow velocity is obtained by averaging the flux over the patch, and then applying it in the direction normal to the patch faces

6：pressureInletOutletParSlipVelocityFvPatchVectorField： This velocity inlet/outlet boundary condition for pressure boundary where the pressure is specified. A zero-gradient is applied for outflow (as defined by the flux); for inflow, the velocity is obtained from the flux with the specified inlet direction.

7： pressureDirectedInletVelocityFvPatchVectorField： This velocity inlet boundary condition is applied to patches where the pressure is specified. The inflow velocity is obtained from the flux with the specified inlet direction" direction.

8：pressureDirectedInletOutletVelocityFvPatchVectorField： This velocity inlet/outlet boundary condition is applied to pressure boundaries where the pressure is specified. A zero-gradient condtion is applied for outflow (as defined by the flux); for inflow, the velocity is obtained from the flux with the specified inlet direction.

9： flowRateInletVelocityFvPatchVectorField： This boundary condition provides a velocity boundary condition, derived from the flux (volumetric or mass-based), whose direction is assumed to be normal to the patch.

温度进口条件：

1： totalTemperatureFvPatchScalarField： This boundary condition provides a total temperature condition.

耗散率进口条件：epsilon进口条件：

1: turbulentMixingLengthDissipationRateInletFvPatchScalarField：This boundary condition provides a turbulence dissipation, $\epsilon$ (epsilon) inlet condition based on a specified mixing length

2: atmBoundaryLayerInletEpsilonFvPatchScalarField: This boundary condition specifies an inlet value for the turbulence dissipation, $\epsilon$ (epsilon), appropriate for atmospheric boundary layers (ABL), and designed to be used in conjunction with the ABLInletVelocity inlet velocity boundary condition.

3:

湍动能进口条件：

1: turbulentIntensityKineticEnergyInletFvPatchScalarField: This boundary condition provides a turbulent kinetic energy condition, based on user-supplied turbulence intensity, defined as a fraction of the mean velocity:

outlet boundary condition出口边界条件：

freestreamFvPatchField< Type >：This boundary condition provides a free-stream condition. It is a 'mixed' condition derived from the inletOutlet condition, whereby the mode of operation switches between fixed (free stream) value and zero gradient based on the sign of the flux.

inletOutletFvPatchField< Type >：This boundary condition provides a generic outflow condition, with specified inflow for the case of return flow.

uniformInletOutletFvPatchField< Type >：Variant of inletOutlet boundary condition with uniform inletValue.

速度出口条件：

fixedNormalInletOutletVelocityFvPatchVectorField: This velocity inlet/outlet boundary condition combines a fixed normal component obtained from the "normalVelocity" patchField supplied with a fixed or zero-gradiented tangential component depending on the direction of the flow and the setting of "fixTangentialInflow"

fluxCorrectedVelocityFvPatchVectorField： This boundary condition provides a velocity outlet boundary condition for patches where the pressure is specified. The outflow velocity is obtained by "zeroGradient" and then corrected from the flux:

pressureDirectedInletOutletVelocityFvPatchVectorField：This velocity inlet/outlet boundary condition is applied to pressure boundaries where the pressure is specified. A zero-gradient condtion is applied for outflow (as defined by the flux); for inflow, the velocity is obtained from the flux with the specified inlet direction.

pressureInletOutletVelocityFvPatchVectorField：This velocity inlet/outlet boundary condition is applied to pressure boundaries where the pressure is specified. A zero-gradient condition is applied for outflow (as defined by the flux); for inflow, the velocity is obtained from the patch-face normal component of the internal-cell value

pressureNormalInletOutletVelocityFvPatchVectorField：This velocity inlet/outlet boundary condition is applied to patches where the pressure is specified. A zero-gradient condition is applied for outflow (as defined by the flux); for inflow, the velocity is obtained from the flux with a direction normal to the patch faces.

压力出口条件：

uniformTotalPressureFvPatchScalarField：This boundary condition provides a time-varying form of the uniform total pressure boundary condition.

totalPressureFvPatchScalarField：This boundary condition provides a total pressure condition. Four variants are possible

温度出口条件：

inletOutletTotalTemperatureFvPatchScalarField：This boundary condition provides an outflow condition for total temperature for use with supersonic cases, where a user-specified value is applied in the case of reverse flow

totalTemperatureFvPatchScalarField:This boundary condition provides a total temperature condition.

一般性的边界条件： generic boundary conditions：http://foam.sourceforge.net/docs/cpp/a10590.html

calculatedFvPatchField< Type >：This boundary condition is not designed to be evaluated; it is assmued that the value is assigned via field assignment, and not via a call to e.g. updateCoeffs or evaluate.

directionMixedFvPatchField< Type >：Base class for direction-mixed boundary conditions

fixedGradientFvPatchField< Type >：This boundary condition supplies a fixed gradient condition, such that the patch values are calculated using

fixedValueFvPatchField< Type >： This boundary condition supplies a fixed value constraint, and is the base class for a number of other boundary conditions.

zeroGradientFvPatchField< Type >：This boundary condition applies a zero-gradient condition from the patch internal field onto the patch faces.

oscillatingFixedValueFvPatchField< Type >：This boundary condition provides an oscillating condition in terms of amplitude and frequency.

uniformFixedGradientFvPatchField< Type >：This boundary condition provides a uniform fixed gradient condition.

uniformFixedValueFvPatchField< Type >：This boundary condition provides a uniform fixed value condition