CFD中文网

好

@stardust 用过，网格量上去的话，确实算不动。特别慢

好

@shrine 时间太长了，我有些找不到当时调通的文件啦。印象比较深刻的是，天河超算上不好调通，我是在ubuntu上的openfoam里调通的，你可以考虑先在ubuntu上Allwmake一下。

好

各位大佬，我参考文献 3rd High-Lift Workshop Summary Paper - OpenFOAM, STAR-CCM+ & LAVA simulations on Unstructured Grids链接https://arc.aiaa.org/doi/10.2514/6.2018-1253。用rhoPimpleFoam计算飞机，但是算了几步就发散了，希望各位大佬帮忙看看是fvSchemes和fvSolution的问题嘛，还是nut、nuTilda设置的问题。

目前用的是作者上传到会议网站的网格，按照文章里说的，采用的是SA模型，以及动量和湍流项采用的是二阶迎风。

Courant Number mean: 1.4881475e-08 max: 1.8237221
Time = 3.2e-09


* 
* PIMPLE: Iteration 1
diagonal:  Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
DILUPBiCGStab:  Solving for Ux, Initial residual = 2.421018e-11, Final residual = 2.421018e-11, No Iterations 0
DILUPBiCGStab:  Solving for Uy, Initial residual = 2.1538356e-10, Final residual = 2.1538356e-10, No Iterations 0
DILUPBiCGStab:  Solving for Uz, Initial residual = 2.0979924e-10, Final residual = 2.0979924e-10, No Iterations 0
DILUPBiCGStab:  Solving for e, Initial residual = 4.1760218e-08, Final residual = 2.1963671e-13, No Iterations 1
GAMG:  Solving for p, Initial residual = 7.5844072e-07, Final residual = 7.5844072e-07, No Iterations 0
GAMG:  Solving for p, Initial residual = 7.5844072e-07, Final residual = 7.5844072e-07, No Iterations 0
GAMG:  Solving for p, Initial residual = 7.5844072e-07, Final residual = 7.5844072e-07, No Iterations 0
GAMG:  Solving for p, Initial residual = 7.5844072e-07, Final residual = 7.5844072e-07, No Iterations 0
diagonal:  Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
time step continuity errors : sum local = 5.0250634e-11, global = 5.0210593e-11, cumulative = 3.7505015e-09
GAMG:  Solving for p, Initial residual = 7.5844071e-07, Final residual = 7.5844071e-07, No Iterations 0
GAMG:  Solving for p, Initial residual = 7.5844071e-07, Final residual = 7.5844071e-07, No Iterations 0
GAMG:  Solving for p, Initial residual = 7.5844071e-07, Final residual = 7.5844071e-07, No Iterations 0
GAMG:  Solving for p, Initial residual = 7.5844071e-07, Final residual = 7.5844071e-07, No Iterations 0
diagonal:  Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
time step continuity errors : sum local = 5.0250633e-11, global = 5.0210593e-11, cumulative = 3.8007121e-09
PIMPLE: Iteration 2
DILUPBiCGStab:  Solving for Ux, Initial residual = 1.3321419e-08, Final residual = 1.002751e-14, No Iterations 1
DILUPBiCGStab:  Solving for Uy, Initial residual = 6.2117311e-09, Final residual = 6.2117311e-09, No Iterations 0
DILUPBiCGStab:  Solving for Uz, Initial residual = 1.4225061e-10, Final residual = 1.4225061e-10, No Iterations 0
DILUPBiCGStab:  Solving for e, Initial residual = 1.5129023e-07, Final residual = 2.4037934e-13, No Iterations 1
GAMG:  Solving for p, Initial residual = 1.5142431e-06, Final residual = 1.5142431e-06, No Iterations 0
GAMG:  Solving for p, Initial residual = 1.5142431e-06, Final residual = 1.5142431e-06, No Iterations 0
GAMG:  Solving for p, Initial residual = 1.5142431e-06, Final residual = 1.5142431e-06, No Iterations 0
GAMG:  Solving for p, Initial residual = 1.5142431e-06, Final residual = 1.5142431e-06, No Iterations 0
diagonal:  Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
time step continuity errors : sum local = 1.0032647e-10, global = 1.0025895e-10, cumulative = 3.900971e-09
GAMG:  Solving for p, Initial residual = 1.5142431e-06, Final residual = 1.5142431e-06, No Iterations 0
GAMG:  Solving for p, Initial residual = 1.5142431e-06, Final residual = 1.5142431e-06, No Iterations 0
GAMG:  Solving for p, Initial residual = 1.5142431e-06, Final residual = 1.5142431e-06, No Iterations 0
GAMG:  Solving for p, Initial residual = 1.5142431e-06, Final residual = 1.5142431e-06, No Iterations 0
diagonal:  Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
time step continuity errors : sum local = 1.0032647e-10, global = 1.0025895e-10, cumulative = 4.00123e-09
PIMPLE: Iteration 3
DILUPBiCGStab:  Solving for Ux, Initial residual = 2.4133e-08, Final residual = 2.8587861e-14, No Iterations 1
DILUPBiCGStab:  Solving for Uy, Initial residual = 1.1420948e-08, Final residual = 1.9949995e-14, No Iterations 1
DILUPBiCGStab:  Solving for Uz, Initial residual = 1.4399722e-10, Final residual = 1.4399722e-10, No Iterations 0
DILUPBiCGStab:  Solving for e, Initial residual = 1.3956424e-07, Final residual = 2.3199397e-13, No Iterations 1
[73] #0  Foam::error::printStack(Foam::Ostream&)sh: addr2line: command not found
 addr2line failed
