I'm run a richard's case.
And my condition involve recharge(some >=0; some <0). [b](* ps: Because this area didn't have enough well data, so the negative recharge is replace well) [/b]
And my richard's parameters are following :
Property Sr alpha n
Gravel 0.0349 5.79 2.2355
Sand 0.045 14.5 2.68
Clay 0.068 0.8 1.09
Actually i'm not sure the gravel parameters(right/wrong?)
(reference:https://www.researchgate.net/publication/270118948_Evaluation_of_Hydraulic_Parameters_Obtained_by_Different_Measurement_Methods_for_Heterogeneous_Gravel_Soil)
[b]My problem is :
I set two layer for this case(aquifer, aquitard), and get convergence result.
However I discrete first layer to 10 layer(parameter be the same), get a result that no convergence(very high head).[/b]
Is my richard's parameters set wrong?
In my comprehension, the higher 'alpha' will cause water drain easy, so the recharge water will easy to drain to aquitard layer which numerical simulation easy divergence part,
therefore I try to let all layer 'alpha' be 0.0008, and get convergence. But it's very unreasonable.
thanks