PEC Namespace Reference

Functions
AMREX_GPU_DEVICE AMREX_FORCE_INLINE bool	is_boundary_PEC (amrex::GpuArray< int, 3 > const &fboundary, int dir)
	Determines if the field boundary condition stored in fboundary in direction, dir, is PEC. More...
AMREX_GPU_DEVICE AMREX_FORCE_INLINE bool	is_boundary_reflecting (amrex::GpuArray< ParticleBoundaryType, 3 > const &pboundary, int dir)
	Determines if the particle boundary condition stored in pboundary in direction, dir, is reflecting. More...
AMREX_GPU_DEVICE AMREX_FORCE_INLINE int	get_cell_count_to_boundary (const amrex::IntVect &dom_lo, const amrex::IntVect &dom_hi, const amrex::IntVect &ijk_vec, const amrex::IntVect &is_nodal, const int idim, const int iside)
	Calculates the number of grid points the given index is pass the domain boundary i.e. a value of +1 means the current cell is outside of the simulation domain by 1 cell. Note that the high side domain boundary is between cell dom_hi and dom_hi+1 for cell centered grids and on cell dom_hi+1 for nodal grid. This is why (dom_hi[idim] + is_nodal[idim]) is used below. More...
AMREX_GPU_DEVICE AMREX_FORCE_INLINE void	SetEfieldOnPEC (const int icomp, const amrex::IntVect &dom_lo, const amrex::IntVect &dom_hi, const amrex::IntVect &ijk_vec, const int n, amrex::Array4< amrex::Real > const &Efield, const amrex::IntVect &is_nodal, amrex::GpuArray< int, 3 > const &fbndry_lo, amrex::GpuArray< int, 3 > const &fbndry_hi)
	Sets the electric field value tangential to the PEC boundary to zero. The tangential Efield components in the guard cells outside the domain boundary are set equal and opposite to the field in the valid cells at their mirrored locations. The normal Efield components in the guard cells are set equal to the field in the valid cells at their mirrored locations. The number or depth of guard cells updated is equal to the shape factor of particles in each dimension. For corner cells with mixed boundaries, the mirror location could be outside valid region, while still ensuring PEC condition is maintained across the PEC boundary, and the necessary sign change is accounted for depending on if the component, icomp, is tangential or normal to the PEC boundary. More...
AMREX_GPU_DEVICE AMREX_FORCE_INLINE void	SetBfieldOnPEC (const int icomp, const amrex::IntVect &dom_lo, const amrex::IntVect &dom_hi, const amrex::IntVect &ijk_vec, const int n, amrex::Array4< amrex::Real > const &Bfield, const amrex::IntVect &is_nodal, amrex::GpuArray< int, 3 > const &fbndry_lo, amrex::GpuArray< int, 3 > const &fbndry_hi)
	Sets the magnetic field value normal to the PEC boundary to zero. The tangential (and normal) field value of the guard cells outside the domain boundary are set equal (and opposite) to the respective field components in the valid cells at their mirrored locations. The number or depth of guard cells updated is equal to the shape factor of particles in each dimension. More...
AMREX_GPU_DEVICE AMREX_FORCE_INLINE void	SetRhoOrJfieldFromPEC (const int n, const amrex::IntVect &ijk_vec, amrex::Array4< amrex::Real > const &field, amrex::GpuArray< GpuArray< int, 2 >, AMREX_SPACEDIM > const &mirrorfac, amrex::GpuArray< GpuArray< amrex::Real, 2 >, AMREX_SPACEDIM > const &psign, amrex::GpuArray< GpuArray< bool, 2 >, AMREX_SPACEDIM > const &is_pec, amrex::GpuArray< bool, AMREX_SPACEDIM > const &tangent_to_bndy, amrex::Box const &fabbox)
	Sets the rho or J field value in cells close to and on a PEC boundary. The charge/current density deposited in the guard cells are either reflected back into the simulation domain (if a reflecting particle boundary is used), or the opposite charge/current density is deposited back in the domain to capture the effect of an image charge. The charge/current density on the PEC boundary is set to 0 while values in the guard cells are set equal (and opposite) to their mirror location inside the domain - representing image charges - in the normal (tangential) direction. More...
