gautelinga / bernaise Goto Github PK

Create a small module where dimensionless parameters can be converted into simulation units.

Constant velocity at inlet and outlet.

Hello, I am trying to run some of the included cases (which run fine) but have a problem visualising the results with Paraview (and with the included postprocessing script).

The xdmf file contains lines like

<DataItem Dimensions="2145 2" Format="HDF">V_from_tstep_0.h5</DataItem>

for every reference to DataItem, while a reference xdmf file should contain a semicolon and then a path, like this example (from https://github.com/MattH688/BERNAISE/blob/master/meshes/mf_flowfocus.xdmf#L6)

<DataItem DataType="Float" Dimensions="9617 3" Format="HDF" Precision="8">mf_flowfocus.h5:/data0</DataItem>

Opening with paraview thus leads to errors like

XDMF Error in /build/paraview/src/ParaView-v5.10.1/VTK/ThirdParty/xdmf2/vtkxdmf2/libsrc/XdmfValuesHDF.cxx line 76 (Can't Open Dataset V_from_tstep_0.h5)
XDMF Error in /build/paraview/src/ParaView-v5.10.1/VTK/ThirdParty/xdmf2/vtkxdmf2/libsrc/XdmfDataItem.cxx line 534 (Reading Values Failed)
XDMF Error in /build/paraview/src/ParaView-v5.10.1/VTK/ThirdParty/xdmf2/vtkxdmf2/libsrc/XdmfGrid.cxx line 755 (Error in Update() of Topology)
(   4.603s) [paraview        ]      vtkXdmfReader.cxx:488    ERR| vtkXdmfReader (0x55b7cb7e40a0): Failed to read data.
(   4.675s) [paraview        ]       vtkExecutive.cxx:752    ERR| vtkPVCompositeDataPipeline (0x55b7cf2f9920): Algorithm vtkFileSeriesReader(0x55b7d40b83a0) returned failure for request: vtkInformation (0x55b7cf36bd20)

Postprocess.py

Errors related to the same problem:

python postprocess.py folder=results_charged_droplet/1/ method=plot save=True                Tue 04 Oct 2022 07:56:27 PM UTC
BERNAISE: Post-processing tool
Sought fields: All
Opening folder: results_charged_droplet/1/
Traceback (most recent call last):
  File "/home/<snip>/BERNAISE/postprocess.py", line 157, in <module>
    main()
  File "/home/<snip>/BERNAISE/postprocess.py", line 145, in main
    ts = TimeSeries(folder, sought_fields=sought_fields)
  File "/home/<snip>/BERNAISE/utilities/TimeSeries.py", line 62, in __init__
    self._load_timeseries(sought_fields)
  File "/home/<snip>/BERNAISE/utilities/TimeSeries.py", line 136, in _load_timeseries
    = parse_xdmf(xml_file, get_mesh_address=True)
  File "/home/<snip>/BERNAISE/common/io.py", line 241, in parse_xdmf
    dset_address = prop[0].text.split(":")[1]
IndexError: list index out of range

Software

how do i make the mesh for rounde barbell periodic, and way is Snausen/Snoevsen periodic?

After a simulation is killed, there are no checkpoint files in the Checkpoints folder.

The basic solver Crashes when enable_NS, enable_EC and enable_PF is set to False (One at a time) due to variabols that are not defined! see setup function especilaly line 119 -138 (solvers/basic.py)
Futher there is a problem in line 284-294 (solvers/basic.py) when calling the solver

We need to write a post-processing tool to analyze simulations after they have been run. Therefore we should write a small piece of code to read meshes and data from xdmf files and analyze the simulations a posteriori.

(I have some code for this lying around, so I will assign myself to this issue.)