Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “open source model”

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 1,585 records · Page 88

High-Energy Emission From Millisecond Pulsars

The X-ray and gamma-ray spectrum of rotation-powered millisecond pulsars is investigated in a model for acceleration and pair cascades on open field lines above the polar caps. Although these pulsars have low surface magnetic fields, their short periods allow them to have large magnetospheric potential drops, but the majority do not produce sufficient pairs to completely screen the accelerating electric field. In these sources, the primary and secondary electrons continue to accelerate to high altitude and their Lorentz factors are limited by curvature and synchrotron radiation reaction. The accelerating particles maintain high Lorentz factors and undergo cyclotron resonant absorption of radio emission, that produces and maintains a large pitch angle, resulting in a strong synchrotron component. The resulting spectra consist of several distinct components: curvature radiation from primary electrons dominating from 1 - 100 GeV, synchrotron radiation from primary and secondary electrons dominating up to about 100 MeV, and much weaker inverse-Compton radiation from primary electrons a t 0.1 - 1 TeV. We find that the relative size of these components depends on pulsar period, period derivative, and neutron star mass and radius with the level of the synchrotron component also depending sensitively on the radio emission properties. This model is successful in describing the observed X-ray and gamma-ray spectrum of PSR J0218+4232 as synchrotron radiation, peaking around 100 MeV and extending up to a turnover around several GeV. The predicted curvature radiation components from a number of millisecond pulsars, as well as the collective emission from the millisecond pulsars in globular clusters, should be detectable with AGILE and GLAST. We also discuss a hidden population of X-ray-quiet and radio-quiet millisecond pulsars which have evolved below the pair death line, some of which may be detectable by telescopes sensitive above 1 GeV. Subject headings: pulsars: general - radiation mechanisms: nonthermal - stars: neutron - gamma rays: theory

Harding, Alice K.↗

Physics-Informed Learning Machines for Multiscale and Multiphysics Problems (PHILMS) (Technical Report)

