Skip to content

Commit

Permalink
Adding multi-aquifer well energy transport (MWE) package. Includes ne…
Browse files Browse the repository at this point in the history
…w autotest
  • Loading branch information
emorway-usgs committed Feb 9, 2024
1 parent a1415f2 commit 3a6b1ef
Show file tree
Hide file tree
Showing 22 changed files with 2,338 additions and 16 deletions.
570 changes: 570 additions & 0 deletions autotest/test_gwe_mwe_conduction.py

Large diffs are not rendered by default.

16 changes: 8 additions & 8 deletions doc/Common/gwe-mweobs.tex
Original file line number Diff line number Diff line change
@@ -1,13 +1,13 @@
% general APT observations
MWE & temperature & ifno or boundname & -- & Well temperature. If boundname is specified, boundname must be unique for each well. \\
MWE & temperature & mawno or boundname & -- & Well temperature. If boundname is specified, boundname must be unique for each well. \\
%flowjaface not included
MWE & storage & ifno or boundname & -- & Simulated energy storage flow rate for a well or group of wells. \\
MWE & constant & ifno or boundname & -- & Simulated energy constant-flow rate for a well or group of wells. \\
MWE & from-mvr & ifno or boundname & -- & Simulated energy inflow into a well or group of wells from the MVE package. Energy inflow is calculated as the product of provider temperature and the mover flow rate. \\
MWE & mwt & ifno or boundname & \texttt{iconn} or -- & Energy flow rate for a well or group of wells and its aquifer connection(s). If boundname is not specified for ID, then the simulated well-aquifer flow rate at a specific well connection is observed. In this case, ID2 must be specified and is the connection number \texttt{iconn} for well \texttt{ifno}. \\
MWE & storage & mawno or boundname & -- & Simulated energy storage flow rate for a well or group of wells. \\
MWE & constant & mawno or boundname & -- & Simulated energy constant-flow rate for a well or group of wells. \\
MWE & from-mvr & mawno or boundname & -- & Simulated energy inflow into a well or group of wells from the MVE package. Energy inflow is calculated as the product of provider temperature and the mover flow rate. \\
MWE & mwe & mawno or boundname & \texttt{iconn} or -- & Energy flow rate for a well or group of wells and its aquifer connection(s). If boundname is not specified for ID, then the simulated well-aquifer flow rate at a specific well connection is observed. In this case, ID2 must be specified and is the connection number \texttt{iconn} for well \texttt{mawno}. \\

% observations specific to the mwt package
MWE & rate & ifno or boundname & -- & Simulated energy flow rate for a well or group of wells. \\
MWE & fw-rate & ifno or boundname & -- & Simulated energy flow rate for a flowing well or group of flowing wells. \\
% observations specific to the mwe package
MWE & rate & mawno or boundname & -- & Simulated energy flow rate for a well or group of wells. \\
MWE & fw-rate & mawno or boundname & -- & Simulated energy flow rate for a flowing well or group of flowing wells. \\
MWE & rate-to-mvr & well or boundname & -- & Simulated energy flow rate that is sent to the MVE Package for a well or group of wells.\\
MWE & fw-rate-to-mvr & well or boundname & -- & Simulated energy flow rate that is sent to the MVE Package from a flowing well or group of flowing wells. \\
4 changes: 4 additions & 0 deletions doc/mf6io/gwe/gwe.tex
Original file line number Diff line number Diff line change
Expand Up @@ -133,6 +133,10 @@ \subsection{Streamflow Energy Transport (SFE) Package}
\subsection{Lake Energy Transport (LKE) Package}
\input{gwe/lke}

\newpage
\subsection{Multi-Aquifer Well Energy Transport (MWE) Package}
\input{gwe/mwe}

