Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “software differences”

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,315 records · Page 73

Tycho: Realistically Simulating Exoplanets within Stellar Clusters. I. Improving the Monte Carlo Approach

To fully understand the diverse population of exoplanets, we must study their early lives within open clusters, the birthplace of most stars with masses >0.5M {sub ⊙} (including those currently in the field). Indeed, when we observe planets within clustered environments, we notice highly eccentric and odd systems that suggest the importance of dynamical pathways created by interactions with additional bodies (as in the case of HD 285507b). However, it has proven difficult to investigate these effects, as many current numerical solvers for the multi-scale N-body problem are simplified and limited in scope. To remedy this, we aim to create a physically complete computational solution to explore the role of stellar close encounters and interplanetary interactions in producing the observed exoplanet populations for both open cluster stars and field stars. We present a new code, Tycho, which employs a variety of different computational techniques, including multiple N-body integration methods, close encounter handling, modified Monte Carlo scattering experiments, and a variety of empirically informed initial conditions. We discuss the methodology in detail, and its implementation within the AMUSE software framework. Approximately 1% of our systems are promptly disrupted by star-star encounters contributing to the rogue planets occurrence rate. Additionally, we find that close encounters which that perturb long-period planets lead to 38.3% of solar-system-like planetary systems becoming long-term unstable.

79 ASTRONOMY AND ASTROPHYSICS↗

Towards Common Components and a Standard Model for Workflow Systems

The role of scalable high-performance workflows and flexible workflow management systems that can support multiple simulations will continue to increase in importance. For example, with the end of Dennard scaling, there is a need to substitute a single long running simulation with multiple repeats of shorter simulations, or concurrent replicas. Further, many scientific problems involve ensembles of simulations in order to solve a higher-level problem or produce statistically meaningful results. However most supercomputing software development and performance enhancements have focused on optimizing single- simulation performance. On the other hand, there is a strong inconsistency in the definition and practice of workflows and workflow management systems. This inconsistency often centers around the difference between several different types of workfiows, including modeling and simulation, grid, uncertainty quantification, and purely conceptual workflows. This work explores this phenomenon by examining the different types of workflows and workflow management systems, reviewing the perspective of a large supercomputing facility, examining the common features and problems of workflow management systems, and finally presenting a proposed solution based on the concept of common building blocks. The implications of the continuing proliferation of workflow management systems and the lack of interoperability between these systems are discussed from a practical perspective. In doing so, we have begun an investigation of the design and implementation of open workflow systems for supercomputers based upon common components.

Billings, Jay Jay↗

Performance Evaluation of Distributed Energy Resource Management Systems Using Software- and Hardware-in-the-Loop

Distributed energy resource management systems (DERMS) are an important emerging technology for coordinating distributed energy resources (DERs) integrated in power distribution networks. These solutions implement a software platform to monitor, control, and coordinate DERs in real time to provide aggregate grid services to the electric utility or distribution system operator. To move DERMS technology forward, research and performance validation in a laboratory environment are needed to ensure proper function of DERMS prior to field deployment. This presentation presents a generic hardware-in-the-loop (HIL) platform to evaluate performance of DERMS using software-in-the-loop (SIL) and power HIL. The platform provides realistic testing in a laboratory environment, including accurate modeling of full-scale distribution systems, DERMS software controllers, and at-power distributed energy resources (DERs), such as solar and battery inverters. Two case studies are presented to demonstrate how the platform can be used independently as an integrated development and validation approach for different projects. Laboratory testing results are included to showcase the voltage regulation performance of each study and the benefits of implementing such a platform for DERMS evaluation.

DERMS↗

Using a Large Language Model for Accurate Technical Language Generation in the Predictive Maintenance of Circulating Water Systems in Nuclear Power Plants

Machine learning (ML) methods for predictive maintenance (PdM) are emerging as effective proactive strategies for diagnosing equipment degradation and enabling effective decision-making. However, explainability and trustworthiness of artificial intelligence are two salient challenges that need to be addressed for wider deployment of these technologies in nuclear power plants (NPPs). Large language models (LLMs) offer a unique approach to tackle these challenges by explaining PdM, work orders, diagnosis results, and ML algorithms to users, who may not be familiar with ML and PdM in general. Moreover, by dynamically retrieving relevant information from technical documents and evaluating factuality of LLM generation, the accuracy and relevance of LLM generations can be improved. This work demonstrates using LLMs to explain the causes and consequences of circulating water system failures based on multiyear NPP work orders. This work tests the capability of multimodal LLM approaches in explaining the differences in the circulating water system from both the Salem and Hope Creek NPPs using both text and image resources. This work also demonstrates the use of multimodal LLMs in describing the diagnosis tab of a predictive maintenance software named VIsualization for PrEdictive maintenance Recommendation (VIPER) to users.

