OpenFOAM编程指南

Op
O p enFOA M编程指
指南南
OpenFOAM User Guide Chinese Edition
Christopher J. Greenshields 著
李东岳译
2
Copyright©2011-2014 OpenFOAM Foundation.
This work is licensed under a Creative Commons Attribution - NonCommercial

NoDerivs 3.0 Unported License.
版权声明©2011-2014 OpenFOAM Foundation
本指南遵守：Creative Commons Attribution - NonCommercial NoDerivs 3.0

Unported License.
License
THE WORK (AS DEFINED BELOW) IS PROVIDED UNDER THE TERMS OF THIS
CREATIVE COMMONS PUBLIC LICENSE ("CCPL" OR "LICENSE"). THE WORK IS
PROTECTED BY COPYRIGHT AND/OR OTHER APPLICABLE LAW. ANY USE OF THE
WORK OTHER THAN AS AUTHORIZED UNDER THIS LICENSE OR COPYRIGHT LAW IS
PROHIBITED.
BY EXERCISING ANY RIGHTS TO THE WORK PROVIDED HERE, YOU ACCEPT AND
AGREE TO BE BOUND BY THE TERMS OF THIS LICENSE. TO THE EXTENT THIS
LICENSE MAY BE CONSIDERED TO BE A CONTRACT, THE LICENSOR GRANTS YOU
THE RIGHTS CONTAINED HERE IN CONSIDERATION OF YOUR ACCEPTANCE OF
SUCH TERMS AND CONDITIONS.
1. Definitions
a. "Adaptation" means a work based upon the Work, or upon the Work and other pre-
existing works, such as a translation, adaptation, derivative work, arrangement of
music or other alterations of a literary or artistic work, or phonogram or performance
and includes cinematographic adaptations or any other form in which the Work may
be recast, transformed, or adapted including in any form recognizably derived from
the original, except that a work that constitutes a Collection will not be considered an
Adaptation for the purpose of this License. For the avoidance of doubt, where the
Work is a musical work, performance or phonogram, the synchronization of the Work
in timed-relation with a moving image ("synching") will be considered an Adaptation
for the purpose of this License.
b. "Collection" means a collection of literary or artistic works, such as encyclopedias

and anthologies, or performances, phonograms or broadcasts, or other works or
subject matter other than works listed in Section 1(f) below, which, by reason of the
3
selection and arrangement of their contents, constitute intellectual creations, in which

the Work is included in its entirety in unmodified form along with one or more other
contributions, each constituting separate and independent works in themselves,
which together are assembled into a collective whole. A work that constitutes a
Collection will not be considered an Adaptation (as defined above) for the purposes of
this License.
c. "Distribute" means to make available to the public the original and copies of the
Work through sale or other transfer of ownership.
d. "Licensor" means the individual, individuals, entity or entities that offer(s) the Work
under the terms of this License.
e. "Original Author" means, in the case of a literary or artistic work, the individual,
individuals, entity or entities who created the Work or if no individual or entity can be
identified, the publisher; and in addition (i) in the case of a performance the actors,
singers, musicians, dancers, and other persons who act, sing, deliver, declaim, play
in, interpret or otherwise perform literary or artistic works or expressions of folklore;
(ii) in the case of a phonogram the producer being the person or legal entity who first
fixes the sounds of a performance or other sounds; and, (iii) in the case of
broadcasts, the organization that transmits the broadcast.
f. "Work" means the literary and/or artistic work offered under the terms of this License
including without limitation any production in the literary, scientific and artistic domain,
whatever may be the mode or form of its expression including digital form, such as a
book, pamphlet and other writing; a lecture, address, sermon or other work of the
same nature; a dramatic or dramatico-musical work; a choreographic work or
entertainment in dumb show; a musical composition with or without words; a
cinematographic work to which are assimilated works expressed by a process
analogous to cinematography; a work of drawing, painting, architecture, sculpture,
engraving or lithography; a photographic work to which are assimilated works
expressed by a process analogous to photography; a work of applied art; an
illustration, map, plan, sketch or three-dimensional work relative to geography,
topography, architecture or science; a performance; a broadcast; a phonogram; a
compilation of data to the extent it is protected as a copyrightable work; or a work
performed by a variety or circus performer to the extent it is not otherwise considered
a literary or artistic work.
g. "You" means an individual or entity exercising rights under this License who has not
previously violated the terms of this License with respect to the Work, or who has
received express permission from the Licensor to exercise rights under this License
despite a previous violation.
h. "Publicly Perform" means to perform public recitations of the Work and to

communicate to the public those public recitations, by any means or process,
including by wire or wireless means or public digital performances; to make available
to the public Works in such a way that members of the public may access these
Works from a place and at a place individually chosen by them; to perform the Work
4
to the public by any means or process and the communication to the public of the
performances of the Work, including by public digital performance; to broadcast and
rebroadcast the Work by any means including signs, sounds or images.
i. "Reproduce" means to make copies of the Work by any means including without
limitation by sound or visual recordings and the right of fixation and reproducing
fixations of the Work, including storage of a protected performance or phonogram in
digital form or other electronic medium.
2. Fair Dealing Rights. Nothing in this License is intended to reduce, limit, or restrict any
uses free from copyright or rights arising from limitations or exceptions that are provided for in
connection with the copyright protection under copyright law or other applicable laws.
3. License Grant. Subject to the terms and conditions of this License, Licensor hereby grants
You a worldwide, royalty-free, non-exclusive, perpetual (for the duration of the applicable
copyright) license to exercise the rights in the Work as stated below:
a. to Reproduce the Work, to incorporate the Work into one or more Collections, and to
Reproduce the Work as incorporated in the Collections; and,
b. to Distribute and Publicly Perform the Work including as incorporated in Collections.
The above rights may be exercised in all media and formats whether now known or hereafter
devised. The above rights include the right to make such modifications as are technically
necessary to exercise the rights in other media and formats, but otherwise you have no rights
to make Adaptations. Subject to 8(f), all rights not expressly granted by Licensor are hereby
reserved, including but not limited to the rights set forth in Section 4(d).
4. Restrictions. The license granted in Section 3 above is expressly made subject to and
limited by the following restrictions:
a. You may Distribute or Publicly Perform the Work only under the terms of this License.
You must include a copy of, or the Uniform Resource Identifier (URI) for, this License
with every copy of the Work You Distribute or Publicly Perform. You may not offer or
impose any terms on the Work that restrict the terms of this License or the ability of
the recipient of the Work to exercise the rights granted to that recipient under the
terms of the License. You may not sublicense the Work. You must keep intact all
notices that refer to this License and to the disclaimer of warranties with every copy of
the Work You Distribute or Publicly Perform. When You Distribute or Publicly Perform
the Work, You may not impose any effective technological measures on the Work
that restrict the ability of a recipient of the Work from You to exercise the rights
5
granted to that recipient under the terms of the License. This Section 4(a) applies to
the Work as incorporated in a Collection, but this does not require the Collection apart
from the Work itself to be made subject to the terms of this License. If You create a
Collection, upon notice from any Licensor You must, to the extent practicable, remove
from the Collection any credit as required by Section 4(c), as requested.
b. You may not exercise any of the rights granted to You in Section 3 above in any
manner that is primarily intended for or directed toward commercial advantage or
private monetary compensation. The exchange of the Work for other copyrighted
works by means of digital file-sharing or otherwise shall not be considered to be
intended for or directed toward commercial advantage or private monetary
compensation, provided there is no payment of any monetary compensation in
connection with the exchange of copyrighted works.
c. If You Distribute, or Publicly Perform the Work or Collections, You must, unless a
request has been made pursuant to Section 4(a), keep intact all copyright notices for
the Work and provide, reasonable to the medium or means You are utilizing: (i) the
name of the Original Author (or pseudonym, if applicable) if supplied, and/or if the
Original Author and/or Licensor designate another party or parties (e.g., a sponsor
institute, publishing entity, journal) for attribution ("Attribution Parties") in Licensor's
copyright notice, terms of service or by other reasonable means, the name of such
party or parties; (ii) the title of the Work if supplied; (iii) to the extent reasonably
practicable, the URI, if any, that Licensor specifies to be associated with the Work,
unless such URI does not refer to the copyright notice or licensing information for the
Work. The credit required by this Section 4(c) may be implemented in any reasonable
manner; provided, however, that in the case of a Collection, at a minimum such credit
will appear, if a credit for all contributing authors of Collection appears, then as part of
these credits and in a manner at least as prominent as the credits for the other
contributing authors. For the avoidance of doubt, You may only use the credit
required by this Section for the purpose of attribution in the manner set out above
and, by exercising Your rights under this License, You may not implicitly or explicitly
assert or imply any connection with, sponsorship or endorsement by the Original
Author, Licensor and/or Attribution Parties, as appropriate, of You or Your use of the
Work, without the separate, express prior written permission of the Original Author,
Licensor and/or Attribution Parties.
d. For the avoidance of doubt:

i. Non-waivable Compulsory License Schemes. In those jurisdictions in
which the right to collect royalties through any statutory or compulsory
licensing scheme cannot be waived, the Licensor reserves the exclusive right
to collect such royalties for any exercise by You of the rights granted under
this License;
ii. Waivable Compulsory License Schemes. In those jurisdictions in which the

right to collect royalties through any statutory or compulsory licensing
scheme can be waived, the Licensor reserves the exclusive right to collect
such royalties for any exercise by You of the rights granted under this
6
License if Your exercise of such rights is for a purpose or use which is

otherwise than noncommercial as permitted under Section 4(b) and
otherwise waives the right to collect royalties through any statutory or
compulsory licensing scheme; and,
iii. Voluntary License Schemes. The Licensor reserves the right to collect
royalties, whether individually or, in the event that the Licensor is a member
of a collecting society that administers voluntary licensing schemes, via that
society, from any exercise by You of the rights granted under this License
that is for a purpose or use which is otherwise than noncommercial as
permitted under Section 4(b).
e. Except as otherwise agreed in writing by the Licensor or as may be otherwise

permitted by applicable law, if You Reproduce, Distribute or Publicly Perform the
Work either by itself or as part of any Collections, You must not distort, mutilate,
modify or take other derogatory action in relation to the Work which would be
prejudicial to the Original Author's honor or reputation.
5. Representations, Warranties and Disclaimer
UNLESS OTHERWISE MUTUALLY AGREED BY THE PARTIES IN WRITING, LICENSOR

OFFERS THE WORK AS-IS AND MAKES NO REPRESENTATIONS OR WARRANTIES OF
ANY KIND CONCERNING THE WORK, EXPRESS, IMPLIED, STATUTORY OR
OTHERWISE, INCLUDING, WITHOUT LIMITATION, WARRANTIES OF TITLE,
MERCHANTIBILITY, FITNESS FOR A PARTICULAR PURPOSE, NONINFRINGEMENT, OR
THE ABSENCE OF LATENT OR OTHER DEFECTS, ACCURACY, OR THE PRESENCE OF
ABSENCE OF ERRORS, WHETHER OR NOT DISCOVERABLE. SOME JURISDICTIONS
DO NOT ALLOW THE EXCLUSION OF IMPLIED WARRANTIES, SO SUCH EXCLUSION
MAY NOT APPLY TO YOU.
6. Limitation on Liability. EXCEPT TO THE EXTENT REQUIRED BY APPLICABLE LAW,

IN NO EVENT WILL LICENSOR BE LIABLE TO YOU ON ANY LEGAL THEORY FOR ANY
SPECIAL, INCIDENTAL, CONSEQUENTIAL, PUNITIVE OR EXEMPLARY DAMAGES
ARISING OUT OF THIS LICENSE OR THE USE OF THE WORK, EVEN IF LICENSOR HAS
BEEN ADVISED OF THE POSSIBILITY OF SUCH DAMAGES.
7. Termination
a. This License and the rights granted hereunder will terminate automatically upon any
breach by You of the terms of this License. Individuals or entities who have received
Collections from You under this License, however, will not have their licenses
7
terminated provided such individuals or entities remain in full compliance with those
licenses. Sections 1, 2, 5, 6, 7, and 8 will survive any termination of this License.
b. Subject to the above terms and conditions, the license granted here is perpetual (for
the duration of the applicable copyright in the Work). Notwithstanding the above,
Licensor reserves the right to release the Work under different license terms or to
stop distributing the Work at any time; provided, however that any such election will
not serve to withdraw this License (or any other license that has been, or is required
to be, granted under the terms of this License), and this License will continue in full
force and effect unless terminated as stated above.
8. Miscellaneous
a. Each time You Distribute or Publicly Perform the Work or a Collection, the Licensor
offers to the recipient a license to the Work on the same terms and conditions as the
license granted to You under this License.
b. If any provision of this License is invalid or unenforceable under applicable law, it

shall not affect the validity or enforceability of the remainder of the terms of this
License, and without further action by the parties to this agreement, such provision
shall be reformed to the minimum extent necessary to make such provision valid and
enforceable.
c. No term or provision of this License shall be deemed waived and no breach

consented to unless such waiver or consent shall be in writing and signed by the
party to be charged with such waiver or consent.
d. This License constitutes the entire agreement between the parties with respect to the
Work licensed here. There are no understandings, agreements or representations
with respect to the Work not specified here. Licensor shall not be bound by any
additional provisions that may appear in any communication from You. This License
may not be modified without the mutual written agreement of the Licensor and You.
e. The rights granted under, and the subject matter referenced, in this License were
drafted utilizing the terminology of the Berne Convention for the Protection of Literary
and Artistic Works (as amended on September 28, 1979), the Rome Convention of
1961, the WIPO Copyright Treaty of 1996, the WIPO Performances and Phonograms
Treaty of 1996 and the Universal Copyright Convention (as revised on July 24, 1971).
These rights and subject matter take effect in the relevant jurisdiction in which the
License terms are sought to be enforced according to the corresponding provisions of
the implementation of those treaty provisions in the applicable national law. If the
standard suite of rights granted under applicable copyright law includes additional
rights not granted under this License, such additional rights are deemed to be
included in the License; this License is not intended to restrict the license of any
rights under applicable law.
8
This OpenFOAM Programmer’s Guide, Chinese Edition is translated into Chinese by Dongyue
Li1(Chapter 1, 2), Xiangyu Gao2 (Chaper 3).
1
DYFLUID Studio, Torino, Piedmont, Italy, 10143. 2 Purdue University, West Lafayette, IN, USA, 47907
《OpenFOAM 编程指南中文版》由李东岳 1 翻译第一、二章，高翔宇 2 翻译第三章

1
东岳流体, 意大利皮埃蒙特区都灵市, 10143. 2 普渡大学, 美国印第安纳州西拉法叶, 47907.
9
目录
License............................................................................................................................................ 2
第一章 ............................................................................................................................................ 11
1.1 坐标系统.......................................................................................................................... 11
1.2 张量.................................................................................................................................. 11
1.2.1 张量表示法 ........................................................................................................... 13
1.3 张量算术.......................................................................................................................... 13
1.3.1 内积....................................................................................................................... 14
1.3.2 双内积................................................................................................................... 15
1.3.3 两个三阶张量的三内积 ....................................................................................... 15
1.3.4 外积....................................................................................................................... 15
1.3.5 矢量的叉乘 ........................................................................................................... 15
1.3.6 其他张量操作 ....................................................................................................... 16
1.3.7 几何变形和单位张量 ........................................................................................... 16
1.3.8 张量等式............................................................................................................... 17
1.3.9 二阶张量算术 ....................................................................................................... 17
1.3.10 标量算术............................................................................................................. 18
1.4 OpenFOAM 张量类 ......................................................................................................... 19
1.4.1 OpenFOAM 张量操作........................................................................................... 19
1.5 量纲单位.......................................................................................................................... 21
第二章 ............................................................................................................................................ 23
2.1 微分算符.......................................................................................................................... 23
2.1.1 梯度....................................................................................................................... 23
2.1.2 散度....................................................................................................................... 24
2.1.3 旋度....................................................................................................................... 24
2.1.4 拉普拉斯............................................................................................................... 24
2.1.5 时间导数............................................................................................................... 24
2.2 离散简述.......................................................................................................................... 25
2.2.1 OpenFOAM 链表和场........................................................................................... 25
2.3 计算域离散...................................................................................................................... 25
2.3.1 定义 OpenFOAM 网格......................................................................................... 27
2.3.2 定义 geometricField.............................................................................................. 28
2.4 方程离散.......................................................................................................................... 29
2.4.1 拉普拉斯项 ........................................................................................................... 34
2.4.2 对流项................................................................................................................... 34
2.4.3 一阶时间导数 ....................................................................................................... 35
2.4.4 二阶时间导数项 ................................................................................................... 35
2.4.5 散度项................................................................................................................... 35
2.4.6 梯度项................................................................................................................... 36
2.4.7 梯度-梯度的平方 ................................................................................................. 37
2.4.8 旋度....................................................................................................................... 37
2.4.9 源项....................................................................................................................... 37
2.4.10 其他显性离散 ..................................................................................................... 37
10
2.5 时间离散.......................................................................................................................... 38
2.5.1 OpenFOAM 中时间离散的代码........................................................................... 39
2.6 边界条件.......................................................................................................................... 39
2.6.1 物理边界条件 ....................................................................................................... 40
第三章 ............................................................................................................................................ 41
3.1 圆柱绕流.......................................................................................................................... 41
3.1.1 问题阐述............................................................................................................... 42
3.1.2 potentialFoam 求解器的注意事项 ........................................................................ 43
3.1.3 网格的生成 ........................................................................................................... 43
3.1.4 边界条件和初始条件 ........................................................................................... 45
3.1.5 运行算例............................................................................................................... 46
3.2 稳态后向台阶湍流模拟 .................................................................................................. 49
3.2.1 问题阐述............................................................................................................... 49
3.2.2 网格的生成 ........................................................................................................... 50
3.2.3 边界条件和初始场 ............................................................................................... 53
3.2.4 求解器的控制参数 ............................................................................................... 53
3.2.5 运行算例和后处理 ............................................................................................... 54
3.3 前向台阶的超音速绕流 .................................................................................................. 54
3.3.1 问题阐述............................................................................................................... 54
3.3.2 网格生成............................................................................................................... 56
3.3.3 运行算例............................................................................................................... 57
3.3.4 练习....................................................................................................................... 57
3.4 高压水箱放水问题 .......................................................................................................... 58
3.4.1 问题阐述............................................................................................................... 58
3.4.2 网格的生成 ........................................................................................................... 59
3.4.3 准备运行算例 ....................................................................................................... 61
3.4.4 运行算例............................................................................................................... 62
3.4.5 细化网格来获得更精确解 ................................................................................... 63
3.5 磁流体问题...................................................................................................................... 63
3.5.1 问题阐述............................................................................................................... 64
3.5.2 网格的生成 ........................................................................................................... 65
3.5.3 运行算例............................................................................................................... 66
11
第一章
张量
这一章我们介绍张量和张量算术，以及他们在本手册中是如何进行表达的。然后我们解释张
量和张量算术如何在 OpenFOAM 中的实现。
1.1 坐标系统
OpenFOAM 主要用来解决连续介质问题，例如：固体应力分析、液体气体流动问题以及物
理材料的变形问题。因此，OpenFOAM 在三维空间和时间中处理用张量来进行描述的物理
问题。参见图 1.1，OpenFOAM 使用的坐标系统符合右手规则。建立于右手规则的坐标系符
合下述特征：三个坐标轴起源于 O 点，他们的名字为 Ox、Oy、Oz，当从 Oz 向上观测的时
候（O 点离观测者较近），从 Ox 上某点指向 Oy 上某点的圆弧为顺时针方向1。
图 1.1 右手规则坐标轴
1.2 张量
张量用来描述从属某一空间的实体，其遵循一定的数学法则。简要来说，张量就是和一系列
单位矢量相关的一系列的值；在 OpenFOAM 中，单位矢量𝒊𝒙 , 𝒊𝒚 , 𝒊𝒛和上文提及的符合右手规
1 本手册相对于 OpenFOAM ProgrammersGuide 对行距有所调整以保证页数匹配

12
则的坐标系相关。单位矢量是正交的，每个张量具有如下特征：
 维度 d：我们用 d 来表示，在 OpenFOAM 中 d = 3；
 阶（秩）r：其为一个整数，表示分量的个数，例如我们可以这样表示：𝑑 𝑟 ；
在 OpenFOAM 1.x 中，我们默认设置为 3 维（d = 3）
，并提供 0 到 3 阶张量作为标
准张量，对张量进行扩展也非常容易。我们熟知 0 阶张量和 1 阶张量代表标量和矢
量，2 阶张量和 3 阶张量有些读者可能会感到有些陌生。为了文章的完整性我们把
他们都进行了描述；
 0 阶标量：我们用 0 阶张量表示任何一个可以用单个数值表示的物理量，例如：质
量 m，体积 V，压力 p，粘度𝜇；
 1 阶矢量：它用来表示具有大小和方向的物理量。他们通常的形式为： a = (𝑎1 , 𝑎2 , 𝑎3 )，

其和笛卡尔坐标系的 x，y，z 相对应。有时我们用索引符 i 表示这个张量，例如
𝑎𝑖 , 𝑖 = 1,2,3，但是本书中我们将其省略，因为我们处理的问题均为 3 维的；
 2 阶张量：我们用 T 来表示，T 具有 9 个元素，参见如下：
𝑇11 𝑇12 𝑇13

T = 𝑇𝑖𝑗 = (𝑇21 𝑇22 𝑇23 ) (1.1)
𝑇31 𝑇32 𝑇33
由于其为 2 阶张量，因此分量𝑇𝑖𝑗 使用两个索引符来表示，如上述，我们在下文省略

ij 索引符号。当𝑖 = 𝑗的时候，𝑇𝑖𝑗 为对角元素。那些𝑖 ≠ 𝑗的元素为非对角元素。T 的
转置矩阵如下：
𝑇11 𝑇21 𝑇31

𝑇
T = 𝑇𝑗𝑖 = ( 21 𝑇22 𝑇32 ) (1.2)
𝑇31 𝑇23 𝑇33
需要注意的是，我们通常将二阶张量称为张量，因为再往上的高阶张量非常罕见；
 例如𝑇𝑖𝑗 = 𝑇𝑗𝑖 ，
对称 2 阶张量：对称张量指得是非对角元素关于对角元素是对称的。
在这种情况下，只有 6 个独立元素（ 𝑇12 = 𝑇21 ，𝑇13 = 𝑇31 ，𝑇23 = 𝑇32 ）。在
OpenFOAM 中我们区分对称张量和非对称张量来节省内存。因为对称张量只需要
存储 6 个元素，非对称张量需要存储 9 个元素。大多数连续介质问题中的张量为对
称张量；
 3 阶张量：具有 27 个元素，我们用𝑃𝑖𝑗𝑘 来表示，由于其太复杂在这里我们不表示；
 对称 3 阶张量：OpenFOAM 将其定义为下述元素相等的 3 阶张量：𝑃𝑖𝑗𝑘 = 𝑃𝑖𝑘𝑗 =

𝑃𝑗𝑖𝑘 = 𝑃𝑗𝑘𝑖 = 𝑃𝑘𝑖𝑗 = 𝑃𝑘𝑗𝑖 ，因此它具有 10 个独立分量。更确切来说，3 个相同的矢
量进行外积即可得到对称 3 阶张量，我们在 1.3.4 节中介绍外积；
13
1.2.1 张量表示法
本书是一本关于连续介质计算力学的书籍，主要处理 3 维空间和时间的复杂偏微分方程组。
在开始后续章节之前我们很有必要了解一下张量的表示方法以防止在后续章节发生混乱。为
了很好的理解方程组，我们希望读者接受这样一个概念：张量是一个整体而不是一系列标量
的组合。另外，任何张量操作都是在这个整体上进行操作而不是对某个元素进行操作。
在本书下面的章节中，所有的张量（除了标量）我们都是用粗体，a。我们使用单一字
符来表示张量，这更加的表明张量本身是一个实体，并且这样表示张量非常简洁。这种表达
方式潜在的缺点就是张量的阶数并没有有效的表达。很明显，张量的阶数不能为 0。实际上，
由于我们很清楚的了解各个张量代表的物理量的实际意义，缺少张量阶数的这种表达方式也
不会对我们造成太大的麻烦。例如我们知道速度 U 的阶数为 1。
另外一个关于张量操作符的基本概念是张量和坐标系统的无关性。例如矢量 a 不管我
们从哪个角度来观测它都是矢量 a。这意味着张量和坐标系统无关。然而其他符号，例如𝑎𝑖 ，
它则和坐标系统有关，因为它表示的是张量的某个分量。这样的结果就是张量由一系列和坐
标系统有关的元素值来表示，这些元素值依赖于坐标系统。
正如 1.2 节所说，对于张量的索引符表达方式将贯彻全书，它的主要目的是为了将张量
运算拓展到其中的每个分量上。当使用索引符表示张量的时候，我们定义，如果在一项中某
个索引符出现了两次，则这一项表示元素（例如第 1、2、3 元素）的相加。例如：
𝑎𝑖 𝑏𝑖 = ∑3𝑖=1 𝑎𝑖 𝑏𝑖 = 𝑎1 𝑏1 + 𝑎2 𝑏2 + 𝑎3 𝑏3 (1.3)
利用简化的上述定义，在下面的章节中，∑ 将被省略。
1.3 张量算术
在这一章节我们描述 OpenFOAM 可用的张量算术。我们首先来看一下最常见的张量操作：