[73] #1  Foam::sigFpe::sigHandler(int)sh: addr2line: command not found
 addr2line failed
[73] #2  ?sh: addr2line: command not found
 addr2line failed
[73] #3  Foam::hePsiThermo<Foam::psiThermo, Foam::pureMixture<Foam::sutherlandTransport<Foam::species::thermo<Foam::hConstThermo<Foam::perfectGas<Foam::specie> >, Foam::sensibleInternalEnergy> > > >::calculate()sh: addr2line: command not found
 addr2line failed
[73] #4  Foam::hePsiThermo<Foam::psiThermo, Foam::pureMixture<Foam::sutherlandTransport<Foam::species::thermo<Foam::hConstThermo<Foam::perfectGas<Foam::specie> >, Foam::sensibleInternalEnergy> > > >::correct()sh: addr2line: command not found
 addr2line failed
[73] #5  ?sh: addr2line: command not found
 addr2line failed
[73] #6  __libc_start_mainsh: addr2line: command not found
 addr2line failed
[73] #7  ?sh: addr2line: command not found
 addr2line failed
Fatal error in PMPI_Waitall: Unknown error class, error stack:
PMPI_Waitall(405).....................: MPI_Waitall(count=40, req_array=0x3c75cd0, status_array=0x1) failed
MPIR_Waitall_impl(181)................: fail failed
MPIDU_Complete_posted_with_error(1710): Process failed
Fatal error in PMPI_Waitall: Unknown error class, error stack:
PMPI_Waitall(405).....................: MPI_Waitall(count=22, req_array=0x3966230, status_array=0x1) failed
MPIR_Waitall_impl(181)................: fail failed
MPIDU_Complete_posted_with_error(1710): Process failed
Fatal error in PMPI_Waitall: Unknown error class, error stack:
PMPI_Waitall(405).....................: MPI_Waitall(count=52, req_array=0x304c450, status_array=0x1) failed
MPIR_Waitall_impl(181)................: fail failed
MPIDU_Complete_posted_with_error(1710): Process failed
Fatal error in PMPI_Waitall: Unknown error class, error stack:
PMPI_Waitall(405).....................: MPI_Waitall(count=46, req_array=0x355b680, status_array=0x1) failed
MPIR_Waitall_impl(181)................: fail failed
MPIDU_Complete_posted_with_error(1710): Process failed
Fatal error in PMPI_Waitall: Unknown error class, error stack:
PMPI_Waitall(405).....................: MPI_Waitall(count=36, req_array=0x3dfabb0, status_array=0x1) failed
MPIR_Waitall_impl(181)................: fail failed
MPIDU_Complete_posted_with_error(1710): Process failed
Fatal error in PMPI_Waitall: Unknown error class, error stack:
PMPI_Waitall(405).....................: MPI_Waitall(count=40, req_array=0x349e540, status_array=0x1) failed
MPIR_Waitall_impl(181)................: fail failed
MPIDU_Complete_posted_with_error(1710): Process failed
Fatal error in PMPI_Waitall: Unknown error class, error stack:
PMPI_Waitall(405).....................: MPI_Waitall(count=26, req_array=0x3bd5ce0, status_array=0x1) failed
MPIR_Waitall_impl(181)................: fail failed
MPIDU_Complete_posted_with_error(1710): Process failed
Fatal error in PMPI_Waitall: Unknown error class, error stack:
PMPI_Waitall(405).....................: MPI_Waitall(count=12, req_array=0x3f18790, status_array=0x1) failed
MPIR_Waitall_impl(181)................: fail failed
MPIDU_Complete_posted_with_error(1710): Process failed
Fatal error in PMPI_Waitall: Unknown error class, error stack:
PMPI_Waitall(405).....................: MPI_Waitall(count=18, req_array=0x3f97b30, status_array=0x1) failed
MPIR_Waitall_impl(181)................: fail failed
MPIDU_Complete_posted_with_error(1710): Process failed
Fatal error in PMPI_Waitall: Unknown error class, error stack:
PMPI_Waitall(405).....................: MPI_Waitall(count=14, req_array=0x345ac20, status_array=0x1) failed
MPIR_Waitall_impl(181)................: fail failed
MPIDU_Complete_posted_with_error(1710): Process failed
Fatal error in PMPI_Waitall: Unknown error class, error stack:
PMPI_Waitall(405).....................: MPI_Waitall(count=32, req_array=0x3df8960, status_array=0x1) failed
MPIR_Waitall_impl(181)................: fail failed
MPIDU_Complete_posted_with_error(1710): Process failed
Fatal error in PMPI_Waitall: Unknown error class, error stack:
PMPI_Waitall(405).....................: MPI_Waitall(count=28, req_array=0x4316400, status_array=0x1) failed
MPIR_Waitall_impl(181)................: fail failed
MPIDU_Complete_posted_with_error(1710): Process failed
Fatal error in PMPI_Waitall: Unknown error class, error stack:
PMPI_Waitall(405).....................: MPI_Waitall(count=28, req_array=0x2fe0f30, status_array=0x1) failed
MPIR_Waitall_impl(181)................: fail failed
MPIDU_Complete_posted_with_error(1710): Process failed
Fatal error in PMPI_Waitall: Unknown error class, error stack:
PMPI_Waitall(405).....................: MPI_Waitall(count=50, req_array=0x3a36c60, status_array=0x1) failed
MPIR_Waitall_impl(181)................: fail failed
MPIDU_Complete_posted_with_error(1710): Process failed
Fatal error in PMPI_Waitall: Unknown error class, error stack:
PMPI_Waitall(405).....................: MPI_Waitall(count=24, req_array=0x3f6c5f0, status_array=0x1) failed
MPIR_Waitall_impl(181)................: fail failed
MPIDU_Complete_posted_with_error(1710): Process failed
Fatal error in PMPI_Waitall: Unknown error class, error stack:
PMPI_Waitall(405).....................: MPI_Waitall(count=24, req_array=0x520edf0, status_array=0x1) failed
MPIR_Waitall_impl(181)................: fail failed
MPIDU_Complete_posted_with_error(1710): Process failed
slurmd[cn10696]: *** STEP 17962804.0 KILLED AT 2021-12-21T17:17:15 WITH SIGNAL 9 ***
slurmd[cn10696]: *** STEP 17962804.0 KILLED AT 2021-12-21T17:17:15 WITH SIGNAL 9 ***
slurmd[cn10696]: *** STEP 17962804.0 KILLED AT 2021-12-21T17:17:15 WITH SIGNAL 9 ***
slurmd[cn10696]: *** STEP 17962804.0 KILLED AT 2021-12-21T17:17:15 WITH SIGNAL 9 ***
slurmd[cn10696]: *** STEP 17962804.0 KILLED AT 2021-12-21T17:17:15 WITH SIGNAL 9 ***
slurmd[cn10696]: *** STEP 17962804.0 KILLED AT 2021-12-21T17:17:15 WITH SIGNAL 9 ***
slurmd[cn10696]: *** STEP 17962804.0 KILLED AT 2021-12-21T17:17:15 WITH SIGNAL 9 ***
slurmd[cn10696]: *** STEP 17962804.0 KILLED AT 2021-12-21T17:17:15 WITH SIGNAL 9 ***
slurmd[cn10696]: *** STEP 17962804.0 KILLED AT 2021-12-21T17:17:15 WITH SIGNAL 9 ***
slurmd[cn10696]: *** STEP 17962804.0 KILLED AT 2021-12-21T17:17:15 WITH SIGNAL 9 ***
slurmd[cn10696]: *** STEP 17962804.0 KILLED AT 2021-12-21T17:17:15 WITH SIGNAL 9 ***
slurmd[cn10696]: *** STEP 17962804.0 KILLED AT 2021-12-21T17:17:15 WITH SIGNAL 9 ***
slurmd[cn10696]: *** STEP 17962804.0 KILLED AT 2021-12-21T17:17:15 WITH SIGNAL 9 ***
slurmd[cn10696]: *** STEP 17962804.0 KILLED AT 2021-12-21T17:17:15 WITH SIGNAL 9 ***
slurmd[cn10696]: *** STEP 17962804.0 KILLED AT 2021-12-21T17:17:16 WITH SIGNAL 9 ***
slurmd[cn10696]: *** STEP 17962804.0 KILLED AT 2021-12-21T17:17:16 WITH SIGNAL 9 ***
slurmd[cn10696]: *** STEP 17962804.0 KILLED AT 2021-12-21T17:17:16 WITH SIGNAL 9 ***