AMREX_GPU_DEVICE AMREX_FORCE_INLINE void	SetNeumannOnPEC (const int n, const amrex::IntVect &ijk_vec, amrex::Array4< amrex::Real > const &field, amrex::GpuArray< GpuArray< int, 2 >, AMREX_SPACEDIM > const &mirrorfac, amrex::GpuArray< GpuArray< bool, 2 >, AMREX_SPACEDIM > const &is_pec, amrex::Box const &fabbox)
	This function sets the given field value on a PEC boundary to enforce a Neumann boundary condition (zero derivative) in the normal direction. More...
bool	isAnyBoundaryPEC ()
void	ApplyPECtoEfield (std::array< amrex::MultiFab *, 3 > Efield, int lev, PatchType patch_type, bool split_pml_field=false)
	Sets the tangential electric field at the PEC boundary to zero. The guard cell values are set equal and opposite to the valid cell field value at the respective mirror locations. More...
void	ApplyPECtoBfield (std::array< amrex::MultiFab *, 3 > Bfield, int lev, PatchType patch_type)
	Sets the normal component of the magnetic field at the PEC boundary to zero. The guard cell values are set equal and opposite to the valid cell field value at the respective mirror locations. More...
void	ApplyPECtoRhofield (amrex::MultiFab *rho, int lev, PatchType patch_type)
	Reflects charge density deposited over the PEC boundary back into the simulation domain. More...
void	ApplyPECtoJfield (amrex::MultiFab *Jx, amrex::MultiFab *Jy, amrex::MultiFab *Jz, int lev, PatchType patch_type)
	Reflects current density deposited over the PEC boundary back into the simulation domain. More...
void	ApplyPECtoElectronPressure (amrex::MultiFab *Pefield, int lev, PatchType patch_type)
	Apply the PEC boundary to the electron pressure field. More...

Function Documentation

ApplyPECtoBfield()

void PEC::ApplyPECtoBfield	(	std::array< amrex::MultiFab *, 3 >	Bfield,
		int	lev,
		PatchType	patch_type
	)

Sets the normal component of the magnetic field at the PEC boundary to zero. The guard cell values are set equal and opposite to the valid cell field value at the respective mirror locations.

Parameters

[in,out]	Bfield	Boundary values of normal Bfield are set to zero.
[in]	lev	level of the Multifab
[in]	patch_type	coarse or fine

ApplyPECtoEfield()

void PEC::ApplyPECtoEfield	(	std::array< amrex::MultiFab *, 3 >	Efield,
		int	lev,
		PatchType	patch_type,
		bool	split_pml_field = false
	)

Sets the tangential electric field at the PEC boundary to zero. The guard cell values are set equal and opposite to the valid cell field value at the respective mirror locations.

Parameters

[in,out]	Efield	Boundary values of tangential Efield are set to zero
[in]	lev	level of the Multifab
[in]	patch_type	coarse or fine
[in]	split_pml_field	whether pml the multifab is the regular Efield or split pml field

ApplyPECtoElectronPressure()

void PEC::ApplyPECtoElectronPressure	(	amrex::MultiFab *	Pefield,
		int	lev,
		PatchType	patch_type
	)

Apply the PEC boundary to the electron pressure field.

Parameters

[in,out]	Pefield	Multifab containing the electron pressure
[in]	lev	level of the Multifab
[in]	patch_type	coarse or fine

ApplyPECtoJfield()

void PEC::ApplyPECtoJfield	(	amrex::MultiFab *	Jx,
		amrex::MultiFab *	Jy,
		amrex::MultiFab *	Jz,
		int	lev,
		PatchType	patch_type
	)

Reflects current density deposited over the PEC boundary back into the simulation domain.

Parameters

[in,out]	Jx,Jy,Jz	Multifabs containing the current density
[in]	lev	level of the Multifab
[in]	patch_type	coarse or fine

ApplyPECtoRhofield()

void PEC::ApplyPECtoRhofield	(	amrex::MultiFab *	rho,
		int	lev,
		PatchType	patch_type
	)

