Search NASA⌕ Search

SEARCH · Search NASA

Results for “model structural error”

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,225 records · Page 68
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Multiscale and Machine Learning Modeling for Process-informed Microstructure Prediction in Additively Manufactured Materials Using MALAMUTE

Advanced Materials and Manufacturing Technologies (AMMT) program under the Department of Energy Office of Nuclear Energy, aims to develop and qualify additively-manufactured materials for nuclear applications. The key challenges to these efforts are the microstructural variabilities observed on the AM products and their impact on the properties and performance of the material in extreme environments. AMMT is using a combination of high-through-put experimental and modeling techniques to accelerate the qualification efforts. Conventionally, in-situ and ex-situ characterizations and testing are performed to correlate different aspects of the AM process to the final product and its performance. However, adopting a trial-and-error approach to experimentally evaluate the vast range of process parameters required to capture the microstructural variabilities is cost-prohibitive. Modeling and simulation provide a comparatively inexpensive way to understand and correlate the microstructural evolution to the processing conditions. The modeling and simulation work-packages within the AMMT program aims to use physics-based and machine learning modeling capabilities to develop a digital twin for AM that can correlate the process conditions to the final product and establish a process-structure-property-performance (PSPP) correlation for AM materials. The melting and subsequent solidification that occurs during the AM process is a complex phenomenon that requires multiscale multiphysics analysis. Idaho National Laboratory’s (INL) Multiphysics Object-Oriented Simulation Environment (MOOSE), specifically the MOOSE Application Library for Advanced Manufacturing UTilitiEs (MALAMUTE) software, provides an ideal platform for developing the multiphysics multiscale model to explore the intricacies of the microstructural evolution during the AM processes within a single framework. Furthermore, given that such full-fidelity simulations can be computationally intensive, reduced order models are necessary to explore the PSPP space for AM materials in an efficient, reliable, and cost-effective way. This work package focuses on understanding the role of process variabilities on the various microstructural characteristics of the AM materials. Microstructures unique to AM materials, such as compositional micro-heterogeneity and dislocation cells, are of particular interest here since they can influence the creep properties and radiation performance. In fiscal year (FY) 24, we significantly advanced upon our work in the last fiscal year, both on physics-based and ML models. The alloy solidification model available in MOOSE has been extended to incorporate the thermodynamic properties and free energy relevant to 316SS. The model demonstrates the Cr segregation that occurs during solidifcation. It is demonstrated that rate of solidification and solute segregation is primarily influence by the cooling rate dictating the level of freezing. This work captures the microstructural variabilities at the subgrain level that are often missing in the part-scale models. With an aim to connect the microstructural evolution model to realistic process conditions, a reduced order model is developed for predicting the thermal conditions around meltpool from high-fidelity process simulations. Furthermore, machine learning approach is used to accelerate the temperature prediction during the AM process. In the following years, MALAMUTE will be used to connect different aspects of the models and quantitatively predict the microstructural evolution. The developed ML-based surrogate model will consider the process conditions as the input to predict the microstructural features in a cost-effective way. The generated microstructures can be used by other work packages under AMMT to evaluate the properties and environmental response of the material at the mesoscale. Thus, this work help identify the key microstructural features at the subgrain level that are significant in property/performance prediction of the AM products. This work will provide inputs to the large-scale process variability models to reevaluate and validate assumptions/simplifications made in the part-scale models. Furthermore, through active learning this work will help identify the data need from both modeling and experimental sides for development of a robust digital twin for AM.

36 MATERIALS SCIENCE↗
DOE OSTI2016

The CMIP6 Sea-Ice Model Intercomparison Project (SIMIP): Understanding sea ice through climate-model simulations

A better understanding of the role of sea ice for the changing climate of our planet is the central aim of the diagnostic Coupled Model Intercomparison Project 6 (CMIP6)-endorsed Sea-Ice Model Intercomparison Project (SIMIP). To reach this aim, SIMIP requests sea-ice-related variables from climate-model simulations that allow for a better understanding and, ultimately, improvement of biases and errors in sea-ice simulations with large-scale climate models. This then allows us to better understand to what degree CMIP6 model simulations relate to reality, thus improving our confidence in answering sea-ice-related questions based on these simulations. Furthermore, the SIMIP protocol provides a standard for sea-ice model output that will streamline and hence simplify the analysis of the simulated sea-ice evolution in research projects independent of CMIP. To reach its aims, SIMIP provides a structured list of model output that allows for an examination of the three main budgets that govern the evolution of sea ice, namely the heat budget, the momentum budget, and the mass budget. Furthermore, we explain the aims of SIMIP in more detail and outline how its design allows us to answer some of the most pressing questions that sea ice still poses to the international climate-research community.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2019

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↗
NASA NTRS1984

