colin/fdtd-coremem
1
0
Fork 0
Code Issues Pull Requests Projects Releases Wiki Activity

Improve the UPML parameters by properly summing the axis conductivities in the corners

Browse Source
This commit is contained in:
Colin
2020-09-19 20:10:10 -07:00
parent f6e4ed8cb9
commit d713c79bbd
2 changed files with 28 additions and 18 deletions
Download Patch File Download Diff File Expand all files Collapse all files

4
examples/boundary_conditions.rs
View File

@@ -3,9 +3,9 @@ use coremem::{Driver, mat};
fn main() { fn main() {
coremem::init_logging(); coremem::init_logging();
let width = 201; let width = 201;
let boundary = 20; let boundary = 50;
let mut driver = Driver::new(width, width, 1e-3 /* feature size */); let mut driver = Driver::new(width, width, 1e-3 /* feature size */);
driver.add_y4m_renderer(&*format!("boundary_conditions_upml-{}.y4m", boundary)); driver.add_y4m_renderer(&*format!("boundary_conditions_upml_round-{}.y4m", boundary));
// driver.add_term_renderer(); // driver.add_term_renderer();
//driver.add_boundary(boundary, 0.1); //driver.add_boundary(boundary, 0.1);

42
src/driver.rs
View File

@@ -76,22 +76,32 @@ impl Driver {
pub fn add_upml_boundary(&mut self, thickness: u32) { pub fn add_upml_boundary(&mut self, thickness: u32) {
// Based on explanation here (slide 63): https://empossible.net/wp-content/uploads/2020/01/Lecture-The-Perfectly-Matched-Layer.pdf // Based on explanation here (slide 63): https://empossible.net/wp-content/uploads/2020/01/Lecture-The-Perfectly-Matched-Layer.pdf
for inset in 0..thickness { let h = self.state.height();
let depth = thickness - inset; let w = self.state.width();
let conductivity = 0.5 * consts::EPS0 / self.state.timestep() * (depth as Flt/thickness as Flt).powf(3.0); for y in 0..h {
for x in inset..self.state.width() - inset { for x in 0..w {
let conductor = mat::Static::anisotropic_conductor(Vec3::new(0.0, conductivity, 0.0)); let depth_x = if x < thickness {
// left thickness - x
*self.state.get_mut(x, inset).mat_mut() = conductor.clone().into(); } else if x >= w - thickness {
// right w - x
*self.state.get_mut(x, self.state.height() - 1 - inset).mat_mut() = conductor.into(); } else {
} 0
for y in inset..self.state.height() - inset { };
let conductor = mat::Static::anisotropic_conductor(Vec3::new(conductivity, 0.0, 0.0)); let depth_y = if y < thickness {
// top thickness - y
*self.state.get_mut(inset, y).mat_mut() = conductor.clone().into(); } else if y >= h - thickness {
// bottom h - y
*self.state.get_mut(self.state.width() - 1 - inset, y).mat_mut() = conductor.into(); } else {
0
};
if depth_x > 0 || depth_y > 0 {
let scale = 0.5 * consts::EPS0 / self.state.timestep();
let cond_x = scale * (depth_x as Flt/thickness as Flt).powf(3.0);
let cond_y = scale * (depth_y as Flt/thickness as Flt).powf(3.0);
let conductor = mat::Static::anisotropic_conductor(Vec3::new(cond_x, cond_y, 0.0));
*self.state.get_mut(x, y).mat_mut() = conductor.into();
}
} }
} }
} }