Uintah and Related C-SAFE Publications


Y. Guo, J.A. Nairn. “Calculation of J-Integral and Stress Intensity Factors using the Material Point Method,” In Computer Modeling in Engineering and Sciences, Vol. 6, No. 3, pp. 295--308. 2004.


The Material Point Method (MPM), which is a particle-based, meshless method that discretizes material bodies into a collection of material points (the particles), is a new method for numerical analysis of dynamic solid mechanics problems. Recently, MPM has been generalized to include dynamic stress analysis of structures with explicit cracks. This paper considers evaluation of crack-tip parameters, such as J-integral and stress intensity factors, from MPM calculations involving explicit cracks. Examples for both static and dynamic problems for pure modes I and II or mixed mode loading show that MPM works well for calculation of fracture parameters. The MPM results agree well with results obtained by other numerical methods and with analytical solutions.

Y. He, T.R. Lutz, M.D. Ediger, C. Ayyagari, D. Bedrov, G.D. Smith. “NMR Experiments and Molecular Dynamics Simulations of the Segmental Dynamics of Polystyrene,” In Macromolecules, Vol. 37, No. 13, pp. 5032--5039. May, 2004.
DOI: 10.1021/ma049843r


We have performed NMR spin-lattice relaxation experiments and molecular dynamics (MD) computer simulations on atactic polystyrene (a-PS). The segmental correlation times of three different molecular weight a-PS (Mn = 1600, 2100, 10 900 g/mol) were extracted from NMR by measuring the 2H spin-lattice relaxation times (T1) over a broad temperature range (390-510 K). MD simulations of an a-PS melt of molecular weight 2200 g/mol were carried out at 475, 500, and 535 K. Comparisons between experiments and simulations show that the MD simulations reproduce both the shape of the P2(t) orientation autocorrelation function and its temperature dependence, while the simulated segmental correlation times are slower than experimental results by a factor of 1.8. If the simulations are rescaled by this factor, they reproduce both the experimental T1 values and the slight difference in dynamics between the backbone and side group of PS.

C.R. Johnson, R.S. MacLeod, S.G. Parker, D.M. Weinstein. “Biomedical Computing and Visualization Software Environments,” In Comm. ACM, Vol. 47, No. 11, pp. 64--71. 2004.

G. Krishnamoorthy, R. Rawat, P.J. Smith. “Parallel Computations of Radiative Heat Transfer Using the Discrete Ordinates Method,” In Numerical Heat Transfer, Part B: Fundamentals, Vol. 47, No. 1, pp. 19--38. 2004.
DOI: 10.1080/10407790490487451


The discrete ordinates method is spatially decomposed to solve the radiative transport equation on parallel computers. Mathematical libraries developed by third parties are used to solve the matrices that result during the solution procedure. The radiation component is verified by comparing computed values against a benchmark. Fixed and scaled problem size efficiencies are examined. Contrary to most previous studies, the parallel efficiencies did not depend strongly on the optical thickness of the medium for our model problem. Timing studies show that GMRES, BiCGSTAB iterative methods with block Jacobi preconditioning perform the best for solving these matrix systems.

A. Violi. “Modeling of Soot Particle Inception in Aromatic and Aliphatic Premixed Flames,” In Combustion and Flame, Vol. 139, No. 4, pp. 279--287. December, 2004.
DOI: 10.1016/j.combustflame.2004.08.013


The growth of hydrocarbon molecules up to sizes of incipient soot is computed in premixed laminar flames using kinetic Monte Carlo and molecular dynamic methodologies (AMPI code). This approach is designed to preserve atomistic scale structure (bonds, bond angles, dihedral angles) as soot precursors evolve into three-dimensional structures. Application of this code to aliphatic (acetylene) and aromatic (benzene) flame environments is able to explain results in the literature on the differences in properties of soot precursors from these two classes of flames, particularly relating to H/C ratio, particle sphericity, and depolarization ratio.

A. Violi, T.N. Truong, A.F. Sarofim. “Kinetics of Hydrogen Abstraction Reactions from Polycyclic Aromatic Hydrocarbons by H Atoms,” In Journal of Physical Chemistry, A, Vol. 108, No. 22, pp. 4846--4852. May, 2004.
DOI: 10.1021/jp026557d


An application of the Reaction Class Transition State Theory/Linear Energy Relationship (RC-TST/LER) is presented for the evaluation of the thermal rate constants of hydrogen abstraction reactions by H atoms from Polycyclic Aromatic Hydrocarbons (PAH). Two classes of reactions have been considered, namely hydrogen bonded to six- and five-membered rings, respectively, and twenty-two reactions have been used to develop the RC-TST/LER parameters. B3LYP and BH&HLYP density functional theory methods were used to calculate necessary potential energy surface information. Detailed analyses of RC-TST/LER reaction factors lead to the conclusion that rate constants for any reaction in these two classes can be approximated by those of its corresponding principal reaction corrected by the reaction symmetry factor. Specifically, for hydrogen abstraction from six-membered rings such as naphthalene and pyrene, k(T) = (σ/σH+C6H6) kH+C6H6 = (σ/6)1.42 × 108T1.77 exp(−6570/T)(cm3/mol·s), and for hydrogen abstraction from five-membered rings such as acenaphthylene and acephenanthrylene, k(T) = (σ/σH+C12H8) kH+C12H8 = (σ/2)3.27 × 108T1.71 exp(−8170/T) (cm3/mol·s), where σ is the reaction symmetry number.


