Hi, can anyone help me to program a code for the equation with initial and boundary conditions using finite difference method
5 views (last 30 days)
Show older comments
Hi, I really need help from anyone help me with how to do a program code for governing equation with initial and boundary condition. Also, I want to plot graph velocity(u) against radius (r). The radius will be the x-axis and the range should be from 0 to 1 only. I will be really appreciated if anyone can help me. Thank you so much. Inside the coding there must be include discrerization of the governing equation. If u need , I wiil send it here right away.
The governing equation is:
![](https://www.mathworks.com/matlabcentral/answers/uploaded_files/848755/image.png)
![](https://www.mathworks.com/matlabcentral/answers/uploaded_files/848760/image.png)
0 Comments
Answers (1)
William Rose
on 31 Dec 2021
@Nur Nadhirah Syed Malik, See attached.
I have left out some parts, because this sounds like a homework assignment, and you will learn by filling in the unknown parts.
4 Comments
William Rose
on 31 Dec 2021
A script is attached that implements the equaitons in the document I posted earlier.
The script generates two figures, shown below:
![](https://www.mathworks.com/matlabcentral/answers/uploaded_files/848875/image.jpeg)
![](https://www.mathworks.com/matlabcentral/answers/uploaded_files/848880/image.jpeg)
William Rose
on 31 Dec 2021
I believe that the constant
in your equation corresponds to kinematic viscosity:
=viscosity/density. I also believe that the units for A0 and A1 are pressure/length, and that the terms A0+A1*cos wt correspond to the pressure difference divided by the length of the tube. Are those beliefs correct? The Womersley number for pulsatile flow is
. If I am correct about
, then the Womersley number for this problem (in which R=1) is
.
![](https://www.mathworks.com/matlabcentral/answers/uploaded_files/848885/image.png)
![](https://www.mathworks.com/matlabcentral/answers/uploaded_files/848890/image.png)
![](https://www.mathworks.com/matlabcentral/answers/uploaded_files/848895/image.png)
![](https://www.mathworks.com/matlabcentral/answers/uploaded_files/848900/image.png)
![](https://www.mathworks.com/matlabcentral/answers/uploaded_files/848905/image.png)
See Also
Categories
Find more on 2-D and 3-D Plots in Help Center and File Exchange
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!