Meshing questions in SC

[garcs2]

Hello,

I am attempting to transfer the solved power density from my heat conduction model to a subchannel model, using a MultiAppProjectionTransfer. I'm finding that I'm getting the following error:

libMesh terminating:
Element type not supported:EDGE2
Stack frames: 15
0: libMesh::print_trace(std::ostream&)
1: libMesh::MacroFunctions::report_error(char const*, int, char const*, char const*, std::ostream&)
2: libMesh::QGauss::init_2D(libMesh::ElemType, unsigned int)
3: libMesh::FE<2u, (libMesh::FEFamily)0>::reinit(libMesh::Elem const*, std::vector<libMesh::Point, std::allocator<libMesh::Point> > const*, std::vector<double, std::allocator<double> > const*)
4: MultiAppProjectionTransfer::execute()
5: FEProblemBase::execMultiAppTransfers(MooseEnumItem, Transfer::DIRECTION)
6: FEProblemBase::execMultiApps(MooseEnumItem, bool)
7: FixedPointSolve::solveStep(double&, double&, std::set<unsigned long, std::less<unsigned long>, std::allocator<unsigned long> > const&)
8: FixedPointSolve::solve()
9: SteadyBase::execute()
10: MooseApp::executeExecutioner()
11: MooseApp::run()
12: /home/garcsamu/moose_exec/snapbench/snapbench-opt(+0x50e8) [0x5555555590e8]
13: __libc_start_main
14: /home/garcsamu/moose_exec/snapbench/snapbench-opt(+0x538a) [0x55555555938a]
[0] ../src/quadrature/quadrature_gauss_2D.C, line 1300, compiled Mar  4 2025 at 21:19:48

Abort(1) on node 0 (rank 0 in comm 0): application called MPI_Abort(MPI_COMM_WORLD, 1) - process 0
libMesh terminating:
Element type not supported:EDGE2
Stack frames: 15
0: libMesh::print_trace(std::ostream&)
1: libMesh::MacroFunctions::report_error(char const*, int, char const*, char const*, std::ostream&)
2: libMesh::QGauss::init_2D(libMesh::ElemType, unsigned int)
3: libMesh::FE<2u, (libMesh::FEFamily)0>::reinit(libMesh::Elem const*, std::vector<libMesh::Point, std::allocator<libMesh::Point> > const*, std::vector<double, std::allocator<double> > const*)
4: MultiAppProjectionTransfer::execute()
5: FEProblemBase::execMultiAppTransfers(MooseEnumItem, Transfer::DIRECTION)
6: FEProblemBase::execMultiApps(MooseEnumItem, bool)
7: FixedPointSolve::solveStep(double&, double&, std::set<unsigned long, std::less<unsigned long>, std::allocator<unsigned long> > const&)
8: FixedPointSolve::solve()
9: SteadyBase::execute()
10: MooseApp::executeExecutioner()
11: MooseApp::run()
12: /home/garcsamu/moose_exec/snapbench/snapbench-opt(+0x50e8) [0x5555555590e8]
13: __libc_start_main
14: /home/garcsamu/moose_exec/snapbench/snapbench-opt(+0x538a) [0x55555555938a]
[0] ../src/quadrature/quadrature_gauss_2D.C, line 1300, compiled Mar  4 2025 at 21:19:48

This suggests to me that this is an issue with the meshing produced by SCMTriSubChannelMeshGenerator. This would make sense to me given that when I view standalone subchannel solutions in Paraview (without using the detailed mesh generator) all I see are thin lines for where the pins and subchannels are located. Would the suggestion then be that I should instead create a user object that records the average power density along the length of the fuel, and send that to the fuel pins in subchannel? I just would like to make sure that there's not a simpler solution or if there's an example that I haven't seen that has gone through a similar multiapps transfer.

[GiudGiud]

Hello

You should use the transfers used in the examples of SCM
Every other transfer is untested in combination with subchannel at this point

I ll make sure the MultiAppGeneralFieldTransfers work. ProjectionTransfer probably will be much later at this point

Guillaume

[kyriv1980]

An example of power transfer would look like this.

[q_prime] # send heat flux from /BISON/heatConduction to subchannel
  type = MultiAppGeneralFieldNearestLocationTransfer
  from_multi_app = sub
  variable = q_prime
  source_variable = q_prime
  from_boundaries = right
  execute_on = 'timestep_end'
[]

Bare in mind that SCM uses linear heat rate W/m not density.