# Help me convert this Equation

2 visualizzazioni (ultimi 30 giorni)
shammas mohamed il 4 Apr 2023 i converted this equation before and got this figure. I forgot how but i lost the file that i wrote the code to get the figure above.
below i am trying to rewrtie it but it is not working. the exponential function is zero or infinity.
Edit! this is the original equation for the R_e.
i checked the parameters with my report and they are the same.
maybe there is a problem here too?
i did what you asked about the final denominator, still the same results.
clc
Contact_angle = 2.4783666516568;
Volume = 5;
ST_L = 72;
Density = 0.973;
Gravity = 9.807;
Shape_factor = 37.1;
T_zero = 0;
Dynamic_viscosity = 8.9e-4;
R_e = ( 4*Volume / (pi*Contact_angle) )^1/3;
T = 0:0.1:6;
R = R_e*(1 - exp(-(2*ST_L/R_e^12 + Density*Gravity/(9*R_e^10))*(24*Shape_factor*Volume^4*(T+T_zero)/(pi^2*Dynamic_viscosity)))).^(1/6);
plot(T,R);
Thank you All for the help
##### 8 CommentiMostra 7 commenti meno recentiNascondi 7 commenti meno recenti
shammas mohamed il 4 Apr 2023
yeah sorry its not as the original equation. check the updated question

Accedi per commentare.

### Risposte (3)

Sam Chak il 4 Apr 2023
Modificato: Sam Chak il 4 Apr 2023
Edit: There are two correction in the equation. Thanks @Walter Roberson for pointing out.
% Parameters
Contact_angle = 2.4783666516568;
Volume = 5;
ST_L = 72;
Density = 0.973;
Gravity = 9.807;
Shape_factor = 37.1;
T_zero = 0;
Dynamic_viscosity = 8.9e-4;
R_e = ( 4*Volume / (pi*Contact_angle) )^1/3;
% suggest to split up the terms
T = linspace(0,6e-12, 601);
f5 = 24*Shape_factor*Volume^4*(T + T_zero)/((pi^2)*Dynamic_viscosity); % <--- correction
f4 = (Density*Gravity)/(9*(R_e^10));
f3 = (2*ST_L)/(R_e^12);
f2 = f3 + f4;
f1 = - f2*f5;
R = R_e*(1 - exp(f1)).^(1/6); % <--- correction
plot(T, R), grid on, xlabel('T') ##### 1 CommentoMostra NessunoNascondi Nessuno
shammas mohamed il 4 Apr 2023
Thank you for the correction. it gave me a good start on the initial rate that the radius will start increasing exponentially from. however there is still problem that the function is not increasing. (the exponent part is 0).
I dont understand why this is happening. i edited my problem with more information if you need it.

Accedi per commentare.

Walter Roberson il 4 Apr 2023
x^1/3 means (x^1) divided by 3 which is x/3 . If you want to raise to the power of 1/3 you need x^(1/3)
Watch out for ^1/6 for the same issue.
##### 1 CommentoMostra NessunoNascondi Nessuno
Dyuman Joshi il 4 Apr 2023
Walter, I did implement this in the code in my comment, and it still doesn't give any similar result to what OP achieved earlier.

Accedi per commentare.

Cris LaPierre il 4 Apr 2023
Modificato: Cris LaPierre il 4 Apr 2023
As written, you need to either include your final denominator in parentheses, or divide by Dynamic_Viscocity.
We have to take the rest of the parameters you have given as correct. One potential error to look into - should your angles be in radians instead of degrees?
Contact_angle = 142;
Volume = 5;
ST_L = 72;
Density = 0.973;
Gravity = 9.807;
Shape_factor = 37.1;
T_zero = 0;
Dynamic_viscosity = 8.9e-4;
R_e = ( 400*Volume / (pi*Contact_angle))^1/3;
T = 0:300;
R = R_e * ( 1 - exp(-(2*ST_L/R_e^12 + (Density*Gravity) / (9*R_e^10))*24*Shape_factor*Volume^4*(T+T_zero)/(pi^2*Dynamic_viscosity))).^(1/6);
plot(T,R); ##### 6 CommentiMostra 5 commenti meno recentiNascondi 5 commenti meno recenti
shammas mohamed il 4 Apr 2023
i understand thanks

Accedi per commentare.

### Categorie

Scopri di più su Computational Fluid Dynamics (CFD) in Help Center e File Exchange

R2022b

### Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!