HP1 SEASONCHAIN
Posted: Thu Aug 01, 2024 3:40 pm
I am trying to repeat HP1 SEASONCHAIN example in Hydrus 1D with HP1. Yet I need to account for stratified medium, say two layers. Each strata has its own initial conditions, sorption properties and first order decay reaction coefficients.
I think that first order decay reaction coefficients, variable per layer, can be accounted for by using two Kinetics blocks with appropriate reference to cells say Kinetics 1-10 and Kinetics 11-100. Each block calls RATEs with parameters appropriate for cells 1-10 and 11-100.
However, I am not sure how to define sorption. In SEASONCHAIN there is only surface with sites Sora, Sorb and Sorc. Should I define (accordingly to Kinetics) Surface 1-10 and Surface 11-100 and equilibrate them with appropriate initial solution?
If this is the case, can sites Sora, Sorb and Sorc be the same in each Surface 1-10 and Surface 11-100? Or each surface should have its “own” Sor(abc) sites?
I think that first order decay reaction coefficients, variable per layer, can be accounted for by using two Kinetics blocks with appropriate reference to cells say Kinetics 1-10 and Kinetics 11-100. Each block calls RATEs with parameters appropriate for cells 1-10 and 11-100.
However, I am not sure how to define sorption. In SEASONCHAIN there is only surface with sites Sora, Sorb and Sorc. Should I define (accordingly to Kinetics) Surface 1-10 and Surface 11-100 and equilibrate them with appropriate initial solution?
If this is the case, can sites Sora, Sorb and Sorc be the same in each Surface 1-10 and Surface 11-100? Or each surface should have its “own” Sor(abc) sites?