22 - GENERAL STUDIES OF NUCLEAR REACTORS↗

Study of 201 Non-Small Cell Lung Cancer Patients Given Stereotactic Ablative Radiation Therapy Shows Local Control Dependence on Dose Calculation Algorithm

Purpose: Pencil beam (PB) and collapsed cone convolution (CCC) dose calculation algorithms differ significantly when used in the thorax. However, such differences have seldom been previously directly correlated with outcomes of lung stereotactic ablative body radiation (SABR). Methods and Materials: Data for 201 non-small cell lung cancer patients treated with SABR were analyzed retrospectively. All patients were treated with 50 Gy in 5 fractions of 10 Gy each. The radiation prescription mandated that 95% of the planning target volume (PTV) receive the prescribed dose. One hundred sixteen patients were planned with BrainLab treatment planning software (TPS) with the PB algorithm and treated on a Novalis unit. The other 85 were planned on the Pinnacle TPS with the CCC algorithm and treated on a Varian linac. Treatment planning objectives were numerically identical for both groups. The median follow-up times were 24 and 17 months for the PB and CCC groups, respectively. The primary endpoint was local/marginal control of the irradiated lesion. Gray's competing risk method was used to determine the statistical differences in local/marginal control rates between the PB and CCC groups. Results: Twenty-five patients planned with PB and 4 patients planned with the CCC algorithms to the same nominal doses experienced local recurrence. There was a statistically significant difference in recurrence rates between the PB and CCC groups (hazard ratio 3.4 [95% confidence interval: 1.18-9.83], Gray's test P=.019). The differences (Δ) between the 2 algorithms for target coverage were as follows: ΔD99{sub GITV} = 7.4 Gy, ΔD99{sub PTV} = 10.4 Gy, ΔV90{sub GITV} = 13.7%, ΔV90{sub PTV} = 37.6%, ΔD95{sub PTV} = 9.8 Gy, and ΔD{sub ISO} = 3.4 Gy. GITV = gross internal tumor volume. Conclusions: Local control in patients receiving who were planned to the same nominal dose with PB and CCC algorithms were statistically significantly different. Possible alternative explanations are described in the report, although they are not thought likely to explain the difference. We conclude that the difference is due to relative dosimetric underdosing of tumors with the PB algorithm.

62 RADIOLOGY AND NUCLEAR MEDICINE↗

MixPI: Mixed-time slicing path integral software for quantized molecular dynamics simulations

We introduce the MixPI software to implement path integral molecular dynamics (PIMD) simulations for the study of condensed phase systems where nuclear quantum effects (NQEs) are important. In contrast to existing PIMD simulation software, MixPI enables the implementation of mixed quantum–classical path integral simulations where only a subset of system degrees of freedom (dofs) are treated quantum mechanically in an extended phase space while the remaining dofs are described classically. We expect this software to be particularly useful for simulations of electron and proton transfer in condensed phase systems, as well as for the study of biological and material systems where only a handful of dofs contribute significantly to the observed NQEs. We demonstrate the use of MixPI in two different systems. The first is a simple water model where we implement a set of mixed quantum–classical simulations to compute average energy and radial distribution functions. We use these simulations to benchmark the effectiveness of MixPI and to demonstrate how it enables systematic investigation into the origin of observed NQEs. We then compute radial distribution functions for a system where MixPI is essential: a solvated metal (M 2+ ) cation described using an explicit quantized electron localized on an M 3+ ion in water.

chemical physics↗

Temperature dependent simulation of diamond depleted Schottky PIN diodes

Diamond is considered as an ideal material for high field and high power devices due to its high breakdown field, high lightly doped carrier mobility, and high thermal conductivity. The modeling and simulation of diamond devices are therefore important to predict the performances of diamond based devices. In this context, we use Silvaco{sup ®} Atlas, a drift-diffusion based commercial software, to model diamond based power devices. The models used in Atlas were modified to account for both variable range and nearest neighbor hopping transport in the impurity bands associated with high activation energies for boron doped and phosphorus doped diamond. The models were fit to experimentally reported resistivity data over a wide range of doping concentrations and temperatures. We compare to recent data on depleted diamond Schottky PIN diodes demonstrating low turn-on voltages and high reverse breakdown voltages, which could be useful for high power rectifying applications due to the low turn-on voltage enabling high forward current densities. Three dimensional simulations of the depleted Schottky PIN diamond devices were performed and the results are verified with experimental data at different operating temperatures.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Exploring the effects of data quality, data worth, and redundancy of CO2 gas pressure and saturation data on reservoir characterization through PEST Inversion

