Kelvin-Helmholtz instability - Discontinuous Galerkin hydrodynamics
Ғылым және технология
2D Simulation of a Kelvin-Helmholtz instability with 4th order discontinuous Galerkin (DG) and adaptive mesh refinement. The simulation starts with 64^2 cells and is refined down to an effective resolution of 4096^2 cells. Shown is the surface density of the fluid. DG offers several advantages over traditional finite volume (FV) methods.It directly solves also for the higher-order moments of the solution, no reconstruction is needed, resulting in an inherent conservation of angular momentum and less advection and diffusion errors compared to a FV method. Furthermore, DG is a higher-order method with a small stencil and many local computations, which renders it highly suitable for high performance computing on massively parallel systems.
You may find the corresponding publication on arXiv: arxiv.org/abs/1506.06140
Пікірлер: 59
Thank you very much for sharing this stunning visualization! I will use this video to explain turbulence scales to my students. Your work is greatly appreciated
Fantastic work! Thanks for making the article public!
Beautiful work!
This is so amazing!
what a beautiful work :D
This is amazing!
Looks beautiful. What was the initial pertubation applied on the system?
oh my god this is literally a fractal
@dwightlumain8645
5 ай бұрын
That's what i'm thinking!
@StojanBarbaric
5 ай бұрын
Hell yeah it is, nature is all about fractal behaviour
@sca4723
5 ай бұрын
Not as that, fractal are between dimentions not in the linear, or in the plane nor the 3d space, something between them, thats because the sum and división fot their sides and nodes aren't integers
@petterlarsson7257
5 ай бұрын
@@sca4723 no. you are incorrect.
@pabilbadoespecial
3 ай бұрын
Yes, this is a feature of turbulence. It's called the Kolmogorov cascade where you will develop a power-law spectrum of modes. This is a lot of talking but it basically means that if you zoom in the image it should look vaguely the same, i.e. fractal :)
Just a question, what type of mesh is this? How can a cell have more than one neighbor for each side?
What is the Reynolds number here?
nice channel! subscribed
Imagine turning this into an awesome wallpaper
Amazing! Love those fractals. What mesh resolution did you use? It's a structured mesh, right?
@m31coding
8 жыл бұрын
Hi Gaurav! It's a structured adaptive mesh with an initial resolution of 64^2 cells. During the simulation the mesh is refined where the hydrodynamic gradients are large, resulting in an effective resolution of 4096^2.
Hi bro I have a question maybe you know how it can be realized in 3d software like blender and etc.
Hey Kevin, Great work! I had struggled with KH myself for a while. :) How did you visualize this? The visualization is as awesome as the calculations themselves!
Beautiful!!! This colour, haha, are you from swden?
Congrats! What software did you used?
Any tips on adaptive mesh refinement? I need to implement an adaptive mesh code and any references would be greatly appreciated.
@m31coding
7 жыл бұрын
Hi! The code paper of the RAMSES code is a great reference: arxiv.org/abs/astro-ph/0111367 best, Kevin
Great visualization! so detailed! Is it better compared to 2D AREPO simulation with similar resolution? and how much memory do you store per one AMR cell?
@m31coding
9 жыл бұрын
Hey Denis, good to hear from you! I haven't done comparisons with the moving mesh code, but it will depend on the problem at hand. In 2D the number of base functions is 1/2*p*(p+1), where p is the order of the scheme. Hence for fourth order DG there are 10 weights stored per conserved variable in every cell. The paper will be online on Tuesday :) See you soon hopefully!
Could you give us a comparison about Discontinuous Galerien Method and a Multidomain Decomposition Chebyshev Method?
@m31coding
3 жыл бұрын
Hi Karl, unfortunately I haven't worked on the method you suggested. Best regards!
Hi, I am interested in using this simulation if you allow me and with due credits of course in a fulldome planetarium movie about chaos. Do you have a 4096x4096 video of it?
@m31coding
3 жыл бұрын
Hi Julien, thank you for your request! Unfortunately I only have a 1920x1080 version, feel free to use it. A fulldome movie about chaos sounds great, have fun! Best regards!
What did you use for the visualization?
@m31coding
8 жыл бұрын
Hi, thank you for your question. The answer is I wrote my own Python scripts using numpy and matplotlib. After creating the single frames I combined them to a movie with avconv.
@kerch00
2 жыл бұрын
@@m31coding How many seconds(about) does it take to calculate one of these frames in Python?
What software did you use to render this?
@m31coding
7 ай бұрын
Thank you for your question. I created the individual frames with Python (Matplotlib) and combined them with ffmpeg to a video.
looks like you are simulating jupiter
Just curious. Why does initially sharp interface become smooth? Do you approximate equations with viscosity or this is due to the numerical viscosity?
@m31coding
4 ай бұрын
Thank you for your question! Yes, I believe this can be attributed to numerical viscosity. The Euler equations, which the code solves, do not account for viscosity.
@mkraposhin
4 ай бұрын
@@m31coding Thank you. Do you have any idea how to decrease it?
And that’s how Jupiters are made
Hey Kevin ! Great Visualization ! Did you use any numerical library/software!!
@m31coding
Жыл бұрын
Hi, thanks for your comment! I implemented DG into the Arepo Code from scratch. In the meantime, Arepo is publicly available, you may check it out. Best regards!
@arthsojitra53
Жыл бұрын
@@m31coding Thanks for the info !
or it could be black holes growing
this is a tsunami thrower that works using a tsunami
BOKEEE
галеркин! 😃
fractals
It looks like fractals
video starts *Ukrainian anthem starts playing but is soon cut off*
extrim
Wtf art
Меандр
heehee ukraine flag
ukraine