[MaCl]

Hello everyone,

I’m new to OpenFOAM and currently have an issue modeling power generation in a reactor using the scalarSemiImplicitSource function.

The simulation is run in steady state with the internal values for both solid and fluid set to 1023 K. Then, in each reactor pin I apply the Su & SpX terms of the scalar semi-implicit source, where I vary Su between 1e7 W/m³ and 1 W/m³, while SpX follows a cosine function as shown in the figure below.

Each fuel pin has a volume of 0.0275 m³, and for the 1:12 reactor model I’m using, the total volume is 0.687 m³.

The result in the 1:12 reactor simulation shows that the temperature keeps rising above 2000 K, whereas in the single-pin simulation the temperature is around 1400 K. Should I reduce the Su value so that the temperature in the 1:12 model can drop to around 1400 K?

However, I don’t think that’s the real solution, because even when I try with Su = 1 W/m³, the temperature still blows up to 2000 K. It seems that this is influenced by the SpX value?

[HLX（CDUT）]

Originally Posted by MaCl

Hello everyone,

I’m new to OpenFOAM and currently have an issue modeling power generation in a reactor using the scalarSemiImplicitSource function.

The simulation is run in steady state with the internal values for both solid and fluid set to 1023 K. Then, in each reactor pin I apply the Su & SpX terms of the scalar semi-implicit source, where I vary Su between 1e7 W/m³ and 1 W/m³, while SpX follows a cosine function as shown in the figure below.

Each fuel pin has a volume of 0.0275 m³, and for the 1:12 reactor model I’m using, the total volume is 0.687 m³.

The result in the 1:12 reactor simulation shows that the temperature keeps rising above 2000 K, whereas in the single-pin simulation the temperature is around 1400 K. Should I reduce the Su value so that the temperature in the 1:12 model can drop to around 1400 K?

However, I don’t think that’s the real solution, because even when I try with Su = 1 W/m³, the temperature still blows up to 2000 K. It seems that this is influenced by the SpX value?

Actually, you can calculate it (heat transfer equation)