Page 1 of 1

imposing pressure wave as inlet using GenBC

Posted: Sat Jul 04, 2020 2:53 pm
by masihbln
Hi SimVascular team,
I am intending to run a simulation featuring pressure BC as inlet using GenBC. there is tutorial on this on website but compared to setting RCR BC,imposing pressure wave BC as inlet is not completely elaborated. it is mentioned that svpre should be modified with neuman pressure wave and to run presolver but its implementation is not covered and editing fortran codes isn't presented.
More specifically, how do svpre file and fortran codes should be accurately modified?
along with that, there is a topic in which SimVascular member prepare a file as steady pressure example but the link doesn't work.
Would absolutely appreciate if the clear corresponding tutorial or example for pressure wave inlet BC would be provided.
Thank you SV great Team

Re: imposing pressure wave as inlet using GenBC

Posted: Tue Jul 07, 2020 7:34 pm
by wgyang
Try the GenBC.f and USER.f attached and see if they work. Basically, GenBC asks users to provide a subroutine to calculate f=dX/dt, where X are ODE unknowns for the lumped parameter model. To impose a known pressure wave, I can set f=0 and update genBC unknown X directly, which is the pressure in this problem. In addition, you can replace the linear interpolation with an analytical Fourier series for the pressure waveform.