Commit dbfed74e authored by Matthieu Constant's avatar Matthieu Constant
Browse files

test avec equations de Stokes

parent 8ddf62b2
......@@ -241,7 +241,7 @@ static void f_u_0(
for (int k = 0; k < D; ++k) {
utau += u[k]*(du[j][k]-u[j]*dq[k]/q);
}
f[j] = rho*((u[j]-uold[j])/dt + u[j]/q*divu + utau/q);
f[j] = 0*rho*((u[j]-uold[j])/dt + 0*u[j]/q*divu + 0*utau/q);
}
}
......@@ -265,7 +265,7 @@ static void f_u_1(double f[D][D], double u[D], const double du[D][D], const doub
double tau[D][D];
for (int i = 0; i < D; ++i) {
for (int j = 0; j < D; ++j) {
tau[i][j] = du[i][j]-u[i]*dq[j]/q;
tau[i][j] = du[i][j]-0*u[i]*dq[j]/q;
}
}
for (int j = 0; j < D; ++j) {
......@@ -289,7 +289,7 @@ static void fluid_problem_assemble_system(FluidProblem *problem, double *rhs, co
HXTLinearSystem *lsys = problem->linear_system;
const Mesh *mesh = problem->mesh;
//a generaliser
double deps = 1e-8;
double deps = 1e-10;
const double *solution = problem->solution;
const double *porosity = problem->porosity;
const double *old_porosity = problem->old_porosity;
......
......@@ -49,7 +49,7 @@ tEnd = 1000 # final time
#numerical parameters
lcmin = 0.1 # mesh size
dt = 5e-2 # time step
alpha = 1e-5 # stabilization coefficient
alpha = 7e-6 # stabilization coefficient
epsilon = alpha*lcmin**2 /nu # stabilization parametre
print('epsilon',epsilon)
......
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