Periodic constant heat flux boundary conditions

18 views

Skip to first unread message

MWTWM M

unread,

Mar 5, 2026, 9:02:17 AM (14 days ago) Mar 5

to Nek5000

Dear community,

I encountered some difficulties when computing a case with periodic constant heat flux boundary conditions. Following the nondimensionalization procedure shown in the Quick Start, I derived the following relations (the superscript “+” denotes nondimensional quantities). Could you please clarify how the nondimensional thermal conductivity, Λ+, is treated in the program? In particular, can the code shown in the red rectangle, param(8), be interpreted as 1/Pe rather than simply the thermal conductivity?

In addition, I modified the turbChannel case based on several references in order to avoid the issue of “unbounded temperature increase” that occurs in periodic constant heat-flux boundary simulations. Specifically:

0205: constant wall temperature boundary (temp = 1.0), qvol = vx(i,j,k,e) / 2330.0 (Pe = 2330)

0206: constant wall temperature boundary (temp = 1.0), qvol = vx(i,j,k,e)

0207: constant heat-flux boundary (flux = 1.0, temp = 0.0), qvol = vx(i,j,k,e)

Although the numerical values of the temperature field differ significantly in the ParaView visualization, they become close to each other after being processed using t_tau. Are such computational results reasonable?

Any suggestions and criticisms are appreciated.

Student W. T. Ma

Reply all

Reply to author

Forward

Message has been deleted

0 new messages