Help me convert this Equation

3 次查看（过去 30 天）
shammas mohamed
shammas mohamed 2023-4-4
评论： shammas mohamed 2023-4-4
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 个评论
显示 6更早的评论
Cris LaPierre
Cris LaPierre 2023-4-4
The 400 is likely because that is what was in your original post before it was edited.
shammas mohamed
shammas mohamed 2023-4-4
yeah sorry its not as the original equation. check the updated question

请先登录，再进行评论。

请先登录，再回答此问题。

回答（3 个）

Sam Chak
Sam Chak 2023-4-4
编辑：Sam Chak 2023-4-4
Ran in:
Hi @shammas mohamed
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 个评论
shammas mohamed
shammas mohamed 2023-4-4
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.

请先登录，再进行评论。

Walter Roberson
Walter Roberson 2023-4-4
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 个评论
Dyuman Joshi
Dyuman Joshi 2023-4-4
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.

请先登录，再进行评论。

Cris LaPierre
Cris LaPierre 2023-4-4
编辑：Cris LaPierre 2023-4-4
Ran in:
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 个评论
显示 4更早的评论
Cris LaPierre
Cris LaPierre 2023-4-4
编辑：Cris LaPierre 2023-4-4
Ran in:
The exponential function becomes 0 because that's what happens when you raise e to a larger and larger negative number.
exp(0)
ans = 1
exp(-0.1)
ans = 0.9048
exp(-1)
ans = 0.3679
exp(-10)
ans = 4.5400e-05
Because the entire expression 1-exp(...) is multiplied by , as time increases, your result will approach .
If it happens to be doing that quicker or slower than expected, check your units.

请先登录，再进行评论。

类别

在 Help Center 和 File Exchange 中查找有关 Computational Fluid Dynamics (CFD) 的更多信息

标签

产品


版本

R2022b

Community Treasure Hunt

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

Start Hunting!

Translated by