Difference between revisions of "Small Scale Heterogeneities Validation"

From SCECpedia
Jump to navigationJump to search
Line 81: Line 81:
 
==== Plots for bottom slice ====
 
==== Plots for bottom slice ====
  
{|
 
| '''Matlab Code 0.0 Hurst'''
 
| '''C SSH Code 0.0 Hurst'''
 
| '''Difference 0.0 Hurst'''
 
|-
 
| [[File:matlab_H0_0.png|600px]]
 
| [[File:myssh_H0_0.png|600px]]
 
| [[File:compare_H0_0.png|600px]]
 
|}
 
 
{|
 
| '''Matlab Code 0.2 Hurst'''
 
| '''C SSH Code 0.2 Hurst'''
 
| '''Difference 0.2 Hurst'''
 
|-
 
| [[File:matlab_H2_0.png|600px]]
 
| [[File:myssh_H2_0.png|600px]]
 
| [[File:compare_H2_0.png|600px]]
 
|}
 
 
{|
 
| '''Matlab Code 0.4 Hurst'''
 
| '''C SSH Code 0.4 Hurst'''
 
| '''Difference 0.4 Hurst'''
 
|-
 
| [[File:matlab_H4_0.png|600px]]
 
| [[File:myssh_H4_0.png|600px]]
 
| [[File:compare_H4_0.png|600px]]
 
|}
 
