schwaral
Subscriber
Could it even possibly be the mesh quality, if it converges with a Diffusion Doefficient (D_c) =0? For small and large D_c the UDS values go towards infinite. I did find people are using linearization of source term. Might this be a problem? When I disable the FLUX Term for both UDS (I think the D_c should not be used then??) the solution diverges as well. This last point confuses me a lot.nnMy Source Term for UDS1/2 are: nC is scalar from UDS 1, M is scalar of UDS2 --- ast. adyn, gst, gdy, C_sat, M_sat are coefficients --- c_vel = velocity.nnsource (C)= -c_density * ((ast + adyn * c_vel) * (1 - M / M_sat) * C + (gst + gdyn * c_vel) * (1 - C / C_sat) * M);ndS[eqn] = -c_density * ((ast + adyn * c_vel) * (1 - M / M_sat) - (gst + gdyn * c_vel)*M/C_sat);nnsource (M)= c_density * ((ast + adyn * c_vel) * (1 - M / M_sat) * C + (gst + gdyn * c_vel) * (1 - C / C_sat) * M);ndS[eqn2] = -c_density*((gst + gdyn * c_vel) * (1 - C / C_sat)+(ast + adyn * c_vel)*C/M_sat);n