CFD软件：OpenFOAM二次开发_（21）.OpenFOAM高级应用与研究.docx

PAGE1

PAGE1

OpenFOAM高级应用与研究

1.高级求解器开发

1.1.自定义求解器的基本步骤

在OpenFOAM中，开发自定义求解器是一个重要的高级应用。通过自定义求解器，可以针对特定的流体动力学问题进行优化和扩展。以下是开发自定义求解器的基本步骤：

创建新的求解器目录：

在OpenFOAM的$FOAM_TUTORIALS目录中，创建一个新的求解器目录。例如，创建一个名为customSolver的目录。

mkdir-p$FOAM_TUTORIALS/incompressible/customSolver

复制标准求解器模板：

从OpenFOAM的标准求解器中复制一个模板到新的求解器目录。例如，复制icoFoam求解器的模板。

cp-r$FOAM_APP/incompressible/icoFoam$FOAM_TUTORIALS/incompressible/customSolver

重命名求解器文件：

将模板中的求解器文件重命名为自定义求解器的名称。例如，将icoFoam.C重命名为customSolver.C。

mv$FOAM_TUTORIALS/incompressible/customSolver/icoFoam.C$FOAM_TUTORIALS/incompressible/customSolver/customSolver.C

修改求解器文件：

打开customSolver.C文件，根据需要进行修改。例如，添加新的物理模型、边界条件或数值方法。

//customSolver.C

#includefvCFD.H

intmain(intargc,char*argv[])

{

#includesetRootCaseLists.H

#includecreateTime.H

#includecreateMesh.H

#includecreateFields.H

#includeinitContinuityErrs.H

//添加新的物理模型

volScalarFieldalpha

(

IOobject

(

alpha,

runTime.timeName(),

mesh,

IOobject::MUST_READ,

IOobject::AUTO_WRITE

),

mesh

);

//主时间循环

while(runTime.loop())

{

#includereadPISOControls.H

#includesolveMomentum.H

#includesolvePressure.H

#includecontinuityErrs.H

#includewrite.H

//添加新的求解步骤

solve(fvm::ddt(alpha)+fvm::div(phi,alpha)==0);

}

InfoEnd\nendl;

return0;

}

创建Make文件：

在customSolver目录中创建Make文件，指定求解器的源文件和依赖库。

#Make/FILES

customSolver.o

#Make/options

EXE=$(FOAM_USER_APPBIN)/customSolver

LIB=\

-lincompressibleTransportModels\

-lincompressibleTurbulenceModel

编译求解器：

使用wmake命令编译新的求解器。

cd$FOAM_TUTORIALS/incompressible/customSolver

wmake

1.2.添加新的物理模型

在OpenFOAM中，可以添加新的物理模型来扩展求解器的功能。例如，添加一个新的标量场alpha来模拟某种物理现象。

定义新的标量场：

在求解器文件中定义新的标量场alpha。

volScalarFieldalpha

(

IOobject

(

alpha,

runTime.timeName(),