search

Exawind / hfm-fsi

Input files and data for HFM FSI Work
0 stars 0 forks source link

(FY23 Q3) Run IEA 15-MW mesh in rigid mode #30

Open gantech opened 9 months ago

gantech commented 9 months ago 
Realm::initialize() End 
Filter Size 0.00733094
 1/1         myNDTW
        WallDistEQS           21         0.0144574       0.565271             1
Realm shall provide output files at : currentTime/timeStepCount: 0/0 (realm_1)
Realm shall provide restart files at: currentTime/timeStepCount: 0/0 (realm_1)
*******************************************************
Time Step Count: 1 Current Time: 0.00733094
 dtN: 0.00733094 dtNm1: 0.00733094 gammas: 1 -1 0
 DualNodalVolume min: 2.3919e-12 max: 9.07357 total: 2.0783e+06
 Max Courant: 46064.7 Max Reynolds: 5.21918e+06 (realm_1)
realm_1::advance_time_step() 
NLI    Name           Linear Iter      Linear Res     NLinear Res    Scaled NLR
---    ----           -----------      ----------     -----------    ----------
1/2    Equation System Iteration
 1/1      myLowMach
        MomentumEQS_X          8        0.00320046       0.363179             1
        MomentumEQS_Y          6        0.00384854       0.433726             1
        MomentumEQS_Z          6        0.00379047       0.433726             1
        MomentumEQS           20        0.00627869       0.712836             1
        ContinuityEQS         16          0.153119        4.90823             1
 1/1ShearStressTransportWrap
        TurbKineticEnergyEQS   4         0.0487879        8.99762             1
        SpecDissRateEQS        5           39663.7    2.06785e+06             1
2/2    Equation System Iteration
 1/1      myLowMach
        MomentumEQS_X          6          0.154319        4.76515       13.1207
        MomentumEQS_Y          7         0.0789501        3.56085       8.20991
        MomentumEQS_Z          7         0.0773811        3.56085        8.2099
        MomentumEQS           20          0.189829        6.93296       9.72588
        ContinuityEQS         14           1.28877         44.151       8.99529
 1/1ShearStressTransportWrap
        TurbKineticEnergyEQS  14         0.0450503        1.21939      0.135524
        SpecDissRateEQS        2           40152.8    3.87718e+07       18.7499
Mass Balance Review:
Density accumulation: 0
Integrated inflow:    0
Integrated open:      0
Total mass closure:   0
Filter Size 0.01466188737703458
Mean System Norm: 7754379.130866034 1 0.00733094
*******************************************************
Time Step Count: 2 Current Time: 0.0146619
 dtN: 0.00733094 dtNm1: 0.00733094 gammas: 1.5 -2 0.5
 DualNodalVolume min: 2.3919e-12 max: 9.07357 total: 2.0783e+06
 Max Courant: 24482 Max Reynolds: 5.37048e+06 (realm_1)
realm_1::advance_time_step() 
NLI    Name           Linear Iter      Linear Res     NLinear Res    Scaled NLR
---    ----           -----------      ----------     -----------    ----------
1/2    Equation System Iteration
 1/1      myLowMach
        MomentumEQS_X        100           3.28317        8.00367             1
        MomentumEQS_Y        100           4.56993        5.78416             1
        MomentumEQS_Z        100            2.8222        5.78418             1
        MomentumEQS          300           6.29509        11.4443             1
        ContinuityEQS         15          0.585038         27.127             1
 1/1ShearStressTransportWrap
        TurbKineticEnergyEQS   4         0.0299349        4.01112             1
        SpecDissRateEQS        1       4.26301e+06    1.54987e+09             1
2/2    Equation System Iteration
 1/1      myLowMach
        MomentumEQS_X          9          0.180345        4.83528      0.604132
        MomentumEQS_Y         19          0.129354        3.70795      0.641052
        MomentumEQS_Z         18          0.129844        3.70793      0.641046
        MomentumEQS           46          0.257131        7.13285      0.623266
        ContinuityEQS         16          0.276007        9.97026       0.36754
 1/1ShearStressTransportWrap
        TurbKineticEnergyEQS   4         0.0219501       0.663959       0.16553
        SpecDissRateEQS        1       3.11692e+06    5.97738e+09       3.85669
