Reply To: Advection-Diffusion Dirichlet Temperature BC
OpenLB – Open Source Lattice Boltzmann Code › Forums › on OpenLB › General Topics › Advection-Diffusion Dirichlet Temperature BC › Reply To: Advection-Diffusion Dirichlet Temperature BC
Dear Mathias,
At the inlet I define bulkdynamics, use setAdvectionDiffusionTemperatureBoundary and finally ADlattice.defineRho(superGeometry,inlet,T_hot)
; where T_hot is an AnalyticalConst3D
. For the NS lattice, I set the velocity using setInterpolatedVelocityBoundary
.
Following your suggestion of resetting the temperature, I added the following to the function setBoundaryValues
:
void setBoundaryValues( ThermalUnitConverter<T, NSDESCRIPTOR, TDESCRIPTOR> const& converter,
SuperLattice3D<T, NSDESCRIPTOR>& NSlattice,
SuperLattice3D<T, TDESCRIPTOR>& ADlattice,
int iT, SuperGeometry3D<T>& superGeometry)
{
int iTmaxStart = converter.getLatticeTime( maxPhysT*0.2 );
int iTupdate = 100;
if ( iT%iTupdate == 0 && iT<= iTmaxStart ) {
AnalyticalConst3D<T,T> T_hot(Thot);
ADlattice.defineRho(superGeometry, 3 , T_hot);
}
}
With this routine I haven’t noticed any changes, I still see a decrease in temperature at the inlet from timestep 2 onwards. However, I may not be doing the reset appropriately.
Thank you for your time,
Ramiro