Monday, June 29, 2015

Recent Results

This is a plot of 30 runs for the 4mm cube, poisson dist from 1 to 25um. The variation occurs in the later stages once failure has already set in. 1.528 GPa
2619035[].hhpc-hn1
Bypass inclusion on
4mm/Atest.uda.0032-->62
...

2619260.hhpc-hn1
4mm/Atest.uda.063
4mm with 3 grid
Inclusion calculation off
B4C.xml
dk2:13e12
eta:6e-12
poisson, 1-25^2.6

2619695.hhpc-hn1
4mm/Atest.uda.065
4mm with 3 grid
B4C.xml
dk2:6e-5
eta:832000
poisson, 1-25^2.6
1.45GPa
takes about 15 minutes irl


 Turned out that $RANDOM doesn't work in tcsh, use bash instead.
cat >
cat >>

2619704.hhpc-hn1
4mm/Atest.uda.067
testing the built B4C.xml
1.48GPa!!

Testing the built addric.pbs
2619706.hhpc-hn1
4mm/Atest.uda.068
 1.48GPa- this one isnt right.

2619765.hhpc-hn1

Lets try a new setup that allows us to retain input decks so we know the random seed. We will use jobno as the .xml input into the .ups file.  Lets try it out.

2619765.hhpc-hn1- should make a new $PBS_JOBID.ups and .xml files and run off them...
4mm/Atest.uda.069
1.58GPa

2619766[].hhpc-hn1
4mm/Atest.uda.070-30?
eta=13e12*4e-3^3=832000
-1.5137e+09 GPa
stdev:1.1212e+08
30 runs
lets get more runs on this after we get data for .125mm



2619769[].hhpc-hn1
1mm/Atest.uda.000-99?
Running into a problem with simultaneous calls to sus I think.
Only got to 72 runs.
eta=(13e12*1e-3^3)=13000
-1.4877e+09
9.2777e+07


2619770[].hhpc-hn1
p125mm/Atest.uda.003-102?
did not work, trying again.


eta=13e12*.125e-3^3=25.39

2619843.hhpc-hn1
p125mm/Atest.uda.103
Dropped 25um max to 5um...is it possible the flaw size is too big for the cell size??

This is behavior we get at .p125 mm with 1um to 25um flaw size:
Why does stress critically increase??

Wednesday, June 17, 2015

Post June 17th Update

2617527.hhpc-hn1
HuRamCode2/Atest.uda.003
flawdensity: 13e12
velocity defined
 distribution: delta6um
interval:1e-8
dkappa2:6e-12
chi:0

As an aside: the diluteness can be characterized as eta*s^3 / volume. In a delta distribution of 6um, that comes out to (13e12* 6e-6^3)/(4e-3^3)=43875.  As discussed, we need to drop eta appropriately.

If eta is 13e12 #/m^3, then for a 4mm edge length, #=832000.

832000*6e-6^3/4e-3^3=.0028, which is properly dilute.  We will need to do this going forward.


poisson 1um to 25um exp2.6
 possible sus or pdTdt error.
2617783.hhpc-hn1

recopied sus, same error
fixed pdTdt, recompiling.
2617784.hhpc-hn1

 visit

2.7.3 doesnt work
2.6.0 seems to work
/data1/hhpc-hn8/Programs/visit/visit_install
hhpc-hn8.idies.jhu.edu

Looks like particle id 281483566776320 is the central particle.

Questions for Ravi- why do we need ghost cells?
Why does ghost cells generate 27 particles? I guess 26 are ghost particles....check this with ravi.

Looks like a non-uniaxial stress state. will modify to drop other stresses.

2618383.hhpc-hn1
def grad: .88 1 1
looks like it went the wrong direction. lets try 1.1
2618389.hhpc-hn1
def grad: .88 1.1 1.1
Atest.uda.009
around 100th timestep- looks like too much tension on the non central axis:




Atest.uda.010
def grad .88 1.05 1.05
2618401.hhpc-hn1
B4C.xml
dk2:6e-6
2618413.hhpc-hn1
 Atest.uda.012
 def grad .88 1.05 1.05
B4C.xml
dk2:6e-12

328th timestep, doesnt make curve better.  1.57 GPa. I thnk 1.04 is best.


Atest.uda.013
2618479.hhpc-hn1
def grad .88 1.04 1.04
B4C.xml
dk2:6e-12
eta:10e12

comes to 2.38e9 GPa

3e12 = .81 GPa... 2.38-3.75=1.37/.81--> subtract 5e12.

Atest.uda.014
def grad .88 1.04 1.04
B4C.xml
dk2:6e-12
eta:5e12
2618559.hhpc-hn1
2.42 GPa


Atest.uda.016
def grad .88 1.04 1.04
B4C.xml
dk2:6e-10
eta:5e10
2618563.hhpc-hn1

2.43 seems to be the saturated strength.

Lets try a run with the expected flaw density for a 4mm cube.
2618569.hhpc-hn1
Atest.uda.017
def grad .88 1.04 1.04
B4C.xml
dk2:6e-5
eta:832000

2.52e9, beautiful curve.

Atest.uda.018
def grad .88 1.04 1.04
B4C.xml
dk2:6e-5
eta:832000
poisson, 1-25^2.6
2618581.hhpc-hn1

1.85 GPa!!!
Go back to 6um delta.  try different volume sizes.
Atest.uda.019
def grad .88 1.04 1.04
.4mm edge length
B4C.xml
dk2:6e-5
eta:832000
delta,6um
2618647.hhpc-hn1
2.52...... wtf


2618652.hhpc-hn1
Atest.uda.020
def grad .88 1.04 1.04
.4mm edge length
B4C.xml
dk2:6e-2
eta:832
delta,6um

3.54


2618678.hhpc-hn1
2618679.hhpc-hn1
Atest.uda.021
def grad .88 1.04 1.04
1mm edge length
B4C.xml
dk2:6e-3
eta:832
delta,6um

3.54GPa again....makes sense, didnt adapt eta.

 2618681.hhpc-hn1
Atest.uda.022
def grad .88 1.04 1.04
1mm edge length
B4C.xml
dk2:6e-3
eta:13000
delta,6um


Array with poisson:
2618682[].hhpc-h - didnt run far enough, i think im under an order of magnitude
2618687[].hhpc-h  -1.52, -1.52, -1.52....

2618695.hhpc-hn1
Atest.uda.000
1mm, new grid.
B4C.xml
dk2:6e-3
eta:13000
poisson, 1-25^2.6

1.53


2618696.hhpc-hn1
4mm/Atest.uda.000
4mm, new grid.
B4C.xml
dk2:6e-5
eta:832000
poisson, 1-25^2.6
1.53

2618697.hhpc-hn1
p125mm/Atest.uda.000
.125mm with 3 grid
B4C.xml
dk2: 25.39
eta:6e0
poisson, 1-25^2.6

1.53??
curve wasn't as pretty, no dropoff.

2618698.hhpc-hn1
p125mm/Atest.uda.001
.125mm with 3 grid
B4C.xml
dk2: 25.39
eta:6e-1
poisson, 1-25^2.6