Vamos considerar uma condição inicial suave e a equação do calor em uma dimensão: no intervalo aberto , e vamos assumir que queremos resolvê-lo numericamente com diferenças finitas.∂tu=∂xxu∂tu=∂xxu \partial_t u = \partial_{xx} u]0,1[]0,1[]0,1[ Eu sei que para o meu problema a ser bem colocada,...
Estou tentando aprender sobre resolver numericamente o PDE sozinho. Estou começando com o método da diferença finita (FDM) há algum tempo, porque ouvi dizer que o FDM é o fundamento de vários métodos numéricos para o PDE. Até agora, eu tenho um entendimento básico do FDM e fui capaz de escrever...
Eu implementei um solucionador de Euler para trás no python 3 (usando numpy). Para minha própria conveniência e como exercício, também escrevi uma pequena função que calcula uma aproximação de diferença finita do gradiente para que nem sempre seja necessário determinar analiticamente o jacobiano...
Vamos começar com um problema do formulário (L+k2)u=0(L+k2)u=0(\mathcal{L} + k^2) u=0 com um conjunto de determinadas condições de contorno ( Dirichlet , Neumann , Robin , Periodic , Bloch-Periodic ). Isso corresponde à localização dos valores próprios e dos vetores próprios para algum operador ,...
Gostaria de escrever meu próprio solucionador para equações de Euler compressíveis e, mais importante, quero que ele funcione de maneira robusta em todas as situações. Eu gostaria que fosse baseado em FE (DG está ok). Quais são os métodos possíveis? Estou ciente de fazer a DG de ordem 0 (volumes...
Estou trabalhando na resolução das equações de poroelasticidade unidimensional acopladas (modelo de biot), dadas como: −(λ+2μ)∂2u∂x2+∂p∂x=0−(λ+2μ)∂2u∂x2+∂p∂x=0-(\lambda+ 2\mu) \frac{\partial^2 u}{\partial x^2} + \frac{\partial p}{\partial x} =
Se empregarmos o Método de Linhas para discretização (discretização separada de tempo e espaço) de EDPs hiperbólicas, obtemos após a discretização espacial por nosso método numérico favorito (fx. Método do Volume Finito) na prática importa qual solucionador de ODE empregamos para a discretização...
Termos de origem, como os devidos à batimetria nas equações de águas rasas, precisam ser integrados de uma maneira especial para preservar os estados físicos estáveis. Existe uma maneira geral de construir métodos bem equilibrados ou requer técnicas especiais para cada
A forma forte de um PDE exige que a solução desconhecida pertença a . Mas a forma fraca requer apenas que a solução desconhecida pertença a .H2H2H^2H1H1H^1 Como você reconcilia isso?
Estou tentando simular modelos básicos de semicondutores para fins pedagógicos - a partir do modelo de difusão à deriva. Embora eu não queira usar um simulador de semicondutores pronto para uso - aprenderei outros modelos (comuns, recentes ou obscuros), quero usar um solucionador de PDE pronto para...
Para usar a Fast Fourier Transform (FFT) em dados de amostra uniforme, por exemplo, em conexão com os solucionadores de PDE, é sabido que a FFT é um algoritmo ). Quão bem a escala FFT é processada em paralelo para (isto é, muito grande)?n → ∞O (nlog( N )O(nlog(n)\mathcal{O}(n\log(n)n →...
Suponha que eu tenha um problema de valor limite: d2vocêdx2+ dvdx= f em Ωd2udx2+dvdx=f in Ω\frac{d^2u}{dx^2} + \frac{dv}{dx}=f \text{ in } \Omega dvocêdx+ d2vdx2= g em Ωdudx+d2vdx2=g in Ω\frac{du}{dx} +\frac{d^2v}{dx^2} =g \text{ in } \Omega u = h em ∂Ωu=h in ∂Ωu=h \text{ in }...
Os altos fluxos de número de Reynolds produzem camadas de limite muito finas. Se a resolução de parede for usada na simulação de redemoinho grande, a proporção da imagem pode ser da ordem de 10610610^6 . Muitos métodos se tornam instáveis nesse regime porque a constante inf-sup se degrada como...
Para algum domínio convexo simples em 2D, temos alguns u ( x ) que satisfazem a seguinte equação: - d i v ( A ∇ u ) + c u n = f com certas condições de contorno de Dirichlet e / ou Neumann. Que eu saiba, aplicar o método de Newton em um espaço de elementos finitos seria uma maneira relativamente...
Estou tentando entender como o método de otimização baseado em adjuntos funciona para uma otimização restrita do PDE. Particularmente, estou tentando entender por que o método adjunta é mais eficiente para problemas em que o número de variáveis de design é grande, mas o "número de equações é...
Estou implementando o artigo " Transporte de massa ideal para registro e entortamento ", meu objetivo é colocá-lo online, pois não consigo encontrar nenhum código de transporte de massa euleriano on-line e isso seria interessante pelo menos para a comunidade de pesquisa em processamento de...
Eu gostaria de resolver alguns PDEs em variedades, digamos, por exemplo, uma equação elíptica em uma esfera. Por onde começo? Gostaria de encontrar algo que use bibliotecas / códigos preexistentes em 2d, nada tão sofisticado (no momento) Adicionado mais tarde: Artigos e relatórios são...
Qual é o estado da arte atual para resolver EDPs parabólicas de maior dimensão (3-10) no domínio complexo com pólos simples (da forma ) e absorvendo condições de contorno?1|r⃗ 1−r⃗ 2|1|r→1−r→2| \frac{1}{|\vec{r}_1 - \vec{r}_2|} Especificamente, estou interessado em resolver a equação de...
Alguém pode me ajudar a encontrar os livros sobre soluções numéricas (diferença finita e métodos Crank-Nicolson) das equações de Poisson e difusão, incluindo exemplos de geometria irregular, como um domínio que consiste na área entre um retângulo e um círculo (especialmente livros ou links em...
Sou iniciante na FE. Minha aplicação é a precificação de derivativos financeiros em que o espaço é tridimensional. Então, adicionando tempo, o problema tem seis dimensões. Eu tentei olhar em volta (Fenics, escript, deal.II, ...), mas meu entendimento é que esses softwares são limitados a 3 + 1...