Output variables of the DualSPHysics5.0_NNewtonian

Hello DualSPHysics community.

I'm trying to simulate Non-Newtonian fluid using DualSPHysics5.0_NNewtonian.

The simulation has been successfully applied. We wanna evaluate the shear stress, the magnitude shear rate and the viscosity by the Hershel-Bulkley-Papanastasiou model using "measuretool", but these variables are not included in bi4 files.

Although we tried to modify the source, function JSphGPU::ParticlesDataDown in JSphGpu.cpp, we have no idea which variables are which...

The choice of reology treatment in Def.xml file are followings:

  1. RheologyTreatment: Single and multi-phase
  2. VelocityGradientType: SPH
  3. ViscoTreatment: Laminar+SPS

I would appreciate it if you could give me advice.



  • Hi Sato,

    Thank you for your message.

    In the non-Newtonian code (CPU and GPU) we have intentionally included an array named auxnn which can be used to output any (scalar) variable you like. You will need to uncomment the array in the code (assign the value you like) and pass it through to savedata().

    There is also a commented example with pressure in the code (please follow the savedata() function).

    if not clear please email the dualsphysics team and I might be able to help you out (however, I am bit on the busy side, so you might have to wait a few days).



  • Dear Dr.George,

    Thank you for your reply.

    I am trying to build the source uncommenting the variables.

    I might send email the team if it doesn't work.


    K. Sato,

Sign In or Register to comment.