Transformação de coordenadas por diferenças finitas para coordenadas polares esféricas

8

Eu tenho parte de um problema descrito pela equação de conservação do momento:

ρt+1sinθθ(ρusinθ)=0

Onde e ρ = f ( θ , t ) (velocidade constante).u=f(θ)ρ=f(θ,t)

Ingenuamente, pode-se aplicar uma das soluções listadas aqui . O problema em questão é melhor descrito em coordenadas polares esféricas (casca esférica fina) e nessas soluções em cartesiano. devo fazer algum tipo de transformação de coordenadas antes de discretizar essa equação ou posso discretizá-la diretamente?

Em segundo lugar, existe alguma razão para primeiro expandirmos a derivada em θ e depois tentarmos discretizar?

Como observação - eu fiz várias das opções acima e obtive soluções que não parecem consistentes (fisicamente, algumas parecem fazer sentido). Estou interessado em saber se há uma transformação de coordenadas adequada que deve ser feita ou se algum dos métodos mencionados anteriormente será suficiente.

EDITAR:

Eu defino fluxo como:

Φi+1/2=ui+1/2+|ui+1/2|2ρisinθi+ui+1/2|ui+1/2|2ρi+1sinθi+1

Φi1/2=ui1/2+|ui1/2|2ρi1sinθi1+ui1/2|ui1/2|2ρisinθi

sinθ±12

θ=0

Marm0t
fonte
2
hsinθsinθ

Respostas:

3

Φ

θuρ

Godric Seer
fonte
Obrigado pela resposta - investiguei possíveis soluções para esse problema, mas ainda não encontrei uma solução que pareça razoável. Eu tenho o BC's em baixa - mas não há benchmarks para testar meu esquema.
Marm0t
Φρsin(θ)u