Spherical shell geometry

[Bhargav Mantravadi]

Dear Team,

I hope you are doing well. I recently began using Nek5000, and I am working with spherical shell geometry. I am currently using the MATLAB code that Prof. Fischer provided previously in this group. 

After generating the .rea file with the MATLAB script, I run reatore2 to obtain the .re2 file and then compile the code. However, when I run the command nekbmpi test 6, I get the error shown below. Could you please advise on how to resolve this issue?

Also, do you recommend any alternative approaches for building a spherical shell mesh for Nek5000 simulations?

Thanks,

Bhargav

--------------------------------------

 Number of MPI ranks :           6

 Reading test.par
general = [(null)]
general:stopat = [numSteps]
general:numsteps = [10]
general:dt = [5.00E-02]
general:variabledt = [yes]
general:targetcfl = [0.5]
general:writecontrol = [timeStep]
general:writeinterval = [1]
velocity = [(null)]
velocity:viscosity = [1e-03]
velocity:density = [1.0]
pressure = [(null)]
pressure:solver = [GMRES]
pressure:residualtol = [1e-06]
temperature = [(null)]
temperature:conductivity = [1e-03]
temperature:rhocp = [1.0]
 reading test.re2
 hdr:#v002     5120  3     5120 this is the hdr
 partioning elements to MPI ranks
 reading mesh
   readp_re2_mesh:pack/cr/unpack : 0.13E-02 0.19E-03 0.47E-04
   readp_re2_mesh:byte_read_mpi  : 0.69E-01
 reading boundary faces                 1024 for ifield   1
   readp_re2_bc:pack/cr/unpack : 0.54E-05 0.62E-03 0.24E-04
   readp_re2_bc:byte_read_mpi  : 0.45E-01
ERROR: Error reading .re2 boundary data  ierr=  6

[YuHsiang Lan]

Hi Bhargav,

Have you specified the number of BC in par file?

  https://nek5000.github.io/NekDoc/faq.html#troubleshooting

ChatGPT response (https://chatgpt.com/g/g-68a731e21f7c8191bcc684cc265a1d79-nek5000-assistant)

  https://chatgpt.com/share/6929cbdf-ca24-8009-aae9-895763ad16a5

See the 4.1 for prenek that makes spherical shell. It might not be up-to-date but the functionality is still there.

Hope this helps,

Yu-Hsiang

--

[Bhargav Mantravadi]

Hi Yu-Hsiang,

Thank you very much for your response. I checked the BCs, and the .rea file generated from the MATLAB script does not have thermal BCs. I have modified the script file, and with the updated .rea file, I run reatore2, genmap, and compile the code. I can run the simulation now.

Thank you for the inputs. 

Regards,

Bhargav

--
You received this message because you are subscribed to the Google Groups "Nek5000" group.
To unsubscribe from this group and stop receiving emails from it, send an email to nek5000+u...@googlegroups.com.
To view this discussion visit https://groups.google.com/d/msgid/nek5000/0940c70b-76d2-48a6-a720-69fb32f8bb23n%40googlegroups.com.