Shortcuts nix

Alegra: /usr/cta/restricted/MHD

From JHU:
Have to use ossh
Have to use wired ethernet for ./pkinit dmallic1@HPCMP.HPC.MIL


SSH:
ssh -X -Y dmallic1@hhpc-hn8.idies.jhu.edu
ssh -X -Y dmallic1@hhpc-hn1.pha.jhu.edu
Moving files
scp works great for hn8

Emacs:
ctrl x ctrl s - save
ctrl x ctrl c- quit
ctrl v- page down

visit
/data1/hhpc-hn8/Programs/visit 2.6.0

find . -name "foo*"

Aliases/Path:
emacs -nw ~/.bash_profile

alias xyz='cd /xyz'
PATH=$xyz/

source ~/.bash_profile

Uintah-
To install:
config log
make cleanreally
 make -j 6
After making changes:
/data1/hhpc-hn8/dmallic1/JHUUintah_SVN_ravi/src/CCA/Components/MPM/ConstitutiveModel
move HuRam.cc and .h, then make cleanreally then make -j 6IN THE UINTAH SVN DIR

and constitutive model.xml into standalone/inputs/UPS_SPEC
 HuRamesh.cc
remove pdTdt "_preReloc"


Test

4mm.uda.010-My input deck with Ravi's dimensions gives same answer as Ravi's input deck with my dimensions. This indicates that its not the input deck thats causing the problem.

Ravi suggests increasing flaw density and checking dkappa2 parameter on Jun3

2613095.hhpc-hn1 4mm.uda.008?

2613096.hhpc-hn1? 4mm.uda.009

 

This run isn't complete yet but i suspect the same outcome.
2613164.hhpc-hn1  4mm.uda.010? i think this was dkappa e-6

New calculation:
for 4mm^3 volume, at 22e12 #/m^3, 4mm^3 = 4e-9 m^3, so #=88000  This makes eta way too low to make sense.

2613218.hhpc-hn1
4mm.uda.011
flawdensity:13e14
distribution:delta6um
interval:1e-8
dkappa2:6e-14
chi:2

What is that sub peak? I thought moving dkappa2 would fix it?
Max peak occurs at about 40um.... this is 15 um past what experiments say it should be. Perhaps increase flaw density?
3.75e9 strength!

2613226.hhpc-hn1
4mm.uda.012
flawdensity:13e14
distribution:delta6um
interval:1e-8
dkappa2:6e-6
chi:2

Increasing dkappa clearly makes the curve cleaner.
3.83GPa
SR: 1242


2613228.hhpc-hn1
4mmHigherEta/4mm.uda.000
flawdensity:13e16
distribution:delta6um
interval:1e-8
dkappa2:6e-16
chi:2

Increasing the flaw density to 13e16 puts us far past a normal strength for the material. It looks like the inflection point at about 40us was missed, so the actual time to failure may not have changed.

2613296.hhpc-hn1

4mmHigherEta/4mm.uda.001
flawdensity:13e15
distribution:delta6um
interval:1e-8
dkappa2:6e-15
chi:2

2613297.hhpc-hn1
HuRamCode2/p125mm.uda.001
flawdensity:13e14
distribution:delta6um
interval:1e-8
dkappa2:6e-6
chi:2

5.5Gpa at smallest size (.125mm^3) at 1101 1e3 SR

2613314.hhpc-hn1
HuRamCode2/p125mm.uda.002
flawdensity:13e14
velocity .18
distribution:delta6um
interval:1e-8
dkappa2:6e-6
chi:2
check this at correct SR- came out to 1029 from 0 to 4300 step.

3.51GPa

Try .19 next

2613324.hhpc-hn1
HuRamCode2/p125mm.uda.003
flawdensity:13e14
velocity .19
distribution:delta6um
interval:1e-8
dkappa2:6e-6
chi:2
CANCELLED THIS, ABOVE SR IS FINE AT 1000



2613316.hhpc-hn1
HuRamCode2/p25mm.uda.001
flawdensity:13e14
velocity .4
distribution:delta6um
interval:1e-8
dkappa2:6e-6
chi:2

SR Came to 1138
3.69GPa

2613334.hhpc-hn1
HuRamCode2/p25mm.uda.002
flawdensity:13e14
velocity .35
distribution:delta6um
interval:1e-8
dkappa2:6e-6
chi:2

STR; 3.69GPa
SR 994.62

2613317.hhpc-hn1
HuRamCode2/p5mm.uda.001
flawdensity:13e14
velocity .8
distribution:delta6um
interval:1e-8
dkappa2:6e-6
chi:2

SR:1141-drop to .7
3.51 GPa

2613332.hhpc-hn1
HuRamCode2/p5mm.uda.002
flawdensity:13e14
velocity .7
distribution:delta6um
interval:1e-8
dkappa2:6e-6
chi:2

This result is not converging!!!!

I looked at dStress:

Looks like the two peaks are there but simply hidden.  They occur at ~6050 and ~6700.

SR is 964 through 6000
4.4GPa

