modify examples/3D_turb_mixing

Author: Hyeoksu Lee

examples/3D_turb_mixing/case.py

@@ -31,11 +31,11 @@
 
 # Time advancement
 cfl = 0.5
-T = 200.0
+T = 150.0
 dt = cfl * dx / (u0 / M0 + 1)
 Ntfinal = int(T / dt)
 Ntstart = 0
-Nfiles = 40
+Nfiles = 30
 t_save = int(math.ceil((Ntfinal - 0) / float(Nfiles)))
 Nt = t_save * Nfiles
 t_step_start = Ntstart
@@ -74,7 +74,7 @@
             "weno_order": 5,
             "weno_eps": 1.0e-40,
             "weno_Re_flux": "F",
-            "mapped_weno": "T",
+            "wenoz": "T",
             "riemann_solver": 2,
             "wave_speeds": 1,
             "avg_state": 2,

examples/3D_turb_mixing/turbulence_stat/average_tke_over_self_similar.m

@@ -1,6 +1,8 @@
 close all; clear all;
 
+
 % Setup
+disp("Start average_tke_over_self_similar ..."); tic;
 set_user_inputs(); load variables/user_inputs.mat;
 
 ybeg = -5; yend = 5; ny = 101;
@@ -12,7 +14,7 @@
 D_averaged = zeros(ny,1);
 
 % Compute averaged TKE budget
-for q = 1:navg
+for q = 1:Nfiles
     load("results/tke_budget_data/tstep_"+string(timesteps(q))+".mat");
 
     % Normalization
@@ -40,6 +42,7 @@
 % Plot
 plot_tke_budget(T_averaged, P_averaged, D_averaged, y, mth);
 
+disp("End of program"); toc;
 
 %% FUNCTIONS
 % Compute the wall-normal derivative of a discretized function, fun(y)