Structure and evolution of the large scale solar and heliospheric magnetic fields

Structure and evolution of large scale photospheric and coronal magnetic fields in the interval 1976-1983 were studied using observations from the Stanford Solar Observatory and a potential field model. The solar wind in the heliosphere is organized into large regions in which the magnetic field has a componenet either toward or away from the sun. The model predicts the location of the current sheet separating these regions. Near solar minimum, in 1976, the current sheet lay within a few degrees of the solar equator having two extensions north and south of the equator. Soon after minimum the latitudinal extent began to increase. The sheet reached to at least 50 deg from 1978 through 1983. The complex structure near maximum occasionally included multiple current sheets. Large scale structures persist for up to two years during the entire interval. To minimize errors in determining the structure of the heliospheric field particular attention was paid to decreasing the distorting effects of rapid field evolution, finding the optimum source surface radius, determining the correction to the sun's polar field, and handling missing data. The predicted structure agrees with direct interplanetary field measurements taken near the ecliptic and with coronameter and interplanetary scintillation measurements which infer the three dimensional interplanetary magnetic structure. During most of the solar cycle the heliospheric field cannot be adequately described as a dipole.

Hoeksema, J. T.↗
NASA NTRS1991

Finite element structural redesign by large admissible perturbations

In structural redesign, two structural states are involved; the baseline (known) State S1 with unacceptable performance, and the objective (unknown) State S2 with given performance specifications. The difference between the two states in performance and design variables may be as high as 100 percent or more depending on the scale of the structure. A Perturbation Approach to Redesign (PAR) is presented to relate any two structural states S1 and S2 that are modeled by the same finite element model and represented by different values of the design variables. General perturbation equations are derived expressing implicitly the natural frequencies, dynamic modes, static deflections, static stresses, Euler buckling loads, and buckling modes of the objective S2 in terms of its performance specifications, and S1 data and Finite Element Analysis (FEA) results. Large Admissible Perturbation (LEAP) algorithms are implemented in code RESTRUCT to define the objective S2 incrementally without trial and error by postprocessing FEA results of S1 with no additional FEAs. Systematic numerical applications in redesign of a 10 element 48 degree of freedom (dof) beam, a 104 element 192 dof offshore tower, a 64 element 216 dof plate, and a 144 element 896 dof cylindrical shell show the accuracy, efficiency, and potential of PAR to find an objective state that may differ 100 percent from the baseline design.

Bernitsas, Michael M.↗
NASA NTRS1994

General purpose architecture for intelligent computer-aided training

An intelligent computer-aided training system having a general modular architecture is provided for use in a wide variety of training tasks and environments. It is comprised of a user interface which permits the trainee to access the same information available in the task environment and serves as a means for the trainee to assert actions to the system; a domain expert which is sufficiently intelligent to use the same information available to the trainee and carry out the task assigned to the trainee; a training session manager for examining the assertions made by the domain expert and by the trainee for evaluating such trainee assertions and providing guidance to the trainee which are appropriate to his acquired skill level; a trainee model which contains a history of the trainee interactions with the system together with summary evaluative data; an intelligent training scenario generator for designing increasingly complex training exercises based on the current skill level contained in the trainee model and on any weaknesses or deficiencies that the trainee has exhibited in previous interactions; and a blackboard that provides a common fact base for communication between the other components of the system. Preferably, the domain expert contains a list of 'mal-rules' which typifies errors that are usually made by novice trainees. Also preferably, the training session manager comprises an intelligent error detection means and an intelligent error handling means. The present invention utilizes a rule-based language having a control structure whereby a specific message passing protocol is utilized with respect to tasks which are procedural or step-by-step in structure. The rules can be activated by the trainee in any order to reach the solution by any valid or correct path.

Loftin, R. Bowen↗
NASA NTRS2004

Toward a Global Map of Raindrop Size Distributions: Rain-Type Classification and Its Implications for Validating Global Rainfall Products - Part 1