\newpage
\subsection{Flow Model Interface (FMI) Package}
\input{gwe/fmi}
Expand Down
55 changes: 55 additions & 0 deletions doc/mf6io/gwe/mwe.tex
Original file line number Diff line number Diff line change
@@ -0,0 +1,55 @@
Multi-Aquifer Well Energy Transport (MWE) Package information is read from the file that is specified by ``MWE6'' as the file type. There can be as many MWE Packages as necessary for a GWE model. Each MWE Package is designed to work with flows from a corresponding GWF MAW Package. By default \mf uses the MWE package name to determine which MAW Package corresponds to the MWE Package. Therefore, the package name of the MWE Package (as specified in the GWE name file) must match with the name of the corresponding MAW Package (as specified in the GWF name file). Alternatively, the name of the flow package can be specified using the FLOW\_PACKAGE\_NAME keyword in the options block. The GWE MWE Package cannot be used without a corresponding GWF MAW Package.

The MWE Package does not have a dimensions block; instead, dimensions for the MWE Package are set using the dimensions from the corresponding MAW Package. For example, the MAW Package requires specification of the number of wells (NMAWWELLS). MWE sets the number of wells equal to NMAWWELLS. Therefore, the PACKAGEDATA block below must have NMAWWELLS entries in it.

\vspace{5mm}
\subsubsection{Structure of Blocks}
\lstinputlisting[style=blockdefinition]{./mf6ivar/tex/gwe-mwe-options.dat}
\lstinputlisting[style=blockdefinition]{./mf6ivar/tex/gwe-mwe-packagedata.dat}
\lstinputlisting[style=blockdefinition]{./mf6ivar/tex/gwe-mwe-period.dat}

\vspace{5mm}
\subsubsection{Explanation of Variables}
\begin{description}
\input{./mf6ivar/tex/gwe-mwe-desc.tex}
\end{description}

\vspace{5mm}
\subsubsection{Example Input File}
\lstinputlisting[style=inputfile]{./mf6ivar/examples/gwe-mwe-example.dat}

\vspace{5mm}
\subsubsection{Available observation types}
Multi-Aquifer Well Energy Transport Package observations include well temperature and all of the terms that contribute to the continuity equation for each well. Additional MWE Package observations include energy flow rates for individual wells, or groups of wells; the well volume (\texttt{volume}); and the conductance for a well-aquifer connection conductance (\texttt{conductance}). The data required for each MWE Package observation type is defined in table~\ref{table:gwe-mweobstype}. Negative and positive values for \texttt{mwe} observations represent a loss from and gain to the GWE model, respectively. For all other flow terms, negative and positive values represent a loss from and gain from the MWE package, respectively.

\begin{longtable}{p{2cm} p{2.75cm} p{2cm} p{1.25cm} p{7cm}}
\caption{Available MWE Package observation types} \tabularnewline

\hline
\hline
\textbf{Stress Package} & \textbf{Observation type} & \textbf{ID} & \textbf{ID2} & \textbf{Description} \\
\hline
\endfirsthead

\captionsetup{textformat=simple}
\caption*{\textbf{Table \arabic{table}.}{\quad}Available MWE Package observation types.---Continued} \tabularnewline

\hline
\hline
\textbf{Stress Package} & \textbf{Observation type} & \textbf{ID} & \textbf{ID2} & \textbf{Description} \\
\hline
\endhead


\hline
\endfoot

\input{../Common/gwe-mweobs.tex}
\label{table:gwe-mweobstype}
\end{longtable}

\vspace{5mm}
\subsubsection{Example Observation Input File}
\lstinputlisting[style=inputfile]{./mf6ivar/examples/gwe-mwe-example-obs.dat}


1 change: 1 addition & 0 deletions doc/mf6io/gwe/namefile.tex
Original file line number Diff line number Diff line change
Expand Up @@ -36,6 +36,7 @@ \subsubsection{Explanation of Variables}
ESL6 & Energy Source Loading Package & * \\
SFE6 & Streamflow Energy Transport Package & * \\
LKE6 & Lake Energy Transport Package & * \\
MWE6 & Multi-Aquifer Well Energy Transport Package & * \\
OBS6 & Observations Option \\
\hline
\end{tabular*}
Expand Down
Loading

0 comments on commit 3a6b1ef

Please sign in to comment.