Get Help

Turn on suggestions

Auto-suggest helps you quickly narrow down your search results by suggesting possible matches as you type.

Showing results for

- Community
- :
- PTC Mathcad
- :
- PTC Mathcad
- :
- rate constants determination: from ODEs solving v....

Options

- Subscribe to RSS Feed
- Mark Topic as New
- Mark Topic as Read
- Float this Topic for Current User
- Bookmark
- Subscribe
- Mute
- Printer Friendly Page

Auto-suggest helps you quickly narrow down your search results by suggesting possible matches as you type.

Showing results for

Highlighted

- Mark as New
- Bookmark
- Subscribe
- Mute
- Subscribe to RSS Feed
- Permalink
- Email to a Friend
- Notify Moderator

07-03-2013
09:01 AM

07-03-2013
09:01 AM

Again in the book "chemical kinetics with maple and mathcad", chapter 4 introduces different ways for rate constants determination.

Take figure 4.8 as an example (corresponding attached file is "Fig_4_08_PineneIsomerization_InvTask_MC14" ), numerical integration method is used searching rate constants:

but i think that solving ODEs is more direct and general than integration method, so I solve the problem with ODEsolve method:

Both methods looks good enough, but differential method results are closer to that of DYANFIT than for integral method.

But differential method seems much more initial-condition sensitive than integral method, e.g. if guess values of k1--k4 are all set to 1, the calculation will be very time-consuming. and I dont know how to improve this situation.

Does anybody give any suggestions on "lowering/ reducing" the initial - condition sensitivity of differential method?

Solved! Go to Solution.

Labels:

1 ACCEPTED SOLUTION

Accepted Solutions

Highlighted

- Mark as New
- Bookmark
- Subscribe
- Mute
- Subscribe to RSS Feed
- Permalink
- Email to a Friend
- Notify Moderator

07-03-2013
02:34 PM

07-03-2013
02:34 PM

Try basing your initial guesses on the information available to you from the basic data and, where possible, the form of the equations. See attached.

Alan

Improved file attached

4 REPLIES 4

Highlighted
##

- Mark as New
- Bookmark
- Subscribe
- Mute
- Subscribe to RSS Feed
- Permalink
- Email to a Friend
- Notify Moderator

07-03-2013
01:38 PM

07-03-2013
01:38 PM

Re: rate constants determination: from ODEs solving v.s. numerical integration

For parameter estimation, all parameters should be scaled to the same order of magnitude and then you try guesses for two orders of magnitude surrounding the expected values.

Highlighted

- Mark as New
- Bookmark
- Subscribe
- Mute
- Subscribe to RSS Feed
- Permalink
- Email to a Friend
- Notify Moderator

07-03-2013
02:34 PM

07-03-2013
02:34 PM

Try basing your initial guesses on the information available to you from the basic data and, where possible, the form of the equations. See attached.

Alan

Improved file attached

Highlighted
##

- Mark as New
- Bookmark
- Subscribe
- Mute
- Subscribe to RSS Feed
- Permalink
- Email to a Friend
- Notify Moderator

07-03-2013
09:03 PM

07-03-2013
09:03 PM

Re: rate constants determination: from ODEs solving v.s. numerical integration

Thanks. Alan

your method is very interesting, time-saving and skillful than mine. I appreciate very much. Now I got three ways of solving the problem, lol.

by the way, I think your answer is more educational than practical, coz the reaction orders are all integers, and this is not very often in engineering practice.

thanks a lot again for your very detailed answer to my question

Highlighted
##

- Mark as New
- Bookmark
- Subscribe
- Mute
- Subscribe to RSS Feed
- Permalink
- Email to a Friend
- Notify Moderator

07-03-2013
03:46 PM

07-03-2013
03:46 PM

Re: rate constants determination: from ODEs solving v.s. numerical integration

Thanks to all for the interesting to the book http://twt.mpei.ac.ru/TTHB/New-Chem-Kin/En-Ru-book.html

Now we try to solve book's problems in Prime with units - mole, L, sec etc.

This discussion is very usefull for the second edition of the book.