Hello, It depends of the value of the Reynolds number and the gradients of K with respect to x or t, but generally the last two terms do not generally pose problems. The first additional term to the right leads to a mass matrix, which is well conditioned. The second term itself is trickier. If you want to use an analytical Jacobian formulation and Newton's method, you will need to calculate the Frechet derivative of the velocity magnitude. You can also use a Picard iteration for this term, which will greatly simplify the expression of the Jacobian at the cost more Newton iteration. Issues generally arise when you have jumps the value of K in space. Then generally it is better to use some sort of upwinding to prevent oscillations in the velocity field (i.e SUPG).
On Friday, 27 September 2019 22:41:22 UTC-4, FU wrote: > > Hi， > I want to solve the problem about Darcy-Brinkman-Forchheimer equations, > but don't know how to discretizate this equation. > > [image: Darcy equation.png] > > This equation has a similar N-S equation. But the discretization of the > last item of the equation and the programming statements are somewhat > unclear. > > > Yours， > > FU > -- The deal.II project is located at http://www.dealii.org/ For mailing list/forum options, see https://groups.google.com/d/forum/dealii?hl=en --- You received this message because you are subscribed to the Google Groups "deal.II User Group" group. To unsubscribe from this group and stop receiving emails from it, send an email to dealii+unsubscr...@googlegroups.com. To view this discussion on the web visit https://groups.google.com/d/msgid/dealii/c5bd3188-e6b8-4338-92f3-3a0ec422bb03%40googlegroups.com.