Variability in the global distribution of precipitation is recognized as a key element in assessing the impact of climate change for life on earth. The response of precipitation to climate forcings is, however, poorly understood because of discrepancies in the magnitude and sign of climatic trends in satellite-based rainfall estimates. Quantifying and ultimately removing these biases is critical for studying the response of the hydrologic cycle to climate change. In addition, estimates of random errors owing to variability in algorithm assumptions on local spatial and temporal scales are critical for establishing how strongly their products should be weighted in data assimilation or model validation applications and for assigning a level of confidence to climate trends diagnosed from the data. This paper explores the potential for refining assumed drop size distributions (DSDs) in global radar rainfall algorithms by establishing a link between satellite observables and information gleaned from regional validation experiments where polarimetric radar, Doppler radar, and disdrometer measurements can be used to infer raindrop size distributions. By virtue of the limited information available in the satellite retrieval framework, the current method deviates from approaches adopted in the ground-based radar community that attempt to relate microphysical processes and resultant DSDs to local meteorological conditions. Instead, the technique exploits the fact that different microphysical pathways for rainfall production are likely to lead to differences in both the DSD of the resulting raindrops and the three-dimensional structure of associated radar reflectivity profiles. Objective rain-type classification based on the complete three-dimensional structure of observed reflectivity profiles is found to partially mitigate random and systematic errors in DSDs implied by differential reflectivity measurements. In particular, it is shown that vertical and horizontal reflectivity structure obtained from spaceborne radar can be used to reproduce significant differences in Z(sub dr) between the easterly and westerly climate regimes observed in the Tropical Rainfall Measuring Mission Large-scale Biosphere-Atmosphere (TRMM-LBA) field experiment as well as the even larger differences between Amazonian rainfall and that observed in eastern Colorado. As such, the technique offers a potential methodology for placing locally observed DSD information into a global framework.

L'Ecuyer, Tristan S.↗
DOE OSTI2017

Development of performance specifications for hybrid modeling of floating wind turbines in wave basin tests

Hybrid modeling—combining physical testing and numerical simulation in real time$-$opens new opportunities in floating wind turbine research. Wave basin testing is an important validation step for floating support structure design, but the conventional approaches that use physical wind above the basin are limited by scaling problems in the aerodynamics. Applying wind turbine loads with an actuation system that is controlled by a simulation responding to the basin test in real time offers a way to avoid scaling problems and reduce cost barriers for floating wind turbine design validation in realistic coupled wind and wave conditions. This paper demonstrates the development of performance specifications for a system that couples a wave basin experiment with a wind turbine simulation. Two different points for the hybrid coupling are considered: the tower-base interface and the aero-rotor interface (the boundary between aerodynamics and the rotor structure). Analyzing simulations of three floating wind turbine designs across seven load cases reveals the motion and force requirements of the coupling system. By simulating errors in the hybrid coupling system, the sensitivity of the floating wind turbine response to coupling quality can be quantified. The sensitivity results can then be used to determine tolerances for motion tracking errors, force actuation errors, bandwidth limitations, and latency in the hybrid coupling system. These tolerances can guide the design of hybrid coupling systems to achieve desired levels of accuracy. An example demonstrates how the developed methods can be used to generate performance specifications for a system at 1:50 scale. Results show that sensitivities vary significantly between support structure designs and that coupling at the aero-rotor interface has less stringent requirements than those for coupling at the tower base. As a result, the methods and results presented here can inform design of future hybrid coupling systems and enhance understanding of how test results are affected by hybrid coupling quality.

16 TIDAL AND WAVE POWER↗
NASA NTRS2011

Effective Tree Scattering at L-Band

