Mirror 3 d #70
Open
Mirror 3 d #70
Conversation
This file contains hidden or bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
…ion from electric field term. Code is running but on much slower time scale than simulation done by Yuri.
…parameters in cgs units, and changed params to params_si, which has the parameters in SI units. Fixed typo of maxwellian function in mirror_parameters.m
… v_max to 10*params_si.a.vth to be consistent with Franz.
…lts for small energies. plot_mirror_output.m plots normalized conductivities. Instability still found for high energies of target distribution.
…riate condition for x, i.e. x>pi/2 and x<pi/2.
…ode appears to be relatively stable up to dt=1e-5. Also created PDE/mirror2_advection_Efield_div.m to test out advection terms from electric field acceleration.
…r2_velocity_div.m. Changed term names in mirror2_velocity_div.m to termE's and termC's. Fixed sign error in termC1 for collisional advection, and now case2 gives a deformed Maxwellian.
… off of fixing-fp2-stefan2 branch.
… rates to linear functions of velocity. Test passes
…is an additional option for the magnetic field. Modified mirror3.m to have a case where the magnetic field is set to this function.
…nd temperature in eV. Set up parameters for upcoming simulation with Juan, with exception of boundary conditions. Initial condition have been modified to be composed of maxwell and gaussian distributions.
…ed boundary condition for s-dimension.
…case 3 inside mirror3 to have initial condition of beam in physical space that can be modified.
…-dimensional moment and to plot temperatures over space inside mirror
…icient also has cosine in denominator because of math correction.
…ction coefficient
…ected with cos(z) in numerator.
…ave to clean it up
…s now the jacobian is an md_func, which means we need to pull it out of dimension, similar to the bc _funcs
…testing to see what effect advection term has in mirror3.m
…the space-time evolution of the solution. Changed the pitch angle distribution to a function of energy and temperature in mirror_parameters.m. Removed collisional terms from mirror3 to perform comparison with particle code.
…m included and sin(z) denominator eliminated from RHS in mirror1_pitch.m
…ode appears to be relatively stable up to dt=1e-5. Also created PDE/mirror2_advection_Efield_div.m to test out advection terms from electric field acceleration.
…ffusion term within the FP collision operator and the mass term.
…ients from mirror_FokkerPlanckCoeffs.m for the nonlinear solver. Still working on the solver. Also adjusting asgard_run_pde.m to update the params object correctly.
…ror_3D off fixing-fp2-complete-stefan2.
…pts into mirror_parameters.m for their params objects. asgard_run_pde.m has the FP coefficients commented out, and runaway_electron_plots.m has been refitted to run mirror2_velocity_div.m
…rd/DG-SparseGrid into fixing-fp2-complete-stefan2
…s of -1, +1 within DIV terms) of mirror2_velocity_div.m. This stabilizes the electric field terms. Modified runaway_electron_plots.m and spitzer_test.m to produce E field scans of mirror2_velocity_div.m.
…terms. All tests expect for mirror_* pass
…ide moment_integral_test.m. Merged fixing-fp2-complete-stefan2 with mirror_3D.
…rd/DG-SparseGrid into fixing-fp2-complete-stefan2
… in mirror_parameters.m from command line. Changed scaling of this value to 0.005*ratio2 in runaway_electron_plots.m which now produces data consistent with analytic values down to 0.05*E_D.
added 2 commits
December 14, 2021 14:33
…ts.m after rebase.
…allows mirror1_collision_div.m to pass test. Moved test parameters for mirror2_advection_Efield_div_test.m in case 4 inside mirror_parameters.m and now this test pasts. All tests in mirror_test.m pass after rebasing off reference.
Sign up for free
to join this conversation on GitHub.
Already have an account?
Sign in to comment
Add this suggestion to a batch that can be applied as a single commit.
This suggestion is invalid because no changes were made to the code.
Suggestions cannot be applied while the pull request is closed.
Suggestions cannot be applied while viewing a subset of changes.
Only one suggestion per line can be applied in a batch.
Add this suggestion to a batch that can be applied as a single commit.
Applying suggestions on deleted lines is not supported.
You must change the existing code in this line in order to create a valid suggestion.
Outdated suggestions cannot be applied.
This suggestion has been applied or marked resolved.
Suggestions cannot be applied from pending reviews.
Suggestions cannot be applied on multi-line comments.
Suggestions cannot be applied while the pull request is queued to merge.
Suggestion cannot be applied right now. Please check back later.
No description provided.