Search NASA⌕ Search

SEARCH · Search NASA

Results for “open source codes”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 55 records · Page 3
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Optical second harmonic generation in anisotropic multilayers with complete multireflection of linear and nonlinear waves using ♯SHAARP.ml package

Abstract Optical second harmonic generation (SHG) is a nonlinear optical effect widely used for nonlinear optical microscopy and laser frequency conversion. Closed-form analytical solution of the nonlinear optical responses is essential for evaluating materials whose optical properties are unknown a priori. A recent open-source code, ♯SHAARP. si , can provide such closed form solutions for crystals with arbitrary symmetries, orientations, and anisotropic properties at a single interface. However, optical components are often in the form of slabs, thin films on substrates, and multilayer heterostructures with multiple reflections of both the fundamental and up to ten different SHG waves at each interface, adding significant complexity. Many approximations have therefore been employed in the existing analytical approaches, such as slowly varying approximation, weak reflection of the nonlinear polarization, transparent medium, high crystallographic symmetry, Kleinman symmetry, easy crystal orientation along a high-symmetry direction, phase matching conditions and negligible interference among nonlinear waves, which may lead to large errors in the reported material properties. To avoid these approximations, we have developed an open-source package named Second Harmonic Analysis of Anisotropic Rotational Polarimetry in Multilayers (♯SHAARP. ml ). The reliability and accuracy are established by experimentally benchmarking with both the SHG polarimetry and Maker fringes using standard and commonly used nonlinear optical materials as well as twisted 2-dimensional heterostructures.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Enhancing ChatPORT with CUDA-to-SYCL Kernel Translation Capability

Large Language Models (LLMs) have shown strong capabilities in general code translation. However, code translation involving parallel programming models remains largely unexplored. This work enhances the capabilities of code LLMs in CUDA-to-SYCL kernel translation with parameter-efficient fine-tuning. The resultant fine-tuned LLM, called ChatPORT, is an effort to provide high-fidelity translations from one programming model to another. We describe the preparation of datasets from heterogeneous computing benchmarks for model fine-tuning and testing, the parameter-efficient fine-tuning of 19 open-source code models ranging in size from 0.5 to 34 billion parameters and evaluate the correctness rates of the SYCL kernels by the fine-tuned models. The experimental results show that most code models fail to translate CUDA codes to SYCL correctly. However, fine-tuning these models using a small set of CUDA and SYCL kernels can enhance the capabilities of these models in kernel translation. Depending on the sizes of the models, the correctness rate ranges from 19.9% to 81.7% for a test dataset of 62 CUDA kernels.

Jin, Zheming [ORNL] (ORCID:000000027197780X)↗
DOE OSTI2020

CTF Theory Manual (V.4.2)

Coolant-Boiling in Rod Arrays—Two Fluids (COBRA-TF) is a thermal/hydraulic (T/H) simulation code designed for light water reactor (LWR) vessel analysis. It uses a two-fluid, three-field (i.e. fluid film, fluid drops, and vapor) modeling approach. Both sub-channel and 3D Cartesian forms of 9 conservation equations are available for LWR modeling. The code was originally developed by Pacific Northwest Laboratory in 1980 and had been used and modified by several institutions over the last few decades. COBRA-TF also found use at the Pennsylvania State University (PSU) by the Reactor Dynamics and Fuel Modeling Group (RDFMG) and has been improved, updated, and subsequently re-branded as CTF. As part of the improvement process, it was necessary to generate sufficient documentation for the open-source code which had lacked such material upon being adopted by RDFMG. This document serves mainly as a theory manual for CTF, detailing the many two-phase heat transfer, drag, and important accident scenario models contained in the code as well as the numerical solution process utilized. Coding of the models is also discussed, all with consideration for updates that have been made when transitioning from COBRA-TF to CTF. Further documentation outside of this manual is also available at RDFMG which focus on code input deck generation and source code global variable and module listings

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2022

HiPACE++: A portable, 3D quasi-static particle-in-cell code

Modeling plasma accelerators is a computationally challenging task and the quasi-static particle-in-cell algorithm is a method of choice in a wide range of situations. In this work, we present the first performance-portable, quasi-static, three-dimensional particle-in-cell code HiPACE++. By decomposing all the computation of a 3D domain in successive 2D transverse operations and choosing appropriate memory management, HiPACE++ demonstrates orders-of-magnitude speedups on modern scientific GPUs over CPU-only implementations. The 2D transverse operations are performed on a single GPU, avoiding time-consuming communications. The longitudinal parallelization is done through temporal domain decomposition, enabling near-optimal strong scaling from 1 to 512 GPUs. HiPACE++ is a modular, open-source code enabling efficient modeling of plasma accelerators from laptops to state-of-the-art supercomputers.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

