Small Scale Heterogeneities Validation: Difference between revisions

From SCECpedia
Jump to navigationJump to search
← Older edit
Davidgil (talk | contribs)
1,669 edits
Davidgil (talk | contribs)
1,669 edits
No edit summary
 
(10 intermediate revisions by the same user not shown)
Line 3: Line 3:
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.
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.


=== Mesh 0.0 Hurst ===
Common parameters:


Specifications:
Hurst = 0.0<br />
L1 = 200<br />
d = 20<br />
d = 20<br />
st23 = 5<br />
st23 = 5<br />
Line 14: Line 10:
n2 = 1000<br />
n2 = 1000<br />
n3 = 1500<br />
n3 = 1500<br />
== 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 ====
==== Plots at surface level ====


{|
{|
| '''Matlab Code'''
| '''Matlab Code 0.0 Hurst'''
| '''C SSH Code'''
| '''C SSH Code 0.0 Hurst'''
| '''Difference 0.0 Hurst'''
|-
|-
| [[File:matlab_H0_0.png|600px]]
| [[File:matlab_H0_0.png|600px]]
| [[File:myssh_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'''
|-
|-
| '''Difference'''
| [[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:compare_H0_0.png|600px]]
| [[File:matlab_H4_0.png|600px]]
|
| [[File:myssh_H4_0.png|600px]]
| [[File:compare_H4_0.png|600px]]
|}
|}


Line 34: Line 50:


{|
{|
| '''Matlab Code'''
| '''Matlab Code 0.0 Hurst'''
| '''C SSH Code'''
| '''C SSH Code 0.0 Hurst'''
| '''Difference 0.0 Hurst'''
|-
|-
| [[File:matlab_H0_250.png|600px]]
| [[File:matlab_H0_250.png|600px]]
| [[File:myssh_H0_250.png|600px]]
| [[File:myssh_H0_250.png|600px]]
| [[File:compare_H0_250.png|600px]]
|}
{|
| '''Matlab Code 0.2 Hurst'''
| '''C SSH Code 0.2 Hurst'''
| '''Difference 0.2 Hurst'''
|-
|-
| '''Difference'''
| [[File:matlab_H2_250.png|600px]]
|
| [[File:myssh_H2_250.png|600px]]
| [[File:compare_H2_250.png|600px]]
|}
 
{|
| '''Matlab Code 0.4 Hurst'''
| '''C SSH Code 0.4 Hurst'''
| '''Difference 0.4 Hurst'''
|-
|-
| [[File:compare_H0_250.png|600px]]
| [[File:matlab_H4_250.png|600px]]
|
| [[File:myssh_H4_250.png|600px]]
| [[File:compare_H4_250.png|600px]]
|}
|}


Line 50: Line 82:


{|
{|
| '''Matlab Code'''
| '''Matlab Code 0.0 Hurst'''
| '''C SSH Code'''
| '''C SSH Code 0.0 Hurst'''
| '''Difference 0.0 Hurst'''
|-
|-
| [[File:matlab_H0_500.png|600px]]
| [[File:matlab_H0_500.png|600px]]
| [[File:myssh_H0_500.png|600px]]
| [[File:myssh_H0_500.png|600px]]
| [[File:compare_H0_500.png|600px]]
|}
{|
| '''Matlab Code 0.2 Hurst'''
| '''C SSH Code 0.2 Hurst'''
| '''Difference 0.2 Hurst'''
|-
|-
| '''Difference'''
| [[File:matlab_H2_500.png|600px]]
|
| [[File:myssh_H2_500.png|600px]]
| [[File:compare_H2_500.png|600px]]
|}
 
{|
| '''Matlab Code 0.4 Hurst'''
| '''C SSH Code 0.4 Hurst'''
| '''Difference 0.4 Hurst'''
|-
|-
| [[File:compare_H0_499.png|600px]]
| [[File:matlab_H4_500.png|600px]]
|
| [[File:myssh_H4_500.png|600px]]
| [[File:compare_H4_500.png|600px]]
|}
|}


=== Mesh 0.2 Hurst ===
== Validation by Correlation Length ==


Specifications:
In this condition, we hold hurst to be a constant 0.2, but vary correlation length to be either 100m, 200m, or 300m.


Hurst = 0.2<br />
==== Plots at surface level ====
L1 = 200<br />
d = 20<br />
st23 = 5<br />
n1 = 500<br />
n2 = 1000<br />
n3 = 1500<br />


==== Plots at surface level ====
{|
| '''Matlab Code 100m L1'''
| '''C SSH Code 100m L1'''
| '''Difference 100m L1'''
|-
| [[File:matlab_H0_0.png|600px]]
| [[File:myssh_H0_0.png|600px]]
| [[File:compare_H0_0.png|600px]]
|}


{|
{|
| '''Matlab Code'''
| '''Matlab Code 200m L1'''
| '''C SSH Code'''
| '''C SSH Code 200m L1'''
| '''Difference 200m L1'''
|-
|-
| [[File:matlab_H2_0.png|600px]]
| [[File:matlab_H2_0.png|600px]]
| [[File:myssh_H2_0.png|600px]]
| [[File:myssh_H2_0.png|600px]]
| [[File:compare_H2_0.png|600px]]
|}
{|
| '''Matlab Code 300m L1'''
| '''C SSH Code 300m L1'''
| '''Difference 300m L1'''
|-
|-
| '''Difference'''
| [[File:matlab_L300_0.png|600px]]
|
| [[File:myssh_L300_0.png|600px]]
|-
| [[File:compare_L300_0.png|600px]]
| [[File:compare_H2_0.png|600px]]
|
|}
|}


Line 94: Line 150:


{|
{|
| '''Matlab Code'''
| '''Matlab Code 100m L1'''
| '''C SSH Code'''
| '''C SSH Code 100m L1'''
| '''Difference 100m L1'''
|-
| [[File:matlab_L100_250.png|600px]]
| [[File:myssh_L100_250.png|600px]]
| [[File:compare_L100_250.png|600px]]
|}
 
{|
| '''Matlab Code 200m L1'''
| '''C SSH Code 200m L1'''
| '''Difference 200m L1'''
|-
|-
| [[File:matlab_H2_250.png|600px]]
| [[File:matlab_H2_250.png|600px]]
| [[File:myssh_H2_250.png|600px]]
| [[File:myssh_H2_250.png|600px]]
|-
| '''Difference'''
|
|-
| [[File:compare_H2_250.png|600px]]
| [[File:compare_H2_250.png|600px]]
|
|}
|}
==== Plots for bottom slice ====


{|
{|
| '''Matlab Code'''
| '''Matlab Code 300m L1'''
| '''C SSH Code'''
| '''C SSH Code 300m L1'''
| '''Difference 300m L1'''
|-
|-
| [[File:matlab_H2_500.png|600px]]
| [[File:matlab_L300_250.png|600px]]
| [[File:myssh_H2_500.png|600px]]
| [[File:myssh_L300_250.png|600px]]
|-
| [[File:compare_L300_250.png|600px]]
| '''Difference'''
|
|-
| [[File:compare_H2_500.png|600px]]
|
|}
|}


=== Mesh 0.4 Hurst ===
==== Plots for bottom slice ====
 
Specifications:
 
Hurst = 0.4<br />
L1 = 200<br />
d = 20<br />
st23 = 5<br />
n1 = 500<br />
n2 = 1000<br />
n3 = 1500<br />
 
==== Plots at surface level ====


{|
{|
| '''Matlab Code'''
| '''Matlab Code 100m L1'''
| '''C SSH Code'''
| '''C SSH Code 100m L1'''
| '''Difference 100m L1'''
|-
|-
| [[File:matlab_H4_0.png|600px]]
| [[File:matlab_L100_500.png|600px]]
| [[File:myssh_H4_0.png|600px]]
| [[File:myssh_L100_500.png|600px]]
|-
| [[File:compare_L100_500.png|600px]]
| '''Difference'''
|
|-
| [[File:compare_H4_0.png|600px]]
|
|}
|}
==== Plots at 50% depth ====


{|
{|
| '''Matlab Code'''
| '''Matlab Code 200m L1'''
| '''C SSH Code'''
| '''C SSH Code 200m L1'''
| '''Difference 200m L1'''
|-
|-
| [[File:matlab_H4_250.png|600px]]
| [[File:matlab_H2_500.png|600px]]
| [[File:myssh_H4_250.png|600px]]
| [[File:myssh_H2_500.png|600px]]
|-
| [[File:compare_H2_500.png|600px]]
| '''Difference'''
|
|-
| [[File:compare_H4_250.png|600px]]
|
|}
|}
==== Plots for bottom slice ====


{|
{|
| '''Matlab Code'''
| '''Matlab Code 300m L1'''
| '''C SSH Code'''
| '''C SSH Code 300m L1'''
|-
| '''Difference 300m L1'''
| [[File:matlab_H4_500.png|600px]]
| [[File:myssh_H4_500.png|600px]]
|-
|-
| '''Difference'''
| [[File:matlab_L300_500.png|600px]]
|
| [[File:myssh_L300_500.png|600px]]
|-
| [[File:compare_L300_500.png|600px]]
| [[File:compare_H4_500.png|600px]]
|
|}
|}

Latest revision as of 23:40, 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 Code 0.2 Hurst C SSH Code 0.2 Hurst Difference 0.2 Hurst
Matlab Code 0.4 Hurst C SSH Code 0.4 Hurst Difference 0.4 Hurst

Plots at 50% depth

Matlab Code 0.0 Hurst C SSH Code 0.0 Hurst Difference 0.0 Hurst
Matlab Code 0.2 Hurst C SSH Code 0.2 Hurst Difference 0.2 Hurst
Matlab Code 0.4 Hurst C SSH Code 0.4 Hurst Difference 0.4 Hurst

Plots for bottom slice

Matlab Code 0.0 Hurst C SSH Code 0.0 Hurst Difference 0.0 Hurst
Matlab Code 0.2 Hurst C SSH Code 0.2 Hurst Difference 0.2 Hurst
Matlab Code 0.4 Hurst C SSH Code 0.4 Hurst Difference 0.4 Hurst

Validation by Correlation Length

In this condition, we hold hurst to be a constant 0.2, but vary correlation length to be either 100m, 200m, or 300m.

Plots at surface level

Matlab Code 100m L1 C SSH Code 100m L1 Difference 100m L1
Matlab Code 200m L1 C SSH Code 200m L1 Difference 200m L1
Matlab Code 300m L1 C SSH Code 300m L1 Difference 300m L1

Plots at 50% depth

Matlab Code 100m L1 C SSH Code 100m L1 Difference 100m L1
Matlab Code 200m L1 C SSH Code 200m L1 Difference 200m L1
Matlab Code 300m L1 C SSH Code 300m L1 Difference 300m L1

Plots for bottom slice

Matlab Code 100m L1 C SSH Code 100m L1 Difference 100m L1
Matlab Code 200m L1 C SSH Code 200m L1 Difference 200m L1
Matlab Code 300m L1 C SSH Code 300m L1 Difference 300m L1
Retrieved from "https://strike.scec.org/index.php?title=Small_Scale_Heterogeneities_Validation&oldid=9799"