{|
 
| '''Matlab Code 0.0 Hurst'''
 
| '''C SSH Code 0.0 Hurst'''
 
| '''Difference 0.0 Hurst'''
 
|-
 
| [[File:matlab_H0_0.png|600px]]
 
| [[File:myssh_H0_0.png|600px]]
 
| [[File:compare_H0_0.png|600px]]
 
|}
 
 
{|
 
| '''Matlab Code 0.2 Hurst'''
 
| '''C SSH Code 0.2 Hurst'''
 
| '''Difference 0.2 Hurst'''
 
|-
 
| [[File:matlab_H2_0.png|600px]]
 
| [[File:myssh_H2_0.png|600px]]
 
| [[File:compare_H2_0.png|600px]]
 
|}
 
 
{|
 
| '''Matlab Code 0.4 Hurst'''
 
| '''C SSH Code 0.4 Hurst'''
 
| '''Difference 0.4 Hurst'''
 
|-
 
| [[File:matlab_H4_0.png|600px]]
 
| [[File:myssh_H4_0.png|600px]]
 
| [[File:compare_H4_0.png|600px]]
 
|}
 
 
{|
 
{|
 
| '''Matlab Code 0.0 Hurst'''
 
| '''Matlab Code 0.0 Hurst'''

Revision as of 23:35, 5 March 2014

Introduction

We would like to validate that the media with the small-scale heterogeneities that are being generated with the new C code are correct. One way to do this would be to take a look at some plots to see if the generated heterogeneities look correct. We present a few different meshes and plots showing the heterogeneities at various depths, as well as difference plots.

Common parameters:

d = 20
st23 = 5
n1 = 500
n2 = 1000
n3 = 1500

Validation by Hurst Value

In this test, we vary hurst to be either 0.0, 0.2, or 0.4. L1 is held constant at 200m.

Plots at surface level

Matlab Code 0.0 Hurst C SSH Code 0.0 Hurst Difference 0.0 Hurst
Matlab H0 0.png Myssh H0 0.png Compare H0 0.png
Matlab Code 0.2 Hurst C SSH Code 0.2 Hurst Difference 0.2 Hurst
Matlab H2 0.png Myssh H2 0.png Compare H2 0.png
Matlab Code 0.4 Hurst C SSH Code 0.4 Hurst Difference 0.4 Hurst
Matlab H4 0.png Myssh H4 0.png Compare H4 0.png

Plots at 50% depth

Matlab Code 0.0 Hurst C SSH Code 0.0 Hurst Difference 0.0 Hurst
Matlab H0 250.png Myssh H0 250.png Compare H0 250.png
Matlab Code 0.2 Hurst C SSH Code 0.2 Hurst Difference 0.2 Hurst
Matlab H2 250.png Myssh H2 250.png Compare H2 250.png
Matlab Code 0.4 Hurst C SSH Code 0.4 Hurst Difference 0.4 Hurst
Matlab H4 250.png Myssh H4 250.png Compare H4 250.png

Plots for bottom slice

Matlab Code 0.0 Hurst C SSH Code 0.0 Hurst Difference 0.0 Hurst
Matlab H0 500.png Myssh H0 500.png Compare H0 500.png
Matlab Code 0.2 Hurst C SSH Code 0.2 Hurst Difference 0.2 Hurst
Matlab H2 500.png Myssh H2 500.png Compare H2 500.png
Matlab Code 0.4 Hurst C SSH Code 0.4 Hurst Difference 0.4 Hurst
Matlab H4 500.png Myssh H4 500.png Compare H4 500.png
Matlab Code C SSH Code
Matlab H0 500.png Myssh H0 500.png
Difference
Compare H0 499.png

Mesh 0.2 Hurst

Specifications:

Hurst = 0.2
L1 = 200
d = 20
st23 = 5
n1 = 500
n2 = 1000
n3 = 1500

Plots at surface level

Matlab Code C SSH Code
Matlab H2 0.png Myssh H2 0.png
Difference
Compare H2 0.png

Plots at 50% depth

Matlab Code C SSH Code
Matlab H2 250.png Myssh H2 250.png
Difference
Compare H2 250.png

Plots for bottom slice

Matlab Code C SSH Code
Matlab H2 500.png Myssh H2 500.png
Difference
Compare H2 500.png

Mesh 0.4 Hurst

Specifications:

Hurst = 0.4
L1 = 200
d = 20
st23 = 5
n1 = 500
n2 = 1000
n3 = 1500

Plots at surface level

Matlab Code C SSH Code
Matlab H4 0.png Myssh H4 0.png
Difference
Compare H4 0.png

Plots at 50% depth

Matlab Code C SSH Code
Matlab H4 250.png Myssh H4 250.png
Difference
Compare H4 250.png

Plots for bottom slice

Matlab Code C SSH Code
Matlab H4 500.png Myssh H4 500.png
Difference
Compare H4 500.png

Validation by Correlation Length

100m Correlation Length

Specifications:

Hurst = 0.2
L1 = 100
d = 20
st23 = 5
n1 = 500
n2 = 1000
n3 = 1500

Plots at surface level

Matlab Code C SSH Code
Matlab L100 0.png Myssh L100 0.png
Difference
Compare L100 0.png

Plots at 50% depth

Matlab Code C SSH Code
Matlab L100 250.png Myssh L100 250.png
Difference
Compare L100 250.png

Plots for bottom slice

Matlab Code C SSH Code
Matlab L100 500.png Myssh L100 500.png
Difference
Compare L100 500.png

200m Correlation Length

Specifications:

Hurst = 0.2
L1 = 200
d = 20
st23 = 5
n1 = 500
n2 = 1000
n3 = 1500

Plots at surface level

Matlab Code C SSH Code
Matlab H2 0.png Myssh H2 0.png
Difference
Compare H2 0.png

Plots at 50% depth

Matlab Code C SSH Code
Matlab H2 250.png Myssh H2 250.png
Difference
Compare H2 250.png

Plots for bottom slice

Matlab Code C SSH Code
Matlab H2 500.png Myssh H2 500.png
Difference
Compare H2 500.png

300m Correlation Length

Specifications:

Hurst = 0.2
L1 = 300
d = 20
st23 = 5
n1 = 500
n2 = 1000
n3 = 1500

Plots at surface level

Matlab Code C SSH Code
Matlab L300 0.png Myssh L300 0.png
Difference
Compare L300 0.png

Plots at 50% depth

Matlab Code C SSH Code
Matlab L300 250.png Myssh L300 250.png
Difference
Compare L300 250.png

Plots for bottom slice

Matlab Code C SSH Code
Matlab L300 500.png Myssh L300 500.png
Difference
Compare L300 500.png