PeleC Channel Flow

[fahadrhisat]

I was trying to simulate the channel flow using PeleC. There are a few things I want to modify,

1. Rather than using a pressure gradient (dpdx), which is given as default I want to use a fixed value for the flow velocity.
2. The default code computes u, v, and w. I want to include a few more derived variables such as wall shear stress, skin friction coefficient, turbulent kinetic energy, and Reynolds Stress components.
3. Is there any way to estimate the near wall resolution (y+) by modifying any part of the code?

[baperry2]

1. Running a periodic channel, the average velocity at statistically stationary state will be a function of dpdx. You can adjust dpdx to achieve the average velocity you want. If you want to run a spatially developing channel case, you can set inflow/outflow boundaries and use the bcnormal function to specify your desired velocity at the inflow.
2. PeleC supports user-defined derived variables, defined as part of the ProblemSpecificFunctions struct. See prob.H in Exec/Regtests/EB-C9 for an example. You can code up anything that's a function of the the state in a cell or adjacent cells.
3. We don't have a built in way to compute y+. You should be able to post-process it from information in plt files though.