Yes, I also realized that it could be explained by the conservation of energy. For a short diverging pipe, the energy loss can be neglected. The flow slows down when it enters a region with a larger cross-sectional area. The kinetic energy is converted into a potential. So the dynamic pressure is re...

Regarding solver.inp, I don't notice anything particularly wrong. Could you please specify what parameters you set in the svpre file? I may try to downsample your inflow.flow from 781 to approximately 100 data points while applying some smoothing. For the svpre file, consider reducing the number of ...

The model is quite short relative to its diameter. I am wondering if you have tried an expansion vessel with a larger length. My intuition is that the prescribed inflow (125 ml/s) has a sufficient momentum to reach outlet even though outlet pressure is higher by 150 dyn/cm^2 which is only 0.1 mmHg. ...

Could you grant me the access to your files so that I can take a look?
Thanks!

Impedance BC cannot be prescribed via the SV GUI. It is not used as often as RCR. To use impedance BC, you need to prepare two files impt.data and Qhistor.dat in order to compute the impedance convolution. impt.dat has a format similar to rcrt.dat. Each block represent a set of impedance for an outl...

I expect to see a pressure difference of 20 Pa instead of 40 Pa. Can you check if the two rcrt.dat files differ in the last two rows only. Take the rcrt file used in the Simvascular cylinder demo project as an example. The original Pd was prescribed as follows, 0.0 0.0 1.0 0.0 To change to Pd=20Pa, ...

The rules for distributing resistance are based on the assumption that flow is proportional to A^alpha. When alpha equals 1, the outflow is proportional to the outlet area. Similarly, the capacitance is proportional to the outlet area. The total capacitance is the sum of each outlet's capacitance wh...

To my knowledge, the resistance split doesn't take account of the 3D model resistance of the 3D model. It just distributes the total resistance you provide following a power law. If we assume the pressure across the resistors is the same (i.e. in parallel), the resistance is split such that Q is pro...

Did you also use the resistance BC for other CFD software packages? You can start with a zero pressure BC for all software packages you want to compare as the baseline. Then consider comparing resistance or RCR. For commercial CFD software packages, did you use a UDF function to implement the resist...

Maybe you can do it in the post-processing. Since the wall nodal velocity is available, apply the membrane linear constitutive relationship to get the wall stress tensor: \sigma=C:\epsilon, where C is the material stiffness tensor and \epsilon is the strain rate tensor. See eqns 32-35 in Figueroa et...