Cyclic 'U' and Non-cyclic 'Scalar' on cyclic boundary
Dear FOAMers:
I am using OpenFOAM 7.0 to simulate the scalar transport in curved open channel, with interFoam and smagorinsky LES model. In order to get the fully developed flow field, I impose cyclic BC for U and p_rgh at INLET and OUTLET. For the variable scalar, I impose Dirichlet BC at INLET and Neumann BC at OUTLET. It turns out not working. OpenFOAM requires the same cyclic BC for scalar. BTW, the scalar is released in the form of point source near the INLET. How could I realize the above request?:D Any ideas are appreciated. Thanks in advance! Here I post my boundary file, U file and scalar file: Code:
FoamFile Code:
FoamFile Code:
FoamFile |
May be you can use the following appoach:
Calculate the flux (U and p) with your periodic b.c. but without the scaler. Uses scalaTransportFoam afterwards. You cannot get a transient result for the scalar, however. |
Hi, Pilz, thanks for your reply!
Your idea sounds reasonable, I gonna adopt the advice and apply it in my simulation. However, I am still wondering there are any simpler ways for cyclic cases. |
It is always valuable de-coupling effects in a simulation. This way, scalarTransportFoam is easier not more complicated. And the results are more reliable in comparison ton al-in-all simulation.
Cyclic b.c. often arise problems and lead to nonphysical results. It may be you have to look for another approach anyway. |
Cyclic 'U' and Non-cyclic 'Scalar' on cyclic boundary
Dear hhu lulu,
Were you able to solve your problem? I faced the same problem when I tried to simulate 2D periodic hills with passive scalar transport. Using scalarTransportFoam solver, as mentioned by piuy58, seems more perspective. |
Have you tried to add a heat sink using fvOptions/fvModels? In the paste several channel flows with similar scalar conditions were performed.
|
Quote:
|
Quote:
I have been searching but I only found constant wall heat flux, I thought I saw something for fixed scalar value. I have some ideas to deal with this problem
Regards |
Quote:
Sorry for my slow reply. I find a workaround with my issue. The periodic B.C. has still been used for the scalar in my case. To prevent the scalar from recycling, the scalarFixedValueConstraint has been used to set the value of scalar to 0 at the end of the computational domain. That is, the scalar is released from the source and it vanishes near the outlet, which would not influence the scalar transport in a single period. I don't know whether this method matches your problem. Maybe you can elaborate on your cases. Thanks! |
All times are GMT -4. The time now is 13:26. |