using Plasma
e = species.e
function TwoStream(vs1, vs2)
function P(x,v)
if !(v isa Array)
v = [v]
end
if !(x isa Array)
x = [x]
end
v = sqrt(sum(v .^2))
x = sqrt(sum(x .^2))
1/2*(1/(sqrt(2*π)*vs1)*exp(-(v-vs2)^2)/(2*vs1^2) + 1/(sqrt(2*π)*vs1)*exp(-(v+vs2)^2)/(2*vs1^2))
end
end
D_e = Distribution(TwoStream(1.6, -1.4), e)
G = Geometry()
plasma = ElectrostaticPlasma([D_e], G)
sol = Plasma.solve(plasma, dim=1, GPU=false, time_ub = 4.0, ub=4.0, lb=-4.0)
-
Notifications
You must be signed in to change notification settings - Fork 0
Benchmarks for Plasma.jl and neural Vlasov differential equation solvers
License
killah-t-cell/PlasmaBenchmarks.jl
Folders and files
Name | Name | Last commit message | Last commit date | |
---|---|---|---|---|
Repository files navigation
About
Benchmarks for Plasma.jl and neural Vlasov differential equation solvers
Resources
License
Stars
Watchers
Forks
Releases
No releases published
Packages 0
No packages published