Journal of Computational Physics

Videos of 1)Transonic flow past cylinder using WENO 5th-order scheme 2) Kelvin-Helmholtz instability with three stencil types (T1, T2, T3) using WENO 3rd and WENO 6th-order schemes 3) Blown-up version Kelvin-Helmholtz instability with three stencil types (T1, T2, T3) using WENO 3rd and WENO 6th-order schemes

TopImorphSmall is a stl triangulation. remesh_TopImorphSmall is a Python script performing remeshing (3D) of TopImorphSmall. convergenceMVC is a Python script performing the convergence test of section 3.3. mvcP2 is a Python script performing derivations of section 5.2.

This is the original data used in our paper "A robust and effcient segregated algorithm for fluid flow: The EPPL method", including the three test cases and the post-processing scripts to generate the figures. All the cases are in the format of OpenFOAM , any text editors are enough to view the settings, and paraview , tecplot and gnuplot are recommanded to view the fields. For more information about the settings, please have a look at our article.

These files contain numerical results for the scalar flux and current obtained from several orders of approximation in the post-processed P_N method described in the manuscript. They intend to show how the method converged for the three test cases studied in the manuscript.

Surface and volume mesh files for all computations shown in the manuscript

Fortran implementation of the perturbed truncated and shifted (PeTS) equation of state (Heier et al., Mol. Phys., 2018) for the Lennard-Jones truncated and shifted potential. The implementation is based on the reduced Helmholtz energy. It is possible to choose from a variety of input variables, e.g. density and temperature, density and pressure, etc. Only for density and temperature as input variables, the PeTS EOS can be directly evaluated. Otherwise, Newton algorithms are used to invert the EOS.

Data for the paper in the title submitted to JCP

See Appendix C.

This dataset contains code examples for different symplectic integrators with non-canonical quadrature points described in [1]. Here guiding-center motion is implemented in its axisymmetric variant for tokamak magnetic fields in canonicalized flux coordinates. [1] C. G. Albert, S. V. Kasilov, and W. Kernbichler, Symplectic integration with non-canonical quadrature for guiding-center orbits in magnetic confinement devices, Mar. 2019, arXiv:1903.06885. Submitted to J. Comp. Phys

frkg2.tgz is a C implementation of FRKG numerical integration at second order. The code consists of two files: the file frkgcore.c contains the code necessary to obtain a solution to the problem prescribed in frkguser.c. The files Tables_N1-N2/arks_<N>_<a>_<b>.csv contain the coefficents for FRKG schemes at orders N=1...2 and, at each order, over a range of 18 Gegenbauer coefficient values running from 0 to 2N. The coefficients for each scheme is provided for up to Nx257 stages. The default problem is similar to the Brusselator presented by Hairer & Wanner (see (1.6) of Sect. IV.1): a stiff nonlinear diffusion-reaction problem describing chemical kinetics of a tri-molecular chemical reaction subject to periodic boundary conditions. Advection is added following Abdulle & Vilmart (2013). The problem is specified via three user defined functions: deriv() - populates f(W)=W_t; spec() - returns spectral radius; and Wexact(). User parameters are: RTOL - error tolerence; NE - number of grid points (including boundary) assuming a uniform square mesh; TEND; XMAX; XMIN. To apply the FRKG scheme to other problems it should only be necessary to supply a modified version of frkguser.c. Note that complex variables are used via complex.h internally to the time-marching function trystep(). Above trystep() no complex variables are visible. Usage: to compile executable frkg using gcc gcc -O3 frkgcore.c frkguser.c -lm -o frkg Solution is written to sol.dat Reference solution is sol.ref To plot first species via gnuplot: splot "sol.dat" every 10:10 u 1:2:3 w pm3d ----------------------------- The additional files Tables_N3-N8/arks_<N>_<a>_<b>.csv contain the coefficents for FRKG schemes at orders N=3...8 and, at each order, over a range of 18 Gegenbauer coefficient values running from 0 to 2N. The coefficients for each scheme is provided for up to Nx257 stages.