High-order finite element method for atomic structure calculations

We introduce featom, an open source code that implements a high-order finite element solver for the radial Schrödinger, Dirac, and Kohn-Sham equations. The formulation accommodates various mesh types, such as uniform or exponential, and the convergence can be systematically controlled by increasing the number and/or polynomial order of the finite element basis functions. The Dirac equation is solved using a squared Hamiltonian approach to eliminate spurious states. Here, to address the slow convergence of the $κ=±1$ states due to divergent derivatives at the origin, we incorporate known asymptotic forms into the solutions. We achieve a high level of accuracy (10 -8 Hartree) for total energies and eigenvalues of heavy atoms such as uranium in both Schrödinger and Dirac Kohn-Sham solutions. We provide detailed convergence studies and computational parameters required to attain commonly required accuracies. Finally, we compare our results with known analytic results as well as the results of other methods. In particular, we calculate benchmark results for atomic numbers (Z) from 1 to 92, verifying current benchmarks. We demonstrate significant speedup compared to the state-of-the-art shooting solver dftatom. An efficient, modular Fortran 2008 implementation, is provided under an open source, permissive license, including examples and tests, wherein particular emphasis is placed on the independence (no global variables), reusability, and generality of the individual routines.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2023

Combustion of petroleum-based transportation fuels and their blends with biofuels: a new approach for developing surrogates and understanding the effects of blending

Accurate and versatile reaction mechanisms are necessary to simulate combustion of transportation fuels blended with biofuels for improving fuel efficiency, reducing fuel consumption and mitigating the formation and emissions of combustion products harmful to the atmosphere. This project involved a combined effort of simulation and experimentation of liquid fuel burning using a combustion configuration amenable to numerical modeling with a level of detail not previously achieved for blends of biofuels with gasoline certification and surrogate fuels. The configuration of an isolated droplet burning under conditions where gas transport arises solely from fuel evaporation was selected as the platform for experiments and numerical modeling. The spherical symmetry and one-dimensional gas transport that results enabled simulating the droplet burning process with an precedented level of detail. Processes associated with unsteady gas and liquid transport, formation of particulate and gaseous products, radiation, multicomponent phase equilibrium at the droplet surface and moving boundary effects from droplet evaporation, were incorporated in a single numerical framework. The numerical model was based on the open source code OPENSmoke++ (OS) adapted to incorporate the effects noted above. In addition to simulation, an experimental design was developed for the isolated droplet configuration to acquire the data needed for validating the numerical model. Several broad accomplishments of the project were the following: demonstrating generally excellent agreement between measured and simulated combustion parameters for the fuel blends examined that included mixtures of heptane and isobutanol as a model system, and surrogates comprised of up to seven miscible components mixed with ethanol or isobutanol as representative biofuel additives; using a new approach to validate reaction mechanisms of biofuel blends which incorporated fuel evaporation into the process along with developing an experimental design for acquiring data to compare with simulations; and showing that a wealth of information could be obtained on the combustion physics of biofuels from experiments requiring volumes on the order of only nanoliters at a time thus opening the way to evaluating biofuels synthesized by new processes early in development.

02 PETROLEUM↗
DOE OSTI2021

Coupling particle scale model and SuperDEM–CFD for multiscale simulation of biomass pyrolysis in a packed bed pyrolyzer

An efficient biomass pyrolysis process requires a comprehensive understanding of the chemical and physical phenomena that occur at multi-length and time scales. In this study, a multiscale computational approach was developed and validated for biomass pyrolysis in a packed-bed reactor by integrating pyrolysis kinetics, a particle scale model, and SuperDEM-CFD in open-source code MFiX. A one-dimensional particle-scale model that discretizes the characteristic length of biomass particle into layers was developed to predict the intraparticle phenomena inside a single particle. The 1D model was validated by comparing it with a single biomass particle pyrolysis experiment. Here, a recently developed SuperDEM-CFD model was employed to simulate the non-spherical particle-particle contact and fluid-particle interaction. The coupled model was applied to simulate the pyrolysis of cubic biomass particles in a packed bed and validated by comparing with experimental data. Simulation with and without particle-scale model was compared, and the effect of the gas-solid heat transfer models was also investigated.

09 BIOMASS FUELS↗
DOE OSTI2021

The ezSpectra suite: An easy‐to‐use toolkit for spectroscopy modeling

Abstract A molecule's spectrum encodes information about its structure and electronic properties. It is a unique fingerprint that can serve as a molecular ID. Quantum chemistry calculations provide key ingredients for interpreting spectra, but modeling the spectra rarely ends there; it requires additional steps that entail combined treatments of electronic and nuclear degrees of freedom and account for specifics of the experimental setup (light energy, polarization, averaging over molecular orientations, temperature, etc.). This Software Focus article describes the ezSpectra suite, which currently comprises two stand‐alone open‐source codes: ezFCF and ezDyson . ezFCF calculates Franck–Condon factors, which yield vibrational progressions for polyatomic molecules, within the double‐harmonic approximation. ezDyson calculates absolute cross‐sections for photodetachment/photoionization processes and photoelectron angular distributions using Dyson orbitals computed by a quantum chemistry program. This article is categorized under: Electronic Structure Theory > Ab Initio Electronic Structure Methods Theoretical and Physical Chemistry > Spectroscopy Software > Simulation Methods

Gozem, Samer↗
DOE OSTI2022

Generalized thermo-mechanical framework for heterogeneous materials through asymptotic homogenization

Abstract A fundamental understanding of the interaction between microstructure and underlying physical mechanisms is essential, especially for developing more accurate multi-physics models for heterogeneous materials. Effects of microstructure on the material response at the macroscale are modeled by using the generalized thermomechanics. In this study, strain gradient theory is employed as a higher-order theory on the macroscale with thermodynamics modeled as a first-order theory on the microscale. Hence, energy depends only on the temperature such that we circumvent an extension of Fourier’s law and analyze the “simplest” thermo-mechanical model in strain gradient elasticity. Developing multiphysics models for heterogeneous materials is indeed a challenge and even this “simplest” model in generalized thermomechanics creates dozens of parameters to be determined. We develop a thermo-mechanical framework, in which microstructure is modeled as a periodic structure and through asymptotic homogenization approach, higher-order parameters at macroscopic scale are calculated. To illustrate the importance of higher-order parameters in overall thermo-mechanical response of a heterogeneous materials, finite element method (FEM) is employed with the aid of open-source codes (FEniCS). Verification example of a bulk system and several case studies of porous structures demonstrate how such numerical framework can be beneficial in the design of materials with tailored microstructures.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

Fluid Dynamics Effects on Microstructure Prediction in Single-Laser Tracks for Additive Manufacturing of IN625

Single-track laser fusion were simulated using a heat-transfer-solidification-only (HTS) model and its extension with fluid dynamics (HTS_FD) model using a parallel open-source code, which included laminar fluid dynamics, flat-free surface of the molten alloy, heat transfer, phase-change, evaporation, and surface tension phenomena. The results illustrate that the fluid dynamics affects the solidification and ensuing microstructure. For the HTS_FD simulations, thermal gradient, G was found to exhibit a maximum at the extremity of the solidified pool ( i.e. , at the free surface), while for HTS simulations, G exhibited a maximum around the entire edge of the solidified pool. HTS_FD simulations predicted a wider range of cooling rates than the HTS simulations, exhibited an increased spread in the solidification speed, V variation within the melt-pool with respect to the HTS model results. Primary dendrite arm spacing (PDAS) were evaluated based on power law correlations and marginal stability theory models using the ( G , V ) from HTS and HTS_FD simulations to quantify the effect of the fluid dynamics on the microstructure. At low-laser powers and low-scan speeds, the PDAS obtained with the fluid dynamics model (HTS_FD) was larger by more than 30 pct with respect to the PDAS calculated with the simple HTS model. A new PDAS correlation, i.e. , \( \lambda_{1} \left[ {\mu {\text{m}}} \right] = 832\;G\left[ {\text{K/m}} \right]^{ - 0.5} V\left[ {\text{m/s}} \right]^{ - 0.25} \) , which uses the ( G , V ) results from the HTS_FD model was developed and validated against experimental results.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Validation and application of a multiphase CFD model for hydrodynamics, temperature field and RTD simulation in a pilot-scale biomass pyrolysis vapor phase upgrading reactor

