Newer
Older
% Generated by roxygen2: do not edit by hand
% Please edit documentation in R/Rphree_ReactTrans.R
\name{AdvectionPQC}
\alias{AdvectionPQC}
\title{Simple explicit 1D forward-Euler advection, multispecies}
\usage{
AdvectionPQC(conc, inflow = rep(0, ncol(conc)), dx, dt, U)
\item{conc}{either a vector or a matrix containing the chemical
state to be updated}
\item{inflow}{vector of concentrations entering the inlet}
\item{dx}{the grid spacing in m}
\item{dt}{the required time step in s}
\item{U}{Darcy velocity in m^3/s}
}
\value{
a matrix containing the state after the advection
}
\description{
This function takes the current state of a chemical system in form
of a matrix, a vector of concentrations representing boundary
conditions (injection), the required time step (dt/s), grid spcing
(dx/m) and U (Darcy velocity/m^3/s). Note that it is
responsibility of user to check that the required dt does not
trespass the CFL condition.