- Notifications
You must be signed in to change notification settings - Fork269
Open
Description
I noticed a severe bug in the PML when using steel boxes that makes it reflect the EM waves. The bug happens consistently when using the same geometry, but not all geometries using the material show the bug, in my tests it happens almost half of the times. The steel is implemented as#material: 1.0 1400000.0 1000.0 0.0 steel_sleepers. The reflections seem to be stronger when a Peplinski soil is also present in the geometry, didn't try with other materials.
Some field snapshots:
With steel boxes and a Peplinski soil:
Only steel boxes:
This bug makes entire B-scans completely wrong:
With bug:
Similar geometry without bug:
Code to reproduce:
## General commands:#title: scan_00185 #domain: 1.5 1.7 0.002 #dx_dy_dz: 0.002 0.002 0.002 #time_window: 2.5e-08 #output_dir: pml_bug#messages: n ## Source and receiver:#waveform: ricker 1 1000000000.0 source_wave #hertzian_dipole: z 0.65 1.5 0.001 source_wave #rx: 0.85 1.5 0.001 #src_steps: 0.005 0.0 0.0 #rx_steps: 0.005 0.0 0.0 ## PSS:#soil_peplinski: 0.42592250416541744 0.13647460641294687 2.3 2.65 0.06364996446113151 0.06494363641117931 pss #fractal_box: 0 0.4814529318664053 0 1.5 0.7935274720150615 0.002 1.5 1 1 1 15 pss pss_fractal_box 921236401 #add_surface_roughness: 0 0.7935274720150615 0 1.5 0.7935274720150615 0.002 1.5 1 1 0.7835274720150615 0.7935274720150615 pss_fractal_box 2626582198 ## Sleepers:#material: 1.0 1400000.0 1000.0 0.0 steel_sleepers #box: 0.0 0.9349999999999999 0.0 0.28 1.035 0.002 steel_sleepers #box: 0.61 0.9349999999999999 0.0 0.9299999999999999 1.035 0.002 steel_sleepers #box: 1.26 0.9349999999999999 0.0 1.5 1.035 0.002 steel_sleepers ## Snapshots:#python:snapshot_times = [i*1e-10 for i in range(1, 50)]for t in snapshot_times: print(f"#snapshot: 0 0 0 1.5 1.7 0.002 0.002 0.002 0.002 {t} {'pml_bug_' + str(t)}")#end_python:Metadata
Metadata
Assignees
Labels
No labels