Now I want to impose a value for scalar transport at inlet face by cs_user_boundary_conditions.f90, not GUI.
Can we impose a value of scalar transport at inlet by cs_user_boundary_conditions.f90?
At the same time I am imposing pressure value following the code.
However, errors happend when computing.
Just wIth the code for pressure, it works, but with the code of scalar it does not work.
Can we impose the specific value of scalar transpot at inlet by usersubroutine?
Once I succeed, I will edit the code to impose time-dependent value for scalar transport.
Code: Select all
call getfbr('inlet', nlelt, lstelt)
do ilelt = 1, nlelt
ifac = lstelt(ilelt)
itypfb(ifac) = ientre
!Dirichlet B.C for pressure:
icodcl(ifac,ipr) = 1
rcodcl(ifac,ipr,1) = 101325.633
!Dirichlet B.C for scalar:
icodcl(ifac,isca(1))=1
rcodcl(ifac,isca(1),1)=10 ! desired value
!Neumann B.C for velocity:
icodcl(ifac,iu)= 3
icodcl(ifac,iv)= 3
icodcl(ifac,iw)= 3
enddo
What is "i"? and Can I find "i" in setup log?
Would someone help me?
Tsubasa