orionwp.inp

pic 模拟程序！面向对象

OrionWP
{

 - An electromagnetic wave with a Gaussian profile is launched from the 
   left boundary in Cartesian geometry for a short pulse in an 
   uniform plasma.

 - The wave is defined by the variables:
     - gauss_w      = spot size at beam waist. 
     - gauss_rl     = Rayleigh length.
     - omegag       = frequency.
     - zf           = position of focus. 
     - Agw          = wave amplitude.    
}

Variables
{
// Grid variables

      Nx = 1600 
      Ny = 150 
      Lx = 1.6e-04 
      Ly = 4.0e-05 

// General parameters

      electronMass = 9.1094e-31 
      electronCharge = -1.6022e-19
      permit = 8.8542e-12 
      speedLight = 2.9979e8
      speedLight2 = speedLight*speedLight 
      electronCharge2 = electronCharge*electronCharge 
      qOverm = electronCharge/electronMass

// Plasma parameters

      elecPlasmaDensity =  6.9e23 

      elecPlasmaFreq = sqrt(electronCharge*qOverm*elecPlasmaDensity*electronCharge/permit)
 
//    Gaussian wave parameters

      waist = 2.0e-05 
      angFreq = 1.88e15 
      angFreq2 = angFreq*angFreq

      waveVector = sqrt((angFreq2 - elecPlasmaFreq*elecPlasmaFreq)/speedLight2)
      rayleighLength = waist*waist*waveVector/2.
      focus = Lx*0.5
}
Region
{
Grid
{
	J = Nx 
	x1s = 0.0
	x1f = Lx 
	n1 = 1.0
	K = Ny 
	x2s = 0.0
	x2f = Ly 
	n2 = 1.0
	Geometry = 1
	Rule
	{
	 Limit
	 n1 < 0.25
	 Fatal -- n1 < 0.25 grid spacing too nonuniform to ensure accuracy
	}
	Rule
	{
	 Algebra
	 J * K > 10000
	 Warning -- J*K >= 10000 may mean memory problems!
	}
}
Control
{
	dt = 1.6e-16 
	ElectrostaticFlag = 0
}
Species
{
        name = electrons
        m = electronMass 
        q = electronCharge 
        collisionModel = 1
}
Species
{
        name = ions
        m = 1.67e-27
        q = 1.6e-19
}
Load
{
        x1MinMKS = 0.000
        x1MaxMKS = Lx
        x2MinMKS = 0.00
        x2MaxMKS = Ly
        speciesName = electrons
        density = elecPlasmaDensity
        np2c = 2.0e11
        temperature = 0.0
}
Load
{
        x1MinMKS = 0.000
        x1MaxMKS = Lx
        x2MinMKS = 0.00
        x2MaxMKS = Ly
        speciesName = ions
        density = elecPlasmaDensity
        np2c = 0 
        temperature = 0.
}

PortGauss
{
	j1 = 0 
	k1 = 0 
	j2 = 0 
	k2 = Ny 
	normal = 1

	C = 1.0 
	a1 = 1.0 
	a0 = 0.0
	tdelay = 0
	trise = 1.0e-14 
	tpulse = 1.2e-13 
	tfall = 1.0e-14 

        gauss_w = waist 
        omegag = angFreq
        gauss_rl = rayleighLength 
        zf = focus 
        Agw = 6.0e12 
	
	EFFlag = 1 
	name = PortGauss
}
Conductor
{
	j1 = 0
	k1 = 0 
	j2 = Nx 
	k2 = 0 
	normal = 1
        reflection = 1
}
Conductor
{
	j1 = 0
	k1 = Ny 
	j2 = Nx 
	k2 = Ny 
	normal = -1
        reflection = 1
}
Conductor
{
        j1 = 0
        k1 = 0
        j2 = 0
        k2 = Ny
        normal = 1
        reflection = 1
}
Conductor
{
        j1 = Nx
        k1 = 0
        j2 = Nx
        k2 = Ny
        normal = -1
        reflection = 1
}

ExitPort
{
	j1 = Nx 
	k1 = 0 
	j2 = Nx 
	k2 = Ny 
	EFFlag = 1
	normal = -1
	name = ExitPort
	C = 0
	A = 0
}
}