www.geongrid.org CYBERINFRASTRUCTURE FOR THE GEOSCIENCES A community modeling A community modeling environment: geodynamic environment: geodynamic integration of multi-scale integration of multi-scale geoscience data geoscience data Mian Liu Mian Liu 1 , Huai Zhang , Huai Zhang 1,2 1,2 , , Youqing Yang Youqing Yang 1 , Qingsong Li , Qingsong Li 1 , , Yaolin Shi Yaolin Shi 2 1-University of Missouri-Columbia 2-Computational Geodynamic Lab, CAS, China
35
Embed
CYBERINFRASTRUCTURE FOR THE GEOSCIENCES A community modeling environment: geodynamic integration of multi-scale geoscience data Mian Liu.
This document is posted to help you gain knowledge. Please leave a comment to let me know what you think about it! Share it to your friends and learn new things together.
Transcript
www.geongrid.orgCYBERINFRASTRUCTURE FOR THE GEOSCIENCES
A community modeling environment: A community modeling environment: geodynamic integration of multi-scale geodynamic integration of multi-scale
DO J=1,NMATEPRMT(J) = EMATE((IMATE-1)*NMATE+J)End doPRMT(NMATE+1)=TIMEPRMT(NMATE+2)=DTprmt(nmate+3)=imateprmt(nmate+4)=num
Other element matrix computing SubsPDE expressionContains information of the physical model, such as variables and equations for generating element stiffness matrix.
Fortran Segmentscodes that realize the physical model at element level.
variables
equation
Automated Code Generator
Step 1: From PDE expression to Fortran segments
Segment 1
Segment 2
Segment 3
Segment 4
www.geongrid.orgCYBERINFRASTRUCTURE FOR THE GEOSCIENCES
Step 2: From algorithm expression to Fortran segments
do i=1,k do j=1,k estifn(i,j)=0.0 end do end do do i=1,k estifn(i,i)=estifn(i,i) do j=1,k estifn(i,j)=estifn(i,j)+es(i,j) end do end do
U(IDGF,NODI)=U(IDGF,NODI) *+ef(i)
defistif Smass Mload Ftype emdty lstep 0
equationmatrix = [S]FORC=[F]
SOLUTION Uwrite(s,unod) U
end
Algorithm ExpressionContains information for forming global
stiffness matrix for the model.
Fortran Segmentscodes that realize the physical model at global level.
Stiffness matrix Segment 5
Segment 6
www.geongrid.orgCYBERINFRASTRUCTURE FOR THE GEOSCIENCES
SUBROUTINE ETSUB(KNODE,KDGOF,IT,KCOOR,KELEM,K,KK, *NUMEL,ITYP,NCOOR,NUM,TIME,DT,NODVAR,COOR,NODE,#SUBET.sub *U) implicit double precision (a-h,o-z) DIMENSION NODVAR(KDGOF,KNODE),COOR(KCOOR,KNODE), *U(KDGOF,KNODE),EMATE(300),#SUBDIM.sub *R(500),PRMT(500),COEF(500),LM(500)#SUBFORT.sub#ELEM.subC WRITE(*,*) 'ES EM EF ='C WRITE(*,18) (EF(I),I=1,K)#MATRIX.sub L=0 M=0 I=0 DO 700 INOD=1,NNE ……… U(IDGF,NODI)=U(IDGF,NODI)#LVL.sub DO 500 JNOD=1,NNE ………500 CONTINUE700 CONTINUE ……… return end
Program StencilFortran Segments generated
Step 3: Plug Fortran segments into a stencil, forming final FE program
Segment 1
Segment 2
Segment 4
Segment 3
Segment 5
Segment 6
…………..
www.geongrid.orgCYBERINFRASTRUCTURE FOR THE GEOSCIENCES
Examples
• Western US tectonics
• Deformation of Asian continent
• Stress evolution and strain localization in the San Andreas Fault
www.geongrid.orgCYBERINFRASTRUCTURE FOR THE GEOSCIENCES
A Preliminary Finite Element Model of Active A Preliminary Finite Element Model of Active Crustal Deformation in the Western USCrustal Deformation in the Western US
www.geongrid.orgCYBERINFRASTRUCTURE FOR THE GEOSCIENCES
www.geongrid.orgCYBERINFRASTRUCTURE FOR THE GEOSCIENCES
The Power of GEON Cluster NodeThe Power of GEON Cluster NodeOriginal series model (single CPU)
•Less than 3000 elements
•Three layers in R-direction
•2 min for per time step
Preliminary parallel model (16-nodes, 32 CPUs)
•More than 800,000 unstructured elements
•Major Faults and more deformation zones
•Subduction of Juan de Fuca slab
•21 layers in R-direction
(x 40vertical topographic exaggeration)
www.geongrid.orgCYBERINFRASTRUCTURE FOR THE GEOSCIENCES
www.geongrid.orgCYBERINFRASTRUCTURE FOR THE GEOSCIENCES
www.geongrid.orgCYBERINFRASTRUCTURE FOR THE GEOSCIENCES
The model now allows simulation of large scale The model now allows simulation of large scale continental deformation with unprecedented detailcontinental deformation with unprecedented detail
www.geongrid.orgCYBERINFRASTRUCTURE FOR THE GEOSCIENCES
The model now allows simulation of large scale The model now allows simulation of large scale continental deformation with unprecedented detailcontinental deformation with unprecedented detail
www.geongrid.orgCYBERINFRASTRUCTURE FOR THE GEOSCIENCES
Ongoing Effort:Ongoing Effort: Toward a new 3D model of Toward a new 3D model of continental deformation in Asiacontinental deformation in Asia
www.geongrid.orgCYBERINFRASTRUCTURE FOR THE GEOSCIENCES
Predicted surface shear stress
www.geongrid.orgCYBERINFRASTRUCTURE FOR THE GEOSCIENCES
Loading the San Andreas Fault by relative PA-NA motion
•Fully 3D•Dynamic•Plastic-viscoelastic•Co-seismic/interseismic cyclesfrom seconds to 104 years•Parallel computing on PC clusters
www.geongrid.orgCYBERINFRASTRUCTURE FOR THE GEOSCIENCES
Comparison of predicted surface velocity and GPS data
www.geongrid.orgCYBERINFRASTRUCTURE FOR THE GEOSCIENCES
Predicted maximum shear stress
www.geongrid.orgCYBERINFRASTRUCTURE FOR THE GEOSCIENCES
www.geongrid.orgCYBERINFRASTRUCTURE FOR THE GEOSCIENCES
Predicted rate of plastic strain energy release outside the SAF
www.geongrid.orgCYBERINFRASTRUCTURE FOR THE GEOSCIENCES
Dream on…
• Integrating the community modeling environment with the geoscience data cyberinfrastructure;
• Grid computation and data integration;
• Automated (optimized?) work flow management (the Kepler system?)
www.geongrid.orgCYBERINFRASTRUCTURE FOR THE GEOSCIENCES