- 136
- 208 272
Rahmat Sunarya
Приєднався 15 вер 2010
It's all my lecture notes.
Please enjoy it.
Please enjoy it.
Friction Factor Colebrook Equation (MATLAB) using fzero syntax
#matlab #colebrook #frictionfactor #fluidmechanics
clc; clear all; close all;
%initial guess
f0=0.1;
[f, yval]=fzero(@cb,f0)
%Colebrook Equation
function y=cb(f)
NRe=6.5e4; %Reynold Number
eD = 0.00013; % e/D
y=1/sqrt(f)+0.86*log(eD/3.7+2.51/NRe/sqrt(f));
end
clc; clear all; close all;
%initial guess
f0=0.1;
[f, yval]=fzero(@cb,f0)
%Colebrook Equation
function y=cb(f)
NRe=6.5e4; %Reynold Number
eD = 0.00013; % e/D
y=1/sqrt(f)+0.86*log(eD/3.7+2.51/NRe/sqrt(f));
end
Переглядів: 48
Відео
Dimensionless Generic Breakthrough Curve
Переглядів 116Місяць тому
#matlab #adsorption #curvefitting Inspired by : Saikat Roy & Arun S. Moharir (2019). Modeling the Generic Breakthrough Curve for Adsorption Process; arxiv.org/pdf/1907.00195 clc;clear;close all; %parameter L=0.3; %m eps=0.4; %- u=0.005; %m/s LDF=0.5; %1/s LDF=15*De/Rp^2 (LDF = for simplification) b=0.1; %m3/kg qm=0.02; %kg adsd/kg adst c0=1.5; %kg/m3 rhop=1100; %kg/m3 %time-mesh t=linspace(0,50...
Ordinary Differential Equations (ODEs) Second-Order Finite Difference Solved by MATLAB
Переглядів 1412 місяці тому
28.8 Compound A diffuses through a 4-cm-long tube and reacts as it diffuses. The equation governing diffusion with reaction is: 𝐷 (𝑑^2 𝐴)/(𝑑𝑥^2 )−𝑘𝐴=0 At one end of the tube, there is a large source of A at a concentration of 0.1 M. At the other end of the tube there is an adsorbent material that quickly adsorbs any A, making the concentration 0 M. If D = 1.5 x 10-6 cm2/s and k = 5 x 10-6 1/s, ...
System Ordinary Differential Equations (ODEs) Solved by MATLAB
Переглядів 902 місяці тому
28.7 The following equations define the concentrations of three reactants: (𝑑𝑐_𝑎)/𝑑𝑡=−10𝑐_𝑎 𝑐_𝑐 𝑐_𝑏 (𝑑𝑐_𝑏)/𝑑𝑡=10𝑐_𝑎 𝑐_𝑐−𝑐_𝑏 (𝑑𝑐_𝑐)/𝑑𝑡=−10𝑐_𝑎 𝑐_𝑐 𝑐_𝑏−2𝑐_𝑐 If the initial conditions are ca = 50, cb = 0, and cc = 40, find the concentrations for the times from 0 to 3 s. #odesolver #chapra #canale #matlab
Non Linear Curve Fitting : Signal-to-Noise Ratio (SNR). Acoustic Concentration Measurement
Переглядів 943 місяці тому
#curvefitting #acoustic #matlab Paper : www.researchgate.net/publication/329773388_Applicability_of_Acoustic_Concentration_Measurements_in_Suspensions_of_Artificial_and_Natural_Sediments_Using_an_Acoustic_Doppler_Velocimeter ta in comment section Cdata=[0.001 0.002 0.002 0.004 0.005 0.010 0.020 0.042 0.107 0.202 0.206 0.399 0.600 0.855 1.161 1.699 2.078 2.751 4.569 7.401 9.524 13.574 20.902 26....
Breakthrough Adsorption using Dimensionless Number. Solved byt Method of Line
Переглядів 2403 місяці тому
The Advection PDE-ODE Equation Partial Differential Equation (PDE) Ordinary Differential Equation (ODE) #ode15s #adsorption #maths #mathematics #dimensionless #sherwoodnumber
The Simplest Way to Make an Adsorption Breakthrough Curve
Переглядів 2894 місяці тому
The Advection PDE-ODE Equation Partial Differential Equation (PDE) Ordinary Differential Equation (ODE) #ode15s #adsorption #maths #mathematics
Partial Differential Equation (PDE) Burgers - Huxley Equation | Method of Line (MATLAB)
Переглядів 3286 місяців тому
#Burgers #Huxley #waveequation #pde #matlab
Partial Differential Equation (PDE) Fix Bed Reactor with The Danckwerts boundary conditions
Переглядів 3167 місяців тому
#PDE #danckwerts #matlab #methodofline Partial Differential Equation (PDE) Fix Bed Reactor with The Danckwerts boundary conditions Method of Line Matlab Pe=200; %Peclet Number Da=2; mkohler number v=2; Nx=151; L=1; x=linspace(0,L,Nx); dx=x(2)-x(1); tf=2; Nt=101; t=linspace(0,tf,101); IC=zeros(1,Nx);
Flash Tank Calculation Real System Scatchard-Hildebrand and Second Virial
Переглядів 1657 місяців тому
#VLE #thermodynamics #flash Flash Tank Calculation Real System Scatchard-Hildebrand and Second Virial z=[0.2 0.3 0.4 0.1]; ed mole fraction t=40; lcius T=t 273.15; %Kelvin Pt=3; %atm R=82.057; %mL atm/mol/K %Antoine Eq for saturated vap pressure %Psat=10^(A-B/(C T)) A=[6.804 6.809 6.876 6.870]; B=[803.81 935.86 1075.78 1168.72]; C=[246.99 238.73 233.21 224.21]; PmmHg=10.^(A-B./(C t)); %in mmHg...
Molar Volume gas mixture using Peng-Robinson Equation of State
Переглядів 3017 місяців тому
#matlab #pengrobinson #thermodynamics #fugacity Calculate the molar volume of a gas mixture consisting of 50 mol % propane (1), 20% n-butane (2), and 30% n-pentane (3) at 40 bar and 480 K using the Peng-Robinson equation of state. Ismail Tosun (2021) , the thermodynamic of phase and reaction equilibria p.227 clc;clear;close all; R=83.14; %cm3.bar/mol/K P=40; r T=480; %K %propane(1)-nbutane(2)-...
Fugacity coefficient mixture (3 components) using Peng-Robinson Equation of state
Переглядів 4527 місяців тому
#matlab #pengrobinson #thermodynamics #fugacity R=8.314; %J/mol/K ideal gas constant T=233.2; %K system temp P=0.1; %Mpa system pressure x=[0.1 0.3 0.6]; %- fraction mol Tc=[190.56 305.32 369.83]; %Critical Temp K Pc=[4.599 4.872 4.248]; %Critical Press MPa b=[26.803 40.538 56.316]; %Peng-Robinson covolume w=[0.011 0.099 0.152]; entric factor kij=[0 0.00340 0.01070; 0.00340 0 0.00900; 0.01070 ...
Van Laar Model | Curve Fitting Vapor-Liquid Equilibrium of Binary System.
Переглядів 1527 місяців тому
Reference : Lo, Teh Cheng., Bieber, H. Harold., Karr, Andrew E. (1962). Vapor-Liquid Equilibrium of n-Pentane-Acetone. Journal of Chemical & Engineering Data 1962 7 (3), 327-331. DOI: 10.1021/je60014a002 #ConfidenceInterval #CurveFitting #SOLVER #ExcessGibbsEnergy #Redlich #Kister #thermodynamics #activitycoefficient
Redlich - Kister Expansion One Parameter | Curve Fitting Vapor-Liquid Equilibrium of Binary System.
Переглядів 4077 місяців тому
Reference : Lo, Teh Cheng., Bieber, H. Harold., Karr, Andrew E. (1962). Vapor-Liquid Equilibrium of n-Pentane-Acetone. Journal of Chemical & Engineering Data 1962 7 (3), 327-331. DOI: 10.1021/je60014a002 #ConfidenceInterval #CurveFitting #SOLVER #ExcessGibbsEnergy #Redlich #Kister #thermodynamics #activitycoefficient
IAST Selectivity Using Excel Solver Combined Langmuir-Freundlich Model
Переглядів 2537 місяців тому
#IAST #excel #adsorption Ideal Adsorbed Solution Theory (IAST) Selectivity Combined Langmuir Freundlich Model. Using Microsoft Excel Solver
Breakthrough Curve Desorption Solved by Method of Line GNU Octave (PDE-ODE)
Переглядів 2428 місяців тому
Breakthrough Curve Desorption Solved by Method of Line GNU Octave (PDE-ODE)
PDE Burgers’ Equation GNU Octave Method of Line
Переглядів 2058 місяців тому
PDE Burgers’ Equation GNU Octave Method of Line
Breakthrough Curve Adsorption Solved by pde1dm GNU Octave. Solved PDE-ODE (coupled system)
Переглядів 3448 місяців тому
Breakthrough Curve Adsorption Solved by pde1dm GNU Octave. Solved PDE-ODE (coupled system)
Breakthrough Curve of Adsorption and IAST Selectivity in Binary Liquid System. Method of Line MATLAB
Переглядів 2818 місяців тому
Breakthrough Curve of Adsorption and IAST Selectivity in Binary Liquid System. Method of Line MATLAB
Boundary Value Problem (BVP) | Steady State Mass Transfer Diffusion-Advection-Reaction
Переглядів 1729 місяців тому
Boundary Value Problem (BVP) | Steady State Mass Transfer Diffusion-Advection-Reaction
Boundary Value Problem (BVP) | Steady State Heat Transfer in Rod
Переглядів 2289 місяців тому
Boundary Value Problem (BVP) | Steady State Heat Transfer in Rod
Method of Line - PDE Burgers’ Equation
Переглядів 3589 місяців тому
Method of Line - PDE Burgers’ Equation
Non-Linear Curve Fitting - Indirect Fitting : Bacterial Growth Kinetic Rate
Переглядів 1819 місяців тому
Non-Linear Curve Fitting - Indirect Fitting : Bacterial Growth Kinetic Rate
Non-Linear Curve Fitting using MATLAB : Finding Parameter in Virial Equation
Переглядів 3659 місяців тому
Non-Linear Curve Fitting using MATLAB : Finding Parameter in Virial Equation
Vapor Liquid Equilibrium (VLE) Binary System Formic acid-Acetic acid Using UNIQUAC (MATLAB)
Переглядів 4139 місяців тому
Vapor Liquid Equilibrium (VLE) Binary System Formic acid-Acetic acid Using UNIQUAC (MATLAB)
Boundary Value Problem (BVP) | Analytical vs Numerical Solution
Переглядів 1219 місяців тому
Boundary Value Problem (BVP) | Analytical vs Numerical Solution
UNIQUAC Activity Coefficient & Excess Gibbs Energy Binary Mixture Benzene and 2,2,4-Trimethylpentane
Переглядів 3449 місяців тому
UNIQUAC Activity Coefficient & Excess Gibbs Energy Binary Mixture Benzene and 2,2,4-Trimethylpentane
Elliptic PDE-Laplace Equation. Concentration Distribution in the Tank
Переглядів 2409 місяців тому
Elliptic PDE-Laplace Equation. Concentration Distribution in the Tank
Elliptic PDE-Laplace Equation. Heated in the center of the square plate
Переглядів 1709 місяців тому
Elliptic PDE-Laplace Equation. Heated in the center of the square plate
Elliptic PDE-Laplace Equation Heat Transfer in Square Plate insulated. Neumann Boundary Condition
Переглядів 7059 місяців тому
Elliptic PDE-Laplace Equation Heat Transfer in Square Plate insulated. Neumann Boundary Condition
hi, did you post the code somewhere ?
Sorry, I did not. If you need code, please let me know. Thx
great video, but the background music is very annoying.
I am so sorry 🙏🏻
Hi, can you helpme? please i can´t solve a system of PDE Cylindrical are six PDE of reaction
Thank you for sharing useful knowledge
my pleasure sir
my pleasure sir
Thank you for sharing your knowledge. You're code helped me a lot with my thesis
Really? I am happy to hear that 😁
hello sir, may I ask. I'm currently casting the membrane, but I'm having a problem. After I cast it and then put it in the coagulation bath, my membrane sticks to the glass plate. Is there a solution that I can do to prevent the membrane from sticking to the glass plate? thank you..
You should sink the casting membrane about 3 minute before you pick up your membrane. Did you let the membrane totally under coagulation bath in 3 minutes?
@@rahmatsunarya yes, i even sink the membrane more than 3 minutes because it was still sticking to the glass plate. anw i am using cellulose acetate and NMP solvent to make the membrane
Please tell me the detail composition your membrane solution at rahmatsunarya@gmail.com
@@rahmatsunarya i've emailed you, thank you so much for helping me
I just answered your email. Kindly check
Hi thanks for the video, it's a very basic question but s it distilled water in the bath you drop the glass? if it isn't, which liquid is it?
that's right. I used demin water
hello sir, i need thé code of the parametrical study of this same situation do u havé it please?? is about study thé effect of different parametrs on the breakthrough curve, is in the article source u put in this video. can u help me please i havé a presentation next week and i am really struggling 😢
Please get me at rahmatsunarya@gmail.com
Hi, i like your videos they are very usuful, i have a question about the breackthrough curve of adsorption can u help me please? Please answer me as fast as u can, i have a presentation in two weeks and i struggle to solve the problem in matlab, i am a biginner. Thank u so much
Hi If you have question, please tell me at rahmatsunarya@gmail.com thanks
@@rahmatsunarya I sent you an email. Thank you so much
Can you please check your email? I sent you a private message. Thank you.
Hi, at 6:25, why T(j-1,i) not T(j,i). You explained at 2-1=1 but I don't understand why. As you specified, T(1,:)=28 so j=2 will result T(1,i)=28 =B(i).. I am so confused.
because of j starting from 2. in other hand we need as j=1, so, we put j-1.
I have a code for my research problem.I am struggling to get a plot...can you help me. it's a fdm code.
Fdm code? It's about PDE or ODE?
Could you pass me those books?
I'm sorry. It's an illegal way
Hello, this video is very helpful. I have a question, How would you write the boundary condition matrix A and B if there are nonconstant flux terms such as radiation and convection which depend on the length(t)? Thank you!
Please let me know your case in my email
@@rahmatsunaryaYes, I will. Thank you, sir.
Sorry, could you share me your script, please? I had problems with the funtion @fres
Surely, please get me in gmail 🙏🏻
Very good video! How can i include a non-linear ecuation in a boundary condition ?
Yes surely Sir
Can Orthogonal Collocation method solve this problem with axial diffusion items？
I think it can't be solved using orthogonal collocation method.
Thanks this video is very good
My pleasure
if I have a question about the sheet, can I contact you?
Please tell me
Halo kak kenapa nilai Nt bernilai 101 dan nilai Nz 151 bagaimana cara mendapatkan kedua nilai tersebut?
Halo. Itu bisa disesuaikan. Tidak harus segitu. Bole Nt= 50, Nz= 20. Nilainya bebas. Itu hanya menunjukkan pembagian range waktu dan jarak. Trims
@@rahmatsunarya Peran Nt dan Nz pada program itu sebagai ap ya kak?
Jika nilainya makin besar, makan hasilnya lebih akurat. Jika nilainya makin kecil, hasilnya tidak begitu akurat. Dan juga apabila nilainya besar tentu akan memakan waktu perhitungan yg lama.
@@rahmatsunarya siap terimakasih kak
Terima kasih kembali 😊
Why Nt=101 and Nz=151 sir?
It's a number for meshing time and space
At 3.18 how are you doing ? Am not getting. Please explain
Make a function, turn PDE to ODEs function
Dear, just a question, what kind of Matlab License you need to upload these beautiful videos? Regards.
Student license sir
excelent video really enjoy you’re content, could next video be about volume control solved by elliptic equations? Maybe solve problem 29.9 on chapras book 7th edition
Is it PDE in rectangular case?
@@rahmatsunarya yes!
Hmm it is a irreguler mesh point. What do you think ?
@@rahmatsunarya i have the equations for the desired nodes, the thing im struggling with is making the matlab code 😓
@arathcamacho1247 would u please to tell me?
Hy.. I tried putting the same code on matlab and it's not working.. I am actually trying to solve a system of 4PDE in matlab. ( Absolutely new to this) I feel like this video is very similar to the problem..so I was trying to draw parallels.. unfortunately I couldn't run the code in matlab.. Can you help me out too?..
Sure. Please let me know
@@rahmatsunarya Error using vertcat Dimensions of arrays being concatenated are not consistent. Error in self3>f (line 124) dydt=[dc1dt;dc2dt;dc3dt;dq1dt;dq2dt;dq3dt]; Error in odearguments (line 92) f0 = ode(t0,y0,args{:}); % ODE15I sets args{1} to yp0. Error in ode15s (line 148) odearguments(odeIsFuncHandle, odeTreatAsMFile, solver_name, ode, tspan, y0, options, varargin); Error in self3 (line 36) [t, y] =ode15s(@f,t,IC,[],eps,u,LDF,b1,b2,b3,qm1,qm2,qm3,c10,c20,c30,rhop,Nz,dz); >> these kinda errors i am getting .can i contact you via gmail?
@@rahmatsunarya i am getting erros like Error using vertcat Dimensions of arrays being concatenated are not consistent. Error in self3>f (line 124) dydt=[dc1dt;dc2dt;dc3dt;dq1dt;dq2dt;dq3dt]; Error in odearguments (line 92) f0 = ode(t0,y0,args{:}); % ODE15I sets args{1} to yp0. Error in ode15s (line 148) odearguments(odeIsFuncHandle, odeTreatAsMFile, solver_name, ode, tspan, y0, options, varargin); Error in self3 (line 36) [t, y] =ode15s(@f,t,IC,[],eps,u,LDF,b1,b2,b3,qm1,qm2,qm3,c10,c20,c30,rhop,Nz,dz); >> for the same code youre using. Can i contact you via gmail ?
Surely, rahmatsunarya@gmail.com
@@rahmatsunarya I have mailed you
Hello sir. Can I solve FTFSCS method by using this formula. Please respose.
What is FTFSCS?
Awesome video
thank you sir
Pls send your mail id and ph no
Thanks a lot. Excellent work
My pleasure
I see that the units in the model does not match, when q is in kg/kg..
I am sorry for the mistake. But in this video, I want to share how to solve PDE-ODE. That's all. I am not focused on unit or unsubstantial content except solving the transport equation.
Can you link me the paper wher you get the data from or explain what is yd, xd and ntotd?
Please find the book for the data : Duong D Do (1988) Adsorption Analysis: Equilibria and Kinetics.
Thank you, it helped in my class presentation. I kept your video on references.
Wow. Amazing. I glad to hear that 😊😊
hi sir how can we make this and what material use for this kindly share details with weight of material used in it for making
Made by 20% polysulfone and 80% NMP in 50 gram solution
@@rahmatsunarya sir still confused kindly explain simply and where i got these chemicals
10 gram PS mixed well in 40 gram NMP. Then you have a polymer solution. Now, the solution is casted like my video
@@rahmatsunarya thanks sir and can i learn from you to make a membrane
Sure. If you have some questuions again, let me know at rahmatsunarya@gmail.com
Great work
Thanks
Indeed, I have the same error. Error in ode15s (line 153) odearguments(odeIsFuncHandle, odeTreatAsMFile, solver_name, ode, tspan, y0, options, varargin); Error in Sips (line 40) [t y]=ode15s(@fsips,t,IC,[],eps,u,LDF1,LDF2,c10,c20,rhop,rhog,Nz,dz,Cp,Cpg,delH1,delH2,Qw)
Please let me know at my gmail. Thanks
how would you do this for multicomponent ?
Maybe you can use another method, such as method of line.
Hello, nice work, could You share the paper where you found the binary interaction parameters?
I found it from internet.
it depends on what kind of mixture you want to evaluate. The basics are in ismail Tosun, 2021. There are special consideration of binary interaction parameter calculation in Gao, 2008, included in ismail Tosun, 2021. or choose the older book of ismail Tosun, 2013 for manual calculation. You can search other formula for binary interaction parameter (a mix, b mix, or k_ij) elsewhere on the internet.
Thank you for explanation 👍
Which paper is the equation referenced from?
Hi. I found out from my viewer, who emailed me.
MINE GROUP 4
Thanks
MINE GROUP 4
Thanks
Jacinth Castro
Thanks
Thanks
¡Hola! Me gustaría saber si tienes algún blog o página donde subas las presentaciones, para si es posible descargarlas ¡Gracias!
I am sorry, could you please use proper English? I can't understand, sorry
@@rahmatsunarya Hi! I would like to know if you have any blog or page where you upload the presentations, so if it is possible to download them Thank you!
@@gabrielomarperezmeza3079 I will send it to you, please send me email at Rahmatsunarya@gmail.com
Where did you get the binary interaction parameters of the mixture?
I found it from paper, sir.
@@rahmatsunarya Thank You for your answer
My pleasure bro
@@rahmatsunarya Hello, can you share the name of the paper?. Thanks in advance!
@@JOSEANTONIO198207 oh sorry, after I checked again, I found it from internet and I download as excel file.
Awesome video ! Is it working for chromatography column ?
Hello. I am not familiar with chromatography. If the process is adsorption, I think it will work in chromatography. Thanks
Can i get code ?
surely 😄
Please brother Share
Send me an email at rahmatsunarya@gmail.com
Thank you so much your videos have been so helpful you saved my degree
My pleasure Nouran
I need your help , how to contact you
This is my email rahmatsunarya@gmail.com
Good work
My pleasure 🙏🏻
Hi Sir, really thanks for this video. I tried the same code, but got a different graph. I sent you the code in gmail. Would you please help me to find out where I went wrong? Again I really appreciate all your videos.
Surely, which on your email. I cant see your name JollyG1rl in my inbox 😄
@@rahmatsunarya hi Sir I sent an email again just now
The email subject would be "regarding an issue in MATLAB code"
What is z?
Hello. z is UNIQUAC coordination number