fvSchmes

// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

ddtSchemes
{
    default         Euler;
}

gradSchemes
{
    default         Gauss linear;

    limited         cellLimited Gauss linear 1;
    grad(U)         $limited;
    grad(k)         $limited;
    grad(omega)     $limited;
}

divSchemes
{
    default         none;

    div(phi,U)      Gauss linearUpwind grad(U);//Gauss linearUpwind limited;

    turbulence      Gauss linearUpwind default;//Gauss linearUpwind limited;
    energy          Gauss linearUpwind default;//Gauss linearUpwind limited;

    div(phi,k)      $turbulence;
    div(phi,omega)  $turbulence;

    div(phi,e)      $energy;
    div(phi,K)      $energy;
    div(phi,Ekp)    $energy;
    div(phi,nuTilda)  $energy;

    div(phiv,p)     Gauss upwind;
    div((phi|interpolate(rho)),p) Gauss upwind;

    div(((rho*nuEff)*dev2(T(grad(U)))))    Gauss linear;
}

laplacianSchemes
{
    default         Gauss linear limited 0.5;//Gauss linear corrected;
}

interpolationSchemes
{
    default         linear;
}

snGradSchemes
{
    default         limited 0.5;//corrected;
}

wallDist
{
    method meshWave;
}

fvSolution


// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

solvers
{
    p
    {
        solver          GAMG;
        smoother        GaussSeidel;
        tolerance       1e-5;
        relTol          0.01;
		nPreSweeps      0;
		nPostSweeps     2;
		cacheAgglomeration  on;
		agglomerator      faceAreaPair;
		ncellsInCoarsetLevel      1000;
		mergeLevels          1;
    }

    pFinal
    {
        $p;
        relTol          0;
    }

    "(rho|U|k|omega|e)"
    {
        solver          PBiCGStab;
        preconditioner  DILU;
        tolerance       1e-8;
        relTol          0.0;
    }

    "(rho|U|k|omega|e|nuTilda)Final"
    {
        $U;
        relTol          0;
    }
}

SIMPLE
{
    residualControl
    {
        p               1e-5;
        U               1e-5;
        "(k|omega|e)"   1e-5;
    }

    nNonOrthogonalCorrectors 3;
    pMinFactor      0.1;
    pMaxFactor      2;
}

PIMPLE
{
    nCorrectors              2;
    nNonOrthogonalCorrectors 3;
    nOuterCorrectors         3;
    pMinFactor      0.1;
    pMaxFactor      2;
}

relaxationFactors
{
	p       0.3;
	U       0.7;
	k       0.7;
	omega   0.7;
 
}

// ************************************************************************* //

nut


// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

dimensions      [0 2 -1 0 0 0 0];

internalField   uniform 0;

boundaryField
{
    freestream
    {
        type            calculated;
        value           uniform 0;
    }

    wall
    {
        type            nutkWallFunction;
        type            nutkWallFunction;
        value           uniform 0;
    }

    #includeEtc "caseDicts/setConstraintTypes"
}