加、减、乘、除。加法和减法只能同阶操作。他们的操作即为将张量对应的元素分别进行加
减法。例如，矢量 a 和矢量 b 的减法即为：
a − b = 𝑎𝑖 − 𝑏𝑖 = (𝑎1 − 𝑏1 , 𝑎2 − 𝑏2 , 𝑎3 − 𝑏3 ) (1.4)
张量 a 和标量 s 的乘法也是可交换的，即为将张量的每个元素都乘以标量，例如：
s𝐚 = 𝑠𝑎𝑖 = (𝑠𝑎1 , 𝑠𝑎2 , 𝑠𝑎3 ) (1.5)
在做张量 a 和标量 s 的除法的时候，标量 s 必须为分母，例如：
𝐚/𝑠 = 𝒂𝑖 /𝑠 = (𝑎1 /𝑠, 𝑎2 /𝑠, 𝑎3 /𝑠) (1.6)
依据上述基本定义，在接下来章节中，我们讨论一阶张量的更为复杂的乘积操作。
14
1.3.1 内积
张量的内积可以在阶数为𝑟1 的张量和阶数为𝑟2 的张量上进行操作。结果张量的阶数为𝑟 =

𝑟1 + 𝑟2 − 2，下面我们描述至高 3 阶的张量内积：
 矢量 a 和矢量 b 的内积是可交换的，并且结果为一个标量 s：
s = 𝑎𝑖 𝑏𝑖 = 𝑎1 𝑏1 + 𝑎2 𝑏2 + 𝑎3 𝑏3 (1.7)
 张量 T 和矢量 a 的内积结果为矢量 b，我们用列向量来表示𝐛 = 𝐓 ∙ 𝐚（看起来更加

方便）：
𝑇11 𝑎1 + 𝑇12 𝑎2 + 𝑇13 𝑎3

𝑏𝑖 = 𝑇𝑖𝑗 𝑎𝑗 = (𝑇21 𝑎1 + 𝑇22 𝑎2 + 𝑇23 𝑎3 ) (1.8)
𝑇31 𝑎1 + 𝑇32 𝑎2 + 𝑇33 𝑎3
如果 T 是非对称张量，那么这个内积操作是非对称的，即𝐛 = 𝐚 ∙ 𝐓 = 𝐓 𝐓 ∙ 𝐚：
𝑇11 𝑎1 + 𝑇21 𝑎2 + 𝑇31 𝑎3

𝑏𝑖 = 𝑎𝑗 𝑇𝑗𝑖 = (𝑇12 𝑎1 + 𝑇22 𝑎2 + 𝑇32 𝑎3 ) (1.9)
𝑇13 𝑎1 + 𝑇23 𝑎2 + 𝑇33 𝑎3
 两个张量 T 和 S 的内积也是一个张量，即 P = T ∙ S，P 元素的值这样计算：
𝑃𝑖𝑗 = 𝑇𝑖𝑘 𝑆𝑘𝑗 (1.10)
并且这种运算是不可交换的，例如：T ∙ S = (𝑆 𝑇 ∙ 𝑇 𝑇 )𝑇
 矢量 a 和三阶张量 P 的内积结果为一个二阶张量 T，T = a ∙ P，T 的元素为：
𝑇𝑖𝑗 = 𝑎𝑘 𝑃𝑘𝑖𝑗 (1.11)
同样，这个运算也是不可交换的，因此 T = P ∙ a 的元素为：
𝑇𝑖𝑗 = 𝑃𝑖𝑗𝑘 𝑎𝑘 (1.12)
 二阶张量 T 和三阶张量 P 的内积为三阶张量 Q，Q = T ∙ P，Q 的元素为：
𝑄𝑖𝑗𝑘 = 𝑇𝑖𝑙 𝑃𝑙𝑗𝑘 (1.13)
同样，这个运算也是不可交换的，因此 Q = P ∙ T 的元素为：
𝑄𝑖𝑗𝑘 = 𝑃𝑖𝑗𝑙 𝑇𝑙𝑘 (1.14)

15
1.3.2 双内积
两个二阶张量 T 和 S 的双内积结果为一个标量 s（s = T : S），s 这样计算：
𝑠 = 𝑇𝑖𝑗 𝑆𝑖𝑗 = 𝑇11 𝑆11 + 𝑇12 𝑆12 + 𝑇13 𝑆13 + (1.15)

𝑇21 𝑆21 + 𝑇22 𝑆22 + 𝑇23 𝑆23
𝑇31 𝑆31 + 𝑇32 𝑆32 + 𝑇33 𝑆33
二阶张量 T 和三阶张量 P 的双内积结果为一个矢量 a（a = T : P）:
𝑎𝑖 = 𝑇𝑗𝑘 𝑃𝑗𝑘𝑖 (1.16)
这个运算是不可交换的，因此 a = P : T 为：
𝑎𝑖 = 𝑃𝑖𝑗𝑘 𝑇𝑗𝑘 (1.17)
1.3.3 两个三阶张量的三内积
两个三阶张量 P 和 Q 的三内积结果为一个标量 s（s = P ⋮ Q），s 这样计算：
𝑠 = 𝑃𝑖𝑗𝑘 𝑄𝑖𝑗𝑘 (1.18)
1.3.4 外积
张量以及矢量的外积这样定义：
 两个矢量 a 和 b 的外积是不可交换的，结果为张量T = ab = (ba)T ，其这样计算：
𝑎1 𝑏1 𝑎1 𝑏2 𝑎1 𝑏3
𝑇𝑖𝑗 = 𝑎𝑖 𝑏𝑗 = (𝑎2 𝑏1 𝑎2 𝑏2 𝑎2 𝑏3 ) (1.19)
𝑎3 𝑏1 𝑎3 𝑏2 𝑎3 𝑏3
 一个矢量 a 和一个张量 T 的外积为一个三阶张量 P（P = a T），其这样计算：
𝑃𝑖𝑗𝑘 = 𝑎𝑖 𝑇𝑗𝑘 (1.20)
同样这个运算为不可交换的，因此 P = T a 的结果为：
𝑃𝑖𝑗𝑘 = 𝑇𝑗𝑘 𝑎𝑖 (1.21)
1.3.5 矢量的叉乘
叉乘只能用于矢量，对于两个矢量 a 和 b，他们的叉乘为矢量 c，其这样计算：