Mass Balance Review:
Density accumulation: 0
Integrated inflow:    0
Integrated open:      0
Total mass closure:   0
Filter Size 0.02199283106555188
Mean System Norm: 1195475231.010823 2 0.0146619
*******************************************************
Time Step Count: 3 Current Time: 0.0219928
 dtN: 0.00733094 dtNm1: 0.00733094 gammas: 1.5 -2 0.5
 DualNodalVolume min: 2.3919e-12 max: 9.07357 total: 2.0783e+06
 Max Courant: 14416.1 Max Reynolds: 5.57586e+06 (realm_1)
realm_1::advance_time_step() 
NLI    Name           Linear Iter      Linear Res     NLinear Res    Scaled NLR
---    ----           -----------      ----------     -----------    ----------
1/2    Equation System Iteration
 1/1      myLowMach
        MomentumEQS_X        100          0.249208       0.867707             1
        MomentumEQS_Y        100         0.0529252        0.67987             1
        MomentumEQS_Z        100          0.216324         0.6801             1
        MomentumEQS          300          0.334218        1.29525             1
        ContinuityEQS         13          0.168122         4.9971             1
 1/1ShearStressTransportWrap
        TurbKineticEnergyEQS   6         0.0527322        1.47715             1
        SpecDissRateEQS        1       2.05895e+06    6.36843e+09             1
2/2    Equation System Iteration
 1/1      myLowMach
        MomentumEQS_X        100           1.18481        1.32108       1.52249
        MomentumEQS_Y        100          0.919294        1.41266       2.07784
        MomentumEQS_Z        100            1.2369        1.40779       2.06997
        MomentumEQS          300           1.94391        2.39222       1.84692
        ContinuityEQS         13         0.0692238        2.86999       0.57433
 1/1ShearStressTransportWrap
        TurbKineticEnergyEQS  32         0.0147594       0.415858      0.281527
        SpecDissRateEQS        1            772594     1.2624e+10       1.98228
Mass Balance Review:
Density accumulation: 0
Integrated inflow:    0
Integrated open:      0
Total mass closure:   0
Filter Size 0.02932377475406917
Mean System Norm: 2524796361.623323 3 0.0219928
*******************************************************
Time Step Count: 4 Current Time: 0.0293238
 dtN: 0.00733094 dtNm1: 0.00733094 gammas: 1.5 -2 0.5
 DualNodalVolume min: 2.3919e-12 max: 9.07357 total: 2.0783e+06
 Max Courant: 10712.4 Max Reynolds: 6.04041e+06 (realm_1)
realm_1::advance_time_step() 
NLI    Name           Linear Iter      Linear Res     NLinear Res    Scaled NLR
---    ----           -----------      ----------     -----------    ----------
1/2    Equation System Iteration
 1/1      myLowMach
        MomentumEQS_X        100           1.27955        1.00262             1
        MomentumEQS_Y        100           1.21623        1.03022             1
        MomentumEQS_Z        100           1.54366        1.02991             1
        MomentumEQS          300           2.34507        1.76842             1
        ContinuityEQS         13         0.0495999        2.38697             1
 1/1ShearStressTransportWrap
        TurbKineticEnergyEQS  11         0.0255042       0.722766             1
        SpecDissRateEQS        1            622428    3.19726e+08             1
2/2    Equation System Iteration
 1/1      myLowMach
        MomentumEQS_X        100           0.47212       0.903584      0.901225
        MomentumEQS_Y        100          0.365379       0.663218      0.643761
        MomentumEQS_Z        100          0.386519       0.664835      0.645525
        MomentumEQS          300          0.711194         1.3032      0.736926
        ContinuityEQS         13          0.035942         1.5839      0.663562
 1/1ShearStressTransportWrap
        TurbKineticEnergyEQS  14        0.00957105       0.269549      0.372941
        SpecDissRateEQS        1            204289    2.27132e+08      0.710395
Mass Balance Review:
Density accumulation: 0
Integrated inflow:    0
Integrated open:      0
Total mass closure:   0
Filter Size 0.03665471844258646
Mean System Norm: 45426313.54060858 4 0.0293238
*******************************************************
Time Step Count: 5 Current Time: 0.0366547
 dtN: 0.00733094 dtNm1: 0.00733094 gammas: 1.5 -2 0.5
 DualNodalVolume min: 2.3919e-12 max: 9.07357 total: 2.0783e+06
 Max Courant: 7928.66 Max Reynolds: 6.61084e+06 (realm_1)
