How to resolve the error 'Cutting back alpha due to evaluation error'? #775
Unanswered
sento1suzu8
asked this question in
Bug
Replies: 1 comment 1 reply
-
|
Check this post: #765 |
Beta Was this translation helpful? Give feedback.
1 reply
-
|
Thank you, Dr. He! Increasing the scale of the design variables has resolved most of the issues, but there are still a few cases where errors occur no matter how much the scale is increased. |
Beta Was this translation helpful? Give feedback.
Sign up for free
to join this conversation on GitHub.
Already have an account?
Sign in to comment
-
When I run airfoil optimization, I encounter the error: 'Cutting back alpha due to evaluation error'. It appears after solving the adjoint equation, where the mesh check failed, leading to this error. The process keeps looping until the mesh check is OK.The detailed context in which the error occurs is as follows:
Calculating [dRdW]^T * Psi using reverse-mode AD
Solving Linear Equation... 893 s
Main iteration 0 KSP Residual norm 8.271317064397e+01 893 s
Main iteration 100 KSP Residual norm 3.437881200480e+01 919 s
Main iteration 200 KSP Residual norm 1.166244343949e+01 947 s
Main iteration 300 KSP Residual norm 6.872823262659e+00 977 s
Main iteration 400 KSP Residual norm 5.419108065393e+00 1010 s
Main iteration 500 KSP Residual norm 4.612946765591e+00 1046 s
Main iteration 600 KSP Residual norm 3.148501385521e-01 1085 s
Main iteration 700 KSP Residual norm 3.323126128814e-02 1127 s
Main iteration 800 KSP Residual norm 1.336409716420e-02 1171 s
Main iteration 900 KSP Residual norm 1.023194960931e-02 1216 s
Main iteration 1000 KSP Residual norm 5.312335204715e-04 1266 s
Main iteration 1100 KSP Residual norm 1.616053211267e-04 1291 s
Main iteration 1198 KSP Residual norm 8.231506822385e-05 1318 s
Total iterations 1198
Solving Linear Equation... Completed! 1318 s
Residual tolerance satisfied, solution finished!
Calculating [dRdAOA]^T * Psi using reverse-mode AD
Calculating [dRdXv]^T * Psi using reverse-mode AD
Total number of variables............................: 129
variables with only lower bounds: 0
variables with lower and upper bounds: 129
variables with only upper bounds: 0
Total number of equality constraints.................: 72
Total number of inequality constraints...............: 42
inequality constraints with only lower bounds: 2
inequality constraints with lower and upper bounds: 40
inequality constraints with only upper bounds: 0
iter objective inf_pr inf_du lg(mu) ||d|| lg(rg) alpha_du alpha_pr ls
0 -6.2098204e+01 0.00e+00 7.67e+00 0.0 0.00e+00 - 0.00e+00 0.00e+00 0
Driver debug print for iter coord: rank0:pyOptSparse_IPOPT|2
Design Vars
{'dvs.aoa': array([2.006402956750587]),
'dvs.decalage': array([0.]),
'dvs.shape_low': array([ 0. , 0. , 0. ,
0. , 0.0592309205931247, 0.0592309205931247,
0.1089897309805795, 0.1089897309805795, 0.036081727683251 ,
0.036081727683251 , 0.0613595346707307, 0.0613595346707307,
0.0231740827505762, 0.0231740827505762, 0.1188000099 ,
0.1188000099 , 0.0094489236910009, 0.0094489236910009,
0.0340606311780038, 0.0340606311780038, 0.0055003955333538,
0.0055003955333538, 0.0481221582088257, 0.0481221582088257,
-0.0014036751432956, -0.0014036751432956, 0.0324944427098825,
0.0324944427098825, -0.0085340146917926, -0.0085340146917926,
0.017471512979505 , 0.017471512979505 , -0.0068474705191121,
-0.0068474705191121, 0.0051954299239382, 0.0051954299239382,
0.0326128650543985, 0.0326128650543985, 0.0434096741359753,
0.0434096741359753, 0. , 0. ,
0. , 0. ]),
'dvs.shape_up': array([ 0. , 0. , 0. ,
0. , 0.0681084428826185, 0.0681084428826185,
0.2135197356126622, 0.2135197356126622, 0.0583130180639486,
0.0583130180639486, 0.1499636659232295, 0.1499636659232295,
0.0403427478162987, 0.0403427478162987, 0.0863641812810199,
0.0863641812810199, -0.0048861858875737, -0.0048861858875737,
0.0630977930037021, 0.0630977930037021, -0.0293464721225661,
-0.0293464721225661, 0.0514873534675934, 0.0514873534675934,
-0.0141853376542055, -0.0141853376542055, 0.0449136599512155,
0.0449136599512155, -0.0410650153801825, -0.0410650153801825,
0.0336125884446796, 0.0336125884446796, -0.0088999485840658,
-0.0088999485840658, 0.0074242766236691, 0.0074242766236691,
-0.0401374379146979, -0.0401374379146979, 0.0488687719920184,
0.0488687719920184, -0.007436995498791 , -0.007436995498791 ,
0.0196902864096293, 0.0196902864096293, -0.0367144176604189,
-0.0367144176604189, 0.0455672742611457, 0.0455672742611457,
-0.0144046273258031, -0.0144046273258031, 0.0114905001028183,
0.0114905001028183, -0.035099104635883 , -0.035099104635883 ,
0.0313615665267206, 0.0313615665267206, -0.0227142600796746,
-0.0227142600796746, 0.0141990850251321, 0.0141990850251321,
-0.0241534007837717, -0.0241534007837717, 0.0208275843055889,
0.0208275843055889, -0.0301026562878543, -0.0301026562878543,
0.0129447013671949, 0.0129447013671949, -0.0372899729078434,
-0.0372899729078434, -0.01435094050669 , -0.01435094050669 ,
-0.046211364361681 , -0.046211364361681 , 0.0009555733942524,
0.0009555733942524, 0.0956543282040274, 0.0956543282040274,
0.114582304252481 , 0.114582304252481 , 0. ,
0. , 0. , 0. ]),
'dvs.translate': array([0., 0.])}
Setting up primal boundary conditions based on pyOptions:
Setting primal boundary conditions...
Setting U = (23.96529835758002 0.8395681224133503 0) at inlet
Setting primal boundary conditions...
Setting U = (23.96529835758002 0.8395681224133503 0) at bottomlet
Setting nut wall BC for upwing. BCType=nutLowReWallFunction
Setting nut wall BC for lowwing. BCType=nutLowReWallFunction
Setting k wall BC for upwing. BCType=fixedValue
Setting k wall BC for lowwing. BCType=fixedValue
Setting omega wall BC for upwing. BCType=omegaWallFunction
Setting omega wall BC for lowwing. BCType=omegaWallFunction
Running Primal Solver 003
Previous solution time 3000 not found and nothing deleted.
No finite volume options present
Starting time loop
Checking mesh quality for time = 0
Overall domain bounding box (-200.4110445027096 -200.3836137111711 0) (200.5416894763988 200.5069468205042 0.1000402611574791)
Mesh has 3 geometric (non-empty/wedge) directions 3{1}
Mesh has 3 solution (non-empty) directions 3{1}
Boundary openness (5.457718130858038e-20 -1.121410591104332e-19 3.300317064448473e-16) OK.
***High aspect ratio cells found, Max aspect ratio: 2.302979857130167e+97, number of cells 331
<<Writing 331 cells with high aspect ratio to set highAspectRatioCells
Minimum face area = 5.580738475076541e-10. Maximum face area = 61.00619349745944. Face area magnitudes OK.
***Zero or negative cell volume detected. Minimum negative volume: -1.421770034323735e-06, Number of negative volume cells: 331
<<Writing 331 zero volume cells to set zeroVolumeCells
Mesh non-orthogonality Max: 178.2931482739922 average: 21.20777616167234
*Number of severely non-orthogonal (> 75 degrees) faces: 627.
***Number of non-orthogonality errors: 686.
<<Writing 1313 non-orthogonal faces to set nonOrthoFaces
***Error in face pyramids: 2106 faces are incorrectly oriented.
<<Writing 1555 faces with incorrect orientation to set wrongOrientedFaces
***Max skewness = 44.13720311526393, 20 highly skew faces detected which may impair the quality of the results
<<Writing 20 skew faces to set skewFaces
Coupled point location match (average 0) OK.
Failed 5 mesh checks.
Nonlinear constraints
{'geometry.linearcon_low': array([0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0.,
0., 0., 0., 0., 0.]),
'geometry.linearcon_up': array([0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0.,
0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0.,
0., 0., 0., 0., 0., 0., 0., 0.]),
'geometry.thickcon_low': array([1.066142444830694 , 1.1871321808623083, 1.238803800693251 ,
1.1891742724133143, 1.1421375928359707, 1.135669744896856 ,
1.0985294921139213, 1.0848225366277606, 1.047234482155241 ,
0.9446011063063113, 1.1602228583030005, 1.1871321808623072,
1.2388038006932507, 1.1891742724133139, 1.14213759283597 ,
1.135669744896853 , 1.0985294921139221, 1.0848225366277586,
1.0472344821552355, 0.9446011063063113]),
'geometry.thickcon_up': array([1.0571631464417859, 1.1568701753503274, 1.1544146029154017,
1.1367586548474078, 1.1375927085355826, 1.1307852105853577,
1.1116869413080803, 1.0904766573882994, 1.075080034134348 ,
1.0051290390019074, 1.0571631464417852, 1.1568701753503277,
1.1544146029154017, 1.1367586548474078, 1.1375927085355826,
1.1307852105853577, 1.1116869413080799, 1.0904766573882994,
1.0750800341343474, 1.005129039001896 ]),
'geometry.volcon_low': array([1.1661400694942512]),
'geometry.volcon_up': array([1.1309221641689904])}
Objectives
{'LoD.val': array([62.09820375830306])}
Warning: Cutting back alpha due to evaluation error
Driver debug print for iter coord: rank0:pyOptSparse_IPOPT|3
Design Vars
{'dvs.aoa': array([2.0032014783752934]),
'dvs.decalage': array([0.]),
'dvs.shape_low': array([ 0. , 0. , 0. ,
0. , 0.0296154602965624, 0.0296154602965624,
0.0544948654902897, 0.0544948654902897, 0.0180408638416255,
0.0180408638416255, 0.0306797673353654, 0.0306797673353654,
0.0115870413752881, 0.0115870413752881, 0.05940000495 ,
0.05940000495 , 0.0047244618455004, 0.0047244618455004,
0.0170303155890019, 0.0170303155890019, 0.0027501977666769,
0.0027501977666769, 0.0240610791044129, 0.0240610791044129,
-0.0007018375716478, -0.0007018375716478, 0.0162472213549413,
0.0162472213549413, -0.0042670073458963, -0.0042670073458963,
0.0087357564897525, 0.0087357564897525, -0.0034237352595561,
-0.0034237352595561, 0.0025977149619691, 0.0025977149619691,
0.0163064325271993, 0.0163064325271993, 0.0217048370679876,
0.0217048370679876, 0. , 0. ,
0. , 0. ]),
'dvs.shape_up': array([ 0. , 0. , 0. ,
0. , 0.0340542214413093, 0.0340542214413093,
0.1067598678063311, 0.1067598678063311, 0.0291565090319743,
0.0291565090319743, 0.0749818329616147, 0.0749818329616147,
0.0201713739081493, 0.0201713739081493, 0.0431820906405099,
0.0431820906405099, -0.0024430929437868, -0.0024430929437868,
0.0315488965018511, 0.0315488965018511, -0.0146732360612831,
-0.0146732360612831, 0.0257436767337967, 0.0257436767337967,
-0.0070926688271028, -0.0070926688271028, 0.0224568299756077,
0.0224568299756077, -0.0205325076900913, -0.0205325076900913,
0.0168062942223398, 0.0168062942223398, -0.0044499742920329,
-0.0044499742920329, 0.0037121383118346, 0.0037121383118346,
-0.020068718957349 , -0.020068718957349 , 0.0244343859960092,
0.0244343859960092, -0.0037184977493955, -0.0037184977493955,
0.0098451432048147, 0.0098451432048147, -0.0183572088302094,
-0.0183572088302094, 0.0227836371305728, 0.0227836371305728,
-0.0072023136629016, -0.0072023136629016, 0.0057452500514091,
0.0057452500514091, -0.0175495523179415, -0.0175495523179415,
0.0156807832633603, 0.0156807832633603, -0.0113571300398373,
-0.0113571300398373, 0.0070995425125661, 0.0070995425125661,
-0.0120767003918859, -0.0120767003918859, 0.0104137921527944,
0.0104137921527944, -0.0150513281439272, -0.0150513281439272,
0.0064723506835975, 0.0064723506835975, -0.0186449864539217,
-0.0186449864539217, -0.007175470253345 , -0.007175470253345 ,
-0.0231056821808405, -0.0231056821808405, 0.0004777866971262,
0.0004777866971262, 0.0478271641020137, 0.0478271641020137,
0.0572911521262405, 0.0572911521262405, 0. ,
0. , 0. , 0. ]),
'dvs.translate': array([0., 0.])}
Setting up primal boundary conditions based on pyOptions:
Setting primal boundary conditions...
Setting U = (23.96534523216085 0.8382290279782065 0) at inlet
Setting primal boundary conditions...
Setting U = (23.96534523216085 0.8382290279782065 0) at bottomlet
Setting nut wall BC for upwing. BCType=nutLowReWallFunction
Setting nut wall BC for lowwing. BCType=nutLowReWallFunction
Setting k wall BC for upwing. BCType=fixedValue
Setting k wall BC for lowwing. BCType=fixedValue
Setting omega wall BC for upwing. BCType=omegaWallFunction
Setting omega wall BC for lowwing. BCType=omegaWallFunction
Running Primal Solver 004
Previous solution time 3000 not found and nothing deleted.
No finite volume options present
Starting time loop
Checking mesh quality for time = 0
Overall domain bounding box (-201.1019107348662 -201.0882625812153 0) (201.1556081639635 201.1502360502998 0.1000302194724538)
Mesh has 3 geometric (non-empty/wedge) directions 3{1}
Mesh has 3 solution (non-empty) directions 3{1}
Boundary openness (6.38066780420224e-21 9.223368789357613e-20 2.627129359312861e-16) OK.
***High aspect ratio cells found, Max aspect ratio: 1.118985155010599e+97, number of cells 81
<<Writing 81 cells with high aspect ratio to set highAspectRatioCells
Minimum face area = 1.332149265143189e-09. Maximum face area = 63.67598684828134. Face area magnitudes OK.
***Zero or negative cell volume detected. Minimum negative volume: -7.532366570454725e-07, Number of negative volume cells: 81
<<Writing 81 zero volume cells to set zeroVolumeCells
Mesh non-orthogonality Max: 177.4075208472619 average: 18.75126633577062
*Number of severely non-orthogonal (> 75 degrees) faces: 231.
***Number of non-orthogonality errors: 158.
<<Writing 389 non-orthogonal faces to set nonOrthoFaces
***Error in face pyramids: 516 faces are incorrectly oriented.
<<Writing 404 faces with incorrect orientation to set wrongOrientedFaces
***Max skewness = 2791.270908145276, 11 highly skew faces detected which may impair the quality of the results
<<Writing 13 skew faces to set skewFaces
Coupled point location match (average 0) OK.
Failed 5 mesh checks.
Nonlinear constraints
{'geometry.linearcon_low': array([0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0.,
0., 0., 0., 0., 0.]),
'geometry.linearcon_up': array([0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0.,
0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0.,
0., 0., 0., 0., 0., 0., 0., 0.]),
'geometry.thickcon_low': array([1.0330557999377392, 1.0934616438376736, 1.119242820421181 ,
1.0944806845970905, 1.0710047890241605, 1.06777603960721 ,
1.0492320707246217, 1.042386576340598 , 1.023609155768332 ,
0.972287529831544 , 1.0800320660550435, 1.0934616438376725,
1.1192428204211806, 1.0944806845970902, 1.0710047890241594,
1.0677760396072071, 1.049232070724622 , 1.042386576340597 ,
1.0236091557683282, 0.972287529831544 ]),
'geometry.thickcon_up': array([1.028581573220893 , 1.078435087675164 , 1.0772073014577006,
1.068379327423704 , 1.0687963542677914, 1.0653926052926788,
1.0558434706540403, 1.0452383286941496, 1.0375400170671738,
1.0025645195009578, 1.0285815732208923, 1.0784350876751643,
1.0772073014577006, 1.0683793274237041, 1.0687963542677912,
1.0653926052926788, 1.0558434706540398, 1.0452383286941496,
1.0375400170671736, 1.0025645195009445]),
'geometry.volcon_low': array([1.0830700347472004]),
'geometry.volcon_up': array([1.0654610820845292])}
Objectives
{'LoD.val': array([62.09820375830306])}
Warning: Cutting back alpha due to evaluation error
Driver debug print for iter coord: rank0:pyOptSparse_IPOPT|4
Design Vars
{'dvs.aoa': array([2.001600739187647]),
'dvs.decalage': array([0.]),
'dvs.shape_low': array([ 0. , 0. , 0. ,
0. , 0.0148077301482812, 0.0148077301482812,
0.0272474327451449, 0.0272474327451449, 0.0090204319208128,
0.0090204319208128, 0.0153398836676827, 0.0153398836676827,
0.0057935206876441, 0.0057935206876441, 0.029700002475 ,
0.029700002475 , 0.0023622309227502, 0.0023622309227502,
0.008515157794501 , 0.008515157794501 , 0.0013750988833385,
0.0013750988833385, 0.0120305395522064, 0.0120305395522064,
-0.0003509187858239, -0.0003509187858239, 0.0081236106774706,
0.0081236106774706, -0.0021335036729482, -0.0021335036729482,
0.0043678782448763, 0.0043678782448763, -0.001711867629778 ,
-0.001711867629778 , 0.0012988574809845, 0.0012988574809845,
0.0081532162635996, 0.0081532162635996, 0.0108524185339938,
0.0108524185339938, 0. , 0. ,
0. , 0. ]),
'dvs.shape_up': array([ 0. , 0. , 0. ,
0. , 0.0170271107206546, 0.0170271107206546,
0.0533799339031655, 0.0533799339031655, 0.0145782545159872,
0.0145782545159872, 0.0374909164808074, 0.0374909164808074,
0.0100856869540747, 0.0100856869540747, 0.021591045320255 ,
0.021591045320255 , -0.0012215464718934, -0.0012215464718934,
0.0157744482509255, 0.0157744482509255, -0.0073366180306415,
-0.0073366180306415, 0.0128718383668983, 0.0128718383668983,
-0.0035463344135514, -0.0035463344135514, 0.0112284149878039,
0.0112284149878039, -0.0102662538450456, -0.0102662538450456,
0.0084031471111699, 0.0084031471111699, -0.0022249871460165,
-0.0022249871460165, 0.0018560691559173, 0.0018560691559173,
-0.0100343594786745, -0.0100343594786745, 0.0122171929980046,
0.0122171929980046, -0.0018592488746977, -0.0018592488746977,
0.0049225716024073, 0.0049225716024073, -0.0091786044151047,
-0.0091786044151047, 0.0113918185652864, 0.0113918185652864,
-0.0036011568314508, -0.0036011568314508, 0.0028726250257046,
0.0028726250257046, -0.0087747761589707, -0.0087747761589707,
0.0078403916316802, 0.0078403916316802, -0.0056785650199187,
-0.0056785650199187, 0.003549771256283 , 0.003549771256283 ,
-0.0060383501959429, -0.0060383501959429, 0.0052068960763972,
0.0052068960763972, -0.0075256640719636, -0.0075256640719636,
0.0032361753417987, 0.0032361753417987, -0.0093224932269608,
-0.0093224932269608, -0.0035877351266725, -0.0035877351266725,
-0.0115528410904203, -0.0115528410904203, 0.0002388933485631,
0.0002388933485631, 0.0239135820510069, 0.0239135820510069,
0.0286455760631202, 0.0286455760631202, 0. ,
0. , 0. , 0. ]),
'dvs.translate': array([0., 0.])}
Setting up primal boundary conditions based on pyOptions:
Setting primal boundary conditions...
Setting U = (23.96536864139236 0.8375594797789651 0) at inlet
Setting primal boundary conditions...
Setting U = (23.96536864139236 0.8375594797789651 0) at bottomlet
Setting nut wall BC for upwing. BCType=nutLowReWallFunction
Setting nut wall BC for lowwing. BCType=nutLowReWallFunction
Setting k wall BC for upwing. BCType=fixedValue
Setting k wall BC for lowwing. BCType=fixedValue
Setting omega wall BC for upwing. BCType=omegaWallFunction
Setting omega wall BC for lowwing. BCType=omegaWallFunction
Running Primal Solver 005
Previous solution time 3000 not found and nothing deleted.
No finite volume options present
Starting time loop
Checking mesh quality for time = 0
Overall domain bounding box (-200.7582020660039 -200.7514790921915 0) (200.7808042469127 200.782317596295 0.1000179770089662)
Mesh has 3 geometric (non-empty/wedge) directions 3{1}
Mesh has 3 solution (non-empty) directions 3{1}
Boundary openness (-2.637556292590489e-20 1.239482112167082e-19 6.128800849192251e-16) OK.
***High aspect ratio cells found, Max aspect ratio: 5.29093994399472e+96, number of cells 47
<<Writing 47 cells with high aspect ratio to set highAspectRatioCells
Minimum face area = 2.225532648569384e-09. Maximum face area = 64.51757722620928. Face area magnitudes OK.
***Zero or negative cell volume detected. Minimum negative volume: -3.623958377317874e-07, Number of negative volume cells: 47
<<Writing 47 zero volume cells to set zeroVolumeCells
Mesh non-orthogonality Max: 170.9867118055276 average: 18.06504939828748
*Number of severely non-orthogonal (> 75 degrees) faces: 115.
***Number of non-orthogonality errors: 95.
<<Writing 210 non-orthogonal faces to set nonOrthoFaces
***Error in face pyramids: 282 faces are incorrectly oriented.
<<Writing 219 faces with incorrect orientation to set wrongOrientedFaces
***Max skewness = 445.7206294589248, 5 highly skew faces detected which may impair the quality of the results
<<Writing 5 skew faces to set skewFaces
Coupled point location match (average 0) OK.
Failed 5 mesh checks.
Nonlinear constraints
{'geometry.linearcon_low': array([0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0.,
0., 0., 0., 0., 0.]),
'geometry.linearcon_up': array([0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0.,
0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0.,
0., 0., 0., 0., 0., 0., 0., 0.]),
'geometry.thickcon_low': array([1.0165238562154117, 1.0467012100180748, 1.0595749200607585,
1.04721012478197 , 1.035484745644942 , 1.0338718643454723,
1.024607277471374 , 1.0211867306738616, 1.0118024857909238,
0.9861406462797405, 1.0399939003014653, 1.046701210018073 ,
1.0595749200607583, 1.0472101247819705, 1.0354847456449416,
1.0338718643454687, 1.0246072774713744, 1.0211867306738598,
1.0118024857909216, 0.9861406462797405]),
'geometry.thickcon_up': array([1.0142907866104471, 1.039217543837582 , 1.0386036507288503,
1.034189663711852 , 1.0343981771338957, 1.0326963026463394,
1.02792173532702 , 1.0226191643470746, 1.0187700085335865,
1.001282259750482 , 1.0142907866104465, 1.0392175438375826,
1.03860365072885 , 1.0341896637118517, 1.0343981771338957,
1.0326963026463396, 1.02792173532702 , 1.0226191643470743,
1.0187700085335865, 1.0012822597504687]),
'geometry.volcon_low': array([1.0415350173736753]),
'geometry.volcon_up': array([1.0327305410422984])}
Objectives
{'LoD.val': array([62.09820375830306])}
Warning: Cutting back alpha due to evaluation error
Driver debug print for iter coord: rank0:pyOptSparse_IPOPT|5
Design Vars
{'dvs.aoa': array([2.0008003695938235]),
'dvs.decalage': array([0.]),
'dvs.shape_low': array([ 0. , 0. , 0. ,
0. , 0.0074038650741406, 0.0074038650741406,
0.0136237163725724, 0.0136237163725724, 0.0045102159604064,
0.0045102159604064, 0.0076699418338413, 0.0076699418338413,
0.002896760343822 , 0.002896760343822 , 0.0148500012375 ,
0.0148500012375 , 0.0011811154613751, 0.0011811154613751,
0.0042575788972505, 0.0042575788972505, 0.0006875494416692,
0.0006875494416692, 0.0060152697761032, 0.0060152697761032,
-0.000175459392912 , -0.000175459392912 , 0.0040618053387353,
0.0040618053387353, -0.0010667518364741, -0.0010667518364741,
0.0021839391224381, 0.0021839391224381, -0.000855933814889 ,
-0.000855933814889 , 0.0006494287404923, 0.0006494287404923,
0.0040766081317998, 0.0040766081317998, 0.0054262092669969,
0.0054262092669969, 0. , 0. ,
0. , 0. ]),
'dvs.shape_up': array([ 0. , 0. , 0. ,
0. , 0.0085135553603273, 0.0085135553603273,
0.0266899669515828, 0.0266899669515828, 0.0072891272579936,
0.0072891272579936, 0.0187454582404037, 0.0187454582404037,
0.0050428434770373, 0.0050428434770373, 0.0107955226601275,
0.0107955226601275, -0.0006107732359467, -0.0006107732359467,
0.0078872241254628, 0.0078872241254628, -0.0036683090153208,
-0.0036683090153208, 0.0064359191834492, 0.0064359191834492,
-0.0017731672067757, -0.0017731672067757, 0.0056142074939019,
0.0056142074939019, -0.0051331269225228, -0.0051331269225228,
0.0042015735555849, 0.0042015735555849, -0.0011124935730082,
-0.0011124935730082, 0.0009280345779586, 0.0009280345779586,
-0.0050171797393372, -0.0050171797393372, 0.0061085964990023,
0.0061085964990023, -0.0009296244373489, -0.0009296244373489,
0.0024612858012037, 0.0024612858012037, -0.0045893022075524,
-0.0045893022075524, 0.0056959092826432, 0.0056959092826432,
-0.0018005784157254, -0.0018005784157254, 0.0014363125128523,
0.0014363125128523, -0.0043873880794854, -0.0043873880794854,
0.0039201958158401, 0.0039201958158401, -0.0028392825099593,
-0.0028392825099593, 0.0017748856281415, 0.0017748856281415,
-0.0030191750979715, -0.0030191750979715, 0.0026034480381986,
0.0026034480381986, -0.0037628320359818, -0.0037628320359818,
0.0016180876708994, 0.0016180876708994, -0.0046612466134804,
-0.0046612466134804, -0.0017938675633362, -0.0017938675633362,
-0.0057764205452101, -0.0057764205452101, 0.0001194466742816,
0.0001194466742816, 0.0119567910255034, 0.0119567910255034,
0.0143227880315601, 0.0143227880315601, 0. ,
0. , 0. , 0. ]),
'dvs.translate': array([0., 0.])}
Setting up primal boundary conditions based on pyOptions:
Setting primal boundary conditions...
Setting U = (23.96538033899339 0.8372247054341552 0) at inlet
Setting primal boundary conditions...
Setting U = (23.96538033899339 0.8372247054341552 0) at bottomlet
Setting nut wall BC for upwing. BCType=nutLowReWallFunction
Setting nut wall BC for lowwing. BCType=nutLowReWallFunction
Setting k wall BC for upwing. BCType=fixedValue
Setting k wall BC for lowwing. BCType=fixedValue
Setting omega wall BC for upwing. BCType=omegaWallFunction
Setting omega wall BC for lowwing. BCType=omegaWallFunction
Running Primal Solver 006
Previous solution time 3000 not found and nothing deleted.
No finite volume options present
Starting time loop
Checking mesh quality for time = 0
Overall domain bounding box (-200.3749847484872 -200.3717997489057 0) (200.3849392351411 200.3867704049215 0.1000102009996678)
Mesh has 3 geometric (non-empty/wedge) directions 3{1}
Mesh has 3 solution (non-empty) directions 3{1}
Boundary openness (-1.464365984757529e-19 -4.46035006257846e-20 5.898705284396428e-17) OK.
***High aspect ratio cells found, Max aspect ratio: 2.341505292859145e+96, number of cells 15
<<Writing 15 cells with high aspect ratio to set highAspectRatioCells
Minimum face area = 7.914312637348374e-11. Maximum face area = 64.77797616223728. Face area magnitudes OK.
***Zero or negative cell volume detected. Minimum negative volume: -1.558246623838535e-07, Number of negative volume cells: 15
<<Writing 15 zero volume cells to set zeroVolumeCells
Mesh non-orthogonality Max: 179.1976726779579 average: 17.78991055052065
*Number of severely non-orthogonal (> 75 degrees) faces: 82.
***Number of non-orthogonality errors: 34.
<<Writing 116 non-orthogonal faces to set nonOrthoFaces
***Error in face pyramids: 111 faces are incorrectly oriented.
<<Writing 89 faces with incorrect orientation to set wrongOrientedFaces
***Max skewness = 234.0707992292716, 9 highly skew faces detected which may impair the quality of the results
<<Writing 9 skew faces to set skewFaces
Coupled point location match (average 0) OK.
Failed 5 mesh checks.
Nonlinear constraints
{'geometry.linearcon_low': array([0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0.,
0., 0., 0., 0., 0.]),
'geometry.linearcon_up': array([0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0.,
0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0.,
0., 0., 0., 0., 0., 0., 0., 0.]),
'geometry.thickcon_low': array([1.008260892313176 , 1.0233426950303477, 1.0297748283160137,
1.0235969851288578, 1.0177377297631456, 1.0169316914468642,
1.012301369408262 , 1.0105916743974275, 1.0059007106657818,
0.9930695597995646, 1.0199910913997194, 1.0233426950303464,
1.0297748283160135, 1.023596985128858 , 1.0177377297631447,
1.0169316914468607, 1.0123013694082625, 1.0105916743974257,
1.0059007106657762, 0.9930695597995852]),
'geometry.thickcon_up': array([1.0071453933052235, 1.0196087719187912, 1.019301825364425 ,
1.0170948318559259, 1.017199088566948 , 1.0163481513231696,
1.0139608676635103, 1.0113095821735372, 1.009385004266793 ,
1.0006411298752407, 1.0071453933052228, 1.0196087719187914,
1.0193018253644248, 1.017094831855926 , 1.017199088566948 ,
1.0163481513231698, 1.0139608676635103, 1.0113095821735372,
1.0093850042667931, 1.0006411298752296]),
'geometry.volcon_low': array([1.0207675086869126]),
'geometry.volcon_up': array([1.0163652705211832])}
Objectives
{'LoD.val': array([62.09820375830306])}
Warning: Cutting back alpha due to evaluation error
Driver debug print for iter coord: rank0:pyOptSparse_IPOPT|6
Design Vars
{'dvs.aoa': array([2.0004001847969115]),
'dvs.decalage': array([0.]),
'dvs.shape_low': array([ 0. , 0. , 0. ,
0. , 0.0037019325370703, 0.0037019325370703,
0.0068118581862862, 0.0068118581862862, 0.0022551079802032,
0.0022551079802032, 0.0038349709169207, 0.0038349709169207,
0.001448380171911 , 0.001448380171911 , 0.00742500061875 ,
0.00742500061875 , 0.0005905577306876, 0.0005905577306876,
0.0021287894486252, 0.0021287894486252, 0.0003437747208346,
0.0003437747208346, 0.0030076348880516, 0.0030076348880516,
-0.000087729696456 , -0.000087729696456 , 0.0020309026693677,
0.0020309026693677, -0.000533375918237 , -0.000533375918237 ,
0.0010919695612191, 0.0010919695612191, -0.0004279669074445,
-0.0004279669074445, 0.0003247143702461, 0.0003247143702461,
0.0020383040658999, 0.0020383040658999, 0.0027131046334985,
0.0027131046334985, 0. , 0. ,
0. , 0. ]),
'dvs.shape_up': array([ 0. , 0. , 0. ,
0. , 0.0042567776801637, 0.0042567776801637,
0.0133449834757914, 0.0133449834757914, 0.0036445636289968,
0.0036445636289968, 0.0093727291202018, 0.0093727291202018,
0.0025214217385187, 0.0025214217385187, 0.0053977613300637,
0.0053977613300637, -0.0003053866179734, -0.0003053866179734,
0.0039436120627314, 0.0039436120627314, -0.0018341545076604,
-0.0018341545076604, 0.0032179595917246, 0.0032179595917246,
-0.0008865836033878, -0.0008865836033878, 0.002807103746951 ,
0.002807103746951 , -0.0025665634612614, -0.0025665634612614,
0.0021007867777925, 0.0021007867777925, -0.0005562467865041,
-0.0005562467865041, 0.0004640172889793, 0.0004640172889793,
-0.0025085898696686, -0.0025085898696686, 0.0030542982495011,
0.0030542982495011, -0.0004648122186744, -0.0004648122186744,
0.0012306429006018, 0.0012306429006018, -0.0022946511037762,
-0.0022946511037762, 0.0028479546413216, 0.0028479546413216,
-0.0009002892078627, -0.0009002892078627, 0.0007181562564261,
0.0007181562564261, -0.0021936940397427, -0.0021936940397427,
0.00196009790792 , 0.00196009790792 , -0.0014196412549797,
-0.0014196412549797, 0.0008874428140708, 0.0008874428140708,
-0.0015095875489857, -0.0015095875489857, 0.0013017240190993,
0.0013017240190993, -0.0018814160179909, -0.0018814160179909,
0.0008090438354497, 0.0008090438354497, -0.0023306233067402,
-0.0023306233067402, -0.0008969337816681, -0.0008969337816681,
-0.0028882102726051, -0.0028882102726051, 0.0000597233371408,
0.0000597233371408, 0.0059783955127517, 0.0059783955127517,
0.0071613940157801, 0.0071613940157801, 0. ,
0. , 0. , 0. ]),
'dvs.translate': array([0., 0.])}
Setting up primal boundary conditions based on pyOptions:
Setting primal boundary conditions...
Setting U = (23.96538618604021 0.8370573182004816 0) at inlet
Setting primal boundary conditions...
Setting U = (23.96538618604021 0.8370573182004816 0) at bottomlet
Setting nut wall BC for upwing. BCType=nutLowReWallFunction
Setting nut wall BC for lowwing. BCType=nutLowReWallFunction
Setting k wall BC for upwing. BCType=fixedValue
Setting k wall BC for lowwing. BCType=fixedValue
Setting omega wall BC for upwing. BCType=omegaWallFunction
Setting omega wall BC for lowwing. BCType=omegaWallFunction
Running Primal Solver 007
Previous solution time 3000 not found and nothing deleted.
No finite volume options present
Starting time loop
Checking mesh quality for time = 0
Overall domain bounding box (-200.1783841258557 -200.1768379846372 0) (200.1826909633366 200.1838894513989 0.100005391570451)
Mesh has 3 geometric (non-empty/wedge) directions 3{1}
Mesh has 3 solution (non-empty) directions 3{1}
Boundary openness (-1.10355134116653e-19 1.81109991986004e-20 -1.365880582418776e-16) OK.
***High aspect ratio cells found, Max aspect ratio: 1.474485647289277e+96, number of cells 7
<<Writing 7 cells with high aspect ratio to set highAspectRatioCells
Minimum face area = 3.327111080994577e-08. Maximum face area = 64.88133911119255. Face area magnitudes OK.
***Zero or negative cell volume detected. Minimum negative volume: -5.077337642218166e-08, Number of negative volume cells: 7
<<Writing 7 zero volume cells to set zeroVolumeCells
Mesh non-orthogonality Max: 158.8543809331002 average: 17.63615380861624
*Number of severely non-orthogonal (> 75 degrees) faces: 34.
***Number of non-orthogonality errors: 10.
<<Writing 44 non-orthogonal faces to set nonOrthoFaces
***Error in face pyramids: 41 faces are incorrectly oriented.
<<Writing 33 faces with incorrect orientation to set wrongOrientedFaces
***Max skewness = 40.67666462481496, 3 highly skew faces detected which may impair the quality of the results
<<Writing 3 skew faces to set skewFaces
Coupled point location match (average 0) OK.
Failed 5 mesh checks.
Nonlinear constraints
{'geometry.linearcon_low': array([0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0.,
0., 0., 0., 0., 0.]),
'geometry.linearcon_up': array([0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0.,
0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0.,
0., 0., 0., 0., 0., 0., 0., 0.]),
'geometry.thickcon_low': array([1.0041301840121308, 1.0116693015617022, 1.0148841172135046,
1.0117964025894395, 1.008867673498194 , 1.0084647588369602,
1.006150106966387 , 1.0052954077822098, 1.0029502211011552,
0.9965345910552835, 1.0099940375160967, 1.0116693015617009,
1.0148841172135041, 1.01179640258944 , 1.0088676734981932,
1.0084647588369564, 1.0061501069663872, 1.0052954077822092,
1.002950221101152 , 0.9965345910552835]),
'geometry.thickcon_up': array([1.0035726966526124, 1.0098043859593957, 1.0096509126822124,
1.0085474159279628, 1.0085995442834743, 1.0081740756615847,
1.0069804338317554, 1.0056547910867684, 1.0046925021333961,
1.0003205649376234, 1.0035726966526117, 1.009804385959396 ,
1.0096509126822126, 1.0085474159279628, 1.008599544283474 ,
1.0081740756615847, 1.0069804338317554, 1.0056547910867684,
1.0046925021333966, 1.0003205649376123]),
'geometry.volcon_low': array([1.0103837543435314]),
'geometry.volcon_up': array([1.0081826352606256])}
Objectives
{'LoD.val': array([62.09820375830306])}
Warning: Cutting back alpha due to evaluation error
Driver debug print for iter coord: rank0:pyOptSparse_IPOPT|7
Design Vars
{'dvs.aoa': array([2.000200092398456]),
'dvs.decalage': array([0.]),
'dvs.shape_low': array([ 0. , 0. , 0. ,
0. , 0.0018509662685351, 0.0018509662685351,
0.0034059290931431, 0.0034059290931431, 0.0011275539901016,
0.0011275539901016, 0.0019174854584603, 0.0019174854584603,
0.0007241900859555, 0.0007241900859555, 0.003712500309375 ,
0.003712500309375 , 0.0002952788653438, 0.0002952788653438,
0.0010643947243126, 0.0010643947243126, 0.0001718873604173,
0.0001718873604173, 0.0015038174440258, 0.0015038174440258,
-0.000043864848228 , -0.000043864848228 , 0.0010154513346838,
0.0010154513346838, -0.0002666879591185, -0.0002666879591185,
0.0005459847806095, 0.0005459847806095, -0.0002139834537223,
-0.0002139834537223, 0.0001623571851231, 0.0001623571851231,
0.00101915203295 , 0.00101915203295 , 0.0013565523167492,
0.0013565523167492, 0. , 0. ,
0. , 0. ]),
'dvs.shape_up': array([ 0. , 0. , 0. ,
0. , 0.0021283888400818, 0.0021283888400818,
0.0066724917378957, 0.0066724917378957, 0.0018222818144984,
0.0018222818144984, 0.0046863645601009, 0.0046863645601009,
0.0012607108692593, 0.0012607108692593, 0.0026988806650319,
0.0026988806650319, -0.0001526933089867, -0.0001526933089867,
0.0019718060313657, 0.0019718060313657, -0.0009170772538302,
-0.0009170772538302, 0.0016089797958623, 0.0016089797958623,
-0.0004432918016939, -0.0004432918016939, 0.0014035518734755,
0.0014035518734755, -0.0012832817306307, -0.0012832817306307,
0.0010503933888962, 0.0010503933888962, -0.0002781233932521,
-0.0002781233932521, 0.0002320086444897, 0.0002320086444897,
-0.0012542949348343, -0.0012542949348343, 0.0015271491247506,
0.0015271491247506, -0.0002324061093372, -0.0002324061093372,
0.0006153214503009, 0.0006153214503009, -0.0011473255518881,
-0.0011473255518881, 0.0014239773206608, 0.0014239773206608,
-0.0004501446039313, -0.0004501446039313, 0.0003590781282131,
0.0003590781282131, -0.0010968470198713, -0.0010968470198713,
0.00098004895396 , 0.00098004895396 , -0.0007098206274898,
-0.0007098206274898, 0.0004437214070354, 0.0004437214070354,
-0.0007547937744929, -0.0007547937744929, 0.0006508620095497,
0.0006508620095497, -0.0009407080089954, -0.0009407080089954,
0.0004045219177248, 0.0004045219177248, -0.0011653116533701,
-0.0011653116533701, -0.0004484668908341, -0.0004484668908341,
-0.0014441051363025, -0.0014441051363025, 0.0000298616685704,
0.0000298616685704, 0.0029891977563759, 0.0029891977563759,
0.00358069700789 , 0.00358069700789 , 0. ,
0. , 0. , 0. ]),
'dvs.translate': array([0., 0.])}
Setting up primal boundary conditions based on pyOptions:
Setting primal boundary conditions...
Setting U = (23.96538910912521 0.8369736245683314 0) at inlet
Setting primal boundary conditions...
Setting U = (23.96538910912521 0.8369736245683314 0) at bottomlet
Setting nut wall BC for upwing. BCType=nutLowReWallFunction
Setting nut wall BC for lowwing. BCType=nutLowReWallFunction
Setting k wall BC for upwing. BCType=fixedValue
Setting k wall BC for lowwing. BCType=fixedValue
Setting omega wall BC for upwing. BCType=omegaWallFunction
Setting omega wall BC for lowwing. BCType=omegaWallFunction
Running Primal Solver 008
Previous solution time 3000 not found and nothing deleted.
No finite volume options present
Starting time loop
Checking mesh quality for time = 0
Overall domain bounding box (-200.0999120254473 -200.0991149018927 0) (200.1017493064568 200.102456242731 0.1000027644020987)
Mesh has 3 geometric (non-empty/wedge) directions 3{1}
Mesh has 3 solution (non-empty) directions 3{1}
Boundary openness (1.416144782950207e-19 -8.419486351775251e-20 3.796571340996087e-17) OK.
***High aspect ratio cells found, Max aspect ratio: 2.913138629025231e+95, number of cells 1
<<Writing 1 cells with high aspect ratio to set highAspectRatioCells
Minimum face area = 3.315151661345229e-10. Maximum face area = 64.92572085232982. Face area magnitudes OK.
***Zero or negative cell volume detected. Minimum negative volume: -3.315982040266679e-11, Number of negative volume cells: 1
<<Writing 1 zero volume cells to set zeroVolumeCells
Mesh non-orthogonality Max: 156.5005135437737 average: 17.56171051188666
*Number of severely non-orthogonal (> 75 degrees) faces: 11.
***Number of non-orthogonality errors: 2.
<<Writing 13 non-orthogonal faces to set nonOrthoFaces
***Error in face pyramids: 7 faces are incorrectly oriented.
<<Writing 7 faces with incorrect orientation to set wrongOrientedFaces
Max skewness = 3.155094101388551 OK.
Coupled point location match (average 0) OK.
Failed 4 mesh checks.
Nonlinear constraints
{'geometry.linearcon_low': array([0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0.,
0., 0., 0., 0., 0.]),
'geometry.linearcon_up': array([0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0.,
0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0.,
0., 0., 0., 0., 0., 0., 0., 0.]),
'geometry.thickcon_low': array([1.002065026064755 , 1.0058341303901865, 1.0074412161014732,
1.0058976696086375, 1.0044335349585423, 1.0042321042609839,
1.0030749077203178, 1.0026475956864522, 1.0014750768443532,
0.9982672485623284, 1.0049966360872538, 1.005834130390185 ,
1.0074412161014725, 1.0058976696086375, 1.0044335349585418,
1.0042321042609803, 1.0030749077203174, 1.002647595686452 ,
1.0014750768443483, 0.9982672485623492]),
'geometry.thickcon_up': array([1.0017863483263063, 1.004902192979698 , 1.004825456341106 ,
1.0042737079639814, 1.0042997721417373, 1.0040870378307927,
1.0034902169158777, 1.0028273955433837, 1.0023462510666972,
1.000160282468815 , 1.0017863483263056, 1.0049021929796982,
1.0048254563411063, 1.0042737079639812, 1.0042997721417373,
1.0040870378307927, 1.0034902169158781, 1.0028273955433837,
1.0023462510666978, 1.0001602824688014]),
'geometry.volcon_low': array([1.0051918771718404]),
'geometry.volcon_up': array([1.0040913176303468])}
Objectives
{'LoD.val': array([62.09820375830306])}
Warning: Cutting back alpha due to evaluation error
Driver debug print for iter coord: rank0:pyOptSparse_IPOPT|8
Design Vars
{'dvs.aoa': array([2.0001000461992278]),
'dvs.decalage': array([0.]),
'dvs.shape_low': array([ 0. , 0. , 0. ,
0. , 0.0009254831342676, 0.0009254831342676,
0.0017029645465716, 0.0017029645465716, 0.0005637769950508,
0.0005637769950508, 0.0009587427292302, 0.0009587427292302,
0.0003620950429778, 0.0003620950429778, 0.0018562501546875,
0.0018562501546875, 0.0001476394326719, 0.0001476394326719,
0.0005321973621563, 0.0005321973621563, 0.0000859436802087,
0.0000859436802087, 0.0007519087220129, 0.0007519087220129,
-0.000021932424114 , -0.000021932424114 , 0.0005077256673419,
0.0005077256673419, -0.0001333439795593, -0.0001333439795593,
0.0002729923903048, 0.0002729923903048, -0.0001069917268611,
-0.0001069917268611, 0.0000811785925615, 0.0000811785925615,
0.000509576016475 , 0.000509576016475 , 0.0006782761583746,
0.0006782761583746, 0. , 0. ,
0. , 0. ]),
'dvs.shape_up': array([ 0. , 0. , 0. ,
0. , 0.0010641944200409, 0.0010641944200409,
0.0033362458689478, 0.0033362458689478, 0.0009111409072492,
0.0009111409072492, 0.0023431822800505, 0.0023431822800505,
0.0006303554346297, 0.0006303554346297, 0.0013494403325159,
0.0013494403325159, -0.0000763466544933, -0.0000763466544933,
0.0009859030156828, 0.0009859030156828, -0.0004585386269151,
-0.0004585386269151, 0.0008044898979311, 0.0008044898979311,
-0.000221645900847 , -0.000221645900847 , 0.0007017759367377,
0.0007017759367377, -0.0006416408653154, -0.0006416408653154,
0.0005251966944481, 0.0005251966944481, -0.000139061696626 ,
-0.000139061696626 , 0.0001160043222448, 0.0001160043222448,
-0.0006271474674172, -0.0006271474674172, 0.0007635745623753,
0.0007635745623753, -0.0001162030546686, -0.0001162030546686,
0.0003076607251505, 0.0003076607251505, -0.000573662775944 ,
-0.000573662775944 , 0.0007119886603304, 0.0007119886603304,
-0.0002250723019657, -0.0002250723019657, 0.0001795390641065,
0.0001795390641065, -0.0005484235099357, -0.0005484235099357,
0.00049002447698 , 0.00049002447698 , -0.0003549103137449,
-0.0003549103137449, 0.0002218607035177, 0.0002218607035177,
-0.0003773968872464, -0.0003773968872464, 0.0003254310047748,
0.0003254310047748, -0.0004703540044977, -0.0004703540044977,
0.0002022609588624, 0.0002022609588624, -0.0005826558266851,
-0.0005826558266851, -0.000224233445417 , -0.000224233445417 ,
-0.0007220525681513, -0.0007220525681513, 0.0000149308342852,
0.0000149308342852, 0.0014945988781879, 0.0014945988781879,
0.001790348503945 , 0.001790348503945 , 0. ,
0. , 0. , 0. ]),
'dvs.translate': array([0., 0.])}
Setting up primal boundary conditions based on pyOptions:
Setting primal boundary conditions...
Setting U = (23.9653905705581 0.836931777748428 0) at inlet
Setting primal boundary conditions...
Setting U = (23.9653905705581 0.836931777748428 0) at bottomlet
Setting nut wall BC for upwing. BCType=nutLowReWallFunction
Setting nut wall BC for lowwing. BCType=nutLowReWallFunction
Setting k wall BC for upwing. BCType=fixedValue
Setting k wall BC for lowwing. BCType=fixedValue
Setting omega wall BC for upwing. BCType=omegaWallFunction
Setting omega wall BC for lowwing. BCType=omegaWallFunction
Running Primal Solver 009
Previous solution time 3000 not found and nothing deleted.
No finite volume options present
Starting time loop
Checking mesh quality for time = 0
Overall domain bounding box (-200.05708988685 -200.0566753764777 0) (200.0579167557994 200.0583070359965 0.1000013987329265)
Mesh has 3 geometric (non-empty/wedge) directions 3{1}
Mesh has 3 solution (non-empty) directions 3{1}
Boundary openness (1.535515940087082e-19 -1.128681359749453e-20 -5.577546879439494e-16) OK.
Max cell openness = 5.687698894642443e-15 OK.
Max aspect ratio = 845.6910038781288 OK.
Minimum face area = 1.091587689619168e-07. Maximum face area = 64.93997465563967. Face area magnitudes OK.
Min volume = 1.095704285682162e-08. Max volume = 6.493996705161612. Total volume = 15998.65793360898. Cell volumes OK.
Mesh non-orthogonality Max: 72.01880967127045 average: 17.530355500006
Non-orthogonality check OK.
Face pyramids OK.
Max skewness = 2.73244002952865 OK.
Coupled point location match (average 0) OK.
Mesh OK.
Beta Was this translation helpful? Give feedback.
All reactions