TLDR: 2D RVE (brick matrix + circular air holes), ~650 K, radiation matters. S2S doesn't support periodic BCs, so I used adiabatic top/bottom + fixed T on left/right instead. The half/quarter circles cut at domain boundaries ended up in a default adiabatic zone (not in my named selections), and I'm unsure if radiation is handled correctly there. Main question: is this adiabatic corner treatment acceptable for effective conductivity, or should I switch to DO (which supports periodic)?
Edit: i added new images to the comments
Full details: I'm running a 2D steady-state conjugate heat transfer analysis on an RVE (representative volume element) for a thermal energy storage material for my master's thesis. The geometry is a rectangular brick matrix (alumina firebrick, k=1.8 W/m·K) with circular air holes arranged in a regular pattern. Operating temperature is around 650 K, so radiation is significant.
Setup:
Fluent 2022 R2, 2D planar, pressure-based, steady
Energy equation on, Viscous: Laminar
Radiation: S2S (surface-to-surface)
Left edge: fixed T = 655 K, Right edge: fixed T = 645 K
Top/bottom edges: adiabatic (I originally wanted periodic, but S2S is incompatible with periodic BCs, so I switched to adiabatic)
Air holes are Fluid zones (no gravity, so essentially stagnant — pure conduction + radiation)
Surface emissivity = 0.85 for all walls
The problem:
The circular holes are cut at the domain boundaries (half-circles at edges, quarter-circles at corners). During meshing, some of these boundary arc segments were not assigned to my named selections, so Fluent grouped them into a default zone ("part_3-air") with a shadow pair, set as adiabatic.
In my temperature contour, most of the domain looks physically correct (smooth gradient from hot to cold side, isotherms bending around the low-conductivity holes). However, I'm concerned about how radiation is being handled at these corner/edge arc zones, since they are not part of a proper periodic or continuous boundary.
My questions:
Since S2S doesn't support periodic BCs, what is the recommended way to handle radiation in an RVE where the physical situation should be periodic?
How should the boundary-cut arc segments (half/quarter circles at domain edges) be treated for radiation? Should they participate in S2S or be excluded?
Would switching to the DO (Discrete Ordinates) model — which supports periodic BCs — be more appropriate for this type of RVE radiation problem? Any drawbacks for a 2D case with many radiating surfaces?
Is the adiabatic assumption at these corner zones acceptable given that the net heat flux through them appears negligible, or does it introduce significant error in the effective conductivity calculation?
I can share screenshots of the mesh, contour, and boundary condition setup if helpful. Thanks in advance.