hydro:brooks-corey
Unterschiede
Hier werden die Unterschiede zwischen zwei Versionen angezeigt.
Beide Seiten der vorigen RevisionVorhergehende ÜberarbeitungNächste Überarbeitung | Vorhergehende Überarbeitung | ||
hydro:brooks-corey [2017/06/16 18:56] – ckuells | hydro:brooks-corey [2024/04/10 10:02] (aktuell) – Externe Bearbeitung 127.0.0.1 | ||
---|---|---|---|
Zeile 1: | Zeile 1: | ||
- | ~~NOTRANS~~ | + | |
====== Brooks-Corey Funktion ====== | ====== Brooks-Corey Funktion ====== | ||
Die Brooks Corey Funktion stellt eine Beziehung zwischen Bodenfeuchte oder Sättigung einerseits und der Saugspannung andererseits her: | Die Brooks Corey Funktion stellt eine Beziehung zwischen Bodenfeuchte oder Sättigung einerseits und der Saugspannung andererseits her: | ||
- | $$ | + | $$ \frac{\theta-\theta_r}{\theta_s-\theta_r} = {\left[\frac{\psi_0}{\psi}\right]}^\lambda$$ |
- | \frac{\theta-\theta_r}{\theta_s-\theta_r} = {\left[\frac{\psi_0}{\psi}\right]}^\lambda | + | |
- | mit den Parametern: | + | mit den Parametern |
- | + | ||
- | \begin{tabular}{ll} | + | |
- | | + | |
- | | + | |
- | | + | |
- | | + | |
- | | + | |
- | \end{tabular} | + | |
- | $$ | + | |
Die ungesättigte Leitfähigkeit wird nach Brooks-Corey durch folgende Beziehung dargestellt: | Die ungesättigte Leitfähigkeit wird nach Brooks-Corey durch folgende Beziehung dargestellt: | ||
- | $$ | + | $$ \frac{K({\theta})}{K_s} = {\left[{\frac{\theta-\theta_r}{\theta_s-\theta_r}}\right]}^{3+2*\lambda}$$ |
- | \frac{K({\theta})}{K_s} = {\left[{\frac{\theta-\theta_r}{\theta_s-\theta_r}}\right]}^{3+2*\lambda} | + | |
+ | mit den Parametern $K_\theta$ ungesättigte Bodenfeuchte, | ||
+ | |||
+ | <code S+ | Brooks-Corey.R> | ||
+ | |||
+ | # Brooks-Corey (1966) | ||
+ | brookscorey <- function(theta, | ||
+ | # Standard values | ||
+ | ths <- seq(0.05, | ||
+ | trs <- 0.035 | ||
+ | tms <- 0.45 # maximale Feuchte | ||
+ | kss <- 1E-5 # hydraul. Leitfähigkeit in m/s | ||
+ | ips <- 0.2 # pore size index | ||
+ | # Define values | ||
+ | if (is.null(theta)) { | ||
+ | theta <- ths # Default Bodenfeuchte | ||
+ | } | ||
+ | if (is.null(tred)) | ||
+ | tred <- trs # Default Restfeuchte | ||
+ | } | ||
+ | if (is.null(tmax)) | ||
+ | tmax <- tms # Default max. Feuchte | ||
+ | } | ||
+ | if (is.null(ksat)) | ||
+ | ksat <- kss # Default Durchlaessigkeit | ||
+ | } | ||
+ | if (is.null(pi)) | ||
+ | pi <- ips # Default Porenindex | ||
+ | } | ||
+ | # Calculate unsaturated hydraulic conductivity | ||
+ | Se <- (theta-tred)/ | ||
+ | if ((Se < 0) || (Se > 1)) { | ||
+ | | ||
+ | } | ||
+ | else{ | ||
+ | kunsat <- Se^(2+3*pi) | ||
+ | if (is.null(kunsat)) { | ||
+ | |||
+ | result <- "Error in calculation of unsaturated conductivity." | ||
+ | } | ||
+ | else{ | ||
+ | ratio <- Se^(2+3*pi) | ||
+ | kunsat <- ratio * ksat | ||
+ | result <- list(Se, | ||
+ | } | ||
+ | return(result) | ||
+ | } | ||
+ | |||
+ | #Test | ||
+ | unsat <- brookscorey(0.25, | ||
- | mit den Parametern: | + | # Values |
+ | ths <- seq(0.05, | ||
+ | trs <- 0.045 | ||
+ | tms <- 0.45 # maximale Feuchte | ||
+ | kss <- 1E-5 # hydraul. Leitfähigkeit in m/s | ||
+ | ips <- 0.2 # pore size index | ||
- | \begin{tabular}{ll} | + | unsat <- brookscorey(ths, |
- | | + | Se <- unsat[[1]] |
- | K_s & gesättigte Bodenfeuchte \\ | + | ku <- unsat[[3]] |
- | \lambda & Porengrößenindex \\ | + | |
- | \end{tabular} | + | |
- | $$ | + | |
+ | plot(Se,ku) | ||
+ | </ |
/usr/www/users/uhydro/doku/data/attic/hydro/brooks-corey.1497632161.txt.gz · Zuletzt geändert: 2024/04/10 10:12 (Externe Bearbeitung)