realm_1::advance_time_step() 
NLI    Name           Linear Iter      Linear Res     NLinear Res    Scaled NLR
---    ----           -----------      ----------     -----------    ----------
1/2    Equation System Iteration
 1/1      myLowMach
        MomentumEQS_X        100           1.34536       0.735358             1
        MomentumEQS_Y        100          0.536883       0.808004             1
        MomentumEQS_Z        100          0.595074       0.803582             1
        MomentumEQS          300             1.566        1.35623             1
        ContinuityEQS         12          0.050842        1.36611             1
 1/1ShearStressTransportWrap
        TurbKineticEnergyEQS   8         0.0159425        0.52783             1
        SpecDissRateEQS        1            364627    1.08334e+08             1
2/2    Equation System Iteration
 1/1      myLowMach
        MomentumEQS_X         40         0.0133732       0.470913      0.640387
        MomentumEQS_Y         50         0.0171178       0.461501      0.571162
        MomentumEQS_Z         41         0.0130318        0.44561      0.554529
        MomentumEQS          131         0.0253316       0.795808      0.586779
        ContinuityEQS         13         0.0186929       0.922798      0.675495
 1/1ShearStressTransportWrap
        TurbKineticEnergyEQS   5        0.00542401       0.180274      0.341538
        SpecDissRateEQS        1            143847    3.99399e+07      0.368674
Mass Balance Review:
Density accumulation: 0
Integrated inflow:    0
Integrated open:      0
Total mass closure:   0
Filter Size 0.04398566213110375
Mean System Norm: 7987984.063651924 5 0.0366547
*******************************************************
Time Step Count: 6 Current Time: 0.0439857
 dtN: 0.00733094 dtNm1: 0.00733094 gammas: 1.5 -2 0.5
 DualNodalVolume min: 2.3919e-12 max: 9.07357 total: 2.0783e+06
 Max Courant: 7702.47 Max Reynolds: 7.0898e+06 (realm_1)
realm_1::advance_time_step() 
NLI    Name           Linear Iter      Linear Res     NLinear Res    Scaled NLR
---    ----           -----------      ----------     -----------    ----------
1/2    Equation System Iteration
 1/1      myLowMach
        MomentumEQS_X         17         0.0203029       0.729528             1
        MomentumEQS_Y         17         0.0242875       0.665948             1
        MomentumEQS_Z         17          0.017488       0.657295             1
        MomentumEQS           51         0.0361652        1.18648             1
        ContinuityEQS         12         0.0285371       0.955776             1
 1/1ShearStressTransportWrap
        TurbKineticEnergyEQS   4        0.00969194        0.42986             1
        SpecDissRateEQS        1            245338    3.69397e+07             1
2/2    Equation System Iteration
 1/1      myLowMach
        MomentumEQS_X          9          0.011989        0.33666      0.461476
        MomentumEQS_Y          9         0.0083409       0.319629      0.479962
        MomentumEQS_Z          7         0.0100818       0.312948      0.476116
        MomentumEQS           25         0.0177468       0.559857      0.471864
        ContinuityEQS         13         0.0136414        0.67811      0.709486
 1/1ShearStressTransportWrap
        TurbKineticEnergyEQS   5        0.00166459        0.12665      0.294631
        SpecDissRateEQS        1            107375     5.0637e+07        1.3708
Mass Balance Review:
Density accumulation: 0
Integrated inflow:    0
Integrated open:      0
Total mass closure:   0
Filter Size 0.05131660581962105
Mean System Norm: 10127394.88214785 6 0.0439857
*******************************************************
Time Step Count: 7 Current Time: 0.0513166
 dtN: 0.00733094 dtNm1: 0.00733094 gammas: 1.5 -2 0.5
 DualNodalVolume min: 2.3919e-12 max: 9.07357 total: 2.0783e+06
 Max Courant: 7458.02 Max Reynolds: 7.41285e+06 (realm_1)
realm_1::advance_time_step() 
NLI    Name           Linear Iter      Linear Res     NLinear Res    Scaled NLR
---    ----           -----------      ----------     -----------    ----------
1/2    Equation System Iteration
 1/1      myLowMach
        MomentumEQS_X          6         0.0199112       0.683524             1
        MomentumEQS_Y          6         0.0137452       0.602462             1
        MomentumEQS_Z          6         0.0136654       0.570455             1
        MomentumEQS           18         0.0277872        1.07498             1
        ContinuityEQS         12         0.0201928       0.754423             1
 1/1ShearStressTransportWrap
        TurbKineticEnergyEQS   4        0.00701925       0.348656             1
        SpecDissRateEQS        1            177264    5.59184e+07             1
