|
[Sponsors] |
![]() |
![]() |
#1 |
New Member
Pouria Taghikhani
Join Date: Sep 2016
Posts: 14
Rep Power: 9 ![]() |
Hi all
I am trying to simulate a ducted propeller with pimpleDyMFoam solver. There is no problem with propeller alone. But when I try to add the duct and solve it again, solver stops. I am now trying to figure the problem by using: Code:
runApplication moveDynamicMesh -checkAMI but the following message apears: Code:
/*---------------------------------------------------------------------------*\ | ========= | | | \\ / F ield | OpenFOAM: The Open Source CFD Toolbox | | \\ / O peration | Version: 4.x | | \\ / A nd | Web: www.OpenFOAM.org | | \\/ M anipulation | | \*---------------------------------------------------------------------------*/ /* Windows 32 and 64 bit porting by blueCAPE: http://www.bluecape.com.pt *\ | Based on Windows porting (2.0.x v4) by Symscape: http://www.symscape.com | \*---------------------------------------------------------------------------*/ Build : 4.x-ed69f631ce88 Exec : C:/PROGRA~1/BLUECF~1/OpenFOAM-4.x/platforms/mingw_w64GccDPInt32Opt/bin/moveDynamicMesh.exe -checkAMI Date : May 14 2017 Time : 21:54:08 Host : "DESKTOP-PU3N81D" PID : 5112 Case : C:/Users/Pouria/Desktop/propeller400 nProcs : 1 SigFpe : Enabling floating point exception trapping (FOAM_SIGFPE). fileModificationChecking : Monitoring run-time modified files using timeStampMaster allowSystemOperations : Allowing user-supplied system call operations // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // Create time Create mesh for time = 0 Selecting dynamicFvMesh solidBodyMotionFvMesh Selecting solid-body motion function rotatingMotion Applying solid body motion to cellZone innerCylinderSmall Writing VTK files with weights of AMI patches. PIMPLE: no residual control data found. Calculations will employ 2 corrector loops forces Rotor: --> FOAM Warning : From function Foam::labelHashSet Foam::polyBoundaryMesh::patchSet(const Foam::UList<Foam::wordRe>&, bool, bool) const in file meshes/polyMesh/polyBoundaryMesh/polyBoundaryMesh.C at line 794 Cannot find any patch or group names matching propellerRotor Not including porosity effects forces DAS: --> FOAM Warning : From function Foam::labelHashSet Foam::polyBoundaryMesh::patchSet(const Foam::UList<Foam::wordRe>&, bool, bool) const in file meshes/polyMesh/polyBoundaryMesh/polyBoundaryMesh.C at line 794 Cannot find any patch or group names matching propellerDuct Not including porosity effects Time = 5e-006 PIMPLE: iteration 1 --> FOAM Warning : From function virtual bool Foam::solidBodyMotionFvMesh::update() in file solidBodyMotionFvMesh/solidBodyMotionFvMesh.C at line 228 Did not find volVectorField U Not updating Uboundary conditions. PIMPLE: iteration 2 Point usage OK. Upper triangular ordering OK. Topological cell zip-up check OK. Face vertices OK. <<Number of duplicate (not baffle) faces found: 2. This might indicate a problem. <<Number of faces with non-consecutive shared points: 3. This might indicate a problem. Mesh topology OK. Boundary openness (3.72253e-016 -9.13045e-017 9.01232e-017) OK. Max cell openness = 4.35445e-016 OK. Max aspect ratio = 11.8457 OK. Minimum face area = 3.57658e-009. Maximum face area = 0.00294544. Face area magnitudes OK. Min volume = 9.05772e-011. Max volume = 0.000110349. Total volume = 1.31244. Cell volumes OK. Mesh non-orthogonality Max: 64.9606 average: 9.12233 Non-orthogonality check OK. Face pyramids OK. Max skewness = 3.94105 OK. Mesh geometry OK. Mesh OK. Calculating AMI weights between owner patch: AMI1 and neighbour patch: AMI2 AMI: Creating addressing and weights between 45130 source faces and 46232 target faces AMI: Patch source sum(weights) min/max/average = 0.0734142, 2.13295, 0.999935 AMI: Patch target sum(weights) min/max/average = 0, 1.84503, 0.999232 ExecutionTime = 39.753 s ClockTime = 40 s --> FOAM FATAL ERROR: Could not find U, p From function void Foam::functionObjects::forces::initialise() in file forces/forces.C at line 186. FOAM exiting Code:
/*---------------------------------------------------------------------------*\ | ========= | | | \\ / F ield | OpenFOAM: The Open Source CFD Toolbox | | \\ / O peration | Version: 4.1 | | \\ / A nd | Web: www.OpenFOAM.org | | \\/ M anipulation | | \*---------------------------------------------------------------------------*/ Build : 4.1 Exec : pimpleDyMFoam -parallel Date : May 15 2017 Time : 12:40:46 Host : "pouria" PID : 3167 Case : /home/pouria/Desktop/propeller400 nProcs : 4 Slaves : 3 ( "pouria.3168" "pouria.3169" "pouria.3170" ) Pstream initialized with: floatTransfer : 0 nProcsSimpleSum : 0 commsType : nonBlocking polling iterations : 0 sigFpe : Enabling floating point exception trapping (FOAM_SIGFPE). fileModificationChecking : Monitoring run-time modified files using timeStampMaster allowSystemOperations : Allowing user-supplied system call operations // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // Create time Create mesh for time = 0 Selecting dynamicFvMesh solidBodyMotionFvMesh Selecting solid-body motion function rotatingMotion Applying solid body motion to cellZone innerCylinderSmall PIMPLE: no residual control data found. Calculations will employ 2 corrector loops Reading field p Reading field U Reading/calculating face flux field phi AMI: Creating addressing and weights between 45146 source faces and 46392 target faces AMI: Patch source sum(weights) min/max/average = 0.565692, 2.10459, 1.00072 AMI: Patch target sum(weights) min/max/average = 0, 1.84453, 0.998807 Selecting incompressible transport model Newtonian Selecting turbulence model type RAS Selecting RAS turbulence model kEpsilon bounding k, min: 0 max: 0.06 average: 0.06 bounding epsilon, min: 0 max: 0.0495 average: 0.0495 kEpsilonCoeffs { Cmu 0.09; C1 1.44; C2 1.92; C3 -0.33; sigmak 1; sigmaEps 1.3; } No MRF models present Reading/calculating face velocity Uf No finite volume options present Courant Number mean: 6.98835e-07 max: 0.00055208 forces propellerTip: Not including porosity effects forces DAS: Not including porosity effects Starting time loop Courant Number mean: 6.98835e-07 max: 0.00055208 deltaT = 1e-05 Time = 1e-05 AMI: Creating addressing and weights between 45146 source faces and 46392 target faces AMI: Patch source sum(weights) min/max/average = 0.00638659, 2.10717, 1.00071 AMI: Patch target sum(weights) min/max/average = 0, 1.84552, 0.998866 PIMPLE: iteration 1 smoothSolver: Solving for Ux, Initial residual = 1, Final residual = 0.0302785, No Iterations 1 smoothSolver: Solving for Uy, Initial residual = 1, Final residual = 0.000126242, No Iterations 1 smoothSolver: Solving for Uz, Initial residual = 1, Final residual = 0.0305505, No Iterations 1 [1] #0 Foam::error::printStack(Foam::Ostream&) at ??:? [1] #1 Foam::sigFpe::sigHandler(int) at ??:? [1] #2 ? in "/lib/x86_64-linux-gnu/libc.so.6" [1] #3 Foam::divide(Foam::Field<double>&, double const&, Foam::UList<double> const&) at ??:? [1] #4 Foam::tmp<Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> > Foam::operator/<Foam::fvPatchField, Foam::volMesh>(Foam::dimensioned<double> const&, Foam::tmp<Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> > const&) at ??:? [1] #5 ? at ??:? [1] #6 __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6" [1] #7 ? at ??:? [pouria:03168] *** Process received signal *** [pouria:03168] Signal: Floating point exception (8) [pouria:03168] Signal code: (-6) [pouria:03168] Failing at address: 0x3e800000c60 [pouria:03168] [ 0] /lib/x86_64-linux-gnu/libc.so.6(+0x354b0)[0x7f15358094b0] [pouria:03168] [ 1] /lib/x86_64-linux-gnu/libc.so.6(gsignal+0x38)[0x7f1535809428] [pouria:03168] [ 2] /lib/x86_64-linux-gnu/libc.so.6(+0x354b0)[0x7f15358094b0] [pouria:03168] [ 3] /opt/openfoam4/platforms/linux64GccDPInt32Opt/lib/libOpenFOAM.so(_ZN4Foam6divideERNS_5FieldIdEERKdRKNS_5UListIdEE+0x114)[0x7f1536c46264] [pouria:03168] [ 4] pimpleDyMFoam(_ZN4FoamdvINS_12fvPatchFieldENS_7volMeshEEENS_3tmpINS_14GeometricFieldIdT_T0_EEEERKNS_11dimensionedIdEERKS8_+0x1d7)[0x474ba7] [pouria:03168] [ 5] pimpleDyMFoam[0x4259ac] [pouria:03168] [ 6] /lib/x86_64-linux-gnu/libc.so.6(__libc_start_main+0xf0)[0x7f15357f4830] [pouria:03168] [ 7] pimpleDyMFoam[0x4286c9] [pouria:03168] *** End of error message *** -------------------------------------------------------------------------- mpirun noticed that process rank 1 with PID 3168 on node pouria exited on signal 8 (Floating point exception). -------------------------------------------------------------------------- by the way. I am uploading a section of my project geometry. The purple block shows AMI region. thank you. Last edited by Taghi; May 15, 2017 at 04:25. |
|
![]() |
![]() |
![]() |
![]() |
#2 |
Super Moderator
Tobias Holzmann
Join Date: Oct 2010
Location: Bad Wörishofen
Posts: 2,711
Blog Entries: 6
Rep Power: 51 ![]() ![]() ![]() |
Hi,
for your moveDynamicMesh application you should call the following: Code:
moveDynamicMesh -noFunctionObjects -checkAMI
__________________
Keep foaming, Tobias Holzmann |
|
![]() |
![]() |
![]() |
![]() |
#3 |
New Member
Pouria Taghikhani
Join Date: Sep 2016
Posts: 14
Rep Power: 9 ![]() |
Thanks. It worked. But I still have problems with pimplyDyMFoam
|
|
![]() |
![]() |
![]() |
![]() |
#4 |
Super Moderator
Tobias Holzmann
Join Date: Oct 2010
Location: Bad Wörishofen
Posts: 2,711
Blog Entries: 6
Rep Power: 51 ![]() ![]() ![]() |
Is your AMI moving correctly and are the BC correct? You have a floating point exception (division by zero). Maybe you missed some
Code:
value uniform xxx;
__________________
Keep foaming, Tobias Holzmann |
|
![]() |
![]() |
![]() |
![]() |
#5 |
New Member
Pouria Taghikhani
Join Date: Sep 2016
Posts: 14
Rep Power: 9 ![]() |
Hi again. It took a long time but I think I figured were the problem is. The AMI is too close to Boundaries. Is there any way that I can rotate my propeller without AMI. MRF of SRF?
|
|
![]() |
![]() |
![]() |
Tags |
moving mesh, pimpledymfoam, propellers |
|
|
![]() |
||||
Thread | Thread Starter | Forum | Replies | Last Post |
[swak4Foam] GroovyBC the dynamic cousin of funkySetFields that lives on the suburb of the mesh | gschaider | OpenFOAM Community Contributions | 300 | October 29, 2014 18:00 |
How to install CGNS under windows xp? | lzgwhy | Main CFD Forum | 1 | January 11, 2011 18:44 |
[swak4Foam] groovyBC: problems compiling: "flex: not found" and "undefined reference to ..." | sega | OpenFOAM Community Contributions | 12 | February 17, 2010 09:30 |
How to get the max value of the whole field | waynezw0618 | OpenFOAM Running, Solving & CFD | 4 | June 17, 2008 05:07 |
Compiling problems with hello worldC | fw407 | OpenFOAM Installation | 21 | January 6, 2008 17:38 |