On 3/13/23 07:09, Tom Li wrote:
I used the matrixfree to solve a solidificatio phase field problem. I think I
have no way to change the value in the local matrix and right hand side. And
I'm trying to use a variable substitution. ie .the old numerical
solutiuon(phi) is in the range of [-1,1]. I used psi(x) = sqrt(2.0) *
atanh(phi(x)) to replace the old value phi. After this change the solution
will become linear in the hole domain , psi will be very large when the point
is far away from liquid/solid interface . So I want to limit psi in a range
.for example[-10,10]. What should I do in dealii.
Tom:
you don't say how you want that to happen. There are many ways in which you
can limit psi to that range:
* You can change the equations in such a way that they only produce psi values
in that range
* You can just cut off values that are larger
* You can transform values, for example you can transform
psi -> atan(psi)*10/(pi/2)
In other works, your question describes what you want to be (a state) rather
than how you want to get there (an action). Software describes actions, not
states, and you first need to think of what you want that action to be.
Best
W.
--
------------------------------------------------------------------------
Wolfgang Bangerth email: bange...@colostate.edu
www: http://www.math.colostate.edu/~bangerth/
--
The deal.II project is located at http://www.dealii.org/
For mailing list/forum options, see
https://groups.google.com/d/forum/dealii?hl=en
---
You received this message because you are subscribed to the Google Groups "deal.II User Group" group.
To unsubscribe from this group and stop receiving emails from it, send an email
to dealii+unsubscr...@googlegroups.com.
To view this discussion on the web visit
https://groups.google.com/d/msgid/dealii/d92dcfd4-6893-e59c-279e-b2bf4f91e149%40colostate.edu.
