## WSS calculations directly from SimVascular

Provides a system for patient-specific cardiovascular modeling and simulation.
Atefeh Razavi
Posts: 29
Joined: Sun Mar 04, 2018 10:47 am

### WSS calculations directly from SimVascular

Hello,

I was wondering to know how I may extract (calculate) WSS values directly from SimVascular not Paraview.

Thank you and looking forward to your reply,
Ati

David Parker
Posts: 225
Joined: Tue Aug 23, 2005 2:43 pm

### Re: WSS calculations directly from SimVascular

Hi Ati,

WSS is calculated from the simulation results using the svPost program, which is part of the svSolver programs. You could have a look at the source code to see how the values are calculated (see https://github.com/SimVascular/svSolver ... ver/svPost).

Cheers,
Dave

Atefeh Razavi
Posts: 29
Joined: Sun Mar 04, 2018 10:47 am

### Re: WSS calculations directly from SimVascular

Hi Dave!

I looked at the code, however due to high number of variables, I did not figure out which variables were used to calculate WSS.
As defined, WSS is the velocity derived tangential stress. I was interested to see if I can calculate WSS from the velocity as well as the normal to the cell and then compare it with the WSS values given in Paraview.

Thanks,
Atefeh

David Parker
Posts: 225
Joined: Tue Aug 23, 2005 2:43 pm

### Re: WSS calculations directly from SimVascular

Hi Atefeh,

Sorry, I mislead you! The WSS calculation is actually done in the solver in solver_subroutines.f at line 1667:

Code: Select all

``````
rtmp = zero
wallssVec  = zero

do i=1, numnp
if (invflx(i) .ne. 0) then
rtmp(i,2:4) = flxres(i,1:3) !viscous flux
c     calculate the WSS
tn = flxres(i,1) * flxnrm(i,1)
&            + flxres(i,2) * flxnrm(i,2)
&            + flxres(i,3) * flxnrm(i,3)

wallssVec(i,1) = flxres(i,1) - tn * flxnrm(i,1)
wallssVec(i,2) = flxres(i,2) - tn * flxnrm(i,2)
wallssVec(i,3) = flxres(i,3) - tn * flxnrm(i,3)
endif
enddo
``````
where flxnrm is the outward normal and flxres is the viscous flux.

Cheers,
Dave

Atefeh Razavi
Posts: 29
Joined: Sun Mar 04, 2018 10:47 am

### Re: WSS calculations directly from SimVascular

Hi Dave,

Thanks for the info.
I was wondering to know if it is possible and how I may use this code directly (probably in SV console) to extract WSS.
Another question is when using Paraview, if a line is intersected with a radial cross-section, instead of having presumably two WSS values (two points on the wall), a series of WSS values along that line are reported. I appreciate your answer.

Thank you,
Ati

David Parker
Posts: 225
Joined: Tue Aug 23, 2005 2:43 pm

### Re: WSS calculations directly from SimVascular

Hi Ati,

If you want to calculate wss yourself then you would need to incorporate the code I posted into a program that could read in SV results and mesh files, calculate surface normals, viscous fluxes and then wss. This would be quite involved.

If you would like to just look at the wss values then you can write them to a text file using ParaView. If you could explain to me what you would like to do with the wss values then maybe I could write a Python script to do what you want.

I believe the ParaView line probe samples the mesh at several points along the line.

Cheers,
Dave

Atefeh Razavi
Posts: 29
Joined: Sun Mar 04, 2018 10:47 am

### Re: WSS calculations directly from SimVascular

Hi Dave,

I would like to check WSS results read in Paraview with a step- by step extracted one,as you explained, in a program (that could be what?).
The reason is that I have seen regions of high WSS (compared to previously reported data for the coronary simulations) and try to make sure they are accurate.
I need to mention that I have done mesh study as well (with the flow and pressure matching clinical data) and the values do not change (less than 3%) for the consecutive meshes.

Thank you,
Ati

Atefeh Razavi
Posts: 29
Joined: Sun Mar 04, 2018 10:47 am

### Re: WSS calculations directly from SimVascular

About the wss along a line, I still don't understand why more than one point on each end of the line should be reported.
The image attached showed the velocity and wss along a segment of the line (all the pink dots are WSS values which increase dramatically along the line). I would expect to see one value for WSS and then zero for the rest of the mesh. So I am not sure if the highest value of these is the one showed in WSS contours.

Thank you,
Ati
Attachments vel_WSS_2.JPG (47.37 KiB) Viewed 253 times

David Parker
Posts: 225
Joined: Tue Aug 23, 2005 2:43 pm

### Re: WSS calculations directly from SimVascular

Hi Ati,

I will write a program to extract the wss from the SV results files, this is something I've been wanting to do anyway.

If you are reading in.vtu files to visualize wss then the wss results are calculated for all nodes in the mesh but are only non-zero on the surface. So the line probe will sample wss values in the mesh volume. The wss will drop from a non-zero value on the surface to a zero value within an element. You can see this in the attached image, where the wss fades from red to blue within a single element.

Cheers,
Dave
Attachments wss-probe.png (56.68 KiB) Viewed 229 times

Atefeh Razavi
Posts: 29
Joined: Sun Mar 04, 2018 10:47 am

### Re: WSS calculations directly from SimVascular

Hi Dave,

Thanks for the answer and it would be great to extract wss from SV results.
So there is a range of WSS values not just WSS on the surface. So which value of this range is reported if a vtp file is used? just the maximum on the surface? or the mean of the range within the element?

Thank you,
Ati