2/2    Equation System Iteration
 1/1      myLowMach
        MomentumEQS_X          7        0.00501686       0.328457      0.480535
        MomentumEQS_Y          7        0.00339508       0.245419      0.407361
        MomentumEQS_Z          7        0.00338034       0.245855      0.430981
        MomentumEQS           21        0.00693701       0.478079      0.444732
        ContinuityEQS         13         0.0188009       0.507423      0.672598
 1/1ShearStressTransportWrap
        TurbKineticEnergyEQS   5        0.00131756      0.0935769      0.268393
        SpecDissRateEQS        1           83173.6    5.90224e+07       1.05551
Mass Balance Review:
Density accumulation: 0
Integrated inflow:    0
Integrated open:      0
Total mass closure:   0
Filter Size 0.05864754950813834
Mean System Norm: 11804487.27186151 7 0.0513166
*******************************************************
Time Step Count: 8 Current Time: 0.0586475
 dtN: 0.00733094 dtNm1: 0.00733094 gammas: 1.5 -2 0.5
 DualNodalVolume min: 2.3919e-12 max: 9.07357 total: 2.0783e+06
 Max Courant: 7158.42 Max Reynolds: 7.60703e+06 (realm_1)
realm_1::advance_time_step() 
NLI    Name           Linear Iter      Linear Res     NLinear Res    Scaled NLR
---    ----           -----------      ----------     -----------    ----------
1/2    Equation System Iteration
 1/1      myLowMach
        MomentumEQS_X          6         0.0164497       0.650597             1
        MomentumEQS_Y          6         0.0105427       0.571354             1
        MomentumEQS_Z          6         0.0104713       0.541268             1
        MomentumEQS           18         0.0221673        1.02112             1
        ContinuityEQS         13           0.01007       0.484602             1
 1/1ShearStressTransportWrap
        TurbKineticEnergyEQS   4        0.00656031        0.28515             1
        SpecDissRateEQS        1            135826    5.15558e+07             1
2/2    Equation System Iteration
 1/1      myLowMach
        MomentumEQS_X          7        0.00491754       0.289844      0.445504
        MomentumEQS_Y          7        0.00337972       0.208492      0.364909
        MomentumEQS_Z          7        0.00283671       0.210151      0.388257
        MomentumEQS           21        0.00660694       0.414297      0.405726
        ContinuityEQS         13          0.012944       0.343768      0.709381
 1/1ShearStressTransportWrap
        TurbKineticEnergyEQS   5        0.00144341      0.0721788      0.253126
        SpecDissRateEQS        1           67356.5    4.08397e+07      0.792146
Mass Balance Review:
Density accumulation: 0
Integrated inflow:    0
Integrated open:      0
Total mass closure:   0
Filter Size 0.06597849319665564
Mean System Norm: 8167949.69294741 8 0.0586475
*******************************************************
Time Step Count: 9 Current Time: 0.0659785
 dtN: 0.00733094 dtNm1: 0.00733094 gammas: 1.5 -2 0.5
 DualNodalVolume min: 2.3919e-12 max: 9.07357 total: 2.0783e+06
 Max Courant: 6878.45 Max Reynolds: 7.72969e+06 (realm_1)
realm_1::advance_time_step() 
NLI    Name           Linear Iter      Linear Res     NLinear Res    Scaled NLR
---    ----           -----------      ----------     -----------    ----------
1/2    Equation System Iteration
 1/1      myLowMach
        MomentumEQS_X          6         0.0218313       0.649685             1
        MomentumEQS_Y          6         0.0140587         0.5693             1
        MomentumEQS_Z          6         0.0102757       0.548024             1
        MomentumEQS           18         0.0279257          1.023             1
        ContinuityEQS         14        0.00842321       0.363634             1
 1/1ShearStressTransportWrap
        TurbKineticEnergyEQS   4        0.00826946       0.242212             1
        SpecDissRateEQS        1            108948    2.87948e+07             1
2/2    Equation System Iteration
 1/1      myLowMach
        MomentumEQS_X         10        0.00775003       0.283809      0.436841
        MomentumEQS_Y          7        0.00607667       0.206836      0.363317
        MomentumEQS_Z         10        0.00492943       0.204505      0.373168
        MomentumEQS           27         0.0110131       0.406388      0.397252
        ContinuityEQS         14            0.0104       0.283314      0.779119
 1/1ShearStressTransportWrap
        TurbKineticEnergyEQS   6        0.00100094      0.0577367      0.238373
        SpecDissRateEQS        1           56944.6    1.84951e+07      0.642307
