abarret / multiphase-stokes

Solver a mixture of fluids based on IBAMR
1 stars 0 forks source link

Variable coefficients #71

Closed abarret closed 5 months ago

abarret commented 6 months ago

I want to be able to specify a spatially (and temporally) varying drag coefficient. This can be set via a CartGridFunction for xi in the hierarchy integrator, but this information needs to be propagated further along the line to the operators and preconditioner. We could make VCTwoPhaseStaggeredStokesOperator take flag to specify either variable or constant drag coefficients and then have an if statement in the apply function. This seems not ideal, as what happens if we want to, e.g. specify variable viscosity later as well? An alternative is to have a separate operator for each type of variable coefficient.

We first need to break down the operators and preconditioner to smaller chunks that are easier to replace.

abarret commented 6 months ago

Also, we should replace the names TwoFluid with Multiphase and remove VC and INS.