2613318.hhpc-hn1
HuRamCode2/1mm.uda.001
flawdensity:13e14
velocity 1.6
distribution:delta6um
interval:1e-8
dkappa2:6e-6
chi:2

SR:1141- lets drop to 1.4
3.39GPa
Probably need to drop dkappa2 or change chi

2613330.hhpc-hn1
HuRamCode2/1mm.uda.002
flawdensity:13e14
velocity 1.4
distribution:delta6um
interval:1e-8
dkappa2:6e-6
chi:2

4.54GPa
984 SR

2613319.hhpc-hn1
HuRamCode2/2mm.uda.001
flawdensity:13e14
velocity 3
distribution:delta6um
interval:1e-8
dkappa2:6e-6
chi:2

SR:1073- a little high.  Should decrease velocity to 2.85
3.54 GPa?? Higher SR explains this

2613325.hhpc-hn1
HuRamCode2/2mm.uda.002
flawdensity:13e14
velocity 2.825
distribution:delta6um
interval:1e-8
dkappa2:6e-6
chi:2
SR: 1005
3.68GPa

2613321.hhpc-hn1
4mmHigherEta/4mm.uda.002
flawdensity:13e14
velocity 7
distribution:poisson 4um to 8um exp3
interval:1e-8
dkappa2:6e-6
chi:2

what the hell is going on with the poisson distribution??? gonna zero out variability and hit it with 6um to 6um.

2613350.hhpc-hn1
4mmHigherEta/4mm.uda.003
flawdensity:13e14
velocity 5.6
distribution:poisson 6um to 7um exp3
interval:1e-8
dkappa2:6e-6
chi:2


2613327.hhpc-hn1
HuRamCode2/4mm.uda.013
flawdensity:13e14
velocity 5.75
distribution:delta6um
interval:1e-8
dkappa2:6e-6
chi:2
3.61GPa
SR: 1025

Final Results for delta 6 distribution
4mm:3.61
2mm:3.68
1mm:4.54
.5mm: 4.4
.25mm: 3.69
.125mm: 3.51

2613498.hhpc-hn1
4mmHigherEta/4mm.uda.004
flawdensity:13e14
velocity 5.6
distribution:poisson 6um to 7um exp3
interval:1e-8
dkappa2:6e-14
chi:2

2613499.hhpc-hn1
HuRamCode2/4mm.uda.014
flawdensity:13e13
velocity 5.6
distribution:poisson 6um to 7um exp3
interval:1e-8
dkappa2:6e-6
chi:2

2613504.hhpc-hn1
4mmHigherEta/4mm.uda.005
flawdensity:13e12
velocity 5.6
distribution:poisson 6um to 7um exp3
interval:1e-8
dkappa2:6e-6
chi:2

We should be able to follow one of ravi's poisson runs to get results. Lets do that.
The run with zeros changed nothing.

2613912.hhpc-hn1
HuRamCode2/4mm.uda.015
flawdensity:13e12
velocity 5.75
distribution:poisson 1um to 7um exp3
interval:1e-8
dkappa2:6e-12
chi:2
Bottom of first peak: 1.2e-1

2613957.hhpc-hn1
HuRamCode2/4mm.uda.016
flawdensity:13e12
velocity 5.75
distribution:poisson 1um to 7um exp3
interval:1e-8
dkappa2:6e-12
chi:0
Looks like it went earlier.... not what u want

2613958.hhpc-hn1

4mmHigherEta/4mm.uda.006
flawdensity:13e11
velocity 5.75
distribution:poisson 1um to 7um exp3
interval:1e-8
dkappa2:6e-11
chi:2
looks a little higher....

2613960.hhpc-hn1
HuRamCode2/4mm.uda.017
flawdensity:13e10
velocity 5.75
distribution:poisson 1um to 7um exp3
interval:1e-8
dkappa2:6e-10
chi:0
e4.

2614091.hhpc-hn1
 HuRamCode2/4mm.uda.018
flawdensity:13e9
velocity 5.75
distribution:poisson 1um to 7um exp3
interval:1e-8
dkappa2:6e-9
chi:0

Trying addric.ups
2614756.hhpc-hn1
flawdensity:22e12
velocity time def grad
distribution:poisson 1um to 25um exp2.6
interval:1e-8
dkappa2:6e-12
chi:0
will need a try with chi if this doesnt work.
2614842.hhpc-hn1 again
HuRamCode2/Atest.uda.001
2615163.hhpc-hn1 again
HuRamCode2/Atest.uda.001
Is it possible that the boundary conditions are causing non uniaxial stress??
Make the box on both sides.


2615370.hhpc-hn1
HuRamCode2/4mm.uda.019 possibly 20 its 19
flawdensity:22e12
velocity 5.75
distribution:poisson 1um to 25um exp2.6
interval:1e-8
dkappa2:6e-12
chi:0

2615371.hhpc-hn1
HuRamCode2/Atest.uda.002
flawdensity:22e12
velocity defined
distribution:poisson 1um to 25um exp2.6
interval:1e-8
dkappa2:6e-12
chi:0

=.5*(Y2^2+AB2^2+AE2^2-1)

Code updated on June 17th