From 674de755788ba2c7f1556c0d4d8fb65125167408 Mon Sep 17 00:00:00 2001
From: Lars von Wolff <lars.von-wolff@ians.uni-stuttgart.de>
Date: Mon, 14 Feb 2022 15:54:31 +0100
Subject: [PATCH] Fixed some details in a fem formulation

---
 .../phasefield/localoperator/fem_ffs_chns_r_monolithic.hh | 8 ++++----
 1 file changed, 4 insertions(+), 4 deletions(-)

diff --git a/dune/phasefield/localoperator/fem_ffs_chns_r_monolithic.hh b/dune/phasefield/localoperator/fem_ffs_chns_r_monolithic.hh
index 3adb2b0..417b051 100644
--- a/dune/phasefield/localoperator/fem_ffs_chns_r_monolithic.hh
+++ b/dune/phasefield/localoperator/fem_ffs_chns_r_monolithic.hh
@@ -120,7 +120,7 @@ public:
 
           (pf_c_reg-pf_c_reg_Old)*(chns.u-param.phys.uAst)*chns.basis.phi[i]/param.time.dt
           + pf_c_reg_Old*(chns.u-chnsOld.u)*chns.basis.phi[i]/param.time.dt
-          + param.phys.D*pf_c_reg_Old*(chns.gradu*chns.basis.gradphi[i][0])
+          + param.phys.D*pf_c_reg*(chns.gradu*chns.basis.gradphi[i][0])
           + J1*(chns.u-param.phys.uAst)
           //- pf_c*(u-param.phys.uAst)*(v*gradphi[i][0])
           + (pf_c*(chns.u-param.phys.uAst)*chns.div_v
@@ -147,8 +147,8 @@ public:
 
           RF gradxi3k = -param.tw.Sigma3*(chns.gradxi[k]/param.tw.Sigma1 + chns.gradxi2[k]/param.tw.Sigma2);
           RF S = 0;
-          S += - chns.xi2*(chns.gradpf[k]*pf_f  - chns.pf  * gradpf_f[k])/(pf_f*pf_f);
-          S += - chns.xi *(chns.gradpf2[k]*pf_f - chns.pf2 * gradpf_f[k])/(pf_f*pf_f);
+          S += - chns.xi2*pf_f*(chns.gradpf[k]* pf_f  - chns.pf  * gradpf_f[k])/(pf_f*pf_f);
+          S += - chns.xi *pf_f*(chns.gradpf2[k]*pf_f - chns.pf2 * gradpf_f[k])/(pf_f*pf_f);
           S += - 2*param.delta*pf3*(gradxi3k-chns.gradxi2[k]-chns.gradxi[k]);
 
           //Stronger form - flow in interfaces possible
@@ -162,7 +162,7 @@ public:
           r.accumulate(child(chns.vspace,k),i, (   rho_f_reg*(chns.v[k]-chnsOld.v[k])*chns.vbasis.phi[i]/param.time.dt          //TimeEvolution
                                             //- ns.p * (pf_f*ns.vbasis.gradphi[i][0][k] + ns.vbasis.phi[i]*gradpf_f[k])            //Pressure in weak form
                                             + pf_f * (chns.gradp[k] * chns.vbasis.phi[i]) //Pressure in strong form
-                                            + param.phys.eval_mu(pf_f) * (chns.jacv[k]*chns.vbasis.gradphi[i][0]) //+vphi[i]*(jacv[k]*gradpf)   )  //Viscosity
+                                            + param.phys.eval_mu(pf_f) * (chns.jacv[k]*chns.vbasis.gradphi[i][0]) //+symmetry   )  //Viscosity
                                             + param.vDis.eval(pf_f)*chns.v[k]*chns.vbasis.phi[i]              //Velocity Dissipation in Solid
                                             //- rho*g[k]*vphi[i]                //Gravitation
                                             - S                                //Surface Tension
-- 
GitLab