port.cpp

pic 模拟程序！面向对象

/*
====================================================================

PORT.CPP

0.99	(NTG 12-29-93) Separated into individual module from pic.h.
0.991	(JohnV 01-03-93) Aesthetics, compile.
0.992	(JohnV 03-23-94) Streamline includes, remove pic.h.

====================================================================
*/

#include "port.h"
#include "fields.h"
#include "ptclgrp.h"
#include "misc.h"

#ifdef _MSC_VER
using std::cout;
#endif
//#include "cylwvgd.h"

Scalar AnalyticEz(Scalar z, Scalar r, Scalar f, Scalar t, Scalar a)
{Scalar omega = 2.0*3.14159*f, c = SPEED_OF_LIGHT, gamma = 2.4048;
// Scalar k = sqrt((omega*omega/(c*c)) - gamma*gamma/(a*a));
//Scalar Ez = bessj0(gamma*r/a)*cos(k*z - omega*t);
 Scalar Ez = bessj0(gamma*r/a);
 return Ez;
}

Scalar AnalyticEr(Scalar z, Scalar r, Scalar f, Scalar t, Scalar a)
{Scalar omega = 2.0*3.14159*f, c = SPEED_OF_LIGHT, gamma = 3.832;
// Scalar k = sqrt((omega*omega/(c*c)) - gamma*gamma/(a*a));
// Scalar Er = (k*a/gamma)*bessj1(gamma*r/a)*sin(k*z - omega*t);
 Scalar Ez = bessj0(gamma*r/a);
 return Ez;
}

Scalar AnalyticBphi(Scalar z, Scalar r, Scalar f, Scalar t, Scalar a)
{Scalar omega = 2.0*3.14159*f, c = SPEED_OF_LIGHT, gamma = 2.4048;
 Scalar k = sqrt((omega*omega/(c*c)) - gamma*gamma/(a*a));
 Scalar coeff = sqrt(gamma*gamma + k*k*a*a)/gamma;
// Scalar Bphi = coeff*bessj1(gamma*r/a)*sin(k*z - omega*t);
 Scalar Bphi = coeff*bessj1(gamma*r/a);
 return Bphi;
}

Scalar AnalyticEphi(Scalar z, Scalar r, Scalar f, Scalar t, Scalar a)
{Scalar omega = 2.0*3.14159*f, c = SPEED_OF_LIGHT, gamma = 3.832;
 Scalar k = sqrt((omega*omega/(c*c)) - gamma*gamma/(a*a));
 Scalar Ephi = bessj1(gamma*r/a)*sin(k*z - omega*t);
 return Ephi;
}

Port::Port(oopicList <LineSegment> *segments) : Boundary(segments)
{

//	no code yet
	if(segments->nItems() > 1) {
	  cout << "Warning, Port-derived boundaries can only have 1 segment.\n";
	  cout << "Check your input file.\n";
	}
	BCType=DIELECTRIC_BOUNDARY;

}

//--------------------------------------------------------------------
//	Apply boundary conditions to fields locally.

void Port::applyFields(Scalar t,Scalar dt)
{
   if(dt==0) return;
	Grid*	grid = fields->get_grid();
	Scalar	freq = 10E9;
	Scalar	zforEz = (grid->getMKS(j1 + 1, k1).e1()
		+ grid->getMKS(j1, k1).e1())/2.0;
	Scalar	zforEr = grid->getMKS(j1, k1).e1();
	Scalar   a = grid->getMKS(0, grid->getK()).e2();
	for (int k = k1; k <= k2; k++)
	{
		Scalar	rforEz = grid->getMKS(j1, k).e2();
   	Scalar	intEz = grid->dl1(j1, k)
			*AnalyticEz(zforEz, rforEz, freq, t, a);
   	fields->setIntEdl(j1, k, 1, intEz);
   	if (k < k2)
    	{
			Scalar	rforEr = (grid->getMKS(j1, k + 1).e2()
				+ grid->getMKS(j1, k).e2())/2.0;
	   	Scalar	intEr = grid->dl2(j1, k)*
				AnalyticEr(zforEr, rforEr, freq, t, a);
   		fields->setIntEdl(j1, k, 2, intEr);
		}
	}
}

//--------------------------------------------------------------------
//	Set the passive bc for fields at the port.  Currently just a
//	copy of those for conductor.

void Port::setPassives()
{
	if (alongx2())								//	vertical
	{
		for (int k=MIN(k1,k2); k<MAX(k1,k2); k++)
		{
			fields->set_iC2(j1, k, 0);       	
			fields->set_iC3(j1, k, 0);
		}
		fields->set_iC3(j1, k2, 0);
	}
	else											//	horizontal
	{
		for (int j=MIN(j1,j2); j<MAX(j1,j2); j++)
      {
			fields->set_iC1(j, k1, 0);
			fields->set_iC3(j, k1, 0);
		}
		fields->set_iC3(j2, k1, 0);
	}
}

//--------------------------------------------------------------------
//	Port::emit() simply deletes Particles in its stack.  May add some
//	diagnostics for particles collected in the future.
#ifndef __linux__
#ifndef _MSC_VER
#pragma argsused
#endif
#endif

ParticleList& Port::emit(Scalar t,Scalar dt, Species* species)
{
	while(!particleList.isEmpty())
	{
		Particle* p = particleList.pop();
		delete p;
	}
   return particleList;
}


void Port::collect(Particle& p, Vector3& dxMKS)
{
	Boundary::collect(p, dxMKS); // do statistics
	delete &p;
}