VectorDepositor

VectorDepositor

Flexible way of depositing current from charged particles in a simulation into occurrences of depField, instead of the classical way of depositing current into the last three components of a SumRhoJ field. With a VectorDepositor you can simply employ a VectorDepositor and a three-vector J in an EM PIC simulation instead of using the four-vector SumRhoJ, as has been classically used, thus saving on memory.

Note

The VectorDepositor should not be used in MultiField block of kind = emMultiField, or errors will occur in the simulation.

Note

Higher order particles are incompatible with the polar/cylindrical axis. If particles in the simulation pass near the polar/cylindrical coordinate axis at r=0, only areaWeightingCP or esirk1stOrderCP can be used.

VectorDepositor Parameters

kind (string, default = areaWeighting)

Type of deposition algorithm; choices are:

• areaWeighting

• esirk1stOrder

• esirk2ndOrder

• esirk7thOrder

• areaWeightingCP (for coordProdGrid)

• esirk1stOrderCP (for coordProdGrid) When working with CoordProd grids these work slighty different for cartesian and cylindrical grids. If working with cylindrical grids the deposited quantity is Charge, and not ChargeDensity. To compensate for this it is necessary to use a scaling field that will apply the function 1/(PI*DZ*( (y+DR/2)^2 - (max(0, y-DR/2) )^2)) For an example of this it is recommended to review a visual setup example input file using the “maximum stabilty” particle deposition.

• esirk2ndOrderCP (for coordProdGrid)

• esirk7thOrderCP (for coordProdGrid)

depField (string)

Vector depfield defined in the MultiField block that will contain the deposited current.

Note

When using higher order deposition (esirk2ndOrder, …, esirk7thOrder), maxIntDepHalfWidth must be set accordingly in the Grid block (see Additional Attributes for Particle Simulations).

Example VectorDepositor Block

<VectorDepositor~currDep>
    kind = areaWeighting
    depField = myEmField.J
</VectorDepositor>

Example Reference to VectorDepositor in a Species Block

<Species  electrons>
    kind = relBorisDF
    currDeps = [ currDep ]
    .
    .
    .
</Species>