Я решаю множество PDE численно, но прикладная математика не является моей областью. Я не понял, какие журналы по прикладной математике я должен читать, чтобы не отставать от последних событий в этой области. Какие хорошие журналы нужно почитать, чтобы не отставать от последних достижений в...
Прямо сейчас у меня есть код, который использует алгоритм Кранка-Николсона, но я думаю, что я хотел бы перейти к алгоритму более высокого порядка для временного перехода. Я знаю, что алгоритм Кранка-Николсона стабилен в той области, в которой я хочу работать, но я обеспокоен тем, что некоторые...
Я знаю, что мы можем использовать методы математического анализа, чтобы доказать, есть ли у IVP или BVP решение, оно уникально и постоянно зависит от граничных / начальных значений. Для некоторых PDE, особенно нелинейных pde, очень трудно, если не невозможно, доказать правильность. Есть ли...
Пусть - выпуклая полигонально ограниченная липшицева область в , пусть .R 2 f ∈ L 2 ( Ω )ΩΩ\OmegaR2R2\mathbb R^2f∈L2(Ω)f∈L2(Ω)f \in L^2(\Omega) Тогда решение задачи Дирихле в , в имеет единственное решение в и корректно, т. Е. Для некоторой константы имеем .Δu=fΔu=f\Delta u =...
Я хотел бы узнать, как работает элемент Raviart-Thomas (RT). С этой целью я хотел бы аналитически описать, как базисные функции выглядят на контрольном квадрате. Цель здесь не в том, чтобы реализовать это самостоятельно, а в том, чтобы просто получить интуитивное понимание элемента. Я в...
Я пытаюсь решить двумерное уравнение Пуассона, используя метод прерывистого Галеркина (DG) и следующую дискретизацию (у меня есть файл png, но мне не разрешено загружать его, извините): Уравнение: ∇ ⋅ ( х ∇ Т) + f= 0∇⋅(κ∇T)+f=0\nabla \cdot( \kappa \nabla T) + f = 0 Новые уравнения: Q= κ ∇ T∇ ⋅ кв=...
Я решаю систему реакции-диффузии Тьюринга с помощью следующего кода C ++. Это слишком медленно: для текстуры 128x128 пикселей допустимое количество итераций составляет 200, что приводит к задержке в 2,5 секунды. Мне нужно 400 итераций, чтобы получить интересное изображение, но 5 секунд ожидания -...
Рассмотрим следующую задачу где форсирующий член может зависеть от (см. Формулировку 1 ниже для формулировки), а также от и его первых производных. Это 1 + 1 мерное волновое уравнение. У нас есть начальные данные, прописанные в .Wuv=FWuv=F W_{uv} = F u,vu,vu,vWWW{u+v=0}{u+v=0}\{u+v = 0\} Меня...
Я использую конечно-разностную схему Кранка-Николсона для решения одномерного уравнения теплопроводности. Мне интересно, если принцип максимума / минимума уравнения теплопроводности (то есть, что максимум / минимум возникает при начальном условии или на границах) также имеет место для дискретного...
Я знаю, что кусочно-линейное приближение конечных элементов uhuhu_h из Δu(x)=f(x)in Uu(x)=0on ∂UΔu(x)=f(x)in Uu(x)=0on ∂U \Delta u(x)=f(x)\quad\text{in }U\\ u(x)=0\quad\text{on }\partial U удовлетворяет ∥u−uh∥H10(U)≤Ch∥f∥L2(U)‖u−uh‖H01(U)≤Ch‖f‖L2(U) \|u-u_h\|_{H^1_0(U)}\leq Ch\|f\|_{L^2(U)} при...
Я слышал, что быстрое преобразование Фурье можно использовать для решения проблемы Пуассона, когда все граничные условия одного типа ... ряды синусов для дирихле, косинусы для Неймана и оба для периодических. Рассматривая двумерную прямоугольную область, предположим, что две противоположные стороны...
У меня есть пара вопросов относительно следующего: Я пытаюсь решить уравнение Шредингера в 1D, используя дискретизацию кривошипа Николсона с последующим инвертированием получающейся трехдиагональной матрицы. Моя проблема теперь превратилась в проблему с периодическими граничными условиями, и...
Я читаю статью [1], где они решают следующее нелинейное уравнение UT+UИкс+ тыUИкс-Uх х т= 0ut+ux+uux−uxxt=0\begin{equation} u_t + u_x + uu_x - u_{xxt} = 0 \end{equation} используя методы конечных разностей. Они также анализируют устойчивость схем с помощью анализа устойчивости фон Неймана. Однако,...
Я заинтересован в предложениях для книжных ссылок на предмет числовых PDE и ODE, в частности, строгий анализ таких методов в манере, написанной для профессиональных математиков. Он не должен быть чрезвычайно всеобъемлющим в смысле перечисления сотен или тысяч различных методов, но мне было бы...
В литературе FEM полувариационные методы обычно используются при решении зависящих от времени PDE. Я не видел полностью вариационного подхода, то есть, когда FEM дискретизирует пространство и время, возможно, позволяя использовать неструктурированные сетки пространства-времени. Хотя методы...
Я обнаружил, что метод линий - очень естественный способ думать о дискретизации PDE. Поэтому я всегда по умолчанию придерживаюсь этого подхода, когда мне предлагают новый набор уравнений. Я никогда не видел PDE, где это не будет работать. Мне интересно, есть ли методы дискретизации (или типы PDE),...
Я сделал странное наблюдение, вычисляя решение простого одномерного уравнения диффузии реакции: ∂∂Tа =∂2∂Икс2- б∂∂ta=∂2∂x2a−ab\frac{\partial}{\partial t}a=\frac{\partial^2}{\partial x^2}a-ab ∂∂Tb = - a b∂∂tb=−ab\frac{\partial}{\partial t}b=-ab ∂∂Tс = а∂∂tc=a\frac{\partial}{\partial t}c = a...
Предположим, у меня была следующая периодическая проблема 1D адвекции: ∂u∂t+c∂u∂x=0∂u∂t+c∂u∂x=0\frac{\partial u}{\partial t} + c\frac{\partial u}{\partial x} = 0 в Ω=[0,1]Ω=[0,1]\Omega=[0,1] u(0,t)=u(1,t)u(0,t)=u(1,t)u(0,t)=u(1,t) u(x,0)=g(x)u(x,0)=g(x)u(x,0)=g(x) где g(x)g(x)g(x) имеет разрыв...
Я погружаюсь в увлекательный мир анализа методом конечных элементов и хотел бы решить большую термомеханическую проблему (только термическая механическая, без обратной связи).→→\rightarrow Для механической задачи я уже понял из ответа Джеффа , что мне нужно будет использовать итеративный решатель...
Я использую набор инструментов Matlab pde для решения определенного эллиптического уравнения в 2D. Решение в порядке, хотя мне нужно построить его вдоль заданной линии, т.е. вырезать плоский срез из трехмерной сетки, представляющей решение. Я не могу придумать, как разумно задействовать функции...