colpul.inp

pic 模拟程序！面向对象

colpul.inp
{

File:     colpul.inp

Purpose:  Colliding pulses simulation.


Basic parameters:
    Plasma of density 6.9e23, into which is injected a laser
    pulse that is 132.7 fs half-sine pulses.
    Cartesian geometry.

Plasma variables:
  elecPlasmaDensity = 6.9e23 

Laser variables:
    Half sine pulse
  pulShp_p0 = 2
    Time delay to start pulse
  tdelay_p0 = 0.

====================================================================

 - The wave is defined by the variables:

     - spotSize_p0      = Half width at amplitude of central amplitde over
                          e at the beam waist. 
     - waveLeng_p0      = wavelenght.
     - focus_p0         = position of focus. 
     - amp_p0           = peak wave amplitude.    
     - chirp_p0         = change of wavelength from begining of pulse
                          to end.
     - pulShp_p0        = pulse shape (0=trapezoid; 1=Gaussian; 2=half-sine)
     - tdelay_p0        = time delay to start pulse. 
     - pulLeng_p0       = pulse length (meters). (half-width for trapezoid and
                          Gaussian full for half-sine).
}

Variables
{
// Grid variables

      Nx = 4000 
      Ny = 150 
      Lx = 160.0e-06 
      Ly = 80.0e-06 
      DX = Lx/Nx

// 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
      Pi = 3.14159265358979323846 

// Plasma parameters

      elecPlasmaDensity =  6.9e23 
      elecPlasmaFreq = sqrt(electronCharge*qOverm*elecPlasmaDensity/permit)
 
//=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-


}
Region
{
Grid
{
	J = Nx 
	x1s = 0.0
	x1f = Lx 
	n1 = 1.0
	K = Ny 
	x2s = 0.0
	x2f = Ly 
	n2 = 1.0
	Geometry = 1
}
Control
{
	dt = 8.4e-17 
	ElectrostaticFlag = 0
	movingWindow = 1 
	shiftDelayTime = 5.17e-13 
}
Species
{
        name = electrons
        m = electronMass 
        q = electronCharge 
}
Species
{
        name = ions
        m = 1.67e-27
        q = 1.6e-19
}
Load
{
        x1MinMKS = 0.0
        x1MaxMKS = 20.0e-06
        x2MinMKS = 0.00
        x2MaxMKS = Ly
        speciesName = electrons
        density = 0.0
        np2c = 0.0
        LoadMethodFlag = 1
}
Load
{
        x1MinMKS = 0.0
        x1MaxMKS = 20.0e-06
        x2MinMKS = 0.00
        x2MaxMKS = Ly
        speciesName = ions
        density = 0.0
        np2c = 0.0
}
Load
{
        x1MinMKS = 20.0e-06
        x1MaxMKS = 100.0e-06
        x2MinMKS = 0.00
        x2MaxMKS = Ly
        speciesName = electrons
        density = elecPlasmaDensity
        analyticF = 0.5*elecPlasmaDensity*(1.0 - cos(PI*(x1-20.0e-06)/80.0e-06))
        np2c = 3.68e09
        LoadMethodFlag = 1
}
Load
{
        x1MinMKS = 20.0e-06
        x1MaxMKS = 100.0e-06
        x2MinMKS = 0.00
        x2MaxMKS = Ly
        speciesName = ions
        density = elecPlasmaDensity
        analyticF = 0.5*elecPlasmaDensity*(1.0 - cos(PI*(x1-20.0e-06)/80.0e-06))
        np2c = 0.0
        LoadMethodFlag = 1
}
Load
{
        x1MinMKS = 100.0e-06
        x1MaxMKS = Lx - DX
        x2MinMKS = 0.00
        x2MaxMKS = Ly
        speciesName = electrons
        density = elecPlasmaDensity
        np2c = 3.68e09
        LoadMethodFlag = 1
}
Load
{
        x1MinMKS = 100.0e-06
        x1MaxMKS = Lx - DX
        x2MinMKS = 0.00
        x2MaxMKS = Ly
        speciesName = ions
        density = elecPlasmaDensity
        np2c = 0
        LoadMethodFlag = 1
}
Load
{
        Name = shiftLoad
        x1MinMKS = Lx - DX
        x1MaxMKS = Lx
        x2MinMKS = 0.00
        x2MaxMKS = Ly
        speciesName = electrons
        density = elecPlasmaDensity
        np2c = 3.68e09
        LoadMethodFlag = 1
}
Load
{
        Name = shiftLoad
        x1MinMKS = Lx - DX
        x1MaxMKS = Lx
        x2MinMKS = 0.00
        x2MaxMKS = Ly
        speciesName = ions
        density = elecPlasmaDensity
        np2c = 0
        LoadMethodFlag = 1
}

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

// Wave (0)
	A = 0
	C = 1.0 
	pulShp_p0 = 2 
	tdelay_p0 = 0.0 
	pulLeng_p0 = 40.0e-06 
        chirp_p0 = 0
        spotSize_p0 = 1.5e-05  
        waveLeng_p0 = 0.8e-06 
        focus_p0 = 3.0e-05 
        amp_p0 = 5.0e12 

// Wave (1)
	A = 0
	C = 1.0 
        pulShp_p1 = 2 
        tdelay_p1 = 2.1681e-13
        pulLeng_p1 = 20.0e-06 
        chirp_p1 = 0
        spotSize_p1 = 15.0e-06
        waveLeng_p1 = 0.83e-06
        focus_p1 = 2.0e-05 
        amp_p1 = 3.6e12 

	EFFlag = 0 
	name = PortGauss
}
PortGauss
{
        j1 = Nx 
        k1 = 0
        j2 = Nx 
        k2 = Ny
        normal = -1

// Wave (0)
        A = 0
        C = 1.0
        pulShp_p0 = 2
        tdelay_p0 = 10.0
        pulLeng_p0 = 40.0e-06
        chirp_p0 = 0
        spotSize_p0 = 1.5e-05
        waveLeng_p0 = 0.8e-06
        focus_p0 = 3.0e-05 
        amp_p0 = 0.0 

// Wave (1)
        A = 0
        C = 1.0
        pulShp_p1 = 2
        tdelay_p1 = 4.336e-13
        pulLeng_p1 = 20.0e-06
        chirp_p1 = 0
        spotSize_p1 = 15.0e-06
        waveLeng_p1 = 0.8e-06
        focus_p1 = 2.0e-05
        amp_p1 = 3.84e12 

        EFFlag = 0
        name = PortGauss
}

Conductor
{
	j1 = 0
	k1 = Ny 
	j2 = Nx 
	k2 = Ny 
	normal = -1
}
Conductor
{
	j1 = 0
	k1 = 0 
	j2 = Nx 
	k2 = 0 
	normal = 1
}
}