Evaluator for an individual transport component

[Sam Shaheen]

Hi all,

In a multi-component transport simulation, I am wondering if it's possible to create an evaluator that considers only a single component. For example, the first-order_tracer.xml problem from the most recent short course includes first-order decay of a tracer through the source term tracer2_rate:

<ParameterList name="tracer2_rate" type="ParameterList"> 

  <Parameter name="evaluator type" type="string" value="multiplicative evaluator"/>

  <Parameter name="coefficient" type="double" value="-1.0e-3"/> 

  <Parameter name="dependencies" type="Array(string)" value="{mole_fraction}"/> 

</ParameterList>

In this case, only a single tracer was included in the problem. If there were two tracers in the problem (e.g., Tracer1 and Tracer2) is it possible to call only one of the transport components (e.g., Tracer2) into an evaluator? I tried testing dependencies mole_fraction.Tracer2 and mole_fraction.1 in the evaluator above, but was getting error messages that "Evaluator "mole_fraction.Tracer2" @ "transport_coupler_next" cannot be created in State". Just curious if it's possible to target a single component as such, and if so what I should call the dependency.

Thanks,

Sam

[Sergi Molins Rafa]

Hi Sam

I created an Issue after the course to remind us we may want more flexibility in evaluators for multicomponent fields.

https://github.com/amanzi/ats/issues/326

I do realize the issue description is not too clear but in any case I think we have not done anything about it unless I missed something.

Sergi

--
You received this message because you are subscribed to the Google Groups "Amanzi-ATS Users" group.
To unsubscribe from this group and stop receiving emails from it, send an email to ats-users+...@googlegroups.com.
To view this discussion visit https://groups.google.com/d/msgid/ats-users/79652970-b7a3-480c-b2da-1cee2ec37693n%40googlegroups.com.

[Sergi Molins Rafa]

Hi Sam

I updated https://github.com/amanzi/ats/issues/326 with information about a branch that I started to attempt to make this work for the short course example. Do you think something like this would work for your case? Would you be available to test it?

Thanks

Sergi

[Sam Shaheen]

Hi Sergi,

I built this new version and tested out both a more sophisticated reactive transport problem then a two-tracer version of the short course demo. In both cases, I am getting an error message 

"terminate called after throwing an instance of 'Errors::Message' 

what(): Requested components are not consistent with previous request."

I attached my input script for the tracer demo in case I overlooked something in making it a multi-component problem.

Thanks,
Sam