Commit 86b9d551 authored by Alexander Nikolas Breuer's avatar Alexander Nikolas Breuer
Browse files

Updated layer of halfspace 3 setup.

parent 32d2d95f
Layer Over Halfspace Benchmark 3 (LOH.3)
========================================
This directory contains the data for the LOH.3 benchmark.
The benchmark has a planar free surface and internal material contrast.
The benchmarks assesses the accuracy of constant Q.
Layer Over Halfspace 3 (LOH.3)
==============================
* Contact: Alexander Breuer
The meshing was done on a local desktop.
The simulations on a single socket of the Frontera machine.
Installation and Meshing
------------------------
.. code-block:: bash
mkdir logs
CC=gcc CXX=g++ ./scripts/install.sh 2>&1 | tee logs/install.log
./scripts/meshes.sh 2>&1 | tee logs/meshes.log
Runs
----
.. code-block:: bash
sbatch scripts/run.slurm
Visualization
-------------
.. code-block:: bash
mkdir vis
cd vis
wget http://www.sismowine.org/WP2_LOH3/ref.zip
unzip ref.zip
ln -s ../runs/loh3_uns_* .
python ../scripts/vis.py 2>&1 | tee ../logs/vis.log
Config for the LOH.3 Benchmark
==============================
This directory contains the config for the LOH.3 benchmark.
* Contact: Alex Breuer
<?xml version="1.0" encoding="UTF-8"?>
<!--
@file This file is part of EDGE.
@author Alexander Breuer (anbreuer AT ucsd.edu)
@section LICENSE
Copyright (c) 2021, Friedrich Schiller University Jena
Copyright (c) 2016-2017, Regents of the University of California
All rights reserved.
Redistribution and use in source and binary forms, with or without modification, are permitted provided that the following conditions are met:
1. Redistributions of source code must retain the above copyright notice, this list of conditions and the following disclaimer.
2. Redistributions in binary form must reproduce the above copyright notice, this list of conditions and the following disclaimer in the documentation and/or other materials provided with the distribution.
3. Neither the name of the copyright holder nor the names of its contributors may be used to endorse or promote products derived from this software without specific prior written permission.
THIS SOFTWARE IS PROVIDED BY THE COPYRIGHT HOLDERS AND CONTRIBUTORS "AS IS" AND ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED TO, THE IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR PURPOSE ARE DISCLAIMED. IN NO EVENT SHALL THE COPYRIGHT HOLDER OR CONTRIBUTORS BE LIABLE FOR ANY DIRECT, INDIRECT, INCIDENTAL, SPECIAL, EXEMPLARY, OR CONSEQUENTIAL DAMAGES (INCLUDING, BUT NOT LIMITED TO, PROCUREMENT OF SUBSTITUTE GOODS OR SERVICES; LOSS OF USE, DATA, OR PROFITS; OR BUSINESS INTERRUPTION) HOWEVER CAUSED AND ON ANY THEORY OF LIABILITY, WHETHER IN CONTRACT, STRICT LIABILITY, OR TORT (INCLUDING NEGLIGENCE OR OTHERWISE) ARISING IN ANY WAY OUT OF THE USE OF THIS SOFTWARE, EVEN IF ADVISED OF THE POSSIBILITY OF SUCH DAMAGE.
@section DESCRIPTION
Runtime config of the Layer Over Halfspace benchmark 3 (LOH3).
-->
<edge>
<cfr>
<mesh>
<in>
<base>meshes/gen/edge_v/loh3_100</base>
<extension>.msh</extension>
</in>
<boundary>
<free_surface>101</free_surface>
<outflow>105</outflow>
</boundary>
</mesh>
<velocity_model>
<domain>
<half_space>
<origin>
<x>0</x>
<y>0</y>
<z>0</z>
</origin>
<normal>
<x>0</x>
<y>0</y>
<z>1</z>
</normal>
</half_space>
<half_space>
<origin>
<x> 0</x>
<y> 0</y>
<z>1000</z>
</origin>
<normal>
<x> 0</x>
<y> 0</y>
<z>-1</z>
</normal>
</half_space>
<rho>2600</rho>
<lambda>20.8E9</lambda>
<mu>10.4E9</mu>
<qp>40</qp>
<qs>120</qs>
</domain>
<domain>
<half_space>
<origin>
<x> 0</x>
<y> 0</y>
<z>100000</z>
</origin>
<normal>
<x> 0</x>
<y> 0</y>
<z>-1</z>
</normal>
</half_space>
<rho>2700</rho>
<lambda>32.4038016E9</lambda>
<mu>32.3980992E9</mu>
<qp>69.3</qp>
<qs>155.9</qs>
</domain>
</velocity_model>
<setups>
<attenuation>
<central_frequency>2.5</central_frequency>
<frequency_ratio>100</frequency_ratio>
</attenuation>
<point_sources>
<file>sources/loh3_uns_src.h5</file>
</point_sources>
<end_time>9.0</end_time>
</setups>
<output>
<receivers>
<path_to_dir>loh3_uns_100</path_to_dir>
<freq>0.005</freq>
<receiver>
<name>loh3_1</name>
<coords>
<x> 0</x>
<y> 693</y>
<z> 0</z>
</coords>
</receiver>
<receiver>
<name>loh3_2</name>
<coords>
<x> 0</x>
<y>5543</y>
<z> 0</z>
</coords>
</receiver>
<receiver>
<name>loh3_3</name>
<coords>
<x> 0</x>
<y>10392</y>
<z> 0</z>
</coords>
</receiver>
<receiver>
<name>loh3_4</name>
<coords>
<x> 490</x>
<y> 490</y>
<z> 0</z>
</coords>
</receiver>
<receiver>
<name>loh3_5</name>
<coords>
<x>3919</x>
<y>3919</y>
<z> 0</z>
</coords>
</receiver>
<receiver>
<name>loh3_6</name>
<coords>
<x>7348</x>
<y>7348</y>
<z> 0</z>
</coords>
</receiver>
<receiver>
<name>loh3_7</name>
<coords>
<x> 577</x>
<y> 384</y>
<z> 0</z>
</coords>
</receiver>
<receiver>
<name>loh3_8</name>
<coords>
<x>4612</x>
<y>3075</y>
<z> 0</z>
</coords>
</receiver>
<receiver>
<name>loh3_9</name>
<coords>
<x>8647</x>
<y>5764</y>
<z> 0</z>
</coords>
</receiver>
</receivers>
<wave_field>
<type>vtk_binary</type>
<sparse_type>101</sparse_type>
<file>loh3_uns_100/wf</file>
<int>1.0</int>
</wave_field>
</output>
</cfr>
</edge>
......@@ -5,7 +5,7 @@
@author Alexander Breuer (anbreuer AT ucsd.edu)
@section LICENSE
Copyright (c) 2019, Alexander Breuer
Copyright (c) 2021, Friedrich Schiller University Jena
Copyright (c) 2016-2017, Regents of the University of California
All rights reserved.
......@@ -20,22 +20,15 @@
THIS SOFTWARE IS PROVIDED BY THE COPYRIGHT HOLDERS AND CONTRIBUTORS "AS IS" AND ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED TO, THE IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR PURPOSE ARE DISCLAIMED. IN NO EVENT SHALL THE COPYRIGHT HOLDER OR CONTRIBUTORS BE LIABLE FOR ANY DIRECT, INDIRECT, INCIDENTAL, SPECIAL, EXEMPLARY, OR CONSEQUENTIAL DAMAGES (INCLUDING, BUT NOT LIMITED TO, PROCUREMENT OF SUBSTITUTE GOODS OR SERVICES; LOSS OF USE, DATA, OR PROFITS; OR BUSINESS INTERRUPTION) HOWEVER CAUSED AND ON ANY THEORY OF LIABILITY, WHETHER IN CONTRACT, STRICT LIABILITY, OR TORT (INCLUDING NEGLIGENCE OR OTHERWISE) ARISING IN ANY WAY OUT OF THE USE OF THIS SOFTWARE, EVEN IF ADVISED OF THE POSSIBILITY OF SUCH DAMAGE.
@section DESCRIPTION
Runtime config for the Layer Over Halfspace benchmark 3 (LOH3).
Runtime config of the Layer Over Halfspace benchmark 3 (LOH3).
-->
<edge>
<build>
<element_type>tet4</element_type>
<moab>yes</moab>
</build>
<cfr>
<mesh>
<options>
<read>;</read>
</options>
<files>
<in>meshes/tet4_200.msh</in>
</files>
<in>
<base>meshes/gen/edge_v/loh3_200</base>
<extension>.msh</extension>
</in>
<boundary>
<free_surface>101</free_surface>
<outflow>105</outflow>
......@@ -91,8 +84,8 @@
</half_space>
<rho>2700</rho>
<lambda>32.4E9</lambda>
<mu>32.4E9</mu>
<lambda>32.4038016E9</lambda>
<mu>32.3980992E9</mu>
<qp>69.3</qp>
<qs>155.9</qs>
</domain>
......@@ -111,6 +104,7 @@
<end_time>9.0</end_time>
</setups>
<output>
<receivers>
<path_to_dir>loh3_uns_200</path_to_dir>
......@@ -121,7 +115,7 @@
<coords>
<x> 0</x>
<y> 693</y>
<z>0</z>
<z> 0</z>
</coords>
</receiver>
......@@ -130,7 +124,7 @@
<coords>
<x> 0</x>
<y>5543</y>
<z>0</z>
<z> 0</z>
</coords>
</receiver>
......@@ -139,7 +133,7 @@
<coords>
<x> 0</x>
<y>10392</y>
<z> 0.01</z>
<z> 0</z>
</coords>
</receiver>
......@@ -148,7 +142,7 @@
<coords>
<x> 490</x>
<y> 490</y>
<z>0</z>
<z> 0</z>
</coords>
</receiver>
......@@ -157,7 +151,7 @@
<coords>
<x>3919</x>
<y>3919</y>
<z>0</z>
<z> 0</z>
</coords>
</receiver>
......@@ -166,7 +160,7 @@
<coords>
<x>7348</x>
<y>7348</y>
<z>0</z>
<z> 0</z>
</coords>
</receiver>
......@@ -175,7 +169,7 @@
<coords>
<x> 577</x>
<y> 384</y>
<z>0</z>
<z> 0</z>
</coords>
</receiver>
......@@ -184,7 +178,7 @@
<coords>
<x>4612</x>
<y>3075</y>
<z>0</z>
<z> 0</z>
</coords>
</receiver>
......@@ -193,7 +187,7 @@
<coords>
<x>8647</x>
<y>5764</y>
<z>0</z>
<z> 0</z>
</coords>
</receiver>
......@@ -205,7 +199,6 @@
<file>loh3_uns_200/wf</file>
<int>1.0</int>
</wave_field>
</output>
</cfr>
</edge>
This source diff could not be displayed because it is stored in LFS. You can view the blob instead.
This source diff could not be displayed because it is stored in LFS. You can view the blob instead.
......@@ -4,6 +4,6 @@
To the extent possible under law, `Alexander Nikolas Breuer <https://dial3343.org/>`_ has waived all copyright and related or neighboring rights to the following work:
* `tet4_200.msh.tar.xz`
* `gen.tar.xz`
This work is published from: Germany.
Mesh for the LOH.3 Benchmark
============================
This directory contains the meshes for the LOH.3 benchmark.
* Contact: Alex Breuer
Used Commands
=============
.. code-block:: bash
gmsh -format msh2 -3 tet4_200.geo 2>&1 | tee tet4_200.msh.log
gmsh -3 tet4_200.msh -o /tmp/tet4_200.vtk
mesh_probe /tmp/tet4_200.vtk > tet4_200.csv
<?xml version="1.0" encoding="UTF-8"?>
<!--
@file This file is part of EDGE.
@author Alexander Breuer (anbreuer AT ucsd.edu)
@section LICENSE
Copyright (c) 2021, Friedrich Schiller University Jena
All rights reserved.
Redistribution and use in source and binary forms, with or without modification, are permitted provided that the following conditions are met:
1. Redistributions of source code must retain the above copyright notice, this list of conditions and the following disclaimer.
2. Redistributions in binary form must reproduce the above copyright notice, this list of conditions and the following disclaimer in the documentation and/or other materials provided with the distribution.
3. Neither the name of the copyright holder nor the names of its contributors may be used to endorse or promote products derived from this software without specific prior written permission.
THIS SOFTWARE IS PROVIDED BY THE COPYRIGHT HOLDERS AND CONTRIBUTORS "AS IS" AND ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED TO, THE IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR PURPOSE ARE DISCLAIMED. IN NO EVENT SHALL THE COPYRIGHT HOLDER OR CONTRIBUTORS BE LIABLE FOR ANY DIRECT, INDIRECT, INCIDENTAL, SPECIAL, EXEMPLARY, OR CONSEQUENTIAL DAMAGES (INCLUDING, BUT NOT LIMITED TO, PROCUREMENT OF SUBSTITUTE GOODS OR SERVICES; LOSS OF USE, DATA, OR PROFITS; OR BUSINESS INTERRUPTION) HOWEVER CAUSED AND ON ANY THEORY OF LIABILITY, WHETHER IN CONTRACT, STRICT LIABILITY, OR TORT (INCLUDING NEGLIGENCE OR OTHERWISE) ARISING IN ANY WAY OUT OF THE USE OF THIS SOFTWARE, EVEN IF ADVISED OF THE POSSIBILITY OF SUCH DAMAGE.
@section DESCRIPTION
EDGE-V config for the Layer Over Halfspace benchmark 3 (LOH1).
-->
<edge_v>
<mesh>
<periodic>no</periodic>
<write_element_annotations>1</write_element_annotations>
<n_partitions>1</n_partitions>
<files>
<in>gen/gmsh/tet4_100.msh</in>
<out>
<base>gen/edge_v/loh3_100</base>
<extension>.msh</extension>
</out>
</files>
</mesh>
<velocity_model>
<seismic_expression>
<![CDATA[
if( z < 1000 ) {
vp := 4000;
vs := 2000;
qp := 40;
qs := 120;
}
else {
vp := 6000;
vs := 3464;
qp := 69.3;
qs := 155.9;
};
]]>
</seismic_expression>
</velocity_model>
<time>
<n_groups>5</n_groups>
<centroid_based_cfl>true</centroid_based_cfl>
<files>
<out>
<time_steps>gen/edge_v/time_steps_100.csv</time_steps>
</out>
</files>
</time>
</edge_v>
<?xml version="1.0" encoding="UTF-8"?>
<!--
@file This file is part of EDGE.
@author Alexander Breuer (anbreuer AT ucsd.edu)
@section LICENSE
Copyright (c) 2021, Friedrich Schiller University Jena
All rights reserved.
Redistribution and use in source and binary forms, with or without modification, are permitted provided that the following conditions are met:
1. Redistributions of source code must retain the above copyright notice, this list of conditions and the following disclaimer.
2. Redistributions in binary form must reproduce the above copyright notice, this list of conditions and the following disclaimer in the documentation and/or other materials provided with the distribution.
3. Neither the name of the copyright holder nor the names of its contributors may be used to endorse or promote products derived from this software without specific prior written permission.
THIS SOFTWARE IS PROVIDED BY THE COPYRIGHT HOLDERS AND CONTRIBUTORS "AS IS" AND ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED TO, THE IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR PURPOSE ARE DISCLAIMED. IN NO EVENT SHALL THE COPYRIGHT HOLDER OR CONTRIBUTORS BE LIABLE FOR ANY DIRECT, INDIRECT, INCIDENTAL, SPECIAL, EXEMPLARY, OR CONSEQUENTIAL DAMAGES (INCLUDING, BUT NOT LIMITED TO, PROCUREMENT OF SUBSTITUTE GOODS OR SERVICES; LOSS OF USE, DATA, OR PROFITS; OR BUSINESS INTERRUPTION) HOWEVER CAUSED AND ON ANY THEORY OF LIABILITY, WHETHER IN CONTRACT, STRICT LIABILITY, OR TORT (INCLUDING NEGLIGENCE OR OTHERWISE) ARISING IN ANY WAY OUT OF THE USE OF THIS SOFTWARE, EVEN IF ADVISED OF THE POSSIBILITY OF SUCH DAMAGE.
@section DESCRIPTION
EDGE-V config for the Layer Over Halfspace benchmark 3 (LOH1).
-->
<edge_v>
<mesh>
<periodic>no</periodic>
<write_element_annotations>1</write_element_annotations>
<n_partitions>1</n_partitions>
<files>
<in>gen/gmsh/tet4_200.msh</in>
<out>
<base>gen/edge_v/loh3_200</base>
<extension>.msh</extension>
</out>
</files>
</mesh>
<velocity_model>
<seismic_expression>
<![CDATA[
if( z < 1000 ) {
vp := 4000;
vs := 2000;
qp := 40;
qs := 120;
}
else {
vp := 6000;
vs := 3464;
qp := 69.3;
qs := 155.9;
};
]]>
</seismic_expression>
</velocity_model>
<time>
<n_groups>5</n_groups>
<centroid_based_cfl>true</centroid_based_cfl>
<files>
<out>
<time_steps>gen/edge_v/time_steps_200.csv</time_steps>
</out>
</files>
</time>
</edge_v>
/**
* @file This file is part of EDGE.
*
* @author Alexander Breuer (alex.breuer AT uni-jena.de)
*
* @section LICENSE
* Copyright (c) 2021, Friedrich Schiller University Jena
* Copyright (c) 2017, Regents of the University of California
* All rights reserved.
*
* Redistribution and use in source and binary forms, with or without modification, are permitted provided that the following conditions are met:
*
* 1. Redistributions of source code must retain the above copyright notice, this list of conditions and the following disclaimer.
*
* 2. Redistributions in binary form must reproduce the above copyright notice, this list of conditions and the following disclaimer in the documentation and/or other materials provided with the distribution.
*
* 3. Neither the name of the copyright holder nor the names of its contributors may be used to endorse or promote products derived from this software without specific prior written permission.
*
* THIS SOFTWARE IS PROVIDED BY THE COPYRIGHT HOLDERS AND CONTRIBUTORS "AS IS" AND ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED TO, THE IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR PURPOSE ARE DISCLAIMED. IN NO EVENT SHALL THE COPYRIGHT HOLDER OR CONTRIBUTORS BE LIABLE FOR ANY DIRECT, INDIRECT, INCIDENTAL, SPECIAL, EXEMPLARY, OR CONSEQUENTIAL DAMAGES (INCLUDING, BUT NOT LIMITED TO, PROCUREMENT OF SUBSTITUTE GOODS OR SERVICES; LOSS OF USE, DATA, OR PROFITS; OR BUSINESS INTERRUPTION) HOWEVER CAUSED AND ON ANY THEORY OF LIABILITY, WHETHER IN CONTRACT, STRICT LIABILITY, OR TORT (INCLUDING NEGLIGENCE OR OTHERWISE) ARISING IN ANY WAY OUT OF THE USE OF THIS SOFTWARE, EVEN IF ADVISED OF THE POSSIBILITY OF SUCH DAMAGE.
*
* @section DESCRIPTION
* Tet4 setup of the LOH3 benchmark with a char length of 100 in the region of interest of the layer.
**/
// characteristic length in the
cl1 = 100;
// remaining roi (scaled by 3464/2000)
cl2 = cl1*1.732;
// remaining region (limited to 900m to avoid ill-shaped elements in the layer)
cl3 = (900 < cl2*10) ? 900 : cl2*10;
Include "tet4_base.geo";
/**
* @file This file is part of EDGE.
*
* @author Alexander Breuer (anbreuer AT ucsd.edu)
* @author Alexander Breuer (alex.breuer AT uni-jena.de)
*
* @section LICENSE
* Copyright (c) 2019, Alexander Breuer
* Copyright (c) 2021, Friedrich Schiller University Jena
* Copyright (c) 2017, Regents of the University of California
* All rights reserved.
*
......@@ -29,6 +29,6 @@ cl1 = 200;
cl2 = cl1*1.732;
// remaining region (limited to 900m to avoid ill-shaped elements in the layer)
cl3 = (900 < cl2*3) ? 900 : cl2*3;
cl3 = (900 < cl2*10) ? 900 : cl2*10;
Include "tet4_base.geo";
/**
* @file This file is part of EDGE.
*
* @author Alexander Breuer (anbreuer AT ucsd.edu)
* @author Alexander Breuer (alex.breuer AT uni-jena.de)
*
* @section LICENSE
* Copyright (c) 2019, Alexander Breuer.
* Copyright (c) 2021, Friedrich Schiller University Jena
* Copyright (c) 2016, Regents of the University of California
* All rights reserved.
*
......@@ -19,14 +19,11 @@
* THIS SOFTWARE IS PROVIDED BY THE COPYRIGHT HOLDERS AND CONTRIBUTORS "AS IS" AND ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED TO, THE IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR PURPOSE ARE DISCLAIMED. IN NO EVENT SHALL THE COPYRIGHT HOLDER OR CONTRIBUTORS BE LIABLE FOR ANY DIRECT, INDIRECT, INCIDENTAL, SPECIAL, EXEMPLARY, OR CONSEQUENTIAL DAMAGES (INCLUDING, BUT NOT LIMITED TO, PROCUREMENT OF SUBSTITUTE GOODS OR SERVICES; LOSS OF USE, DATA, OR PROFITS; OR BUSINESS INTERRUPTION) HOWEVER CAUSED AND ON ANY THEORY OF LIABILITY, WHETHER IN CONTRACT, STRICT LIABILITY, OR TORT (INCLUDING NEGLIGENCE OR OTHERWISE) ARISING IN ANY WAY OUT OF THE USE OF THIS SOFTWARE, EVEN IF ADVISED OF THE POSSIBILITY OF SUCH DAMAGE.
*
* @section DESCRIPTION
* Tet4 setup of Layer Over Halfspace (LOH) benchmark.
* Tet4 setup of Layer Over Halfspace (LOH) benchmark 3.