Reflects charge density deposited over the PEC boundary back into the simulation domain.

Sets the rho field value in cells close to and inside a PEC boundary. The charge density deposited in the guard cells are either reflected back into the simulation domain (if a reflecting particle boundary is used), or the opposite charge density is deposited back in the domain to capture the effect of an image charge. The charge density on the PEC boundary is set to 0 while values in the guard cells are set equal and opposite to their mirror location inside the domain - representing image charges.

Parameters

[in,out]	rho	Multifab containing the charge density
[in]	lev	level of the Multifab
[in]	patch_type	coarse or fine

get_cell_count_to_boundary()

AMREX_GPU_DEVICE AMREX_FORCE_INLINE int PEC::get_cell_count_to_boundary	(	const amrex::IntVect &	dom_lo,
		const amrex::IntVect &	dom_hi,
		const amrex::IntVect &	ijk_vec,
		const amrex::IntVect &	is_nodal,
		const int	idim,
		const int	iside
	)

Calculates the number of grid points the given index is pass the domain boundary i.e. a value of +1 means the current cell is outside of the simulation domain by 1 cell. Note that the high side domain boundary is between cell dom_hi and dom_hi+1 for cell centered grids and on cell dom_hi+1 for nodal grid. This is why (dom_hi[idim] + is_nodal[idim]) is used below.

Parameters

[in]	dom_lo,dom_hi	Domain boundaries
[in]	ijk_vec	Cell coordinates
[in]	is_nodal	Whether the field of interest is nodal
[in]	idim	Dimension of interest
[in]	iside	0 for low and 1 for high

Returns

number of grid points to the boundary

is_boundary_PEC()

AMREX_GPU_DEVICE AMREX_FORCE_INLINE bool PEC::is_boundary_PEC	(	amrex::GpuArray< int, 3 > const &	fboundary,
		int	dir
	)

Determines if the field boundary condition stored in fboundary in direction, dir, is PEC.

Parameters

[in]	fboundary	Value containing boundary type
[in]	dir	direction

Returns

1 if the boundary type is PEC else 0

is_boundary_reflecting()

AMREX_GPU_DEVICE AMREX_FORCE_INLINE bool PEC::is_boundary_reflecting	(	amrex::GpuArray< ParticleBoundaryType, 3 > const &	pboundary,
		int	dir
	)

Determines if the particle boundary condition stored in pboundary in direction, dir, is reflecting.

Parameters

[in]	pboundary	Value containing boundary type
[in]	dir	direction

Returns

1 if the boundary type is reflecting else 0

isAnyBoundaryPEC()

bool PEC::isAnyBoundaryPEC

(

)

Returns 1 if any domain boundary is set to PEC, else returns 0.

SetBfieldOnPEC()

AMREX_GPU_DEVICE AMREX_FORCE_INLINE void PEC::SetBfieldOnPEC	(	const int	icomp,
		const amrex::IntVect &	dom_lo,
		const amrex::IntVect &	dom_hi,
		const amrex::IntVect &	ijk_vec,
		const int	n,
		amrex::Array4< amrex::Real > const &	Bfield,
		const amrex::IntVect &	is_nodal,
		amrex::GpuArray< int, 3 > const &	fbndry_lo,
		amrex::GpuArray< int, 3 > const &	fbndry_hi
	)

Sets the magnetic field value normal to the PEC boundary to zero. The tangential (and normal) field value of the guard cells outside the domain boundary are set equal (and opposite) to the respective field components in the valid cells at their mirrored locations. The number or depth of guard cells updated is equal to the shape factor of particles in each dimension.

For 3D : x component is tangential to the y-boundary and z-boundary y component is tangential to the x-boundary and z-boundary z component is tangential to the x-boundary and y-boundary x component is normal to the x-boundary y component is normal to the y-boundary z component is normal to the z-boundary where, x-boundary is the yz-plane at x=xmin and x=xmax y-boundary is the xz-plane at y=ymin and y=ymax z-boundary is the xy-plane at z=zmin and z=zmax

