diff options
Diffstat (limited to 'Examples')
| -rwxr-xr-x | Examples/Tests/ElectrostaticSphereEB/analysis_rz.py | 67 | ||||
| -rw-r--r-- | Examples/Tests/ElectrostaticSphereEB/inputs_rz | 29 |
2 files changed, 96 insertions, 0 deletions
diff --git a/Examples/Tests/ElectrostaticSphereEB/analysis_rz.py b/Examples/Tests/ElectrostaticSphereEB/analysis_rz.py new file mode 100755 index 000000000..bb170f921 --- /dev/null +++ b/Examples/Tests/ElectrostaticSphereEB/analysis_rz.py @@ -0,0 +1,67 @@ +#!/usr/bin/env python + +# Copyright 2019-2021 Yinjian Zhao +# +# This file is part of WarpX. +# +# License: BSD-3-Clause-LBNL + +# This script tests the embedded boundary in RZ. +# A cylindrical surface (r=0.1) has a fixed potential 1 V. +# The outer surface has 0 V fixed. +# Thus the analytical solution has the form: +# phi(r) = A+B*log(r), Er(r) = -B/r. + +# Possible errors: +# tolerance: 0.004 +# Possible running time: < 1 s + +import os +import sys + +import numpy as np +import yt + +sys.path.insert(1, '../../../../warpx/Regression/Checksum/') +import checksumAPI + +tolerance = 0.004 + +fn = sys.argv[1] +ds = yt.load( fn ) + +all_data_level_0 = ds.covering_grid(level=0,left_edge=ds.domain_left_edge, dims=ds.domain_dimensions) +phi = all_data_level_0['boxlib', 'phi'].v.squeeze() +Er = all_data_level_0['boxlib', 'Ex'].v.squeeze() + +Dx = ds.domain_width/ds.domain_dimensions +dr = Dx[0] +rmin = ds.domain_left_edge[0] +rmax = ds.domain_right_edge[0] +nr = phi.shape[0] +r = np.linspace(rmin+dr/2.,rmax-dr/2.,nr) +B = 1.0/np.log(0.1/0.5) +A = -B*np.log(0.5) + +err = 0.0 +errmax_phi = 0.0 +errmax_Er = 0.0 +for i in range(len(r)): + if r[i]>=0.1: + phi_theory = A+B*np.log(r[i]) + Er_theory = -B/float(r[i]) + err = abs( phi_theory - phi[i,:] ).max() / phi_theory + if err>errmax_phi: + errmax_phi = err + err = abs( Er_theory - Er[i,:] ).max() / Er_theory + # Exclude the last inaccurate interpolation. + if err>errmax_Er and i<len(r)-1: + errmax_Er = err + +print('max error of phi = ', errmax_phi) +print('max error of Er = ', errmax_Er) +print('tolerance = ', tolerance) +assert(errmax_phi<tolerance and errmax_Er<tolerance) + +test_name = os.path.split(os.getcwd())[1] +checksumAPI.evaluate_checksum(test_name, fn, do_particles=False) diff --git a/Examples/Tests/ElectrostaticSphereEB/inputs_rz b/Examples/Tests/ElectrostaticSphereEB/inputs_rz new file mode 100644 index 000000000..af5f7fcbf --- /dev/null +++ b/Examples/Tests/ElectrostaticSphereEB/inputs_rz @@ -0,0 +1,29 @@ +max_step = 1 +amr.n_cell = 64 64 +amr.max_level = 0 +amr.blocking_factor = 8 +amr.max_grid_size = 128 +boundary.field_lo = none periodic +boundary.field_hi = pec periodic +boundary.potential_lo_x = 0 +boundary.potential_hi_x = 0 +boundary.potential_lo_y = 0 +boundary.potential_hi_y = 0 +boundary.potential_lo_z = 0 +boundary.potential_hi_z = 0 +geometry.dims = RZ +geometry.prob_lo = 0.0 -0.5 +geometry.prob_hi = 0.5 0.5 +warpx.const_dt = 1e-6 + +warpx.do_electrostatic = labframe +warpx.eb_implicit_function = "-(x**2-0.1**2)" +warpx.eb_potential(x,y,z,t) = "1." +warpx.self_fields_required_precision = 1.e-7 + +algo.field_gathering = momentum-conserving + +diagnostics.diags_names = diag1 +diag1.intervals = 1 +diag1.diag_type = Full +diag1.fields_to_plot = Ex phi |