Implementing JJ BC for solid phase and NSW for gas phase

[mohsenclick]

@cgw HI Charles,
YOU are right. I am trying to implement the NSW for the fluid phase and JJ BC for the solid phase. I used the Mixed wall type in GUI as you mentioned. However, the *.mfix file is like the 2nd format. According to the:

Wall BC for TFM simulation

I did not find any issue. Please see the definition of wall boundary condition (bold sentence), that is why you see all wall boundary conditions are set to PSW no matter what boundary condition you chose. |FREE_SLIP_WALL| |Velocity gradients at the wall vanish. If BC_JJ_PS is equal to 1, the Johnson-Jackson boundary condition is used for solids. A FSW is equivalent to using a PSW with hw=0.| |NO_SLIP_WALL| |All components of the velocity vanish at the wall. If BC_JJ_PS is equal to 1, …

, it should be correct, but I think the BC does not implement correctly. What do you think?
Also, when we use the JJ BC for the solid phase, and the energy equation is not solved, the following terms should not be used in the solution:
bc_c_t_g(3) = 0.0 <
bc_c_t_s(3,1) = 0.0 <
bc_hw_t_g(3) = 0.0 <
bc_hw_t_s(3,1) = 0.0 <
bc_thetaw_m(3,1) = 0.0 <

                             >  bc_c_theta_m(3,1) = 0.0
                             >  bc_hw_s(3,1) = 0.0

6.mfx (18.9 KB)

Thanks,
Mohsen