For 2D : WarpX uses X-Z as the two dimensions x component is tangential to the z-boundary y component is tangential to the x-boundary and z-boundary z component is tangential to the x-boundary x component is normal to the x-boundary y component is not normal to any boundary (Only xz dimensions in 2D) z component is normal to the z-boundary where, x-boundary is along the line z at x=xmin and x=xmax z-boundary is along the line x at z=zmin and z=zmax

For 1D : WarpX uses Z as the only dimension x component is tangential to the z-boundary y component is tangential to the z-boundary z component is not tangential to the z-boundary x component is not normal to any boundary (Only z dimension in 1D) y component is not normal to any boundary (Only z dimension in 1D) z component is normal to the z-boundary where, z-boundary is a point at z=zmin and z=zmax

For RZ : WarpX uses R-Z as the two dimensions r component is tangential to the z-boundary theta_component is tangential to the r-boundary and z-boundary z component is tangential to the r-boundary r component is normal to the r-boundary theta_component is not normal to any boundary (on RZ dimensions are modeled) z component is normal to the z-boundary where, r-boundary is along the line z at r=rmin and r=rmax z-boundary is along the line r at z=zmin and z=zmax

Parameters

[in]	icomp	component of the Bfield being updated (0=x, 1=y, 2=z in Cartesian) (0=r, 1=theta, 2=z in RZ)
[in]	dom_lo	index value of the lower domain boundary (cell-centered)
[in]	dom_hi	index value of the higher domain boundary (cell-centered)
[in]	ijk_vec	indices along the x(i), y(j), z(k) of Efield Array4
[in]	n	index of the MultiFab component being updated
[in]	Bfield	field data to be updated if (ijk) is at the boundary or a guard cell
[in]	is_nodal	staggering of the field data being updated.
[in]	fbndry_lo	Field boundary type at the lower boundaries
[in]	fbndry_hi	Field boundary type at the upper boundaries

SetEfieldOnPEC()

AMREX_GPU_DEVICE AMREX_FORCE_INLINE void PEC::SetEfieldOnPEC	(	const int	icomp,
		const amrex::IntVect &	dom_lo,
		const amrex::IntVect &	dom_hi,
		const amrex::IntVect &	ijk_vec,
		const int	n,
		amrex::Array4< amrex::Real > const &	Efield,
		const amrex::IntVect &	is_nodal,
		amrex::GpuArray< int, 3 > const &	fbndry_lo,
		amrex::GpuArray< int, 3 > const &	fbndry_hi
	)

Sets the electric field value tangential to the PEC boundary to zero. The tangential Efield components in the guard cells outside the domain boundary are set equal and opposite to the field in the valid cells at their mirrored locations. The normal Efield components in the guard cells are set equal to the field in the valid cells at their mirrored locations. The number or depth of guard cells updated is equal to the shape factor of particles in each dimension. For corner cells with mixed boundaries, the mirror location could be outside valid region, while still ensuring PEC condition is maintained across the PEC boundary, and the necessary sign change is accounted for depending on if the component, icomp, is tangential or normal to the PEC boundary.

For 3D : x component is tangential to the y-boundary and z-boundary y component is tangential to the x-boundary and z-boundary z component is tangential to the x-boundary and y-boundary x component is normal to the x-boundary y component is normal to the y-boundary z component is normal to the z-boundary where, x-boundary is the yz-plane at x=xmin and x=xmax y-boundary is the xz-plane at y=ymin and y=ymax z-boundary is the xy-plane at z=zmin and z=zmax

For 2D : WarpX uses X-Z as the two dimensions x component is tangential to the z-boundary y component is tangential to the x-boundary and z-boundary z component is tangential to the x-boundary x component is normal to the x-boundary y component is not normal to any boundary (Only xz dimensions in 2D) z component is normal to the z-boundary where, x-boundary is along the line z at x=xmin and x=xmax z-boundary is along the line x at z=zmin and z=zmax

For 1D : WarpX uses Z as the only dimension x component is tangential to the z-boundary y component is tangential to the z-boundary z component is not tangential to the z-boundary x component is not normal to any boundary (Only z dimension in 1D) y component is not normal to any boundary (Only z dimension in 1D) z component is normal to the z-boundary where, z-boundary is a point at z=zmin and z=zmax

