Considere o seguinte problema que o termo forçado pode depender de (consulte o Edit 1 abaixo para a formulação) e e seus primeiros derivados. Esta é uma equação de onda dimensional 1 + 1. Temos dados iniciais prescritos em .Wuv=FWuv=F W_{uv} = F u,vu,vu,vWWW{u+v=0}{u+v=0}\{u+v = 0\} Estou...
Estou usando o esquema de diferenças finitas da Crank-Nicolson para resolver uma equação de calor 1D. Gostaria de saber se o princípio máximo / mínimo da equação do calor (ou seja, que o máximo / mínimo ocorre na condição inicial ou nos limites) também vale para a solução discretizada. Isso...
Estou lendo um papel [1] onde eles resolver o seguinte não-linear equação usando métodos de diferenças finitas. Eles também analisam a estabilidade dos esquemas usando a análise de estabilidade de Von Neumann. No entanto, como os autores percebem, isso é aplicável apenas aos PDE lineares....
Antecedentes: Eu criei apenas uma solução numérica para o 2d Navier-Stokes, para um curso. Era uma solução para o fluxo da cavidade acionado pela tampa. O curso, no entanto, discutiu vários esquemas para discretizações espaciais e discretizações de tempo. Também fiz mais cursos de manipulação de...
Estou aplicando o método das diferenças finitas em um sistema de 3 equações acopladas. Duas das equações não são acopladas, no entanto, a terceira equação se une às outras duas. Notei que, alterando a ordem das equações, digamos de para que a matriz do coeficiente se torna simétrica.( x , z , y )(...
Como aproximar o número de condição de uma matriz grande GGG , se GGG é uma combinação de transformadas de Fourier (não uniforme ou uniforme), diferenças finitas e matrizes diagonais ?FFFRRRSSS As matrizes são muito grandes e não são armazenadas na memória e estão disponíveis apenas como...
Suponha que eu tenha o seguinte problema periódico de advecção 1D: ∂você∂t+ c ∂você∂x= 0∂você∂t+c∂você∂x=0 0\frac{\partial u}{\partial t} + c\frac{\partial u}{\partial x} = 0 in onde g (x) tem uma descontinuidade de salto em x ^ * \ in (0,1) . Ω = [ 0 , 1 ]Ω=[0 0,1 1]\Omega=[0,1] u ( 0 , t ) = u...
Quando é preferível usar os polinômios de Bernstein para aproximar uma função contínua, em vez de usar os únicos métodos preliminares de Análise Numérica: "Polinômios de Lagrange", "Operadores de diferenças finitas simples". A questão é sobre comparar esses
você pode me dar um link para uma explicação boa e simples do esquema de diferenças finitas de Shortley-Weller? Tentei pesquisar no Google, mas tudo o que recebo são publicações acadêmicas (inacessíveis). Também tentei ler o capítulo dedicado (4.8) do livro de Wolfgang Hackbusch "Equações...
Suponha que você tenha um sistema de PDEs para resolver. Pelo menos por simplicidade, vamos supor que seja independente do tempo, quase linear (linear em suas derivadas) resolvido em uma grade retangular no espaço (x, y) e com condições de contorno especificadas ao redor. Minha pergunta é mais...
Estou tentando modelar a condução de calor dentro de um cilindro de madeira usando métodos implícitos de diferenças finitas. A equação geral de calor que estou usando para formas cilíndricas e esféricas é: Onde p é o fator de forma, p = 1 para cilindro ep = 2 para esfera. As condições de...
Meu desafio é resolver o seguinte sistema de equações, que descreve a combustão de gases em meios porosos: 1) Continuidade ε∂ρg∂t+∂∂x(ρgux)=0ε∂ρg∂t+∂∂x(ρgux)=0\varepsilon \frac{\partial \rho_g}{\partial t} +\frac{\partial}{\partial x} \left(\rho_g u_x\right)=0 2) Lei de Darcy...
Eu tenho parte de um problema descrito pela equação de conservação do momento: ∂ρ∂t+ 1pecadoθ∂∂θ( ρ u sinθ ) = 0∂ρ∂t+1sinθ∂∂θ(ρusinθ)=0\frac{\partial \rho}{\partial t} + \frac{1}{\sin\theta} \frac{\partial}{\partial \theta}(\rho u \sin \theta) =0 Onde e ρ = f ( θ , t ) (velocidade constante).u...
O cálculo dos modos próprios de uma membrana semicircular se reduz ao seguinte problema de valor próprio ∇2u=k2u,∇2u=k2u,\nabla^2u=k^2u\;, onde a região de interesse é um semicírculo definido por e φ ∈ [ 0 , π ] .r ∈ [ 0 , 1 ]r∈[0,1]r\in[0,1]& Phi; ∈ [ 0 , π]φ∈[0,π]\varphi\in[0,\pi] É...
As equações de movimento para um sólido elástico são dadas por ∇⋅σ+f=ρu¨σ=Cεε=12(∇u+[∇u]T)∇⋅σ+f=ρu¨σ=Cεε=12(∇u+[∇u]T)\begin{align} &\nabla \cdot \boldsymbol{\sigma} + \mathbf{f} = \rho \ddot{\mathbf{u}}\\ &\boldsymbol{\sigma} = \mathbb{C}\boldsymbol{\varepsilon}\\ &\boldsymbol{\varepsilon} =...
Eu gostaria de saber o que é mais vantajoso quando se trata de resolver equações hiperbólicas não lineares, métodos de elementos finitos ou diferenças finitas? Qual método será melhor na captura de choques? É possível fornecer uma resposta / referências detalhadas? Além disso, quero resolver...