For routine microwave Soil Moisture (SM) retrieval through vegetation, the tau-omega [1] model [zero-order Radiative Transfer (RT) solution] is attractive due to its simplicity and eases of inversion and implementation. It is the model used in baseline retrieval algorithms for several planned microwave space missions, such as ESA's Soil Moisture Ocean Salinity (SMOS) mission (launched November 2009) and NASA's Soil Moisture Active Passive (SMAP) mission (to be launched 2014/2015) [2 and 3]. These approaches are adapted for vegetated landscapes with effective vegetation parameters tau and omega by fitting experimental data or simulation outputs of a multiple scattering model [4-7]. The model has been validated over grasslands, agricultural crops, and generally light to moderate vegetation. As the density of vegetation increases, sensitivity to the underlying SM begins to degrade significantly and errors in the retrieved SM increase accordingly. The zero-order model also loses its validity when dense vegetation (i.e. forest, mature corn, etc.) includes scatterers, such as branches and trunks (or stalks in the case of corn), which are large with respect to the wavelength. The tau-omega model (when applied over moderately to densely vegetated landscapes) will need modification (in terms of form or effective parameterization) to enable accurate characterization of vegetation parameters with respect to specific tree types, anisotropic canopy structure, presence of leaves and/or understory. More scattering terms (at least up to first-order at L-band) should be included in the RT solutions for forest canopies [8]. Although not really suitable to forests, a zero-order tau-omega model might be applied to such vegetation canopies with large scatterers, but that equivalent or effective parameters would have to be used [4]. This requires that the effective values (vegetation opacity and single scattering albedo) need to be evaluated (compared) with theoretical definitions of these parameters. In a recent study [9], effective vegetation opacity of coniferous trees was compared with two independent estimates of the same parameter. First, a zero-order RT model was fitted to multiangular microwave emissivity data in a least-square sense to provide effective vegetation optical depth as done in spaceborne retrieval algorithms. Second, a ratio between radar backscatter measurements with a corner reflector under trees and in an open area was calculated to obtain measured tree propagation characteristics. Finally, the theoretical propagation constant was determined by forward scattering theorem using detailed measurements of size/angle distributions and dielectric constants of the tree constituents (trunk, branches, and needles). Results indicated that the effective attenuation values are smaller than but of similar magnitude to both the theoretical and measured values. This study will complement the previous work [9] and will focus on characterization of effective scattering albedo by assuming that effective vegetation opacity is same as theoretical opacity. The resultant effective albedo will not be the albedo of single forest canopy element anymore, but it becomes a global parameter, which depends on all the processes taking place within the canopy including multiple scattering as described.

Kurum, Mehmet↗
DOE OSTI2024

Multi-resolution partial differential equations preserved learning framework for spatiotemporal dynamics

Traditional data-driven deep learning models often struggle with high training costs, error accumulation, and poor generalizability in complex physical processes. Physics-informed deep learning (PiDL) addresses these challenges by incorporating physical principles into the model. Most PiDL approaches regularize training by embedding governing equations into the loss function, yet this depends heavily on extensive hyperparameter tuning to weigh each loss term. To this end, we propose to leverage physics prior knowledge by “baking” the discretized governing equations into the neural network architecture via the connection between the partial differential equations (PDE) operators and network structures, resulting in a PDE-preserved neural network (PPNN). This method, embedding discretized PDEs through convolutional residual networks in a multi-resolution setting, largely improves the generalizability and long-term prediction accuracy, outperforming conventional black-box models. The effectiveness and merit of the proposed methods have been demonstrated across various spatiotemporal dynamical systems governed by spatiotemporal PDEs, including reaction-diffusion, Burgers’, and Navier-Stokes equations.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1995

(abstract) Characterization of Structural Response for Systems with Loose Joints

This paper describes a technique to locate and characterize loose joints in a large truss structure to generate an accurate structural model. The joint looseness is modeled as a gap in the member that opens and closes depending on the loading. Arbitrarily placed actuators are used to prestress the structure to first linearize the response. Next the actuator displacements are systematically reduced while monitoring the displacement response. The gap locations are determined by comparing the measured displacements with sets of calculated displacements and the sizes are estimated by monitoring the gap member length changes using the appropriate linear force-displacement relationship for the load level. The effect of measurement error in the truss displacements and the actuator length changes are investigated.

truss structures loose joints displacements actuat↗
NASA NTRS2008

Application of FUN3D and CFL3D to the Third Workshop on CFD Uncertainty Analysis

Two Reynolds-averaged Navier-Stokes computer codes - one unstructured and one structured - are applied to two workshop cases (for the 3rd Workshop on CFD Uncertainty Analysis, held at Instituto Superior Tecnico, Lisbon, in October 2008) for the purpose of uncertainty analysis. The Spalart-Allmaras turbulence model is employed. The first case uses the method of manufactured solution and is intended as a verification case. In other words, the CFD solution is expected to approach the exact solution as the grid is refined. The second case is a validation case (comparison against experiment), for which modeling errors inherent in the turbulence model and errors/uncertainty in the experiment may prevent close agreement. The results from the two computer codes are also compared. This exercise verifies that the codes are consistent both with the exact manufactured solution and with each other. In terms of order property, both codes behave as expected for the manufactured solution. For the backward facing step, CFD uncertainty on the finest grid is computed and is generally very low for both codes (whose results are nearly identical). Agreement with experiment is good at some locations for particular variables, but there are also many areas where the CFD and experimental uncertainties do not overlap.

