I found that a 3D model of viscous flow in a rectangular channel with jump section implemented on the basis of FEM demonstrates a transition to turbulence. We use an algorithm for unsteady flows described on

https://community.wolfram.com/groups/-/m/t/1433064 Solver for unsteady flow with the use of the Navier-Stokes and Mathematica FEM

Effectively this algorithm works as LES. At one step, the average like values {UX[i], VY[i], WZ[i], P0[i]} are calculated here, and at another step it is used to calculate {u, v, w, p}. Here is the code and visualization of the flow:

H = 1/2; L = 2;
\[CapitalOmega] = 
   0 <= z <= H && 0 <= x <= L && 
    0 <= y <= 
     H && ! (x >= 1 && y <= 0.1) && ! (x >= 1 && 
       y >= 0.4) && ! (x >= 1 && z <= 0.1) && ! (x >= 1 && 
       z >= 0.4), {x, y, z}];



Figure 1

We see that this is similar to turbulence, although we did not use any model of turbulence. There is a question about the algorithm for calculating the interpolation function in 3D in the case of NDSolveValue[] with FEM. How is the interpolation function calculated in NDSolveValue[] with 3D FEM? Is it possible to make this algorithm faster?

Update 1. This algorithm can be used up to Reynolds number $Re=4500$ (it corresponds to $\nu =5\times 10^{-5}$ in the code above). Flow visualisation in the cross-section $x=0.95$ with k=22

frames = Table[
   StreamDensityPlot[{VY[i][.95, y, z], WZ[i][.95, y, z]}, {y, 0, 
     H}, {z, 0, H}, ColorFunction -> "Rainbow", Frame -> False], {i, 
    1, 22}];

Figure 2

  • $\begingroup$ I don't understand: This looks like you are solving the steady flow equations. Where does time come in? There is no "turbulence" in steady flows. $\endgroup$
    – Pirx
    Mar 27 '20 at 22:22
  • 1
    $\begingroup$ Do you see type (u[x,y,z]-UX[i-1][x,y,z])/t0 terms there? This is approximation of $\frac {\partial u}{\partial t}$ in a numerical scheme that I developed for transient flows. $\endgroup$ Mar 27 '20 at 23:33
  • 1
    $\begingroup$ O.k., got you. We're looking at one long piece of "spaghetti code", so I missed that. $\endgroup$
    – Pirx
    Mar 27 '20 at 23:39
  • 1
    $\begingroup$ This does not work in another layout. $\endgroup$ Mar 27 '20 at 23:42

Your Answer

By clicking “Post Your Answer”, you agree to our terms of service, privacy policy and cookie policy

Browse other questions tagged or ask your own question.