MRI evp in cylindarical coordinate

[Luis]

Hi Dedalus Community!

I want to implement an eigenvalue problem to solve the linear growth rate of MRI in cylindrical coordinate. I saw a recent post regarding using disk basis for 3Dsimulation: https://groups.google.com/g/dedalus-users/c/a9MJHPjDFt8. I just begin implementing this and my starting point is to create a coordinate that is a direct product of d3.Coordinate('z') and d3.PolarCoordinates('phi', 'r'). I soon found that I cannot extract unit vector of this system using coords.unit_vector_fields(dist). I wonder does dedalus support periodic cylinders for now? If so, how should I initialized the normal vectors in cylindrical coordinate. I attached my code below for your reference.

Thank you!

Best

Hongke

[Keaton Burns]

Hi Hongke,

This is still under development, but is going well and we expect it to release working examples of periodic cylinder problems soon.  In the meantime, if you just want to do an EVP with some axial wavenumber kz (as opposed to a true 3d simulation), you can follow the pipe flow EVP example here:

https://dedalus-project.readthedocs.io/en/latest/pages/examples/evp_disk_pipe_flow.html

Best,

-Keaton

--
You received this message because you are subscribed to the Google Groups "Dedalus Users" group.
To unsubscribe from this group and stop receiving emails from it, send an email to dedalus-user...@googlegroups.com.
To view this discussion on the web visit https://groups.google.com/d/msgid/dedalus-users/61f0ae3a-8a37-47f9-99f6-0312b3ed6a79n%40googlegroups.com.