Mass Balance Review:
Density accumulation: 0
Integrated inflow:    0
Integrated open:      0
Total mass closure:   0
Filter Size 0.07330943688517293
Mean System Norm: 3699015.83349829 9 0.0659785
*******************************************************
Time Step Count: 10 Current Time: 0.0733094
 dtN: 0.00733094 dtNm1: 0.00733094 gammas: 1.5 -2 0.5
 DualNodalVolume min: 2.3919e-12 max: 9.07357 total: 2.0783e+06
 Max Courant: 6659.83 Max Reynolds: 7.81602e+06 (realm_1)
realm_1::advance_time_step() 
NLI    Name           Linear Iter      Linear Res     NLinear Res    Scaled NLR
---    ----           -----------      ----------     -----------    ----------
1/2    Equation System Iteration
 1/1      myLowMach
        MomentumEQS_X          6         0.0159898       0.669508             1
        MomentumEQS_Y          6         0.0130236       0.591143             1
        MomentumEQS_Z          6         0.0123115       0.578439             1
        MomentumEQS           18         0.0240179        1.06409             1
        ContinuityEQS         14        0.00643768       0.332387             1
 1/1ShearStressTransportWrap
        TurbKineticEnergyEQS   5        0.00426045       0.217077             1
        SpecDissRateEQS        1           90532.3    1.08216e+07             1
2/2    Equation System Iteration
 1/1      myLowMach
        MomentumEQS_X         31         0.0108549       0.287177      0.428937
        MomentumEQS_Y         12        0.00802065       0.214564      0.362964
        MomentumEQS_Z         12        0.00779919       0.210938      0.364668
        MomentumEQS           55          0.015588       0.415937      0.390886
        ContinuityEQS         14        0.00969405       0.283917      0.854178
 1/1ShearStressTransportWrap
        TurbKineticEnergyEQS   6        0.00155653      0.0478579      0.220465
        SpecDissRateEQS        1             49920     5.8919e+06      0.544455
Mass Balance Review:
Density accumulation: 0
Integrated inflow:    0
Integrated open:      0
Total mass closure:   0
Filter Size 0.08064038057369022
Mean System Norm: 1178379.383812134 10 0.0733094
gantech commented 9 months ago

Something's not right with the hole cut. There's more than 1 nalu-wind block that is a field in the same location of the overlap region between the blades. @ashesh2512 Could you please take a look and verify that it's ok?

Image

pscrozi commented 9 months ago

Overset has a problem with this case. Ashesh suggests we change the overset a bit on this. Ganesh is re-exporting the mesh and trying that approach out.

gantech commented 9 months ago

Changing the overset boundary to this seems to give the correct hole cut.

Image


- overset_boundary_condition: bc_overset
    overset_connectivity_type: tioga
    overset_user_data:
      tioga_populate_inactive_part: false
      mesh_group:
      - overset_name: blade1
        mesh_parts:
        - blade1-HEX
        wall_parts: [blade1]
        ovset_parts: [outerBCblade1]
      - overset_name: blade2
        mesh_parts:
        - blade2-HEX
        wall_parts: [blade2]
        ovset_parts: [outerBCblade2]
      - overset_name: blade3
        mesh_parts:
        - blade3-HEX
        wall_parts: [blade3]
        ovset_parts: [outerBCblade3]```
gantech commented 9 months ago

This one blows up after 70 time steps as well. Need to fix tip mesh region

pscrozi commented 9 months ago

Problem with the tip near the trailing edge. Poor quality mesh correlated with the tip turbulence etc. Now done to the point where there are failures on simplified examples.

Working on NM80 right now (10/9).

pscrozi commented 9 months ago

Ganesh made progress on this. Ganesh found a way to get a smooth curve on the trailing edge. He's making a beveled trailing edge. Working on the CAD and the mesh soon.

gantech commented 9 months ago

Continuing to make progress here. Integrating tip into mesh.

pscrozi commented 9 months ago

Ganesh resubmitting this run to Eagle's high priority queue.

pscrozi commented 9 months ago

Neil volunteered to help debug this issue as needed.

pscrozi commented 8 months ago

Ganesh built the mesh, Nate was building another one. The one that Nate built did run for 800 timesteps or so. Nate tried to viz it last night, but it wasn't working right. All the solves looked ok and it seemed fine. We have to submit the FSI job, maybe today. Seems that the mesh is good.