Accurate prediction of transport phenomena is critical for VPU reactor design, optimization, and scale-up. The current study focused on the validation and application of a multiphase CFD model within an open-source code MFiX for hydrodynamics, temperature field, and residence time distribution (RTD) simulation in a non-reacting circulating fluidized bed riser for biomass pyrolysis vapor phase upgrading (VPU). First, an Eulerian-Eulerian approach three-dimensional CFD model was employed to simulate the pilot-scale VPU riser on the supercomputer Joule. Excellent quantitative agreement between experimental and simulated results was achieved for pressure drops and temperature field in a range of operating conditions. Then the validated multiphase CFD model was applied to predict gas and solid residence time distributions (RTDs) since prediction and analysis of RTD is an important tool to study the complex multiphase flow behavior and mixing inside chemical reactors. The predictions show that solid mean residence time is 3.5 times the gas residence time; the solid RTD is more sensitive to the process gas flow rate than the solids circulation rate.

09 BIOMASS FUELS↗
DOE OSTI2022

Fast estimation of reaction rates in spherical and non-spherical porous catalysts

Here we present a methodology for modeling multi-step reaction rates in porous catalyst particles for use in CFD-DEM and two fluid models. Single-step effectiveness factors based on a Thiele modulus, while useful, cannot accurately capture the cascading reaction systems common in high temperature vapor-phase chemical reactors like fluidized catalytic cracking units and catalytic biomass fast pyrolysis systems. Instead, multi-step effectiveness vectors derived from steady-state solutions to the governing reaction-diffusion equations are needed. Solutions for various catalyst shapes are presented, including spheres, cylinders, and prisms. Computational challenges inherent in repeated evaluation of reaction rates with diffusion limitations are discussed, and an efficient implementation based on pre-computed lookup tables is proposed and demonstrated on a simulation of a fluidized bed reactor. Open-source code is provided for the compilation of reaction rate tables for use in ODE, DEM, and two-fluid models.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2020

A computationally-efficient method for flamelet calculations

A new open-source code for the simulation of the diffusion flamelet equations is proposed. Emphasis is placed on using an approximate Jacobian to reduce the computational cost of the matrix operations. Performance of the proposed solvers is tested by performing flamelet calculations with kinetic mechanisms of varying sizes. For the unity Lewis number equations, the present iterative Newton solver using an approximate Jacobian greatly outperforms direct Newton solvers using exact Jacobians. The computation cost scales linearly with the number of species, leading to a reduction in solution times by two orders of magnitude for mechanisms containing thousands of species. The applicability of the Jacobian approximations to the solution of the non-unity Lewis number flamelet equations is assessed. The approximations are generally inadequate to solve the full non-unity Lewis number equations but can be used in some applications depending on the balance of terms in the flamelet equations. As an example, the flamelet solver is applied to the study of sooting tendencies in laminar co-flow diffusion flames where modified non-unity Lewis number flamelet equations, previously shown to accurately reproduce experimentally-measured Yield Sooting Indices (YSI), are solved. Here, the accelerated flamelet solver is well suited for sensitivity analysis and uncertainty quantification with large detailed kinetic mechanisms, tasks for which the computational cost was previously prohibitive.

42 ENGINEERING↗
DOE OSTI2025

A generalized and adaptable tensor-contraction-based cluster expansion formalism for multicomponent solids

Density functional theory (DFT)-based simulations of materials have first-principles accuracy, but are very computationally expensive. For simulating various properties of multi-component alloys, the cluster expansion (CE) technique has served as the standard workaround to improve computational efficiency. However, the standard CE technique is difficult to extend to exotic and/or low-symmetry lattices, often implemented via iteration over particular cluster types, which must be enumerated per lattice structure. In this work, we introduce the tensor cluster expansion (TCE), implemented in the open-source code tce-lib, which maps correlation functions to mixed tensor contractions, eliminating the need to iterate over cluster types and additionally making the calculation of correlation functions well-suited for massively parallel architectures like GPUs. We show that local interaction energies are an immediate consequence of the TCE formalism, yielding nearly $\mathcal{O}$(1) energy difference calculations. We then use this formalism to fit CE models for the TaW and CoNiCrFeMn systems, and use these models to respectively compute the enthalpy of mixing curve and Cowley short-range order parameters, showing excellent agreement with ground truth data.

Cluster expansion↗
DOE OSTI2020

A cross-sectional aeroelastic analysis and structural optimization tool for slender composite structures

