scipy.integrate.odeint with pulsed time dependent args
I try to implement an chemical rate solver, where reaction rate constants change in time corresponding to a pulsed discharge.
The problem is that odeint increases the step size too fast.
So my question is if there is a facy way how to influence the step size change.
Images: Plot of density(t) and the Puls (discharge ON or OFF)
E.g. Solution of rate equation x: time, y: concentration of species
E.g. The discharge pulses x: time, y: electric field
So for some reactions the reaction rate (ki) is high within the puls and low in between.
My code is to long and confusing so I hope this will show what I mean.
def rate_func_odenit(Ci,time,ki,puls_funcion...):
if puls_func==1:
....ki(t)=function(t)...
else:
ki(t)=constant...
def ki(time,...):
odeint(rate_func_odenit, Ci0, time,args=(ki,...))
thanks for your ideas!
python scipy args odeint chemistry
asked Jan 18 at 17:36
Matthias KMatthias K
33
add a comment |
I try to implement an chemical rate solver, where reaction rate constants change in time corresponding to a pulsed discharge.
The problem is that odeint increases the step size too fast.
So my question is if there is a facy way how to influence the step size change.
Images: Plot of density(t) and the Puls (discharge ON or OFF)
E.g. Solution of rate equation x: time, y: concentration of species
E.g. The discharge pulses x: time, y: electric field
So for some reactions the reaction rate (ki) is high within the puls and low in between.
My code is to long and confusing so I hope this will show what I mean.
def rate_func_odenit(Ci,time,ki,puls_funcion...):
if puls_func==1:
....ki(t)=function(t)...
else:
ki(t)=constant...
def ki(time,...):
odeint(rate_func_odenit, Ci0, time,args=(ki,...))
thanks for your ideas!
python scipy args odeint chemistry
asked Jan 18 at 17:36
Matthias KMatthias K
33
"... my question is [...] how to influence the step size change..." Take a look at thehmaxargument ofodeint. It is described in the "Other Parameters" section of the odeint docstring.
– Warren Weckesser
Jan 18 at 18:28
add a comment |
I try to implement an chemical rate solver, where reaction rate constants change in time corresponding to a pulsed discharge.
The problem is that odeint increases the step size too fast.
So my question is if there is a facy way how to influence the step size change.
Images: Plot of density(t) and the Puls (discharge ON or OFF)
E.g. Solution of rate equation x: time, y: concentration of species
E.g. The discharge pulses x: time, y: electric field
So for some reactions the reaction rate (ki) is high within the puls and low in between.
My code is to long and confusing so I hope this will show what I mean.
def rate_func_odenit(Ci,time,ki,puls_funcion...):
if puls_func==1:
....ki(t)=function(t)...
else:
ki(t)=constant...
def ki(time,...):
odeint(rate_func_odenit, Ci0, time,args=(ki,...))
thanks for your ideas!
python scipy args odeint chemistry
asked Jan 18 at 17:36
Matthias KMatthias K
33
I try to implement an chemical rate solver, where reaction rate constants change in time corresponding to a pulsed discharge.
The problem is that odeint increases the step size too fast.
So my question is if there is a facy way how to influence the step size change.
Images: Plot of density(t) and the Puls (discharge ON or OFF)
E.g. Solution of rate equation x: time, y: concentration of species
E.g. The discharge pulses x: time, y: electric field
So for some reactions the reaction rate (ki) is high within the puls and low in between.
My code is to long and confusing so I hope this will show what I mean.
def rate_func_odenit(Ci,time,ki,puls_funcion...):
if puls_func==1:
....ki(t)=function(t)...
else:
ki(t)=constant...
def ki(time,...):
odeint(rate_func_odenit, Ci0, time,args=(ki,...))
thanks for your ideas!
python scipy args odeint chemistry
python scipy args odeint chemistry
asked Jan 18 at 17:36
Matthias KMatthias K
33
asked Jan 18 at 17:36
Matthias KMatthias K
33
asked Jan 18 at 17:36
Matthias KMatthias K
33
asked Jan 18 at 17:36
Matthias KMatthias K
33
asked Jan 18 at 17:36
Matthias KMatthias K
33
33
"... my question is [...] how to influence the step size change..." Take a look at thehmaxargument ofodeint. It is described in the "Other Parameters" section of the odeint docstring.
– Warren Weckesser
Jan 18 at 18:28
add a comment |
"... my question is [...] how to influence the step size change..." Take a look at thehmaxargument ofodeint. It is described in the "Other Parameters" section of the odeint docstring.
– Warren Weckesser
Jan 18 at 18:28
"... my question is [...] how to influence the step size change..." Take a look at the
hmax argument of odeint. It is described in the "Other Parameters" section of the odeint docstring.– Warren Weckesser
Jan 18 at 18:28
"... my question is [...] how to influence the step size change..." Take a look at the
hmax argument of odeint. It is described in the "Other Parameters" section of the odeint docstring.– Warren Weckesser
Jan 18 at 18:28
add a comment |
0
active
oldest
votes
Your Answer
StackExchange.ifUsing("editor", function () {
StackExchange.using("externalEditor", function () {
StackExchange.using("snippets", function () {
StackExchange.snippets.init();
});
});
}, "code-snippets");
StackExchange.ready(function() {
var channelOptions = {
tags: "".split(" "),
id: "1"
};
initTagRenderer("".split(" "), "".split(" "), channelOptions);
StackExchange.using("externalEditor", function() {
// Have to fire editor after snippets, if snippets enabled
if (StackExchange.settings.snippets.snippetsEnabled) {
StackExchange.using("snippets", function() {
createEditor();
});
}
else {
createEditor();
}
});
function createEditor() {
StackExchange.prepareEditor({
heartbeatType: 'answer',
autoActivateHeartbeat: false,
convertImagesToLinks: true,
noModals: true,
showLowRepImageUploadWarning: true,
reputationToPostImages: 10,
bindNavPrevention: true,
postfix: "",
imageUploader: {
brandingHtml: "Powered by u003ca class="icon-imgur-white" href="https://imgur.com/"u003eu003c/au003e",
contentPolicyHtml: "User contributions licensed under u003ca href="https://creativecommons.org/licenses/by-sa/3.0/"u003ecc by-sa 3.0 with attribution requiredu003c/au003e u003ca href="https://stackoverflow.com/legal/content-policy"u003e(content policy)u003c/au003e",
allowUrls: true
},
onDemand: true,
discardSelector: ".discard-answer"
,immediatelyShowMarkdownHelp:true
});
}
});
Sign up or log in
StackExchange.ready(function () {
StackExchange.helpers.onClickDraftSave('#login-link');
});
Sign up using Google
Sign up using Facebook
Sign up using Email and Password
Post as a guest
Required, but never shown
StackExchange.ready(
function () {
StackExchange.openid.initPostLogin('.new-post-login', 'https%3a%2f%2fstackoverflow.com%2fquestions%2f54258910%2fscipy-integrate-odeint-with-pulsed-time-dependent-args%23new-answer', 'question_page');
}
);
Post as a guest
Required, but never shown
0
active
oldest
votes
0
active
oldest
votes
active
oldest
votes
active
oldest
votes
Thanks for contributing an answer to Stack Overflow!
- Please be sure to answer the question. Provide details and share your research!
But avoid …
- Asking for help, clarification, or responding to other answers.
- Making statements based on opinion; back them up with references or personal experience.
To learn more, see our tips on writing great answers.
Sign up or log in
StackExchange.ready(function () {
StackExchange.helpers.onClickDraftSave('#login-link');
});
Sign up using Google
Sign up using Facebook
Sign up using Email and Password
Post as a guest
Required, but never shown
StackExchange.ready(
function () {
StackExchange.openid.initPostLogin('.new-post-login', 'https%3a%2f%2fstackoverflow.com%2fquestions%2f54258910%2fscipy-integrate-odeint-with-pulsed-time-dependent-args%23new-answer', 'question_page');
}
);
Post as a guest
Required, but never shown
Sign up or log in
StackExchange.ready(function () {
StackExchange.helpers.onClickDraftSave('#login-link');
});
Sign up using Google
Sign up using Facebook
Sign up using Email and Password
Post as a guest
Required, but never shown
Sign up or log in
StackExchange.ready(function () {
StackExchange.helpers.onClickDraftSave('#login-link');
});
Sign up using Google
Sign up using Facebook
Sign up using Email and Password
Post as a guest
Required, but never shown
Sign up or log in
StackExchange.ready(function () {
StackExchange.helpers.onClickDraftSave('#login-link');
});
Sign up using Google
Sign up using Facebook
Sign up using Email and Password
Sign up using Google
Sign up using Facebook
Sign up using Email and Password
Post as a guest
Required, but never shown
Required, but never shown
Required, but never shown
Required, but never shown
Required, but never shown
Required, but never shown
Required, but never shown
Required, but never shown
Required, but never shown
"... my question is [...] how to influence the step size change..." Take a look at the
hmaxargument ofodeint. It is described in the "Other Parameters" section of the odeint docstring.– Warren Weckesser
Jan 18 at 18:28