Rumsey, C. L.↗
NASA NTRS2018

Image Guider Subsystem Analysis for the GHAPS Project

The Gondola for High-Altitude Planetary Science (GHAPS) project is a balloon-borne astronomical observatory designed operate in the UV, Visible, and near-mid IR spectral region. The GHAPS Optical Telescope Assembly (OTA) is designed around a one meter aperture narrow field-of-view telescope with near diffraction limited performance. GHAPS will utilize Wallops Arc-Second Pointing System (WASP) for pointing the OTA with an accuracy of 1 arc second or better. WASP relies heavily on a self-contained star tracker assembly to determine the OTA line of sight. Preliminary structural analysis indicated that potential misalignments could be present between the OTA line of sight and the star tracker FOV center during the expected flight conditions that could compromise GHAPS pointing accuracy. An analysis is presented that estimates the worst case long-term (or bias) pointing misalignments for the GHAPS opto-mechanical configuration, examines the need for additional hardware to correct pointing errors, and determines the best hardware and software implementation to do so. The analysis includes the development of a photometric model to predict the limiting magnitude of a GHAPS Facility Guide System (FGS) for a variety of image sensor technologies from visible to mid IR. FGS system performance is estimated using star magnitude and spectral data from actual star catalogs for more realistic predictions. Based on expected performance and other factors, such as size, weight, power, and cost, a final recommendation is made for the GHAPS configuration.

Telescope Guider↗
NASA NTRS1995

Task Analytic Models to Guide Analysis and Design: Use of the Operator Function Model to Represent Pilot-Autoflight System Mode Problems

Task-analytic models structure essential information about operator interaction with complex systems, in this case pilot interaction with the autoflight system. Such models serve two purposes: (1) they allow researchers and practitioners to understand pilots' actions; and (2) they provide a compact, computational representation needed to design 'intelligent' aids, e.g., displays, assistants, and training systems. This paper demonstrates the use of the operator function model to trace the process of mode engagements while a pilot is controlling an aircraft via the, autoflight system. The operator function model is a normative and nondeterministic model of how a well-trained, well-motivated operator manages multiple concurrent activities for effective real-time control. For each function, the model links the pilot's actions with the required information. Using the operator function model, this paper describes several mode engagement scenarios. These scenarios were observed and documented during a field study that focused on mode engagements and mode transitions during normal line operations. Data including time, ATC clearances, altitude, system states, and active modes and sub-modes, engagement of modes, were recorded during sixty-six flights. Using these data, seven prototypical mode engagement scenarios were extracted. One scenario details the decision of the crew to disengage a fully automatic mode in favor of a semi-automatic mode, and the consequences of this action. Another describes a mode error involving updating aircraft speed following the engagement of a speed submode. Other scenarios detail mode confusion at various phases of the flight. This analysis uses the operator function model to identify three aspects of mode engagement: (1) the progress of pilot-aircraft-autoflight system interaction; (2) control/display information required to perform mode management activities; and (3) the potential cause(s) of mode confusion. The goal of this paper is twofold: (1) to demonstrate the use of the operator functio model methodology to describe pilot-system interaction while engaging modes And monitoring the system, and (2) to initiate a discussion of how task-analytic models might inform design processes. While the operator function model is only one type of task-analytic representation, the hypothesis of this paper is that some type of task analytic structure is a prerequisite for the design of effective human-automation interaction.

Degani, Asaf↗
DOE OSTI2019

Viscosity of Iron Oxide Aluminosilicate Melts

