Skip to content
Snippets Groups Projects
Commit 674de755 authored by Lars von Wolff's avatar Lars von Wolff
Browse files

Fixed some details in a fem formulation

parent 4009c940
No related branches found
No related tags found
No related merge requests found
...@@ -120,7 +120,7 @@ public: ...@@ -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-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 + 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) + J1*(chns.u-param.phys.uAst)
//- pf_c*(u-param.phys.uAst)*(v*gradphi[i][0]) //- pf_c*(u-param.phys.uAst)*(v*gradphi[i][0])
+ (pf_c*(chns.u-param.phys.uAst)*chns.div_v + (pf_c*(chns.u-param.phys.uAst)*chns.div_v
...@@ -147,8 +147,8 @@ public: ...@@ -147,8 +147,8 @@ public:
RF gradxi3k = -param.tw.Sigma3*(chns.gradxi[k]/param.tw.Sigma1 + chns.gradxi2[k]/param.tw.Sigma2); RF gradxi3k = -param.tw.Sigma3*(chns.gradxi[k]/param.tw.Sigma1 + chns.gradxi2[k]/param.tw.Sigma2);
RF S = 0; RF S = 0;
S += - chns.xi2*(chns.gradpf[k]*pf_f - chns.pf * 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 *(chns.gradpf2[k]*pf_f - chns.pf2 * 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]); S += - 2*param.delta*pf3*(gradxi3k-chns.gradxi2[k]-chns.gradxi[k]);
//Stronger form - flow in interfaces possible //Stronger form - flow in interfaces possible
...@@ -162,7 +162,7 @@ public: ...@@ -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 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 //- 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 + 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 + param.vDis.eval(pf_f)*chns.v[k]*chns.vbasis.phi[i] //Velocity Dissipation in Solid
//- rho*g[k]*vphi[i] //Gravitation //- rho*g[k]*vphi[i] //Gravitation
- S //Surface Tension - S //Surface Tension
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Please register or to comment