12.07.2015 Views

Formacion de estructura con campos escalares: la ... - Cinvestav

Formacion de estructura con campos escalares: la ... - Cinvestav

Formacion de estructura con campos escalares: la ... - Cinvestav

SHOW MORE
SHOW LESS
  • No tags were found...

You also want an ePaper? Increase the reach of your titles

YUMPU automatically turns print PDFs into web optimized ePapers that Google loves.

Capítulo 5Hidrodinámica: solución numéricaEn este capítulo presentamos los resultados numéricos <strong>de</strong> <strong>la</strong> forma hidrodinámica <strong>de</strong><strong>la</strong> ecuación <strong>de</strong> tipo Schrödinger, en <strong>la</strong> cual no se <strong>con</strong>si<strong>de</strong>ró el potencial cuántico. La razón<strong>de</strong> ello es por su <strong>de</strong>pen<strong>de</strong>ncia cuadrática <strong>de</strong> <strong>la</strong> <strong>con</strong>stante <strong>de</strong> P<strong>la</strong>nck, lo cual hace razonabledicha <strong>con</strong>si<strong>de</strong>ración. En esta primera aproximación, se resolvieron tres casos <strong>de</strong> interés:<strong>la</strong> ecuación <strong>de</strong> advección (U = 0), <strong>la</strong> ecuación <strong>de</strong> Burgers (U = 0) y por último el casoU = kx 2 /2. Para el primero, a<strong>de</strong>más <strong>de</strong> ser <strong>la</strong> ecuación diferencial lineal más simple <strong>de</strong> tipohiperbólico, ésta tiene solución analítica y figura a <strong>la</strong> partícu<strong>la</strong> libre en mecánica cuántica.La ecuación <strong>de</strong> Burgers es <strong>la</strong> ecuación diferencial no lineal más sencil<strong>la</strong> <strong>de</strong> naturalezahiperbólica, <strong>la</strong> cual se obtiene <strong>de</strong> <strong>la</strong> ecuación <strong>de</strong> momentos <strong>de</strong>l sistema hidrodinámicopara U = 0. En el caso U = kx 2 /2, nuestro interés es ver qué suce<strong>de</strong> cuando le damos alsistema <strong>con</strong>diciones iniciales cercanas al <strong>de</strong>l estado base <strong>de</strong>l osci<strong>la</strong>dor armónico cuántico,es <strong>de</strong>cir, queremos ver qué pasa cuando <strong>de</strong>spreciamos el término cuántico. Finalmente, sehicieron <strong>la</strong>s pruebas <strong>de</strong> <strong>con</strong>vergencia que nos darán credibilidad en los resultados.5.1. Código numéricoPara realizar los casos mencionados, hemos e<strong>la</strong>borado un código <strong>de</strong> primer or<strong>de</strong>n tantoen <strong>la</strong> discretización espacial como en <strong>la</strong> temporal que resuelve el sistema∂ t ρ + ∂ x (ρv) = 0 , (5.1)∂ t v + ∂ x (v 2 /2) = −αkx , (5.2)don<strong>de</strong> es evi<strong>de</strong>nte que ambas expresiones están en forma <strong>con</strong>servativa, siendo f ρ = ρv elflujo en <strong>la</strong> ecuación <strong>de</strong> <strong>con</strong>tinuidad y f v = v 2 /2 y S v = kx el flujo y fuente en <strong>la</strong> ecuación35

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!