From 0d63936a1cc4c2a7f0b48d3aa34af8d90ff16564 Mon Sep 17 00:00:00 2001
From: Antoine Cyril David Hoffmann <ahoffman@spcpc606.epfl.ch>
Date: Fri, 30 Oct 2020 15:49:22 +0100
Subject: [PATCH] Notes about value found on stationary particle flux
---
wk/flux_results.m | 24 ++++++++++++++++++++++++
1 file changed, 24 insertions(+)
create mode 100644 wk/flux_results.m
diff --git a/wk/flux_results.m b/wk/flux_results.m
new file mode 100644
index 00000000..651d4b34
--- /dev/null
+++ b/wk/flux_results.m
@@ -0,0 +1,24 @@
+t0 = 100; t1 = 140; [~,it0] = min(abs(t0-Ts2D)); [~,it1] = min(abs(t1-Ts2D));
+range = it0:it1;
+avg = mean(Flux_ri(range))
+stdev = std(Flux_ri(range))^(.5)
+figure
+hist(Flux_ri(range),20)
+%%
+Gamma_ = [0.62, 0.29];
+N_ = [ 256, 128];
+L_ = [ 50, 25];
+P_ = [ 2, 2];
+J_ = [ 1, 1];
+NU_ = [0.01, 0.01];
+etaB_ = [ 0.5, 0.5];
+
+
+if 0
+%% Fig 3 of Ricci Rogers 2006
+fig = figure;
+for i = 1:numel(Gamma_)
+ semilogy(etaB_(i),Gamma_(i),'o'); hold on;
+end
+ xlabel('$\eta_B$'); ylabel('$\Gamma^\infty_{part}$')
+end
\ No newline at end of file
--
GitLab