Versions Compared

Key

  • This line was added.
  • This line was removed.
  • Formatting was changed.

...

Panel

Table of content:

Table of Contents

 

 

Title:

...

 Variably Saturated Flow Model (VSFM): A Unified treatment of hydrologic processes in the unsaturated-saturated zone

Requirements and Design

ACME

...

Land Group

Date:  

Summary

The purpose of this section is to summarize what capability is to be added to the ACME Ocean and Ice system Land Model through this design process.
It should be clear what new code will do that the current code does not. Summarizing the primary challenges with respect to software design and implementation is also appropriate for this section. Finally, this statement should contain general statement with regard to what is “success.”

 

Requirements

Requirement: name-of-requirement-here

Date last modified:  
Contributors:  (add your name to this list if it does not appear)


Each requirement is to be listed under a ”section” heading, as there will be a one-to-one correspondence between requirements, design, proposed imple- mentation and testing. Requirements should not discuss technical software issues, but rather focus on model capability. To the extent possible, require- ments should be relatively independent of each other, thus allowing a clean design solution, implementation and testing plan.

 

Algorithmic Formulations

Design solution: short-description-of-proposed-solution-here

Date last modified: 10/20/2015
Contributors: Gautam Bisht

 

For each requirement, there is a design solution that is intended to meet that requirement. Design solutions can include detailed technical discussions of PDEs, algorithms, solvers and similar, as well as technical discussion of performance issues. In general, this section should steer away from a detailed discussion of low-level software issues such as variable declarations, interfaces and sequencing.

 

  • The governing PDE for flow in porous media is given as:

d/dt (rho * phi * sat) = - Div (rho * q) + Q  ... (1)

where 

rho = density of water, which is a function of pressure and temperature,

phi = porosity,

sat = liquid saturation

q = darcy flux

Q = source/sink of water (e.g. infiltration, evapotranspiration, bare soil evaporation)

Div = divergence operator

 

  • The Darcy flux is given by q = - k kr / mu * Del (P)   ... (2)

where 

k = permeability,

kr = relative permeability.

mu = viscosity,

P = liquid pressure

 

  • Following discretization schemes are used to solve Eq(2):
    • Spatial discretization: Finite volume,
    • Temporal discretization: Backward Euler scheme 

Design and Implementation

Implementation: short-desciption-of-implementation-here

Date last modified: 10/20/2015
Contributors: Gautam Bisht

 

Numerical implementation of the VSFM would use the following capabilities of PETScSNES: To
This section should detail the plan for implementing the design solution for requirement XXX. In general, this section is software-centric with a focus on software implementation. Pseudo code is appropriate in this section. Links to actual source code are appropriate. Project management items, such as svn branches, timelines and staffing are also appropriate. How do we typeset pseudo code?
  • Initial prototype of the VSFM model will be developed in MATLAB.
  • The VSFM implementation in ACME would use PETSc.
    • SNES solver of PETSc will be used to solve the system of nonlinear equations that is obtained after applying spatial and temporal discretization.
    • DMComposite: The future goal of ACME is to solve tightly coupled transport of water in the soil with other relevant physics (e.g. transport of water in roots [ACME V2]). Thus, the current VSFM will use DMComposite feature of PETSc. DMComposite allows application to 
This section should detail the plan for implementing the design solution for requirement XXX. In general, this section is software-centric with a focus on software implementation. Pseudo code is appropriate in this section. Links to actual source code are appropriate. Project management items, such as svn branches, timelines and staffing are also appropriate. How do we typeset pseudo code?
    • an application to brea

 

Planned Verification and Unit Testing 

Verification and Unit Testing: short-desciption-of-testing-here

Date last modified:

  

 10/20/2015 

Contributors: 

(add your name to this list if it does not appear)

Gautam Bisht

 

How will XXX be tested? i.e. how will be we know when we have met requirement XXX. Will these unit tests be included in the ongoing going forward?

Planned Validation Testing 

Validation Testing: short-desciption-of-testing-here

Date last modified: 10/20/2015 

Contributors: 

(add your name to this list if it does not appear)

Gautam Bisht

 

How will XXX be tested? What observational or other dataset will be used?  i.e. how will be we know when we have met requirement XXX. Will these unit tests be included in the ongoing going forward?

  • The VSFM will be tested against previously benchmark problems published.
  • Additionally, VSFM will be tested against PFLOTRAN for certain benchmark problems.

Planned Performance Testing 

Performance Testing: short-desciption-of-testing-here

Date last modified:
Contributors: (add your name to this list if it does not appear)

 

How will XXX be tested? i.e. how will be we know when we have met requirement XXX. Will these unit tests be included in the ongoing going forward?

 

...