OpenFOAM计入曲率影响的速度边界条件

pengdi

在partialSlip的基础上创建计入曲率影响的速度边界，但是计算fvc::div(n)后，边界条件可以编译，但是无法进行计算。具体边界条件如下：
tmp<vectorField> nHat = this->patch().nf();

const fvMesh& mesh = this->internalField().mesh();
const surfaceVectorField& Sf = mesh.Sf();  
const surfaceScalarField& magSf = mesh.magSf();             // Face areas
const surfaceVectorField& N = Sf/(magSf); 

Field<Type>::operator=
(
    (1/(1+((1/(this->patch().deltaCoeffs()))/valueFraction_)+ ((mag(fvc::div(N)))*(1/(this->patch().deltaCoeffs()))))) 
      *transform(I - sqr(nHat), this->patchInternalField())
);

transformFvPatchField<Type>::evaluate();

bestucan

倒是贴下报错什么的

pengdi

此回复已被删除！

pengdi

@bestucan
错误如下：
in file D:/blue-CFD/blueCFD-Core-2017/OpenFOAM-5.x/src/OpenFOAM/fields/Fields/Field/Field.T.C at line 324
Reading "E:/jiyi/pspsps/0/U.boundaryField.FOIL" from line 42 to line 45
expected keyword 'uniform' or 'nonuniform', assuming deprecated Field format from Foam version 2.0.
Generating stack trace...

bestucan

@pengdi
http://bugs.openfoam.org/view.php?id=471
看看是不是因为这个

CFD中文网

OpenFOAM计入曲率影响的速度边界条件