subroutine problem
Posted: Wed May 08, 2019 6:40 pm
Hello
i have a problem in the the boundary conditions subroutine cs_user_boundary_conditions.f90, when i set the temperature change by time in the code as :
call getfbr('16',nlelt,lstelt)
if (ttcabs.lt.3.8d0) then
do ilelt = 16, nlelt
ifac = lstelt(ilelt)
rcodcl(ifac,isca(iscalt),16) = 20.d0 + 100.d0*ttcabs
enddo
else
do ilelt = 16 , nlelt
ifac = lstelt(ilelt)
rcodcl(ifac,isca(iscalt),16) = 400.d0
enddo
endif
call getfbr('18',nlelt,lstelt)
if (ttcabs.lt.3.8d0) then
do ilelt = 18, nlelt
ifac = lstelt(ilelt)
rcodcl(ifac,isca(iscalt),18) = 20.d0 + 100.d0*ttcabs
enddo
else
do ilelt = 18 , nlelt
ifac = lstelt(ilelt)
rcodcl(ifac,isca(iscalt),18) = 400.d0
enddo
endif
in two inlets, which are above the cylinder , i notice in paraview that they have same variation (in red )as the other two inlets which are in the side of cylinder) , as shown in the picture attached does that means the temperature variation by time are not working properly ?
Thank you
i have a problem in the the boundary conditions subroutine cs_user_boundary_conditions.f90, when i set the temperature change by time in the code as :
call getfbr('16',nlelt,lstelt)
if (ttcabs.lt.3.8d0) then
do ilelt = 16, nlelt
ifac = lstelt(ilelt)
rcodcl(ifac,isca(iscalt),16) = 20.d0 + 100.d0*ttcabs
enddo
else
do ilelt = 16 , nlelt
ifac = lstelt(ilelt)
rcodcl(ifac,isca(iscalt),16) = 400.d0
enddo
endif
call getfbr('18',nlelt,lstelt)
if (ttcabs.lt.3.8d0) then
do ilelt = 18, nlelt
ifac = lstelt(ilelt)
rcodcl(ifac,isca(iscalt),18) = 20.d0 + 100.d0*ttcabs
enddo
else
do ilelt = 18 , nlelt
ifac = lstelt(ilelt)
rcodcl(ifac,isca(iscalt),18) = 400.d0
enddo
endif
in two inlets, which are above the cylinder , i notice in paraview that they have same variation (in red )as the other two inlets which are in the side of cylinder) , as shown in the picture attached does that means the temperature variation by time are not working properly ?
Thank you