// ************************************************************************* //

nuTilda


// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

dimensions      [0 2 -1 0 0 0 0];

internalField   uniform 0;

boundaryField
{
    freestream
    {
        type            inletOutlet;
        inletValue      $internalField;
        value           $internalField;
    }

    wall
    {
        type            fixedValue;
        value           uniform 0;
    }

    symmetry
    {
        type            symmetry;
    }

    #includeEtc "caseDicts/setConstraintTypes"
}

// ************************************************************************* //

好

请教一下，我在icem保存用acsii格式保存了msh文件，用fluent3DMeshToFoam转化网格时出现如下报错，不知道是什么原因呢，操作是在天河二号上进行的。

patch 15 from Fluent indices: 59793665 to: 59902377 type: wall
#0  Foam::error::printStack(Foam::Ostream&)sh: addr2line: command not found
 addr2line failed
#1  Foam::sigSegv::sigHandler(int)sh: addr2line: command not found
 addr2line failed
#2  ?sh: addr2line: command not found
 addr2line failed
#3  Foam::polyTopoChange::renumberCompact(Foam::List<int> const&, Foam::List<int>&)sh: addr2line: command not found
 addr2line failed
#4  Foam::polyTopoChange::compact(bool, bool, int&, Foam::List<int>&, Foam::List<int>&)sh: addr2line: command not found
 addr2line failed
#5  Foam::polyTopoChange::compactAndReorder(Foam::polyMesh const&, bool, bool, bool, int&, Foam::Field<Foam::Vector<double> >&, Foam::List<int>&, Foam::List<int>&, Foam::List<Foam::objectMap>&, Foam::List<Foam::objectMap>&, Foam::List<Foam::objectMap>&, Foam::List<Foam::objectMap>&, Foam::List<Foam::objectMap>&, Foam::List<Foam::objectMap>&, Foam::List<Foam::objectMap>&, Foam::List<Foam::objectMap>&, Foam::List<Foam::Map<int> >&, Foam::List<int>&, Foam::List<int>&, Foam::List<Foam::Map<int> >&)sh: addr2line: command not found
 addr2line failed
#6  Foam::polyTopoChange::changeMesh(Foam::polyMesh&, bool, bool, bool, bool)sh: addr2line: command not found
 addr2line failed
#7  ?sh: addr2line: command not found
 addr2line failed
#8  __libc_start_mainsh: addr2line: command not found
 addr2line failed
#9  ?sh: addr2line: command not found
 addr2line failed

好

@number44 我这个p和U的tolerance分别设置的是1e-8和1e-9，感觉挺小的，还要调的更小嘛

好

各位大佬，想问一下我计算一个飞机模型的外流，设置的dt为2e-8，但是不知道为什么从最开始就是No Iterations 0，这是没在计算嘛？刚开始设置的dt为5e-8，计算一段时间就发散了。

Courant Number mean: 7.5311001e-06 max: 0.56432736
Time = 1.2e-07

diagonal:  Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
DILUPBiCGStab:  Solving for Ux, Initial residual = 6.9178663e-11, Final residual = 6.9178663e-11, No Iterations 0
DILUPBiCGStab:  Solving for Uy, Initial residual = 2.5448919e-10, Final residual = 2.5448919e-10, No Iterations 0
DILUPBiCGStab:  Solving for Uz, Initial residual = 4.2785621e-10, Final residual = 4.2785621e-10, No Iterations 0
DILUPBiCGStab:  Solving for e, Initial residual = 9.7719937e-09, Final residual = 1.0109782e-12, No Iterations 1
GAMG:  Solving for p, Initial residual = 7.8947472e-09, Final residual = 7.8947472e-09, No Iterations 0
GAMG:  Solving for p, Initial residual = 7.8947472e-09, Final residual = 7.8947472e-09, No Iterations 0
diagonal:  Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
time step continuity errors : sum local = 1.0995375e-11, global = -1.6176582e-13, cumulative = -8.4645739e-13
GAMG:  Solving for p, Initial residual = 7.8947534e-09, Final residual = 7.8947534e-09, No Iterations 0
GAMG:  Solving for p, Initial residual = 7.8947534e-09, Final residual = 7.8947534e-09, No Iterations 0
diagonal:  Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
time step continuity errors : sum local = 1.0995384e-11, global = -1.6176582e-13, cumulative = -1.0082232e-12
DILUPBiCGStab:  Solving for nuTilda, Initial residual = 9.7105217e-10, Final residual = 9.7105217e-10, No Iterations 0

