I want to calculate the orientation distribution of non spherical particles in a packed bed from the coordinates x, y and z and radii of the particle.

I want to calculate the orientation distribution of non spherical particles in a packed bed from the coordinates x, y and z and radii of the particle. The results can be ploted in histogram or at the different planes X-O-Y, X-O-Z , Y-O-Z plane. The sample plots are attached along this question. I am new to coding so dont know how to proceed with this problem.

Thanks in advance

Answer 1

[Anurag]

Hi Aashna,

I understand that you want to calculate the orientation distribution of non-spherical particles in a packed bed from coordinates and radii of the particles.

It can be done in the following way:

• Load your data and calculate orientation for each particle based on its coordinates. Orientation might be represented by angles or vectors.
• Create histograms to represent the orientation distribution.
• Plot the particle orientations on different planes using scatter plot or visualization techniques.

Refer to the code snippet below to implement the steps above:

load('particle_data.mat');
theta = rand(size(x)) * 2*pi; % Replace this with your actual calculation 
% Plot histograms 
figure;

% Histogram of orientation angles
subplot(2, 2, 1);
histogram(theta, 30, 'Normalization', 'probability');
title('Orientation Distribution');

% Scatter plot on X-O-Y plane
subplot(2, 2, 2);
scatter(x, y, [], theta, 'filled');
xlabel('X');
ylabel('Y');
title('X-O-Y Plane');

% Scatter plot on X-O-Z plane
subplot(2, 2, 3);
scatter(x, z, [], theta, 'filled');
xlabel('X');
ylabel('Z');
title('X-O-Z Plane');

% Scatter plot on Y-O-Z plane
subplot(2, 2, 4);
scatter(y, z, [], theta, 'filled');
xlabel('Y');
ylabel('Z');
title('Y-O-Z Plane');

Further on plotting histograms:

https://in.mathworks.com/help/matlab/ref/matlab.graphics.chart.primitive.histogram.html

Hope this helps,

Regards,

Anurag