This study examined the impacts of reservoir properties on CO2 migration after subsurface injection and evaluated the possibility of characterizing reservoir properties using CO2 monitoring data such as saturation distribution. The injection reservoir was assumed to be located 1400-1500 m below the ground surface such that CO2 remained in the supercritical state. The reservoir was assumed to contain layers with alternating conductive and resistive properties, which is analogous to actual geological formations such as the Mount Simon Sandstone unit. The CO2 injection simulation used a cylindrical grid setting in which the injection well was situated at the center of the domain, which extended up to 8000 m from the injection well. The CO2 migration was simulated using the PNNL-developed simulator STOMP-CO2e (the water-salt-CO2 module). We adopted a nonlinear parameter estimation and optimization modeling software package, PEST, for automated reservoir parameter estimation. We explored the effects of data quality, data worth, and data redundancy on the detectability of reservoir parameters using CO2 saturation monitoring data, by comparing PEST inversion results using data with different levels of noises, various numbers of monitoring wells and locations, and different data collection spacing and temporal sampling intervals. This study yielded insight into the use of CO2 saturation monitoring data for reservoir characterization and how to design the monitoring system to optimize data worth and reduce data redundancy.

Carbon sequestration↗

Dynamic calibration of differential equations using machine learning, with application to turbulence models

We present a methodology for calibration of parametric ordinary and partial differential equation models, using off-the-shelf software for back-propagation in Neural Networks (NN). As a prototypical example, we consider calibration of a Reynolds-averaged Navier-Stokes (RANS) turbulence closure model, against ground truth data from direct numerical simulations (DNS) of two different turbulent flows. Numerical time integration is represented as a custom NN, where only the RANS model parameters are trainable. A loss function is defined to quantify the mismatch between the NN prediction and the ground truth over a predefined, finite time integration window. This loss function is then minimized using a gradient descent method utilizing the back-propagation algorithm. Furthermore, this dynamic approach to training is to be contrasted with a static approach, wherein a least square regression estimate for parameters is obtained in the limit of an infinitesimal time integration window. In a first test of static and dynamic approaches against ground truth data generated by the model, the former proves to be significantly faster and more accurate than the latter at recovering the parameters. When both calibration approaches are tested against DNS data, for which it is known that the model cannot achieve a perfect fit, the static approach yields a good prediction only for short times, while the dynamic approach results in physical and stable predictions over the entire integration window. After optimization of the dynamic approach for time step, spatial resolution, stability, and physics-based constraints, we obtain a 50% improvement of outcomes over those obtained from the existing, manually calibrated set of parameters, demonstrating the merits of this systematic and automated procedure.

97 MATHEMATICS AND COMPUTING↗

Implicit solvation model for density-functional study of nanocrystal surfaces and reaction pathways

Solid-liquid interfaces are at the heart of many modern-day technologies and provide a challenge to many materials simulation methods. A realistic first-principles computational study of such systems entails the inclusion of solvent effects. In this work, we implement an implicit solvation model that has a firm theoretical foundation into the widely used density-functional code Vienna ab initio Software Package. The implicit solvation model follows the framework of joint density functional theory. We describe the framework, our algorithm and implementation, and benchmarks for small molecular systems. We apply the solvation model to study the surface energies of different facets of semiconducting and metallic nanocrystals and the S N 2 reaction pathway. We find that solvation reduces the surface energies of the nanocrystals, especially for the semiconducting ones and increases the energy barrier of the S N 2 reaction.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Towards scanning nanostructure X-ray microscopy

This article demonstrates spatial mapping of the local and nanoscale structure of thin film objects using spatially resolved pair distribution function (PDF) analysis of synchrotron X-ray diffraction data. This is exemplified in a lab-on-chip combinatorial array of sample spots containing catalytically interesting nanoparticles deposited from liquid precursors using an ink-jet liquid-handling system. A software implementation is presented of the whole protocol, including an approach for automated data acquisition and analysis using the atomic PDF method. The protocol software can handle semi-automated data reduction, normalization and modeling, with user-defined recipes generating a comprehensive collection of metadata and analysis results. By slicing the collection using included functions, it is possible to build images of different contrast features chosen by the user, giving insights into different aspects of the local structure.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

MO-C-17A-05: A Three-Dimensional Head-And-Neck Phantom for Validation of Kilovoltage- and Megavoltage-Based Deformable Image Registration

