zhyang-dev/CFDEMcoupling-PFM
1
0
Fork 0
mirror of https://github.com/ParticulateFlow/CFDEMcoupling-PFM.git synced 2025-12-08 06:37:44 +00:00
Code Issues Actions Packages Projects Releases Wiki Activity

change data plotting of falling_sphere_two_way_coupling case

Browse Source 
simply plot particle's omega-z and pos-y of actual case instead of
comparing omega-z of actual case to a case without torque coupling
(requiring a 2nd simulation with modified setup)
This commit is contained in:
danielque
2023-01-10 12:23:29 +01:00
parent 29410326d1
commit 4283ff09e7
1 changed files with 27 additions and 21 deletions
Download Patch File Download Diff File Expand all files Collapse all files

48
tutorials/cfdemSolverIB/falling_sphere_two_way_coupling/CFD/octave/plot_data.m
View File

@ -1,33 +1,39 @@
clear all;
clc;
%% Script to plot the angular velocities and the particle positions
% plot the angular velocity and the particle position
A = importdata('../../DEM/post/angular_velocity_no_coupling.txt',' ',1);
B = importdata('../../DEM/post/position_no_coupling.txt',' ',1);
C = importdata('../../DEM/post/angular_velocity.txt',' ',1);
D = importdata('../../DEM/post/position.txt',' ',1);
A = importdata('../../DEM/post/angular_velocity.txt',' ',1);
B = importdata('../../DEM/post/position.txt',' ',1);
pos1 = B.data();
omega1 = A.data();
pos2 = D.data();
omega2 = C.data();
omega = A.data();
pos = B.data();
time = omega1(:,1);
omegax1 = omega1(:,2);
omegay1 = omega1(:,3);
omegaz1 = omega1(:,4);
time1 = omega(:,1);
omegax = omega(:,2);
omegay = omega(:,3);
omegaz = omega(:,4);
time2 = omega2(:,1);
omegax2 = omega2(:,2);
omegay2 = omega2(:,3);
omegaz2 = omega2(:,4);
time2 = pos(:,1);
posx = pos(:,2);
posy = pos(:,3);
posz = pos(:,4);
figure
plot(time,omegaz1,'-.-',time2,omegaz2,'Linewidth',1.5)
figure(1)
plot(time1,omegaz,'-.-','Linewidth',1.5)
xlabel('Time (s)')
ylabel('Angular velocity (1/s)')
legend('One-way coupling','Two-way coupling')
legend('Two-way torque coupling')
axis([0 0.5 0 11])
set(gca,'FontSize',12)
print('angular_velocity_compare.eps')
print('angular_velocity.eps')
figure(2)
plot(time2,posy,'-.-','Linewidth',1.5)
xlabel('Time (s)')
ylabel('Y-position (m)')
legend('Two-way torque coupling')
axis([0 0.5 0 0.1])
set(gca,'FontSize',12)
print('position.eps')