A fully open-source available framework for the parametric cross-sectional analysis and design optimization of slender composite structures, such as helicopter or wind turbine blades, is presented. The framework—Structural Optimization and Aeroelastic Analysis (SONATA)—incorporates two structural solvers, the commercial tool VABS, and the novel open-source code ANBA4. SONATA also parameterizes the design inputs, postprocesses and visualizes the results, and generates the structural inputs to a variety of aeroelastic analysis tools. It is linked to the optimization library OpenMDAO. This work presents the methodology and explains the fundamental approaches of SONATA. Structural characteristics were successfully verified for both VABS and ANBA4 using box beam examples from literature, thereby verifying the parametric approach to generating the topology and mesh in a cross section as well as the solver integration. The framework was furthermore exercised by analyzing and evaluating a fully resolved highly flexible wind turbine blade. Computed structural characteristics correlated between VABS and ANBA4, including off-diagonal terms. Stresses, strains, and deformations were recovered from loads derived through coupling with aeroelastic analysis. The framework, therefore, proves effective in accurately analyzing and optimizing slender composite structures on a high-fidelity level that is close to a three-dimensional finite element model.

42 ENGINEERING↗
DOE OSTI2024

Software stewardship and advancement of a high-performance computing scientific application: QMCPACK

Here, we provide an overview of the software engineering efforts and their impact in QMCPACK, a production-level ab-initio Quantum Monte Carlo open-source code targeting high-performance computing (HPC) systems. Aspects included are: (i) strategic expansion of continuous integration (CI) targeting CPUs, using GitHub Actions own runners, and NVIDIA and AMD GPUs used in pre-exascale systems, (ii) incremental reduction of memory leaks using sanitizers, (iii) incorporation of Docker containers for CI and reproducibility, and (iv) refactoring efforts to improve maintainability, testing coverage, and memory lifetime management. We quantify the value of these improvements by providing metrics to illustrate the shift towards a predictive, rather than reactive, maintenance approach. Our goal, in documenting the impact of these efforts on QMCPACK, is to contribute to the body of knowledge on the importance of research software engineering (RSE) for the stewardship and advancement of community HPC codes to enable scientific discovery at scale.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Self-consistent solution of the Frank–Bilby equation for interfaces containing disconnections

The quantized Frank–Bilby equation can be used to identify interfacial line defect array configurations which relax the misorientation and/or misfit of a coherent crystalline interface. These line defect arrays may be comprised of dislocations and/or disconnections, which are interfacial steps with dislocation character. When an interface contains disconnections, solution of the quantized Frank–Bilby equation is complicated by the fact that the habit plane orientation is not known in advance because it depends on the unknown spacing of the disconnection array. We present a root-finding-based method for addressing this issue, enabling a self-consistent solution for arbitrary defect content. Our method has been implemented in an open-source code which enumerates all possible solutions given a list of candidate line defects. Two cases are presented employing the code: a misoriented FCC twin boundary and an FCC/BCC phase boundary with the Nishiyama-Wasserman orientation relationship. Both cases exhibit more than 10,000 solutions to the Frank–Bilby equation, with several hundred solutions categorized as ‘‘low energy’’ and thus plausible configurations for the actual interface. The resulting set of solutions can be utilized to predict and understand the properties of a given interface.

42 ENGINEERING↗
DOE OSTI2022

Simulation of dilatancy-controlled gas migration processes in saturated bentonite using a coupled multiphase flow and elastoplastic H 2 M model

Dilatancy-controlled gas flow in preferential pathways plays a key role in the safety analysis of radioactive waste repositories. This is particularly the case for bentonite, an often-preferred barrier material. Gas flow in preferential pathways is characterized by localization and spontaneous behavior, which is challenging to simulate in numerical models due to strong hydro-mechanical coupling. To analyze a laboratory experiment in the framework of the DECOVALEX-2023 project, this study introduced a new approach of combining continuous modelling methods with spatial material properties derived from material heterogeneities and experimental observations. The proposed model utilized hydro-mechanical spatial distributions, namely Young's modulus and gas entry pressure, and elastoplasticity combined with a linear swelling model. A conceptual strain-dependent permeability approach simulated dilatancy-controlled gas flow based on hydro-mechanical coupling. To test the effectiveness of the presented approach, a gas injection test in a compacted, saturated bentonite sample was simulated using the open-source code OpenGeoSys 5.8 and compared with experimental observations. The presented methodology is capable of simulating localized gas flow in preferential pathways. The spatial distributions of Young's modulus and gas entry pressure affect the swelling pressure, relative permeability and, in combination with the strain-dependent permeability model, also the intrinsic permeability.

42 ENGINEERING↗