The research work at University of California Santa Barbara (UCSB) resulted in several new developments in the areas of scientific machine learning, numerical analysis, and practical methods for data-driven modeling, prediction, reductions, and simulation. Many of the projects were carried out in collaboration with members of the national laboratories at Sandia National Laboratories (SNL), Pacific Northwestern National Laboratories (PNNL), and other institutions. Results included developing new scientific machine learning methods, related theory and mathematical frameworks for analysis and training, data-driven numerical solvers, and related tools and software for scientific computation. During the support period, over 16+ papers were submitted for publication, and 4 open-source software packages were developed and released (available at http://atzberger.org/). In addition, 7+ students and 2 post-docs were mentored in collaboration with the laboratory staff for future careers in academia, government labs, and industry.

97 MATHEMATICS AND COMPUTING↗

Improving and Automating Building Model Data Exchange

There are many instances throughout a project’s lifecycle where there arises a need for quick and accurate risk assessment of building designs. For example, an unexpected design change during construction may necessitate structural engineers to perform a seismic risk assessment on analytical models of the updated building design using high fidelity structural analysis software, such as ANSYS or Abaqus. However, the efficiency of such workflows often depends upon the interoperability of architectural design software and structural analysis software. When the quality of this interoperability is lacking or even non-existent, the efficiency of virtual engineering workflows is hampered, which increases project costs. A McGraw Hill industry survey of professional users of Building Information Modeling (BIM) technologies found that there is high demand for BIM interoperability for structural analysis, but that the value/difficulty ratio is currently too low for practical use. There have been efforts by the academic community to facilitate model data exchange between the architectural design and structural analysis domains, but such solutions have not been widely adopted by industry, face technical challenges, and oftentimes are limited in applicability for users of various BIM software. Therefore, INL is developing capabilities to improve, automate, and generalize model data exchange between architectural BIM software (e.g., Revit) and structural analysis software (e.g., SAP2000, ANSYS). The goal is to help expedite and automate as much of the pre-processing step for creating analytical models in finite element analysis software as reasonably as possible. Such a "BIM-to-FEA" conversion tool should provide direct benefit to end-users through accuracy, automation, quick turn-around, and wide applicability. To generalize the application of this BIM-to-FEA conversion tool and increase its useability among the many different commercial BIM software currently used by industry, the program is being developed with the concept of openBIM. OpenBIM is the application of non-proprietary, open data standards that allow for BIM model data exchange in a format that is accessible, retainable, and useable for all users. The most widely used open, non-proprietary data exchange format for BIM is the Industry Foundation Classes (IFC) schema. IFC is developed by buildingSMART international and is ISO certified (ISO 16739-1:2018). The BIM-to-FEA conversion tool is being developed for compatibility with typical commercial building designs of steel framed structures. The tool is currently capable of importing architectural BIM data of framed building structures, recognizing and extracting the aspects of the model that are required for structural analysis, adjusting the connectivity of frame members, and finally exporting to an analytical model stored in the IFC format. The exported IFC analytical model can then be imported into various openBIM compliant software, such as SAP2000. Such capabilities have already been tested on commercial software, as shown above, and continue to be improved. Work is underway to test the conversion on various commercial BIM software, develop a user-friendly interface, incorporate the program into the broader DeepLynx data warehouse project being developed by INL, and to eventually open-source the tool for the benefit of the community. Future development of the tool envisions the ability for efficient iterative risk assessment of generative building designs, all within a workflow utilizing open-source tools. One such open-source tool will be MOOSE, an advanced finite element analysis tool developed at INL. The conversion tool will also branch out from typical commercial building designs and will aim to incorporate nuclear construction. The aim will be to convert both structural and non-structural components of nuclear facilities, such as curved concrete containment structures and piping systems, respectively.

97 MATHEMATICS AND COMPUTING↗

GLDAS-2 Land Surface Model Data and Data Services at NASA GES DISC

The goal of the NASA Global Land Data Assimilation System (GLDAS, https://ldas.gsfc.nasa.gov/gldas(https://ldas.gsfc.nasa.gov/gldas)) is to generate optimal fields of land surface states and fluxes by ingesting satellite- and ground-based observational data products, using advanced land surface modeling and data assimilation techniques (Rodell et al., 2004).The GLDAS dataset currently archived at and distributed by the NASA Goddard Earth Sciences Data and Information Services Center (GES DISC, https://disc.gsfc.nasa.gov/ (https://disc.gsfc.nasa.gov/)) is GLDAS Version 2 (GLDAS-2). It contains a series of output fields from the upgraded Noah-3.6, Catchment-F2.5, and VIC-4.1.2 Land Surface Models (LSMs) in the Land Information System (LIS-V7, https://lis.gsfc.nasa.gov/ (https://lis.gsfc.nasa.gov/)). GLDAS-2 has three components:GLDAS-2.0, GLDAS-2.1, and GLDAS-2.2. GLDAS-2.0 is forced entirely with the upgraded Princeton Meteorological ForcingV2.2 Dataset and provides a temporally consistent series from 1948 through 2014. GLDAS-2.1 is forced with a combination of model and observation data, with data spanning from 2000 to the present. The GLDAS-2.2 product suite uses data assimilation(DA), whereas the GLDAS-2.0 and GLDAS-2.1 products are "open-loop" (i.e., no data assimilation). The choice of forcing data, as well as DA observation source, variable, and scheme, varies for different GLDAS-2.2 products. The currently availableGLDAS-2.2 data contain a daily 0.25-degree output from the Catchment-F2.5 LSM in LIS-V7. The data are forced with the meteorological analysis fields from the operational European Centre for Medium-Range Weather Forecasts Integrated Forecasting System (ECMWF-IFS) and assimilated with GRACE and GRACE-FO data, ranging from February 1, 2003 to the present. The current GLDAS-2.0 and 2.1 Noah LSM data were reprocessed in November 2019 and January 2020 respectively and their data from Catchment and VIC LSMs are new to the GLDAS-2 collection. This presentation provides a summary of theGLDAS-2 data products, their land surface fields, and their related data services at the GES DISC; and a description of the majorGLDAS-2 climatological characteristics as well as the intercomparison with the data of the previous version.

Hydrology↗

MFC 5.0: An exascale many-physics flow solver

Many problems of interest in engineering, medicine, and the fundamental sciences rely on high-fidelity flow simulation, making performant computational fluid dynamics solvers a mainstay of the open-source software community. Previous work MFC 3.0 was made a published, documented, and open-source solver via Bryngelson et al. Comp. Phys. Comm. (2021) with numerous physical features, numerical methods, and scalable infrastructure. MFC 5.0 is a significant update to MFC 3.0, featuring a broad set of well-established and novel physical models and numerical methods, as well as the introduction of GPU and APU (or superchip) acceleration. Here, we exhibit state-of-the-art performance and ideal scaling on the first two exascale supercomputers, OLCF Frontier and LLNL El Capitan. Combined with MFC’s single-accelerator performance, MFC achieves exascale computation in practice, and achieved the largest-to-date public CFD simulation at 200 trillion grid points as a 2025 ACM Gordon Bell Prize finalist. New physical features include the immersed boundary method, N-fluid phase change, Euler–Euler and Euler–Lagrange sub-grid bubble models, fluid-structure interaction, hypo- and hyper-elastic materials, chemically reacting flow, two-material surface tension, magnetohydrodynamics (MHD), and more. Numerical techniques now represent the current state-of-the-art, including general relaxation characteristic boundary conditions, WENO variants, Strang splitting for stiff sub-grid flow features, and low Mach number treatments. Weak scaling to tens of thousands of GPUs on OLCF Summit and Frontier and LLNL El Capitan achieves efficiencies within 5% of ideal to over 90% of their respective system sizes. Strong scaling results for a 16-times increase in device count show parallel efficiencies over 90% on OLCF Frontier. MFC’s software stack has undergone further improvements, including continuous integration, which ensures code resilience and correctness through over 300 regression tests; metaprogramming, which reduces code length while maintaining performance portability; and code generation for computing chemical reactions

Computational fluid dynamics↗

Deep Electric Field Predictions by Drift-Reduced Braginskii Theory with Plasma-Neutral Interactions Based on Experimental Images of Boundary Turbulence

We present 2-dimensional turbulent electric field calculations via physics-informed deep learning consistent with (i) drift-reduced Braginskii theory under the framework of an axisymmetric fusion plasma with purely toroidal field and (ii) experimental estimates of the fluctuating electron density and temperature on open field lines obtained from analysis of gas puff imaging of a discharge on the Alcator C-Mod tokamak. The inclusion of effects from the locally puffed atomic helium on particle and energy sources within the reduced plasma turbulence model are found to strengthen correlations between the electric field and electron pressure. The neutrals are also directly associated with broadening the distribution of turbulent field amplitudes and increasing ${\bf E \times B}$ shearing rates. As a result, this demonstrates a novel approach in plasma experiments by solving for nonlinear dynamics consistent with partial differential equations and data without encoding explicit boundary nor initial conditions.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

PVade (PV Aerodynamic Design Engineering) [SWR-23-49]

PV Aerodynamic Design Engineering (PVade) is an open-source fluid-structure interaction (FSI) solver which can accurately simulate wind loads and aerodynamic stability in solar-tracking PV arrays. PVade’s pressure profiles and inertial load time series can be used as inputs into mechanical module models to study degradation mechanisms including cracking of cells, weathering of cracked cells, and glass breakage. This software enables PV plant owners to predict optimal high-wind stow strategies under a variety of weather and site conditions.

Young, Ethan↗

Acoustic Research under the Source Physics Experiment

The Source Physics Experiment series is a long-term research and development (R&D) effort under the U.S. Department of Energy’s National Nuclear Security Administration focused on improving the physical understanding of how chemical explosions generate seismoacoustic signals. Beginning in 2011, a series of subsurface chemical explosions in two different and highly contrasting geologies were conducted at the Nevada National Security Site in Nevada, USA with the objective of improving simulation and modeling approaches to explosion identification, yield estimation and other monitoring applications. The two executed phases of the series provide new explosion signature source data from a wide range of geophysical diagnostic equipment; recorded data from the test series is now openly available to the broader seismoacoustic community. This manuscript details the executed test series, deployed seismoacoustic networks, and summarizes major scientific achievements utilizing recorded signatures from the explosive tests.

45 MILITARY TECHNOLOGY, WEAPONRY, AND NATIONAL DEF↗

Open Building Operating System: An Open-Source Grid Responsive Control Platform for Buildings

Grid-interactive efficient buildings (GEBs) with flexible loads are a promising method to decarbonize buildings, shift loads during peak hours, and lower energy use and electricity costs. Despite the promising benefits of GEBs, automation systems that manage flexible loads in response to energy prices or other grid signals are still uncommon in small and medium commercial buildings. Recent literature demonstrates such control solutions, but they often rely on custom integrations lacking the tools and drivers needed for scalability. To address these gaps, our team has created a fully open-source software stack capable of integrating heterogeneous flexible building loads and implementing integrated portable control applications called the Open Building Operating System (OpenBOS). The software can be deployed over existing control architecture with a small capital cost. OpenBOS leverages semantic models, which have been the subject of recent investigations to facilitate application portability. The use of semantic data reduces the labor and expense required to deploy and update smart control applications, increasing scalability. In this paper, the semantic modeling schema "Brick" was used, but the proposed approach can also be applied to ASHRAE standard 223P, when released. This paper describes the methodology and software components of OpenBOS and demonstrates its functionality with a rule-based demand flexibility control application configured using a semantic model. This application was tested at a real building in NY that uses a dual-fuel heating system made up of five ductless heat pump mini-splits and a central furnace serving a single zone. The demonstration reduced electricity costs at the site by 27%, demand during a shed event by 49%, and furnace usage by 35%.

24 POWER TRANSMISSION AND DISTRIBUTION↗

A Partitioned -Task Parallel Implementation of the NASA Multiscale Analysis Tool for High Performance Computing

The NASA Multiscale Analysis Tool (NASMAT) is a “plug and play” software package that allows users to conduct massively multiscale modeling of hierarchical and nonlinear materials. This work extends the scalability and improves the High Performance Computing friendliness of NASMAT by adopting a Partitioned Task-Parallel approach. Interoperability of NASMAT with external software is enhanced through preCICE, a open source library for multiphysics coupling in a partitioned manner. Enhancement through preCICE allows for easy integration of NASMAT to other macro solvers and dissociates the parallelization strategy adopted within NASMAT from the macro solver. The task-parallel framework based on Master-Worker approach is implemented as the parallelization scheme. The scheme accounts for hierarchy of multiple scales (task-dependence) and heterogeneous nature (dynamic load balancing) of computations. The applicability and scalability of the framework will be evaluated by analyzing large scale engineering problems through massively multiscale methods.

NASMAT↗

Exascale Multiphysics Nuclear Reactor Simulations for Advanced Designs

ENRICO is a coupled application developed under the U.S. Department of Energy's Exascale Computing Project (ECP) targeting the modeling of advanced nuclear reactors. It couples radiation transport with heat and fluid simulation, including the high-fidelity, highresolution Monte-Carlo code Shift and the Computational fluid dynamics code NekRS. NekRS is a highly-performant open-source code for simulation of incompressible and low-Mach fluid flow, heat transfer, and combustion with a particular focus on turbulent flows in complex domains. It is based on rapidly convergent high-order spectral element discretizations that feature minimal numerical dissipation and dispersion. State-of-the-art multilevel preconditioners, efficient high-order time-splitting methods, and runtime-adaptive communication strategies are built on a fast OCCA-based kernel library, libParanumal, to provide scalability and portability across the spectrum of current and future high-performance computing platforms. On Frontier, Nek5000/RS has recently achieved an unprecedented milestone in breaching over 1 billion spectral elements and 350 billion degrees of freedom. Shift has demonstrated the capability to transport upwards of 1 billion particles per second in full core nuclear reactor simulations featuring complete temperature-dependent, continuous-energy physics on Frontier. Shift achieved a weak-scaling efficiency of 97.8% on 8192 nodes of Frontier and calculated 6 reactions in 214,896 fuel pin regions below 1% statistical error yielding first-of-a-kind resolution for a Monte Carlo transport application.

Hamilton, Steven P.↗

Coupling the Solar-Wind/IMF to the Ionosphere through the High Latitude Cusps

Magnetic merging is a primary means for coupling energy from the solar wind into the magnetosphere-ionosphere system. The location and nature of the process remain as open questions. By correlating measurements form diverse locations and using large-scale MHD models to put the measurements in context, it is possible to constrain out interpretations of the global and meso-scale dynamics of magnetic merging. Recent evidence demonstrates that merging often occurs at high latitudes in the vicinity of the cusps. The location is in part controlled by the clock angle in the interplanetary magnetic field (IMF) Y-Z plane. In fact, B(sub Y) bifurcated the cusp relative to source regions. The newly opened field lines may couple to the ionosphere at MLT locations of as much as 3 hr away from local noon. On the other side of noon the cusp may be connected to merging sites in the opposite hemisphere. In face, the small convection cell is generally driven by opposite hemisphere merging. B(sub X) controls the timing of the interaction and merging sites in each hemisphere, which may respond to planar features in the IMF at different times. Correlation times are variable and are controlled by the dynamics of the tilt of the interplanetary electric field phase plane. The orientation of the phase plane may change significantly on time scales of tens of minutes. Merging is temporally variable and may be occurring at multiple sites simultaneously. Accelerated electrons from the merging process excite optical signatures at the foot of the newly opened field lines. All-sky photometer observations of 557.7 nm emissions in the cusp region provide a "television picture" of the merging process and may be used to infer the temporal and spatial variability of merging, tied to variations in the IMF.

Maynard, Nelson C.↗

Radionuclide Screening Analysis and Transport Parameters for Pahute Mesa Detonations, Nevada National Security Site

A “screening analysis” is implemented at Pahute Mesa (PM), Nevada National Security Site (NNSS), to determine subsets of the 43 radionuclides listed in the radionuclide inventory of Finnegan et al. (2016) that are of relevance or potential relevance to the hydrologic source term (HST) for assessment of radionuclide transport in groundwater. Consideration is added to how levels of contamination are defined and whether any other radionuclides not in the inventory are relevant to the HST. A model is developed to estimate a range of possible source concentrations in groundwater that account for uncertainty in partitioning into melt glass and sorption into surrounding rock of the exchange zone for the 82 PM underground nuclear tests detonated in vertical shafts. The model also accounts for the varied hydrogeochemical settings. Transport parameters needed for the screening model calculations are developed from databases for hydrogeologic units, chemistry, mineralogy, fracture spacing, fracture aperture, fracture openness, matrix porosity, bulk density, saturation, and alteration. Of key importance is consideration of diagenetic zonation of the mineralogy related to hydrothermal alteration. The screening analysis compares model results with available groundwater radiochemistry data. Consistency between screening model results and available data helps validate the model for application to all 43 radionuclides, most of which have no measurements or only non-detect measurements of concentration in groundwater. The screening analysis determines that ten radionuclides – tritium, Sr-90, I-129, Cs-137, U-232, U-233, U-234, U-238, and the total of Pu-239 and Pu-240 are relevant to the HST. Determination of relevance is primarily based on data and/or model results indicating source concentration exceeding a maximum contaminant level (MCL) out to 100 years since the time of detonation. U-238 is relevant to assessment of contamination relative to MCLs for uranium and gross alpha particle activity. Five other radionuclides – C-14, Cl-36, Tc-99, Np-237, and Pu-238 – are determined potentially relevant to the HST based on a 0.1 MCL threshold. The screening analysis recommends additional attention to the natural daughter radionuclides of the uranium series (U-238) and thorium series (Th-232) decay chains, particularly Ra-226 and Ra-228, that were not included the inventory.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗

Hydrodynamic conditions in laser irradiated buried layer experiments

The calculation of open shell ionization level and radiative properties of materials in Non-Local Thermal Equilibrium (NLTE) is currently still a major challenge for any atomic model. The predictions of various NLTE atomic codes at these conditions still differ significantly. In recent years, a new buried layer platform was developed at the Lawrence Livermore National Laboratory and the Laboratory for Laser Energetics. This platform is used to measure ionization distribution and emission of open L-shell, mid-Z ions and open M-shell, high-Z ions at NLTE conditions that are relevant in many laser plasma applications. These experiments offer a unique chance for benchmarking the atomic models. In order to perform these experiments, a uniform well characterized plasma source is required. In this work, we present one-dimensional (1D) and two-dimensional simulations of the experimental platform. These simulations were used for both the design and the analysis of the experiments. The simulations demonstrate the different phases of hydrodynamic evolution of the target and identify the time windows in which uniform conditions can be achieved. A 1D expansion of the target was found to be adequate to describe the target's evolution for most of the experiment duration. The fast 1D simulations were compared with recent experimental results from the Omega laser facility. The sensitivity of the results to several modeling parameters such as the electron flux limiter and laser resonant absorption is reported.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Optimal Energy Dispatch of Distributed PVs for the Next Generation of Distribution Management Systems

Advanced Distribution Management Systems (ADMS) are being widely adopted by electric utilities for managing and optimizing the operations of their distribution systems. Distributed photovoltaic (DPV) systems with smart inverters can be controlled to adjust active power and reactive power outputs, and they are envisioned to become a part of (centrally or distributed) controllable assets managed by the ADMS for optimizing grid operations. This paper proposes an optimal energy dispatch strategy controlling DPV systems for regulating distribution voltages and achieving conservation voltage reduction. A convex optimization model is proposed with the use of linearized power flow, and the gradient projection algorithm is used to solve the optimal active power and reactive power outputs of smart inverters. The proposed optimal energy dispatch is implemented using an open-source ADMS platform, and simulation results have demonstrated the effectiveness of the proposed approach on improving distribution grid operations.

42 ENGINEERING↗

ASCR@40

One page flyers on: A Revolution in Modeling and Simulation; Applying Equations to Complex Problems; Uncertainty Quantification; Moving Big Data; Grid Computing for High-Speed Collaboration; Making Sense of Big Data; Overcoming Scaling Challenges; Building Better Computers; World-Leading Computing Facilities; Supporting Science Through Open-Source Software; and Building the Computational Workforce. Content provided by Department of Energy multiprogram laboratory researchers. Prepared by the Krell Institute for the ASCAC Subcommittee on the 40-year History of ASCR.

97 MATHEMATICS AND COMPUTING↗

Elevating SolTrace's Capabilities for the Next Generation of Concentrating Solar Analysis

SolTrace is an open-source Monte Carlo ray tracing software developed at NREL. SolTrace can characterize concentrating solar thermal (CST) collector optical performance and is CST technology agnostic. Shown in Fig. 1, SolTrace is a foundational tool in NREL's CST system and component modeling suite. SolTrace's generic surface elements can flexibly model novel collector and receiver designs to predict spatial and temporal flux distributions - critical to understand for CST component design, performance prediction, and system integration. Since its initial development, SolTrace has over 1,650 references on Google Scholar, over 9,800 downloads since 2017, and has served the CST research and development community as a benchmark of 3rd party verification. SolTrace provides users with many options for defining surface shape and boundaries. However, SolTrace provides limited documentation which can result in a steep learning curve for new users. Additionally, SolTrace lacks the computational performance required to evaluate optical performance of a CST system over the course of a year and/or iteratively over design parameters in a timely manner. To address this, we are working towards a new release of SolTrace that enables increased computational throughput by implementing ray tracing acceleration structures and enabling GPU parallelization. Additionally, we are working to improve SolTrace's usability, accessibility, and maintainability by (1) automating solar position time-dependent simulation processes, (2) creating general CST collector templates of grouped elements, (3) updating the user interface to better visualize model inputs and outputs, and (4) creating a user support network through forums, "how to" videos, and documentation.

14 SOLAR ENERGY↗

Parabolic harmonics in magnetospheric modeling - The main dipole and the ring current

A new representation is developed for that portion of the magnetospheric magnetic field which is produced by currents on the magnetopause or outside it. For a magnetopause which is a paraboloid of rotation and which confines all field lines of magnetospheric sources, the fields due to magnetopause currents are derived for two sources, a dipole and a two-parameter model of the ring current, each arbitrarily inclined to the earth-sun line. Methods are also developed for accomodating changes in the size of the magnetospheric cavity and for using nonparaboloidal boundaries, open at the nightside. Neutral points and polar cusps are derived, and in agreement with observations it is found that a growth of the ring current shifts them equatorward. By extrapolation from these results, it is argued that the observation of such shifts during magnetic storms may help resolve the relative importance of two processes proposed for magnetic storms.

Stern, D. P.↗