Purpose: To develop a three-dimensional (3D) deformable head-and-neck (H and N) phantom with realistic tissue contrast for both kilovoltage and megavoltage computed tomography and use it to objectively evaluate deformable image registration (DIR) algorithms. Methods: The phantom represents H and N patient anatomy. It is constructed from thermoplastic, which becomes pliable in boiling water, and hardened epoxy resin. Using a system of additives, the Hounsfield unit (HU) values of these materials were tuned to mimic anatomy for both kilovoltage (kV) and megavoltage (MV) imaging. The phantom opened along a sagittal midsection to reveal nonradiopaque markers, which were used to characterize the phantom deformation. The deformed and undeformed phantom was scanned with kV and MV computed tomography. Additionally, a calibration curve was created to change the HUs of the MV scans to be similar to kV HUs, (MC). The extracted ground-truth deformation was then compared to the results of two commercially available DIR algorithms, from Velocity Medical Solutions and MIM Software. Results: The phantom produced a 3D deformation, representing neck flexion, with a magnitude of up to 8 mm and was able represent tissue HUs for both kV and MV imaging modalities. The two tested deformation algorithms yielded vastly different results. For kV-kV registration, MIM made the lowest mean error, and Velocity made the lowest maximum error. For MV-MV, kV-MV, and kV-MC Velocity produced both the lowest mean and lowest maximum errors. Conclusion: The application of DIR across different imaging modalities is particularly difficult, due to differences in tissue HUs and the presence of imaging artifacts. For this reason, DIR algorithms must be validated specifically for this purpose. The developed H and N phantom is an effective tool for this purpose.

62 RADIOLOGY AND NUCLEAR MEDICINE↗

Qualitative Analysis, Structure-Preserving Discretizations and Scalable Solvers for Multiphysics Problems. Final Technical Report

Many problems in physics, engineering, meteorology, computer science, and other fields involve a certain level of coupling between different physical quantities. With the development of advanced computational modeling, simulation, and hardware, simulations for many important multiphysics/multiscale coupling effects are now more tractable. In this project, we have developed robust, efficient, and practical numerical algorithms for such simulation software that are capable of tackling these complex mathematical models in many fundamental ways such as stability, wellposedness, and consistency. We expect that such an integrated theoretical and algorithmic framework will have a substantial impact on numerous multiscale and multiphysics modeling and simulations.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

The Dark Energy Spectrographic Instrument (DESI) Guide-Focus Assembly

The specific research of the DESI project is to study in detail the expansion history of the Universe over the past 10 billion years. In order to do this, the project designed, fabricated, tested and commissioned the DESI instrument which has been deployed at the Mayall Telescope at the Kitt Peak National Observatory near Tucson, Arizona. The DESI project was conducted by approximately 30 US and foreign national laboratories. Our CRADA with the The Institut de Física d’Altes Energies (IFAE) in Barcelona (Spain) involved both research planning for eventual observing on the DESI telescope; and in providing hardware jointly designed by LBNL and IFAE; ten guide and focus sensor systems (GFAs) that either provide feedback to the Mayall telescope tracking system or provide focus position information to the DESI corrector hexapod. The key major goals of this CRADA were achieved.For survey planning, Drs. Fernandez, Fosalba, Garcia-Bellido, and Gaztañaga were all important contributors to the working groups to which they were assigned. Their application of redshift techniques, development of simulations to support target selection, exploration of the use of alternative observables to extract cosmological information like AP and void statistics, and the use of PAU data to assist in the optimization of targeting selections have been vital in the early planning for targets during commissioning. The major role in the project involved Drs. Castander, Miquel, Padilla, and Sànchez. They were all directly involved with the development of the Guide-Focus Assembly (GFA) hardware, software and filter validation. In order to accomplish the goals of the DESI Project, the Mayall telescope must track DESI targeted galaxies for up to 15 minutes and the DESI corrector must be repositioned on its hexapod system for different telescope orientations. The guide and focus assemblies (GFA) provide in this information. There are ten GFA modules on the focal plane, six providing tracking information and four providing focus information. A common GFA design supports both functions by a simple switching of its operating mode. IFAE designed, fabricated, tested and delivered twelve GFAs. Once received at LBNL, the GFAs were integrated with the focal plane. This was an iterative process as the GFA’s were delivered as they were made. Our working relationship with the IFAE scientists and engineers was remarkable given the distance and time zone issues. The GFAs have performed extremely well during the commissioning phase of the instrument.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE WPTO FY2021 Seedling End Report: Coupling Soil-Foundation Models to MHK Device Dynamic Models to Tighten the Design Envelope

