Skip to content

Commit

Permalink
Long Domain 1um Resolution (AMReX-Microelectronics#32)
Browse files Browse the repository at this point in the history
* change plot.period to plot.intervals

* changed plot.period to plot.intervals

* updated Examples/Tests/Macroscopic_Maxwell/inputs_3d_LLG_noMs

* some tweaks Examples/Tests/Macroscopic_Maxwell/inputs_3d_LLG_noMs

* fixed bug in getting mag_exchange_init_style and mag_anisotropy_init_style

* EOL whitespace

* saving Ricky's test

* more tests

* B_nodal_flag was declared inside of a #ifdef WARPX_MAG_LLG, preventing compilation

Removed references to Bfield_fp_old as it is not used

* remove reference to Bfield_fp_old

* remove geometry.is_periodic and particle species stuff

* fix typo

* remove print statement

* more inputs file cleanup

* more description

* 4 walls

* do_pml not needed

* remove do_pml flag

* verbose=1

* cleanup

* infinite walls in y

* circuit updates

* more inputs

* tweak to prevent leakage

* EOL whitespace

* renname file

* cleanup

* renname stuff

* more comments

* cleanup

* remove bottom PEC wall

* Apply suggestions from code review

Co-authored-by: Revathi  Jambunathan <[email protected]>

* Apply suggestions from code review

Co-authored-by: Revathi  Jambunathan <[email protected]>

* Apply suggestions from code review

* long domain with 1um resolution ant ~5GHz current source

* Update Examples/Tests/circuits/inputs_longdomain

* Apply suggestions from code review

* fix unused inputs

* update README

Co-authored-by: rzamora2 <[email protected]>
Co-authored-by: jackieyao0114 <[email protected]>
Co-authored-by: RevathiJambunathan <[email protected]>
Co-authored-by: Revathi  Jambunathan <[email protected]>
  • Loading branch information
5 people authored Aug 17, 2021
1 parent e37fdbb commit 0b03a3f
Show file tree
Hide file tree
Showing 2 changed files with 62 additions and 1 deletion.
4 changes: 3 additions & 1 deletion Examples/Tests/circuits/README
Original file line number Diff line number Diff line change
Expand Up @@ -22,4 +22,6 @@ inputs_current - provide a current source at the lo-y boundary inside the tube
inputs_current_large - provide a current source at the lo-y boundary in a region larger than the tube
(and zero the current source inside the tube)

inputs_silicon - same as inputs_2 but with the tube wall with -z normal replaced with silicon
inputs_silicon - same as inputs_2 but with the tube wall with -z normal replaced with silicon

inputs_longdomain - This is a 4um x 1.6cm x 4um domain with a current source with wavelength 1.6cm
59 changes: 59 additions & 0 deletions Examples/Tests/circuits/inputs_longdomain
Original file line number Diff line number Diff line change
@@ -0,0 +1,59 @@
################################
####### GENERAL PARAMETERS ######
#################################
max_step = 100000
amr.n_cell = 8 16000 8
amr.max_grid_size = 16000
amr.blocking_factor = 8
geometry.coord_sys = 0
geometry.prob_lo = -4.e-6 -0.8e-2 -4.e-6
geometry.prob_hi = 4.e-6 0.8e-2 4.e-6
amr.max_level = 0
# use pec instead of pml overlaying current source so you don't get a reflection
boundary.field_lo = periodic pec periodic
boundary.field_hi = periodic pml periodic

#################################
############ NUMERICS ###########
#################################
warpx.verbose = 1
warpx.use_filter = 0
warpx.cfl = 0.9

algo.em_solver_medium = macroscopic # vacuum/macroscopic

algo.macroscopic_sigma_method = backwardeuler # laxwendroff or backwardeuler

macroscopic.sigma_function(x,y,z) = "sigma_0"

macroscopic.epsilon_function(x,y,z) = "epsilon_0"

macroscopic.mu_function(x,y,z) = "mu_0 "

#################################
############ FIELDS #############
#################################

my_constants.pi = 3.14159265359
my_constants.c = 299792458.
my_constants.wavelength = 1.6e-2

my_constants.sigma_0 = 0.
my_constants.epsilon_0 = 8.8541878128e-12
my_constants.mu_0 = 1.25663706212e-06

warpx.B_excitation_on_grid_style = parse_B_excitation_grid_function

warpx.Bx_excitation_flag_function(x,y,z) = "(y <= -0.8e-2 + 1.e-6) * (y >= -0.8e-2 - 1.e-6)"
warpx.By_excitation_flag_function(x,y,z) = "0."
warpx.Bz_excitation_flag_function(x,y,z) = "0."

warpx.Bx_excitation_grid_function(x,y,z,t) = "cos((2*pi*x)/0.000016) * sin(2*pi*c/wavelength*t)"
warpx.By_excitation_grid_function(x,y,z,t) = "0."
warpx.Bz_excitation_grid_function(x,y,z,t) = "0."

# Diagnostics
diagnostics.diags_names = plt
plt.intervals = 100
plt.fields_to_plot = Ex Ey Ez Bx By Bz
plt.diag_type = Full

0 comments on commit 0b03a3f

Please sign in to comment.