Variance Reduced DSMC in 3D with 1D Boundary conditions.
This is code that was used in my PhD Thesis
There is sample input in the input.txt file.
I'm not sure if this is SOTA anymore but including for reference.
This java is very old and designed to run as fast as possible on machines circa 2008.
If you are interested in this topic in any serious way I suggest you look here