Skip to content

Commit

Permalink
update mesh and add exodus
Browse files Browse the repository at this point in the history
  • Loading branch information
lin-yang-ly committed Oct 25, 2024
1 parent 5fa5925 commit 5e0dce1
Show file tree
Hide file tree
Showing 4 changed files with 3 additions and 37 deletions.
Binary file modified test/tests/val-2a/gold/val-2a_limited_out.e
Binary file not shown.
Binary file added test/tests/val-2a/gold/val-2a_out.e
Binary file not shown.
4 changes: 2 additions & 2 deletions test/tests/val-2a/tests
Original file line number Diff line number Diff line change
Expand Up @@ -5,14 +5,14 @@
[val-2a_csvdiff]
type = CSVDiff
input = 'val-2a.i'
cli_args = "simulation_time='${units 3e3 s}' Mesh/active='cartesian' inactive='Adaptivity' Executioner/nl_rel_tol=8e-1 Executioner/dt=2.5 Executioner/inactive=TimeStepper Outputs/file_base='val-2a_limited_out'"
cli_args = "simulation_time='${units 3e3 s}' nx_scale=1 Executioner/nl_rel_tol=8e-1 Executioner/dt=2.5 Executioner/inactive=TimeStepper Outputs/file_base='val-2a_limited_out'"
csvdiff = val-2a_limited_out.csv
requirement = 'The system shall be able to model Ion implantation based on the requirements and conditions listed in the TMAP4 V&V report'
[]
[val-2a_exodiff]
type = Exodiff
input = 'val-2a.i'
cli_args = "simulation_time='${units 3e3 s}' Mesh/active='cartesian' inactive='Adaptivity' Executioner/nl_rel_tol=8e-1 Executioner/dt=2.5 Executioner/inactive=TimeStepper Outputs/file_base='val-2a_limited_out'"
cli_args = "simulation_time='${units 3e3 s}' nx_scale=1 Executioner/nl_rel_tol=8e-1 Executioner/dt=2.5 Executioner/inactive=TimeStepper Outputs/file_base='val-2a_limited_out'"
exodiff = val-2a_limited_out.e
prereq = val-2a_csvdiff
should_execute = false # this test relies on the output files from val-2a_csvdiff_heavy, so it shouldn't be run twice
Expand Down
36 changes: 1 addition & 35 deletions test/tests/val-2a/val-2a.i
Original file line number Diff line number Diff line change
@@ -1,4 +1,4 @@
nx_scale = 1
nx_scale = 5
high_dt_max = 100
low_dt_max = 1
simulation_time = '${units 2e4 s}'
Expand All @@ -16,15 +16,6 @@ time_3 = '${units 12062 s}'
time_4 = '${units 14572 s}'
time_5 = '${units 17678 s}'

## Modeling parameters
sample_thickness = '${units 5e-4 m -> nm}'
desired_mesh_size_surface = '${units 1e-9 m -> nm}' # defined by the refinement of the implantation profile and the need to have a fine mesh at BC
refinement_depth_left = '${units 1e-5 m -> nm}'
largest_mesh_size = ${refinement_depth_left}
refinement_steps = ${fparse floor((log(largest_mesh_size)-log(desired_mesh_size_surface)) / log(2))} # design to reach desired_mesh_size_surface and refine refinement_depth_left from the start (otherwise refinment does not happen)
node_size = ${fparse desired_mesh_size_surface*2^refinement_steps} # defined by the refinement of the implantation profile
num_nodes = ${fparse floor(sample_thickness/node_size)}

[Variables]
[concentration]
order = FIRST
Expand All @@ -33,7 +24,6 @@ num_nodes = ${fparse floor(sample_thickness/node_size)}
[]

[Mesh]
active = 'gen'
[cartesian]
type = CartesianMeshGenerator
dim = 1
Expand All @@ -43,29 +33,6 @@ num_nodes = ${fparse floor(sample_thickness/node_size)}
ix = '${fparse 5 * ${nx_scale}} ${nx_scale} ${nx_scale}
${nx_scale} ${nx_scale} ${fparse 10 * ${nx_scale}}'
[]

[gen]
type = GeneratedMeshGenerator
dim = 1
nx = ${num_nodes}
xmax = ${sample_thickness}
[]
[]

[Adaptivity]
steps = ${refinement_steps}
initial_steps = ${refinement_steps}
max_h_level = ${refinement_steps}
marker = box
[Markers]
[box]
type = BoxMarker
bottom_left = '${refinement_depth_left} -1 -1'
top_right = '${sample_thickness} 1 1'
inside = DO_NOTHING
outside = REFINE
[]
[]
[]

[Kernels]
Expand Down Expand Up @@ -142,7 +109,6 @@ num_nodes = ${fparse floor(sample_thickness/node_size)}
[]

[Functions]
################# TMAP4
[Kd_left_func]
type = ParsedFunction
expression = '${dissociation_coefficient_parameter_enclos1} * (1 - 0.9999 * exp(-6e-5 * t))'
Expand Down

0 comments on commit 5e0dce1

Please sign in to comment.