B. Banerjee, D.O. Adams. “Micromechanics-Based Prediction of Effective Elastic Properties of Polymer Bonded Explosives,” In Physica, B, Vol. 338, No. 1-4, pp. 8--15. October, 2003.
DOI: 10.1016/S0921-4526(03)00451-4


Polymer bonded explosives are particulate composites containing a high volume fraction of stiff elastic explosive particles in a compliant viscoelastic binder. Since the volume fraction of particles can be greater than 0.9 and the modulus contrast greater than 20 000, rigorous bounds on the elastic moduli of the composite are an order of magnitude different from experimentally determined values. Analytical solutions are also observed to provide inaccurate estimates of effective elastic properties. Direct finite element approximations of effective properties require large computational resources because of the complexity of the microstructure of these composites. An alternative approach, the recursive cells method (RCM) is also explored in this work. Results show that the degree of discretization and the microstructures used in finite element models of PBXs can significantly affect the estimated Young's moduli.

B. Banerjee, C.M. Cady, D.O. Adams. “Micromechanics Simulations of Glass--estane Mock Polymer Bonded Explosives,” In Modeling and Simulation in Materials Science and Engineering, Vol. 11, No. 4, pp. 457--475. 2003.
DOI: 10.1088/0965-0393/11/4/304


Polymer bonded explosives (PBXs) are particulate composites containing explosive particles and a continuous binder. The elastic modulus of the particles, at room temperature and higher, is often three to four orders of magnitude higher than that of the binder. Additionally, the explosive particles occupy high volume fractions, often greater than 90%. Both experimental and numerical determination of macroscopic properties of these composites is difficult. High modulus contrast mock PBXs provide a means of relatively inexpensive experimentation and validation of numerical approaches to determine properties of these materials. The goal of this investigation is to determine whether the effective elastic properties of monodisperse glass–estane mock PBXs can be predicted from two-dimensional micromechanics simulations using the finite element (FEM) method. In this study, the effect of representative volume element (RVE) size on the prediction of two-dimensional properties is explored. Two-dimensional estimates of elastic properties are compared with predictions from three-dimensional computations and with experimental data on glass–estane composites containing three different volume fractions of spherical glass beads. The effect of particle debonding on the effective elastic properties is also investigated using contact analyses. Results show that two-dimensional unit cells containing 10–20 circular particles are adequate for modelling glass–estane composites containing less than 60% glass particles by volume. No significant difference is observed between properties predicted by the two- and three-dimensional models. FEM simulations of RVEs, containing particles that are perfectly bonded to the binder, produce estimates of Young's modulus that are higher than the experimental data. Incorporation of debonding between particles and the binder causes the effective Young's modulus to decrease. However, the results suggest that cracks in the composite may play a significant role in determining the effective properties of mock polymer bonder explosives composed of glass and estane. The FEM simulations indicate that two-dimensional models that incorporate debonds and cracks can be used to obtain accurate estimates of the effective properties of glass–estane composites and possibly of PBXs.

W. Ciro, E.G. Eddings, A.F. Sarofim. “Fast Cookoff Tests Report,” C-SAFE Internal Report, C-SAFE, October, 2003.

L.E. Fried, R. Manaa, J.P. Lewis. “Modeling the Reactions of Energetic Materials in the Condensed Phase,” In Overviews of Recent Research on Energetic Materials, Edited by D. Thompson and T. Brill and R. Shaw, 2003.

J.D. de St. Germain, S.G. Parker. “Software Integration in an Academic Environment,” In Software Quality Forum (SQF), Arlington, Virginia, pp. (published on CD). March, 2003.

J.D. de St. Germain, A. Morris, S.G. Parker, A.D. Malony, S. Shende. “Performance Analysis Integration in the Uintah Software Development Cycle,” In International Journal of Parallel Programming, Vol. 31, No. 1, pp. 35--53. 2003.


The increasing complexity of high-performance computing environments and programming methodologies presents challenges for empirical performance evaluation. Evolving parallel and distributed systems require performance technology that can be flexibly configured to observe different events and associated performance data of interest. It must also be possible to integrate performance evaluation techniques with the programming paradigms and software engineering methods. This is particularly important for tracking performance on parallel software projects involving many code teams over many stages of development. This paper describes the integration of the TAU and XPARE tools in the Uintah Computational Framework (UCF). Discussed is the use of performance mapping techniques to associate low-level performance data to higher levels of abstraction in UCF and the use of performance regression testing to provide a historical portfolio of the evolution of application performance. A scalability study shows the benefits of integrating performance technology in building large-scale parallel applications.