好

@五好青年用tecplot，先slice后extract，就可以了。具体参考这个https://www.bilibili.com/video/BV1Mr4y1w71c/

好

@五好青年 hello，你解决了这个问题吗

好

@static 你好，我最近也在做三维飞机的相关数值模拟，可否给个联系方式，请教一些问题。谢谢，这是我的qq：3310792254

好

希望各位大佬指点一下，就是我想用双时间步，应该怎么改呢。
如下图代码，Info出来的deltaT还是controldict里的deltaT

runTime.deltaT() =0.00001 * demensiont;
    Info<< "deltaTime = " << runTime.deltaT() << nl << endl;
fvScalarMatrix FEqn
(
     	  fvm::ddt(H)
	+ fvm::div(phi,H)
	- fvm::laplacian(DH,H)
);

好

算出来的等值面结果如图

好

计算一个对称的平板绕流，速度场和涡量大小的云图都对称，利用如下公式计算另一个变量的时候，计算出来的等值面不对称（应该是轻微不对称，但是我现在十分不对称），请问各位大佬，可能有什么原因呢。

	surfaceScalarField phiw
	(
		"phiw",
                fvc::flux(omega) * LL
	);
	

	for (int ntau = 1; ntau <= 5; ntau++)
	{
		H.oldTime() = H;
		fvScalarMatrix HEqn
		(
	             	  fvm::ddt(H)
			+ fvm::div(phiw,H)
			- fvm::laplacian(DH,H)
		);
		HEqn.solve();
	}

好

@random_ran 您好！我最近也在做三维圆柱扰流的相关工作，我先开始做的是长方形流域的二维圆柱绕流，感觉还能差不多吻合升阻力系数，但是三维数值模拟很慢，看到您的帖子，开始转移到圆形流域，发现升阻力系数会随着流域大小的变化而变化，三维流域的话，就算的很慢很慢，差不多要12个小时才要算1个时间单位。如果方便的话，可否把你的算例或者网格让我参考一下~x谢谢！！！这是我的邮箱~3310792254@qq.com~qq.com~祝好！

好

@波流力可否推荐一下上交那篇植入其他涡量识别方法的文章。我目前也在搞相关的工作，想借鉴一下。谢谢！

好

好

在超算上跑$64^{3}$的网格，可以并行，但是跑$256^{3}$的网格的时候不能并行，出现如下报错，请问是什么原因呢，各位大佬。

Create mesh for time = 0

Create mesh for time = 0

Create mesh for time = 0

Create mesh for time = 0

slurmd[cn10253]: *** STEP 15265985.0 KILLED AT 2020-07-29T00:45:34 WITH SIGNAL 9 ***
slurmd[cn10253]: *** STEP 15265985.0 KILLED AT 2020-07-29T00:45:34 WITH SIGNAL 9 ***

好

@东岳是的，东岳老师，的确是挺困难的，我水平现在还不够。我这两天先把那个github上的WENO程序调通了，打算先用这个weno程序跑一下。抱歉回复那么晚~~谢谢东岳老师。

好

谢谢东岳老师！我还想问一下，那我现在想写一个weno格式的话，我是不是要在src/finiteVolume/interpolation/surfaceInterpolation/limitedSchemes目录下新建一个weno的文件，然后在文件里写C文件和H文件。写完以后在src/finiteVolume/Make/files里面添加编译文件，最后在finiteVolume里编译一下，是不是就可以像openfoam自带的格式一样使用了呢。主要之前完全没搞过，不太确定这个步骤对不对。谢谢东岳老师！！！

好

大家好，我最近要把两个方程加入到icoFoam中去，现在需要在openfoam里加入新的差分格式方便和dns程序做对比。但是我目前只了解如何在求解器里加新方程的步骤，不了解加入新格式的步骤。请问各位大佬能稍微讲一下在哪个目录下添加新格式吗，或者有什么参考文献吗？谢谢！！！

好

hello，不知道你后来有设置这个双时间步吗？可能给看一下你的设置吗？谢谢！

CFD中文网

好好科研

帖子