A structure-based viscosity model has been developed for the CaO-MgO-FeO-SiO{sub 2}-Al{sub 2}O{sub 3} system and its subsystems. A critical analysis of the experimental data for iron oxide aluminosilicate systems has been conducted and used in the optimization of a model, which provides a link between the viscosity and the structures of the melts. The structural characteristics of aluminosilicate melt can be described by four types of oxygen: bridging oxygen O{sub i}{sup 0}, non-bridging oxygen O{sub j}{sup −}, free oxygen O{sup 2−}, and excess bridging oxygen O{sup *}. The present model is capable of predicting the viscosities in the CaO-MgO-FeO-SiO{sub 2}-Al{sub 2}O{sub 3} system and its subsystems over wide ranges of composition and temperature above the liquidus within experimental uncertainties, and the average of relative errors for this model was found to be 19.94 pct. This model includes FeO-SiO{sub 2}, CaO-SiO{sub 2}-FeO, Al{sub 2}O{sub 3}-SiO{sub 2}-FeO, MgO-SiO{sub 2}-FeO, CaO-SiO{sub 2}-MgO-FeO, CaO-SiO{sub 2}-Al{sub 2}O{sub 3}-FeO, and CaO-MgO-SiO{sub 2}-Al{sub 2}O{sub 3}-FeO. It was found that the effects of different metal cations on the viscosity can be determined, and the ability to reduce viscosity in the silicate melts follows the order: FeO > CaO > MgO. In the present model, the predicted viscosity decreased when CaO was initially initial added to the Al{sub 2}O{sub 3}-MgO-FeO-SiO{sub 2} system at a fixed SiO{sub 2} content because Ca{sup 2+} has the highest priority for charge compensation and FeO is the strongest modifier for this slag.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Data-Centric Development of Lignin Structure–Solubility Relationships in Deep Eutectic Solvents Using Molecular Simulations

Lignin is a natural source of aromatic chemicals with significant potential as an abundant, renewable feedstock for value-added products. Deep eutectic solvents (DES)–solvents composed of a hydrogen bond donor (HBD) and acceptor (HBA) in varying ratios–have emerged as a highly tunable class of solvents for lignin solubilization. However, the variety of possible DES compositions and limited molecular-scale understanding of lignin solubility makes solvent selection a challenge without laborious trial-and-error experimentation. To address these challenges, we use classical molecular dynamics (MD) simulations to study the interactions of lignin model compounds with various DES–water systems. Quantitative parameters (descriptors) were calculated by postprocessing the MD results and used to train a regression model that predicts experimentally determined solubilities of lignin model compounds. This approach revealed that the most important descriptors of solubility are the system temperature, solute hydrophilicity, and metrics quantifying hydrogen bonding. Maximizing the interactions between solute–HBD (hydrophobic group), water–HBD (hydrophilic group), and water–HBA molecules led to the highest model compound solubility. Our results support a hydrotropic mechanism in which extensive DES–water hydrogen bonding and favorable HBD interactions with the solute promote high solubility. We applied the regression model derived using model compounds to predict the solubility of representative lignin oligomers. The model predicted lignin oligomers’ solubilities in good agreement with experiments, indicating that the simulations of model compounds can be extended to predict the solubility of larger lignin compounds across a range of solvent compositions and temperatures. Furthermore, these findings provide new molecular-scale insight into lignin solubilization mechanisms and a new method for computationally screening potential solvent systems for lignin valorization.

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

Revisiting surface core-level shifts for ionic compounds

The shifts of core-level binding energies can provide powerful information about the electronic structure of a material. Understanding the physical origin of these shifts for catalytically relevant oxides may provide important insight into their properties. This requires reliable theoretical methods which are able to relate the binding energy shifts to the electronic structure. In order to establish such a methodology, the CaO(100) surface to bulk core-level binding energy shifts have been studied with Hartree-Fock and density-functional theory methods using both cluster and periodic slab models. The shifts obtained from the different theoretical methods are compared with each other and with data from synchrotron-based x-ray photoelectron spectroscopy (XPS) measurements. With a common approximation for the slab model treatment of XPS, the predicted binding energy shifts are seriously in error. In this work, the origin of the error is identified as arising from a flawed treatment of the surface atom binding energies, and a method for correcting the failure is presented.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

Tiling Framework for Heterogeneous Computing of Matrix based Tiled Algorithms

Tiling matrix operations can improve the load balancing and performance of applications on heterogeneous computing resources. Writing a tile-based algorithm for each operation with a traditional, hand-tuned tiling approach that uses for loops in C/C++ is cumbersome and error prone. Moreover, it must enable and support the heterogeneous memory management of data objects and also explore architecture-supported, native, tiled-data transfer APIs instead of copying the tiled data to continuous memory before the data transfer. The tiling framework provides a tiled data structure for heterogeneous memory mapping and parameterization to a heterogeneous task specification API. We have integrated our tiled framework into MatRIS (Math kernels library using IRIS). IRIS is a heterogeneous run-time framework with a heterogeneous programming model, memory model, and task execution model. Experiments reveal that the tiled framework for BLAS operations has improved the programmability of tiled BLAS and improved performance by ~20% when compared against the traditional method that copies the data to continuous memory locations for heterogeneous computing.

Miniskar, Narasinga Rao↗