# T106c Low-Pressure Turbine Cascade

**URL:** <https://pyfr.discourse.group/t/t106c-low-pressure-turbine-cascade/1072>\
**Category:** Cases\
**Tags:** config\
**Created:** [2 April 2024 13:30 UTC](https://pyfr.discourse.group/t/t106c-low-pressure-turbine-cascade/1072 "2024-04-02T13:30:01Z")\
**Posts on this page:** 9\
**Page:** 1

<div class="post-metadata">

**Author:** ![fei](https://avatars.discourse-cdn.com/v4/letter/f/59ef9b/32.png) [@fei](https://pyfr.discourse.group/u/fei)\
**Post date:** [2 April 2024 13:30 UTC](https://pyfr.discourse.group/t/t106c-low-pressure-turbine-cascade/1072/1 "2024-04-02T13:30:01Z")

</div>

Hello, I’m replicating the T106c Low-Pressure Turbine Cascade case in the paper ‘Impact of Number Representation for High-Order Implicit Large-Eddy Simulations’. The mesh I used is in [https://how4.cenaero.be/](https://how4.cenaero.be/), the order of mesh is third. But the error comes out: NaNs. Then, I used the anti-alias, the error is AttributeError: ‘NavierStokesElements’ object has no attribute ‘\_grad\_upts’. How can I solve it? Thanks a lot.

The ini file and mesh are attached.  
[https://drive.google.com/drive/folders/1SKZb6hfqh0wLHkcVWptz2KshrginuwXN?usp=sharing](https://drive.google.com/drive/folders/1SKZb6hfqh0wLHkcVWptz2KshrginuwXN?usp=sharing)

With best regards

---

<div class="post-metadata">

**Author:** ![fdw](https://avatars.discourse-cdn.com/v4/letter/f/a587f6/32.png) [@fdw](https://pyfr.discourse.group/u/fdw)\
**Post date:** [3 April 2024 13:06 UTC](https://pyfr.discourse.group/t/t106c-low-pressure-turbine-cascade/1072/2 "2024-04-03T13:06:52Z")

</div>

This could be a bug. Let me get back to you in a couple of days.

Regards, Freddie.

---

<div class="post-metadata">

**Author:** ![fei](https://avatars.discourse-cdn.com/v4/letter/f/59ef9b/32.png) [@fei](https://pyfr.discourse.group/u/fei)\
**Post date:** [3 April 2024 13:20 UTC](https://pyfr.discourse.group/t/t106c-low-pressure-turbine-cascade/1072/3 "2024-04-03T13:20:04Z")

</div>

Thanks for your reply. I simulated this case with anti-alias by pyfr-1.15, the error ‘NaNs’ also comes out. Could you help me check my ini file or mesh?

Thanks a lot.

---

<div class="post-metadata">

**Author:** ![fdw](https://avatars.discourse-cdn.com/v4/letter/f/a587f6/32.png) [@fdw](https://pyfr.discourse.group/u/fdw)\
**Post date:** [3 April 2024 15:59 UTC](https://pyfr.discourse.group/t/t106c-low-pressure-turbine-cascade/1072/4 "2024-04-03T15:59:50Z")

</div>

A large part of running these cases is getting the mesh, start-up procedure, and boundary conditions correct. I would therefore suggest you work to get the simulation stable at p = 1 and go from there.

Regards, Freddie.

---

<div class="post-metadata">

**Author:** ![fei](https://avatars.discourse-cdn.com/v4/letter/f/59ef9b/32.png) [@fei](https://pyfr.discourse.group/u/fei)\
**Post date:** [4 April 2024 15:30 UTC](https://pyfr.discourse.group/t/t106c-low-pressure-turbine-cascade/1072/5 "2024-04-04T15:30:39Z")

</div>

Thaks a lot. I have start at p = 1, but the simulation also breaks up. The span direction is periodic, and the up and bottom of the computational field are also periodic. The boundary condition and initial condition are presented in ini file, as the [https://acdl.mit.edu/HOW5/WorkshopPresentations/CS2\_T106\_LPT\_Cascades/CS2\_0\_T106\_LPT\_summary.pdf](https://acdl.mit.edu/HOW5/WorkshopPresentations/CS2_T106_LPT_Cascades/CS2_0_T106_LPT_summary.pdf) said. I don’t know which part is wrong.

---

<div class="post-metadata">

**Author:** ![fdw](https://avatars.discourse-cdn.com/v4/letter/f/a587f6/32.png) [@fdw](https://pyfr.discourse.group/u/fdw)\
**Post date:** [4 April 2024 18:58 UTC](https://pyfr.discourse.group/t/t106c-low-pressure-turbine-cascade/1072/6 "2024-04-04T18:58:02Z")

</div>

You will first want to check that the boundary conditions are being enforced correctly. Next, you will want to look at where exactly in the domain it is starting to diverge. Then, in response to this you may want to either increase mesh resolution in this area or adapt your start-up procedure.

Regards, Freddie.

---

<div class="post-metadata">

**Author:** ![fei](https://avatars.discourse-cdn.com/v4/letter/f/59ef9b/32.png) [@fei](https://pyfr.discourse.group/u/fei)\
**Post date:** [7 April 2024 15:09 UTC](https://pyfr.discourse.group/t/t106c-low-pressure-turbine-cascade/1072/7 "2024-04-07T15:09:50Z")

</div>

> [@fdw](#):
>
> you will want to look at where exactly in the domain it is starting to diverge

How to know where exactly in the domain it is starting to diverge?

By the way, what’s the aspect ratio in your mesh? Is there a specific numerical limit to the aspect ratio of the grid in the FR simulation?

---

<div class="post-metadata">

**Author:** ![fdw](https://avatars.discourse-cdn.com/v4/letter/f/a587f6/32.png) [@fdw](https://pyfr.discourse.group/u/fdw)\
**Post date:** [7 April 2024 17:36 UTC](https://pyfr.discourse.group/t/t106c-low-pressure-turbine-cascade/1072/8 "2024-04-07T17:36:43Z")

</div>

You will need to look at the solution as it develops and see if there are any regions of low/high pressure and or velocity.

With regards to the aspect ratio this is something you can determine yourself based off of the information given in the paper. Generally, for explicit FR an aspect ratio \< 20 is suggested.

Regards, Freddie.

---

<div class="post-metadata">

**Author:** ![fei](https://avatars.discourse-cdn.com/v4/letter/f/59ef9b/32.png) [@fei](https://pyfr.discourse.group/u/fei)\
**Post date:** [21 May 2024 08:52 UTC](https://pyfr.discourse.group/t/t106c-low-pressure-turbine-cascade/1072/9 "2024-05-21T08:52:34Z")

</div>

The isentropic Mach number is along the chordwise distribution. How to get the isentropic Mach numberalong the chordwise distribution? I know how to get the isentropic Mach number along x direction, but in the chordwise distribution, could tell me how to get it?
