Discussion Closed This discussion was created more than 6 months ago and has been closed. To start a new discussion with a link back to this one, click here.

solve stiff ODE with femlab 3.5

Please login with a confirmed email address before reporting spam

I have made the set of stiff ordinary differential equations system of densities of species in form :dn/dt = source - n*perte, for this my question, for do this,i use the interface "Global equations" , but with step adaptatif wich change at each novel solution.

how can specify this adaptatif step in box of time-dependant of software Femlab 3.5 ???

hoping for a response full with your scientific enthusiasm.

1 Reply Last Post 2010年3月14日 GMT-4 17:22
Ivar KJELBERG COMSOL Multiphysics(r) fan, retired, former "Senior Expert" at CSEM SA (CH)

Please login with a confirmed email address before reporting spam

Posted: 1 decade ago 2010年3月14日 GMT-4 17:22
Hello Isabelle

Sorry but I have some problems to follow your explanations:

1) are you in the PDE application mode coefficients "c" or global "g" mode or in an other module ?
2) you add a "Gobal Equation" of the type: dn/dt = source - n*loss, where n is probably the species density, but is this one or several equations (dn1/dt = source1 - n*loss1 & dn2/dt = source2 - n*loss2 ... ) ?
3) as dn/dt is a time derivative, I understand that you are in the transient solver mode, so you must have defined a time of the type t = range(start, step, end), or perhaps a t = log10(range(start_power, step_power, end_power)) if you need a logaritmic stepping
4) I do not understand your adaptive step issue, do you mean that "source" and "loss" are dependent on the time too, and you have a non-linear problem ?, or is it the time stepping you would want to have adaptive ?

I understand that it is not easy to write all this in English, but it is a good exercice, for all of us ;),
in the worst case write a couple of sentences in French too, there will always be someone to translate. But the idea of keeping it in English is to give an opportunity for the majority to follow

Good luck
Ivar
Hello Isabelle Sorry but I have some problems to follow your explanations: 1) are you in the PDE application mode coefficients "c" or global "g" mode or in an other module ? 2) you add a "Gobal Equation" of the type: dn/dt = source - n*loss, where n is probably the species density, but is this one or several equations (dn1/dt = source1 - n*loss1 & dn2/dt = source2 - n*loss2 ... ) ? 3) as dn/dt is a time derivative, I understand that you are in the transient solver mode, so you must have defined a time of the type t = range(start, step, end), or perhaps a t = log10(range(start_power, step_power, end_power)) if you need a logaritmic stepping 4) I do not understand your adaptive step issue, do you mean that "source" and "loss" are dependent on the time too, and you have a non-linear problem ?, or is it the time stepping you would want to have adaptive ? I understand that it is not easy to write all this in English, but it is a good exercice, for all of us ;), in the worst case write a couple of sentences in French too, there will always be someone to translate. But the idea of keeping it in English is to give an opportunity for the majority to follow Good luck Ivar

Note that while COMSOL employees may participate in the discussion forum, COMSOL® software users who are on-subscription should submit their questions via the Support Center for a more comprehensive response from the Technical Support team.