Если возникла потребность в имитации бурной деятельности, то можно сделать так: настаиваем PCG solver на 500 outer и 500 inner итераций, критерий сходимости — 1e-5. Ну и модельку в нестационар на пять сотен шагов по времени. Все, комп на пару часов будет занят, можно халявить.
А еще можно PEST запустить с этой же моделькой, но тут можно и на недельку «отдохнуть».
Прикладная гидрогеология, геофильтрационное моделирование, инженерная геология. Блог специалиста по гидрогеологическому моделированию Копылова Александра.
27/07/2011
Для халявщиков
20/07/2011
Processing MODFLOW. Введение.
В мои должностные обязанности в свое свое время входило обучение основам моделирования новых сотрудников и коллег из сторонних организаций. Многие из них в впоследствии вполне успешно продолжили деятельность на этом поприще, так что мою методу вполне можно считать состоявшейся, хотя, несомненно, не лишенной недостатков. Комментарии, а особенно конструктивная критика, категорически приветствуются.
Описание работы с препроцессором (о сути этого названия я писал раньше) у меня практически полностью совпадает со структурой программного меню, которая, в свою очередь, соответствует последовательности создания модели. За эту особенность я и люблю использовать эту программу в качестве учебного пособия. Собственно говоря, англоязычная инструкция к программе написана так-же – именно как последовательное описание пунктов меню.
Скачать бесплатную версию Processing Modflow - там, правда, доступна достаточно древняя версия 5.3.1, но она хотя бы бесплатная. Более новая 8-ка вообще неработоспособна (там в коде перепутаны идентификаторы полей ввода в некоторых диалогах, что приводит к совершенно невообразимым глюкам). Лично мне нравится семерка, но она платная. Я буду описывать именно седьмую версию – мне этот подход кажется наиболее универсальным. Если что-то будет непонятно, то я всегда рад общению.
05/07/2011
Профильные модели
Все-таки есть небольшая польза от тупого просиживания штанов в офисе. Придумал способ, как избежать длительного и мучительного итеративного рисования свободной поверхности на профильных моделях в MODFLOW.
Нужно в качестве подошвы ячеек задать координату y их нижней границы. Следовательно, когда уровень в блоке становится ниже подошвы, то ячейка сама отключается, что нам и нужно. Не стоит забывать о сохранении единичной мощности пласта.
Остаются проблемы с инфильтрацией: не совсем понятно, куда ее задать (на всю «площадь» некорректно, да и приводит к тому же к необходимости итеративного же ее уточнения), а при неизвестной верхней границе — невозможно задать на ней инфильтрацию.
Ну и обязательно надо включить модуль wetting capability, что традиционно приводит к проблемам со сходимостью.