16
𝑐 = 𝑎 × 𝑏 = 𝑒𝑖𝑗𝑘 𝑎𝑗 𝑏𝑘 = (𝑎2𝑏3 − 𝑎3 𝑏2 , 𝑎3 𝑏1 − 𝑎1 𝑏3 , 𝑎1 𝑏2 − 𝑎2 𝑏1 ) (1.22)
其中𝑒𝑖𝑗𝑘 这样定义：
0 当两个索引符相同的时候
eijk = {+1 ijk 索引符为关于 123 的偶次序 (1.23)
−1 ijk 索引符为关于 123 的奇次序
其中偶次序为偶数次交换相邻检索符得到 123，231，312，奇次序为奇数次交换交换相邻检
索符得到 132，213，321。
1.3.6 其他张量操作
下面我们描述 OpenFOAM 提供的一些不太常用的张量操作：
 平方：某个张量和自己本身的外积，例如对于矢量 a 的平方为：𝒂2 = 𝒂𝒂；

 n 次幂： 𝒂3 = 𝒂𝒂𝒂；
某个张量和自己本身的 n 次外积，例如对于矢量 a 的 3 次幂为：
 模的平方：某个张量和自己本身的阶次内积，例如对于 2 阶张量 T 有|𝑻|2 = 𝑻 ∶ 𝑻；
 模：为模的平方的平方根，例如对于张量 T，|𝑇|2 = √𝑻 ∶ 𝑻。单位大小的矢量我们
称为单位矢量；
 最大元素：张量的最大元素（包括符号），注意：非最大的模；
 最小元素：张量的最小元素；
 平均元素：张量元素的平均；
 缩放：正如名称所言，缩放用于使用某个张量的元素对另一个张量的元素进行缩放，
这两个张量的阶数需要一致，例如：对 b 矢量使用 a 矢量进行缩放后为矢量 c，计
算如下：
𝑐𝑖 = scale(𝑎, 𝑏) = (𝑎1 𝑏1 , 𝑎2 𝑏2 , 𝑎3 𝑏3 ) (1.24)
1.3.7 几何变形和单位张量
2 阶张量可以定义为矢量的线性函数，例如：向量 a 和向量 b 可以通过 T 联系起来𝐚 = 𝐓 ⋅ 𝐛。

我们可以通过改变 T 中的分量来使 T 能够对张量实现不同的坐标变换，例如从 x、y、z 的
坐标系统变换为 x*、y*、z*坐标系统。这种变换称之为张量转换。然而，标量在坐标系变换
时不发生变化。矢量可实现下述变换：
𝐚∗ = 𝐓 ∙ 𝐚 (1.25)
二阶张量 S 可以通过下述公式变换成为 S*：
𝐒∗ = 𝐓 ∙ 𝐒 ∙ 𝐓 T (1.26)
单位张量定义为：它将某个张量变换为自己本身，即：
17
𝐚=𝐈∙𝐚 (1.27)
因此，单位张量为：
1 0 0
𝐈 = δij = (0 1 0) (1.28)
0 0 1
其中δij 为 Kronecker 标示符。
1.3.8 张量等式
下面我们来列举一下常用的张量等式，我们假定相关的导数都是存在且连续的。其中 s 表示
标量，a 表示矢量：
∇ ⋅ (∇ × 𝐚) = 0
∇ × (∇s) = 0
∇ ⋅ (s𝐚) = s∇ ⋅ 𝐚 + 𝐚 ⋅ ∇s
∇ × (s𝐚) = s∇ × 𝐚 + 𝐚 × ∇s
∇(𝐚 ⋅ 𝐛) = 𝐚 × (∇ × 𝐛) + 𝐛 × (∇ × 𝐚) + (𝐚 ⋅ ∇)𝐛 + (𝐛 ⋅ ∇)𝐚 (1.29)
∇ ⋅ (𝐚 × 𝐛) = 𝐛 ⋅ (∇ × 𝐚) − 𝐚 ⋅ (∇ × 𝐛)
∇ × (𝐚 × 𝐛) = 𝐚(∇ ⋅ 𝐛) − 𝐛(∇ ⋅ 𝐚) + (𝐛 ⋅ ∇)𝐚 − (𝐚 ⋅ ∇)𝐛
∇ × (∇ × 𝐚) = ∇(∇ ⋅ 𝐚) − ∇2 𝐚
(∇ × 𝐚) × 𝐚 = 𝐚 ⋅ (∇𝐚) − ∇(𝐚 ⋅ 𝐚)
了解下面这个𝑒 − 𝛿等式有助于用户对张量在标识符下的操作有更好的理解：
𝑒𝑖𝑗𝑘 𝑒𝑖𝑟𝑠 = 𝛿𝑗𝑟 𝛿𝑘𝑠 − 𝛿𝑗𝑠 𝛿𝑘𝑟 (1.30)
1.3.9 二阶张量算术
下面我们列举几个 2 阶张量的算法：
 转置：正如方程 2.1 所说，𝐓 = 𝑇𝑖𝑗 ，𝐓 T = 𝑇𝑗𝑖 ；

 symm、skew：正如 1.2 节所说，我们有对称张量即其关于对角线是对称的：𝐓 =
𝐓 T 。非对称张量（skew 张量）符合这种特性：𝐓 = −𝐓 T，这意味着𝑇11 = 𝑇22 = 𝑇33 =
0，每个二阶张量都可以分为 symm 部分和 skew 部分即：
⏟ + 𝐓 T ) + 0.5 (𝐓
𝐓 = 0.5 (𝐓 ⏟ − 𝐓 T ) = symm 𝐓 + skew 𝐓 (1.31)
symm skew
 迹：张量的迹为标量（对角线元素的和）：
tr 𝐓 = 𝑇11 + 𝑇22 + 𝑇33 (1.32)
 对角：由张量的对角线元素构成的矢量：
diag 𝐓 = (𝑇11 , 𝑇22 , 𝑇33 ) (1.33)
 dev、hyd：每个二阶张量都可分解为 dev 部分（其迹为 0）和 hyd 部分（𝐓 = 𝑠𝐈，
s 为标量）：
18
1 1
⏟ − (tr 𝐓)𝐈 + ⏟(tr 𝐓)𝐈 = dev 𝐓 + hyd 𝐓
𝐓=𝐓
3 3
(1.34)
dev hyd
 行列式：二阶张量的 det 这样计算：
𝑇11 𝑇12 𝑇13 𝑇11 (𝑇22 𝑇33 − 𝑇23 𝑇32 ) −
1
det 𝐓 = |𝑇21 𝑇22 𝑇23 | = 𝑇12 (𝑇21 𝑇33 − 𝑇23 𝑇31 ) + = 𝑒𝑖𝑗𝑘 𝑒𝑝𝑞𝑟 𝑇𝑖𝑝 𝑇𝑗𝑞 𝑇𝑘𝑟
6
(1.35)
𝑇31 𝑇32 𝑇33 𝑇13 (𝑇21 𝑇32 − 𝑇22 𝑇31 )
 代数余子式：张量的余子式通常这样计算：对于公式 1.35 中的 3×3 行列式，删除
元素所在的行和列，剩下元素取2 × 2行列式，例如𝑇12 的余子式为：
𝑇 𝑇23
| 21 | = T21 𝑇33 − 𝑇23 𝑇31 (1.36)
𝑇31 𝑇33
代数余子式是带符号的余子式，符号这样计算：
+, 𝑖 + 𝑗为偶数
( (1.37)
−, 𝑖 + 𝑗为奇数
T 的代数余子式即为：
1
cof 𝐓 = 𝑒𝑗𝑘𝑟 𝑒𝑖𝑠𝑡 𝑇𝑠𝑘 𝑇𝑡𝑟 (1.38)
2
 逆：张量的逆这样计算：
cof 𝐓T
inv 𝐓 = det 𝐓
(1.39)
 霍奇对偶：张量的霍奇对偶这样计算：
∗ 𝐓 = (𝑇23 , −𝑇13 , 𝑇12 ) (1.40)
1.3.10 标量算术
OpenFOAM 支持大部分的标量算术例如平方根、指数、对数、正弦余弦等，它们在表 1.2 中

列出。OpenFOAM 还提供了 3 个针对标量的附加函数即：
 Sign：
1 如果𝑠 ≥ 0
sgn (𝑠) = { (1.41)
−1 如果𝑠 < 0
 Positive：
1 如果𝑠 ≥ 0
pos (𝑠) = { (1.42)
0 如果𝑠 < 0
 Limit：
19
𝑠 如果𝑠 < n
limit (𝑠, 𝑛) = { (1.43)
0 如果𝑠 ≥ 0
1.4 OpenFOAM 张量类
在 OpenFOAM 中，primitive 类包含了距今为止所有的有关张量算术的计算。表 1.1 列举了

OpenFOAM 提供的基本的张量类型，表中也提供了访问张量元素的函数：
阶数名字基本类访问函数
0 标量 scalar
1 矢量 vector x()，y()，z()
2 张量 tensor xx()，xy()，xz()
表 1.1 OpenFOAM 张量基本类
在 OpenFOAM 中我们可以把张量：
1 2 3
T = (4 5 6) (1.44)
7 8 9
这样定义：
tensor T(1, 2, 3, 4, 5, 6, 7, 8, 9);
我们可以利用 xz()函数来访问𝑇13 或者𝑇𝑥𝑧 ，例如：
Info << “Txz = ” << T.xz() << endl;
这会输出：
Txz = 3
1.4.1 OpenFOAM 张量操作
1.3 节描述的所有张量操作在 OpenFOAM 中都有特定的语法来实现，并且尽可能的和相对

应的算符名称或图形相近。一些函数的名字直接使用操作符的名称，例如 symm()，其他的
一些可以用操作符直接实现，例如 *，所有的函数参见表 1.2：
操作描述数学形式 OpenFOAM 中的实现
加法 𝐚+𝐛 a+b
减法 𝐚−𝐛 a-b
20
乘以标量 𝒔𝐚 s*a
除以标量 𝐚∕𝒔 a/s
外积 a 和 b 的阶数≥ 1 𝐚𝐛 a*b
内积 a 和 b 的阶数≥ 1 𝐚⋅𝐛 a&b
双内积 a 和 b 的阶数≥ 2 𝐚∶𝐛 a && b
叉乘 a 和 b 的阶数= 1 𝐚×𝐛 a^b
𝟐
平方 𝐚 sqr(a)
模的平方 |𝐚|𝟐 magSqr(a)
模 |𝐚| mag(a)
𝒏
幂 n = 0,1,…,4 𝐚 pow(a,n)
元素均值 i = 0,1,…,N 𝑎̅𝒾 cmptAv(a)
元素最大值 i = 0,1,…,N max(𝑎𝒾 ) max(a)
元素最小值 i = 0,1,…,N min(𝑎𝒾 ) min(a)
缩放 scale(𝐚, 𝐛) scale(a,b)
2 阶张量操作
转置 𝐓T T.T()
对角 diag𝐓 diag(T)
迹 tr𝐓 tr(T)
Dev dev𝐓 dev(T)
Sym symm𝐓 symm(T)
Skew-sym skew𝐓 skew(T)
行列式 det𝐓 det(T)
代数余子式 cof𝐓 cof(T)
逆 inv𝐓 inv(T)
∗
霍奇对偶 𝐓 *T
标量算术
取正负 sgn(𝑠) sign(s)

取正 𝑠 >= 0 pos(s)
取负 𝑠<0 neg(s)
极限 lim(𝑠, 𝑛) limit(s,n)
平方根 √𝑠 sqrt(s)
e exp 𝑠 exp(s)
对数 ln 𝑠 log(s)
10 为基的对数 log10 𝑠 log10(s)
sin asin 𝑠 sin(s)
cos cos 𝑠 cos(s)
tan tan 𝑠 tan(s)
arcsin asin 𝑠 asin(s)
arccos acos 𝑠 acos(s)
21
arctan atan 𝑠 atan(s)

sinh sinh 𝑠 sinh(s)
cosh cosh 𝑠 cosh(s)
tanh tanh 𝑠 tanh(s)
arcsinh asinh 𝑠 asinh(s)
arccosh acosh 𝑠 acosh(s)
arctanh atanh 𝑠 atanh(s)
误差函数 erf 𝑠 erf(s)
C 误差函数 erfc 𝑠 erfc(s)
Gamma 对数函数 ln Γ𝑠 lgamma(s)
0 阶 1 型 Bessel 函数 J0 𝑠 j0(s)
1 阶 1 型 Bessel 函数 J1 𝑠 j1(s)
0 阶 2 型 Bessel 函数 Y0 𝑠 y0(s)
1 阶 2 型 Bessel 函数 Y1 𝑠 y1(s)
a 和 b 为任意阶数的张量（除非进行指定）
，s 为标量，N 为元素数
表 1.2 OpenFOAM 中的张量操作
1.5 量纲单位
在连续介质里，物理量都具有量纲。例如质量：kg, 体积：m³，压力：Pa。这些物理量的代
数几何操作要求量纲保持一致。对于某些同样量纲的物理量，只有加、减、等于是具有物理
意义的。为了防止对物理量进行无意义或者错误的运算，OpenFOAM 推荐用户将量纲依附
于场数据以及某个物理量，这样在对其进行张量操作的时候，OpenFOAM 会执行量纲检查。
量纲我们用 dimensionSet 类来定义，例如：
dimensionSet pressureDims(1, -1 ,-2, 0, 0, 0, 0);
编号名称 SI 单位 USCS 单位
1 质量 kg lbm
2 米 m ft
3 时间 s s
4 开尔文 K 。
𝑅
5 摩尔质量 kgmol lbmol
6 电流 A A
7 光强 cd cd
表 1.3 S.I. 量纲
表 1.3 列举了 OpenFOAM 植入的 S.I.量纲制。我们上个范例中的 pressureDims 使用 dimension-

Set 把单位设置为 kg m-1 s-2，因为我们 pressureDims 序列中第一个值为 1，这对应 k1，第二
22
个值为-1，这对应 m-1，以此类推。附带单位的张量使用 dimensioned<Type>模板类来定义，

其中 type 是 scalar，vector，tensor 等。dimensioned<Type>将存储一个 word 类、一种<Type>
类型、和 dimensionSet：
dimensionedTensor sigma
(
"sigma",
dimensionSet (1, -1, -2, 0, 0, 0, 0),
tensor (1e6,0,0,0,1e6,0,0,0,1e6),
);
上述这段代码将会定义这样一个张量（应力单位）
：
106 0 0
σ=( 0 106 0 ) (1.45)
0 0 106
23
第二章
离散
目前我们处理的仅仅是张量算术。但是我们需要求解的是关于时间和空间的偏微分方程组。
因此我们将拓展我们的张量为张量场，例如：一个关于时间和空间变化的张量。在这一章节
我们首先对我们可能会遇到的偏微分算符进行描述。然后我们说明张量场在 OpenFOAM 中
是如何构建的，以及如何将张量场离散成为方程组。
2.1 微分算符
在我们处理导数之前，我们首先通过索引符𝜕𝑖 引入 nabla 算符∇：

𝜕 𝜕 𝜕 𝜕
∇= ∂i = =( , , ) (2.1)
𝜕𝑥𝑖 𝜕𝑥1 𝜕𝑥2 𝜕𝑥3
nabla 算符在使用的时候遵循以下原则：
 当 nabla 算符作用在乘积上的时候，其和乘积的导数相同。例如：∂i 𝑎𝑏 = (𝜕𝑖 𝑎)𝑏 +

𝑎(𝜕𝑖 𝑏)；
 在其它情况下 nabla 算符的运算和一般的向量相同；
2.1.1 梯度
如果定义一个标量场 s 并且 s 为连续可导的，我们有矢量场∇𝑠：
𝜕𝑠 𝜕𝑠 𝜕𝑠
∇𝑠 = ∂i 𝑠 = (𝜕𝑥 , 𝜕𝑥 , 𝜕𝑥 ) (2.2)
1 2 3
梯度我们可以用在任何张量场上，得到的是更高一阶的张量场，例如，矢量场 a 的梯度就是
二阶张量场：
𝜕𝑎1 /𝜕𝑥1 𝜕𝑎1 /𝜕𝑥2 𝜕𝑎1 /𝜕𝑥3

∇𝐚 = ∂i 𝑎𝑖 = (𝜕𝑎2 /𝜕𝑥1 𝜕𝑎2 /𝜕𝑥2 𝜕𝑎2 /𝜕𝑥3 ) (2.3)
𝜕𝑎3 /𝜕𝑥1 𝜕𝑎3 /𝜕𝑥2 𝜕𝑎3 /𝜕𝑥3
24
2.1.2 散度
如果矢量场 a 是连续可导，那么 a 的散度场为标量场：

𝜕𝑎1 𝜕𝑎2 𝜕𝑎3
∇ ⋅ 𝐚 = ∂i 𝑎𝑖 = + + (2.4)
𝜕𝑥1 𝜕𝑥2 𝜕𝑥3
散度可以在任何一阶以及以上阶张量场上操作，结果为低一阶的张量场。例如，二阶张量场
T 的散度为矢量场：
𝜕𝑇11 /𝜕𝑥3 + 𝜕𝑇12 /𝜕𝑥3 + 𝜕𝑇13 /𝜕𝑥3
∇ ⋅ 𝐓 = ∂i 𝑇𝑖𝑗 = ( 𝜕𝑇21 /𝜕𝑥3 + 𝜕𝑇22 /𝜕𝑥3 + 𝜕𝑇23 /𝜕𝑥3 ) (2.5)
𝜕𝑇31 /𝜕𝑥3 + 𝜕𝑇32 /𝜕𝑥3 + 𝜕𝑇33 /𝜕𝑥3
2.1.3 旋度
如果矢量场 a 是连续可导，那么 a 的旋度场为矢量场：

𝜕𝑎3 𝜕𝑎2 𝜕𝑎1 𝜕𝑎3 𝜕𝑎2 𝜕𝑎1
∇ × 𝐚 = eijk 𝜕𝑗 𝑎𝑘 = ( − , − , − ) (2.6)
𝜕𝑥2 𝜕𝑥3 𝜕𝑥3 𝜕𝑥1 𝜕𝑥1 𝜕𝑥2
旋度和散度的关系为：
∇ × 𝐚 = 2(∗ skew∇𝐚) (2.7)
2.1.4 拉普拉斯
拉普拉斯操作在数学上可以看做是散度和梯度的组合，即：∇2 = ∇ ⋅ ∇。然而，拉普拉斯应该
看作为是单步操作：将张量场转化为同阶张量场的操作，而不是两步张量场操作的结合（先
增加一阶然后再减去一阶）。
实际上，拉普拉斯最好定义为标量操作符的形式（正如我们定义∇操作符一样）：
𝜕2 𝜕2 𝜕2
∇2 = ∂2 = + + (2.8)
𝜕𝑥12 𝜕𝑥22 𝜕𝑥32
例如，对于一个标量场 s 的拉普拉斯场同样也为一个标量场：
𝜕2𝑠 𝜕2 𝑠 𝜕2 𝑠
∇2 𝑠 = 𝜕 2 𝑠 = + + (2.9)
𝜕𝑥12 𝜕𝑥22 𝜕𝑥32
2.1.5 时间导数
存在不止一种有关张量场的时间导数定义。为了准确的描述时间导数我们必须首先知道张量
是和一个移动的某个体积下的物理量有关的。如果我们追踪这个无限小的微元（或者微粒）
，
它是即时移动的并且物理张量𝜙会跟随时间而发生变化，我们定义物质导数为：
2
2 拉格朗日方法
25
𝐷𝜙 ∆𝜙
𝐷𝑡
= lim (2.10)
∆𝑡→0 ∆𝑡
在连续介质问题中，特别是流体力学中，我们通常在某个固定点来观测某个物理张量跟随时
间的变化形式3。我们用∂/ ∂𝑡来表示在空间的某个点下的时间导数项，我们用它可以表示物
质导数：
𝐷𝜙 𝜕𝜙
𝐷𝑡
= 𝜕𝑡
+ 𝐔 ⋅ ∇𝜙 (2.11)
其中 U 是某个物理量𝜙的速度，2.11 公式右边的第二项我们称为𝜙的对流率。
2.2 离散简述
离散的意思就是用一系列离散的量来近似代表一系列的连续值。有限体积法以及其他的例如
有限元、有限差分在离散的时候都遵循下面的原则：
 空间离散：通过一系列的点来定义求解域，它填充整个空间并且是有界的；
 时间离散：对于瞬态问题，将计算时间分为有限的小的时间间隔；
 方程离散：针对某个物理问题特定的偏微分方程，生成一系列的代数方程来代表在
空间离散后的物理量；
2.2.1 OpenFOAM 链表和场
OpenFOAM 需要存储大量的数据使用函数对数据进行大量的操作。因此 OpenFOAM 提供模

板类 List<Type>，其可以通过类模板创建一系列的 Type 对象来继承其特性。例如，vector 的
List 即为 List<vector>。
张量链表在 OpenFOAM 中通过 Field<Type>来指定。为了使代码更加的清晰易懂，所有
的 Field<Type>类型，例如 Field<vector>，都使用 typedef 来重新定义为 scalarField，vectorField，
tensorField，symmTensorField，tensorThirdField 以及 symmTensorThirdField。张量场上的代
数算术有一些显而易见的数学规则，例如张量场的元素数量必须一致等。OpenFOAM 也支
持单一场和一个数值的操作。例如，所有的速度场都乘以标量 2，即：𝐔 = 2.0 × 𝐔。
2.3 计算域离散
计算域的离散参考图 2.1 所示。首先我们将整个计算空间离散为网格，在这个计算网格基础

上进行偏微分方程的离散。如果时间也需要进行离散的话是非常简单的。将时间划分为一系
列的时间步∆𝑡（可以随着计算而改变大小）即可。时间步长一般和计算中的某些条件相对应。
更详细来讲，空间离散需要将整个计算域划分为非常多小的网格（有限控制体）。这些
3 欧拉方法
26
网格互相是连续相连的，例如他们不能互相重合，并且他们需要填充整个计算域。典型的空
间离散后的网格如图 2.1 所示：
时间域
表 2.1 计算域离散
表 2.2 有限体积离散的相关参数
27
即使某些物理量以及其他的相关参数需要储存在面或者顶点上，他们大部分都储存在网格中
心 P。网格单元通过一系列的平面来界定，我们用 f 来表示面。在 OpenFOAM 中，网格单
元面的数量是没有限制的，这种网格我们通常称之为任意非结构网格。这种叫法可以和某些
网格（这些网格的网格单元面和坐标轴相依赖）来进行区分。使用这种网格来进行计算的计
算流体力学代码允许非常灵活的网格生成和操作。这在几何体非常复杂以及需要处理动网格
的时候优势明显。
虽然大部分物理量存储在网格中心，但是一些物理量依然存储在网格单元面，在
OpenFOAM 中具有两种网格单元面：
 内部面：连接两个网格单元的面（且仅仅连接两个网格单元）。对于每个内部面，
OpenFOAM 将其连接的两个网格单元一个指定为 owner，一个指定为 neighbor；
 边界面：从属一个网格单元的面，且和计算域边界重合。这些面只有一个 owner 网
格单元；
2.3.1 定义 OpenFOAM 网格
OpenFOAM 中使用不同等级来进行网格描述，最基本的网格类为 polyMesh，之所以叫这个

名字是因为它建立在多面体基础之上。描述 polyMesh 需要使用下面的相关信息来定义几何
和网格，详见图 2.3：
 Points：网格单元点位置矢量链表（vectorField），使用 typedef 定义后的名称为

pointField；
 Faces：一系列网格单元面列表 List<face>，也即 faceList，其中 face 类通过 pointField

的顶点标号相关信息来定义；
 Cells：一系列网格单元列表 List<cell>，也即 cellList，其中 cell 类通过 faceList 的

面标号相关信息来定义；
 Boundary：polyBoundaryMesh 分解为一系列的 patch，polyPatchList 代表不同的边

界区域。边界依靠这种方式来进行区分这样可以使用不同的边界条件。所有的
polyPatch 面都作为一个单一 block 存储在 faceList 中，这样他的面可以使用 slice 类
（存贮 block 的第一个面和最后一个面的引用）来访问。每个 polyPatch 需要从下
面两个信息来组建：
 slice：见上文；
 word：名字；
有限体积离散需要使用存储在 polyMesh 上的几何信息。因此 OpenFOAM 把 polyMesh 拓展

为 fvMesh 类，其存储用于有限体积离散的相关信息。参见图 2.1，fvMesh 从 polyMesh 来组
建，fvMesh 可以在运行的时候进行更新（例如进行网格细化、动网格的算例）。
28
图 2.3 OpenFOAM 的网格描述纲要
2.3.2 定义 geometricField
目前为止我们定义了场（张量链表）和网格。他们可以通过在我们的计算域定义张量场来结
合起来，这通过 OpenFOAM 中的模板类 geometricField<Type>来实现。Field 的值区分为计
算域的内部值（计算域内部网格单元中心的值）和边界值（计算域的边界面的值）。geometric
-Field<Type>存储下列信息：
 Internal field：简单的 Field<Type>类型，详见 2.2.1 节；
 Boundary field：一种 GeometricBoundaryField，需要为每个 patch 的每个网格单元

面定义 Field。于是这就是场的场，存储在 FieldField<Type>里，fvBoundaryMesh 引
用也存储在此处；
 Mesh：fvMesh 的引用，其会添加一下附加信息，例如场信息存储在网格单元中心
还是面上；
 Dimensions：一个 dimensionSet 类型，详见《OpenFOAM 用户指南》4.2.6 节；
 Old values：在时间离散中上一个时间步的场信息。geometricField<Type>会存储上
一个时间步（过去的时间步，甚至上上个时间步）的引用信息；
 Previous iteration values：在使用低松弛求解的时候需要存储之前迭代的信息。因

此，如果需要的话，geometricField<Type>会存储之前迭代的所需信息；
29
类描述符号成员函数
volScalarField 网格体积4 V V()

surfaceVectorField 面矢量 𝐒f Sf()
surfaceScalarField 面大小 |𝐒f | magSf()
volVectorField 网格中心 𝐂 C()
surfaceVectorField 面心 𝐂f Cf()
surfaceScalarField 面通量 ϕg phi()
表 2.1 fvMesh 存储的信息
正如 2.3 节我们讨论的，我们主要将物理量定义在网格单元中心，一部分定义在网格面，很
少的一部分我们定义在网格单元顶点上。geometricField<Type>也使用 typedef 进行了重新命
名来表示定义场的位置信息：
 volField<Type>：网格单元中心场；
 surfaceField<Type>：网格单元面场；
 pointField<Type>：网格单元顶点场；
这些通过 typedef 重新命名的 geometricField<Type>请查阅图 2.4，geometricField<Type>继承

了所有的 Field<Type>数学操作，并且可以使用 dimensionSet 进行量纲检查。它也可以使用
下述章节的有限体积离散的相关函数。geometricField<Type>的类结构请参考图 2.55。
2.4 方程离散
偏微分方程组可以转化为一系列的代数方程组，如下所示：
[𝐴][𝑥] = [𝑏] (2.12)
其中[𝐴]是方阵，[𝑥]是相关物理量的未知矢量，[𝑏]是源项矢量。我们把[𝑥]和[𝑏]称之为矢量
不是因为他们是具有三个元素的矢量，而是他们是一个和矩阵相关的专业词汇。他们是定义
在计算域的一系列的值，例如一个 geometricField<Type>类型，在使用有限体积离散的时候
一般为 volField<Type>类型。
[𝐴]是代数方程组的系数矩阵，它不能用 geometricField<Type>来描述。因此它们使用他
们自己的类进行描述：fvMatrix。fvMatrix<Type>通过对 geometricField<Type>离散进行创建
因此会继承<Type>的特性。他们可以使用常规的代数方程组操作例如加法，减法，乘法之类。
在 OpenFOAM 中，偏微分方程中的每个项都需要使用静态函数类 finiteVolumeMethod
以及 finiteVolumeCalculus 进行表示，简化的写法即为：fvm，fvc。
4此列并不是第一列的翻译
5这个图表并不是 geometricField<Type>的完整描述，在这里只是简单的介绍 geometricField<Type>的主要
衍生类
30
图 2.4 定义在 2 个边界 patch 上的 geometricField<Type>类型

（使用 2 维几何进行简要描述）
31
图 2.5 geometricField<Type>的基本类结构
32
fvm 和 fvc 都相应的包含各自的函数以代表相对应的微分算符，例如∇2 ，∇ ⋅ ， ∂/ ∂t，我们

区分 fvm 和 fvc 的目的在于：
 fvm 对偏微分方程的项进行隐形离散，返回一个 fvMatrix<Type>；
 fvc 的某些函数计算显性导数以及一些其他的显性操作，返回一个 geometricField -

<Type>；
图 2.6 即为定义在网格上的 geometricField<Type>，其有两个边界 patch，为了简化说明如何

将一个类型转化为另一个类型我们使用 2 维图形进行表示：
图 2.6 geometricField<Type>及其数学操作
表 2.2 列举了用于进行偏微分方程离散的 fvm 和 fvc 可用的大部分函数。在进行有限体积离

散的时候，每个偏微分方程的项首先需要在有限控制体 V 内进行积分。大部分空间导数项
都使用高斯定律来将体积分转化为面积分：
∫𝑉 ∇ ⊕ 𝜙 𝑑𝑉 = ∫𝑆 𝑑𝑺 ⊕ 𝜙 (2.13)
其中𝑺是面矢量，𝜙表示任意的张量场，⊕表示任意的张量积，例如内积、外积、叉乘。相对
应的就是散度：∇ ⋅ 𝜙，梯度∇𝜙，以及∇ × 𝜙。在积分之后，我们需要使用合适的格式来进行
线性化，这将在后续章节有所提及。OpenFOAM 提供了很多种格式进行离散，但某些项经
常使用一种格式来离散。离散格式的选择可以在代码中直接指定，或者可以从输入文件中指
定，他们存储在 fvSchemes 类中。
33
离散项显性/隐性数学形式代码实现
拉普拉斯项 Imp/Exp ∇2 𝜙 laplacian(phi)

∇ ⋅ Γ∇𝜙 laplacian(Gamma, phi)
时间导数 Imp/Emp 𝜕𝜙
𝜕𝑡 ddt(phi)
𝜕𝜌𝜙 ddt(rho,phi)
𝜕𝑡
二阶时间导数 Imp/Exp 𝜕 𝜕𝜙
(𝜌 ) d2dt2(rho, phi)
𝜕𝑡 𝜕𝑡
对流项 Imp/Emp ∇ ⋅ (𝜓) div(psi,scheme)*
∇ ⋅ (𝜓𝜓) div(psi, phi, word)*
div(psi, phi)
散度项 Exp ∇⋅χ div(chi)
梯度项 Exp ∇χ grad(chi)
∇𝜙 gGrad(phi)
lsGrad(phi)
snGrad(phi)
snGradCorrection(phi)
梯度梯度项 Exp |∇∇𝜙|2 sqrGradGrad(phi)
旋度项 Exp ∇×𝜙 curl(phi)
源项 Imp 𝜌𝜙 Sp(rho,phi)
Imp/Exp SuSp(rho,phi)
fvm::SuSp 通过判断 rho 的符号进行显性或者隐性离散；

显性项离散可以通过 vol<Type>Field 引入，例如 rho*phi；
函数参数可以为以下类型：
phi：vol<Type>Field；
Gamma：scalar volScalarField, surfaceScalarField, volTensorField, surfaceTensorField；
rho：scalar, volScalarField；
psi：surfaceScalarField；
chi：surface<Type>Field, vol<Type>Field；
表 2.2 OpenFOAM 偏微分方程离散

34
2.4.1 拉普拉斯项
在有限控制体内，拉普拉斯项这样积分以及线性化：
∫𝑉 ∇ ⋅ (Γ∇𝜙) 𝑑𝑉 = ∫𝑆 𝑑𝐒 ⋅ (Γ∇𝜙) = ∑𝑓 Γ𝑓 𝐒𝑓 ⋅ (∇𝜙)𝑓 (2.14)
当网格单元中心 P 和相邻网格单元中心 C 的长度矢量 d 和这两个网格单元的内部面正交的

时候，即 d 平行于𝐒𝐟 的时候，面梯度这样计算6：
𝜙𝑁 −𝜙𝑃
𝐒𝑓 ⋅ (∇𝜙)𝑓 = |𝑆𝑓 | |𝐝|
(2.15)
在网格非正交的情况下，我们需要引入一个显性项7。
2.4.2 对流项
在有限控制体内，对流项这样积分以及线性化：
∫𝑉 ∇ ⋅ (𝜌𝐔𝜙 ) 𝑑𝑉 = ∫𝑆 𝑑𝐒 ⋅ (𝜌𝐔𝜙 ) = ∑𝑓 𝐒𝑓 (𝜌𝐔)𝑓 𝜙𝑓 = ∑𝑓 𝐹𝜙𝑓 (2.16)
面值𝜙f 可以使用一系列的格式进行计算：
 中心差分（CD）：无界二阶精度：
𝜙𝑓 = 𝑓𝑥 𝜙𝑃 + (1 − 𝑓𝑥 )𝜙𝑁 (2.17)
其中𝑓𝑥 = ̅̅̅̅ ̅̅̅̅是面 f 和网格单元中心 N 的距离，𝑃𝑁

̅̅̅̅，𝑓𝑁
𝑓𝑁 /𝑃𝑁 ̅̅̅̅是网格单元中心 P 和
N 的距离；
 迎风格式（UD）：牺牲精度来确保有界，依据流向决定𝜙f 的值：
𝜙𝑃 for 𝐹 ≥ 0
𝜙𝑓 = { (2.18)
𝜙𝑁 for 𝐹 < 0
 混合格式（BD）：即中心差分和迎风格式的混合，力图在合理的精度下保证有界：
𝜙𝑓 = (1 − 𝛾)(𝜙𝑓 ) + 𝛾(𝜙𝑓 ) (2.19)

𝑈𝐷 𝐶𝐷
OpenFOAM 植入了一些 Gamma 差分格式，他们可以选择混合系数γ，同样，Open

-FOAM 也提供了一些其他的著名格式例如 vanLeer，SUPERBEE，MINMOD 等；
6 此公式类似一种格式
7 请参考 Hrv 的博士论文
35
2.4.3 一阶时间导数
在有限控制体内，一阶时间导数项这样积分：
𝜕
∫ 𝜌𝜙 𝑑𝑉
𝜕𝑡 𝑉
(2.20)
这一项需要将时间进行简单的离散，在这里我们引进几个项：
 新值：𝜙 n = 𝜙(𝑡 + ∆𝑡)，即我们在当前时间步需要求的物理量；
 旧值：𝜙 0 = 𝜙(𝑡)，我们在上一个时间步存储的物理量；
 旧-旧值：𝜙 00 = 𝜙(𝑡 − ∆𝑡)，我们在上上一个时间步存储的物理量；
一阶或二阶时间导数项的离散可以在相关的字典文件中通过 timeScheme 关键词来

进行定义，详情请参考《OpenFOAM 用户指南》4.4 节；
 欧拉显性：这样指定：timeScheme EulerImplicit，其为一阶精度：
𝜕 (𝜌𝑃 𝜙𝑃𝑉)𝑛 −(𝜌𝑃𝜙𝑃 𝑉)𝑜

∫ 𝜌𝜙 𝑑𝑉 =
𝜕𝑡 𝑉 ∆𝑡
(2.21)
 向后差分：这样指定：timeScheme BackwardDifferencing，其为需要存储旧-旧值的
二阶精度格式，因此相对于欧拉显性需要更大的内存空间：
𝜕 3(𝜌𝑃 𝜙𝑃 𝑉)𝑛 −4(𝜌𝑃 𝜙𝑃 𝑉)𝑜 +(𝜌𝑃 𝜙𝑃 𝑉)∞

∫ 𝜌𝜙 𝑑𝑉
𝜕𝑡 𝑉
= 2∆𝑡
(2.22)
2.4.4 二阶时间导数项
在有限控制体内，二阶时间导数项这样积分（其为一阶精度）
：
𝜕 𝜕𝜙 (𝜌𝑃 𝜙𝑃 𝑉)𝑛 −2(𝜌𝑃 𝜙𝑃 𝑉)𝑜 +(𝜌𝑃 𝜙𝑃 𝑉)∞

∫ 𝜌 𝑑𝑉 =
𝜕𝑡 𝑉 𝜕𝑡 ∆𝑡 2
(2.23)
2.4.5 散度项
这一节我们所说的散度为显性项（以区别 2.4.2 所提及的对流项），这一节的散度项并不是速

度项和某个物理量的乘积。在有限控制体内，散度项这样积分：
∫𝑉 ∇ ⋅ 𝜙 𝑑𝑉 = ∫𝑆 𝑑𝐒 ⋅ 𝜙 = ∑𝑓 𝐒𝑓 ⋅ 𝜙𝑓 (2.24)
36
其中 fvc::div 的参数可以为 surface<Type>Field（需要指定𝜙f ），或者 vol<Type>Field，其需要

使用 2.4.10 节提及的中心差分格式来进行面插值。
2.4.6 梯度项
梯度项是可以通过很多方式进行计算的显性项。格式可以这样指定：要么通过选择某个格式
的梯度函数来计算，例如 fvc::gGrad，fvc::lsGrad，要么通过 fvc::grad 以及输入文件的 grad -
Scheme 来计算。
 高斯积分：通过 fvc::grad 函数（并在 gradScheme 里面指定 Gauss）指定，或者使

用 fvc::gGrad 函数指定。离散的时候需要在有限控制体内使用高斯定律进行积分：
∫𝑉 ∇𝜙 𝑑𝑉 = ∫𝑆 𝑑𝐒 𝜙 = ∑𝑓 𝐒𝑓 𝜙𝑓 (2.25)
和 fvc::div 函数一样，高斯积分的 fvc::grad 函数的参数可以为 surface<Type>Field

或者 vol<Type>Field。
 最小二乘法：建立在下述规则基础之上：
1. 从点 P 的值，使用 P 点的梯度可以推算 N 点的值；

2. N 点的推断值可以和 N 点的实际值进行比较，其值的不同即为误差；
3. 如果我们使用某种梯度使得 P 节点和临点 N 的梯度的平方和最小，我们认
为这个梯度假定是精准的；
最小二乘法可以通过 fvc::grad 函数（并在 gradScheme 里面指定 leastSquares）

指定，
或者使用 fvc::lsGrad 函数指定。离散首先计算每个 P 点的张量 G
G = ∑𝑁 ω2𝑁 dd (2.26)
d 是 P 点到 N 点的矢量，且权函数wN = 1/|𝐝|。梯度可以这样计算：
(∇𝜙)𝑃 = ∑𝑁 ω2𝑁 G−1 ⋅ d(𝜙𝑁 − 𝜙𝑃 ) (2.27)
 面法相梯度：面法向梯度nf ⋅ (∇𝜙)f 可以使用下面的公式计算：
𝜙𝑁 −𝜙𝑃
(∇𝜙)𝑓 = |𝐝|
(2.28)
梯度项通过 fvc::snGrad 调用并返回一个 surfaceField<Type>，和 2.4.1 节的拉普拉

斯离散一样，∇𝜙f 使用上述公式进行计算。同样的，在非正交网格中，面法相梯度
在计算的时候我们引入一个线性的修正项。在代码中我们可以调用 fvc::snGrad –
Correction。
37
2.4.7 梯度-梯度的平方
梯度-梯度的平方我们这样计算：首先对场取散度；再对这个取散度后的场继续取散度；然
后计算结果的模的平方。在数学上，梯度梯度的平方这样表示：|∇(∇𝜙)2 |。
2.4.8 旋度
旋度可以从梯度计算得来，参见 2.4.6 节。首先我们计算梯度，然后通过下面这个公式（公

式 2.7）来计算旋度：
∇ × 𝜙 = 2 ∗ (skew∇𝜙)
2.4.9 源项
源项的指定有三种方法：
 显性：每个显性的源项为一个 volField<Type>。因此，显性的源项在偏微分方程中
可以简单的看做一个值。例如，如果我们打算求解泊松方程∇2 𝜙 = 𝑓，我们可以定
义 phi 和 f 为一个 volScalarField，然后使用下面的代码：
solve(fvm::laplacian(phi) == f)
 隐性：隐性源项通过下面的公式在有限控制体内进行积分并线性化：
∫𝑉 𝜌𝜙 𝑑𝑉 = 𝜌𝑃 𝑉𝑃 𝜙𝑃 (2.29)
 显性/隐性：隐性源项会改变矩阵的系数。根据源项的正负以及矩阵系数，隐形源
项会使得矩阵的对角占优有所改变。如果降低矩阵的对角占优会在矩阵迭代求解的
时候引入不稳定性。因此，OpenFOAM 提供了混合的源项离散形式：当源项小于 0
的时候进行隐形离散，源项大于 0 的时候进行显性离散8。从数学上来讲，节点 P
的矩阵系数是𝑉𝑝 max(𝜌𝑝 , 0)，源项为：𝑉𝑝 min(𝜌𝑝 , 0)；
2.4.10 其他显性离散
我们还可以使用其他的离散方法在 volField<Type>和 surface<Type>Fields 之间互相转换：
 面积分：fvc::surfaceIntegrate 其将面上的 surface<Type>Field 进行加和，然后除以网

格单元体积：(∑𝑓 𝜙𝑓 )𝑉𝑝 ，其返回一个 volField<Type>；
 面加和：fvc::surfaceSum 其将面上的 surface<Type>Field 进行加和：(∑𝑓 𝜙𝑓 )，其返

回一个 volField<Type>；
8 原文有误
38
 平均：fvc::average 计算 surface<Type>Field 的面平均值：(∑𝑓 𝑆𝑓 𝜙𝑓 )/ ∑𝑓 𝑆𝑓 ，其返回

一个 volField<Type>；
 面插值：faceInterpolate()将 volField<Type>从网格单元中心的值插值到面上并返回
一个 surface<Type>Field（使用中心差分）；
2.5 时间离散
虽然我们在 2.4.3 节和 2.4.4 节讨论了时间项的离散，我们仍然需要讨论一下瞬态问题下的

时间项处理问题。首先我们使用𝒜表示任意的张量算符（例如拉普拉斯），我们假定有下面
的积分形式的偏微分方程：
𝑡+Δ𝑡 𝜕
∫𝑡 [𝜕𝑡 ∫𝑉 𝜌𝜙 𝑑𝑉 + ∫𝑉 𝒜𝜙 𝑑𝑉 ] 𝑑𝑡 = 0 (2.30)
参考公式 2.21 使用欧拉隐形，第一项可以表示为：
𝑡+Δ𝑡 𝜕 𝑡+Δ𝑡 (𝜌𝑃 𝜙𝑃 𝑉)𝑛 −(𝜌𝑃𝜙𝑃 𝑉)𝑜 (𝜌𝑃𝜙𝑃 𝑉)𝑛 −(𝜌𝑃 𝜙𝑃𝑉)𝑜
∫𝑡 [𝜕𝑡 ∫𝑉 𝜌𝜙 𝑑𝑉 ] 𝑑𝑡 = ∫𝑡 ∆𝑡
𝑑𝑡 = ∆𝑡
𝑑𝑡 (2.31)
第二项可以表示为：
𝑡+Δ𝑡 𝑡+Δ𝑡
∫𝑡 [∫𝑉 𝒜𝜙 𝑑𝑉 ] 𝑑𝑡 = ∫𝑡 𝒜∗ 𝜙 𝑑𝑡 (2.32)
其中𝒜∗表示𝒜的离散形式，因此，时间项的积分可以采用以下三种方式：
 欧拉隐性：采用隐形形式离散空间项，因此在计算中取新值𝜙 n：
𝑡+Δ𝑡
∫𝑡 𝒜∗ 𝜙 𝑑𝑡 = 𝒜∗ 𝜙 𝑛 ∆𝑡 (2.33)
其为关于时间的一阶精度，因此为有界且无条件稳定的；
 显性：采用显形形式离散空间项，因此在计算中我们去旧值𝜙 0：
𝑡+Δ𝑡
∫𝑡 𝒜∗ 𝜙 𝑑𝑡 = 𝒜∗ 𝜙 𝑜 ∆𝑡 (2.34)
其为关于时间的一阶精度，如果库郎数大于 1，其会不稳定。库郎数这样定义：
𝐔 ⋅𝐝
𝐶𝑂 = |𝐝|𝑓2∆𝑡 (2.35)
其中𝐔f 为特征速度，例如流体速度等；
39
 C-N 格式：使用显性隐性的混合来离散空间项，因此需要对新值和旧值取平均：
𝑡+Δ𝑡 𝜙𝑛 +𝜙𝑜
∫𝑡 𝒜∗ 𝜙 𝑑𝑡 = 𝒜∗ ( 2
) ∆𝑡 (2.36)
其为关于时间的二阶精度，无条件稳定但是不一定有界；
2.5.1 OpenFOAM 中时间离散的代码
目前，对时间离散的处理要根据空间导数的处理方法（显性隐性）来选定。例如，我们打算
求解下面这个瞬态扩散方程：
𝜕𝜙
𝜕𝑡
= 𝜅∇2 𝜙 (2.37)
欧拉隐性代码如下：
solve(fvm::ddt(phi) == kappa*fvm::laplacian(phi))
其中我们使用 fvm 类将拉普拉斯项隐性离散，如果显性离散的话即为下述：
solve(fvm::ddt(phi) == kappa*fvc::laplacian(phi))
其中我们使用 fvc 类将拉普拉斯项显性离散，C-N 格式的离散可以这样植入：
solve
(
fvm::ddt(phi)
==
kappa*0.5*(fvm::laplacian(phi) + fvc::laplacian(phi))
)
2.6 边界条件
边界条件对于我们需要计算的问题是非常重要的。因此我们在所有的边界面指定边界条件。
边界条件分为两种类型：
 Dirichlet：描述边界上的变量值，因此我们称为固定值边界条件；
 Neumann：描述变量在边界上的法向梯度值，因此我们称之为固定梯度边界条件；
当我们在进行离散的时候，有时会包含这一项：∑，因此我们需要考虑如果某个网格单元面
为边界面的时候会发生什么样的情况：
40
 Fixed value：我们在边界处指定固定值：
 使用边界值𝜙b 代替𝜙f ，可参见方程 2.16；
 如果需要计算面梯度(∇𝜙)𝑓 ，例如拉普拉斯项，我们可以这样用公式表示：
𝜙𝑏 −𝜙𝑃
𝐒𝑓 ⋅ (∇𝜙)𝑓 = |𝑆𝑓 | |𝐝|
(2.38)
 Fixed gradient：固定梯度边界条件需要计算梯度和边界面法向的内积：
𝐒
𝔤𝑏 = (|𝐒| ⋅ ∇𝜙 ) (2.39)
𝑓
 当需要边界的𝜙f 值的时候，我们需要从网格单元中心值进行插值：
𝜙𝑓 = 𝜙𝑃 + 𝐝 ⋅ (∇𝜙)𝑓 = 𝜙𝑃 + |𝐝|𝔤𝑏 (2.40)
 当需要面法向梯度的时候，𝜙b 可以直接进行替换：
𝐒𝑓 ⋅ (∇𝜙)𝑓 = |𝑆𝑓 |𝔤𝑏 (2.41)
2.6.1 物理边界条件
工程师们通常通过实际的物理情况来指定边界条件。相比于前述的各种公式，真实的边界条
件非常具有描述性。在不可压缩缩流中我们有下面几个具有代表性的物理边界条件：
 进口：出口的速度场需要指定，压力应该指定法向梯度为零；
 出口 t：出口的压力场需要指定，速度应该指定法向梯度为零；
 无滑移壁面：壁面处的速度为 0，因此需要指定速度固定值边界条件。压力指定
为法向梯度为零（因为壁面出的通量为 0）
；
有时，计算域以及边界条件会呈对称分布，我们只需要计算其中的一半区域就可以了。其他
区域我们考虑为对称面。其中对称的交界面我们必须指定为：
 对称平面：对称平面边界条件指定平面的法向梯度为零；
41
第三章
算例演示
在这一章中，我们讨论几个 OpenFOAM 所提供的算例。OpenFOAM 对每一个求解器都提供

了对应的算例，有些算例在《OpenFOAM 用户指南》中第二章里进行了介绍。OpenFOAM
提供这些算例的目的在于帮助用户了解 OpenFOAM 求解器和其中的一些工具的特点，比如
前处理和后处理工具以及数值格式与数值算法。这些算例的运行结果也证明了各个求解器的
正确性。
每个算例都从如下几个方面来描述所要求解的问题：几何模型、初始条件、边界条件、
所解控制方程、求解过程中用到的模型和物理参数。这些例子中所选用的计算域有可能只是
原几何模型的一部分（比如当我们引入 symmetryPlane 边界的时候）。每个算例中都有对网
格生成方法的简明介绍（比如 blockMesh）。当然，用户可以直接查看网格，因为 OpenFOAM
当中每一个算例都含有 polyMesh 文件夹，并且在这个文件夹中有描述网格的数据文件。
这些例子都位于 OpenFOAM 安装路径下一个名为 tutorials 的子文件夹下。他们根据相
应的求解器被分配到各个子文件夹中。比如说，每一个 icoFoam 对应的例子都被存储于一个
叫 icoFoam 的子文件夹中。在运行算例之前，我们建议用户将算例复制到其他的文件夹之
中。我们推荐用户将 OpenFOAM 算例存储于$FOAM_RUN 文件夹中。如果这个文件夹还没
有被创建，用户可以用如下命令创建此文件夹：
mkdir –p $FOAM_RUN
教程中的算例可以用如下命令复制到此文件夹中：
cp –r $FOAM_TUTORIALS/* $FOAM_RUN
3.1 圆柱绕流
在这个例子当中，我们将利用 potentialFoam 来研究圆柱势流（potential flow）。这个例子将

会着重介绍 OpenFOAM 如何处理如下问题：
 非正交网格；
 对 OpenFOAM 当中的问题生成解析解；
42
3.1.1 问题阐述
我们将这个圆柱绕流问题描述如下：
计算域：计算域是二维的，由一个正方形和一个与其同心的圆组成，如图 3.1 所示
图 3.1 圆柱绕流几何模型示意图
控制方程：
不可压缩缩流体的质量守恒方程：
∇∙𝐔 = 0 (3.1)
 不可压缩缩无旋场的稳态压力方程：
∇2 𝑝 = 0 (3.2)
边界条件：
 入口（左边界）为固定速度入口 U = (1,0,0) m/s；

 出口（右边界）为固定压力出口 p = 0 pa；
 下边界为无滑移壁面；
 上边界为对称面边界；
初始条件：U = 0 m/s 且 p = 0 pa。虽然对于稳态问题来说，并不需要初始条件，但是对于

OpenFOAM 算例来说，无论求解问题是否稳态都必须给定初始条件；
求解器名称： potentialFoam，其为求解势流的求解器，假设流场不可压缩、无旋、无黏、稳
态、并且忽略重力；
算例名称：cylinder，这个算例位于$FOAM_TUTORIALS/potentialFoam 文件夹中；
43
3.1.2 potentialFoam 求解器的注意事项
potentialFoam 求解器对于 OpenFOAM 正确性的验证有着十分重要的作用，因为势流中的假

设使得许多几何模型比较简单的算例都具有解析解。对于这个圆柱绕流这个算例，我们可以
用已有的解析解来对照我们的数值解。potentialFoam 还可以对其他问题提供符合质量守恒
方程的初始场9。当处理类似问题的时候，此种方法可以避免由于初始场不稳定所带来的计
算不稳定性。总之，potentialFoam 可以用来从用户提供的不符合质量守恒的初始场生成符合
质量守恒的流场。
3.1.3 网格的生成
在《OpenFOAM 用户指南》中，我们已经阐述了如何用 blockMesh 生成网格。在这个算例

中，整个计算域由十个 block 构成，如图 3.2 所示。需要注意的是，在 OpenFOAM 之中，所
有网格都是按照三维网格处理的。如果我们想要解决二维问题，我们必须在第三个维度上加
上一个网格的厚度（虽然我们并没有这第三个维度上解控制方程）。
图 3.2 圆柱绕流中的 block
在图 3.2 中，我们只显示了几何图形的背面，即 z = -0.5 平面，在这个平面上，只有 0 到 18

号节点，而另外 19 个节点在几何图形的正面，z = 0.5 平面。正面 19 个节点的排序和背面相
同，网格节点和 block 的信息被储存在 blockMeshDict 文件夹当中，如下所示：
1 /*--------------------------------*- C++ -*-------------------------------------------------------------------------------*\
2 |========= | |
3 |\\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
4 | \\ / O peration | Version: 2.3.1 |
5 | \\ / A nd | Web: www.OpenFOAM.org |
6 | \\ / M anipulation | |
7 \*--------------------------------------------------------------------------------------------------------------------------*/
8 FoamFile
9 {
10 version 2.0;
11 format ascii;
12 class dictionary;
13 object blockMeshDict;
9 在算例运行之前首先运行势流求解器可以大大加快收敛且增加求解稳定性
44
14 }
15 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
16
17 convertToMeters 1;
18
19 vertices #codeStream
20 {
21 codeInclude
22 #{
23 #include "pointField.H"
24 #};
25
26 code
27 #{
28 pointField points(19);
29 points[0] = point(0.5, 0, -0.5);
30 points[1] = point(1, 0, -0.5);
31 points[2] = point(2, 0, -0.5);
32 points[3] = point(2, 0.707107, -0.5);
33 points[4] = point(0.707107, 0.707107, -0.5);
34 points[5] = point(0.353553, 0.353553, -0.5);
35 points[6] = point(2, 2, -0.5);
36 points[7] = point(0.707107, 2, -0.5);
37 points[8] = point(0, 2, -0.5);
38 points[9] = point(0, 1, -0.5);
39 points[10] = point(0, 0.5, -0.5);
40 points[11] = point(-0.5, 0, -0.5);
41 points[12] = point(-1, 0, -0.5);
42 points[13] = point(-2, 0, -0.5);
43 points[14] = point(-2, 0.707107, -0.5);
44 points[15] = point(-0.707107, 0.707107, -0.5);
45 points[16] = point(-0.353553, 0.353553, -0.5);
46 points[17] = point(-2, 2, -0.5);
47 points[18] = point(-0.707107, 2, -0.5);
48
49 // Duplicate z points
50 label sz = points.size();
51 points.setSize(2*sz);
52 for (label i = 0; i < sz; i++)
53 {
54 const point& pt = points[i];
55 points[i+sz] = point(pt.x(), pt.y(), -pt.z());
56 }
57
58 os << points;
59 #};
60 };
61
62
63 blocks
64 (
65 hex (5 4 9 10 24 23 28 29) (10 10 1) simpleGrading (1 1 1)
75 );
76
77 edges
78 (
79 arc 0 5 (0.469846 0.17101 -0.5)
80 arc 5 10 (0.17101 0.469846 -0.5)
81 arc 1 4 (0.939693 0.34202 -0.5)
82 arc 4 9 (0.34202 0.939693 -0.5)
83 arc 19 24 (0.469846 0.17101 0.5)
84 arc 24 29 (0.17101 0.469846 0.5)
85 arc 20 23 (0.939693 0.34202 0.5)
86 arc 23 28 (0.34202 0.939693 0.5)
87 arc 11 16 (-0.469846 0.17101 -0.5)
88 arc 16 10 (-0.17101 0.469846 -0.5)
89 arc 12 15 (-0.939693 0.34202 -0.5)
90 arc 15 9 (-0.34202 0.939693 -0.5)
91 arc 30 35 (-0.469846 0.17101 0.5)
92 arc 35 29 (-0.17101 0.469846 0.5)
93 arc 31 34 (-0.939693 0.34202 0.5)
45
94 arc 34 28 (-0.34202 0.939693 0.5)

95 );
96
97 boundary
98 (
99 down
100 {
101 type symmetryPlane;
102 faces
103 (
104 (0 1 20 19)
105 (1 2 21 20)
106 (12 11 30 31)
107 (13 12 31 32)
108 );
109 }
110 right
111 {
112 type patch;
113 faces
114 (
115 (2 3 22 21)
116 (3 6 25 22)
117 );
118 }
119 up
120 {
122 faces
123 (
124 (7 8 27 26)
125 (6 7 26 25)
126 (8 18 37 27)
127 (18 17 36 37)
128 );
129 }
130 left
131 {
132 type patch;
133 faces
134 (
135 (14 13 32 33)
136 (17 14 33 36)
137 );
138 }
139 cylinder
140 {
141 type symmetry;
142 faces
143 (
144 (10 5 24 29)
145 (5 0 19 24)
146 (16 10 29 35)
147 (11 16 35 30)
148 );
149 }
150 );
151
152 mergePatchPairs
153 (
154 );
155
156 // ************************************************************************* //
3.1.4 边界条件和初始条件
下面我们根据图 3.1 所示使用 FoamX10或者手动设置边界条件。左边界应该是入口（Inlet），

右边界应该是出口（Outlet），下边界和圆柱边界应该设置为 symmetryPlane。当我们选择上
边界的时候应当遵循以下规则：所选择的边界条件应当允许我们用理论解与数值解进行比
较，因此，我们假设计算域在 y 方向上无穷大。这个假设使得垂直于上边界的速度梯度特别
10 下文提及的所有 FoamX 均停止更新，可采用其他的 OpenFOAM GUI 处理

46
小。因此我们在上边界选择垂直于上边界速度梯度为 0 的边界条件，即 symmetryPlane。这

个边界条件使得理论解和数值解之间的比较变得更加合理。
3.1.5 运行算例
在这个算例中，我们不需要设置流体的性质，因为整个流场是不可压缩无旋无黏的。在 system
子文件夹中，controlDict 文件夹记录着算例运行所需要的参数。需要注意的是，因为我们要
处理的是稳态问题，所以这个算例只应当运行一个时间步：
1 /*--------------------------------*- C++ -*-------------------------------------------------------------------------------*\

2 |========= | |
7 \*--------------------------------------------------------------------------------------------------------------------------*/
8 FoamFile
9 {
10 version 2.0;
11 format ascii;
13 location "system";
14 object controlDict;
15 }
16 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
17
18 application potentialFoam;
19
20 startFrom startTime;
21
22 startTime 0;
23
24 stopAt endTime;
25
26 endTime 1;
27
28 deltaT 1;
29
30 writeControl timeStep;
31
32 writeInterval 1;
33
34 purgeWrite 0;
35
36 writeFormat ascii;
37
38 writePrecision 6;
39
40 writeCompression off;
41
42 timeFormat general;
43
44 timePrecision 6;
45
46 runTimeModifiable true;
47
48 functions
49 {
50 difference
51 {
52 // Load the library containing the 'coded' functionObject
53 functionObjectLibs ("libutilityFunctionObjects.so");
54 type coded;
55 // Name of on-the-fly generated functionObject
56 redirectType error;
57 code
58 #{
59 // Lookup U
60 Info<< "Looking up field U\n" << endl;
61 const volVectorField& U = mesh().lookupObject<volVectorField>("U");
62
63 Info<< "Reading inlet velocity uInfX\n" << endl;
47
64
65 scalar ULeft = 0.0;
66 label leftI = mesh().boundaryMesh().findPatchID("left");
67 const fvPatchVectorField& fvp = U.boundaryField()[leftI];
68 if (fvp.size())
69 {
70 ULeft = fvp[0].x();
71 }
72 reduce(ULeft, maxOp<scalar>());
73
74 dimensionedScalar uInfX
75 (
76 "uInfx",
77 dimensionSet(0, 1, -1, 0, 0),
78 ULeft
79 );
80
81 Info << "U at inlet = " << uInfX.value() << " m/s" << endl;
82
83
84 scalar magCylinder = 0.0;
85 label cylI = mesh().boundaryMesh().findPatchID("cylinder");
86 const fvPatchVectorField& cylFvp = mesh().C().boundaryField()[cylI];
87 if (cylFvp.size())
88 {
89 magCylinder = mag(cylFvp[0]);
90 }
91 reduce(magCylinder, maxOp<scalar>());
92
93 dimensionedScalar radius
94 (
95 "radius",
96 dimensionSet(0, 1, 0, 0, 0),
97 magCylinder
98 );
99
100 Info << "Cylinder radius = " << radius.value() << " m" << endl;
101
102 volVectorField UA
103 (
104 IOobject
105 (
106 "UA",
107 mesh().time().timeName(),
108 U.mesh(),
109 IOobject::NO_READ,
110 IOobject::AUTO_WRITE
111 ),
112 U
113 );
114
115 Info<< "\nEvaluating analytical solution" << endl;
116
117 const volVectorField& centres = UA.mesh().C();
118 volScalarField magCentres(mag(centres));
119 volScalarField theta(acos((centres & vector(1,0,0))/magCentres));
120
121 volVectorField cs2theta
122 (
123 cos(2*theta)*vector(1,0,0)
124 + sin(2*theta)*vector(0,1,0)
125 );
126
127 UA = uInfX*(dimensionedVector(vector(1,0,0))
128 - pow((radius/magCentres),2)*cs2theta);
129
130 // Force writing of UA (since time has not changed)
131 UA.write();
132
133 volScalarField error("error", mag(U-UA)/mag(UA));
134
135 Info<<"Writing relative error in U to " << error.objectPath()
136 << endl;
137
138 error.write();
139 #};
140 }
141 }
142
143
48
144 // ************************************************************************* //
potentialFoam通过迭代来求解压力方程，并且可以对压力方程中的拉普拉斯项进行非正交修
正（压力方程中的显式项会在正交修正迭代之后进行更新）。对压力方程的非正交修正迭代
次数由 controlDict 之中的 nNonOrthogonalCorrectors 来决定。在第一个例子中，我们将
nNonOrthogonalCorrectors设置为 0，也就是说，在求解过程中没有非正交修正迭代步（压力
方程只被求解一次）。与之对应的数值解如图3.3（a）所示（t = 1，当稳态求解结束之时）。
我们期待数值解能够像解析解（如图3.3（c）所示）一样，即穿过整个计算域的流线是光滑
的，但是在网格正交性不好的地方我们却能观察到明显的误差（比如在block 0，1，3的交界
处）。如果我们将nNonOrthogonalCorrectors设置为3，再运行一次算例，其数值解如图3.3（b）
所示。这个解的结果与理论解非常相近（流线非常光滑），因此并没有明显的由非正交性所
引起的误差。
(a) 无非正交修正
(b) 有非正交修正
(c) 解析解
图 3.3 势流中的流线
49
3.2 稳态后向台阶湍流模拟
在这个算例中，我们将会进行稳态后向台阶湍流模拟。算例的原型来自 Pitz 和 Daily 的实验

研究，我们将对算例所得到的数值解与实验值进行比较。这个算例会涉及如下特性：
 利用 blockMesh 当中的网格非均匀分布功能生成网格；
 稳态湍流；
3.2.1 问题阐述
稳态后向台阶湍流模拟定义如下：
计算域：计算域是二维的，由一个窄的入口、一个后向台阶、一个在出口处收缩的喷管组
成，如图 3.4 所示：
图 3.4 后向台阶的几何模型
控制方程：
 不可压缩缩流体的质量守恒方程：
∇∙𝐔 = 0 (3.3)
 稳态动量方程：
∇ ∙ (𝐔𝐔) + ∇ ∙ 𝐑 = −∇𝑝 (3.4)
这里的 p 是动压，𝐑 = ν𝑒𝑓𝑓 ∇𝑼是带有有效动力粘度系数ν𝑒𝑓𝑓 的粘性力项。ν𝑒𝑓𝑓 是通
过传递模型和湍流模型计算而来；
初始条件：U = 0 m/s 且 p = 0 pa。虽然对于稳态问题来说，并不需要初始条件，但是对于

OpenFOAM 算例来说，必须给定初始条件，无论求解问题是否稳态；
边界条件：
 进口为固定速度入口，U = (10, 0, 0)m/s；

 出口为固定压力出口 p = 0 pa；
 其他边界都为非滑移边界；
50
物性：
𝜇
 空气的运动粘度系数𝜈 = 𝜌 = 18.1 × 10−6/1.293 = 14.0𝜇𝑚2 /𝑠 ；
湍流模型：
 标准𝑘 − 𝜀模型；
 常数：𝐶𝜇 = 0.09; 𝐶1 = 1.44; 𝐶2 = 1.92; 𝛼𝑘 = 1; 𝛼𝜖 = 0.76923；
求解器名称：simpleFoam，其为稳态不可压缩流体问题的求解器；
算例名称：pitzDaily，位于$FOAM_TUTORIALS/simpleFoam 文件夹中；
我们利用 simpleFoam 来求解这个问题。simpleFoam 利用 SIMPLE 格式来求解稳态不可压缩

流体的问题。这个求解器可以使用 incompressibleTurbulenceModel 库中所有的湍流模型以及
incompressibleTransportModel 库中的所有牛顿以及非牛顿流体模型。
我么希望这个算例的流场足够的复杂，因此需要利用网格非均匀分布才能得到最优解。一般
来说，剪切应力（速度梯度）大的区域网格细化尤其重要。在这种区域的网格需要比其他剪
切应力较低的区域更细。在进行求解之前，我们可以提前估计一下在哪些地方容易出现大的
剪切应力。在入口处，我们有一个沿着 x 方向的稳态流动，在流体经过台阶时，这个稳态流
动和下方静止的流体之间会产生一个剪切应力，并且在计算域的下半部分产生漩涡。所以，
在这个问题中，剪切应力比较大的区域应该在计算域的中心线附近和靠近壁面的地方。
计算域被分为 12 个 block，如图 3.5 所示。其网格是三维的。图 3.5 所示的是计算域的
背面，z = - 0.5 平面。网格节点和 block 的信息储存在 blockMeshDict 文件夹中，如下所示：
1 /*--------------------------------*- C++ -*-------------------------------------------------------------------------------*\

2 |========= | |
7 \*--------------------------------------------------------------------------------------------------------------------------*/
8 FoamFile
9 {
10 version 2.0;
11 format ascii;
14 }
15 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
16
17 convertToMeters 0.001;
18
19 vertices
20 (
21 (-20.6 0 -0.5)
22 (-20.6 3 -0.5)
23 (-20.6 12.7 -0.5)
24 (-20.6 25.4 -0.5)
25 (0 -25.4 -0.5)
51
图 3.5 后向台阶问题中的 block

26 (0 -5 -0.5)
27 (0 0 -0.5)
28 (0 3 -0.5)
29 (0 12.7 -0.5)
30 (0 25.4 -0.5)
31 (206 -25.4 -0.5)
32 (206 -8.5 -0.5)
33 (206 0 -0.5)
34 (206 6.5 -0.5)
35 (206 17 -0.5)
36 (206 25.4 -0.5)
37 (290 -16.6 -0.5)
38 (290 -6.3 -0.5)
39 (290 0 -0.5)
40 (290 4.5 -0.5)
41 (290 11 -0.5)
42 (290 16.6 -0.5)
43 (-20.6 0 0.5)
44 (-20.6 3 0.5)
45 (-20.6 12.7 0.5)
46 (-20.6 25.4 0.5)
47 (0 -25.4 0.5)
48 (0 -5 0.5)
49 (0 0 0.5)
50 (0 3 0.5)
51 (0 12.7 0.5)
52 (0 25.4 0.5)
53 (206 -25.4 0.5)
54 (206 -8.5 0.5)
55 (206 0 0.5)
56 (206 6.5 0.5)
57 (206 17 0.5)
58 (206 25.4 0.5)
59 (290 -16.6 0.5)
60 (290 -6.3 0.5)
61 (290 0 0.5)
62 (290 4.5 0.5)
63 (290 11 0.5)
64 (290 16.6 0.5)
65 );
66
67 blocks
68 (
69 hex (0 6 7 1 22 28 29 23) (18 7 1) simpleGrading (0.5 1.8 1)
70 hex (1 7 8 2 23 29 30 24) (18 10 1) simpleGrading (0.5 4 1)
73 hex (5 11 12 6 27 33 34 28) (180 9 1) edgeGrading (4 4 4 4 0.5 1 1 0.5 1 1 1 1)
74 hex (6 12 13 7 28 34 35 29) (180 7 1) edgeGrading (4 4 4 4 1.8 1 1 1.8 1 1 1 1)
75 hex (7 13 14 8 29 35 36 30) (180 10 1) edgeGrading (4 4 4 4 4 1 1 4 1 1 1 1)
76 hex (8 14 15 9 30 36 37 31) (180 13 1) simpleGrading (4 0.25 1)
82 );
83
84 edges
85 (
52
86 );
87
88 boundary
89 (
90 inlet
91 {
92 type patch;
93 faces
94 (
95 (0 22 23 1)
96 (1 23 24 2)
97 (2 24 25 3)
98 );
99 }
100 outlet
101 {
102 type patch;
103 faces
104 (
105 (16 17 39 38)
106 (17 18 40 39)
107 (18 19 41 40)
108 (19 20 42 41)
109 (20 21 43 42)
110 );
111 }
112 upperWall
113 {
114 type wall;
115 faces
116 (
117 (3 25 31 9)
118 (9 31 37 15)
119 (15 37 43 21)
120 );
121 }
122 lowerWall
123 {
124 type wall;
125 faces
126 (
127 (0 6 28 22)
128 (6 5 27 28)
129 (5 4 26 27)
130 (4 10 32 26)
131 (10 16 38 32)
132 );
133 }
134 frontAndBack
135 {
136 type empty;
137 faces
138 (
139 (22 28 29 23)
140 (23 29 30 24)
141 (24 30 31 25)
142 (26 32 33 27)
143 (27 33 34 28)
144 (28 34 35 29)
145 (29 35 36 30)
146 (30 36 37 31)
147 (32 38 39 33)
148 (33 39 40 34)
149 (34 40 41 35)
150 (35 41 42 36)
151 (36 42 43 37)
152 (0 1 7 6)
153 (1 2 8 7)
154 (2 3 9 8)
155 (4 5 11 10)
156 (5 6 12 11)
157 (6 7 13 12)
158 (7 8 14 13)
159 (8 9 15 14)
160 (10 11 17 16)
161 (11 12 18 17)
162 (12 13 19 18)
163 (13 14 20 19)
164 (14 15 21 20)
165 );
53
166 }
167 );
168
169 mergePatchPairs
170 (
171 );
172
173 // ************************************************************************* //
这个算例的一大特点，是它大量的使用了《OpenFOAM 用户指南》中 5.3.1 所介绍的 blockMesh

中的网格非均匀功能。用户可以看到，在 block 4、5、6 中，每个边均使用了非均匀分布（每
个膨胀率都对应着 block 当中的一个边）。前四个膨胀率对应与 x1 轴平行的四个边，中间四
个膨胀率对应与 x2 轴平行的四个边，后四个膨胀率对应与 x3 轴平行的四个边。在 block 4、
5、6 中，膨胀率在 x1 和 x3 方向上都是相同的，但是在 x2 方向（y 方向）却是不同的。如果
我们参考《OpenFOAM 用户指南》5.3.1 中对 block 的定义，我们可以发现 block 4，5，6 中
的左边和右边的细化是不一样的。这样做的目的是为了在最关键的地方（台阶的尖角处）生
成细网格。
网格可以在终端由 bockMesh 生成，也可以在 FoamX 之中生成，查看网格的方法再次
不再赘述。
3.2.3 边界条件和初始场
与算例相关的参数文件可以在 FoamX 之中浏览或者修改（当然也可以用手动修改）。在这个

算例中，我们需要为如下参数设置边界条件和初始条件：速度 U，压力 p，湍流动能𝑘，和
耗散率𝜀。边界条件的设置可以通过在 FoamX 之中设置 patch 的类型来实现: 上边界和下边
界被设置为 wall，左边界被设置为 Inlet，右边界被设置为 Outlet。这样的边界条件组合需要
我们在入口处设定 U、𝑘、𝜀的类型为 fixedValue。U 的具体值已经在问题阐述中给出。k 和𝜀
的值应该利用和《OpenFOAM 用户指南》中 2.1.8.1 节一样的方法计算出来。我们假设入口
湍流是各向同性，而且入口处的速度波动为 5%有：
5
𝑈𝑥′ = 𝑈𝑦′ = 𝑈𝑧′ = 10 = 0.5m/s (3.5)
100
且
𝑘 = 1.5(0.5)2 = 0.375𝑚2/𝑠 2 (3.6)
如果我们假定湍流尺度为入口边界长度的十分之一有：
𝐶𝜇0.75 𝑘1.5 0.090.75 0.3751.5
ε= 𝑙
= 0.1×25.4×10−3
= 14855𝑚2 /𝑠 2 (3.7)
在出口处，我们只需要将 p 设置为 0 pa。
3.2.4 求解器的控制参数
fvSchemes 的设置如下：timeScheme 设置为 SteadyState；gradScheme 和 laplacianScheme 设

置为默认值 Gauss；为了确保有界性，divScheme 应该设置为 UD。
我们需要特殊注意 fvTolerances 的设置。虽然顶层的 simpleFoam 代码只包含了压力和
速度的方程，但是在湍流模型中，我们还需要求解𝑘，𝜀和 R 的方程，所以，我们需要为 5 个
方程设置残差（tolerance）。
除了将压力的 solverTolerance 设置为 10-6 且 solverRelativeTolerance
设置为 0.01 之外，对于其他变量，我们将这两个参数分别设置为 10-5 和 0.1。由于这个算例
是一个稳态问题，所以我们需要设置松弛因子（relaxationFactor）。对于压力，为了避免计算
的不稳定，松弛因子被设置为 0.3。对于其他参数（U，𝑘，𝜀，R），松弛因子可以设置为 0.7。
最后，在 controlDict 之中，我们将时间步长设置为 1(因为计算的稳态的)，这里的时间
54
步长只用于稳态计算中的迭代次数进行计数。这个算例需要至少一千次迭代才能收敛，所以，
我们将 endTime 设置为 1000。我们还需要保证存储足够的时间步且不会太多，这样可以设
置为 50 步一写，以保证用户的硬盘空间不会在算例运行过程中被数据填满。
3.2.5 运行算例和后处理
图 3.6 后面台阶问题中的漩涡发展
接下来我们运行程序并且对运行结果进行后处理。在几步迭代之后，比如 50 步，在台阶的
尖角下方会出现一个漩涡，但是漩涡在 x 方向上长度很短，如图 3.6（a）所示。在 1000 个
时间步之后，漩涡沿着 x 方向从台阶处一直延伸到出口处，此时系统已经达到稳态，并且漩
涡已经充分发展，如图 3.6（b-c）所示。
3.3 前向台阶的超音速绕流
在这个算例中，我们将会探究前向台阶的超音速绕流。在这个算例中，入口速度马赫数为 3，
在入口附近有一个台阶，在台阶附近会产生激波。
这个算例主要涉及 OpenFOAM 的如下特性：
 超音速问题的求解
3.3.1 问题阐述
计算域：计算域是二维的，并且由一个入口区域和一个前向台阶组成。台阶的高度占入口长
55
度的 20%，如图 3.7 所示：
图 3.7 前向台阶的几何模型
控制方程：
 质量守恒方程：
∂𝜌
∂t
+ ∇ ∙ (𝜌𝐔) = 0 (3.8)
 理想气体状态方程：
𝑝 = 𝜌𝑅𝑇 (3.9)
 牛顿流体的动量守恒方程：
∂𝜌𝐔
∂t
+ ∇ ∙ (𝜌𝐔𝐔) − ∇ ∙ 𝜇∇𝐔 = −∇𝑝 (3.10)
 ，e = 𝐶𝑣 𝑇，和傅里叶定律，q = −k∇T：
能量方程（忽略粘性项）
∂𝜌e 𝑘
+ ∇ ∙ (𝜌𝐔e) − ∇ ∙ ∇e = −𝑝∇ ∙ 𝐔 (3.11)
∂t 𝐶𝑣
初始条件：𝑈 = 0 m/s, 𝑝 = 1 pa, 𝑇 = 1k；
边界条件：
 入口设置为固定值边界条件 fixedValue，𝑈 = 3 m/s = Mach 3, 𝑝 = 1 pa, 𝑇 = 1k；

 出口对速度，压力，温度（U，p，T）设置为法相梯度为零 zeroGradient；
 下边界被设置为无滑移绝热边界；
 上边界被设置为 symmetryPlane；
物性：
 空气的动力粘度系数为𝜇 = 18.1μPas；
热力学性质：
 定压比热容 𝐶𝑣 = 1.78571 J/kgK；

 气体常数 𝑅 = 0.714286 J/kgK；
56
 热导率 𝑘 = 32.3 μW/mK ；
算例名称：forwardStep，其位于$FOAM_TUTORIALS/sonicFoam 文件夹中；
求解器名称：sonicFoam：为可压缩流体的超音速和跨音速问题设计的求解器；
在本算例中，音速被设计为 1m/s (c = √𝛾𝑅𝑇 = 1𝑚/𝑠)。这样做的结果就是速度的大小直接等

于马赫数。音速的大小可以利用如下关系计算。对于理想气体，𝐶𝑝 − 𝐶𝑣 = 𝑅：
𝐶𝑝 R
γ= 𝐶𝑣
= 𝐶 +1 (3.12)
𝑣
3.3.2 网格生成
本算例中所应用的网格比较简单，用大小相同的长方体网格生成即可。这些长方体网格在 x
方向长 0.06m，y 方向长 0.05m。整个计算域可以被分为三个 block，一个在台阶上沿下方，
两个在台阶上沿上方（分列于台阶两侧）。关于网格的信息被储存在如下所示的
blockMeshDict 文件之中：
1 /*--------------------------------*- C++ -*-------------------------------------------------------------------------------*\

2 |========= | |
7 \*--------------------------------------------------------------------------------------------------------------------------*/
8 FoamFile
9 {
10 version 2.0;
11 format ascii;
14 }
15 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
16
18
19 vertices
20 (
21 (0 0 -0.05)
22 (0.6 0 -0.05)
23 (0 0.2 -0.05)
24 (0.6 0.2 -0.05)
25 (3 0.2 -0.05)
26 (0 1 -0.05)
27 (0.6 1 -0.05)
28 (3 1 -0.05)
29 (0 0 0.05)
30 (0.6 0 0.05)
31 (0 0.2 0.05)
32 (0.6 0.2 0.05)
33 (3 0.2 0.05)
34 (0 1 0.05)
35 (0.6 1 0.05)
36 (3 1 0.05)
37 );
38
39 blocks
40 (
44 );
45
46 edges
47 (
57
48 );
49
50 boundary
51 (
52 inlet
53 {
54 type patch;
55 faces
56 (
57 (0 8 10 2)
58 (2 10 13 5)
59 );
60 }
61 outlet
62 {
63 type patch;
64 faces
65 (
66 (4 7 15 12)
67 );
68 }
69 bottom
70 {
72 faces
73 (
74 (0 1 9 8)
75 );
76 }
77 top
78 {
80 faces
81 (
82 (5 13 14 6)
83 (6 14 15 7)
84 );
85 }
86 obstacle
87 {
88 type patch;
89 faces
90 (
91 (1 3 11 9)
92 (3 4 12 11)
93 );
94 }
95 );
96
97 mergePatchPairs
98 (
99 );
100
101 // ************************************************************************* //
3.3.3 运行算例
本算例将会在 5s 之后达到稳态。10 秒钟时的压力分布如图 3.8 所示。计算结果显示出了压

力的明显不连续性，也就是激波。我们可以看到，激波由从台阶下沿的前方开始发展。
3.3.4 练习
用户可以检验入口速度的升高对于解的影响。
58
图 3.8 前向台阶问题中的激波
3.4 高压水箱放水问题
在这个算例中，我们将会探究一个装满水的高压水箱放水问题。这个问题的一大特点是压力
波的传播，因此，这个算例必须要考虑流体的可压缩性。
这个算例涉及到了 OpenFOAM 的如下特性：
 网格加密；
 流体中的压力波；
3.4.1 问题阐述
计算域：计算域是二维的，由一个水箱和一小截泄流管组成，如图 3.9 所示：
图 3.9 带有泄流管道的水箱示意图
控制方程：在这个问题中，为了能够准确捕捉以有限速度传播的压力波，我们需要考虑流体
的压缩性𝜓。我们将利用正压关系（barotropic relationship）来将压力、密度和压缩性𝜓联系
在一起：
 质量守恒方程：
∂𝜌
+ ∇ ∙ (𝜌𝐔) = 0 (3.13)
∂t
59
 正压关系：
𝜕𝜌 𝜌
= =𝜓 (3.14)
𝜕𝑝 𝐾
K是体积模量；
 公式 3.14 被线性化为：
𝜌 ≈ 𝜌0 + 𝜓(𝑝 − 𝑝0 ) (3.15)
𝜌0 和𝑝0 分别为参考密度和参考压力，𝜌(𝑝0 ) = 𝜌0；
 牛顿流体的动量守恒方程：
∂𝜌𝐔
∂t
+ ∇ ∙ (𝜌𝐔𝐔) − ∇ ∙ 𝜇∇𝐔 = −∇𝑝 (3.16)
边界条件：使用 FoamX，将边界条件设置如下：
 outerWall 设置为 wall；

 axis 设置为 symmetryPlane；
 outlet 设置为 pressureOutlet，出口压力位稳态 p=0 bar；
 front 和 back 设置为 empty；
初始条件：U = 0m/s, 𝑝 = 100 bar；
物性：
 水的动力粘度系数为𝜇 = 1.0mPas；
热力学性质：
 水的密度：𝜌 = 1000kg/m3；
 参考压力：𝑝0 = 1 𝑏𝑎𝑟；
 水的压缩性：𝜓 = 4.54 × 10−7 𝑠 2 /𝑚2 ；
求解器名称：sonicLiquidFoam，用来求解可压缩超音速层流问题；
算例名称：decompressionTank，这个算例位于$FOAM_TUTORIALS/sonicLiquidFoam 文件夹
中；
下面为其 blockMeshDict 字典文件：
1 /*--------------------------------*- C++ -*-------------------------------------------------------------------------------*\

2 |========= | |
60
7 \*--------------------------------------------------------------------------------------------------------------------------*/
8 FoamFile
9 {
10 version 2.0;
11 format ascii;
14 }
15 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
16
17 convertToMeters 0.1;
18
19 vertices
20 (
21 (0 0 -0.1)
22 (1 0 -0.1)
23 (0 0.5 -0.1)
24 (1 0.5 -0.1)
25 (1.5 0.5 -0.1)
26 (0 0.6 -0.1)
27 (1 0.6 -0.1)
28 (1.5 0.6 -0.1)
29 (0 3 -0.1)
30 (1 3 -0.1)
31 (0 0 0.1)
32 (1 0 0.1)
33 (0 0.5 0.1)
34 (1 0.5 0.1)
35 (1.5 0.5 0.1)
36 (0 0.6 0.1)
37 (1 0.6 0.1)
38 (1.5 0.6 0.1)
39 (0 3 0.1)
40 (1 3 0.1)
41 );
42
43 blocks
44 (
49 );
50
51 edges
52 (
53 );
54
55 boundary
56 (
57 outerWall
58 {
59 type wall;
60 faces
61 (
62 (0 1 11 10)
63 (1 3 13 11)
64 (3 4 14 13)
65 (7 6 16 17)
66 (6 9 19 16)
67 (9 8 18 19)
68 );
69 }
70 axis
71 {
73 faces
74 (
75 (0 10 12 2)
76 (2 12 15 5)
77 (5 15 18 8)
78 );
79 }
80 nozzle
81 {
82 type patch;
83 faces
84 (
85 (4 7 17 14)
86 );
61
87 }
88 back
89 {
90 type empty;
91 faces
92 (
93 (0 2 3 1)
94 (2 5 6 3)
95 (3 6 7 4)
96 (5 8 9 6)
97 );
98 }
99 front
100 {
101 type empty;
102 faces
103 (
104 (10 11 13 12)
105 (12 13 16 15)
106 (13 14 17 16)
107 (15 16 19 18)
108 );
109 }
110 );
111
112 mergePatchPairs
113 (
114 );
115
116 // ************************************************************************* //
为了提高计算的准确性，我们将参考压力设置为 1 bar。需要注意的是，初始场和边界条件
中的压力都应该被参考压力值抵消一部分。
3.4.3 准备运行算例
在我们设置 controlDict 的运行参数之前，我们应该考虑一下本算例的特征速度。在本算例

中，流体的速度是很小的，但是压力波却在水中以声速传播。水中的声速计算如下：
1 1
c = √𝜓 = √4.54×10−7 = 1483.2𝑚/𝑠 (3.17)
对于上述网格，特征网格长度越为 2mm(请注意在 blockMeshDict 之中我们把 convertToMeters

被设置为 0.1)。
𝑈∆𝑡
Co = ∆𝑥
(3.18)
本算例合理的时间步应该为∆t = 5 × 10−7 𝑠，在此时间步长下，取水中的声速计算库郎数为

0.35。需要注意的是，因为我们想要捕捉的是压力波在水中的传播，所以程序在运行时所输
出的 Co 数可能会比 0.35 小两个数量级。controlDict 之中的设置如下：
1 /*--------------------------------*- C++ -*-------------------------------------------------------------------------------*\

2 |========= | |
7 \*--------------------------------------------------------------------------------------------------------------------------*/
8 FoamFile
9 {
10 version 2.0;
11 format ascii;
13 location "system";
14 object controlDict;
15 }
16 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
17
18 application sonicLiquidFoam;
62
19
20 startFrom startTime;
21
22 startTime 0;
23
24 stopAt endTime;
25
26 endTime 0.0001;
27
28 deltaT 5e-07;
29
30 writeControl timeStep;
31
32 writeInterval 20;
33
34 purgeWrite 0;
35
36 writeFormat ascii;
37
38 writePrecision 6;
39
40 writeCompression off;
41
42 timeFormat general;
43
44 timePrecision 6;
45
46 runTimeModifiable true;
47
48
49 // ************************************************************************* //
3.4.4 运行算例
图 3.10 压力波的传播
用户可以运行算例，并且在 paraFoam 中观察运行结果。通过喷管流出的液体使得计算域中

形成沿着喷管移动的压力波。当压力波到达喷管和水箱的连接处时，有的压力波通过连接处
传入水箱中，有的被反射回喷管中。由于压力波在喷管中被来回反弹，通过连接处的压力波
会膨胀至整个水箱。参见图 3.10 中的压力云图，这样做的目的是为了比从单纯的等值线更
63
清晰地观察压力波。
如果我们让算例运行足够长的时间，让压力波充分的反射回喷管，我们可以清楚地在喷
管中观察到绝对负压。这个现象非常合理且有物理依据，因为流体能够承受张力（比如说负
压）。然而，在实际情况下，由于液体的非纯净性，液体中为溶解的气体对流体流动会产生
一定的空化效果（由于负压产生的气泡，沸腾和气化的效果一样）。因此，在实际情况下，
我们不会观察到压力下降到气化压力之下（至少在产生气泡之前）。
3.4.5 细化网格来获得更精确解
图 3.11 更加细化的网格中的压力波传播
通过观察压力场随时间的变化，我们可以看到压力波传入水箱中，并且在各个壁面之间进行
反射。我们也可以清楚的看到压力波在某些网格区域有所聚集。现在，为了更好的捕捉压力
波，我们对网格进行细化，并且减小时间步长。我们只需要将 blockMeshDict 里面的网格数
乘以 4，比如 block 0 的网格数从（30，20，1）增加到（120，80，1）。在文件目录下运行
除此之外，
blockMesh。为了将 Co 数控制在 1 以内，我们也应该将时间步长缩小到∆t = 10−7 𝑠。
如图 3.11 所示，加密之后的网格明显能够得到更精确的结果。
3.5 磁流体问题
在这个算例中，我们将会探究流过磁场的导电流体。这类问题被称为磁流体问题，我们用
mhdFoam 来解决类似问题。
64
3.5.1 问题阐述
这个算例被称为 Hartmann 问题，这个问题有解析解，我们可以用解析解来验证 mhdFoam 的

正确性。这个问题的定义如下：
计算域：计算域是二维的，由两个平行的平板组成，流体从两个平板之间流过，如图 3.12 所
示：
图 3.12 Hartmann 问题的几何模型
控制方程：
 不可压缩流体的质量守恒方程：
∇∙𝐔=0 (3.19)
 不可压缩缩流体的动量守恒方程：
∂𝐔
+ ∇ ∙ (𝑼𝑼) + ∇ ∙ (2𝑩Γ𝐵𝑈 𝑩) + ∇ ∙ (𝜈𝑼) + ∇(Γ𝐵𝑈 𝑩: 𝑩) = −∇𝑝 (3.20)
∂t
B 为磁感应强度，Γ𝐵𝑈 = (2𝜇𝜌)−1；
 麦克斯韦方程组：
∂𝐁
∇ × 𝐄 = − ∂t (3.21)
E 为电场强度；
∇∙𝐁 =0 (3.22)
𝜕𝑫
∇×𝑯= 𝑱+ 𝜕𝑡
=𝑱 (3.23)
𝜕𝑫
假设 𝜕𝑡 ≪ 𝑱，这里 H 为磁场强度，J 为电流密度，D 为电通量密度；
 电流连续性：
∇∙𝐉 =0 (3.24)
65
 本构方程：
𝐁 = 𝜇𝐇 (3.25)
 欧姆定律：
𝐉 = σ(𝐄 + 𝐔 × 𝐁) (3.26)
 将方程 3.21，3.23，3.26 联立，并且取旋度运算，可以得到：

∂𝐁
∂t
+ ∇ ∙ (𝑼𝑩) − ∇ ∙ (𝜙𝐵 𝑼) − ∇ ∙ (Γ𝐵𝑈 𝑩) = 0 (3.27)
边界条件：
 入口条件被设置为固定值速度 U = (1,0,0) m/s；

 出口条件被设置为稳态压力出口 p = 0 pa；
 上边界被设置为壁面，B = (0,20,0)T；
 下边界被设置为壁面，B = (0,20,0)T；
 前后面边界被设置为 empty；
初始条件：U = 0 m/s，p = 100 pa，B = (0,20,0)T；
物性：
 运动粘度系数：𝜈 = 1m2 /𝑠；

 密度：𝜌 = 1kg/m3 ；
 电导率：𝜎 = 1(Ωm)−1；
 磁导率：𝜇 = 1H/m；
求解器名称：mhdFoam：用于求解不可压缩层流磁流体问题的求解器；
算例名称：hartmann，这个算例位于$FOAM_TUTORIALS/mhdFoam 文件夹中；
本算例的网格很简单，在 x 方向有 100 个网格，在 y 方向有 40 个网格。blockMeshDict 文件

如下：
1 /*--------------------------------*- C++ -*-------------------------------------------------------------------------------*\
2 |========= | |
7 \*--------------------------------------------------------------------------------------------------------------------------*/
8 FoamFile
9 {
10 version 2.0;
11 format ascii;
14 }
15 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
16
66
18
19 vertices
20 (
21 (0 -1 0)
22 (20 -1 0)
23 (20 1 0)
24 (0 1 0)
25 (0 -1 0.1)
26 (20 -1 0.1)
27 (20 1 0.1)
28 (0 1 0.1)
29 );
30
31 blocks
32 (
34 );
35
36 edges
37 (
38 );
39
40 boundary
41 (
42 inlet
43 {
44 type patch;
45 faces
46 (
47 (0 4 7 3)
48 );
49 }
50 outlet
51 {
52 type patch;
53 faces
54 (
55 (2 6 5 1)
56 );
57 }
58 lowerWall
59 {
60 type patch;
61 faces
62 (
63 (1 5 4 0)
64 );
65 }
66 upperWall
67 {
68 type patch;
69 faces
70 (
71 (3 7 6 2)
72 );
73 }
74 frontAndBack
75 {
76 type empty;
77 faces
78 (
79 (0 3 2 1)
80 (4 5 6 7)
81 );
82 }
83 );
84
85 mergePatchPairs
86 (
87 );
88
89 // ************************************************************************* //
3.5.3 运行算例
用户可以运行算例并且在 paraFoam 中查看运行结果。我们也可以用 Ucomponents 命令将向

67
量场 U 转化为 U 各个分量的标量场。磁流体问题与许多参数相关，这些参数构成了无量纲
的 Hartmann 数。它表征电磁力和粘性力的比值：
𝜎
𝑀 = 𝐵𝐿√𝜌𝜈 (3.28)
这里，L 代表问题的特征长度，𝐵𝑦 = 20𝑇, 𝑀 = 20。在这个算例中，电磁力远大于粘性力。

在计算结果中，流场在 2s 时达到稳态，此时，如果我们在整个计算域的中间，即 x = 10m 的
地方观察会发现，速度曲线几乎是平的。这条曲线可以通过在 paraFoam 中通过绘制 Ux 的
曲线来得到。现在，用户可以将 By 的值调到 1，重新运行算例并运行 Ucomponents。可以看
出，M = 1 的时候电磁力和粘性力大小相近。
图 3.13 在 By=1 和 By=10 时 Hartmann 问题的速度曲线
这时，速度曲线看起来像泊肃叶流中的双曲线，如图 3.13 所示。为了验证数值解的正确性，

解析解也被添加到图 3.13 中以作比较。解析解的表达式如下：
𝑦
𝑈𝑥 (𝑦) cosh 𝑀−cosh 𝑀( )
𝑈𝑥 (0)
= cosh 𝑀−1
𝐿
(3.29)
这里，特征长度 L 为计算域宽度的一半，即 1m。

68
索引
操作符数 A B C D E F G H I J K L M N O P Q R S T U V W X Z
access functions, P-19

操作符 addLayersControls keyword, U-146
* adiabaticFlameT utility, U-95
tensor member function, P-21
+ adiabaticPerfectFluid model, U-99, U-179
tensor member function, P-21 adjointShapeOptimizationFoam solver, U-
- 83
tensor member function, P-21 adjustableRunTime
/ keyword entry, U-60, U-111
tensor member function, P-21 adjustTimeStep keyword, U-60, U-112
/*...*/ agglomerator keyword, U-123
C++ syntax, U-76 algorithms tools, U-96
// alphaContactAngle
C++ syntax, U-76 boundary condition, U-57
OpenFOAM file syntax, U-104 analytical solution, P-41
# include Animations window panel, U-168
C++ syntax, U-70, U-76 anisotropicFilter model, U-100
& Annotation window panel, U-24, U-167
tensor member function, P-21 ansysToFoam utility, U-89
&& APIfunctions model, U-99
tensor member function, P-21 applications, U-67
^ Apply button, U-164, U-168
tensor member function, P-21 applyBoundaryLayer utility, U-88
<LESModel>Coeffs keyword, U-184 applyWallFunctionBoundaryConditions
<RASModel>Coeffs keyword, U-184 utility,
<delta>Coeffs keyword, U-184 U-88
0.000000e+00 directory, U-104 arbitrarily unstructured, P-27
1-dimensional mesh, U-130 arc
1D mesh, U-130 keyword entry, U-139
2-dimensional mesh, U-130 arc keyword, U-138
2D mesh, U-130 As keyword, U-182
ascii
数 keyword entry, U-112
0 directory, U-104 attachMesh utility, U-90
Auto Accept button, U-168
A autoMesh
library, U-96 blocks keyword, U-20, U-30, U-140
autoPatch utility, U-90 boundaries, U-132
autoRefineMesh utility, U-91 boundary, U-132
axes boundary
right-handed, U-136 dictionary, U-129, U-136
right-handed rectangular Cartesian, P- boundary keyword, U-141
11, boundary condition
U-18 alphaContactAngle, U-57
axi-symmetric cases, U-135, U-144 buoyantPressure, U-137
axi-symmetric mesh, U-130 calculated, U-136
cyclic, U-135, U-142
B directionMixed, U-136
background empty, P-59, P-65, U-18, U-130, U-135
process, U-24, U-79 fixedGradient, U-136
backward fixedValue, U-136
keyword entry, U-119 fluxCorrectedVelocity, U-137
Backward differencing, P-35 inlet, P-65
barotropicCompressibilityModels inletOutlet, U-137
library, U-98 mixed, U-136
basicMultiComponentMixture model, U-98, movingWallVelocity, U-137
U-180 outlet, P-65
basicSolidThermo outletInlet, U-137
library, U-99 partialSlip, U-137
basicThermophysicalModels patch, U-134
library, U-98 pressureDirectedInletVelocity, U-137
binary pressureInletVelocity, U-137
keyword entry, U-112 pressureOutlet, P-59
BirdCarreau model, U-101 pressureTransmissive, U-137
blended differencing, P-34 processor, U-135
block setup, U-20
expansion ratio, U-140 slip, U-137
block keyword, U-138 supersonicFreeStream, U-137
blocking surfaceNormalFixedValue, U-137
keyword entry, U-78 symmetryPlane, P-59, U-134
blockMesh totalPressure, U-137
library, U-96 turbulentInlet, U-137
blockMesh solver, P-43 wall, U-40
blockMesh utility, U-37, U-89, U-136 wall, P-59, P-65, U-57, U-134
blockMesh executable wedge, U-130, U-135, U-144
vertex numbering, U-140 zeroGradient, U-136
blockMeshDict boundary conditions, P-39
dictionary, U-18, U-20, U-35, U-48, U- Dirichlet, P-39
136, inlet, P-40
U-144 Neumann, P-39
no-slip impermeable wall, P-40 cases, U-103
outlet, P-40 castellatedMesh keyword, U-146
physical, P-40 castellatedMeshControls
symmetry plane, P-40 dictionary, U-147–U-149
boundaryField keyword, U-21, U-108 castellatedMeshControls keyword, U-146
boundaryFoam solver, U-84 cavitatingDyMFoam solver, U-85
bounded cavitatingFoam solver, U-85
keyword entry, U-117, U-118 cavity flow, U-17
boxToCell keyword, U-58 ccm26ToFoam utility, U-89
boxTurb utility, U-88 CEI ARCH
breaking of a dam, U-55 environment variable, U-173
buoyantBoussinesqPimpleFoam solver, U- CEI HOME
86 environment variable, U-173
buoyantBoussinesqSimpleFoam solver, U- cell
86 expansion ratio, U-140
buoyantPimpleFoam solver, U-86 cell class, P-27
buoyantPressure cell
boundary condition, U-137 keyword entry, U-174
buoyantSimpleFoam solver, U-86 cellLimited
button keyword entry, U-117
Apply, U-164, U-168 cellPoint
Auto Accept, U-168 keyword entry, U-174
Choose Preset, U-166 cellPointFace
Delete, U-164 keyword entry, U-174
Edit Color Map, U-165 cells
Enable Line Series, U-34 dictionary, U-136
Orientation Axes, U-24, U-167 central differencing, P-34
Refresh Times, U-25 cfdTools tools, U-96
Rescale to Data Range, U-25 cfx4ToFoam utility, U-89, U-154
Reset, U-164 changeDictionary utility, U-88
Set Ambient Color, U-166 Charts window panel, U-168
Update GUI, U-165 checkMesh utility, U-90, U-155
Use Parallel Projection, U-24 chemFoam solver, U-86
Use parallel projection, U-167 chemistryModel
library, U-99
C chemistryModel model, U-99
C++ syntax chemistrySolver model, U-99
/*...*/, U-76 chemkinToFoam utility, U-95
//, U-76 Choose Preset button, U-166
# include, U-70, U-76 chtMultiRegionSimpleFoam solver, U-87
cacheAgglomeration keyword, U-123 chtMultiRegionFoam solver, U-87
calculated Chung
boundary condition, U-136 library, U-99
cAlpha keyword, U-61 class
cell, P-27 Color Scale window panel, U-166
dimensionSet, P-21, P-28, P-29 Colors window panel, U-168
face, P-27 compressibleInterDyMFoam solver, U-85
finiteVolumeCalculus, P-29 compressibleInterFoam solver, U-85
finiteVolumeMethod, P-29 compressibleMultiphaseInterFoam solver,
fvMesh, P-27 U-85
fvSchemes, P-32 combinePatchFaces utility, U-91
fvc, P-32 comments, U-76
fvm, P-32 commsType keyword, U-78
pointField, P-27 compressed
polyBoundaryMesh, P-27 keyword entry, U-112
polyMesh, P-27, U-127, U-129 compressibleLESModels
polyPatchList, P-27 library, U-101
polyPatch, P-27 compressibleRASModels
scalarField, P-25 library, U-100
scalar, P-19 constant directory, U-104, U-179
slice, P-27 constant model, U-98
symmTensorField, P-25 constTransport model, U-99, U-180
symmTensorThirdField, P-25 containers tools, U-96
tensorField, P-25 continuum
tensorThirdField, P-25 mechanics, P-11
tensor, P-19 control
vectorField, P-25 of time, U-111
vector, P-19, U-107 controlDict
word, P-21, P-27 dictionary, P-61, U-21, U-30, U-41, U-
class keyword, U-105 50,
clockTime U-60, U-104, U-160
keyword entry, U-111 controlDict file, P-46
cloud keyword, U-175 convection, see divergence, P-34
cloudFunctionObjects convergence, U-38
library, U-96 conversion
cmptAv library, U-97
tensor member function, P-21 convertToMeters keyword, U-138
Co utility, U-92 coordinate
coalChemistryFoam solver, U-87 system, P-11
coalCombustion coordinate system, U-18
library, U-97 corrected
cofactors keyword entry, U-117, U-118
tensor member function, P-21 Courant number, P-38, U-22
coldEngineFoam solver, U-86 Cp keyword, U-181
collapseEdges utility, U-91 cpuTime
Color By menu, U-166 keyword entry, U-111
Color Legend window, U-27 Crank Nicholson
Color Legend window panel, U-166 temporal discretisation, P-38
CrankNicholson decompose model, U-97
keyword entry, U-119 decomposePar utility, U-79, U-80, U-95
createExternalCoupledPatchGeometry decomposeParDict
utility, dictionary, U-79
U-88 decomposition
createBaffles utility, U-90 of field, U-79
createPatch utility, U-90 of mesh, U-79
createTurbulenceFields utility, U-92 decompositionMethods
cross product, see tensor, vector cross library, U-97
product decompression of a tank, P-58
CrossPowerLaw defaultFieldValues keyword, U-58
keyword entry, U-58 deformedGeom utility, U-90
CrossPowerLaw model, U-101 Delete button, U-164
cubeRootVolDelta model, U-101 delta keyword, U-81, U-184
cubicCorrected deltaT keyword, U-111
keyword entry, U-119 dependencies, U-70
cubicCorrection dependency lists, U-70
keyword entry, U-116 det
curl, P-33 tensor member function, P-21
curl determinant, see tensor, determinant
fvc member function, P-33 dev
Current Time Controls menu, U-25, U-165 tensor member function, P-21
curve keyword, U-175 diag
Cv keyword, U-181 tensor member function, P-21
cyclic diagonal
boundary condition, U-135, U-142 keyword entry, U-121, U-122
cyclic DIC
keyword entry, U-135 keyword entry, U-122
cylinder DICGaussSeidel
flow around a, P-41 keyword entry, U-122
dictionary
D LESProperties, U-184
d2dt2 PISO, U-23
fvc member function, P-33 blockMeshDict, U-18, U-20, U-35, U-
fvm member function, P-33 48,
dam U-136, U-144
breaking of a, U-55 boundary, U-129, U-136
datToFoam utility, U-89 castellatedMeshControls, U-147–U-149
db tools, U-96 cells, U-136
ddt controlDict, P-61, U-21, U-30, U-41,
fvc member function, P-33 U-50,
fvm member function, P-33 U-60, U-104, U-160
DeardorffDiffStress model, U-101 decomposeParDict, U-79
debug keyword, U-146 faces, U-129, U-136
fvSchemes, U-61, U-104, U-113 discretisation
fvSolution, U-104, U-120 equation, P-29
mechanicalProperties, U-49 Display window panel, U-23, U-25, U-164,
neighbour, U-129 U-165
owner, U-129 distance
points, U-129, U-136 keyword entry, U-149, U-175
thermalProperties, U-50 distributed model, U-97
thermophysicalProperties, U-179 distributed keyword, U-81, U-82
transportProperties, U-21, U-38, U-41 distributionModels
turbulenceProperties, U-40, U-59, U- library, U-97
183 div
differencing fvc member function, P-33
Backward, P-35 fvm member function, P-33
blended, P-34 divergence, P-33, P-35
central, P-34 divSchemes keyword, U-114
Euler implicit, P-35 dnsFoam solver, U-86
Gamma, P-34 doLayers keyword, U-146
MINMOD, P-34 double inner product, see tensor,double
SUPERBEE, P-34 inner
upwind, P-34 product
van Leer, P-34 DPMFoam solver, U-87
DILU dsmc
keyword entry, U-122 library, U-97
dimension dsmcFieldsCalc utility, U-93
checking in OpenFOAM, P-21, U-107 dsmcFoam solver, U-87
dimensional units, U-107 dsmcInitialise utility, U-88
dimensioned<Type> template class, P-21 dx
dimensionedTypes tools, U-96 keyword entry, U-174
dimensions keyword, U-20, U-108 dynamicFvMesh
dimensionSet class, P-21, P-28, P-29 library, U-97
dimensionSet tools, U-96 dynamicMesh
directionMixed library, U-96
boundary condition, U-136 dynLagrangian model, U-101
directory dynOneEqEddy model, U-101
0.000000e+00, U-104
0, U-104
E
Make, U-71 eConstThermo model, U-99, U-179
constant, U-104, U-179 edgeGrading keyword, U-140
fluentInterface, U-170 edgeMesh
polyMesh, U-104, U-129 library, U-97
processorN, U-80 edges keyword, U-138
run, U-103 Edit menu, U-167, U-168
system, P-46, U-104 Edit Color Map button, U-165
tutorials, P-41, U-17 egrMixture model, U-98, U-180
electrostaticFoam solver, U-88 equationOfState keyword, U-181
empty equilibriumCO utility, U-95
boundary condition, P-59, P-65, U-18, equilibriumFlameT utility, U-95
U-130, U-135 errorReduction keyword, U-153
empty Euler
keyword entry, U-135 keyword entry, U-119
Enable Line Series button, U-34 Euler implicit
endTime keyword, U-22, U-111 differencing, P-35
energy keyword, U-181 temporal discretisation, P-38
engine examples
library, U-97 decompression of a tank, P-58
engineCompRatio utility, U-93 flow around a cylinder, P-41
engineFoam solver, U-86 flow over backward step, P-49
engineSwirl utility, U-88 Hartmann problem, P-63
ensight74FoamExec utility, U-172 supersonic flow over forward step, P-
ENSIGHT7 INPUT 54
environment variable, U-173 execFlowFunctionObjects utility, U-93
ENSIGHT7 READER expandDictionary utility, U-95
environment variable, U-173 expansionRatio keyword, U-152
ensightFoamReader utility, U-91 explicit
enstrophy utility, U-92 temporal discretisation, P-38
environment variable extrude2DMesh utility, U-89
CEI ARCH, U-173 extrudeMesh utility, U-89
CEI HOME, U-173 extrudeToRegionMesh utility, U-89
ENSIGHT7 INPUT, U-173
ENSIGHT7 READER, U-173
F
FOAM RUN, U-103 face class, P-27
WM ARCH OPTION, U-74 face keyword, U-175
WM ARCH, U-74 faceAgglomerate utility, U-88
WM COMPILER BIN, U-74 faceAreaPair
WM COMPILER DIR, U-74 keyword entry, U-123
WM COMPILER LIB, U-74 faceLimited
WM COMPILER, U-74 keyword entry, U-117
WM COMPILE OPTION, U-74 faces
WM DIR, U-74 dictionary, U-129, U-136
WM MPLIB, U-74 FDIC
WM OPTIONS, U-74 keyword entry, U-122
WM PRECISION OPTION, U-74 featureAngle keyword, U-152
WM PROJECT DIR, U-74 features keyword, U-147, U-148
WM PROJECT INST DIR, U-74 field
WM PROJECT USER DIR, U-74 U, U-22
WM PROJECT VERSION, U-74 p, U-22
WM PROJECT, U-74 decomposition, U-79
wmake, U-73 FieldField<Type> template class, P-28
fieldFunctionObjects flow
library, U-96 free surface, U-55
fields, P-25 laminar, U-17
mapping, U-160 steady, turbulent, P-49
fields tools, U-96 supersonic, P-55
fields keyword, U-174 turbulent, U-17
Field<Type> template class, P-25 flow around a cylinder, P-41
fieldValues keyword, U-58 flow over backward step, P-49
file flowType utility, U-92
Make/files, U-72 fluent3DMeshToFoam utility, U-89
controlDict, P-46 fluentInterface directory, U-170
files, U-71 fluentMeshToFoam utility, U-89, U-154
g, U-59 fluxCorrectedVelocity
options, U-71 boundary condition, U-137
snappyHexMeshDict, U-145 fluxRequired keyword, U-114
transportProperties, U-58 OpenFOAM
file format, U-104 cases, U-103
fileFormats FOAM RUN
library, U-97 environment variable, U-103
fileModificationChecking keyword, U-78 foamCalc utility, U-33, U-93
fileModificationSkew keyword, U-78 foamCalcFunctions
files file, U-71 library, U-96
filteredLinear2 foamCorrectVrt script/alias, U-158
keyword entry, U-116 foamDataToFluent utility, U-91, U-170
finalLayerThickness keyword, U-152 foamDebugSwitches utility, U-95
financialFoam solver, U-88 FoamFile keyword, U-105
finite volume foamFile
discretisation, P-23 keyword entry, U-174
mesh, P-27 foamFormatConvert utility, U-95
finiteVolume foamHelp utility, U-95
library, U-96 foamInfoExec utility, U-95
finiteVolume tools, U-96 foamJob script/alias, U-177
finiteVolumeCalculus class, P-29 foamListTimes utility, U-93
finiteVolumeMethod class, P-29 foamLog script/alias, U-177
fireFoam solver, U-86 foamMeshToFluent utility, U-89, U-170
firstTime keyword, U-111 foamToEnsight utility, U-91
fixed foamToEnsightParts utility, U-91
keyword entry, U-112 foamToGMV utility, U-91
fixedGradient foamToStarMesh utility, U-89
boundary condition, U-136 foamToSurface utility, U-89
fixedValue foamToTecplot360 utility, U-91
boundary condition, U-136 foamToVTK utility, U-91
flattenMesh utility, U-90 foamUpgradeCyclics utility, U-88
floatTransfer keyword, U-78 foamUpgradeFvSolution utility, U-88
foamyHexMeshBackgroundMesh utility, U- library, U-97
89 fvSchemes
foamyHexMeshSurfaceSimplify utility, U- dictionary, U-61, U-104, U-113
89 fvSchemes class, P-32
foamyHexMesh utility, U-89 fvSchemes
foamyQuadMesh utility, U-89 menu entry, U-51
forces fvSolution
library, U-96 dictionary, U-104, U-120
foreground
process, U-24
G
format keyword, U-105 g file, U-59
fourth gambitToFoam utility, U-89, U-154
keyword entry, U-117, U-118 GAMG
functions keyword, U-113 keyword entry, U-52, U-121, U-122
fvc class, P-32 Gamma
fvc member function keyword entry, U-116
curl, P-33 Gamma differencing, P-34
d2dt2, P-33 Gauss
ddt, P-33 keyword entry, U-117
div, P-33 Gauss’s theorem, P-32
gGrad, P-33 GaussSeidel
grad, P-33 keyword entry, U-122
laplacian, P-33 General window panel, U-167, U-168
lsGrad, P-33 general
snGrad, P-33 keyword entry, U-112
snGradCorrection, P-33 genericFvPatchField
sqrGradGrad, P-33 library, U-97
fvDOM geometric-algebraic multi-grid, U-122
library, U-98 GeometricBoundaryField template class, P-
FVFunctionObjects 28
library, U-96 geometricField<Type> template class, P-28
fvm class, P-32 geometry keyword, U-146
fvm member function gGrad
d2dt2, P-33 fvc member function, P-33
ddt, P-33 global tools, U-96
div, P-33 gmshToFoam utility, U-89
laplacian, P-33 gnuplot
Su, P-33 keyword entry, U-112, U-174
SuSp, P-33 grad
fvMatrices tools, U-96 fvc member function, P-33
fvMatrix template class, P-29 (Grad Grad) squared, P-33
fvMesh class, P-27 gradient, P-33, P-36
fvMesh tools, U-96 Gauss scheme, P-36
fvMotionSolvers Gauss’s theorem, U-51
least square fit, U-51 identity, see tensor, identity
least squares method, P-36, U-51 incompressibleLESModels
surface normal, P-36 library, U-101
gradSchemes keyword, U-114 incompressiblePerfectGas model, U-99, U-
graph tools, U-96 179
graphFormat keyword, U-112 incompressibleRASModels
GuldersEGRLaminarFlameSpeed model, U- library, U-100
98 incompressibleTransportModels
GuldersLaminarFlameSpeed model, U-98 library, P-50, U-101
incompressibleTurbulenceModels
H library, P-50
hConstThermo model, U-99, U-179 index
heheupsiReactionThermo model, U-98, U- notation, P-12, P-13
180 Information window panel, U-164
Help menu, U-167 inhomogeneousMixture model, U-98, U-
hePsiThermo model, U-98, U-180 180
heRhoThermo model, U-98, U-180 inlet
HerschelBulkley model, U-101 boundary condition, P-65
hExponentialThermo inletOutlet
library, U-99 boundary condition, U-137
Hf keyword, U-181 inner product, see tensor, inner product
hierarchical inotify
keyword entry, U-80, U-81 keyword entry, U-78
highCpCoeffs keyword, U-182 inotifyMaster
homogenousDynOneEqEddy model, U-101 keyword entry, U-78
homogenousDynSmagorinsky model, U- inside
101 keyword entry, U-149
homogeneousMixture model, U-98, U-180 insideCells utility, U-90
hPolynomialThermo model, U-99, U-180 interPhaseChangeDyMFoam solver, U-85
interPhaseChangeFoam solver, U-85
I interDyMFoam solver, U-85
I interfaceProperties
tensor member function, P-21 library, U-102
icoFoam solver, U-17, U-21, U-22, U-24, interfaceProperties model, U-102
U-84 interFoam solver, U-85
icoPolynomial model, U-99, U-179 interMixingFoam solver, U-85
icoUncoupledKinematicParcelDyMFoam internalField keyword, U-21, U-108
solver, interpolation tools, U-96
U-87 interpolationScheme keyword, U-174
icoUncoupledKinematicParcelFoam solver, interpolations tools, U-96
U-87 interpolationSchemes keyword, U-114
ideasToFoam utility, U-154 inv
ideasUnvToFoam utility, U-89 tensor member function, P-21
identities, see tensor, identities iterations
maximum, U-121 deltaT, U-111
delta, U-81, U-184
J dimensions, U-20, U-108
janafThermo model, U-99, U-180 distributed, U-81, U-82
jobControl divSchemes, U-114
library, U-96 doLayers, U-146
jplot edgeGrading, U-140
keyword entry, U-112, U-174 edges, U-138
endTime, U-22, U-111
K energy, U-181
kEpsilon model, U-100 equationOfState, U-181
keyword errorReduction, U-153
As, U-182 expansionRatio, U-152
Cp, U-181 face, U-175
Cv, U-181 featureAngle, U-152
FoamFile, U-105 features, U-147, U-148
Hf, U-181 fieldValues, U-58
LESModel, U-184 fields, U-174
Pr, U-182 fileModificationChecking, U-78
RASModel, U-184 fileModificationSkew, U-78
Tcommon, U-182 finalLayerThickness, U-152
Thigh, U-182 firstTime, U-111
Tlow, U-182 floatTransfer, U-78
Ts, U-182 fluxRequired, U-114
addLayersControls, U-146 format, U-105
adjustTimeStep, U-60, U-112 functions, U-113
agglomerator, U-123 geometry, U-146
arc, U-138 gradSchemes, U-114
blocks, U-20, U-30, U-140 graphFormat, U-112
block, U-138 highCpCoeffs, U-182
boundaryField, U-21, U-108 internalField, U-21, U-108
boundary, U-141 interpolationSchemes, U-114
boxToCell, U-58 interpolationScheme, U-174
cAlpha, U-61 laplacianSchemes, U-114
cacheAgglomeration, U-123 latestTime, U-38
castellatedMeshControls, U-146 layers, U-152
castellatedMesh, U-146 leastSquares, U-51
class, U-105 levels, U-150
cloud, U-175 libs, U-78, U-112
commsType, U-78 locationInMesh, U-148, U-149
convertToMeters, U-138 location, U-105
curve, U-175 lowCpCoeffs, U-182
debug, U-146 manualCoeffs, U-81
defaultFieldValues, U-58 maxAlphaCo, U-60
maxBoundarySkewness, U-153 nSmoothNormals, U-152
maxConcave, U-153 nSmoothPatch, U-150
maxCo, U-60, U-112 nSmoothScale, U-153
maxDeltaT, U-60 nSmoothSurfaceNormals, U-152
maxFaceThicknessRatio, U-152 nSmoothThickness, U-152
maxGlobalCells, U-148 nSolveIter, U-150
maxInternalSkewness, U-153 neighbourPatch, U-142
maxIter, U-121 numberOfSubdomains, U-81
maxLocalCells, U-148 n, U-81
maxNonOrtho, U-153 object, U-105
maxThicknessToMedialRatio, U-152 order, U-81
mergeLevels, U-123 pRefCell, U-23, U-124
mergePatchPairs, U-138 pRefValue, U-23, U-124
mergeTolerance, U-146 p rhgRefCell, U-124
meshQualityControls, U-146 p rhgRefValue, U-124
method, U-81 patchMap, U-160
midPointAndFace, U-175 patches, U-138
midPoint, U-175 preconditioner, U-121, U-122
minArea, U-153 pressure, U-49
minDeterminant, U-153 printCoeffs, U-41, U-184
minFaceWeight, U-153 processorWeights, U-80
minFlatness, U-153 processorWeights, U-81
minMedianAxisAngle, U-152 purgeWrite, U-112
minRefinementCells, U-148 refGradient, U-136
minThickness, U-152 refinementRegions, U-148, U-150
minTriangleTwist, U-153 refinementSurfaces, U-148
minTwist, U-153 refinementRegions, U-149
minVolRatio, U-153 regions, U-58
minVol, U-153 relTol, U-52, U-121
mode, U-149 relativeSizes, U-152
molWeight, U-181 relaxed, U-153
mu, U-182 resolveFeatureAngle, U-148
nAlphaSubCycles, U-61 roots, U-81, U-82
nBufferCellsNoExtrude, U-152 runTimeModifiable, U-112
nCellsBetweenLevels, U-148 scotchCoeffs, U-81
nFaces, U-130 setFormat, U-174
nFinestSweeps, U-123 sets, U-174
nGrow, U-152 simpleGrading, U-140
nLayerIter, U-152 simulationType, U-40, U-59, U-183
nMoles, U-181 smoother, U-123
nPostSweeps, U-123 snGradSchemes, U-114
nPreSweeps, U-123 snapControls, U-146
nRelaxIter, U-150, U-152 snap, U-146
nRelaxedIter, U-152 solvers, U-120
solver, U-52, U-121 GaussSeidel, U-122
specie, U-181 Gauss, U-117
spline, U-138 LESModel, U-40, U-183
startFace, U-130 MGridGen, U-123
startFrom, U-22, U-111 MUSCL, U-116
startTime, U-22, U-111 Newtonian, U-58
stopAt, U-111 PBiCG, U-121
strategy, U-80, U-81 PCG, U-121
surfaceFormat, U-174 QUICK, U-119
surfaces, U-174 RASModel, U-40, U-183
thermoType, U-179 SFCD, U-116, U-119
thermodynamics, U-181 UMIST, U-115
timeFormat, U-112 adjustableRunTime, U-60, U-111
timePrecision, U-112 arc, U-139
timeScheme, U-114 ascii, U-112
tolerance, U-52, U-121, U-150 backward, U-119
topoSetSource, U-58 binary, U-112
traction, U-49 blocking, U-78
transport, U-181 bounded, U-117, U-118
turbulence, U-184 cellLimited, U-117
type, U-132, U-133 cellPointFace, U-174
uniform, U-175 cellPoint, U-174
valueFraction, U-136 cell, U-174
value, U-21, U-136 clockTime, U-111
version, U-105 compressed, U-112
vertices, U-20, U-138, U-139 corrected, U-117, U-118
writeCompression, U-112 cpuTime, U-111
writeControl, U-22, U-60, U-111 cubicCorrected, U-119
writeFormat, U-54, U-112 cubicCorrection, U-116
writeInterval, U-22, U-31, U-111 cyclic, U-135
writePrecision, U-112 diagonal, U-121, U-122
<LESModel>Coeffs, U-184 distance, U-149, U-175
<RASModel>Coeffs, U-184 dx, U-174
<delta>Coeffs, U-184 empty, U-135
keyword entry faceAreaPair, U-123
CrankNicholson, U-119 faceLimited, U-117
CrossPowerLaw, U-58 filteredLinear2, U-116
DICGaussSeidel, U-122 fixed, U-112
DIC, U-122 foamFile, U-174
DILU, U-122 fourth, U-117, U-118
Euler, U-119 general, U-112
FDIC, U-122 gnuplot, U-112, U-174
GAMG, U-52, U-121, U-122 hierarchical, U-80, U-81
Gamma, U-116 inotifyMaster, U-78
inotify, U-78 upwind, U-116, U-119
inside, U-149 vanLeer, U-116
jplot, U-112, U-174 vtk, U-174
laminar, U-40, U-183 wall, U-135
latestTime, U-111 wedge, U-135
leastSquares, U-117 writeControl, U-111
limitedCubic, U-116 writeNow, U-111
limitedLinear, U-116 xmgr, U-112, U-174
limited, U-117, U-118 xyz, U-175
linearUpwind, U-116, U-119 x, U-175
linear, U-116, U-119 y, U-175
line, U-139 z, U-175
localEuler, U-119 kivaToFoam utility, U-89
manual, U-80, U-81 kkLOmega model, U-100
metis, U-81 kOmega model, U-100
midPoint, U-116 kOmegaSST model, U-100
nextWrite, U-111 kOmegaSSTSAS model, U-101
noWriteNow, U-111 Kronecker delta, P-16
nonBlocking, U-78
none, U-115, U-122
L
null, U-174 lagrangian
outside, U-149 library, U-97
patch, U-135, U-176 lagrangianIntermediate
polyLine, U-139 library, U-97
polySpline, U-139 Lambda2 utility, U-92
processor, U-135 LamBremhorstKE model, U-100
raw, U-112, U-174 laminar model, U-100, U-101
runTime, U-31, U-111 laminar
scheduled, U-78 keyword entry, U-40, U-183
scientific, U-112 laminarFlameSpeedModels
scotch, U-80, U-81 library, U-98
simpleSpline, U-139 laplaceFilter model, U-100
simple, U-80, U-81 Laplacian, P-34
skewLinear, U-116, U-119 laplacian, P-33
smoothSolver, U-121 laplacian
startTime, U-22, U-111 fvc member function, P-33
steadyState, U-119 fvm member function, P-33
stl, U-174 laplacianFoam solver, U-83
symmetryPlane, U-135 laplacianSchemes keyword, U-114
timeStampMaster, U-78 latestTime
timeStamp, U-78 keyword entry, U-111
timeStep, U-22, U-31, U-111 latestTime keyword, U-38
uncompressed, U-112 LaunderGibsonRSTM model, U-100
uncorrected, U-117, U-118 LaunderSharmaKE model, U-100
layers keyword, U-152 dsmc, U-97
leastSquares dynamicFvMesh, U-97
keyword entry, U-117 dynamicMesh, U-96
leastSquares keyword, U-51 edgeMesh, U-97
LESdeltas engine, U-97
library, U-100 fieldFunctionObjects, U-96
LESfilters fileFormats, U-97
library, U-100 finiteVolume, U-96
LESModel foamCalcFunctions, U-96
keyword entry, U-40, U-183 forces, U-96
LESModel keyword, U-184 fvDOM, U-98
LESProperties fvMotionSolvers, U-97
dictionary, U-184 genericFvPatchField, U-97
levels keyword, U-150 hExponentialThermo, U-99
libraries, U-67 incompressibleLESModels, U-101
library incompressibleRASModels, U-100
Chung, U-99 incompressibleTransportModels, P-50,
FVFunctionObjects, U-96 U-101
LESdeltas, U-100 incompressibleTurbulenceModels, P-50
LESfilters, U-100 interfaceProperties, U-102
MGridGenGAMGAgglomeration, U- jobControl, U-96
97 lagrangianIntermediate, U-97
ODE, U-97 lagrangian, U-97
OSspecific, U-97 laminarFlameSpeedModels, U-98
OpenFOAM, U-96 linear, U-99
P1, U-98 liquidMixtureProperties, U-99
PV3FoamReader, U-163 liquidProperties, U-99
PVFoamReader, U-163 meshTools, U-97
SLGThermo, U-99 molecularMeasurements, U-97
Wallis, U-99 molecule, U-97
autoMesh, U-96 opaqueSolid, U-98
barotropicCompressibilityModels, U- pairPatchAgglomeration, U-97
98 postCalc, U-96
basicSolidThermo, U-99 potential, U-97
basicThermophysicalModels, U-98 primitive, P-19
blockMesh, U-96 radiationModels, U-98
chemistryModel, U-99 randomProcesses, U-97
cloudFunctionObjects, U-96 reactionThermophysicalModels, U-98
coalCombustion, U-97 sampling, U-96
compressibleLESModels, U-101 solidChemistryModel, U-99
compressibleRASModels, U-100 solidMixtureProperties, U-100
conversion, U-97 solidParticle, U-97
decompositionMethods, U-97 solidProperties, U-99
distributionModels, U-97 solidSpecie, U-100
solidThermo, U-100 List<Type> template class, P-25
specie, U-99 localEuler
spray, U-97 keyword entry, U-119
surfMesh, U-97 location keyword, U-105
surfaceFilmModels, U-102 locationInMesh keyword, U-148, U-149
systemCall, U-96 lowCpCoeffs keyword, U-182
thermophysicalFunctions, U-99 lowReOneEqEddy model, U-101
thermophysical, U-179 LRDDiffStress model, U-101
topoChangerFvMesh, U-97 LRR model, U-100
triSurface, U-97 lsGrad
turbulence, U-97 fvc member function, P-33
twoPhaseProperties, U-102 LTSInterFoam solver, U-85
utilityFunctionObjects, U-96 LTSReactingFoam solver, U-86
viewFactor, U-98 LTSReactingParcelFoam solver, U-87
vtkFoam, U-163
vtkPV3Foam, U-163
M
libs keyword, U-78, U-112 Mach utility, U-92
lid-driven cavity flow, U-17 mag
LienCubicKE model, U-100 tensor member function, P-21
LienCubicKELowRe model, U-100 magneticFoam solver, U-88
LienLeschzinerLowRe model, U-100 magnetohydrodynamics, P-63
Lights window panel, U-167 magSqr
limited tensor member function, P-21
keyword entry, U-117, U-118 Make directory, U-71
limitedCubic make script/alias, U-69
keyword entry, U-116 Make/files file, U-72
limitedLinear manual
keyword entry, U-116 keyword entry, U-80, U-81
line manualCoeffs keyword, U-81
keyword entry, U-139 mapFields utility, U-30, U-37, U-41, U-54,
Line Style menu, U-34 U-89,
linear U-160
library, U-99 mapping
linear fields, U-160
keyword entry, U-116, U-119 Marker Style menu, U-34
linearUpwind matrices tools, U-96
keyword entry, U-116, U-119 max
liquid tensor member function, P-21
electrically-conducting, P-63 maxAlphaCo keyword, U-60
liquidMixtureProperties maxBoundarySkewness keyword, U-153
library, U-99 maxCo keyword, U-60, U-112
liquidProperties maxConcave keyword, U-153
library, U-99 maxDeltaT keyword, U-60
lists, P-25 maxDeltaxyz model, U-101
maxFaceThicknessRatio keyword, U-152 2-dimensional, U-130
maxGlobalCells keyword, U-148 2D, U-130
maximum iterations, U-121 axi-symmetric, U-130
maxInternalSkewness keyword, U-153 basic, P-27
maxIter keyword, U-121 block structured, U-136
maxLocalCells keyword, U-148 decomposition, U-79
maxNonOrtho keyword, U-153 description, U-127
maxThicknessToMedialRatio keyword, U- finite volume, P-27
152 generation, U-136, U-145
mdEquilibrationFoam solver, U-87 grading, U-136, U-140
mdFoam solver, U-87 grading, example of, P-49
mdInitialise utility, U-89 non-orthogonal, P-41
mechanicalProperties refinement, P-58
dictionary, U-49 resolution, U-30
memory tools, U-96 specification, U-127
menu split-hex, U-145
Color By, U-166 Stereolithography (STL), U-145
Current Time Controls, U-25, U-165 surface, U-145
Edit, U-167, U-168 validity constraints, U-127
Help, U-167 Mesh Parts window panel, U-23
Line Style, U-34 meshes tools, U-96
Marker Style, U-34 meshQualityControls keyword, U-146
VCR Controls, U-25, U-165 meshTools
View, U-167 library, U-97
menu entry message passing interface
Plot Over Line, U-34 openMPI, U-81
Save Animation, U-169 method keyword, U-81
Save Screenshot, U-169 metis
Settings, U-168 keyword entry, U-81
Show Color Legend, U-25 metisDecomp model, U-97
Solid Color, U-166 MGridGenGAMGAgglomeration
Toolbars, U-167 library, U-97
View Settings..., U-24 MGridGen
View Settings, U-24, U-167 keyword entry, U-123
Wireframe, U-166 mhdFoam solver, P-65, U-88
fvSchemes, U-51 midPoint
mergeLevels keyword, U-123 keyword entry, U-116
mergeMeshes utility, U-90 midPoint keyword, U-175
mergeOrSplitBaffles utility, U-90 midPointAndFace keyword, U-175
mergePatchPairs keyword, U-138 min
mergeTolerance keyword, U-146 tensor member function, P-21
mesh minArea keyword, U-153
1-dimensional, U-130 minDeterminant keyword, U-153
1D, U-130 minFaceWeight keyword, U-153
minFlatness keyword, U-153 anisotropicFilter, U-100
minMedianAxisAngle keyword, U-152 basicMultiComponentMixture, U-98,
MINMOD differencing, P-34 U-180
minRefinementCells keyword, U-148 chemistryModel, U-99
minThickness keyword, U-152 chemistrySolver, U-99
minTriangleTwist keyword, U-153 constTransport, U-99, U-180
minTwist keyword, U-153 constant, U-98
minVol keyword, U-153 cubeRootVolDelta, U-101
minVolRatio keyword, U-153 decompose, U-97
mirrorMesh utility, U-90 distributed, U-97
mixed dynLagrangian, U-101
boundary condition, U-136 dynOneEqEddy, U-101
mixedSmagorinsky model, U-101 eConstThermo, U-99, U-179
mixtureAdiabaticFlameT utility, U-95 egrMixture, U-98, U-180
mode keyword, U-149 hConstThermo, U-99, U-179
model hPolynomialThermo, U-99, U-180
APIfunctions, U-99 hePsiThermo, U-98, U-180
BirdCarreau, U-101 heRhoThermo, U-98, U-180
CrossPowerLaw, U-101 heheupsiReactionThermo, U-98, U-180
DeardorffDiffStress, U-101 homogenousDynOneEqEddy, U-101
GuldersEGRLaminarFlameSpeed, U- homogenousDynSmagorinsky, U-101
98 homogeneousMixture, U-98, U-180
GuldersLaminarFlameSpeed, U-98 icoPolynomial, U-99, U-179
HerschelBulkley, U-101 incompressiblePerfectGas, U-99, U-
LRDDiffStress, U-101 179
LRR, U-100 inhomogeneousMixture, U-98, U-180
LamBremhorstKE, U-100 interfaceProperties, U-102
LaunderGibsonRSTM, U-100 janafThermo, U-99, U-180
LaunderSharmaKE, U-100 kEpsilon, U-100
LienCubicKELowRe, U-100 kOmegaSSTSAS, U-101
LienCubicKE, U-100 kOmegaSST, U-100
LienLeschzinerLowRe, U-100 kOmega, U-100
NSRDSfunctions, U-99 kkLOmega, U-100
Newtonian, U-101 laminar, U-100, U-101
NonlinearKEShih, U-100 laplaceFilter, U-100
PrandtlDelta, U-101 lowReOneEqEddy, U-101
RNGkEpsilon, U-100 maxDeltaxyz, U-101
RaviPetersen, U-98 metisDecomp, U-97
Smagorinsky2, U-101 mixedSmagorinsky, U-101
Smagorinsky, U-101 multiComponentMixture, U-98, U-180
SpalartAllmarasDDES, U-101 oneEqEddy, U-101
SpalartAllmarasIDDES, U-101 perfectFluid, U-99, U-179
SpalartAllmaras, U-100, U-101 polynomialTransport, U-99, U-180
adiabaticPerfectFluid, U-99, U-179 powerLaw, U-101
psiReactionThermo, U-98, U-180 multiphaseEulerFoam solver, U-85
psiuReactionThermo, U-98, U-180 multiphaseInterFoam solver, U-86
ptsotchDecomp, U-97 MUSCL
pureMixture, U-98, U-180 keyword entry, U-116
qZeta, U-100
reactingMixture, U-98, U-180
N
realizableKE, U-100 n keyword, U-81
reconstruct, U-97 nabla
rhoConst, U-99, U-179 operator, P-23
rhoReactionThermo, U-98, U-180 nAlphaSubCycles keyword, U-61
scaleSimilarity, U-101 nBufferCellsNoExtrude keyword, U-152
scotchDecomp, U-97 nCellsBetweenLevels keyword, U-148
simpleFilter, U-100 neighbour
singleStepReactingMixture, U-98, U- dictionary, U-129
180 neighbourPatch keyword, U-142
smoothDelta, U-101 netgenNeutralToFoam utility, U-89
specieThermo, U-99, U-180 Newtonian
spectEddyVisc, U-101 keyword entry, U-58
sutherlandTransport, U-99, U-180 Newtonian model, U-101
v2f, U-100 nextWrite
vanDriestDelta, U-101 keyword entry, U-111
veryInhomogeneousMixture, U-98, U- nFaces keyword, U-130
180 nFinestSweeps keyword, U-123
modifyMesh utility, U-91 nGrow keyword, U-152
molecularMeasurements nLayerIter keyword, U-152
library, U-97 nMoles keyword, U-181
molecule non-orthogonal mesh, P-41
library, U-97 nonBlocking
molWeight keyword, U-181 keyword entry, U-78
moveDynamicMesh utility, U-90 none
moveEngineMesh utility, U-90 keyword entry, U-115, U-122
moveMesh utility, U-90 NonlinearKEShih model, U-100
movingWallVelocity nonNewtonianIcoFoam solver, U-84
boundary condition, U-137 noWriteNow
MPI keyword entry, U-111
openMPI, U-81 nPostSweeps keyword, U-123
MRFInterFoam solver, U-85 nPreSweeps keyword, U-123
MRFMultiphaseInterFoam solver, U-85 nRelaxedIter keyword, U-152
mshToFoam utility, U-89 nRelaxIter keyword, U-150, U-152
mu keyword, U-182 nSmoothNormals keyword, U-152
multiComponentMixture model, U-98, U- nSmoothPatch keyword, U-150
180 nSmoothScale keyword, U-153
multigrid nSmoothSurfaceNormals keyword, U-152
geometric-algebraic, U-122 nSmoothThickness keyword, U-152
nSolveIter keyword, U-150 dictionary, U-129
NSRDSfunctions model, U-99
null
P
keyword entry, U-174 p field, U-22
numberOfSubdomains keyword, U-81 P1
library, U-98
O p rhgRefCell keyword, U-124
object keyword, U-105 p rhgRefValue keyword, U-124
objToVTK utility, U-90 pairPatchAgglomeration
ODE library, U-97
library, U-97 paraFoam, U-23, U-163
oneEqEddy model, U-101 parallel
Opacity text box, U-167 running, U-79
opaqueSolid partialSlip
library, U-98 boundary condition, U-137
OpenFOAM particleTracks utility, U-93
applications, U-67 patch
file format, U-104 boundary condition, U-134
libraries, U-67 patch
OpenFOAM keyword entry, U-135, U-176
library, U-96 patchAverage utility, U-92
OpenFOAM file syntax patches keyword, U-138
//, U-104 patchIntegrate utility, U-92
openMPI patchMap keyword, U-160
message passing interface, U-81 patchSummary utility, U-95
MPI, U-81 PBiCG
operator keyword entry, U-121
scalar, P-24 PCG
vector, P-23 keyword entry, U-121
Options window, U-168 pdfPlot utility, U-93
options file, U-71 PDRFoam solver, U-86
order keyword, U-81 PDRMesh utility, U-91
Orientation Axes button, U-24, U-167 Pe utility, U-92
orientFaceZone utility, U-90 perfectFluid model, U-99, U-179
OSspecific permutation symbol, P-15
library, U-97 pimpleDyMFoam solver, U-84
outer product, see tensor, outer product pimpleFoam solver, U-84
outlet Pipeline Browser window, U-23, U-164
boundary condition, P-65 PISO
outletInlet dictionary, U-23
boundary condition, U-137 pisoFoam solver, U-17, U-84
outside Plot Over Line
keyword entry, U-149 menu entry, U-34
owner plot3dToFoam utility, U-90
pointField class, P-27 boundary condition, P-59
pointField<Type> template class, P-29 pressureTransmissive
points boundary condition, U-137
dictionary, U-129, U-136 primitive
polyBoundaryMesh class, P-27 library, P-19
polyDualMesh utility, U-90 primitives tools, U-96
polyLine printCoeffs keyword, U-41, U-184
keyword entry, U-139 processorWeights keyword, U-80
polyMesh directory, U-104, U-129 probeLocations utility, U-93
polyMesh class, P-27, U-127, U-129 process
polynomialTransport model, U-99, U-180 background, U-24, U-79
polyPatch class, P-27 foreground, U-24
polyPatchList class, P-27 processor
polySpline boundary condition, U-135
keyword entry, U-139 processor
porousInterFoam solver, U-86 keyword entry, U-135
porousSimpleFoam solver, U-84 processorN directory, U-80
post-processing, U-163 processorWeights keyword, U-81
post-processing Properties window panel, U-25, U-164
paraFoam, U-163 psiReactionThermo model, U-98, U-180
postCalc psiuReactionThermo model, U-98, U-180
library, U-96 ptot utility, U-93
postChannel utility, U-93 ptsotchDecomp model, U-97
potentialFreeSurfaceFoam solver, U-86 pureMixture model, U-98, U-180
potential purgeWrite keyword, U-112
library, U-97 PV3FoamReader
potentialFoam solver, P-42, U-83 library, U-163
pow PVFoamReader
powerLaw model, U-101
pPrime2 utility, U-92
Q
Pr keyword, U-182 Q utility, U-92
PrandtlDelta model, U-101 QUICK
preconditioner keyword, U-121, U-122 keyword entry, U-119
pRefCell keyword, U-23, U-124 qZeta model, U-100
pRefValue keyword, U-23, U-124
pressure keyword, U-49
R
pressure waves R utility, U-92
in liquids, P-58 radiationModels
pressureDirectedInletVelocity library, U-98
boundary condition, U-137 randomProcesses
pressureInletVelocity library, U-97
boundary condition, U-137 RASModel
pressureOutlet keyword entry, U-40, U-183
RASModel keyword, U-184 rhoPimpleFoam solver, U-84
RaviPetersen model, U-98 rhoPimplecFoam solver, U-84
raw rhoReactingFoam solver, U-86
keyword entry, U-112, U-174 rhoReactionThermo model, U-98, U-180
reactingFoam solver, U-86 rhoSimpleFoam solver, U-84
reactingMixture model, U-98, U-180 rhoSimplecFoam solver, U-84
reactingParcelFilmFoam solver, U-87 rmdepall script/alias, U-74
reactingParcelFoam solver, U-87 RNGkEpsilon model, U-100
reactionThermophysicalModels roots keyword, U-81, U-82
library, U-98 rotateMesh utility, U-90
realizableKE model, U-100 run
reconstruct model, U-97 parallel, U-79
reconstructPar utility, U-83 run directory, U-103
reconstructParMesh utility, U-95 runTime
redistributePar utility, U-95 keyword entry, U-31, U-111
refGradient keyword, U-136 runTimeModifiable keyword, U-112
refineHexMesh utility, U-91
refinementRegions keyword, U-149
S
refinementLevel utility, U-91 sammToFoam utility, U-90
refinementRegions keyword, U-148, U-150 sample utility, U-93, U-173
refinementSurfaces keyword, U-148 sampling
refineMesh utility, U-90 library, U-96
refineWallLayer utility, U-91 Save Animation
Refresh Times button, U-25 menu entry, U-169
regions keyword, U-58 Save Screenshot
relative tolerance, U-121 menu entry, U-169
relativeSizes keyword, U-152 scalar, P-12
relaxed keyword, U-153 operator, P-24
relTol keyword, U-52, U-121 scalar class, P-19
removeFaces utility, U-91 scalarField class, P-25
Render View window, U-168 scalarTransportFoam solver, U-83
Render View window panel, U-168 scale
renumberMesh utility, U-90 tensor member function, P-21
Rescale to Data Range button, U-25 scalePoints utility, U-157
Reset button, U-164 scaleSimilarity model, U-101
resolveFeatureAngle keyword, U-148 scheduled
restart, U-38 keyword entry, U-78
Reynolds number, U-17, U-21 scientific
rhoPorousSimpleFoam solver, U-84 keyword entry, U-112
rhoReactingBuoyantFoam solver, U-86 scotch
rhoCentralDyMFoam solver, U-84 keyword entry, U-80, U-81
rhoCentralFoam solver, U-84 scotchCoeffs keyword, U-81
rhoConst model, U-99, U-179 scotchDecomp model, U-97
rhoLTSPimpleFoam solver, U-84 script/alias
foamCorrectVrt, U-158 library, U-99
foamJob, U-177 slice class, P-27
foamLog, U-177 slip
make, U-69 boundary condition, U-137
rmdepall, U-74 Smagorinsky model, U-101
wclean, U-73 Smagorinsky2 model, U-101
wmake, U-69 smapToFoam utility, U-91
second time derivative, P-33 smoothDelta model, U-101
Seed window, U-169 smoother keyword, U-123
selectCells utility, U-91 smoothSolver
Set Ambient Color button, U-166 keyword entry, U-121
setFields utility, U-58, U-89 snap keyword, U-146
setFormat keyword, U-174 snapControls keyword, U-146
sets keyword, U-174 snappyHexMesh utility
setSet utility, U-90 background mesh, U-146
setsToZones utility, U-90 cell removal, U-149
Settings cell splitting, U-147
menu entry, U-168 mesh layers, U-150
settlingFoam solver, U-86 meshing process, U-145
SFCD snapping to surfaces, U-150
keyword entry, U-116, U-119 snappyHexMesh utility, U-89, U-145
shallowWaterFoam solver, U-84 snappyHexMeshDict file, U-145
shape, U-140 snGrad
Show Color Legend fvc member function, P-33
menu entry, U-25 snGradCorrection
SI units, U-107 fvc member function, P-33
simpleReactingParcelFoam solver, U-87 snGradSchemes keyword, U-114
simple Solid Color
keyword entry, U-80, U-81 menu entry, U-166
simpleFilter model, U-100 solidChemistryModel
simpleFoam solver, P-50, U-84 library, U-99
simpleGrading keyword, U-140 solidDisplacementFoam solver, U-88
simpleSpline solidDisplacementFoam solver, U-50
keyword entry, U-139 solidEquilibriumDisplacementFoam solver,
simulationType keyword, U-40, U-59, U- U-88
183 solidMixtureProperties
singleCellMesh utility, U-91 library, U-100
singleStepReactingMixture model, U-98, U- solidParticle
180 library, U-97
skew solidProperties
skewLinear solidSpecie
keyword entry, U-116, U-119 library, U-100
SLGThermo solidThermo
library, U-100 interFoam, U-85
solver interMixingFoam, U-85
DPMFoam, U-87 interPhaseChangeDyMFoam, U-85
LTSInterFoam, U-85 interPhaseChangeFoam, U-85
LTSReactingFoam, U-86 laplacianFoam, U-83
LTSReactingParcelFoam, U-87 magneticFoam, U-88
MRFInterFoam, U-85 mdEquilibrationFoam, U-87
MRFMultiphaseInterFoam, U-85 mdFoam, U-87
PDRFoam, U-86 mhdFoam, P-65, U-88
SRFPimpleFoam, U-84 multiphaseEulerFoam, U-85
SRFSimpleFoam, U-84 multiphaseInterFoam, U-86
XiFoam, U-86 nonNewtonianIcoFoam, U-84
adjointShapeOptimizationFoam, U-83 pimpleDyMFoam, U-84
blockMesh, P-43 pimpleFoam, U-84
boundaryFoam, U-84 pisoFoam, U-17, U-84
buoyantBoussinesqPimpleFoam, U-86 porousInterFoam, U-86
buoyantBoussinesqSimpleFoam, U-86 porousSimpleFoam, U-84
buoyantPimpleFoam, U-86 potentialFreeSurfaceFoam, U-86
buoyantSimpleFoam, U-86 potentialFoam, P-42, U-83
cavitatingDyMFoam, U-85 reactingFoam, U-86
cavitatingFoam, U-85 reactingParcelFilmFoam, U-87
chemFoam, U-86 reactingParcelFoam, U-87
chtMultiRegionFoam, U-87 rhoCentralDyMFoam, U-84
chtMultiRegionSimpleFoam, U-87 rhoCentralFoam, U-84
coalChemistryFoam, U-87 rhoLTSPimpleFoam, U-84
coldEngineFoam, U-86 rhoPimpleFoam, U-84
compressibleInterDyMFoam, U-85 rhoPimplecFoam, U-84
compressibleInterFoam, U-85 rhoReactingFoam, U-86
compressibleMultiphaseInterFoam, U- rhoSimpleFoam, U-84
85 rhoSimplecFoam, U-84
dnsFoam, U-86 rhoPorousSimpleFoam, U-84
dsmcFoam, U-87 rhoReactingBuoyantFoam, U-86
electrostaticFoam, U-88 scalarTransportFoam, U-83
engineFoam, U-86 settlingFoam, U-86
financialFoam, U-88 shallowWaterFoam, U-84
fireFoam, U-86 simpleReactingParcelFoam, U-87
icoFoam, U-17, U-21, U-22, U-24, U- simpleFoam, P-50, U-84
84 solidDisplacementFoam, U-88
icoUncoupledKinematicParcelDyMFoa solidDisplacementFoam, U-50
m, solidEquilibriumDisplacementFoam,
U-87 U-88
icoUncoupledKinematicParcelFoam, sonicDyMFoam, U-84
U-87 sonicFoam, P-56, U-84
interDyMFoam, U-85 sonicLiquidFoam, P-59, U-85
sprayEngineFoam, U-87 steady flow
sprayFoam, U-87 turbulent, P-49
thermoFoam, U-87 steadyParticleTracks utility, U-93
twoLiquidMixingFoam, U-86 steadyState
twoPhaseEulerFoam, U-86 keyword entry, U-119
uncoupledKinematicParcelFoam, U-87 Stereolithography (STL), U-145
solver keyword, U-52, U-121 stitchMesh utility, U-91
solver relative tolerance, U-121 stl
solver tolerance, U-121 keyword entry, U-174
solvers keyword, U-120 stopAt keyword, U-111
sonicDyMFoam solver, U-84 strategy keyword, U-80, U-81
sonicFoam solver, P-56, U-84 streamFunction utility, U-92
sonicLiquidFoam solver, P-59, U-85 stress analysis of plate with hole, U-45
source, P-33 stressComponents utility, U-92
SpalartAllmaras model, U-100, U-101 Style window panel, U-23, U-166
SpalartAllmarasDDES model, U-101 Su
SpalartAllmarasIDDES model, U-101 fvm member function, P-33
specie subsetMesh utility, U-91
library, U-99 summation convention, P-13
specie keyword, U-181 SUPERBEE differencing, P-34
specieThermo model, U-99, U-180 supersonic flow, P-55
spectEddyVisc model, U-101 supersonic flow over forward step, P-54
spline keyword, U-138 supersonicFreeStream
splitCells utility, U-91 boundary condition, U-137
splitMesh utility, U-91 surfaceLambdaMuSmooth utility, U-94
splitMeshRegions utility, U-91 surface mesh, U-145
spray surfaceAdd utility, U-93
library, U-97 surfaceAutoPatch utility, U-93
sprayEngineFoam solver, U-87 surfaceBooleanFeatures utility, U-93
sprayFoam solver, U-87 surfaceCheck utility, U-93
sqr surfaceClean utility, U-93
tensor member function, P-21 surfaceCoarsen utility, U-93
sqrGradGrad surfaceConvert utility, U-93
fvc member function, P-33 surfaceFeatureConvert utility, U-93
SRFPimpleFoam solver, U-84 surfaceFeatureExtract utility, U-93, U-148
SRFSimpleFoam solver, U-84 surfaceField<Type> template class, P-29
star3ToFoam utility, U-90 surfaceFilmModels
star4ToFoam utility, U-90 library, U-102
startFace keyword, U-130 surfaceFind utility, U-93
startFrom keyword, U-22, U-111 surfaceFormat keyword, U-174
starToFoam utility, U-154 surfaceHookUp utility, U-93
startTime surfaceInertia utility, U-94
keyword entry, U-22, U-111 surfaceMesh tools, U-96
startTime keyword, U-22, U-111 surfaceMeshConvert utility, U-94
surfaceMeshConvertTesting utility, U-94 FieldField<Type>, P-28
surfaceMeshExport utility, U-94 Field<Type>, P-25
surfaceMeshImport utility, U-94 geometricField<Type>, P-28
surfaceMeshInfo utility, U-94 List<Type>, P-25
surfaceMeshTriangulate utility, U-94 pointField<Type>, P-29
surfaceNormalFixedValue surfaceField<Type>, P-29
boundary condition, U-137 volField<Type>, P-29
surfaceOrient utility, U-94 temporal discretisation, P-38
surfacePointMerge utility, U-94 Crank Nicholson, P-38
surfaceRedistributePar utility, U-94 Euler implicit, P-38
surfaceRefineRedGreen utility, U-94 explicit, P-38
surfaces keyword, U-174 in OpenFOAM, P-39
surfaceSplitByPatch utility, U-94 temporalInterpolate utility, U-93
surfaceSplitByTopology utility, U-94 tensor, P-11
surfaceSplitNonManifolds utility, U-94 addition, P-13
surfaceSubset utility, U-94 algebraic operations, P-13
surfaceToPatch utility, U-94 algebraic operations in OpenFOAM, P-
surfaceTransformPoints utility, U-94 19
surfMesh antisymmetric, see tensor, skew
library, U- calculus, P-23
97 classes in OpenFOAM, P-19
SuSp cofactors, P-18
fvm member function, P-33 component average, P-16
sutherlandTransport model, U-99, U-180 component maximum, P-16
symm component minimum, P-16
tensor member function, P-21 determinant, P-18
symmetryPlane deviatoric, P-17
boundary condition, P-59, U-134 diagonal, P-17
symmetryPlane dimension, P-12
keyword entry, U-135 double inner product, P-15
symmTensorField class, P-25 geometric transformation, P-16
symmTensorThirdField class, P-25 Hodge dual, P-18
system directory, P-46, U-104 hydrostatic, P-17
systemCall identities, P-17
library, U-96 identity, P-16
inner product, P-14
T inverse, P-18
T() magnitude, P-16
tensor member function, P-21 magnitude squared, P-16
Tcommon keyword, U-182 mathematics, P-11
template class notation, P-13
GeometricBoundaryField, P-28 nth power, P-16
fvMatrix, P-29 outer product, P-15
dimensioned<Type>, P-21 rank, P-12
rank 3, P-12 tensorField class, P-25
scalar division, P-14 tensorThirdField class, P-25
scalar multiplication, P-13 tetgenToFoam utility, U-90
scale function, P-16 text box
second rank, P-12 Opacity, U-167
skew, P-17 thermalProperties
square of, P-16 dictionary, U-50
subtraction, P-13 thermodynamics keyword, U-181
symmetric, P-17 thermoFoam solver, U-87
symmetric rank 2, P-12 thermophysical
symmetric rank 3, P-12 library, U-179
trace, P-17 thermophysicalFunctions
transformation, P-16 library, U-99
transpose, P-12, P-17 thermophysicalProperties
triple inner product, P-15 dictionary, U-179
vector cross product, P-15 thermoType keyword, U-179
tensor class, P-19 Thigh keyword, U-182
tensor member function time
*, P-21 control, U-111
+, P-21 time derivative, P-33
-, P-21 first, P-35
/, P-21 second, P-33, P-35
&, P-21 time step, U-22
&&, P-21 timeFormat keyword, U-112
^, P-21 timePrecision keyword, U-112
cmptAv, P-21 timeScheme keyword, U-114
cofactors, P-21 timeStamp
det, P-21 keyword entry, U-78
dev, P-21 timeStampMaster
diag, P-21 keyword entry, U-78
I, P-21 timeStep
inv, P-21 keyword entry, U-22, U-31, U-111
mag, P-21 Tlow keyword, U-182
magSqr, P-21 tolerance
max, P-21 solver, U-121
min, P-21 solver relative, U-121
pow, P-21 tolerance keyword, U-52, U-121, U-150
scale, P-21 Toolbars
skew, P-21 menu entry, U-167
sqr, P-21 tools
symm, P-21 algorithms, U-96
T(), P-21 cfdTools, U-96
tr, P-21 containers, U-96
transform, P-21 db, U-96
dimensionSet, U-96 turbulence model
dimensionedTypes, U-96 RAS, U-39
fields, U-96 turbulenceProperties
finiteVolume, U-96 dictionary, U-40, U-59, U-183
fvMatrices, U-96 turbulent flow
fvMesh, U-96 steady, P-49
global, U-96 turbulentInlet
graph, U-96 boundary condition, U-137
interpolations, U-96 tutorials
interpolation, U-96 breaking of a dam, U-55
matrices, U-96 lid-driven cavity flow, U-17
memory, U-96 stress analysis of plate with hole, U-45
meshes, U-96 tutorials directory, P-41, U-17
primitives, U-96 twoLiquidMixingFoam solver, U-86
surfaceMesh, U-96 twoPhaseEulerFoam solver, U-86
volMesh, U-96 twoPhaseProperties
topoChangerFvMesh library, U-102
library, U-97 type keyword, U-132, U-133
topoSet utility, U-91
topoSetSource keyword, U-58
U
totalPressure U field, U-22
boundary condition, U-137 Ucomponents utility, P-66
tr UMIST
tensor member function, P-21 keyword entry, U-115
trace, see tensor, trace uncompressed
traction keyword, U-49 keyword entry, U-112
transform uncorrected
tensor member function, P-21 keyword entry, U-117, U-118
transformPoints utility, U-91 uncoupledKinematicParcelFoam solver, U-
transport keyword, U-181 87
transportProperties uniform keyword, U-175
dictionary, U-21, U-38, U-41 units
transportProperties file, U-58 base, U-107
triple inner product, P-15 of measurement, P-21, U-107
triSurface S.I. base, P-21
library, U-97 SI, U-107
Ts keyword, U-182 Syst`eme International, U-107
turbulence United States Customary System, U-
dissipation, U-39 107
kinetic energy, U-39 USCS, U-107
length scale, U-40 Update GUI button, U-165
turbulence uprime utility, U-92
library, U-97 upwind
turbulence keyword, U-184 keyword entry, U-116, U-119
upwind differencing, P-34, U-61 ensightFoamReader, U-91
USCS units, U-107 enstrophy, U-92
Use Parallel Projection button, U-24 equilibriumCO, U-95
Use parallel projection button, U-167 equilibriumFlameT, U-95
utility execFlowFunctionObjects, U-93
Co, U-92 expandDictionary, U-95
Lambda2, U-92 extrude2DMesh, U-89
Mach, U-92 extrudeMesh, U-89
PDRMesh, U-91 extrudeToRegionMesh, U-89
Pe, U-92 faceAgglomerate, U-88
Q, U-92 flattenMesh, U-90
R, U-92 flowType, U-92
Ucomponents, P-66 fluent3DMeshToFoam, U-89
adiabaticFlameT, U-95 fluentMeshToFoam, U-89, U-154
ansysToFoam, U-89 foamCalc, U-33, U-93
applyBoundaryLayer, U-88 foamDataToFluent, U-91, U-170
applyWallFunctionBoundaryCondition foamDebugSwitches, U-95
s, foamFormatConvert, U-95
U-88 foamHelp, U-95
attachMesh, U-90 foamInfoExec, U-95
autoPatch, U-90 foamListTimes, U-93
autoRefineMesh, U-91 foamMeshToFluent, U-89, U-170
blockMesh, U-37, U-89, U-136 foamToEnsightParts, U-91
boxTurb, U-88 foamToEnsight, U-91
ccm26ToFoam, U-89 foamToGMV, U-91
cfx4ToFoam, U-89, U-154 foamToStarMesh, U-89
changeDictionary, U-88 foamToSurface, U-89
checkMesh, U-90, U-155 foamToTecplot360, U-91
chemkinToFoam, U-95 foamToVTK, U-91
collapseEdges, U-91 foamUpgradeCyclics, U-88
combinePatchFaces, U-91 foamUpgradeFvSolution, U-88
createBaffles, U-90 foamyHexMesh, U-89
createPatch, U-90 foamyQuadMesh, U-89
createTurbulenceFields, U-92 foamyHexMeshBackgroundMesh, U-
createExternalCoupledPatchGeometry, 89
U-88 foamyHexMeshSurfaceSimplify, U-89
datToFoam, U-89 gambitToFoam, U-89, U-154
decomposePar, U-79, U-80, U-95 gmshToFoam, U-89
deformedGeom, U-90 ideasToFoam, U-154
dsmcFieldsCalc, U-93 ideasUnvToFoam, U-89
dsmcInitialise, U-88 insideCells, U-90
engineCompRatio, U-93 kivaToFoam, U-89
engineSwirl, U-88 mapFields, U-30, U-37, U-41, U-54, U-
ensight74FoamExec, U-172 89,
U-160 snappyHexMesh, U-89, U-145
mdInitialise, U-89 splitCells, U-91
mergeMeshes, U-90 splitMeshRegions, U-91
mergeOrSplitBaffles, U-90 splitMesh, U-91
mirrorMesh, U-90 star3ToFoam, U-90
mixtureAdiabaticFlameT, U-95 star4ToFoam, U-90
modifyMesh, U-91 starToFoam, U-154
moveDynamicMesh, U-90 steadyParticleTracks, U-93
moveEngineMesh, U-90 stitchMesh, U-91
moveMesh, U-90 streamFunction, U-92
mshToFoam, U-89 stressComponents, U-92
netgenNeutralToFoam, U-89 subsetMesh, U-91
objToVTK, U-90 surfaceLambdaMuSmooth, U-94
orientFaceZone, U-90 surfaceAdd, U-93
pPrime2, U-92 surfaceAutoPatch, U-93
particleTracks, U-93 surfaceBooleanFeatures, U-93
patchAverage, U-92 surfaceCheck, U-93
patchIntegrate, U-92 surfaceClean, U-93
patchSummary, U-95 surfaceCoarsen, U-93
pdfPlot, U-93 surfaceConvert, U-93
plot3dToFoam, U-90 surfaceFeatureConvert, U-93
polyDualMesh, U-90 surfaceFeatureExtract, U-93, U-148
postChannel, U-93 surfaceFind, U-93
probeLocations, U-93 surfaceHookUp, U-93
ptot, U-93 surfaceInertia, U-94
reconstructParMesh, U-95 surfaceMeshConvertTesting, U-94
reconstructPar, U-83 surfaceMeshConvert, U-94
redistributePar, U-95 surfaceMeshExport, U-94
refineHexMesh, U-91 surfaceMeshImport, U-94
refineMesh, U-90 surfaceMeshInfo, U-94
refineWallLayer, U-91 surfaceMeshTriangulate, U-94
refinementLevel, U-91 surfaceOrient, U-94
removeFaces, U-91 surfacePointMerge, U-94
renumberMesh, U-90 surfaceRedistributePar, U-94
rotateMesh, U-90 surfaceRefineRedGreen, U-94
sammToFoam, U-90 surfaceSplitByPatch, U-94
sample, U-93, U-173 surfaceSplitByTopology, U-94
scalePoints, U-157 surfaceSplitNonManifolds, U-94
selectCells, U-91 surfaceSubset, U-94
setFields, U-58, U-89 surfaceToPatch, U-94
setSet, U-90 surfaceTransformPoints, U-94
setsToZones, U-90 temporalInterpolate, U-93
singleCellMesh, U-91 tetgenToFoam, U-90
smapToFoam, U-91 topoSet, U-91
transformPoints, U-91 menu entry, U-24
uprime, U-92 viewFactor
viewFactorsGen, U-89 library, U-98
vorticity, U-92 viewFactorsGen utility, U-89
vtkUnstructuredToFoam, U-90 viscosity
wallFunctionTable, U-89 kinematic, U-21, U-41
wallGradU, U-92 volField<Type> template class, P-29
wallHeatFlux, U-92 volMesh tools, U-96
wallShearStress, U-92 vorticity utility, U-92
wdot, U-93 vtk
writeCellCentres, U-93 keyword entry, U-174
writeMeshObj, U-90 vtkFoam
yPlusLES, U-92 library, U-163
yPlusRAS, U-92 vtkPV3Foam
zipUpMesh, U-91 library, U-163
utilityFunctionObjects vtkUnstructuredToFoam utility, U-90
library, U-96
W
V wall
v2f model, U-100 boundary condition, P-59, P-65, U-57,
value keyword, U-21, U-136 U-134
valueFraction keyword, U-136 wall
van Leer differencing, P-34 keyword entry, U-135
vanDriestDelta model, U-101 wallFunctionTable utility, U-89
vanLeer wallGradU utility, U-92
keyword entry, U-116 wallHeatFlux utility, U-92
VCR Controls menu, U-25, U-165 Wallis
vector, P-12 library, U-99
operator, P-23 wallShearStress utility, U-92
unit, P-16 wclean script/alias, U-73
vector class, P-19, U-107 wdot utility, U-93
vector product, see tensor, vector cross wedge
product boundary condition, U-130, U-135, U-
vectorField class, P-25 144
version keyword, U-105 wedge
vertices keyword, U-20, U-138, U-139 keyword entry, U-135
veryInhomogeneousMixture model, U-98, window
U-180 Color Legend, U-27
View menu, U-167 Options, U-168
View Settings Pipeline Browser, U-23, U-164
menu entry, U-24, U-167 Render View, U-168
View Settings (Render View) window, U- Seed, U-169
167 View Settings (Render View), U-167
View Settings... window panel
Animations, U-168 WM PROJECT USER DIR
Annotation, U-24, U-167 environment variable, U-74
Charts, U-168 WM PROJECT VERSION
Color Legend, U-166 environment variable, U-74
Color Scale, U-166 wmake
Colors, U-168 platforms, U-70
Display, U-23, U-25, U-164, U-165 wmake script/alias, U-69
General, U-167, U-168 word class, P-21, P-27
Information, U-164 writeCellCentres utility, U-93
Lights, U-167 writeCompression keyword, U-112
Mesh Parts, U-23 writeControl
Properties, U-25, U-164 keyword entry, U-111
Render View, U-168 writeControl keyword, U-22, U-60, U-111
Style, U-23, U-166 writeFormat keyword, U-54, U-112
Wireframe writeInterval keyword, U-22, U-31, U-111
menu entry, U-166 writeMeshObj utility, U-90
WM ARCH writeNow
environment variable, U-74 keyword entry, U-111
WM ARCH OPTION writePrecision keyword, U-112
environment variable, U-74
WM COMPILE OPTION
X
environment variable, U-74 x
WM COMPILER keyword entry, U-175
environment variable, U-74 XiFoam solver, U-86
WM COMPILER BIN xmgr
environment variable, U-74 keyword entry, U-112, U-174
WM COMPILER DIR xyz
WM COMPILER LIB
Y
WM DIR y
WM MPLIB yPlusLES utility, U-92
environment variable, U-74 yPlusRAS utility, U-92
WM OPTIONS
Z
WM PRECISION OPTION z
WM PROJECT zeroGradient
environment variable, U-74 boundary condition, U-136
WM PROJECT DIR zipUpMesh utility, U-91
WM PROJECT INST DIR

OpenFOAM编程指南

Uploaded by

Copyright:

Available Formats

OpenFOAM编程指南

Uploaded by

Document Information

Copyright

Available Formats

Share this document

Share or Embed Document

Sharing Options

Did you find this document useful?

Is this content inappropriate?

Copyright:

Available Formats

OpenFOAM编程指南

Uploaded by

Copyright:

Available Formats

Op

Copyright©2011-2014 OpenFOAM Foundation.

This work is licensed under a Creative Commons Attribution - NonCommercial

版权声明©2011-2014 OpenFOAM Foundation

本指南遵守：Creative Commons Attribution - NonCommercial NoDerivs 3.0

b. "Collection" means a collection of literary or artistic works, such as encyclopedias

selection and arrangement of their contents, constitute intellectual creations, in which

h. "Publicly Perform" means to perform public recitations of the Work and to

b. to Distribute and Publicly Perform the Work including as incorporated in Collections.

d. For the avoidance of doubt:

ii. Waivable Compulsory License Schemes. In those jurisdictions in which the

License if Your exercise of such rights is for a purpose or use which is

e. Except as otherwise agreed in writing by the Licensor or as may be otherwise

5. Representations, Warranties and Disclaimer

UNLESS OTHERWISE MUTUALLY AGREED BY THE PARTIES IN WRITING, LICENSOR

6. Limitation on Liability. EXCEPT TO THE EXTENT REQUIRED BY APPLICABLE LAW,

b. If any provision of this License is invalid or unenforceable under applicable law, it

c. No term or provision of this License shall be deemed waived and no breach

《OpenFOAM 编程指南中文版》由李东岳 1 翻译第一、二章，高翔宇 2 翻译第三章

1 本手册相对于 OpenFOAM ProgrammersGuide 对行距有所调整以保证页数匹配

 维度 d：我们用 d 来表示，在 OpenFOAM 中 d = 3；

 1 阶矢量： 它用来表示具有大小和方向的物理量。他们通常的形式为： a = (𝑎1 , 𝑎2 , 𝑎3 )，

 2 阶张量：我们用 T 来表示，T 具有 9 个元素，参见如下：

𝑇11 𝑇12 𝑇13

由于其为 2 阶张量，因此分量𝑇𝑖𝑗 使用两个索引符来表示，如上述，我们在下文省略

𝑇11 𝑇21 𝑇31

 3 阶张量：具有 27 个元素，我们用𝑃𝑖𝑗𝑘 来表示，由于其太复杂在这里我们不表示；

 对称 3 阶张量：OpenFOAM 将其定义为下述元素相等的 3 阶张量：𝑃𝑖𝑗𝑘 = 𝑃𝑖𝑘𝑗 =

在这一章节我们描述 OpenFOAM 可用的张量算术。我们首先来看一下最常见的张量操作：

s𝐚 = 𝑠𝑎𝑖 = (𝑠𝑎1 , 𝑠𝑎2 , 𝑠𝑎3 ) (1.5)

在做张量 a 和标量 s 的除法的时候，标量 s 必须为分母，例如：

𝐚/𝑠 = 𝒂𝑖 /𝑠 = (𝑎1 /𝑠, 𝑎2 /𝑠, 𝑎3 /𝑠) (1.6)

张量的内积可以在阶数为𝑟1 的张量和阶数为𝑟2 的张量上进行操作。结果张量的阶数为𝑟 =

 张量 T 和矢量 a 的内积结果为矢量 b，我们用列向量来表示𝐛 = 𝐓 ∙ 𝐚（看起来更加

𝑇11 𝑎1 + 𝑇12 𝑎2 + 𝑇13 𝑎3

𝑇11 𝑎1 + 𝑇21 𝑎2 + 𝑇31 𝑎3

 两个张量 T 和 S 的内积也是一个张量，即 P = T ∙ S，P 元素的值这样计算：

𝑃𝑖𝑗 = 𝑇𝑖𝑘 𝑆𝑘𝑗 (1.10)

 矢量 a 和三阶张量 P 的内积结果为一个二阶张量 T，T = a ∙ P，T 的元素为：

𝑇𝑖𝑗 = 𝑎𝑘 𝑃𝑘𝑖𝑗 (1.11)

𝑇𝑖𝑗 = 𝑃𝑖𝑗𝑘 𝑎𝑘 (1.12)

 二阶张量 T 和三阶张量 P 的内积为三阶张量 Q，Q = T ∙ P，Q 的元素为：

𝑄𝑖𝑗𝑘 = 𝑇𝑖𝑙 𝑃𝑙𝑗𝑘 (1.13)

𝑄𝑖𝑗𝑘 = 𝑃𝑖𝑗𝑙 𝑇𝑙𝑘 (1.14)

两个二阶张量 T 和 S 的双内积结果为一个标量 s（s = T : S），s 这样计算：

𝑠 = 𝑇𝑖𝑗 𝑆𝑖𝑗 = 𝑇11 𝑆11 + 𝑇12 𝑆12 + 𝑇13 𝑆13 + (1.15)

二阶张量 T 和三阶张量 P 的双内积结果为一个矢量 a（a = T : P）:

𝑎𝑖 = 𝑇𝑗𝑘 𝑃𝑗𝑘𝑖 (1.16)

𝑎𝑖 = 𝑃𝑖𝑗𝑘 𝑇𝑗𝑘 (1.17)

两个三阶张量 P 和 Q 的三内积结果为一个标量 s（s = P ⋮ Q），s 这样计算：

𝑠 = 𝑃𝑖𝑗𝑘 𝑄𝑖𝑗𝑘 (1.18)

 两个矢量 a 和 b 的外积是不可交换的，结果为张量T = ab = (ba)T ，其这样计算：

 一个矢量 a 和一个张量 T 的外积为一个三阶张量 P（P = a T），其这样计算：

𝑃𝑖𝑗𝑘 = 𝑎𝑖 𝑇𝑗𝑘 (1.20)

𝑃𝑖𝑗𝑘 = 𝑇𝑗𝑘 𝑎𝑖 (1.21)

叉乘只能用于矢量，对于两个矢量 a 和 b，他们的叉乘为矢量 c，其这样计算：

𝑐 = 𝑎 × 𝑏 = 𝑒𝑖𝑗𝑘 𝑎𝑗 𝑏𝑘 = (𝑎2𝑏3 − 𝑎3 𝑏2 , 𝑎3 𝑏1 − 𝑎1 𝑏3 , 𝑎1 𝑏2 − 𝑎2 𝑏1 ) (1.22)

下面我们描述 OpenFOAM 提供的一些不太常用的张量操作：

 平方：某个张量和自己本身的外积，例如对于矢量 a 的平方为：𝒂2 = 𝒂𝒂；

𝑐𝑖 = scale(𝑎, 𝑏) = (𝑎1 𝑏1 , 𝑎2 𝑏2 , 𝑎3 𝑏3 ) (1.24)

 1 阶矢量：它用来表示具有大小和方向的物理量。他们通常的形式为： a = (𝑎1 , 𝑎2 , 𝑎3 )，

2 阶张量可以定义为矢量的线性函数，例如：向量 a 和向量 b 可以通过 T 联系起来𝐚 = 𝐓 ⋅ 𝐛。

操作描述数学形式 OpenFOAM 中的实现

离散项显性/隐性数学形式代码实现