Hi, can anyone help me to program a code for the equation with initial and boundary conditions using finite difference method
조회 수: 3 (최근 30일)
이전 댓글 표시
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
답변 (1개)
William Rose
2021년 12월 31일
@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
William Rose
2021년 12월 31일
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
2021년 12월 31일
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)
참고 항목
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!