Upon adding active control to my simulation case, the inlet speed kept going to zero resulting in a flash-back.
I found out by looking at PeleLM.H that m_ctrl_velMax is set to zero by default, resulting in Vnew always being equal to zero:
225 Vnew = std::max(Vnew, 0.0);
227 Vnew = std::min(Vnew, m_ctrl_velMax);
I don't have an issue anymore with using active control as I know I can just supply a value to make my case work, but perhaps something should be changed either in the TripleFlame example:
active_control.velMax = 1.0 # Optional: limit inlet velocity
Or in the code by adding a condition for line 227 of PeleLMFlowControler.cpp (slow ?) or setting an arbitrary large default value in PeleLM.H (not good practice ?).
I would be appy to submit a pull request with the corresponding changes.
Upon adding active control to my simulation case, the inlet speed kept going to zero resulting in a flash-back.
I found out by looking at
PeleLM.Hthatm_ctrl_velMaxis set to zero by default, resulting inVnewalways being equal to zero:I don't have an issue anymore with using active control as I know I can just supply a value to make my case work, but perhaps something should be changed either in the TripleFlame example:
active_control.velMax = 1.0 # Optional: limit inlet velocityOr in the code by adding a condition for line 227 of
PeleLMFlowControler.cpp(slow ?) or setting an arbitrary large default value inPeleLM.H(not good practice ?).I would be appy to submit a pull request with the corresponding changes.