Matt
Petroleum Engineering
Please login with a confirmed email address before reporting spam
Posted:
3 years ago
31 dic 2021, 08:37 GMT-5
Updated:
3 years ago
31 dic 2021, 08:46 GMT-5
** Update, the solution seems to be to define an ODE with the following weak form syntax:

where u is a variable from another PDE/physics interface, and umin is the dependent ODE variable. However, this method only works for Time Dependent solver where derivatives are continuous. This will not work for Time Discrete solver in which time derivatives are replaced by bdf() operator. I wonder if anyone has a solution for time discrete solver?
** Update, the solution seems to be to define an ODE with the following weak form syntax:
\marthrm{0 = -d\left(umin,TIME\right)*test\left(umin\right) + \left(u\leq umin\right)*\left(d\left(u,TIME\right)\leq 0 \right)*d\left(u,TIME\right)*test\left(umin\right)}
where u is a variable from another PDE/physics interface, and umin is the dependent ODE variable. However, this method only works for Time Dependent solver where derivatives are continuous. This will not work for Time Discrete solver in which time derivatives are replaced by bdf() operator. I wonder if anyone has a solution for time discrete solver?
Edgar J. Kaiser
Certified Consultant
Please login with a confirmed email address before reporting spam
Posted:
3 years ago
1 gen 2022, 16:23 GMT-5
Hi Matt,
have a look into this blog contribution by Walter Frei: https://www.comsol.de/blogs/using-the-previous-solution-operator-in-transient-modeling/
I have used this method in many transient and parametric studies. I think it should work in time-discrete as well, but haven't tried that so far.
Cheers
Edgar
-------------------
Edgar J. Kaiser
emPhys Physical Technology
www.emphys.com
Hi Matt,
have a look into this blog contribution by Walter Frei: https://www.comsol.de/blogs/using-the-previous-solution-operator-in-transient-modeling/
I have used this method in many transient and parametric studies. I think it should work in time-discrete as well, but haven't tried that so far.
Cheers
Edgar