定义加速度的phi

danao

重力g的phi可以写成phiG=g & mesh.Sf(),其中g是uniformDimensionedVectorField。如果g变成一个空间上不均匀的加速度场，类型变成volVectorField，这时候的加速度的phi该如何写？？

李东岳

那不就跟速度U一样了么？也是这样写啊？

wwzhao

@东岳 Sf是surfaceVectorField，不能跟volVectorField直接相乘，需要将volVectorField插值到面上再相乘。

李东岳

@wwzhao 主楼知道的吧?

danao

@wwzhao 那应该怎么写啊啊。。。

danao

@东岳 @wwzhao ，抱歉，忘了格式了

volVectorField g
(
IOobject
(
"g",
runTime.timeName(),
mesh,
IOobject::NO_READ,
IOobject::NO_WRITE
),
g=关于网格空间坐标xy的函数
);
dimensionedVectorField g = fvc::interpolate(g,"g");
phig=g & mesh.Sf()

李东岳

volVectorField g
(
IOobject
(
"g",
runTime.timeName(),
mesh,
IOobject::NO_READ,
IOobject::NO_WRITE
),
g=关于网格空间坐标xy的函数
);
//dimensionedVectorField g = fvc::interpolate(g,"g");
phig = fvc::interpolate(g) & mesh.Sf()

你这是做什么？怎么会有非均一重力场

danao

@东岳非惯性坐标系啊，体积力与每个网格的具体位置有关，这样就不用动网格了