The purpose of this Seedling project is to couple a marine renewable energy (MRE) dynamics simulation software with the soil-foundation models in the OC6 Phase II project [Bergua et al., 2021] and evaluate the software’s performance. This is a first step to accurately evaluating soil-foundation impacts on other types of MRE, like wave or current energy converters (WECs, CECs). OC6 Phase II compares offshore wind turbine (OWT) simulations using several different soil-foundation models to identify and fill key gaps in soil-foundation analyses. WEC-Sim was chosen to model the OC6 Phase II offshore wind turbine and various load cases because of its adaptability, accuracy of hydrodynamic loads, and ability to apply an arbitrary wind loading. Of the four methods used in OC6, the apparent fixity soil-foundation method was coupled with WEC-Sim. Technical challenges with flexible hydrodynamic bodies, added mass and external function libraries inhibited the ability to compare the WEC-Sim results to other OC6 participants. These challenges required that the WEC-Sim model of the OC6 OWT use a combination of rigid and flexible bodies to ensure a numerically stable solution. The rigid monopile creates a more stiff system and causes smaller amplitude motion under hydrodynamic loading and higher dominant frequency of motion under wind loading. These discrepancies are expected based on the increased stiffness of the WEC-Sim case.

16 TIDAL AND WAVE POWER↗

From DNA to FBA: How to Build Your Own Genome-Scale Metabolic Model

Microbiological studies are increasingly relying on in silico methods to perform exploration and rapid analysis of genomic data, and functional genomics studies are supplemented by the new perspectives that genome-scale metabolic models offer. A mathematical model consisting of a microbe’s entire metabolic map can be rapidly determined from whole-genome sequencing and annotating the genomic material encoded in its DNA. Flux-balance analysis (FBA), a linear programming technique that uses metabolic models to predict the phenotypic responses imposed by environmental elements and factors, is the leading method to simulate and manipulate cellular growth in silico. However, the process of creating an accurate model to use in FBA consists of a series of steps involving a multitude of connections between bioinformatics databases, enzyme resources, and metabolic pathways. We present the methodology and procedure to obtain a metabolic model using PyFBA, an extensible Pythonbased open-source software package aimed to provide a platform where functional annotations are used to build metabolic models (http://linsalrob.github.io/PyFBA). Backed by the Model SEED biochemistry database, PyFBA contains methods to reconstruct a microbe’s metabolic map, run FBA upon different media conditions, and gap-fill its metabolism. The extensibility of PyFBA facilitates novel techniques in creating accurate genome-scale metabolic models.

genome annotation↗

Quantifying the Opportunity Limits of Automatic Residential Electric Load Shaping

Electric utility residential demand response programs typically reduce load a few times a year during periods of peak energy use. In the future, utilities and consumers may monetarily and environmentally benefit from continuously shaping load by alternatively encouraging or discouraging the use of electricity. One way to shape load and introduce elasticity is to broadcast forecasts of dynamic electricity prices that orchestrate electricity supply and demand in order to maximize the efficiency of conventional generation and the use of renewable resources including wind and solar energy. A binary control algorithm that influences the on and off states of end uses was developed and applied to empirical time series data to estimate price-based instantaneous opportunities for shedding and adding electric load. To overcome the limitations of traditional stochastic methods in quantifying diverse, non-Gaussian, non-stationary distributions of observed appliance behaviour, recent developments in wavelet-based analysis were applied to capture and simulate time-frequency domain behaviour. The performance of autoregressive and spectral reconstruction methods was compared, with phase reconstruction providing the best simulation ensembles. Results show spatiotemporal differences in the amount of load that can be shed and added, which suggest further investigation is warranted in estimating the benefits anticipated from the wide-scale deployment of continuous automatic residential load shaping. Empirical data and documented software code are included to assist in reproducing and extending this work.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗

Quantitative Uncertainty Analysis for Diffraction

Current software uses a least squares statistical method. For each model parameter, current software yields a point estimate and estimated standard error for each point estimate. In contrast, we have employed a Bayesian approach with sampling algorithm to provide a full posterior probability density function on each model parameter. This yields a different type of solution (probabilistic) with better fidelity. The Quantitative Uncertainty Analysis for Diffraction (QUAD) package returns a matrix of affect model parameters as a function of iteration through the MCMC analysis. QUAD is a Python code that allows for analysis of X-ray and neutron diffraction data to infer the structure of materials with quantifiable certainty. The code uses Bayesian statistics and Markov chain sampling algorithms to create posterior probability distributions on all material structure parameters that are modeled by a researcher. The statistical packages use the open source GSASII project as a library to calculate model diffraction patterns.

Fancher, ChristopherM↗