openFoam:报告湍流场时遇到问题
我运行openFoam 4.1,在k-eps模拟中提取湍流应力张量R时遇到问题。4天前我在cfd-online.com上问了同样的问题,但还没有收到回复,所以我现在希望有人能在这里帮助我 为了提取R,我只修改了controlDict,如下所示:openFoam:报告湍流场时遇到问题,openfoam,Openfoam,我运行openFoam 4.1,在k-eps模拟中提取湍流应力张量R时遇到问题。4天前我在cfd-online.com上问了同样的问题,但还没有收到回复,所以我现在希望有人能在这里帮助我 为了提取R,我只修改了controlDict,如下所示: /*--------------------------------*- C++ -*----------------------------------*\ | ========= |
/*--------------------------------*- C++ -*----------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: 4.0 |
| \\ / A nd | Web: www.OpenFOAM.org |
| \\/ M anipulation | |
\*---------------------------------------------------------------------------*/
FoamFile
{
version 2.0;
format ascii;
class dictionary;
object controlDict;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
application simpleFoam;
startTime 0;
stopAt endTime;
endTime 5000;
deltaT 1;
writeControl timeStep;
writeInterval 100;
purgeWrite 2;
writeFormat ascii;
writePrecision 8;
writeCompression uncompressed;
timeFormat general;
timePrecision 6;
runTimeModifiable yes;
functions
{
#includeFunc residuals
turbulenceFields1
{
type turbulenceFields;
libs ("libfieldFunctionObjects.so");
field R;
}
}
// ************************************************************************* //
/*---------------------------------------------------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: 4.1 |
| \\ / A nd | Web: www.OpenFOAM.org |
| \\/ M anipulation | |
\*---------------------------------------------------------------------------*/
Build : 4.1
Exec : simpleFoam
Date : May 18 2017
Time : 07:08:01
Host : "..."
PID : 11073
Case : /home...
nProcs : 1
sigFpe : Enabling floating point exception trapping (FOAM_SIGFPE).
fileModificationChecking : Monitoring run-time modified files using timeStampMaster
allowSystemOperations : Allowing user-supplied system call operations
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
Create time
Create mesh for time = 0
SIMPLE: convergence criteria
field p tolerance 0.001
field U tolerance 0.001
field "(k|epsilon|omega|f|v2)" tolerance 0.001
Reading field p
Reading field U
Reading/calculating face flux field phi
Selecting incompressible transport model Newtonian
Selecting turbulence model type RAS
Selecting RAS turbulence model kEpsilon
kEpsilonCoeffs
{
Cmu 0.09;
C1 1.44;
C2 1.92;
C3 -0.33;
sigmak 1;
sigmaEps 1.3;
}
No MRF models present
No finite volume options present
Starting time loop
turbulenceFields turbulenceFields1: storing fields:
turbulenceProperties:R
Time = 1
smoothSolver: Solving for Ux, Initial residual = 1, Final residual = 0.068919891, No Iterations 6
smoothSolver: Solving for Uy, Initial residual = 1, Final residual = 0.094423261, No Iterations 4
smoothSolver: Solving for Uz, Initial residual = 1, Final residual = 0.065365879, No Iterations 6
GAMG: Solving for p, Initial residual = 1, Final residual = 0.0086113464, No Iterations 26
time step continuity errors : sum local = 0.51812432, global = -0.07783965, cumulative = -0.07783965
smoothSolver: Solving for epsilon, Initial residual = 0.59675441, Final residual = 0.051409907, No Iterations 3
smoothSolver: Solving for k, Initial residual = 1, Final residual = 0.092548244, No Iterations 2
ExecutionTime = 0.81 s ClockTime = 1 s
--> FOAM FATAL ERROR:
object of type N4Foam9Function1INS_10SymmTensorIdEEEE is not allocated
From function T* Foam::autoPtr<T>::operator->() [with T = Foam::Function1<Foam::SymmTensor<double> >]
in file /home/ubuntu/OpenFOAM/OpenFOAM-4.1/src/OpenFOAM/lnInclude/autoPtrI.H at line 176.
FOAM aborting
#0 Foam::error::printStack(Foam::Ostream&) at ??:?
#1 Foam::error::abort() at ??:?
#2 Foam::timeVaryingMappedFixedValueFvPatchField<Foam::SymmTensor<double> >::write(Foam::Ostream&) const at ??:?
#3 Foam::GeometricField<Foam::SymmTensor<double>, Foam::fvPatchField, Foam::volMesh>::Boundary::writeEntry(Foam::word const&, Foam::Ostream&) const at ??:?
#4 Foam::GeometricField<Foam::SymmTensor<double>, Foam::fvPatchField, Foam::volMesh>::writeData(Foam::Ostream&) const at ??:?
#5 Foam::regIOobject::writeObject(Foam::IOstream::streamFormat, Foam::IOstream::versionNumber, Foam::IOstream::compressionType) const at ??:?
#6 Foam::functionObjects::regionFunctionObject::writeObject(Foam::word const&) at ??:?
#7 Foam::functionObjects::turbulenceFields::write() at ??:?
#8 Foam::functionObjectList::execute() at ??:?
#9 Foam::Time::loop() at ??:?
#10 Foam::simpleControl::loop() at ??:?
#11 ? at ??:?
#12 __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6"
#13 ? at ??:?
Aborted (core dumped)
当我尝试使用上述controlDict进行求解时,我得到以下结果:
/*--------------------------------*- C++ -*----------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: 4.0 |
| \\ / A nd | Web: www.OpenFOAM.org |
| \\/ M anipulation | |
\*---------------------------------------------------------------------------*/
FoamFile
{
version 2.0;
format ascii;
class dictionary;
object controlDict;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
application simpleFoam;
startTime 0;
stopAt endTime;
endTime 5000;
deltaT 1;
writeControl timeStep;
writeInterval 100;
purgeWrite 2;
writeFormat ascii;
writePrecision 8;
writeCompression uncompressed;
timeFormat general;
timePrecision 6;
runTimeModifiable yes;
functions
{
#includeFunc residuals
turbulenceFields1
{
type turbulenceFields;
libs ("libfieldFunctionObjects.so");
field R;
}
}
// ************************************************************************* //
/*---------------------------------------------------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: 4.1 |
| \\ / A nd | Web: www.OpenFOAM.org |
| \\/ M anipulation | |
\*---------------------------------------------------------------------------*/
Build : 4.1
Exec : simpleFoam
Date : May 18 2017
Time : 07:08:01
Host : "..."
PID : 11073
Case : /home...
nProcs : 1
sigFpe : Enabling floating point exception trapping (FOAM_SIGFPE).
fileModificationChecking : Monitoring run-time modified files using timeStampMaster
allowSystemOperations : Allowing user-supplied system call operations
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
Create time
Create mesh for time = 0
SIMPLE: convergence criteria
field p tolerance 0.001
field U tolerance 0.001
field "(k|epsilon|omega|f|v2)" tolerance 0.001
Reading field p
Reading field U
Reading/calculating face flux field phi
Selecting incompressible transport model Newtonian
Selecting turbulence model type RAS
Selecting RAS turbulence model kEpsilon
kEpsilonCoeffs
{
Cmu 0.09;
C1 1.44;
C2 1.92;
C3 -0.33;
sigmak 1;
sigmaEps 1.3;
}
No MRF models present
No finite volume options present
Starting time loop
turbulenceFields turbulenceFields1: storing fields:
turbulenceProperties:R
Time = 1
smoothSolver: Solving for Ux, Initial residual = 1, Final residual = 0.068919891, No Iterations 6
smoothSolver: Solving for Uy, Initial residual = 1, Final residual = 0.094423261, No Iterations 4
smoothSolver: Solving for Uz, Initial residual = 1, Final residual = 0.065365879, No Iterations 6
GAMG: Solving for p, Initial residual = 1, Final residual = 0.0086113464, No Iterations 26
time step continuity errors : sum local = 0.51812432, global = -0.07783965, cumulative = -0.07783965
smoothSolver: Solving for epsilon, Initial residual = 0.59675441, Final residual = 0.051409907, No Iterations 3
smoothSolver: Solving for k, Initial residual = 1, Final residual = 0.092548244, No Iterations 2
ExecutionTime = 0.81 s ClockTime = 1 s
--> FOAM FATAL ERROR:
object of type N4Foam9Function1INS_10SymmTensorIdEEEE is not allocated
From function T* Foam::autoPtr<T>::operator->() [with T = Foam::Function1<Foam::SymmTensor<double> >]
in file /home/ubuntu/OpenFOAM/OpenFOAM-4.1/src/OpenFOAM/lnInclude/autoPtrI.H at line 176.
FOAM aborting
#0 Foam::error::printStack(Foam::Ostream&) at ??:?
#1 Foam::error::abort() at ??:?
#2 Foam::timeVaryingMappedFixedValueFvPatchField<Foam::SymmTensor<double> >::write(Foam::Ostream&) const at ??:?
#3 Foam::GeometricField<Foam::SymmTensor<double>, Foam::fvPatchField, Foam::volMesh>::Boundary::writeEntry(Foam::word const&, Foam::Ostream&) const at ??:?
#4 Foam::GeometricField<Foam::SymmTensor<double>, Foam::fvPatchField, Foam::volMesh>::writeData(Foam::Ostream&) const at ??:?
#5 Foam::regIOobject::writeObject(Foam::IOstream::streamFormat, Foam::IOstream::versionNumber, Foam::IOstream::compressionType) const at ??:?
#6 Foam::functionObjects::regionFunctionObject::writeObject(Foam::word const&) at ??:?
#7 Foam::functionObjects::turbulenceFields::write() at ??:?
#8 Foam::functionObjectList::execute() at ??:?
#9 Foam::Time::loop() at ??:?
#10 Foam::simpleControl::loop() at ??:?
#11 ? at ??:?
#12 __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6"
#13 ? at ??:?
Aborted (core dumped)
/*---------------------------------------------------------------------------*\
| ========= | |
|\\/F ield | OpenFOAM:开源CFD工具箱|
|\\/O操作版本:4.1|
|\\/A nd|网址:www.OpenFOAM.org|
|\\/M操作|
\*---------------------------------------------------------------------------*/
构建:4.1
执行官:SimpleFam
日期:2017年5月18日
时间:07:08:01
主持人:“…”
PID:11073
案例:/home。。。
全国人大常委会:1
sigFpe:启用浮点异常捕获(FOAM_sigFpe)。
fileModificationChecking:使用timeStampMaster监视运行时修改的文件
allowSystemOperations:允许用户提供的系统调用操作
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
创造时间
为时间=0创建网格
简单:收敛准则
现场p公差0.001
现场U公差0.001
字段“(k |ε|ω| f | v2)”公差0.001
读取字段p
读域U
读取/计算面通量场φ
牛顿流体不可压缩输运模型的选择
选择湍流模型类型RAS
选择RAS湍流模型kEpsilon
凯普西隆系数
{
债务工具中央结算系统0.09;
C1 1.44;
C2 1.92;
C3-0.33;
sigmak1;
sigmaps1.3;
}
没有MRF型号
不存在有限体积选项
开始时间循环
湍流字段湍流字段1:存储字段:
湍流特性:R
时间=1
smoothSolver:求解Ux,初始残差=1,最终残差=0.068919891,无迭代6
smoothSolver:求解Uy,初始残差=1,最终残差=0.094423261,无迭代4
平滑解算器:求解Uz,初始残差=1,最终残差=0.065365879,无迭代6
GAMG:求解p,初始残差=1,最终残差=0.0086113464,无迭代26
时间步长连续性错误:局部和=0.51812432,全局=-0.07783965,累积=-0.07783965
smoothSolver:求解ε,初始残差=0.59675441,最终残差=0.051409907,无迭代3
smoothSolver:求解k,初始残差=1,最终残差=0.092548244,无迭代2
执行时间=0.81秒时钟时间=1秒
-->泡沫致命错误:
未分配N4Foam9Function1INS\u 10SymmTensorIdeee类型的对象
来自函数T*Foam::autoPtr::operator->()[带T=Foam::Function1]
在第176行的文件/home/ubuntu/OpenFOAM/OpenFOAM-4.1/src/OpenFOAM/lnInclude/autoPtrI.H中。
泡沫流产
#0泡沫::错误::打印堆栈(泡沫::Ostream&)位于?
#1泡沫::错误::中止()位于?
#2泡沫::timeVaryingMappedFixedValueFvPatchField::写入(泡沫::Ostream&)常数在?
#3泡沫::几何字段::边界::写入尝试(泡沫::单词常量和,泡沫::Ostream&)常量在?
#4泡沫::几何字段::写入数据(泡沫::Ostream&)常数在?
#5 Foam::regIOobject::writeObject(Foam::IOstream::streamFormat,Foam::IOstream::versionNumber,Foam::IOstream::compressionType)常量位于?
#6 Foam::functionObjects::regionFunctionObject::writeObject(Foam::word const&)at?
#7 Foam::functionObjects::turbulenceFields::write()位于?
#8 Foam::functionObjectList::execute()位于?
#9泡沫::时间::循环()在?
#10泡沫::simpleControl::loop()位于?
#11 ? 在?
#12“/lib/x86\u 64-linux-gnu/libc.so.6”中的\u libc\u start\u main
#13 ? 在?
中止(堆芯转储)
从上面看,似乎openFoam第一次解决了这个问题,然后出现了一些问题。在目录“1”中,现在只有一个文件“turbulenceProperties:R”,它包含每个点的6个应力张量元素
有人知道怎么解决这个问题吗 我想我在cfd在线论坛的链接中找到了答案