L
@bike-北辰 在 OpenFOAM v2012 waves2foam 造波问题 中说：
@liujm 您好，waves2foam 的版本应该是2148，编译前手动下载了OceanWave3D-Fortran90，然后修改了一下./Allwmake文件，具体可以参考(https://blog.csdn.net/Killig_1/article/details/135386466) ，然后正常编译即可。
嗯嗯，以前编译成功过，可是最近重新编译报错严重，确实主要问题出现在OceanWave3D-Fortran90这里，需要注释掉make中git clone和检查更新部分，感谢~
@李东岳 静态网格也是有这个问题的。按理说不应该啊...我再跑跑其他类似的原生算例
p：
~
#include "../varsSettings" //Include global variables
dimensions [0 2 -2 0 0 0 0];
internalField uniform 0;
boundaryField
{
inlet
{
type zeroGradient;
}
outlet
{
type fixedValue;
value $internalField;
}
frontAndBack
{
type empty;//zeroGradient;
}
topAndBottom
{
type zeroGradient;
}
frame
{
type zeroGradient;//fixedFluxPressure;
}
}
U：
#include "../varsSettings" //Include global variables
dimensions [0 1 -1 0 0 0 0];
internalField uniform $flowVelocity;
boundaryField
{
inlet
{
type fixedValue;
value $internalField;
}
outlet
{
type inletOutlet;
inletValue uniform (0 0 0);
value $internalField;
}
frontAndBack
{
type empty;//slip;
//value uniform (0 0 0);
}
topAndBottom
{
type slip;
value uniform (0 0 0);
}
frame
{
type noSlip;//movingWallVelocity;
//value uniform (0 0 0);
}
}
算例文件 OpenfOAM2012 网格数4w
functionalTest_NoPZRefine2_2D_staticMesh.zip
为了简化算例，我把网格三维转化为二维，简化了结构形式，以便排错。
【算例云图】
【简化版算例下载】
functionalTest_NoPZRefine2_2D.zip
生成网格createMesh，
跑算例runCaseOnly，
全流程runCase
感谢~
各位老师好，
【问题描述】
本人在使用pimpleFoam模拟固定式网框在自由来流（u=0.5m/s）中的响应。整个流场的计算结果符合预期，可是网框表面处的流速竟然高达1e2量级，明显与现实不符，结果云图如下：（侧视图）
【分析】
网框截面原本比本帖中展示的要更小，起初正方形边长仅分布四个网格，为了排除网格数量对结果的影响，特意扩大了截面，增加边长网格数量。但是现在依然存在这个问题。
【初始条件设置】
以下摘抄0文件夹设置：
p：
/*--------------------------------*- C++ -*----------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: v2012 |
| \\ / A nd | Website: www.openfoam.com |
| \\/ M anipulation | |
\*---------------------------------------------------------------------------*/
FoamFile
{
version 2.0;
format ascii;
class volScalarField;
location "0";
object p;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
#include "../varsSettings" //Include global variables
dimensions [0 2 -2 0 0 0 0];
internalField uniform 0;
boundaryField
{
inlet
{
type zeroGradient;
}
outlet
{
type fixedValue;
value $internalField;
}
frontAndBack
{
type zeroGradient;
}
topAndBottom
{
type zeroGradient;
}
frame
{
type fixedFluxPressure;
}
}
// ************************************************************************* //
U:
/*--------------------------------*- C++ -*----------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: v2012 |
| \\ / A nd | Website: www.openfoam.com |
| \\/ M anipulation | |
\*---------------------------------------------------------------------------*/
FoamFile
{
version 2.0;
format ascii;
class volVectorField;
location "0";
object U;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
#include "../varsSettings" //Include global variables
dimensions [0 1 -1 0 0 0 0];
internalField uniform $flowVelocity;
boundaryField
{
inlet
{
type fixedValue;
value $internalField;
}
outlet
{
type inletOutlet;
inletValue uniform (0 0 0);
value $internalField;
}
frontAndBack
{
type slip;
value uniform (0 0 0);
}
topAndBottom
{
type slip;
value uniform (0 0 0);
}
frame
{
/* type noSlip; */
type movingWallVelocity;
value uniform (0 0 0);
}
}
// ************************************************************************* //
pointDisplacement:
/*--------------------------------*- C++ -*----------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: 5 |
| \\ / A nd | Web: www.OpenFOAM.org |
| \\/ M anipulation | |
\*---------------------------------------------------------------------------*/
FoamFile
{
version 2.0;
format ascii;
class pointVectorField;
object pointDisplacement;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
dimensions [0 1 0 0 0 0 0];
internalField uniform (0 0 0);
boundaryField
{
inlet
{
type fixedValue;
value uniform (0 0 0);
}
outlet
{
type fixedValue;
value uniform (0 0 0);
}
frontAndBack
{
type fixedValue;
value uniform (0 0 0);
}
topAndBottom
{
type fixedValue;
value uniform (0 0 0);
}
frame
{
type calculated;
value uniform (0 0 0);
}
}
// ************************************************************************* //
【算例设置】
functionalTest_NoPZRefine1.zip
生成网格createMesh，
跑算例runCaseOnly，
全流程runCase
希望各位老师能给出建议~@李东岳
@Joker
您好，请问您使用的OF是什么版本呢？我在用displacementLaplacian后，使用moveDynamicMesh预览网格运动的时候，系统报错，似乎没有displacementLaplacian这个求解器。请问是我的字典设置有误吗？
Create time
Create mesh for time = 0
Selecting dynamicFvMesh dynamicMotionSolverFvMesh
Selecting motion solver: displacementLaplacian
--> FOAM FATAL IO ERROR: (openfoam-2012)
Unknown solver type displacementLaplacian
Valid solver types :
7(coded displacementInterpolation displacementLayeredMotion multiSolidBodyMotionSolver sixDoFRigidBodyMotion solidBody velocityDisplacement)
file: /OpenFOAM-v2012/myProject/tutorials/floatingPorousPimpleFoam/Condition1_Uin=1.2/constant/dynamicMeshDict at line 18 to 26.
From static Foam::autoPtr<Foam::motionSolver> Foam::motionSolver::New(const Foam::polyMesh&, const Foam::IOdictionary&)
in file motionSolvers/motionSolver/motionSolver.C at line 137.
FOAM exiting
各位前辈大家好，
我在使用const vector cellCentre = mesh_.C()[celli]收集网格中心坐标时，编译器出现报错：
In file included from cfdTools/general/porosityModel/floatingDarcyForchheimer/floatingDarcyForchheimer.H:216,
from cfdTools/general/porosityModel/floatingDarcyForchheimer/floatingDarcyForchheimer.C:30:
cfdTools/general/porosityModel/floatingDarcyForchheimer/floatingDarcyForchheimerTemplates.C: In member function ‘void Foam::porosityModels::floatingDarcyForchheimer::applyInCalcForce(Foam::scalarField&, Foam::vectorField&, const scalarField&, const RhoFieldType&, const scalarField&, const vectorField&, Foam::vectorField&) const’:
cfdTools/general/porosityModel/floatingDarcyForchheimer/floatingDarcyForchheimerTemplates.C:119:48: error: no match for ‘operator[]’ (operand types are ‘const volVectorField’ {aka ‘const Foam::GeometricField<Foam::Vector<double>, Foam::fvPatchField, Foam::volMesh>’} and ‘const label’ {aka ‘const int’})
119 | const vector cellCentre = mesh_.C()[celli];
|
貌似mesh_.C()返回的数据类型不支持使用[]符号，可是为什么在其他代码中就可以这样使用 mesh_.C()[celli]呢？
后期，我通过使用const vector cellCentre = mesh_.C().primitiveField()[celli]通过编译。不过尚不清楚其与mesh_.C()返回的结果是否存在差异。
烦请各位老师解答~
@李东岳 李老师，这也正是我遇到的问题
。
退一步我做一个简化，在图示的myForces.H的结构最外层cell上套一层patch边界，可不可以仅用patch计算cellZone最外层网格的力和力矩[1]，而不阻碍水体穿越被patch包裹的cellZone区域？
这个想法可行吗，如果可行的话老师方便指点一下，如何实现水体不穿越patch的功能？
[1] 只计算外层网格表面的受力是因为网衣通常不会太厚，有文献采用2m厚的多孔介质模型，网格尺寸0.5m来模拟。
在OF2206中，
对于patch，有addToPatchFields函数：
void Foam::functionObjects::forces::addToPatchFields
(
const label patchi,
const vectorField& Md,
const vectorField& fP,
const vectorField& fV
)
{
sumPatchForcesP_ += sum(fP);
sumPatchForcesV_ += sum(fV);
force().boundaryFieldRef()[patchi] += fP + fV;
const vectorField mP(Md^fP);
const vectorField mV(Md^fV);
sumPatchMomentsP_ += sum(mP);
sumPatchMomentsV_ += sum(mV);
moment().boundaryFieldRef()[patchi] += mP + mV;
}
对于多孔介质的cellZone，有addToInternalField函数：
void Foam::functionObjects::forces::addToInternalField
(
const labelList& cellIDs,
const vectorField& Md,
const vectorField& f
)
{
auto& force = this->force();
auto& moment = this->moment();
forAll(cellIDs, i)
{
const label celli = cellIDs[i];
sumInternalForces_ += f[i];
force[celli] += f[i];
const vector m(Md[i]^f[i]);
sumInternalMoments_ += m;
moment[celli] = m;
}
}
