pde - Freefem Fisher's equation -
I'm new to free +++, the problem I am trying to solve is Fisher's equation:
du / dt = d ^ 2u / dx ^ 2 + d ^ 2u / dy ^ 2 + k * u * (1-U) du / dn = 0 - condition of limit
I have a weak problem The problem is with Futur formula (Uh, VH) = int2d (th) (uh * VH / dt + grad (uh)). * Grad (VH) - Int 2D (Th) (K * UH * VH) + Int 2D (Th) (UH 0 * VH / D T) - Int 2D (Gu) (K * Uh * VH * Uh);
Can you please tell me what do I do wrong? There is something wrong with the final conditions.
This is a 2D transient diffusion / conduction equation with a temperature-dependent, non-linear generation duration .
If you skip the non-linear generation period, the equation 2D should look like a weak form for transient diffusion / conduction equation.
How to call free + + non-linear words? How do you plan to handle it?
You know, after all, it means that the solution is a very different animal. You have to walk in time to solve this (eg a Newton-Repson Solver).
The algorithm becomes an iterative, non-linear one. You will not solve for now; You will solve a salary increase and repeat until convergence.
You can repeat the last word like this:
d (k * u (1-u)) = k * du (1-u) - k * U * du = k * (1-2 * u) * du ~ k * du
You still have a product that is * non-linear. What to do? throw it away.
Now you are solving a non-linear transient equation for
du .
Comments
Post a Comment