liangwang0734

-- Gkyl ------------------------------------------------------------------------
--
-- Test for five-moment source updater
--    _______     ___
-- + 6 @ |||| # P ||| +
--------------------------------------------------------------------------------

local ffi  = require "ffi"
local Unit = require "Unit"
local Grid = require "Grid"

liangwang0734

-- GEM-challenge problem

Pi = Lucee.Pi
log = Lucee.logInfo

-- physical parameters
gasGamma = 5./3.
elcCharge = -1.0
ionCharge = 1.0
ionMass = 1.0

liangwang0734

-- Gkyl ------------------------------------------------------------------------
local Moments = require("App.PlasmaOnCartGrid").Moments
local Euler = require "Eq.Euler"

-- physical parameters
gasGamma = 5./3.
elcCharge = -1.0
ionCharge = 1.0
ionMass = 1.0
elcMass = ionMass/25.0

liangwang0734

__all__ = ["gkylFile", "gkylFile_Moment"]

import numpy as np
import adios

import numexpr as ne
import re
import random
import string

liangwang0734

-- XXX testing only; force to use 1st order Lax solver (positivity-conserving but diffusive)
alwaysUseLaxFlux = true
------------------------------------
-- 2D 10-MOMENT SIMULATION OF KHI --
------------------------------------

----------------------
-- HELPER FUNCTIONS --
----------------------

liangwang0734