Timeline for Solving system of PDEs with NDSolve
Current License: CC BY-SA 4.0
13 events
| when toggle format | what | by | license | comment | |
|---|---|---|---|---|---|
| Sep 5, 2019 at 4:42 | comment | added | user21 |
@gumpel have a look at the ref page for NDSolveValue there are plenty examples. Or replace NDSolveValue with NDSolve and use what you have now; I find the NDSolveValue-form the result is returned in more convenient.
|
|
| Sep 5, 2019 at 4:38 | vote | accept | gumpel | ||
| Sep 5, 2019 at 4:37 | comment | added | gumpel | OK, I chose to upgrade to version 12 :) so your code works now. However, could you also help me with plotting this? In 11.3, it was just Plot3D[Evaluate[{z[200, x, y]} /. soln2d], {x, 5, 10}, {y, 5, 10}, but this does not work here. | |
| Sep 4, 2019 at 5:42 | comment | added | user21 |
@gumpel, may a search for NDSolve::mxsst can help, I think there ample examples for this message on this site.
|
|
| Sep 4, 2019 at 5:09 | comment | added | gumpel | OK, thank you. Do you have a suggestion for an alternative? | |
| Sep 4, 2019 at 4:54 | comment | added | user21 | @gumpel, nonlinear finite elements got added in version 12.0 | |
| Sep 4, 2019 at 0:18 | comment | added | gumpel |
Sorry for the late reply. I have version 11.3. It shows the following errors: CoefficientArrays::poly: pp$3209-0.05 (pp$3210+pp$3211)-pp (1-0.1732 pp-z/(1+pp^2)) is not a polynomial. and NDSolveValue::femper: PDE parsing error of {pp$3209-0.05 (pp$3210+pp$3211)-pp (1-0.1732 pp-z/(1+pp^2)),-(-0.3+(1.4434 pp)/(1+pp^2)) z+z$3209-0.05 (z$3212+z$3213)}. Inconsistent equation dimensions.
|
|
| Sep 3, 2019 at 9:10 | comment | added | user21 | @gumpel, the code runs without errors. If it does not for you: what version do you have? what issues does it show? | |
| Sep 3, 2019 at 8:42 | comment | added | gumpel | Or does the code you have given work for you? | |
| Sep 3, 2019 at 8:01 | comment | added | gumpel | OK, thanks. Any idea regarding the second problem? | |
| Sep 3, 2019 at 7:59 | comment | added | user21 | @gumpel, you can store an interpolating function just like any other expression. | |
| Sep 3, 2019 at 7:55 | comment | added | gumpel | 1: OK, then I'll try this. Is there also the possibility to store this in a file and load it later as initial condition? I ask as Mathematica sometimes crashes. 2: Does the "[Element] Rectangle[..]" contain zero flux bc? I copied this code and I get the following error messages: CoefficientArrays::poly: pp$2407-0.05 (pp$2408+pp$2409)-pp (1-0.1732 pp-z/(1+pp^2)) is not a polynomial. NDSolveValue::femper: PDE parsing error of {pp$3017-0.05 (pp$3018+pp$3019)-pp (1-0.1732 pp-z/(1+pp^2)),-(-0.3+(1.4434 pp)/(1+pp^2)) z+z$3017-0.05 (z$3020+z$3021)}. Inconsistent equation dimensions. | |
| Sep 3, 2019 at 7:22 | history | answered | user21 | CC BY-SA 4.0 |