added dsquare_time_averaged.m for figure 13

data/code/3.4_sim_distance_velocity_vectors/dQuadratZeitlichGemittelt.m

+% ** coded by Christian Weiland **
+% ** christian.weiland@tuhh.de **
+%
+% This script "dsquare_timeaveraged.m" creates data for figure 13.
+% It calculates the squared distance between ideal and inertial particles averaged
+% over all particles at each time point. 
+%
+% ADDITIONAL INFORMATION
+% 1. The input files must conform to the standard output structure of the
+% particle data by the option "Moving Probes" of M-Star CFD.
+
+
+clear variables
+close all
+tic
+
+prompt = 'Versatz festgelegter Startpositionen (1) oder Winkel über Radius und Höhe (0)?';
+ANTWORT = input(prompt);
+
+% number of particles in resp. direction, which are set in M-Star
+numX = 5;                               
+numY = 42;
+numZ = 25;
+
+
+R = 0.065;
+H = 0.25;
+
+disp('Lade Daten ein...')
+dataidIdeal = importdata('probeIdeal.txt');
+dataidKL    = importdata('probeKL.txt');
+dataidGS    = importdata('probeGS.txt');
+
+disp('... fertig mit dem Einladen!')
+
+ideal       = dataidIdeal.data;
+KL          = dataidKL.data;
+GS          = dataidGS.data;
+
+%Partikel-ID startet bei 0, soll bei 1 starten
+ideal(:,1)  = ideal(:,1) + 1;
+KL(:,1)     = KL(:,1) + 1;
+GS(:,1)     = GS(:,1) + 1;
+
+%Sortieren nach ID und Zeit
+ideal       = sortrows(ideal);
+KL          = sortrows(KL);
+GS          = sortrows(GS);
+
+%Programm legt (0,0,0) an den äußersten Punkt, soll aber unten auf der Welle liegen
+% -> lineare Verschiebung um halben Durchmesser in x- und z-Richtung
+ideal(:,3)  = ideal(:,3) - 0.065;
+ideal(:,5)  = ideal(:,5) - 0.065;
+KL(:,3)     = KL(:,3) - 0.065;
+KL(:,5)     = KL(:,5) - 0.065;
+GS(:,3)     = GS(:,3) - 0.065;
+GS(:,5)     = GS(:,5) - 0.065;
+
+%Berechnung des Radiusses der Partikel
+rIdeal      = sqrt(ideal(:,3).^2 + ideal(:,5).^2);
+rKL         = sqrt(KL(:,3).^2 + KL(:,5).^2);
+rGS         = sqrt(GS(:,3).^2 + GS(:,5).^2);
+
+%Radius ist der 14. Eintrag
+ideal       = [ideal rIdeal];
+KL          = [KL rKL];
+GS          = [GS rGS];
+
+clear data*
+
+schritte     = find(ideal(:,1)==1, 1, 'last');
+numPart      = max(ideal(:,1));
+zeit = ideal(1:schritte,2)-ideal(1,2);
+
+dXKL          = zeros(numPart,schritte);
+dXGS          = zeros(numPart,schritte);
+
+dYKL          = zeros(numPart,schritte);
+dYGS          = zeros(numPart,schritte);
+
+dZKL          = zeros(numPart,schritte);
+dZGS          = zeros(numPart,schritte);
+
+dKL           = zeros(numPart,schritte);
+dGS           = zeros(numPart,schritte);
+
+for m = 1:numPart
+    for n = 1:schritte
+        
+        dXKL(m,n)  = KL(n + (m-1)*schritte, 3) - ideal(n + (m-1)*schritte, 3);
+        dXGS(m,n)  = GS(n + (m-1)*schritte, 3) - ideal(n + (m-1)*schritte, 3);
+        
+        dYKL(m,n)  = KL(n + (m-1)*schritte, 4) - ideal(n + (m-1)*schritte, 4);
+        dYGS(m,n)  = GS(n + (m-1)*schritte, 4) - ideal(n + (m-1)*schritte, 4);
+         
+        dZKL(m,n)  = KL(n + (m-1)*schritte, 5) - ideal(n + (m-1)*schritte, 5);
+        dZGS(m,n)  = GS(n + (m-1)*schritte, 5) - ideal(n + (m-1)*schritte, 5);
+          
+        dKL(m,n)   = norm([dXKL(m,n) dYKL(m,n) dZKL(m,n)]',2);
+        dGS(m,n)   = norm([dXGS(m,n) dYGS(m,n) dZGS(m,n)]',2);
+    end    
+end
+
+meanGS = mean(dGS.*dGS);
+meanKL = mean(dKL.*dKL);
\ No newline at end of file