Keywords: uintah

J.E. Guilkey, J.A. Weiss. “Implicit Time Integration with the Material Point Method,” In International Journal for Numerical Methods in Engineering, Vol. 57, pp. 1323--1338. 2003.

J.E. Guilkey, T. Harman, A. Xia, B. Kashiwa, P.A. McMurtry. “An Eulerian-Lagrangian Approach for Large Deformation Fluid Structure Interaction Problems, Part 1: Algorithm Development,” In Fluid Structure Interaction II, In Fluid Structure Interaction II: Proceedings of Fluid Structure Interaction , WIT Press, pp. 14. 143--156, 2003.
DOI: 10.2495/FSI030141

T. Harman, J.E. Guilkey, B. Kashiwa, J. Schmidt, P. McMurtry. “An Eulerian-Lagrangian Approach for Large Deformation Fluid Structure Interaction Problems, Part 2: Multi-Physics Simulations within a Modern Computational Framework,” In Fluid Structure Interaction II, WIT Press, pp. 157--166. 2003.


A tightly coupled fluid-structure interaction (FSI) solution technique incorporating fluid and solid mechanics, phase change and chemical reactions is presented. The continuum equations are solved with a cell-centered, multi-material ICE solution method. This formulation is integrated with a Lagrangian, particle based, solid mechanics technique, known as the Material Point Method, as described by Kashiwa et al. [1] and Guilkey et al. [2]. The combined method can handle large deformations and phase change within a single grid, without the need of separate domains for fluids and solids, or the passing of boundary conditions. This paper discusses algorithmic issues involved in accounting for chemical reactions and phase transition among material phases (e.g., solid → gas). Validation is presented as are simulations showing large deformation with phase change. These simulations were performed within a computational framework that contains tools for parallelization, performance analysis, data management, algorithm integration, and data visualization. Features of this framework are described.

G. Krishnamoorthy, J.M. Veranth. “Computational Modeling of CO-CO2 Ratio Inside Single Char Particles during Pulverized Coal Combustion,” In Energy and Fuels, Vol. 17, No. 5, pp. 1367--1371. August, 2003.
DOI: 10.1021/ef030006k


A recently developed model was used to study the CO/CO2 ratio inside a burning pulverized coal particle, to better understand the effect of bulk gas composition on the equilibrium partial pressure of reduced metal species at the surface of ash inclusions. The motivation was to improve the ability to model submicrometer particle formation by ash vaporization, as a function of furnace conditions. Assumptions for the CO/CO2 ratio that have been made in previous studies are compared to predictions from a psuedo-steady-state model for a single porous particle that considers homogeneous and heterogeneous reaction kinetics and mass transfer both in particle pores and in the boundary layer. This is the first publication of model predictions for the CO/CO2 ratio as a function of radius for a coal char particle in a furnace with a bulk gas CO2 concentration in the range of 0%−79%. A method is proposed for summarizing the effects on the CO/CO2 ratio that are due to changes in the bulk furnace gas O2 and CO2 concentration, furnace temperature, and particle size, using an empirical equation that is suitable for incorporation as a submodel into comprehensive computational fluid dynamics-based codes for combustion simulation. Trends from the model simulations show general agreement with experimental data; however, the accuracy of the predictions is limited by the lack of fuel-specific input data.

J.P. Lewis. “Energetics of intermolecular HONO formation in condensed-phase octahydro-1,3,5,7-tetranitro-1,3,5,7-tetrazocine (HMX),” In Chemical Physics Letters, Vol. 371, No. 5-6, pp. 588--593. April, 2003.
DOI: 10.1016/S0009-2614(03)00309-9


We present preliminary work on the formation of HONO in condensed-phase HMX for the three pure polymorphic phases. Our results show that the energetics of the intermolecular hydrogen transfer (to form HONO on an adjacent molecule) is more favorable in δ-HMX than in α-HMX and β-HMX. The energetics of this intermolecular hydrogen transfer process follow the same trends of HMX sensitivity, where δ-HMX is found to be more sensitive to explosion compared to β-HMX.

S.G. Parker. “C-SAFE Uses Linux HPCC in Fire Research,” In Syllabus, Technology for Higher Education, Vol. 16, No. 7, Feburary, 2003.

T.D. Sewell, R. Menikoff, D. Bedrov, G.D. Smith. “A Molecular Dynamics Simulation Study of Elastic Properties of HMX,” In Journal of Chemical Physics, Vol. 119, No. 14, pp. 7471--7426. 2003.
DOI: 10.1063/1.1599273

S. Shellman, K. Sikorski. “A Recursive Algorithm for the Infinity-Norm Fixed Point Problem,” In Journal of Complexity, Vol. 19, No. 6, pp. 799--834. 2003.
DOI: 10.1016/j.jco.2003.06.001