Chemistry. Trying to find dissociation constant with MathCAD

[dtischenko]

Had some problems with calculating the concentration. Couldn't find the way to receive positive solution at the end. How it could be done?

[RichardJ]

Weight the "greater than" constraints very heavily, so rather than PO4>0, use PO4*10^13>0. However, if you add a constraint HPO4*10^13>0 it can't find a solution, so I think you have an error somewhere (perhaps a sign error?)

[ValeryOchkov]

It is a typical case when we must use not the Given-Fine block but programming - bisection method.

I have wrote it in the article (not H3PO4 but H2CO3)

Compare here please

http://twt.mpei.ac.ru/ochkov/CO2/CO2-prog.pdf

picture 1 and 2

[ValeryOchkov]

See please my old message - pH of a buffer water solution

[ViktorKorobov]

Дима, см. картинку. Приведенный мной вариант решения мало чувствителен  начальным приближениям.

[ValeryOchkov]

Я сам об этих логарифмах думал. Но больше о программной реализации.

У нас, Виктор Иванович, по-моему с этих логарифмов и общение началось.

А вот еще про логарифмы - Nostalgia - a slide rule

Пришлите, плз, сам файл с расчётом. Я его хочу в Прайме с единицами измерения прикрутить. Там, кстати, в последнем уравнении, как мне кажется, с.0 нужно убрать.

И второе - вы получили письмо из маплесофт?

[ViktorKorobov]

Да, Валерий Федорович, наше активное сотрудничество стало возможным именно благодаря этим логарифмам!

Прикрепляю расчетные файлы. Я совсем забыл, что химическую символику нельзя использовать в решающих блоках Mathcad 14, 15. Файл работать будет, но при попытке его сохранения появится ошибка:

Если проигнорировать это сообщение, то в дальнейшем такой файл уже открыть не удастся. Поэтому я отдельно прикрепляю файл без химической нотации.

По поводу C.0 в последнем уравнении. У меня ошибки нет. Я специально сверился с подобной задачей из Дж. Батлер. Ионные равновесия. - Л.:Химия, 1973.

Да, только что увидел письмо из маплесофт.

Виктор

[ValeryOchkov]

I hope it will be in Prime 4.

Now we can use S-C-k in Prime thru Mathcad 15 converting

but S-C-j - not - we will get a picture.

But we can use in Prime more than one index.

[ViktorKorobov]

Unfortunately we can not use in Prime 3 superscript to denote the charge of the ions ...

[ValeryOchkov]

Viktor Korobov написал(а):

Да, Валерий Федорович, наше активное сотрудничество стало возможным именно благодаря этим логарифмам!

Прикрепляю расчетные файлы. Я совсем забыл, что химическую символику нельзя использовать в решающих блоках Mathcad 14, 15. Файл работать будет, но при попытке его сохранения появится ошибка:

Если проигнорировать это сообщение, то в дальнейшем такой файл уже открыть не удастся. Поэтому я отдельно прикрепляю файл без химической нотации.

По поводу C.0 в последнем уравнении. У меня ошибки нет. Я специально сверился с подобной задачей из Дж. Батлер. Ионные равновесия. - Л.:Химия, 1973.

Я думаю лишнее - это уравнение электронейтральности, а С.0 - это сумма фосфатов. У Батлера там еще и натрий.

Я эту проблему хорошо знаю - в энергетике фосфаты вводятся в барабан котла для борьбы с накипью. И еще аммиак. Все это держится на уровне рН = 10.2.

[ViktorKorobov]

В таком случае нужно спросить у автора исходной задачи - что он имел ввиду под C.0?

Если аналитическую концентрацию H3PO4 - то тогда уравнение электронейтральности не нужно. Но тогда возникает вопрос корректности расчетов из-за того, что не учитываются коэффициенты активности из-за достаточно большого значения C.0.

Я все-таки думаю, что C.0 - это концентрация добавленного фосфата.