Hello,
I solve phase-field equation, i.e. Cahn-Hilliard equation.
I have the following finite elements for the phase field and the chemical
potential
components of solution
FESystem<dim> fe_phase;
in the initialization list of constructor:
fe_phase (FE_Q<dim> (degree), 2)
(The first component - phase field, the second one - chemical potential)
I am getting the solution in
PetscWrappers::MPI::Vector solution_phase_n;
I am interested in modifying the vector of the solution - I
want to cut the undershoots and overshoots in phase-field component
of the solution. Namely I have values of phase field like -1.005 or 1.005,
however they should be only in <-1,1>.
Is it possible to modify only the relevant parts of the
PetscWrappers::MPI::Vector ?
I want to get afterwards the solution values and gradients using
fe_v_phase.get_function_values
fe_v_phase.get_function_gradients
to get the "repaired" values and gradients of phase field.
Or should I split the system, ie. firstly solve for the phase
field, then modify the phase field solution and afterwards
solve for chemical potential?
Thank You
Marek
--
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 [email protected].
For more options, visit https://groups.google.com/d/optout.
