diff --git a/dune/phasefield/localoperator/fem_ffs_chns_r_monolithic.hh b/dune/phasefield/localoperator/fem_ffs_chns_r_monolithic.hh index 3adb2b014cf2dbcbc591b9537a5d08cbb0d9eefd..417b051a49aa2dc27e397ca2530f4b2041c6b6fc 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