CFD中文网

Lisiye

各CFDer:

最近做了些高精度高分辨率格式的工作，用Shu-Osher、double Mach reflection和 Rayleigh-Taylor做了下验证。

还有其他的一些比较好的算例可以请教一下吗，感谢！

祝好

Lisiye

Dear foamers
在自适应网格运动中，发现了this ->指针的一些运用(如下图)，想请教一下this的成员包含哪些，各成员的具体含义是什么呢，我可以去哪里查询?Good luck!

li
2020.11

Lisiye

@史浩您好，请问对于自适应网格的设置讲解有什么资料可以帮助我学习参数设置的吗

Lisiye

好的，谢谢您

Lisiye

您好感谢您的指导，现在火箭的简单模型可以通过一定速度运动了 0/pointDisplacement 如下

FoamFile
{
version 2.0;
format ascii;
class pointVectorField;
location "0";
object pointDisplacement;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

dimensions [0 1 0 0 0 0 0];

internalField uniform (0 0 0);

boundaryField
{
WALL
{
type fixedValue;
value uniform (0 0 0);
}

GROUND
{
type            fixedValue;
value		uniform (0 0 0);
}

OUTLET
{
type            fixedValue;
value		uniform (0 0 0);
}
    
frontAndBackPlanes
{
    type            fixedValue;
    value		uniform (0 0 0);
}

ROCKET_WALL
{
    type            solidBodyMotionDisplacement;
    solidBodyMotionFunction     linearMotion;
velocity  (0 2 0);

}

INLET
{
    type            solidBodyMotionDisplacement;
    solidBodyMotionFunction     linearMotion;
velocity  (0 2 0);

}

}

其中INLET（火箭喷管）和ROCKET_WALL（火箭壁面）设置为以2m/s运动的边界，但存在以下两个问题

一.用算例运行到2s时提示以下信息

Courant Number mean: 9.51488e+118 max: 1.54754e+123
Time = 0.995

DICPCG: Solving for cellDisplacementx, Initial residual = 0, Final residual = 0, No Iterations 0
DICPCG: Solving for cellDisplacementy, Initial residual = 0.00667525, Final residual = 9.9838e-09, No Iterations 43
smoothSolver: Solving for Ux, Initial residual = 0.749586, Final residual = 9.93052e-07, No Iterations 65
smoothSolver: Solving for Uy, Initial residual = 0.78604, Final residual = 9.80661e-07, No Iterations 62
#0 Foam::error::printStack(Foam::Ostream&) at ??:?
#1 Foam::sigFpe::sigHandler(int) at ??:?
#2 ? in "/lib/x86_64-linux-gnu/libc.so.6"
#3 double Foam::sumProd<double>(Foam::UList<double> const&, Foam::UList<double> const&) at ??:?
#4 Foam::PCG::solve(Foam::Field<double>&, Foam::Field<double> const&, unsigned char) const at ??:?
#5 Foam::GAMGSolver::solveCoarsestLevel(Foam::Field<double>&, Foam::Field<double> const&) const at ??:?
#6 Foam::GAMGSolver::Vcycle(Foam::PtrListFoam::lduMatrix::smoother const&, Foam::Field<double>&, Foam::Field<double> const&, Foam::Field<double>&, Foam::Field<double>&, Foam::Field<double>&, Foam::Field<double>&, Foam::Field<double>&, Foam::PtrList<Foam::Field<double> >&, Foam::PtrList<Foam::Field<double> >&, unsigned char) const at ??:?
#7 Foam::GAMGSolver::solve(Foam::Field<double>&, Foam::Field<double> const&, unsigned char) const at ??:?
#8 Foam::fvMatrix<double>::solveSegregated(Foam::dictionary const&) at ??:?
#9 Foam::fvMatrix<double>::solve(Foam::dictionary const&) in "/home/dyfluid/OpenFOAM/OpenFOAM-7/platforms/linux64GccDPInt32Opt/bin/pimpleFoam"
#10 Foam::fvMatrix<double>::solve() in "/home/dyfluid/OpenFOAM/OpenFOAM-7/platforms/linux64GccDPInt32Opt/bin/pimpleFoam"
#11 ? in "/home/dyfluid/OpenFOAM/OpenFOAM-7/platforms/linux64GccDPInt32Opt/bin/pimpleFoam"
#12 __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6"
#13 ? in "/home/dyfluid/OpenFOAM/OpenFOAM-7/platforms/linux64GccDPInt32Opt/bin/pimpleFoam"
Floating point exception (core dumped)

查看了下网格发现网格变形为负了

请问这种情况用动网格怎么解决？

二我同时想在INLET（火箭喷管口）设置燃气射流环境，目前改了 0/p 0/U 0/T 下的边界条件，并在constant/thermophysicalProperties 设置了燃气参数，但不知道如何指定在INLET上不断喷射燃气。
请问该如何设置呢，十分感谢！

Lisiye

Foamers 大家好

现在我再做一个有关模拟火箭发射的瞬态流场仿真，再研究了不同的方法后有以下几种网格方法可以使用

1.动网格法
2.overset grid（嵌套网格法）
3.自适应网格

之前做过圆柱上下运动绕流的案例，所以打算先选用动网格法进行模拟，但遇到了问题。之前圆柱动网格算例 0文件中 pointDisplacement 圆柱边界设置条件如下：

    cylinder
    {
        type            oscillatingDisplacement;
	omega           0.5; //频率
	amplitude       (0 1 0); //幅度
	value           uniform (0 0 0);
    }

但我对于动网格 0文件中 pointDisplacement 的网格运动类型参数设置不太熟悉，不知道如何去设置一个火箭以一定的加速度向上运动。目前OF提示我的pointDisplacement边界类型如下

{angularOscillatingDisplacement
angularOscillatingVelocity
calculated
codedFixedValue
cyclic
cyclicACMI
cyclicAMI
cyclicRepeatAMI
cyclicSlip
empty
fixedNormalSlip
fixedValue
nonuniformTransformCyclic
oscillatingDisplacement
oscillatingVelocity
processor
processorCyclic
slip
solidBodyMotionDisplacement
surfaceDisplacement
surfaceSlipDisplacement
symmetry
symmetryPlane
timeVaryingMappedFixedValue
timeVaryingUniformFixedValue
uniformFixedValue
uniformInterpolatedDisplacement
value
waveDisplacement
wedge
zeroGradient}

由于我对这些边界类型不太熟，目前我选择的是codedFixedValue或solidBodyMotionDisplacement边界条件。请问Foamers们我该如何去设置一个火箭以一定的加速度向上运动呢？（我在GitHub上查找了关于动网格的算例但大多为overset旋转案例，没有我这种以一定加速度直线运动的 )

不胜感激，相关的资料讲解参数设置也可以！

如果有 overset grid 和自适应网格的FOAM友也欢迎提出一些解决的办法，感谢！（我在youtube上看到了有用CONVERGENCE自适应网格模拟火箭发射的案例，网址https://www.youtube.com/watch?v=JYGJbhRHAzU）

Lisiye

@xpqiu 在 oscillatingFixedValue边界条件的问题中说：

您好请问如何查找0文件夹中pointDisplacement 网格运动边界类型的参数设置说明

CFD中文网

Lisiye

帖子