- Видео 138
- Просмотров 231 063
Rahmat Sunarya
Добавлен 15 сен 2010
It's all my lecture notes.
Please enjoy it.
Please enjoy it.
Heat Transfer with nonconstant initial condition. Method of Line MATLAB
#mathematics #pde #ode #heat #transfer
Heat Transfer with nonconstant initial condition. Method of Line MATLAB
Heat Transfer with nonconstant initial condition. Method of Line MATLAB
Просмотров: 230
Видео
Diffusion - Advection Partial Differential Equation Solved by pdepe MATLAB (PDE)
Просмотров 2894 месяца назад
Diffusion - Advection Partial Differential Equation Solved by pdepe MATLAB (PDE) clc;clear;close all; %meshing time & space t=linspace(0,60,61); %time z=linspace(0,1,100); %space (axial) m=0; %remember m=0 sol=pdepe(m,@pdefun,@icfun,@bcfun,z,t); C=sol(:,:,1); imagesc(z,t,C) colormap jet colorbar title('Concentration profile') xlabel('position (cm)') ylabel('time (s)')
Friction Factor Colebrook Equation (MATLAB) using fzero syntax
Просмотров 1375 месяцев назад
#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
Dimensionless Generic Breakthrough Curve
Просмотров 1946 месяцев назад
#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
Просмотров 1997 месяцев назад
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
Просмотров 1137 месяцев назад
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
Просмотров 1137 месяцев назад
#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
Просмотров 3238 месяцев назад
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
Просмотров 3878 месяцев назад
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)
Просмотров 40011 месяцев назад
#Burgers #Huxley #waveequation #pde #matlab
Partial Differential Equation (PDE) Fix Bed Reactor with The Danckwerts boundary conditions
Просмотров 38611 месяцев назад
#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
Просмотров 21311 месяцев назад
#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
Просмотров 448Год назад
#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
Просмотров 666Год назад
#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.
Просмотров 323Год назад
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.
Просмотров 610Год назад
Redlich - Kister Expansion One Parameter | Curve Fitting Vapor-Liquid Equilibrium of Binary System.
IAST Selectivity Using Excel Solver Combined Langmuir-Freundlich Model
Просмотров 329Год назад
IAST Selectivity Using Excel Solver Combined Langmuir-Freundlich Model
Breakthrough Curve Desorption Solved by Method of Line GNU Octave (PDE-ODE)
Просмотров 305Год назад
Breakthrough Curve Desorption Solved by Method of Line GNU Octave (PDE-ODE)
PDE Burgers’ Equation GNU Octave Method of Line
Просмотров 277Год назад
PDE Burgers’ Equation GNU Octave Method of Line
Breakthrough Curve Adsorption Solved by pde1dm GNU Octave. Solved PDE-ODE (coupled system)
Просмотров 421Год назад
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
Просмотров 353Год назад
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
Просмотров 199Год назад
Boundary Value Problem (BVP) | Steady State Mass Transfer Diffusion-Advection-Reaction
Boundary Value Problem (BVP) | Steady State Heat Transfer in Rod
Просмотров 270Год назад
Boundary Value Problem (BVP) | Steady State Heat Transfer in Rod
Non-Linear Curve Fitting - Indirect Fitting : Bacterial Growth Kinetic Rate
Просмотров 206Год назад
Non-Linear Curve Fitting - Indirect Fitting : Bacterial Growth Kinetic Rate
Non-Linear Curve Fitting using MATLAB : Finding Parameter in Virial Equation
Просмотров 575Год назад
Non-Linear Curve Fitting using MATLAB : Finding Parameter in Virial Equation
Vapor Liquid Equilibrium (VLE) Binary System Formic acid-Acetic acid Using UNIQUAC (MATLAB)
Просмотров 479Год назад
Vapor Liquid Equilibrium (VLE) Binary System Formic acid-Acetic acid Using UNIQUAC (MATLAB)
Boundary Value Problem (BVP) | Analytical vs Numerical Solution
Просмотров 153Год назад
Boundary Value Problem (BVP) | Analytical vs Numerical Solution
UNIQUAC Activity Coefficient & Excess Gibbs Energy Binary Mixture Benzene and 2,2,4-Trimethylpentane
Просмотров 445Год назад
UNIQUAC Activity Coefficient & Excess Gibbs Energy Binary Mixture Benzene and 2,2,4-Trimethylpentane
Elliptic PDE-Laplace Equation. Concentration Distribution in the Tank
Просмотров 273Год назад
Elliptic PDE-Laplace Equation. Concentration Distribution in the Tank
thanks for your explanation
Thanks!
@@RHTQ1 Oh my God. I'm crying for this 😭😭 it's my superthank ever in my YT carier . Thank you so much for your incredible support!
For some reason, I still didn't need fsolve when I switched to a second order reaction, which feels wrong. But that's likely a me problem. TYSM for this awesome video! Seriously, this is better than the terrible resources I have been paying for. Plus there's calm music...
@@RHTQ1 Really? Youur comment really made my day! It gave me such a boost. I'll do my best to make more content and help others 😊😄
Hi, I tried to run this code in version R2024a and R2020b but it doesn't seem to work anymore. Do you know what might be the problem?
I think there's a problem with how lsqnonlin works differently in MATLAB 2024 compared to 2020.
@@rahmatsunarya Thanks for responding. I already discovered a solution by defining the objective function for lsqnonlin as an anonymous function. Line 12 looks like this then: ...= lsqnonlin(@(p) fres(p, t, IC, c), g, lb, ub);
@@KobePeeters-x8i Great. Thank you for sharing
you save the world by making this tutorial series man!
Wow, I’ll let the superheroes know they’re on vacation now! Thanks a lot! 💪🚀
I don't have a Matlab license. Do you know of an open source method for doing this?
Surely. You can use OCTAVE. It's a open source software and Matlab alternative. You can see the tutorial at folowing link ruclips.net/video/cq-XjX7yLD0/видео.htmlsi=giFUA1OamlzN8D7j
@@rahmatsunarya, thanks much. Do you know if there is a Python version?
Why are we interested in finding A12 and A21?
It involved molecular attraction and was useful for calculating the activity coefficient. That's all, sir. Thank you for your comments
sir, can i save the figure result by timestep as video file ?
Hello. Thank you for comment. Sure, you can get the video file. Please see the tutorial in another my videos: ruclips.net/video/XZ0A5qTX5Ko/видео.html
great thanks.............may I get the code?
Sure please send me a email in rahmatsunarya@gmail.com
@@rahmatsunarya great thanks for your help
I am unsure, but at 1:52, the f should be d2Cdx2 instead of dCdx?
Sorry Sir, I think I do same as the MATLAB template for pdepe 🙏
Nice work.. Please simulate this type of advective model from the journal: DEEP-BED FILTRATION WITH MULTISTAGE DEPOSITION KINETICS by Vitaly Gitis etal, 2010.
Sir Please send me your mail for contact.
Hi, are you using Jacobi method for iteration, right? Your videos are very good!
Hi. I'm using simple iteration in matlab. I do not know it was Jacobi method or not 😅
cocok untuk mahasiswa Indonesia
Thank You sir. Appreciate
Can anybody help to find out henry's selectivity from static adsopttion isotherm
did you have a similar problem?
Nice to see your work again
thank you sir for you support
Awesome
thank you sir
Thanks for your video! It's super helpful. Quick question though - I noticed the formula in the slide at 0:52 looks a bit off from the Peng-Robinson equation I'm familiar with. The formula shows (1-\kappa\frac{T}{T_c})^0.5, but I think it should be (1-\kappa(\frac{T}{T_c})^0.5). The code's right, so maybe just a typo in the slide?
OMG, I got typo. That was (1-(T/Tc)^0.5)^2. I am so sorry. Thank you for noticed me
Nice 👍 What about a system of multistage deposition. Where there is more than 1 equations for the rate of particle adsorption and desorption. And also epsilon is time dependent... Very nice and wonderful video Sir. Thank you
My pleasure sir
This video is quite helpful, but i have a question. What if epsilon(porosity) is a function of time, And rate of change of q(particle deposition) has more than one stage time dependent equation. Example in "Deep-bed filtration model with multistage deposition kinetics" by vitaly Gitis etal. How do we code such model. Thank you🎉
Thanks Very much for this video
My pleasure
man, I'm looking for a textbook so I can understand this frfr
Transport phenomena by Bird, Lighfoot
@@rahmatsunarya tysm sir
@Spectrojamz good luck
How about ruclips.net/video/zcDZk1pqyrc/видео.html
can you show the question plz?
Hello. It's not from question, but from paper. Thank you
How about ruclips.net/video/zcDZk1pqyrc/видео.html
Mas Rahmat Sunarya, I am currently trying to develop a direct contact membrane distillation model using an advection-diffusion approach. If you have some time, I would like to ask you a few questions. 😅 The approach I am using is the finite difference method. I am having difficulty because, on one side, there is a boundary condition in the form of a function.The function itself need temperature from interior nodal..
could you please tell me on gmail, Sir?
Thanks for this video. Please can you explain me why in the BC-1 and BC-2 during discretization you divide for 2 (forward and backward 1st order seems not require that)?
BC-1 is known left BC or how we start. So, we use forward. BC-2 is known right BC or how we finish. So, we use backward.
I hope you get the point.
@@rahmatsunarya Sorry the question was ill-posed. I was asking why in forward discretization (or backward) you divide by 2 in BC (at denominator you have the number 2 for both BC).
Omg, it's my mistake. I am sorry. It is not divided by 2. Thanks for your comment
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
Surely sir
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