李老师,我尝试了rhoPimpleFoam没有问题，可以正常运行。但是rhoCentralFoam直接报方向错误。
rhoCentralFoam在不加动网格情况下，可以运行，但最终因为浮点数溢出停止。
Y
yingqing
@yingqing
帖子
-
报错求助“Operator + is undefined for unoriented and oriented types” -
报错求助“Operator + is undefined for unoriented and oriented types”在rhoCentralFoam求解中，使用动网格求解遇到错误。
动网格设置文件/*--------------------------------*- C++ -*----------------------------------*\ | ========= | | | \\ / F ield | OpenFOAM: The Open Source CFD Toolbox | | \\ / O peration | Version: v2206 | | \\ / A nd | Website: www.openfoam.com | | \\/ M anipulation | | \*---------------------------------------------------------------------------*/ FoamFile { version 2.0; format ascii; class dictionary; object dynamicMeshDict; } // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // dynamicFvMesh dynamicMotionSolverFvMesh; motionSolver solidBody; cellZone PASSAGE; // 旋转区域的cellZone名称 solidBodyMotionFunction rotatingMotion; origin (0 0 0); // 旋转中心坐标 axis (0 0 1); // 旋转轴方向 (0 0 1) 表示绕Z轴旋转 omega -1700; // 旋转角速度 [rad/s] // ************************************************************************* //详细的报错结果
/*---------------------------------------------------------------------------*\ | ========= | | | \\ / F ield | OpenFOAM: The Open Source CFD Toolbox | | \\ / O peration | Version: 2206 | | \\ / A nd | Website: www.openfoam.com | | \\/ M anipulation | | \*---------------------------------------------------------------------------*/ Build : _76d719d1e6-20220624 OPENFOAM=2206 version=v2206 Arch : "LSB;label=32;scalar=64" Exec : rhoCentralFoam Date : Aug 27 2025 Time : 20:53:55 Host : sdu-R940xa PID : 3759768 I/O : uncollated Case : /data0/XiaoYongqing/compressSingleBlade/Test10 nProcs : 1 trapFpe: Floating point exception trapping enabled (FOAM_SIGFPE). fileModificationChecking : Monitoring run-time modified files using timeStampMaster (fileModificationSkew 5, maxFileModificationPolls 20) allowSystemOperations : Allowing user-supplied system call operations // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // Create time Create mesh for time = 0 Selecting dynamicFvMesh dynamicMotionSolverFvMesh Selecting motion solver: solidBody Applying motion to cellZone: PASSAGE Selecting solid-body motion function rotatingMotion Reading thermophysical properties Selecting thermodynamics package { type hePsiThermo; mixture pureMixture; transport const; thermo hConst; equationOfState perfectGas; specie specie; energy sensibleInternalEnergy; } AMI: Creating addressing and weights between 4515 source faces and 4515 target faces AMI: Patch source sum(weights) min:1 max:1 average:1 AMI: Patch target sum(weights) min:1 max:1 average:1 AMI: Creating addressing and weights between 5005 source faces and 5005 target faces AMI: Patch source sum(weights) min:1 max:1 average:1 AMI: Patch target sum(weights) min:1 max:1 average:1 AMI: Creating addressing and weights between 1645 source faces and 1645 target faces AMI: Patch source sum(weights) min:1 max:1 average:1 AMI: Patch target sum(weights) min:1 max:1 average:1 Reading field U Creating turbulence model Selecting turbulence model type RAS Selecting RAS turbulence model kOmegaSST Selecting patchDistMethod meshWave RAS { RASModel kOmegaSST; turbulence on; printCoeffs on; alphaK1 0.85; alphaK2 1; alphaOmega1 0.5; alphaOmega2 0.856; gamma1 0.555556; gamma2 0.44; beta1 0.075; beta2 0.0828; betaStar 0.09; a1 0.31; b1 1; c1 10; F3 false; decayControl false; kInf 0; omegaInf 0; } Creating finite volume options from "system/fvOptions" Selecting finite volume options type limitTemperature Source: limitT - selecting all cells - selected 598360 cell(s) with volume 5.59104e-05 Starting time loop surfaceFieldValue outletMassFlow: operation = sum AMI: Creating addressing and weights between 5005 source faces and 5005 target faces AMI: Patch source sum(weights) min:0.999874 max:1 average:0.999999 AMI: Patch target sum(weights) min:0.999898 max:1 average:0.999999 AMI: Creating addressing and weights between 4515 source faces and 4515 target faces AMI: Patch source sum(weights) min:1 max:1 average:1 AMI: Patch target sum(weights) min:1 max:1 average:1 AMI: Creating addressing and weights between 1645 source faces and 1645 target faces AMI: Patch source sum(weights) min:1 max:1 average:1 AMI: Patch target sum(weights) min:1 max:1 average:1 Mean and max Courant Numbers = 0.00114965 0.0538708 Time = 1e-09 diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0 diagonal: Solving for rhoUx, Initial residual = 0, Final residual = 0, No Iterations 0 diagonal: Solving for rhoUy, Initial residual = 0, Final residual = 0, No Iterations 0 diagonal: Solving for rhoUz, Initial residual = 0, Final residual = 0, No Iterations 0 smoothSolver: Solving for Ux, Initial residual = 2.72958e-05, Final residual = 6.36242e-14, No Iterations 1 smoothSolver: Solving for Uy, Initial residual = 1.33319e-05, Final residual = 3.54757e-14, No Iterations 1 smoothSolver: Solving for Uz, Initial residual = 0.00560242, Final residual = 1.89806e-11, No Iterations 1 diagonal: Solving for rhoE, Initial residual = 0, Final residual = 0, No Iterations 0 smoothSolver: Solving for e, Initial residual = 0.00593518, Final residual = 3.01268e-11, No Iterations 1 --> FOAM FATAL ERROR: (openfoam-2206) Operator + is undefined for unoriented and oriented types From Foam::orientedType Foam::operator+(const Foam::orientedType&, const Foam::orientedType&) in file orientedType/orientedType.C at line 488. FOAM aborting #0 Foam::error::printStack(Foam::Ostream&) at ??:? #1 Foam::error::simpleExit(int, bool) at ??:? #2 Foam::error::exiting(int, bool) at ??:? #3 Foam::operator+(Foam::orientedType const&, Foam::orientedType const&) at ??:? #4 void Foam::add<double, double, Foam::fvsPatchField, Foam::surfaceMesh>(Foam::GeometricField<Foam::typeOfSum<double, double>::type, Foam::fvsPatchField, Foam::surfaceMesh>&, Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> const&, Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> const&) in ~/OpenFOAM/OpenFOAM-v2206/platforms/linux64GccDPInt32Opt/bin/rhoCentralFoam #5 Foam::tmp<Foam::GeometricField<Foam::typeOfSum<double, double>::type, Foam::fvsPatchField, Foam::surfaceMesh> > Foam::operator+<double, double, Foam::fvsPatchField, Foam::surfaceMesh>(Foam::tmp<Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> > const&, Foam::tmp<Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> > const&) in ~/OpenFOAM/OpenFOAM-v2206/platforms/linux64GccDPInt32Opt/bin/rhoCentralFoam #6 Foam::fvc::absolute(Foam::tmp<Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> > const&, Foam::GeometricField<Foam::Vector<double>, Foam::fvPatchField, Foam::volMesh> const&) at ??:? #7 Foam::kOmegaSSTBase<Foam::eddyViscosity<Foam::RASModel<Foam::EddyDiffusivity<Foam::ThermalDiffusivity<Foam::CompressibleTurbulenceModel<Foam::fluidThermo> > > > > >::correct() at ??:? #8 ? in ~/OpenFOAM/OpenFOAM-v2206/platforms/linux64GccDPInt32Opt/bin/rhoCentralFoam #9 __libc_start_main in /lib/x86_64-linux-gnu/libc.so.6 #10 ? in ~/OpenFOAM/OpenFOAM-v2206/platforms/linux64GccDPInt32Opt/bin/rhoCentralFoam -
离心式压缩机周期性边界问题问题解决！在使用旋转周期性交界面的时候把 rotationAngle 45去掉
-
离心式压缩机周期性边界问题各位大佬，我在做压缩机并行运算的时候遇到周期性边界条件的问题。
这是我压缩机单流道模型
对于交界面的处理如下：/*--------------------------------*- C++ -*----------------------------------*\ | ========= | | | \\ / F ield | OpenFOAM: The Open Source CFD Toolbox | | \\ / O peration | Version: 2206 | | \\ / A nd | Website: www.openfoam.com | | \\/ M anipulation | | \*---------------------------------------------------------------------------*/ FoamFile { version 2.0; format ascii; arch "LSB;label=32;scalar=64"; class polyBoundaryMesh; location "constant/polyMesh"; object boundary; } // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // 14 ( PASSAGE-HUB { type wall; inGroups 1(wall); nFaces 10375; startFace 1756179; } PASSAGE-SHROUD { type wall; inGroups 1(wall); nFaces 10375; startFace 1766554; } PASSAGE-INFLOW { type patch; nFaces 1330; startFace 1776929; } PASSAGE-OUTFLOW { type cyclicAMI; nFaces 5005; startFace 1778259; matchTolerance 10; transform noOrdering; neighbourPatch OUTBLOCK-INFLOW; AMIMethod faceAreaWeightAMI; restartUncoveredSourceFace 1; } MAIN-BLADE { type wall; inGroups 1(wall); nFaces 9065; startFace 1783264; } SPLITTER-BLADE { type wall; inGroups 1(wall); nFaces 5880; startFace 1792329; } PASSAGE-PER1 { type cyclicAMI; nFaces 4515; startFace 1798209; matchTolerance 10; transform unknown; rotationAxis (0 0 1); rotationAngle 45; rotationCentre (0 0 0); neighbourPatch PASSAGE-PER2; AMIMethod partialFaceAreaWeightAMI; restartUncoveredSourceFace 1; } PASSAGE-PER2 { type cyclicAMI; nFaces 4515; startFace 1802724; matchTolerance 10; transform unknown; rotationAxis (0 0 1); rotationAngle 45; rotationCentre (0 0 0); neighbourPatch PASSAGE-PER1; AMIMethod partialFaceAreaWeightAMI; restartUncoveredSourceFace 1; } OUTBLOCK-HUB { type wall; inGroups 1(wall); nFaces 6721; startFace 1807239; } OUTBLOCK-SHROUD { type wall; inGroups 1(wall); nFaces 6721; startFace 1813960; } OUTBLOCK-INFLOW { type cyclicAMI; nFaces 5005; startFace 1820681; matchTolerance 10; transform noOrdering; neighbourPatch PASSAGE-OUTFLOW; AMIMethod faceAreaWeightAMI; restartUncoveredSourceFace 1; } OUTBLOCK-OUTFLOW { type patch; nFaces 5005; startFace 1825686; } OUTBLOCK-PER1 { type cyclicAMI; nFaces 1645; startFace 1830691; matchTolerance 10; transform unknown; rotationAxis (0 0 1); rotationAngle 45; rotationCentre (0 0 0); neighbourPatch OUTBLOCK-PER2; AMIMethod faceAreaWeightAMI; restartUncoveredSourceFace 1; } OUTBLOCK-PER2 { type cyclicAMI; nFaces 1645; startFace 1832336; matchTolerance 10; transform unknown; rotationAxis (0 0 1); rotationAngle 45; rotationCentre (0 0 0); neighbourPatch OUTBLOCK-PER1; AMIMethod faceAreaWeightAMI; restartUncoveredSourceFace 1; } ) // ************************************************************************* //但是在进行并行运行decomposePar命令时，给出如下错误报告
/*------------------------------------------- // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // Create time Decomposing mesh Create mesh Using [v2012] 'partialFaceAreaWeightAMI' instead of 'faceAreaWeightAMI' in selection table: AMIInterpolation This lookup is 18 months old. Using [v2012] 'partialFaceAreaWeightAMI' instead of 'faceAreaWeightAMI' in selection table: AMIInterpolation This lookup is 18 months old. Calculating distribution of cells Decomposition method scotch [32] (region region0) Finished decomposition in 3.32 s Calculating original mesh data Distributing cells to processors Distributing faces to processors Distributing points to processors Constructing processor meshes Processor 0 Number of cells = 18879 Number of points = 21302 Number of faces shared with processor 7 = 772 Number of faces shared with processor 8 = 140 Number of faces shared with processor 10 = 248 Number of faces shared with processor 16 = 757 Number of faces shared with processor 24 = 772 Number of faces shared with processor 26 = 84 Number of processor patches = 6 Number of processor faces = 2773 Number of boundary faces = 1897 Processor 1 Number of cells = 18685 Number of points = 21243 Number of faces shared with processor 3 = 134 Number of faces shared with processor 9 = 391 Number of faces shared with processor 11 = 175 Number of faces shared with processor 13 = 35 Number of faces shared with processor 15 = 16 Number of faces shared with processor 17 = 1002 Number of faces shared with processor 19 = 429 Number of faces shared with processor 21 = 89 Number of faces shared with processor 25 = 10 Number of faces shared with processor 27 = 34 Number of faces shared with processor 29 = 841 Number of processor patches = 11 Number of processor faces = 3156 Number of boundary faces = 1766 Processor 2 Number of cells = 18514 Number of points = 20845 Number of faces shared with processor 10 = 240 Number of faces shared with processor 18 = 912 Number of faces shared with processor 22 = 640 Number of faces shared with processor 26 = 716 Number of processor patches = 4 Number of processor faces = 2508 Number of boundary faces = 1980 Processor 3 Number of cells = 18648 Number of points = 21261 Number of faces shared with processor 1 = 134 Number of faces shared with processor 15 = 388 Number of faces shared with processor 19 = 734 Number of faces shared with processor 27 = 1075 Number of faces shared with processor 29 = 474 Number of processor patches = 5 Number of processor faces = 2805 Number of boundary faces = 2223 Processor 4 Number of cells = 18884 Number of points = 21280 Number of faces shared with processor 20 = 945 Number of faces shared with processor 30 = 897 Number of processor patches = 2 Number of processor faces = 1842 Number of boundary faces = 2778 Processor 5 Number of cells = 18512 Number of points = 21003 Number of faces shared with processor 13 = 1004 Number of faces shared with processor 17 = 115 Number of faces shared with processor 21 = 657 Number of faces shared with processor 25 = 919 Number of processor patches = 4 Number of processor faces = 2695 Number of boundary faces = 2097 Processor 6 Number of cells = 18884 Number of points = 21264 Number of faces shared with processor 14 = 982 Number of faces shared with processor 22 = 912 Number of processor patches = 2 Number of processor faces = 1894 Number of boundary faces = 2692 Processor 7 Number of cells = 18884 Number of points = 21214 Number of faces shared with processor 0 = 772 Number of faces shared with processor 10 = 386 Number of faces shared with processor 23 = 819 Number of processor patches = 3 Number of processor faces = 1977 Number of boundary faces = 2511 Processor 8 Number of cells = 18523 Number of points = 21152 Number of faces shared with processor 0 = 140 Number of faces shared with processor 12 = 940 Number of faces shared with processor 16 = 834 Number of faces shared with processor 24 = 1180 Number of faces shared with processor 28 = 55 Number of processor patches = 5 Number of processor faces = 3149 Number of boundary faces = 1903 Processor 9 Number of cells = 18802 Number of points = 21444 Number of faces shared with processor 1 = 391 Number of faces shared with processor 13 = 790 Number of faces shared with processor 17 = 136 Number of faces shared with processor 25 = 1159 Number of faces shared with processor 29 = 87 Number of processor patches = 5 Number of processor faces = 2563 Number of boundary faces = 2529 Processor 10 Number of cells = 18878 Number of points = 21365 Number of faces shared with processor 0 = 248 Number of faces shared with processor 2 = 240 Number of faces shared with processor 7 = 386 Number of faces shared with processor 18 = 785 Number of faces shared with processor 24 = 93 Number of faces shared with processor 26 = 992 Number of processor patches = 6 Number of processor faces = 2744 Number of boundary faces = 2046 Processor 11 Number of cells = 18512 Number of points = 20961 Number of faces shared with processor 1 = 175 Number of faces shared with processor 15 = 296 Number of faces shared with processor 17 = 322 Number of faces shared with processor 19 = 904 Number of faces shared with processor 27 = 714 Number of processor patches = 5 Number of processor faces = 2411 Number of boundary faces = 2303 Processor 12 Number of cells = 18520 Number of points = 21507 Number of faces shared with processor 8 = 940 Number of faces shared with processor 16 = 217 Number of faces shared with processor 20 = 271 Number of faces shared with processor 28 = 557 Number of processor patches = 4 Number of processor faces = 1985 Number of boundary faces = 3767 Processor 13 Number of cells = 18690 Number of points = 21228 Number of faces shared with processor 1 = 35 Number of faces shared with processor 5 = 1004 Number of faces shared with processor 9 = 790 Number of faces shared with processor 17 = 22 Number of faces shared with processor 21 = 181 Number of faces shared with processor 25 = 194 Number of faces shared with processor 29 = 498 Number of processor patches = 7 Number of processor faces = 2724 Number of boundary faces = 2158 Processor 14 Number of cells = 18884 Number of points = 21287 Number of faces shared with processor 6 = 982 Number of faces shared with processor 30 = 915 Number of processor patches = 2 Number of processor faces = 1897 Number of boundary faces = 2731 Processor 15 Number of cells = 18742 Number of points = 21822 Number of faces shared with processor 1 = 16 Number of faces shared with processor 3 = 388 Number of faces shared with processor 11 = 296 Number of faces shared with processor 17 = 124 Number of faces shared with processor 21 = 463 Number of faces shared with processor 23 = 408 Number of faces shared with processor 27 = 728 Number of faces shared with processor 29 = 173 Number of faces shared with processor 31 = 759 Number of processor patches = 9 Number of processor faces = 3355 Number of boundary faces = 2521 Processor 16 Number of cells = 18820 Number of points = 21189 Number of faces shared with processor 0 = 757 Number of faces shared with processor 8 = 834 Number of faces shared with processor 12 = 217 Number of faces shared with processor 24 = 117 Number of faces shared with processor 28 = 666 Number of processor patches = 5 Number of processor faces = 2591 Number of boundary faces = 1969 Processor 17 Number of cells = 18514 Number of points = 21000 Number of faces shared with processor 1 = 1002 Number of faces shared with processor 5 = 115 Number of faces shared with processor 9 = 136 Number of faces shared with processor 11 = 322 Number of faces shared with processor 13 = 22 Number of faces shared with processor 15 = 124 Number of faces shared with processor 21 = 719 Number of faces shared with processor 25 = 391 Number of processor patches = 8 Number of processor faces = 2831 Number of boundary faces = 1955 Processor 18 Number of cells = 18779 Number of points = 21403 Number of faces shared with processor 2 = 912 Number of faces shared with processor 10 = 785 Number of faces shared with processor 22 = 296 Number of faces shared with processor 26 = 122 Number of processor patches = 4 Number of processor faces = 2115 Number of boundary faces = 2939 Processor 19 Number of cells = 18512 Number of points = 20978 Number of faces shared with processor 1 = 429 Number of faces shared with processor 3 = 734 Number of faces shared with processor 11 = 904 Number of faces shared with processor 27 = 352 Number of processor patches = 4 Number of processor faces = 2419 Number of boundary faces = 2325 Processor 20 Number of cells = 18884 Number of points = 21151 Number of faces shared with processor 4 = 945 Number of faces shared with processor 12 = 271 Number of processor patches = 2 Number of processor faces = 1216 Number of boundary faces = 3150 Processor 21 Number of cells = 18512 Number of points = 21057 Number of faces shared with processor 1 = 89 Number of faces shared with processor 5 = 657 Number of faces shared with processor 13 = 181 Number of faces shared with processor 15 = 463 Number of faces shared with processor 17 = 719 Number of faces shared with processor 29 = 942 Number of processor patches = 6 Number of processor faces = 3051 Number of boundary faces = 1843 Processor 22 Number of cells = 18884 Number of points = 21297 Number of faces shared with processor 2 = 640 Number of faces shared with processor 6 = 912 Number of faces shared with processor 18 = 296 Number of processor patches = 3 Number of processor faces = 1848 Number of boundary faces = 2802 Processor 23 Number of cells = 18884 Number of points = 21260 Number of faces shared with processor 7 = 819 Number of faces shared with processor 15 = 408 Number of faces shared with processor 31 = 328 Number of processor patches = 3 Number of processor faces = 1555 Number of boundary faces = 3019 Processor 24 Number of cells = 18540 Number of points = 21032 Number of faces shared with processor 0 = 772 Number of faces shared with processor 8 = 1180 Number of faces shared with processor 10 = 93 Number of faces shared with processor 16 = 117 Number of faces shared with processor 26 = 765 Number of processor patches = 5 Number of processor faces = 2927 Number of boundary faces = 1869 Processor 25 Number of cells = 18720 Number of points = 21266 Number of faces shared with processor 1 = 10 Number of faces shared with processor 5 = 919 Number of faces shared with processor 9 = 1159 Number of faces shared with processor 13 = 194 Number of faces shared with processor 17 = 391 Number of processor patches = 5 Number of processor faces = 2673 Number of boundary faces = 2219 Processor 26 Number of cells = 18527 Number of points = 20937 Number of faces shared with processor 0 = 84 Number of faces shared with processor 2 = 716 Number of faces shared with processor 10 = 992 Number of faces shared with processor 18 = 122 Number of faces shared with processor 24 = 765 Number of processor patches = 5 Number of processor faces = 2679 Number of boundary faces = 1959 Processor 27 Number of cells = 18537 Number of points = 21158 Number of faces shared with processor 1 = 34 Number of faces shared with processor 3 = 1075 Number of faces shared with processor 11 = 714 Number of faces shared with processor 15 = 728 Number of faces shared with processor 19 = 352 Number of faces shared with processor 29 = 22 Number of processor patches = 6 Number of processor faces = 2925 Number of boundary faces = 2101 Processor 28 Number of cells = 18526 Number of points = 20963 Number of faces shared with processor 8 = 55 Number of faces shared with processor 12 = 557 Number of faces shared with processor 16 = 666 Number of faces shared with processor 31 = 234 Number of processor patches = 4 Number of processor faces = 1512 Number of boundary faces = 3178 Processor 29 Number of cells = 18512 Number of points = 21051 Number of faces shared with processor 1 = 841 Number of faces shared with processor 3 = 474 Number of faces shared with processor 9 = 87 Number of faces shared with processor 13 = 498 Number of faces shared with processor 15 = 173 Number of faces shared with processor 21 = 942 Number of faces shared with processor 27 = 22 Number of processor patches = 7 Number of processor faces = 3037 Number of boundary faces = 1841 Processor 30 Number of cells = 18884 Number of points = 21229 Number of faces shared with processor 4 = 897 Number of faces shared with processor 14 = 915 Number of processor patches = 2 Number of processor faces = 1812 Number of boundary faces = 2706 Processor 31 Number of cells = 18884 Number of points = 21660 Number of faces shared with processor 15 = 759 Number of faces shared with processor 23 = 328 Number of faces shared with processor 28 = 234 Number of processor patches = 3 Number of processor faces = 1321 Number of boundary faces = 4025 Number of processor faces = 38495 Max number of cells = 18884 (0.990708% above average 18698.8) Max number of processor patches = 11 (131.579% above average 4.75) Max number of faces between processors = 3355 (39.4467% above average 2405.94) Time = 0 AMI: Creating addressing and weights between 4515 source faces and 4515 target faces --> FOAM Warning : From void Foam::advancingFrontAMI::checkPatches() const in file AMIInterpolation/AMIInterpolation/advancingFrontAMI/advancingFrontAMI.C at line 71 Source and target patch bounding boxes are not similar source box span : (0.0420091 0.0511361 0.0349999) target box span : (0.0186894 0.064626 0.0349999) source box : (0.0208063 -0.00369368 0) (0.0628154 0.0474424 0.0349999) target box : (0.0099481 0.0127 0) (0.0286375 0.077326 0.0349999) inflated target box : (0.0061564 0.00890831 -0.0037917) (0.0324292 0.0811177 0.0387916) --> FOAM FATAL ERROR: (openfoam-2206) Unable to set target face for source face 4514 From virtual bool Foam::faceAreaWeightAMI::setNextFaces(Foam::label&, Foam::label&, Foam::label&, const Foam::bitSet&, Foam::labelList&, const Foam::DynamicList<int>&, bool) const in file AMIInterpolation/AMIInterpolation/faceAreaWeightAMI/faceAreaWeightAMI.C at line 363. FOAM aborting #0 Foam::error::printStack(Foam::Ostream&) at ??:? #1 Foam::error::simpleExit(int, bool) at ??:? #2 Foam::error::exiting(int, bool) at ??:? #3 Foam::faceAreaWeightAMI::setNextFaces(int&, int&, int&, Foam::bitSet const&, Foam::List<int>&, Foam::DynamicList<int, 16> const&, bool) const at ??:? #4 Foam::faceAreaWeightAMI::calcAddressing(Foam::List<Foam::DynamicList<int, 16> >&, Foam::List<Foam::DynamicList<double, 16> >&, Foam::List<Foam::DynamicList<Foam::Vector<double>, 16> >&, Foam::List<Foam::DynamicList<int, 16> >&, Foam::List<Foam::DynamicList<double, 16> >&, int, int) at ??:? #5 Foam::faceAreaWeightAMI::calculate(Foam::PrimitivePatch<Foam::SubList<Foam::face>, Foam::Field<Foam::Vector<double> > const&> const&, Foam::PrimitivePatch<Foam::SubList<Foam::face>, Foam::Field<Foam::Vector<double> > const&> const&, Foam::autoPtr<Foam::searchableSurface> const&) at ??:? #6 Foam::cyclicAMIPolyPatch::resetAMI(Foam::UList<Foam::Vector<double> > const&) const at ??:? #7 Foam::cyclicAMIPolyPatch::AMI() const at ??:? #8 Foam::cyclicAMIPolyPatch::applyLowWeightCorrection() const at ??:? #9 Foam::cyclicAMIFvPatchField<double>::patchNeighbourField() const at ??:? #10 Foam::coupledFvPatchField<double>::evaluate(Foam::UPstream::commsTypes) at ??:? #11 Foam::cyclicAMIFvPatchField<double>::cyclicAMIFvPatchField(Foam::fvPatch const&, Foam::DimensionedField<double, Foam::volMesh> const&, Foam::dictionary const&) at ??:? #12 Foam::fvPatchField<double>::adddictionaryConstructorToTable<Foam::cyclicAMIFvPatchField<double> >::New(Foam::fvPatch const&, Foam::DimensionedField<double, Foam::volMesh> const&, Foam::dictionary const&) at ??:? #13 Foam::fvPatchField<double>::New(Foam::fvPatch const&, Foam::DimensionedField<double, Foam::volMesh> const&, Foam::dictionary const&) in ~/OpenFOAM/OpenFOAM-v2206/platforms/linux64GccDPInt32Opt/bin/decomposePar #14 Foam::GeometricBoundaryField<double, Foam::fvPatchField, Foam::volMesh>::readField(Foam::DimensionedField<double, Foam::volMesh> const&, Foam::dictionary const&) in ~/OpenFOAM/OpenFOAM-v2206/platforms/linux64GccDPInt32Opt/bin/decomposePar #15 Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh>::readFields(Foam::dictionary const&) in ~/OpenFOAM/OpenFOAM-v2206/platforms/linux64GccDPInt32Opt/bin/decomposePar #16 Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh>::readFields() in ~/OpenFOAM/OpenFOAM-v2206/platforms/linux64GccDPInt32Opt/bin/decomposePar #17 Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh>::GeometricField(Foam::IOobject const&, Foam::fvMesh const&, bool) in ~/OpenFOAM/OpenFOAM-v2206/platforms/linux64GccDPInt32Opt/bin/decomposePar #18 void Foam::fieldsDistributor::readFields<double, Foam::fvPatchField, Foam::volMesh>(Foam::volMesh::Mesh const&, Foam::IOobjectList const&, Foam::PtrList<Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> >&, bool) at ??:? #19 Foam::fvFieldDecomposer::fieldsCache::readAllFields(Foam::fvMesh const&, Foam::IOobjectList const&) at ??:? #20 ? in ~/OpenFOAM/OpenFOAM-v2206/platforms/linux64GccDPInt32Opt/bin/decomposePar #21 __libc_start_main in /lib/x86_64-linux-gnu/libc.so.6 #22 ? in ~/OpenFOAM/OpenFOAM-v2206/platforms/linux64GccDPInt32Opt/bin/decomposePar请问各位老师有没有什么解决方案🙏
-
请教贴：喷管模拟压力与实验压力值相差较大@李东岳 是的，通过理论计算，出口应该在1.4Ma，fluent模拟的结果也在1.4Ma。但是我模拟的结果出口却达到了4Ma
-
请教贴：喷管模拟压力与实验压力值相差较大各位大佬好，我在模拟喷管的时候，采用改进的rhoCentralFoam求解器，模拟出来的出口压力在0.9MPa，但实验出口压力2.5MPa，相差较大。
我发现喷管出口如果设置fixedValue无法计算。但采用waveTransmissive又无法固定出口压力。以下是我边界条件
想请教一下大家有没有啥更改的建议 -
模拟中的波动问题嗯嗯，我尝试一下
-
模拟中的波动问题李老师，求解器里面的热物性库改了，求解器也做了一些改动，不太能上传。
但我通过第一次波动前后的温度和压力进行计算，发现波动后面的计算取值都是液相的物性值。因为我觉的问题出现在我自定义的液相区物性。
-
模拟中的波动问题请教贴：
各位老师，我在模拟喷管的时候采用的rhoCentralFoam求解器，模拟结果出现这种波动，中心线上的压力分布也是跳跃样式，并没有平滑下降。想请教一下各位老师，出现这个情况的原因可能有什么？喷管压力云图：
中心线压力分布：
-
请教贴：heThermo类中属性he_是内能还是焓的问题和一位师兄交流了一下。是我的判断搞错了，在heThermo类中he_属性取决于energy的选择。如果energy 是sensibleEnthalpy，那么he_=h;如果energy是sensibleInternalEnergy，那么he_=e。对于rhoCentralFoam求解的控制方程是总能，计算的是e。所以energy只能使用sensibleInternalEnergy。
我去openFOAM带的rhoCentralFoam案例库查了一下,确实都是sensibleInternalEnergy！！！ -
请教贴：heThermo类中属性he_是内能还是焓的问题在rhoCentralFoam求解器中，存在以下代码
solve ( fvm::ddt(rhoE) + fvc::div(phiEp) - fvc::div(sigmaDotU) ); dimensionedScalar Emin("Emin", dimEnergy/dimMass, 71800); e =max(rhoE/rho - 0.5*magSqr(U),Emin); e.correctBoundaryConditions();控制方程求解出的是总能，获得内能e。
热物性表constant/thermophysicalProperties中设置的
thermo 是hConst，这意味着he_应该是比焓
energy 设置应该是把he_传递到sensibleEnthalpy类中，使焓来转化为温度感觉这个合理；
但thermo 是hConst
energy设置为sensibleInternalEnergy，这意味着he_应该是焓，传递给sensibleInternalEnergy来转化温度，感觉不是很合理呀。这个地方有理解的老师，能给讲讲嘛，已经懵了 -
请教：可压缩流，压力出口边界条件设置问题各位老师好。在使用rhoCentralFoam模拟喷管的时候遇到一个问题。
压力出口使用waveTransmissive边界条件，保证波能穿过出口边界，但模拟结果显示的出口压力与想要的出口压力相差很大。
即使想要一个高背压的环境，但waveTransmissive依旧会把压力出口计算为低的背压环境。
想请教一下各位老师，在模拟喷管的时候，想要看到由于高背压造成在喷管渐阔段产生的激波，这类模拟出口边界条件应该如何设置。 -
请教：关于openfoam中添加源项的问题在fluent中也存在这种源项线性化，目的是增强求解稳定性。
源项通常表示为S_phi=A+B*phi的形式，其中phi表示因变量，A为源项的显式部分，B*phi为隐式部分A=S_phi-d(S)/d(phi)*phi B=d(S)/d(phi)参考：《ansys fluent 二次开发指南》
rhoU做为因变量，Source_U_linear_explicit 和fvm::SuSp(SourceMomentumFactor, rhoU)大概也是为了形成A和B。 -
OpenFOAM膨胀喷管边界条件设置温度和压力的边界条件与内场压力差非常大。要不要试着从小的入口压力开始跑，稳定之后在增加压力。这种操作方法，我在别人的喷管案例里面看见过
-
请教：关于openfoam中添加源项的问题volVectorField Source_U("Source_U", -U*Source_Y); volVectorField Source_U_linear_explicit("Source_U_linear_explicit", SourceMomentumFactor*rhoU); solve ( fvm::ddt(rhoU) + fvc::div(phiUp) == Source_U - Source_U_linear_explicit + fvm::SuSp(SourceMomentumFactor, rhoU) );在看别人程序的时候发现在源项的添加上多出来两项Source_U_linear_explicit 和fvm::SuSp(SourceMomentumFactor, rhoU)，这两感觉是一样的，不是很理解这样写的原因。哪位老师，能指导一下
-
求教：Laval 喷管模拟产生激波边界条件设置问题李老师，感谢您的修改，对比后我发现我对fvsolutions的设置还不太明确。用您改后的案例，openfoamv9跑出来了，但用openfoam-v2206还是报原来的错误。目前，我准备尝试用clion调试找一下壁面函数中wallDist的问题
-
求教：Laval 喷管模拟产生激波边界条件设置问题网址上的案例看了，使用的是层流模型。我使用层流模型也可以跑通。但使用RAS模型加入壁面函数出问题，麻烦李老师给我看看，debug一下
EosNozzleTestUpdate.zip -
求教：Laval 喷管模拟产生激波边界条件设置问题OF9尝试了，还是避免函数有问题，结果如下：
Selecting thermodynamics package { type heRhoThermo; mixture pureMixture; transport const; thermo hConst; equationOfState perfectGas; specie specie; energy sensibleInternalEnergy; } Reading field U Reading/calculating face flux field phi Creating turbulence model Selecting turbulence model type RAS Selecting RAS turbulence model kEpsilon RAS { model kEpsilon; turbulence on; printCoeffs on; Cmu 0.09; C1 1.44; C2 1.92; C3 0; sigmak 1; sigmaEps 1.3; } Creating thermophysical transport model Selecting thermophysical transport type RAS Selecting default RAS thermophysical transport model unityLewisEddyDiffusivity Creating field dpdt Creating field kinetic energy K No MRF models present No fvModels present No fvConstraints present #0 Foam::error::printStack(Foam::Ostream&) at ??:? #1 Foam::sigFpe::sigHandler(int) at ??:? #2 ? in "/lib/x86_64-linux-gnu/libc.so.6" #3 Foam::nutkWallFunctionFvPatchScalarField::nut() const at ??:? #4 Foam::nutWallFunctionFvPatchScalarField::updateCoeffs() at ??:? #5 Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh>::Boundary::evaluate() in "/home/tsing/OpenFOAM/OpenFOAM-9/platforms/linux64GccDPInt32Opt/bin/rhoPimpleFoam" #6 Foam::RASModels::kEpsilon<Foam::CompressibleMomentumTransportModel<Foam::dynamicTransportModel> >::correctNut() at ??:? #7 ? in "/home/tsing/OpenFOAM/OpenFOAM-9/platforms/linux64GccDPInt32Opt/bin/rhoPimpleFoam" #8 __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6" #9 ? in "/home/tsing/OpenFOAM/OpenFOAM-9/platforms/linux64GccDPInt32Opt/bin/rhoPimpleFoam" Floating point exception (core dumped) -
求教：Laval 喷管模拟产生激波边界条件设置问题换成层流模型可以。但换成kEpsilon模型，也是在Epsilon壁面函数错误。输出结果如下所示。
PIMPLE: Operating solver in PISO mode Reading thermophysical properties Selecting thermodynamics package { type heRhoThermo; mixture pureMixture; transport const; thermo hConst; equationOfState perfectGas; specie specie; energy sensibleInternalEnergy; } Reading field U Reading/calculating face flux field phi Creating turbulence model Selecting turbulence model type RAS Selecting RAS turbulence model kEpsilon RAS { RASModel kEpsilon; turbulence on; printCoeffs on; Cmu 0.09; C1 1.44; C2 1.92; C3 0; sigmak 1; sigmaEps 1.3; } Creating field dpdt Creating field kinetic energy K No MRF models present No finite volume options present Courant Number mean: 4.30558e-05 max: 0.0192008 Starting time loop Courant Number mean: 4.30558e-05 max: 0.0192008 Time = 5e-05 PIMPLE: iteration 1 diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0 smoothSolver: Solving for Ux, Initial residual = 1, Final residual = 4.88632e-08, No Iterations 1 smoothSolver: Solving for Uy, Initial residual = 1, Final residual = 2.83045e-08, No Iterations 1 smoothSolver: Solving for e, Initial residual = 1, Final residual = 3.20771e-07, No Iterations 3 GAMG: Solving for p, Initial residual = 1, Final residual = 0.00612066, No Iterations 3 GAMG: Solving for p, Initial residual = 0.00680187, Final residual = 3.35179e-05, No Iterations 4 diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0 time step continuity errors : sum local = 7.89239e-08, global = -7.08945e-09, cumulative = -7.08945e-09 GAMG: Solving for p, Initial residual = 0.147326, Final residual = 0.000684668, No Iterations 3 GAMG: Solving for p, Initial residual = 0.00064622, Final residual = 4.92284e-07, No Iterations 5 diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0 time step continuity errors : sum local = 1.8365e-09, global = -5.67277e-10, cumulative = -7.65673e-09 #0 Foam::error::printStack(Foam::Ostream&) at ??:? #1 Foam::sigFpe::sigHandler(int) at ??:? #2 ? in /lib/x86_64-linux-gnu/libpthread.so.0 #3 Foam::epsilonWallFunctionFvPatchScalarField::calculate(Foam::turbulenceModel const&, Foam::List<double> const&, Foam::fvPatch const&, Foam::Field<double>&, Foam::Field<double>&) at ??:? #4 Foam::epsilonWallFunctionFvPatchScalarField::calculateTurbulenceFields(Foam::turbulenceModel const&, Foam::Field<double>&, Foam::Field<double>&) at ??:? #5 Foam::epsilonWallFunctionFvPatchScalarField::updateCoeffs() at ??:? #6 Foam::RASModels::kEpsilon<Foam::EddyDiffusivity<Foam::ThermalDiffusivity<Foam::CompressibleTurbulenceModel<Foam::fluidThermo> > > >::correct() at ??:? #7 ? in ~/OpenFOAM/OpenFOAM-v2206/platforms/linux64GccDPInt32Opt/bin/rhoPimpleFoam #8 __libc_start_main in /lib/x86_64-linux-gnu/libc.so.6 #9 ? in ~/OpenFOAM/OpenFOAM-v2206/platforms/linux64GccDPInt32Opt/bin/rhoPimpleFoam Floating point exception (core dumped) -
求教：Laval 喷管模拟产生激波边界条件设置问题想要设置喷管，入口总压102387.1Pa，出口静压42449.07Pa, 入口温度294.45K。采用kOmegaSST湍流模型。但是运行报告湍流边界条件设置有问题。改了好几次都不行，希望各位老师，能帮我看一下，看看是哪里设置存在问题。以下是我的初始条件：
U T p omega nut Ma k alphat
运行报错结果：tsing@tsing:~/OpenFOAM/tsing-v2206/run/EosNozzlePerfect$ rhoPimpleFoam /*---------------------------------------------------------------------------*\ | ========= | | | \\ / F ield | OpenFOAM: The Open Source CFD Toolbox | | \\ / O peration | Version: 2206 | | \\ / A nd | Website: www.openfoam.com | | \\/ M anipulation | | \*---------------------------------------------------------------------------*/ Build : _76d719d1e6-20220624 OPENFOAM=2206 version=v2206 Arch : "LSB;label=32;scalar=64" Exec : rhoPimpleFoam Date : Feb 24 2025 Time : 18:29:00 Host : tsing PID : 218467 I/O : uncollated Case : /home/tsing/OpenFOAM/tsing-v2206/run/EosNozzlePerfect nProcs : 1 trapFpe: Floating point exception trapping enabled (FOAM_SIGFPE). fileModificationChecking : Monitoring run-time modified files using timeStampMaster (fileModificationSkew 5, maxFileModificationPolls 20) allowSystemOperations : Allowing user-supplied system call operations // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // Create time Create mesh for time = 0 PIMPLE: Operating solver in PISO mode Reading thermophysical properties Selecting thermodynamics package { type heRhoThermo; mixture pureMixture; transport const; thermo hConst; equationOfState perfectGas; specie specie; energy sensibleInternalEnergy; } Reading field U Reading/calculating face flux field phi Creating turbulence model Selecting turbulence model type RAS Selecting RAS turbulence model kOmegaSST Selecting patchDistMethod meshWave RAS { RASModel kOmegaSST; turbulence on; printCoeffs on; alphaK1 0.85; alphaK2 1; alphaOmega1 0.5; alphaOmega2 0.856; gamma1 0.555556; gamma2 0.44; beta1 0.075; beta2 0.0828; betaStar 0.09; a1 0.31; b1 1; c1 10; F3 false; decayControl false; kInf 0; omegaInf 0; } Creating field dpdt Creating field kinetic energy K No MRF models present No finite volume options present Courant Number mean: 4.30558e-05 max: 0.0192008 Starting time loop Courant Number mean: 4.30558e-05 max: 0.0192008 Time = 5e-05 PIMPLE: iteration 1 diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0 DILUPBiCGStab: Solving for Ux, Initial residual = 1, Final residual = 8.37036e-10, No Iterations 1 DILUPBiCGStab: Solving for Uy, Initial residual = 1, Final residual = 4.3098e-10, No Iterations 1 DILUPBiCGStab: Solving for e, Initial residual = 1, Final residual = 4.75567e-09, No Iterations 1 GAMG: Solving for p, Initial residual = 1, Final residual = 0.00613491, No Iterations 3 GAMG: Solving for p, Initial residual = 0.00680296, Final residual = 3.5161e-05, No Iterations 4 diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0 time step continuity errors : sum local = 8.28493e-08, global = -7.8312e-09, cumulative = -7.8312e-09 GAMG: Solving for p, Initial residual = 0.147144, Final residual = 0.000685025, No Iterations 3 GAMG: Solving for p, Initial residual = 0.000646506, Final residual = 5.05846e-07, No Iterations 5 diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0 time step continuity errors : sum local = 1.88756e-09, global = -4.94043e-10, cumulative = -8.32524e-09 #0 Foam::error::printStack(Foam::Ostream&) at ??:? #1 Foam::sigFpe::sigHandler(int) at ??:? #2 ? in /lib/x86_64-linux-gnu/libpthread.so.0 #3 Foam::omegaWallFunctionFvPatchScalarField::calculate(Foam::turbulenceModel const&, Foam::List<double> const&, Foam::fvPatch const&, Foam::Field<double>&, Foam::Field<double>&) at ??:? #4 Foam::omegaWallFunctionFvPatchScalarField::calculateTurbulenceFields(Foam::turbulenceModel const&, Foam::Field<double>&, Foam::Field<double>&) at ??:? #5 Foam::omegaWallFunctionFvPatchScalarField::updateCoeffs() at ??:? #6 Foam::kOmegaSSTBase<Foam::eddyViscosity<Foam::RASModel<Foam::EddyDiffusivity<Foam::ThermalDiffusivity<Foam::CompressibleTurbulenceModel<Foam::fluidThermo> > > > > >::correct() at ??:? #7 ? in ~/OpenFOAM/OpenFOAM-v2206/platforms/linux64GccDPInt32Opt/bin/rhoPimpleFoam #8 __libc_start_main in /lib/x86_64-linux-gnu/libc.so.6 #9 ? in ~/OpenFOAM/OpenFOAM-v2206/platforms/linux64GccDPInt32Opt/bin/rhoPimpleFoam Floating point exception (core dumped)