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 19:07] – 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: | + | mit den Parametern |
- | \begin{tabular}{ll} | + | <code S+ | Brooks-Corey.R> |
- | K_\theta & ungesättigte Bodenfeuchte \\ | + | |
- | K_s & gesättigte Bodenfeuchte \\ | + | |
- | \lambda & Porengrößenindex \\ | + | |
- | \end{tabular} | + | |
- | $$ | + | |
- | + | ||
- | <code S+ | Xinanjiang.R> | + | |
# Brooks-Corey (1966) | # Brooks-Corey (1966) | ||
- | th <- seq(0.0,0.45, | + | brookscorey <- function(theta, |
- | tr <- 0.05 # Restfeuchte | + | # Standard values |
- | tm <- 0.45 # maximale Feuchte | + | ths <- seq(0.05,0.4,0.005) |
- | Ks <- 1E-5 # hydraul. Leitfähigkeit in m/s | + | trs <- 0.035 |
- | ip <- 0.05 | + | 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, | ||
- | Se <- (th-tr)/(tm-tr) | + | # 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 | ||
- | ps0<- 0.15 # Lufteintrittspunkt: | + | unsat <- brookscorey(ths,trs,tms,kss,ips) |
- | psa<- seq(0,1E-4,5E-6)# Saugspannung in cm | + | Se <- unsat[[1]] |
+ | ku <- unsat[[3]] | ||
- | S <- (ps0/psa)^ip | + | plot(Se,ku) |
- | plot(S,psa) | ||
</ | </ |
/usr/www/users/uhydro/doku/data/attic/hydro/brooks-corey.1497632865.txt.gz · Zuletzt geändert: 2024/04/10 10:12 (Externe Bearbeitung)