Chorin's method in NSBENCH
In the nsbench framework (test environment for different Navier-Stokes solvers in Fenics) at the implementation of Chorin's projection method for the tentative velocity step the following code is available:
# Tentative velocity step
F1 = (1 / k) * inner(u - u0, v) * dx + inner(grad(u0) * u0, v) * dx + \
nu * inner(grad(u), grad(v)) * dx - inner(f, v) * dx
a1 = lhs(F1)
L1 = rhs(F1)
I am trying to understand the method (and the weak form derivations) and was wondering why the diffusion term in the original formulation nu*laplace(u^n) gets transformed to nu * inner(grad(u), grad(v)) * dx in the dolfin code? Shouldn't this be a weak form for nu*laplace(u^*)? (based on the notation from Wikipedia eqn 1. (http://
Question information
- Language:
- English Edit question
- Status:
- Solved
- For:
- DOLFIN Edit question
- Assignee:
- No assignee Edit question
- Solved by:
- Peter Maday
- Solved:
- Last query:
- Last reply: