Commit 6283226d authored by Alex Breuer's avatar Alex Breuer

Added LOH.3 input.

parent a54b7b89
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.
* Contact: Alexander Breuer
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) 2019, Alexander Breuer
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 for 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>
<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.4E9</lambda>
<mu>32.4E9</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_200</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.01</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_200/wf</file>
<int>1.0</int>
</wave_field>
</output>
</cfr>
</edge>
[![CC0](http://i.creativecommons.org/p/zero/1.0/88x31.png)](http://creativecommons.org/publicdomain/zero/1.0/)
To the extent possible under law, [Alexander Nikolas Breuer](http://dial3343.org) has waived all copyright and related or neighboring rights to the following work:
* `tet4_200.msh.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
/**
* @file This file is part of EDGE.
*
* @author Alexander Breuer (anbreuer AT ucsd.edu)
*
* @section LICENSE
* Copyright (c) 2019, Alexander Breuer
* 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 200 in the region of interest of the layer.
**/
// characteristic length in the
cl1 = 200;
// 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*3) ? 900 : cl2*3;
Include "tet4_base.geo";
This source diff could not be displayed because it is stored in LFS. You can view the blob instead.
/**
* @file This file is part of EDGE.
*
* @author Alexander Breuer (anbreuer AT ucsd.edu)
*
* @section LICENSE
* Copyright (c) 2019, Alexander Breuer.
* Copyright (c) 2016, 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 Layer Over Halfspace (LOH) benchmark.
**/
Mesh.Optimize=1;
/*
* Define Geometry.
*/
// construct surrounding free-surface
Point(1) = { -26000, -26000, 0 };
Point(2) = { 32000, -26000, 0 };
Point(3) = { 32000, 32000, 0 };
Point(4) = { -26000, 32000, 0 };
Line(1) = {1, 2};
Line(2) = {2, 3};
Line(3) = {3, 4};
Line(4) = {4, 1};
// top surface
Line Loop(101) = { 1, 2, 3, 4 };
Plane Surface(201) = {101};
// extrude in depth to get the layer
extList[] = Extrude {0, 0, 1000} { Surface{201}; };
// extrude in depth to get the half-space
extList[] = Extrude {0, 0, 32000} { Surface{extList}; };
/*
* set characteristic lengths
*/
// high-res layer
Field[1] = Box;
Field[1].VIn = cl1;
Field[1].VOut = cl3;
Field[1].XMin = -2000;
Field[1].XMax = 10000;
Field[1].YMin = -2000;
Field[1].YMax = 12000;
Field[1].ZMin = 0;
Field[1].ZMax = 1000;
// high-res half-space
Field[2] = Box;
Field[2].VIn = cl2;
Field[2].VOut = cl3;
Field[2].XMin = -2000;
Field[2].XMax = 10000;
Field[2].YMin = -2000;
Field[2].YMax = 12000;
Field[2].ZMin = 1000;
Field[2].ZMax = 7000;
Field[3] = Min;
Field[3].FieldsList = {1,2};
Background Field = 3;
/*
* Set physical values
*/
// set volumes
Physical Volume(1) = {1};
Physical Volume(2) = {2};
// set free surface boundary
Physical Surface(101) = {201};
// set outflow boundary
Physical Surface(105) = {210, 214, 218, 222, 236, 240, 244, 232, 245};
[![CC0](http://i.creativecommons.org/p/zero/1.0/88x31.png)](http://creativecommons.org/publicdomain/zero/1.0/)
To the extent possible under law, [Alexander Nikolas Breuer](http://dial3343.org) has waived all copyright and related or neighboring rights to the following work:
* `loh3_uns_src.h5`
This work is published from: Germany.
Source for the LOH.3 Benchmark
==============================
This directory contains the point source for the LOH.3 benchmark.
It was converted from the respective LOH.1 source to HDF5.
EDGE-version for the script: `fd980ebfb72f60d86ce503dd47185e85eaa14837`.
Contact: Alex Breuer
Used Commands
-------------
.. code-block: bash
wget https://opt.dial3343.org/3343/edge_opt/raw/master/bench/seismic/wp/loh1/sources/loh1_uns_src.nc
nrfToHdf5.py -l 32.4E9 -m 32.4E9 -i loh1_uns_src.nc -o loh3_uns_src.h5
This source diff could not be displayed because it is stored in LFS. You can view the blob instead.
Markdown is supported
0% or
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment