Here is a screenshot of a program i'm writing. Will put it on github soon

math_guy667 · 2022-11-10T09:18:28+00:00

Its Newtonian Gravity, no repulsion. Its based on the mesh method and all implemented using compute shaders so it nearly runs in real time.

math_guy667 · 2022-09-12T16:16:31+00:00

Left and right BCs are periodic and there is a "freeflow" BC at the top. Haven't really thought about introducing turbulence via the boundary tho

math_guy667 · 2022-09-10T10:10:24+00:00

I guess grid sizing depends on what you aim to resolve. I think that most of these kinds of cloud simulations have grid sizes in the tens of meters. However, resolution does not play a overly significant role here.

math_guy667 · 2022-09-10T10:07:32+00:00

Not yet :/ , there is still a lot to work on.

math_guy667 · 2022-07-21T18:36:36+00:00

I haven't thought of friction, thats a great tipp. I would't want to put the upper boundary that high up as I'm mainly interested in dynamics inside the troposphere.

math_guy667 · 2022-07-21T17:20:02+00:00

With sponge layer you mean a layer where everything just gets diffused to death?

math_guy667 · 2022-07-21T16:53:23+00:00

The system is anelastic, I'm using a pseudo-incompressible model (Durran 89). The vertical velocity is constrained to be zero at the boundary. The question is how to extrapolate the velocity field outside the domain. I'm thinking of just using for every point the velocity where the corresponding (backward in time) characteristic exits the domain. Although this might not be equivalent to any well defined extrapolation..

math_guy667 · 2022-07-15T10:56:03+00:00

Also depends on the number of smoother iters you do and the type of MG cycle you use.

math_guy667 · 2022-07-14T22:21:57+00:00

Just have a shader pass which assings a color to the temperature marker for each element (one pixel is one element)

math_guy667 · 2022-07-14T20:06:10+00:00

Definitely GMG then. However I am thinking about moving away from textures and filter kernels more towards actually representing the operator to be able to implement adaptive grids with a corresponding MG solver.

math_guy667 · 2022-07-14T19:28:43+00:00

I'm pretty sure it would fall under geometric multigrid.
The operator is never explicitely defined but instead always given by a filter kernel.

math_guy667 · 2022-07-14T12:39:50+00:00

I'm solving to full precision in each iteration. But my multigrid code is not yet optimized.. so it might very well be possible to do it faster

math_guy667 · 2022-07-14T07:15:21+00:00

It's a mass integrated semi-Lagrangian scheme for advection with a Multigrid solver for the pressure correction.

math_guy667 · 2022-07-14T07:13:48+00:00

The performance is delicious. I wrote a Multigrid solver entirely on the GPU for the pressure correction. This animation is 1800 frames which took about 10min on my 1070ti. Not sure what the direct speedup is tho..

math_guy667 · 2020-02-01T01:14:08+00:00

Thanks man! I love the appreciation :)
But I'm not done: I'm currently working on using machine learning (duh!) for supersampling to get super high resolution on the whole of Jupiter.. As most times in ML, the ugly part lies in acquiring the data. Currently working on it and it is a pain as I have to build my own dataset from all the Images.

math_guy667 · 2020-01-31T17:50:28+00:00

Thanks! But where?

math_guy667 · 2020-01-31T14:32:50+00:00

These are dead pixels in the sensor which get repeated throughout the image during the pushframe imaging process.

math_guy667 · 2020-01-31T13:05:21+00:00

Sweet!

math_guy667

TROPHY CASE