转载自: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/

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

 

U

p

T

inlet

pressureInletVelocity

totalPressure

inletOutlet

outlet

inletOutlet

totalPressure

inletOutlet

wall

fixedValue

zeroGradient

zeroGradient

k

Epsilon

mut

alpha

inletOutlet

inletOutlet

calculated

calculated

inletOutlet

inletOutlet

calculated

calculated

wallFunction

wallFunction

wallFunction

wallFunction

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

下面是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

OpenFOAM 中边界条件的设定【转载】的更多相关文章

  1. OpenFOAM 中的边界条件(一)【转载】

    链接:http://xiaopingqiu.github.io/2016/04/02/Boundary-conditions-in-OpenFOAM1/ 本系列解读 OpenFOAM 中边界条件的实现 ...

  2. OpenFOAM 中的边界条件(二)【转载】

    转载链接:http://xiaopingqiu.github.io/2016/04/02/Boundary-conditions-in-OpenFOAM2/ 本篇在上一篇的基础上来解读 OpenFOA ...

  3. 根据数据文件自定义边界条件timeVaryingUniformFixedValue【转载】

    转载自:http://blog.sina.com.cn/s/blog_e256415d0101nf9j.html 在OpenFOAM中,可以创建数据文件,自定义边界条件. 下面的例子读取outletP ...

  4. 深度学习---1cycle策略:实践中的学习率设定应该是先增再降

    深度学习---1cycle策略:实践中的学习率设定应该是先增再降 本文转载自机器之心Pro,以作为该段时间的学习记录 深度模型中的学习率及其相关参数是最重要也是最难控制的超参数,本文将介绍 Lesli ...

  5. spring中context:property-placeholder/元素 转载

    spring中context:property-placeholder/元素  转载 1.有些参数在某些阶段中是常量 比如 :a.在开发阶段我们连接数据库时的连接url,username,passwo ...

  6. SqlServer中的merge操作(转载)

    SqlServer中的merge操作(转载)   今天在一个存储过程中看见了merge这个关键字,第一个想法是,这个是配置管理中的概念吗,把相邻两次的更改合并到一起.后来在technet上搜索发现别有 ...

  7. OpenFOAM 中 c++ 基础

    文件布置 在 OpenFOAM 中,所有代码都以注释段开头,使用有限体积的 CFD 类型文件都包括以下头文件 #include "fvCFD.H" 在此头文件种,仅包含类或函数的定 ...

  8. 在OpenFOAM中做用户自定义库——编译library【转载】

    转载自:http://openfoam.blog.sohu.com/22041538.html OpenFOAM自己提供的标准类都是以库的形式提供的,并且利用头文件给出了库的应用接口.这样一来,用户的 ...

  9. OpenFOAM中的基本变量快速认知【转载】

    转载自:http://blog.sina.com.cn/s/blog_a0b4201d0102vsf9.html label 实际上就是整型数据的变体,int,OF对它进行了包装,以适应32或64位系 ...

随机推荐

  1. kvm第四章-- 虚拟化网络管理

  2. VBA Exit For语句

    当想要根据特定标准退出For循环时,就可以使用Exit For语句.当执行Exit For时,控件会立即跳转到For循环之后的下一个语句. 语法 以下是在VBA中Exit For语句的语法. Exit ...

  3. 关于INT_MIN

    来自为知笔记(Wiz)

  4. Dubbo面试

    DUBBO原理.应用与面经总结 SPI 你是否了解SPI,讲一讲什么是SPI,为什么要使用SPI? SPI具体约定:当服务的提供者,提供了服务接口的一种实现之后,在jar包的META-INF/serv ...

  5. C语言面试题目之指针和数组

    说明:所有题目均摘录于网络以及我所见过的面试题目,欢迎补充! 无特殊说明情况下,下面所有题s目都是linux下的32位C程序. 先来几个简单的热热身. 1.计算以下sizeof的值. char str ...

  6. Python字典取键、值对

    1. 取键:keys()方法 #spyder bb={'人才/可怕':23,'伏地魔&波特':'army','哈哈哈,人才,回合':'hhh'} for ii in bb.keys(): pr ...

  7. 第一周助教小结——发布作业&线上答疑

    第一周助教小结 助教博客:https://www.cnblogs.com/jason5689/ 本周点评数目:0份 由于发布的作业还未截至,第一次的作业点评还没开始进行,就描述一下评论博客前的感受吧 ...

  8. 运输层3——传输控制协议TCP概述

    目录 1. TCP最主要的特点 2. TCP的连接 3. socket在不同场景中的含义 写在前面:本文章是针对<计算机网络第七版>的学习笔记 运输层1--运输层协议概述 运输层2--用户 ...

  9. python easygui用法

    一个简单的例子 在 EasyGui 中,所有的 GUI 互动均是通过简单的函数调用,下边一个简单的例子告诉你 EasyGui 确实很 Easy! import easygui as g import ...

  10. js 面向对象之构造器与工厂函数

    字面量模式声明一个对象实例 let m = {name: "lisi", age:15} m.constructor // ƒ Object() { [native code] } ...