aboutsummaryrefslogtreecommitdiff
path: root/Source/FieldSolver/FiniteDifferenceSolver
diff options
context:
space:
mode:
Diffstat (limited to 'Source/FieldSolver/FiniteDifferenceSolver')
-rw-r--r--Source/FieldSolver/FiniteDifferenceSolver/ApplySilverMuellerBoundary.cpp65
-rw-r--r--Source/FieldSolver/FiniteDifferenceSolver/FiniteDifferenceSolver.H4
2 files changed, 45 insertions, 24 deletions
diff --git a/Source/FieldSolver/FiniteDifferenceSolver/ApplySilverMuellerBoundary.cpp b/Source/FieldSolver/FiniteDifferenceSolver/ApplySilverMuellerBoundary.cpp
index 58a382d19..02b712679 100644
--- a/Source/FieldSolver/FiniteDifferenceSolver/ApplySilverMuellerBoundary.cpp
+++ b/Source/FieldSolver/FiniteDifferenceSolver/ApplySilverMuellerBoundary.cpp
@@ -37,7 +37,9 @@ void FiniteDifferenceSolver::ApplySilverMuellerBoundary (
std::array< std::unique_ptr<amrex::MultiFab>, 3 >& Efield,
std::array< std::unique_ptr<amrex::MultiFab>, 3 >& Bfield,
amrex::Box domain_box,
- amrex::Real const dt ) {
+ amrex::Real const dt,
+ amrex::Vector<int> field_boundary_lo,
+ amrex::Vector<int> field_boundary_hi) {
// Ensure that we are using the Yee solver
if (m_fdtd_algo != MaxwellSolverAlgo::Yee) {
@@ -67,6 +69,11 @@ void FiniteDifferenceSolver::ApplySilverMuellerBoundary (
int const nmodes = m_nmodes;
Real const rmin = m_rmin;
+ // Infer whether the Silver-Mueller needs to be applied in each direction
+ bool const apply_hi_r = (field_boundary_hi[0] == FieldBoundaryType::Absorbing_SilverMueller);
+ bool const apply_lo_z = (field_boundary_lo[1] == FieldBoundaryType::Absorbing_SilverMueller);
+ bool const apply_hi_z = (field_boundary_hi[1] == FieldBoundaryType::Absorbing_SilverMueller);
+
// tiling is usually set by TilingIfNotGPU()
// but here, we set it to false because of potential race condition,
// since we grow the tiles by one guard cell after creating them.
@@ -96,12 +103,12 @@ void FiniteDifferenceSolver::ApplySilverMuellerBoundary (
[=] AMREX_GPU_DEVICE (int i, int j, int /*k*/){
// At the +z boundary (innermost guard cell)
- if ( j==domain_box.bigEnd(1)+1 ){
+ if ( apply_hi_z && (j==domain_box.bigEnd(1)+1) ){
for (int m=0; m<2*nmodes-1; m++)
Br(i,j,0,m) = coef1_z*Br(i,j,0,m) - coef2_z*Et(i,j,0,m);
}
// At the -z boundary (innermost guard cell)
- if ( j==domain_box.smallEnd(1)-1 ){
+ if ( apply_lo_z && (j==domain_box.smallEnd(1)-1) ){
for (int m=0; m<2*nmodes-1; m++)
Br(i,j,0,m) = coef1_z*Br(i,j,0,m) + coef2_z*Et(i,j+1,0,m);
}
@@ -110,17 +117,17 @@ void FiniteDifferenceSolver::ApplySilverMuellerBoundary (
[=] AMREX_GPU_DEVICE (int i, int j, int /*k*/){
// At the +z boundary (innermost guard cell)
- if ( j==domain_box.bigEnd(1)+1 ){
+ if ( apply_hi_z && (j==domain_box.bigEnd(1)+1) ){
for (int m=0; m<2*nmodes-1; m++)
Bt(i,j,0,m) = coef1_z*Bt(i,j,0,m) + coef2_z*Er(i,j,0,m);
}
// At the -z boundary (innermost guard cell)
- if ( j==domain_box.smallEnd(1)-1 ){
+ if ( apply_lo_z && (j==domain_box.smallEnd(1)-1) ){
for (int m=0; m<2*nmodes-1; m++)
Bt(i,j,0,m) = coef1_z*Bt(i,j,0,m) - coef2_z*Er(i,j+1,0,m);
}
// At the +r boundary (innermost guard cell)
- if ( i==domain_box.bigEnd(0)+1 ){
+ if ( apply_hi_r && (i==domain_box.bigEnd(0)+1) ){
// Mode 0
Bt(i,j,0,0) = coef1_r*Bt(i,j,0,0) - coef2_r*Ez(i,j,0,0)
+ coef3_r*CylindricalYeeAlgorithm::UpwardDz(Er, coefs_z, n_coefs_z, i, j, 0, 0);
@@ -138,7 +145,7 @@ void FiniteDifferenceSolver::ApplySilverMuellerBoundary (
[=] AMREX_GPU_DEVICE (int i, int j, int /*k*/){
// At the +r boundary (innermost guard cell)
- if ( i==domain_box.bigEnd(0)+1 ){
+ if ( apply_hi_r && (i==domain_box.bigEnd(0)+1) ){
Real const r = rmin + (i + 0.5_rt)*dr; // r on nodal point (Bz is cell-centered in r)
// Mode 0
Bz(i,j,0,0) = coef1_r*Bz(i,j,0,0) + coef2_r*Et(i,j,0,0) - coef3_r*Et(i,j,0,0)/r;
@@ -170,6 +177,18 @@ void FiniteDifferenceSolver::ApplySilverMuellerBoundary (
amrex::Real const coef1_z = (1._rt - cdt_over_dz)/(1._rt + cdt_over_dz);
amrex::Real const coef2_z = 2._rt*cdt_over_dz/(1._rt + cdt_over_dz) / PhysConst::c;
+ bool const apply_lo_x = (field_boundary_lo[0] == FieldBoundaryType::Absorbing_SilverMueller);
+ bool const apply_hi_x = (field_boundary_hi[0] == FieldBoundaryType::Absorbing_SilverMueller);
+#ifdef WARPX_DIM_3D
+ bool const apply_lo_y = (field_boundary_lo[1] == FieldBoundaryType::Absorbing_SilverMueller);
+ bool const apply_hi_y = (field_boundary_hi[1] == FieldBoundaryType::Absorbing_SilverMueller);
+ bool const apply_lo_z = (field_boundary_lo[2] == FieldBoundaryType::Absorbing_SilverMueller);
+ bool const apply_hi_z = (field_boundary_hi[2] == FieldBoundaryType::Absorbing_SilverMueller);
+#else
+ bool const apply_lo_z = (field_boundary_lo[1] == FieldBoundaryType::Absorbing_SilverMueller);
+ bool const apply_hi_z = (field_boundary_hi[1] == FieldBoundaryType::Absorbing_SilverMueller);
+#endif
+
// Loop through the grids, and over the tiles within each grid
#ifdef AMREX_USE_OMP
#pragma omp parallel if (amrex::Gpu::notInLaunchRegion())
@@ -207,23 +226,23 @@ void FiniteDifferenceSolver::ApplySilverMuellerBoundary (
#ifdef WARPX_DIM_3D
// At the +y boundary (innermost guard cell)
- if ( j==domain_box.bigEnd(1)+1 )
+ if ( apply_hi_y && ( j==domain_box.bigEnd(1)+1 ) )
Bx(i,j,k) = coef1_y * Bx(i,j,k) + coef2_y * Ez(i,j,k);
// At the -y boundary (innermost guard cell)
- if ( j==domain_box.smallEnd(1)-1 )
+ if ( apply_lo_y && ( j==domain_box.smallEnd(1)-1 ) )
Bx(i,j,k) = coef1_y * Bx(i,j,k) - coef2_y * Ez(i,j+1,k);
// At the +z boundary (innermost guard cell)
- if ( k==domain_box.bigEnd(2)+1 )
+ if ( apply_hi_z && ( k==domain_box.bigEnd(2)+1 ) )
Bx(i,j,k) = coef1_z * Bx(i,j,k) - coef2_z * Ey(i,j,k);
// At the -z boundary (innermost guard cell)
- if ( k==domain_box.smallEnd(2)-1 )
+ if ( apply_lo_z && ( k==domain_box.smallEnd(2)-1 ) )
Bx(i,j,k) = coef1_z * Bx(i,j,k) + coef2_z * Ey(i,j,k+1);
#elif WARPX_DIM_XZ
// At the +z boundary (innermost guard cell)
- if ( j==domain_box.bigEnd(1)+1 )
+ if ( apply_hi_z && ( j==domain_box.bigEnd(1)+1 ) )
Bx(i,j,k) = coef1_z * Bx(i,j,k) - coef2_z * Ey(i,j,k);
// At the -z boundary (innermost guard cell)
- if ( j==domain_box.smallEnd(1)-1 )
+ if ( apply_lo_z && ( j==domain_box.smallEnd(1)-1 ) )
Bx(i,j,k) = coef1_z * Bx(i,j,k) + coef2_z * Ey(i,j+1,k);
#endif
},
@@ -232,24 +251,24 @@ void FiniteDifferenceSolver::ApplySilverMuellerBoundary (
[=] AMREX_GPU_DEVICE (int i, int j, int k){
// At the +x boundary (innermost guard cell)
- if ( i==domain_box.bigEnd(0)+1 )
+ if ( apply_hi_x && ( i==domain_box.bigEnd(0)+1 ) )
By(i,j,k) = coef1_x * By(i,j,k) - coef2_x * Ez(i,j,k);
// At the -x boundary (innermost guard cell)
- if ( i==domain_box.smallEnd(0)-1 )
+ if ( apply_lo_x && ( i==domain_box.smallEnd(0)-1 ) )
By(i,j,k) = coef1_x * By(i,j,k) + coef2_x * Ez(i+1,j,k);
#ifdef WARPX_DIM_3D
// At the +z boundary (innermost guard cell)
- if ( k==domain_box.bigEnd(2)+1 )
+ if ( apply_hi_z && ( k==domain_box.bigEnd(2)+1 ) )
By(i,j,k) = coef1_z * By(i,j,k) + coef2_z * Ex(i,j,k);
// At the -z boundary (innermost guard cell)
- if ( k==domain_box.smallEnd(2)-1 )
+ if ( apply_lo_z && ( k==domain_box.smallEnd(2)-1 ) )
By(i,j,k) = coef1_z * By(i,j,k) - coef2_z * Ex(i,j,k+1);
#elif WARPX_DIM_XZ
// At the +z boundary (innermost guard cell)
- if ( j==domain_box.bigEnd(1)+1 )
+ if ( apply_hi_z && ( j==domain_box.bigEnd(1)+1 ) )
By(i,j,k) = coef1_z * By(i,j,k) + coef2_z * Ex(i,j,k);
// At the -z boundary (innermost guard cell)
- if ( j==domain_box.smallEnd(1)-1 )
+ if ( apply_lo_z && ( j==domain_box.smallEnd(1)-1 ) )
By(i,j,k) = coef1_z * By(i,j,k) - coef2_z * Ex(i,j+1,k);
#endif
},
@@ -258,17 +277,17 @@ void FiniteDifferenceSolver::ApplySilverMuellerBoundary (
[=] AMREX_GPU_DEVICE (int i, int j, int k){
// At the +x boundary (innermost guard cell)
- if ( i==domain_box.bigEnd(0)+1 )
+ if ( apply_hi_x && ( i==domain_box.bigEnd(0)+1 ) )
Bz(i,j,k) = coef1_x * Bz(i,j,k) + coef2_x * Ey(i,j,k);
// At the -x boundary (innermost guard cell)
- if ( i==domain_box.smallEnd(0)-1 )
+ if ( apply_lo_x && ( i==domain_box.smallEnd(0)-1 ) )
Bz(i,j,k) = coef1_x * Bz(i,j,k) - coef2_x * Ey(i+1,j,k);
#ifdef WARPX_DIM_3D
// At the +y boundary (innermost guard cell)
- if ( j==domain_box.bigEnd(1)+1 )
+ if ( apply_hi_y && ( j==domain_box.bigEnd(1)+1 ) )
Bz(i,j,k) = coef1_y * Bz(i,j,k) - coef2_y * Ex(i,j,k);
// At the -y boundary (innermost guard cell)
- if ( j==domain_box.smallEnd(1)-1 )
+ if ( apply_lo_y && ( j==domain_box.smallEnd(1)-1 ) )
Bz(i,j,k) = coef1_y * Bz(i,j,k) + coef2_y * Ex(i,j+1,k);
#endif
}
diff --git a/Source/FieldSolver/FiniteDifferenceSolver/FiniteDifferenceSolver.H b/Source/FieldSolver/FiniteDifferenceSolver/FiniteDifferenceSolver.H
index bd6add3f0..820f2b876 100644
--- a/Source/FieldSolver/FiniteDifferenceSolver/FiniteDifferenceSolver.H
+++ b/Source/FieldSolver/FiniteDifferenceSolver/FiniteDifferenceSolver.H
@@ -72,7 +72,9 @@ class FiniteDifferenceSolver
std::array< std::unique_ptr<amrex::MultiFab>, 3 >& Efield,
std::array< std::unique_ptr<amrex::MultiFab>, 3 >& Bfield,
amrex::Box domain_box,
- amrex::Real const dt );
+ amrex::Real const dt,
+ amrex::Vector<int> field_boundary_lo,
+ amrex::Vector<int> field_boundary_hi);
void ComputeDivE ( const std::array<std::unique_ptr<amrex::MultiFab>,3>& Efield,
amrex::MultiFab& divE );
generated by cgit v1.2.3 (git 2.46.0) at 2025-08-14 01:00:44 +0000