movingwindow.inp

pic 模拟程序！面向对象

MovingWindow
{

  This input file shows an example of a beam moving in a uniform background.
the geometry is r-z, cylindrical.  All the edges are bounded with conductors:
these boundaries are necessary.

  The initial load consists of stationary deuterium along with protons
moving to the right with v1drift = 3.e9.  This is in units of EV, so I
(JRC) assume that v1drift is actually the momentum in the 1-direction.

  A special load named shiftLoad (its parameter Name = shiftLoad) is executed
whenever a shift of the window occurs.  This special load in this example
loads a uniform plasma into the last cell, which is initially empty when
the moving window shifts.

  The movingWindow has a hard-coded delay time before it starts shifting, 
something like 1/2 or 2/3 of a lightspeed transit time across the system.

  Also, the moving window motion is hard-coded to move to the right, so
rightward-moving particles will stay in the window.  Everything else will
seem to move left.

}
Variables
{
	JMAX = 20
	KMAX = 20
	XMAX = 0.1
	RMAX = 0.02
	DX = XMAX/JMAX
}
Region
{
Grid
{
	J = JMAX
	x1s = 0.0
	x1f = XMAX
	n1 = 1.0
	K = KMAX
	x2s = 0.0
	x2f = RMAX
	n2 = 1.0
}
Control
{
	dt = 2.8E-12
	movingWindow = 1
        //shiftDelayTime =  1e-12 * 10        
        shiftDelayTime = 0
}
Species
{
		  name = proton
		  m = 1.8E-27
		  q = 1.6e-19
}
Species
{
		  name = deuterium
		  m = 3.6E-27
		  q = 1.6e-19
}
Conductor
{
	j1 = 0
	k1 = 0
	j2 = 0
	k2 = 10
	normal = 1
}
Conductor
{
	j1 = 0
	k1 = 10
	j2 = 0
	k2 = 20
	normal = 1
}
Conductor
{
	j1 = 0
	k1 = 20
	j2 = 20
	k2 = 20
	normal = -1
}
Conductor
{
	j1 = 20
	k1 = 20
	j2 = 20
	k2 = 0
	normal = -1
}
Load
{
	x1MinMKS = 0.025
	x1MaxMKS = 0.075
	x2MinMKS = 0.0
	x2MaxMKS = 0.01
	speciesName = proton
	density = 1.0e17
	np2c = 3e8
	units = EV
	v1drift = 3.0e9
}
Load
{
	// initial load
	// i fill everything but the last cell
	x1MinMKS = 0.0
	x1MaxMKS = XMAX - DX
	x2MinMKS = 0.0
	x2MaxMKS = RMAX
	speciesName = deuterium
	density = 1.0e17
	np2c = 3e9
	units = EV
}
Load
{
	// initial load and shift load
	// this fills the last cell when shifts are done, and at t=0
	Name = shiftLoad
	x1MinMKS = XMAX - DX
	x1MaxMKS = XMAX
	x2MinMKS = 0.0
	x2MaxMKS = RMAX
	speciesName = deuterium
	density = 1.0e17
	np2c = 3e9
	units = EV
}

CylindricalAxis
{
	j1 = 0
	k1 = 0
	j2 = 20
	k2 = 0
	normal = 1
}
}