Skip to content

visualization.md

Latest commit

 

History

History
95 lines (72 loc) · 3.6 KB

visualization.md

File metadata and controls

95 lines (72 loc) · 3.6 KB

Visualization

Export VTK Files

You can export particle data as VTK files by using the SolutionSavingCallback. All predefined examples already use this callback to export VTK files to the out directory relative to the current working directory. VTK files can be opened in visualization tools such as ParaView and VisIt.

ParaView

Follow these steps to view the exported VTK files in ParaView:

  1. Click File -> Open.
  2. Navigate to the out directory (relative to the directory that you are running Julia from).
  3. Open both boundary_1.pvd and fluid_1.pvd.
  4. Click "Apply", which by default is on the left pane below the "Pipeline Browser".
  5. Hold the left mouse button to move the solution around.

You will now see the following:

image

It is useful to make the dot size dependent on the actual particle size. For this, first make sure you have "fluid_1.pvd" highlighted in the "Pipeline Browser". Then, in the Properties panel (bottom left), adjust the following settings:

  1. "Representation" to "Point Gaussian".
  2. Choose the right "Shader Preset": "Plain Circle" for 2D and "Sphere" for 3D.
  3. Activate "Scale by Array" and select "particle_spacing" in "Gaussian Scale Array".
  4. Deactivate "Use Scale Function".
  5. Set the "Gaussian Radius" to "0.5".

image

Visualization with a Macro

To simplify visualization of particle data in ParaView, you can use a macro. It reduces the manual steps from the previous section to a single click. Install the macro as follows.

  1. Save the macro code (see below) as a .py file, e.g. PointGaussianMacro.py.
  2. Open ParaView and go to the top menu: Macros → Import New Macro... → Select the saved .py file.
  3. The macro will now appear in the Macros menu and can optionally be pinned to the toolbar.
  4. Load your dataset into ParaView.
  5. Select the dataset in the Pipeline Browser.
  6. Click on the macro name in the Macros menu (or toolbar, if pinned) to run it.
  7. The Point Gaussian representation with particle_spacing scaling will be applied automatically.

Macro Code

# trace generated using paraview version 5.13.1
#from paraview
#paraview.compatibility.major = 5
#paraview.compatibility.minor = 13

from paraview.simple import *
paraview.simple._DisableFirstRenderCameraReset()

# get active source
source = GetActiveSource()

# get active view
renderView1 = GetActiveViewOrCreate('RenderView')

# get display properties
sourceDisplay = GetDisplayProperties(source, view=renderView1)

# change representation type
sourceDisplay.SetRepresentationType('Point Gaussian')

# modified display properties
sourceDisplay.ShaderPreset = 'Plain circle' # for 2D, change to 'Sphere' for 3D
sourceDisplay.ScaleByArray = 1
sourceDisplay.SetScaleArray = ['POINTS', 'particle_spacing']
sourceDisplay.UseScaleFunction = 0
sourceDisplay.GaussianRadius = 0.5

Show Results

To view the result variables, first make sure that "fluid_1.pvd" is highlighted in the "Pipeline Browser", then select a variable in the variable-selection combo box (see the image below). For example, choose "density". To view the time evolution, press the play button (also shown below).

image

API

Modules = [TrixiParticles]
Pages = map(file -> joinpath("io", file), readdir(joinpath("..", "src", "io")))