time dependent solute flux

[Dipankar Dwivedi]

Hello,

How one can assign time dependent solute flux boundary conditions?

Under CONSTRAINT_LIST, how time works (as in below)? Can it be a time series?

 CONSTRAINT_LIST
    0.d0 U_source
    336000.d0 Initial

Thanks,
Dipankar

[Hammond, Glenn E]

Dipankar,

 

As you show below, you should provide a list of times and constraints.  It is essentially a timeseries; however there is no interpolation of geochemistry between times.  Time units can be included as follows:

 

TRANSPORT_CONDITION

  TIME_UNITS h

  CONSTRAINT_LIST

    0. name_of_constraint

    1. name_of_constraint

    2. name_of_constraint

  /

END

 

TRANSPORT_CONDITION

  CONSTRAINT_LIST

    0. h name_of_constraint

    1. h name_of_constraint

    2. h name_of_constraint

  /

END

 

Hope this helps.

 

Glenn

[Peter Lichtner]

Glenn: is the TIME keyword still used:

      case('TIME')

        call InputReadDouble(input,option,default_time)

        call InputErrorMsg(input,option,'TIME','CONDITION') 

      case('TIME_UNITS') 

        call InputReadWord(input,option,word,PETSC_TRUE) 

        call InputErrorMsg(input,option,'UNITS','CONDITION')   

        call StringToLower(word)

        select case(trim(word))     

          case('s','sec','min','m','hr','h','d','day','y','yr')

            default_time_units = trim(word)         

        end select          

      case('CONSTRAINT_LIST')

…Peter

[Dipankar Dwivedi]

Thanks, Glenn. I will try this.

Best,

Dipankar 

[Hammond, Glenn E]

It is still supported, but no one uses it.  Nor does it have any practical use unless the time is 0.  A constraint has to be defined at all times.

 

Glenn