Problem regarding the running of the immersed body.

[APTX-4869-core]

Hello, Incompact3d team. After yesterday's question, I encounter another problem.
I'm trying to run the examples/Sphere case, and I have installed the toolbox successfully with python-3.8. But I didn't find the line must be corrected in
/home/username/anaconda3/envs/x3d/lib/python3.12/site-packages/xcompact3d_toolbox/sandbox.py

So I ignored it and went on. But when I use python3 generate_initial_conditions.py, it reported the following error:

---

Traceback (most recent call last): File "generate_initial_conditions.py", line 11, in <module> prm = x3d.Parameters(loadfile='input.i3d') File "/home/.local/lib/python3.8/site-packages/xcompact3d_toolbox/parameters.py", line 910, in __init__ self.load(raise_warning=raise_warning) File "/home/.local/lib/python3.8/site-packages/xcompact3d_toolbox/parameters.py", line 1261, in load self.from_file(*arg, **kwarg) File "/home/.local/lib/python3.8/site-packages/xcompact3d_toolbox/parameters.py", line 1257, in from_file self.set(raise_warning=raise_warning, **dictionary) File "/home/.local/lib/python3.8/site-packages/xcompact3d_toolbox/parameters.py", line 1171, in set raise KeyError(msg) KeyError: 'nvol is not a valid parameter'

---

If I commented the "nvol" in the input.i3d, it would report the parameters following "i2dsim"(including i2dsim) were all invalid paramenters. So I commented them all, and it reported:

---

epsi Use of not exact is_closed check. This check can lead to misleading results. You could try to use exact=True. See: - false positive: https://github.com/wolph/numpy-stl/issues/198 - false negative: https://github.com/wolph/numpy-stl/pull/213 Use of not exact is_closed check. This check can lead to misleading results. You could try to use exact=True. See: - false positive: https://github.com/wolph/numpy-stl/issues/198 - false negative: https://github.com/wolph/numpy-stl/pull/213 Inflow Plane for Streamwise Velocity Inflow Plane for Vertical Velocity Inflow Plane for Spanwise Velocity Initial Condition for Streamwise Velocity Initial Condition for Vertical Velocity Initial Condition for Spanwise Velocity

---

Anyway it seemed the initial BC was set, but after I ran it, the result would get wrong(the field would get inf).

How can I solve this problem?
Any advice is appreciable and thanks in advance!

Regards,
Arashi

[mathrack]

Hello,

Have you looked at the readme in the examples/Sphere folder ?

Cheers,
Mathrack

[APTX-4869-core]

Hello mathrack!

1. Yes, I did follow the instructions in the README file.
   What's more, besides the instructions, the "iibm = 2" should be changed into "iibm = 1", otherwise it would diverged among few time steps. But after these were all done, the result still diverged as the forces0010000.dat read:
   

2. And I have another question by the way. How can I read the binary files the Channel cases outputed, such as the files in statics/ folder? I tried to open them via Paraview using the "Tecplot Table Reader", but it showed blank in the window.

Thanks in advance!

Regards,
Arashi

[mathrack]

Thank you for the update. I will look at it when I have spare time. With the provided information, it looks like a bug in the toolbox and not in Incompact3d.

For reading 3D binary files from python, I often use the following :

import numpy as np
import argparse

...

default_nx = 1
default_ny = 2
default_nz = 3

...

parser = argparse.ArgumentParser(description='Read and filter binary stats')
parser.add_argument('--nx', type=int, default=default_nx, help='number of cells in x')
parser.add_argument('--ny', type=int, default=default_ny, help='number of cells in y')
parser.add_argument('--nz', type=int, default=default_nz, help='number of cells in z')

...

#
# Read a binary file, store it in a 1D array
def read_bin_file(file, dtype=np.float64):
    return np.fromfile(file, dtype)

#
# Read a binary file, store it in a 3D array
def read_3d_file(file, dtype=np.float64, nx=default_nx, ny=default_ny, nz=default_nz):
    data1d = read_bin_file(file, dtype)
    return data1d.reshape((nx, ny, nz), order='F')