Plotting streamlines uisng axistream.h for axisymmetric domain

[atp...@gmail.com]

Hello all,

I have been trying to plot streamlines of flow past an axisymmetric droplet. The

simulations have already been done and I am plotting streamlines by

post-processing the dump files generated during the simulation.

I have attached my Basilisk script below for reference. I have called restore()

as well as done the field variable assignments in main() itself. This means that

no events are created.

I have also explicitly initialized the axisymmetric metrics since run() is

never called and hence metric event in axi.h is never called.

I am using axistream() in src/axistream.h for getting streamfunction scalar

field. I have followed src/examples/yang.c for this purpose.

I have attached the streamlines generated by basilisk as well as velocity

field and vectors plotted using paraview for a dumpfile.

Velocity field and vectors:

Corresponding streamlines generated using Basilisk:

Streamlines should always be tangential to the velocity field. However, the

streamlines generated using Basilisk are not.

It will be a huge help if someone could look in to this issue of mine.

Many Thanks,

Aditya Parik