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