|
[Sponsors] |
![]() |
![]() |
#1 |
New Member
Anton Afanassiev
Join Date: May 2018
Posts: 1
Rep Power: 0 ![]() |
Hi,
I have been trying to simulate a NREL 5MW reference turbine using rotorDiskSource. I seem to get good results when running on a single processor. But when I decompose and run, the turbine's wake completely disappears. I've noticed that in the decomposed case the rotorDiskSource output shows effective lift = 0, when in the single processor case it is around 100k. Code:
otorDisk output: min/max(AOA) = 1.65132, 1.71517 Effective drag = -5011.09 Effective lift = 0 GAMG: Solving for Ux, Initial residual = 1.00076e-06, Final residual = 1.00076e-06, No Iterations 0 GAMG: Solving for Uy, Initial residual = 7.42757e-06, Final residual = 7.42757e-06, No Iterations 0 GAMG: Solving for Uz, Initial residual = 7.49477e-06, Final residual = 7.49477e-06, No Iterations 0 GAMG: Solving for T, Initial residual = 2.87311e-06, Final residual = 2.87311e-06, No Iterations 0 GAMG: Solving for p_rgh, Initial residual = 0.000266345, Final residual = 2.55409e-06, No Iterations 18 time step continuity errors : sum local = 1.94881e-10, global = 1.05977e-15, cumulative = -1.00503e-09 GAMG: Solving for p_rgh, Initial residual = 6.00008e-06, Final residual = 9.16971e-07, No Iterations 6 time step continuity errors : sum local = 6.99648e-11, global = -5.42417e-15, cumulative = -1.00503e-09 ExecutionTime = 17.06 s ClockTime = 18 s Code:
#include "../meshSetup" numberOfSubdomains $nDomains; // method scotch; // method hierarchical; // method simple; method $method; // method metis; // method manual; // method multiLevel; // method structured; // does 2D decomposition of structured mesh multiLevelCoeffs { // Decomposition methods to apply in turn. This is like hierarchical but // fully general - every method can be used at every level. level0 { numberOfSubdomains 64; //method simple; //simpleCoeffs //{ // n (2 1 1); // delta 0.001; //} method scotch; } level1 { numberOfSubdomains 4; method scotch; } } // Desired output simpleCoeffs { n ($Sx $Sy $Sz); delta 0.001; } hierarchicalCoeffs { n (1 2 1); delta 0.001; order xyz; } metisCoeffs { /* processorWeights ( 1 1 1 1 ); */ } scotchCoeffs { //processorWeights //( // 1 // 1 // 1 // 1 //); //writeGraph true; //strategy "b"; } manualCoeffs { dataFile "decompositionData"; } structuredCoeffs { // Patches to do 2D decomposition on. Structured mesh only; cells have // to be in 'columns' on top of patches. patches (bottomPatch); } Below is the relevant portion of my fvOptions: Code:
disk { type rotorDisk; active on; rotorDiskCoeffs { selectionMode cellSet; cellSet turbine1; fields (U); // Names of fields on which to apply source nBlades 3; // Number of blades tipEffect 0.96; // Normalised radius above which lift = 0 inletFlowType local; // Inlet flow type specification geometryMode specified; origin (0 0 0); axis (1 0 0); refDirection (0 0 1); // Reference direction // - used as reference for psi angle rpm 12.1; trimModel fixedTrim; // fixed || targetForce rhoRef 1.225; rhoInf 1.225; fixedTrimCoeffs { theta0 0; theta1c 0; theta1s 0; } flapCoeffs { beta0 0; // Coning angle [deg] beta1c 0; // Lateral flapping coeff (cos coeff) beta2s 0; // Longitudinal flapping coeff (sin coeff) } |
|
![]() |
![]() |
![]() |
![]() |
#2 |
Member
Kisorthman Vimalakanthan
Join Date: Apr 2011
Posts: 49
Rep Power: 15 ![]() |
Hi Jabbath,
I'm faced with the same issue? were you able to simulate the NREL5MW rotor with rotorDiskSource on parallel.. could you please share your files? many thanks in advance, Kind regards, Kishore
__________________
Kisorthman Vimalakanthan Dept. of Power and Propulsion Cranfield University Email: k.vimalakanthan@gmail.com |
|
![]() |
![]() |
![]() |
![]() |
#3 |
Senior Member
Herpes Free Engineer
Join Date: Sep 2019
Location: The Home Under The Ground with the Lost Boys
Posts: 931
Rep Power: 12 ![]() |
Hi,
Could you please open a bug-issue with a minimal reproducible example in https://develop.openfoam.com/? Many thanks.
__________________
The OpenFOAM community is the biggest contributor to OpenFOAM: User guide/Wiki-1/Wiki-2/Code guide/Code Wiki/Journal Nilsson/Guerrero/Holzinger/Holzmann/Nagy/Santos/Nozaki/Jasak/Primer Governance Bugs/Features: OpenFOAM (ESI-OpenCFD-Trademark) Bugs/Features: FOAM-Extend (Wikki-FSB) Bugs: OpenFOAM.org How to create a MWE New: Forkable OpenFOAM mirror |
|
![]() |
![]() |
![]() |
![]() |
#4 | |
Member
Kisorthman Vimalakanthan
Join Date: Apr 2011
Posts: 49
Rep Power: 15 ![]() |
Quote:
__________________
Kisorthman Vimalakanthan Dept. of Power and Propulsion Cranfield University Email: k.vimalakanthan@gmail.com |
||
![]() |
![]() |
![]() |
Tags |
decomposepar, openfoam 4.1, wind turbine blade, wind turbines |
|
|
![]() |
||||
Thread | Thread Starter | Forum | Replies | Last Post |
Is Playstation 3 cluster suitable for CFD work | hsieh | OpenFOAM | 9 | August 16, 2015 14:53 |
Discrepancy ensight slice and full data for parallel decomposed case | Japsor | OpenFOAM Bugs | 0 | August 27, 2013 04:47 |
Error reading new case | montag dp | FLUENT | 5 | September 15, 2011 06:00 |
decomposed case to 2-cores (Not working) | pkr | OpenFOAM Running, Solving & CFD | 32 | February 27, 2011 18:04 |
Turbulent Flat Plate Validation Case | Jonas Larsson | Main CFD Forum | 0 | April 2, 2004 10:25 |