Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Pulsar debug theory sim #70

Draft
wants to merge 83 commits into
base: Pulsar_BC_CorEInRstar_plus3DR_BdipoleinRstar
Choose a base branch
from

Conversation

RevathiJambunathan
Copy link
Owner

No description provided.

RevathiJambunathan and others added 30 commits May 12, 2021 14:26
…tside, momentum of all particles is zero, fixed function that applied vacuum EB fields on particles
…ramp_omega_time input parameter. If ramp_omega_time < 0, PulsarParm::Omega(t) = omega_star.
…dd input parameter for the timescale for ramping up omega.
… + 1 so the damping factor at r=R_star is 1.0. Add pulsar.damping_scale as an input parameter defaulted to a reasonable value of 10.
RevathiJambunathan and others added 30 commits May 12, 2021 14:28
* Adding pulsar input and reading pulsar parameters using Don's PR

* added external E and B fields and the appropriate input parameters

* EOL whitespace

* seeting up bounds for pulsar surface

* pushing changes for comparison

* particles injected below surface, dipole magnetic field inside and outside, momentum of all particles is zero, fixed function that applied vacuum EB fields on particles

* add particles based on analytical surface charge density

* using local Er and Er_cor to compute particle introduction

* fixing particle injection with local Er that includes Er_external

* changes to pulsar

* michel and li efields inside and outside star. Particle injection rho_GJ

* adding input parameters to toggle rho gj scaling, ninjection, maxden

* removing unwanted insidePulsarBound function

* Damp electric fields internal to the pulsar when compiling with PULSAR on.

* Separate out cartesian ijk -> spherical radius into a utility function so it can be used elsewhere.

* Add PulsarParm::Omega(t) function that returns omega(t) depending on ramp_omega_time input parameter. If ramp_omega_time < 0, PulsarParm::Omega(t) = omega_star.

* Ran the style commands suggested by Travis CI to remove trailing whitespace.

* Upsating GetPosition() based on recent changes in master branch

* Damp electric fields internal to the pulsar when compiling with PULSAR on.

* Undo line removal.

* Calling External E and B for pulsar from GetExternal functions

* a few fixes in the getexternalEB and addplasma to ensure test works

* fixed Ecell

* removing particles if inside R_Star

* recent input file

* eol fix

* Add functor for pulsar particle injection

* fix eol

* Add Functor file

* fix eol

* add updated input file

* push external-EB-on-grid to debug. particle injection off.

* inject fixed number of particles when injection is on

* adding file with external EM field on grid

* fix indent

* fixing for loop to include last dimension

* add notebook with scaling

* fixing E-corotating insdie. E-quad outside. B-dipole everywhere

* gather external pulsar fields, add num_ppc particles, damp Maxwell fields inside pulsar

* include particles only outside r_star

* fix EOL

* particle injection inner ring and particle removal in same location

* damp EB fields inside R_Star, intro and remove part inside R_star

* fix typo in cell cell z coordinate

* add fraction of particles through pcounts of cell.

* Revert "add fraction of particles through pcounts of cell."

This reverts commit f3daa96.

* ensuring that fractional particle injection is uniform in cell

* we dont need the p.id=-1 to remove particles 1-fractional part injection

* updated input file

* terminate ifdef PULSAR

* remove unnecessary call to add external EB on EB Multifab. we call it directly in Fieldgather instead

* clean commented out lines in GetExternalFields

* remove commented out code and old unused PulsarEField and PulsarBField

* remove unused function init

* removing AddPulsarPlasma functor which is not used

* remove commented out code and calls to add plasma functor

* adding missing ;

* FillBoundary for EB after dampingEB

* using gpuarray for staggering used to damp EB fields

* default values in fieldgather doShapeN for problo, probhi, and cur_time

* empty array

* pass by reference

* take user-defined cc to compute r

* Add Dons recommendation to change wt for fractional injection

* modifying particle wt, redefining insidePulsarBoundsCC, and removing unnecessary functions

* delete redundant function definition

* remove redundant function definitions

* remove unused variable

* remove unnecesary empty line

* adding user-inpu for modifying particle weight

* querying input for modiyfing particle weight

* remove header file decleration in Make package

* remove unnecesary header

* remove unused function call

* use input value of star center to allow for shifted centers in the input

* adding a separate file for gather pulsar external EB using grid resolution

* removing Pulsar related calls from the core gather functions

* call pulsar field gather separate from field gather

* deposit charge is modified and therefore not const

* reorder function definitions and add missingsemicolon

* fix typo of variable in pulsarparameters

* add the empty line back

* undo any changes in FieldGather so it is the same as WarpX development

* remove modifications made to fieldgather

* remove unnecessary blank line in at end of file

* remove tab

* EOL

* input parameters for damping EB and external EB to distinguish inside and outside

* adding a few more input parameters

* adding particle absortion radius as input

* remove unnecessary print statement

* add min max bounds for defining particle injection region

* EOL

* fix typo

* fix comments from the review

Co-authored-by: Donald E. Willcox <[email protected]>
Co-authored-by: Revathi Jambunathan <[email protected]>
… plasma sim. new input flags to turn on theory EB on particles/grid
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
None yet
Projects
None yet
Development

Successfully merging this pull request may close these issues.

2 participants