Maximu value for dt

General discussions about how to formulate a script for FlexPDE.

Maximu value for dt

Postby iyaniwura on Wed May 16, 2018 3:38 pm

Dear Moderator,

I'm solving a coupled PDE-ODE model using FlexPDE 6 and observe something that looks like numerical instability around time = 100. I observed that around this time, the dt is large, it is possible for me to specify a maximum value for dt. Maybe I can prevent this instability by keeping the dt as low as possible.

I have a copy of my code incase you will like to see it.

Thanks for your attention.

Last edited by iyaniwura on Wed May 23, 2018 11:59 am, edited 1 time in total.
Posts: 2
Joined: Wed May 09, 2018 1:29 pm

Re: Maximu value for dt

Postby moderator on Thu May 17, 2018 2:08 pm

I ran this script on version 6.51 with the ERRLIM control removed (to let it run faster).
The oscillatory behavior is established from the start, and merely gets more extreme after time 90.
The timestep I observed is about 0.05, and the oscillation period is about 10, so there are 200 timesteps per oscillation.
This is not a timestep size issue.
The oscillation is inherent in your equations.
Posts: 722
Joined: Tue Jan 11, 2011 1:45 pm

Re: Maximu value for dt

Postby iyaniwura on Thu May 17, 2018 5:47 pm

Thanks for the prompt response.

I was also thinking that the value of dt should not determine the stability of the solution since FlexPDE 6 uses BDF. My next guess will be the position of the cells, actually, I just discovered that I made a mistake when computing the positions of the cells.

I will check the code again and let you know if I am able to resolve the issue.

Thanks once again.
Posts: 2
Joined: Wed May 09, 2018 1:29 pm

Return to FlexPDE Discussion

Who is online

Users browsing this forum: No registered users and 1 guest