How to use integrated values in equation.cc file

[Andrew Li]

Hi everyone, 

I am trying to use the integrated value of a field in equation.cc. For example, I would like to use f_tot, which is usually used as a post-processing variable. Another use case would be to use the average value/total sum of the concentration field times a position dependent field. Can anyone give me some recommendations? Thanks!

Best, 

Andrew

[David Montiel Taboada]

Hello, Andrew

Something like what you describe is done in the recent application allenCahn_conserved. I encourage you to look into the documentation. 

Let me know if you have any questions.

Best,

David

--
You received this message because you are subscribed to the Google Groups "PRISMS-PF Users" group.
To unsubscribe from this group and stop receiving emails from it, send an email to prisms-pf-use...@googlegroups.com.
To view this discussion on the web visit https://groups.google.com/d/msgid/prisms-pf-users/034e9b61-e078-45b1-836c-19a9ab57f6c1n%40googlegroups.com.

[Andrew Li]

Hi David, 

I've looked into the allenCahn_conserved application and tried running it, but am getting an error in the customPDE.h file (no member "invM") in the custom solveIncrement function. Is this something that is a real error or am I doing something wrong? I've included the slurm output for your reference.  

Best, 

Andrew

[David Montiel]

Hello, Andrew 

I've replied to this question privately, but please reply to my email within this thread so other people can use the information for reference.

Best,

David

To view this discussion on the web visit https://groups.google.com/d/msgid/prisms-pf-users/6bc3ab74-be5f-496b-8d30-3a8d202d10c8n%40googlegroups.com.

[Andrew Li]

Hi David, 

Thanks! I ran the application twice (1 and 8 processors). Both instances give the same error ("No member "invM"). I am using dea...@9.4.2. 

Best,

Andrew

[Andrew Li]

And this the original application with no modifications

[Jason Landini]

Hi Andrew,

This is a bug on our end. We updated the naming convention on invM.

I'll have a patch out by the end of the day. After that you can update your fork / local repo.

Best,

Jason

To view this discussion on the web visit https://groups.google.com/d/msgid/prisms-pf-users/6bc3ab74-be5f-496b-8d30-3a8d202d10c8n%40googlegroups.com.

[David Montiel]

Thank you for looking into this, Jason.

David

To view this discussion on the web visit https://groups.google.com/d/msgid/prisms-pf-users/CAHaHfUAL6LUPz7AhDUPdhidZsPGFEt3ue40qxJmUhnM0HX0HOQ%40mail.gmail.com.

[Andrew Li]

Hi Jason and David, 

Thanks for the update! On a related note, is this custom solveIncrememnt function something I can port directly to other applications such as CH, CHAC, corrosion and etc? Or are additional modifications required?

Best

Andrew

[David Montiel]

Hello, Andrew

You could port it but modifying it accordingly to obtain the integrated the field(s) you are interested in each case.

David

[Andrew Li]

Thanks! I assume that means I just need to create a new integrated field for each global quantity I’m trying to conserve. That’s very convenient. Really appreciate the help!

[Jason Landini]

Alright the bug should be fixed now. Let me know if you have any other trouble.