For RZ : WarpX uses R-Z as the two dimensions r component is tangential to the z-boundary theta_component is tangential to the r-boundary and z-boundary z component is tangential to the r-boundary r component is normal to the r-boundary theta_component is not normal to any boundary (on RZ dimensions are modeled) z component is normal to the z-boundary where, r-boundary is along the line z at r=rmin and r=rmax z-boundary is along the line r at z=zmin and z=zmax

Parameters

[in]	icomp	component of the Efield being updated (0=x, 1=y, 2=z in Cartesian) (0=r, 1=theta, 2=z in RZ)
[in]	dom_lo	index value of the lower domain boundary (cell-centered)
[in]	dom_hi	index value of the higher domain boundary (cell-centered)
[in]	ijk_vec	indices along the x(i), y(j), z(k) of Efield Array4
[in]	n	index of the MultiFab component being updated
[in]	Efield	field data to be updated if (ijk) is at the boundary or a guard cell
[in]	is_nodal	staggering of the field data being updated.
[in]	fbndry_lo	Field boundary type at the lower boundaries
[in]	fbndry_hi	Field boundary type at the upper boundaries

SetNeumannOnPEC()

AMREX_GPU_DEVICE AMREX_FORCE_INLINE void PEC::SetNeumannOnPEC	(	const int	n,
		const amrex::IntVect &	ijk_vec,
		amrex::Array4< amrex::Real > const &	field,
		amrex::GpuArray< GpuArray< int, 2 >, AMREX_SPACEDIM > const &	mirrorfac,
		amrex::GpuArray< GpuArray< bool, 2 >, AMREX_SPACEDIM > const &	is_pec,
		amrex::Box const &	fabbox
	)

This function sets the given field value on a PEC boundary to enforce a Neumann boundary condition (zero derivative) in the normal direction.

Parameters

[in]	n	index of the MultiFab component being updated
[in]	ijk_vec	indices along the x(i), y(j), z(k) of the rho Array4
[in,out]	field	field data to be updated
[in]	mirrorfac	mirror cell is given by mirrorfac - ijk_vec
[in]	psign	Whether the field value should be flipped across the boundary
[in]	is_pec	Whether the given boundary is PEC
[in]	tangent_to_bndy	Whether a given direction is perpendicular to the boundary
[in]	fabbox	multifab box including ghost cells

SetRhoOrJfieldFromPEC()

AMREX_GPU_DEVICE AMREX_FORCE_INLINE void PEC::SetRhoOrJfieldFromPEC	(	const int	n,
		const amrex::IntVect &	ijk_vec,
		amrex::Array4< amrex::Real > const &	field,
		amrex::GpuArray< GpuArray< int, 2 >, AMREX_SPACEDIM > const &	mirrorfac,
		amrex::GpuArray< GpuArray< amrex::Real, 2 >, AMREX_SPACEDIM > const &	psign,
		amrex::GpuArray< GpuArray< bool, 2 >, AMREX_SPACEDIM > const &	is_pec,
		amrex::GpuArray< bool, AMREX_SPACEDIM > const &	tangent_to_bndy,
		amrex::Box const &	fabbox
	)

Sets the rho or J field value in cells close to and on a PEC boundary. The charge/current density deposited in the guard cells are either reflected back into the simulation domain (if a reflecting particle boundary is used), or the opposite charge/current density is deposited back in the domain to capture the effect of an image charge. The charge/current density on the PEC boundary is set to 0 while values in the guard cells are set equal (and opposite) to their mirror location inside the domain - representing image charges - in the normal (tangential) direction.

Parameters

[in]	n	index of the MultiFab component being updated
[in]	ijk_vec	indices along the x(i), y(j), z(k) of the rho Array4
[in,out]	field	field data to be updated
[in]	mirrorfac	mirror cell is given by mirrorfac - ijk_vec
[in]	psign	Whether the field value should be flipped across the boundary
[in]	is_pec	Whether the given boundary is PEC
[in]	tangent_to_bndy	Whether a given direction is perpendicular to the boundary
[in]	fabbox	multifab box including ghost cells