Commit 82501b31 authored by Matthieu Constant's avatar Matthieu Constant
Browse files

debut q

parent 8556a2ca
......@@ -306,7 +306,7 @@ static void fluid_problem_assemble_system(FluidProblem *problem, double *rhs, co
local_vector[iphi+N_SF*j] += jw*(
dphisigma
//+rho*phii*dudt[j]-rho*utau
+rho*phii*(dudt[j]+u[j]/c*divu+utau/c)
+rho*phii*(dudt[j]+u[j]/q*divu+utau/q)
-p*dphii[j]
);
}
......@@ -335,16 +335,19 @@ static void fluid_problem_assemble_system(FluidProblem *problem, double *rhs, co
int U = j;
LOCAL_MATRIX(U,P) += -jw*dphii[j]*phij;
LOCAL_MATRIX(P,U) += jw*phii*dphij[j];
double utau =0;
for (int k = 0; k < DIMENSION; ++k) {
int V = k;
utau += u[k]*tau[j][k];
//LOCAL_MATRIX(U,V) += -jw*rho*u[j]*dphii[k]*phij/c+jw*2*mu*dphii[k]*dtau[j]*0.5;
LOCAL_MATRIX(U,V) +=
jw*rho*phii*phij*tau[j][k]/c
jw*rho*phii*phij*tau[j][k]/q
+jw*2*mu*dphii[k]*dtau[j]*0.5
+jw*rho*phii*u[j]*dphij[k]/c;
+jw*rho*phii*u[j]*dphij[k]/q;
}
//LOCAL_MATRIX(U,U) += jw*mu*2*0.5*dphiidtau + jw*rho*(phii*phij/dt-phij*dphii[j]*u[j]/c);
LOCAL_MATRIX(U,U) += jw*mu*2*0.5*dphiidtau + jw*rho*phii*(phij/dt+phij/c*divu+udtau/c);
LOCAL_MATRIX(U,U) += jw*mu*2*0.5*dphiidtau + jw*rho*phii*(phij/dt+phij/q*divu+udtau/q);
LOCAL_MATRIX(U,Q) += -jw*rho*phii*(u[j]*divu*phij/(q*q)+utau*phij/(q*q));
}
LOCAL_MATRIX(P,P) += jw*epsilon*dphiidphij;
LOCAL_MATRIX(Q,Q) += -jw*phii*phij;
......
Supports Markdown
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment