Issue importing Plot3D grid

[mburabay]

Hi,

I am using the following in the prep-grid file:

---

fileName = "StructGridExp1.x"

grids = {}

grids[0] = importPlot3DGrid(fileName, 3, scale=1.0)

registerFluidGridArray{
grid = grids[0],
nib=2, njb=8, nkb=4,
fsTag="initial",
bcTags={north="inlet", top="symmetry", bottom="symmetry", east="symmetry", west="symmetry"}
}

---

However, I keep getting the error at the prep-grid stage:

---

lmr.globalconfig.FlowSolverException@./commands/prepgrids.d(127): StructGridExp_grid_v1.lua:7: ')' expected near '='

where line 7 is: grids[0] = importPlot3DGrid(fileName, 3, scale=1.0)

I am not sure why the mistake happens because I looked up this format in the grids-in-eilmer.adoc file

[mburabay]

I could not find an example that imports this mesh format. I know there are some for GridPro, but I use an 800 line script to prepare my mesh in Pointwise, which does not export in the GridPro format.

[rjgollan-on-github]

You are using a function: importPlot3DGrid. You cannot pass a keyword argument (scale=1.0) to a Lua function. Just pass a plain argument.

[mburabay]

Thank you, removing "scale=" fixed the error.

Now I have another error:

lmr.globalconfig.FlowSolverException@./commands/prepgrids.d(127): /myApps/EilmerFiles_v2/gdtkinst/lib/gridarray.lua:51: method 'get_type' is not callable (a nil value)

??:? int lmr.commands.prepgrids.main_(immutable(char)[][]) [0x55809845cd03]
??:? _Dmain [0x5580985bec52]

[mburabay]

The simulation now runs. The prep-grid file looks like this:

---

config.dimensions = 3

grids = {}

fileName = "StructGridExp.x"

grid1 = importPlot3DGrid(fileName, 3, 1.0)

registerFluidGridArray{
grid = grid1[1],
nib=4, njb=8, nkb=8,
fsTag="initial",
bcTags={east="symmetry", west="symmetry", top="symmetry", bottom="symmetry", south="inlet"}
}

registerFluidGridArray{
grid = grid1[2],
nib=4, njb=8, nkb=8,
fsTag="initial",
bcTags={east="wall", west="wall", top="wall", bottom="wall", north="outlet"}
}

identifyGridConnections()

---

The corresponding part of the prep-flow file is:

---

-- Set up flow state
inflow = FlowState:new{p=p_inf, T=T_inf, velx=u_inf, vely=0.0, velz=0.0, nuhat=nuhat_inf}

flowDict = {initial=inflow}
flowDict["inflow"] = inflow

bcDict = {inlet = InFlowBC_Supersonic:new{flowState=inflow},
wall = WallBC_NoSlip_FixedT:new{Twall=T_wall,group="wall"},
outlet=OutFlowBC_Simple:new{},
symmetry=WallBC_WithSlip:new{group="wall"}}

makeFluidBlocks(bcDict, flowDict)
mpiDistributeBlocks{ntasks=512}

---

However, the resulting solution is wrong:

I also tried using identifyBlockConnections(), but that didn't help. What could be the issue?

[rjgollan-on-github]

I would usually inspect the config file at this point to see that my boundary conditions are exactly as I expect them to be.

You'll need to identify your block numbers and your boundaries and inspect that your config file is correct in terms of boundary conditions.

Likely it is not. At that point, you can trace back to your BC specifications in your Lua input file.

(I am NOT suggesting that you hand edit the config file. But it's a good way to track back to your input problem.)