Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Paste Formulation”

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 91 records · Page 5
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Updates to Nek5000: RANS wall functions, NRC support, and documentation

As a key part of the U.S. Department of Energy’s Nuclear Energy Advanced Modeling and Simulation program (NEAMS), the next generation of tools for nuclear reactor design are being developed. A primary metric for success of this program is the adoption of those tools by industry stakeholders. To help meet this metric for Nek5000, development efforts this year have focused on the expansion of the available Reynolds-averaged Navier-Stokes (RANS) turbulence models to include wall function models, continued support of the U.S. Nuclear Regulatory Commission in simulating a hydrogen mitigation benchmark experiment, and expansion of the Nek5000 documentation to enhance the code’s usability. The implementation of wall functions is intended to significantly decrease time to-solution for a wide range of problems and it has been identified as an important feature by stakeholders across the nuclear community. The initial implementation of standard and pressure corrected wall functions has been tested in Nek5000 on the classical problems of channel flow and a backward facing step. It has also been tested on a 2D molten salt fast reactor core and a T-junction. The recommended formulation of choosing y + = 30 as a boundary value combined with a Neumann-Neumann formulation of the boundary condition for k and τ produces the most consistent results compared to a wall resolved approach. The collaboration with the NRC is a continuing exercise that has been ongoing for multiple years. This year we have concluded investigations into appropriate inlet conditions, showing that fully developed turbulent conditions are adequate for the full PANDA domain. We have also concluded simulations for an unobstructed jet and are now refocusing to the original obstructed jet case. The expansion of the code documentation has also been an ongoing effort over the past few years to address feedback from the various stakeholders. The documentation now more fully covers many of the available features and closes identified SQA gaps. Finally we report on three training sessions that were offered over the past year in an effort to expand the Nek5000 user base.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2020

Efficient and robust stability analysis in the internal energy, volume, and moles (UVN) space

Hydrodynamic simulation of fluids in vessels is often based on the internal energy (U), volume (V), and the number of moles of all components ( re N = N 1 , N 2 , ..., N n ) as the working variables to determine the thermodynamic state by finding the global maximum of entropy. The procedure of finding the equilibrium state may be divided into two steps: stability analysis and phase-split computations. Stability analysis is performed first, and the phase-split computation proceeds if the system is unstable. Most of the past work on UVN space has been based on phase-split computations, while the stability analysis is barely investigated. Furthermore, previous studies for the stability analysis use only Newton's method to solve the non-linear algebraic equations. We present an efficient and robust approach for stability analysis where UVN is specified. The successive substitution iteration (SSI) is used to provide good initial guesses for Newton's method. The proposed approach results in a reduced number of unknowns and does not require a large number of iterations to achieve convergence in Newton's method. Overall, our proposed formulation is compatible with different equations of state and is applicable to both pure component and multi-components. The robustness and efficiency of the algorithm for stability analysis are demonstrated in various examples.

02 PETROLEUM↗
DOE OSTI2022

Geometry, Disorder and Phase Transitions in Topological States of Matter

The quantum Hall effect is the birthplace of topological states of matter, a major theme at the forefront of condensed matter physics in the past two decades. The fractional quantum Hall (FQH) effect revolutionized our understanding of phases of electronic matter. FQH states support exotic fractionally charged excitations that obey Abelian or non-Abelian fractional statistics, which are topological excitations that result from the underlying topological order. During this project, our group discovered a previously unrecognized geometric degree of freedom of incompressible FQH states and studied that for a variety of gapped FQH states. We brought this new concept into direct contact with experiments for the first time by generalizing it to Fermi-liquid states of composite fermions. Using the newly formulated powerful infinite Density Matrix Renormalization Group method, our numerical calculations yielded a parameter free prediction that was found to be in excellent agreement with experimental findings on electron systems in semiconductor heterostructures. In parallel, we performed extensive numerical studies on different, competing phases at various Landau level filling factors, and quantum phase transitions that result from such a competition, e.g. Abelian-non-Abelian phase transitions in bilayer systems. We studied geometrical excitations dubbed “gravitons” (because of their analogy with excitations in the theory of gravitation) and ways to excite and detect them, and explored how they couple with topological excitations. In graphene-based chiral materials, we realized the ability to tune through different incompressible and compressible states in a single Landau level, and found appropriate experimental parameters for the exploration of universal Luttinger liquid behavior not obtained in semiconductor-based electron systems. We showed that topological systems had a very different response from nontopological systems to strong disorder (many-body localization) as well as periodic drives.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

Theory Manual for the Fuel Cycle Analysis Code REBUS

REBUS is a system of programs designed for the fuel cycle analysis of fast reactors and accelerator-driven systems. The first version of REBUS was developed in 1970 for calculating the equilibrium operating conditions of liquid metal fast breeder reactors. During the past five decades, various extensions and refinements have been made in the simulation of material motions in fuel cycle facilities and the solution of the neutronics and isotopic transmutation equations that govern the reactor behavior. The latest version complies with the standard code practices and interface dataset specifications of the committee on computer code coordination (CCCC). This report provides a detailed description of the physical processes and their mathematical models and solution methods of the current production version of REBUS. The models of the facilities and physical processes involved in fast reactor fuel cycle are first discussed. Then the mathematical formulation of these physical models is presented, followed by detailed discussions of solution methods of the resulting set of equations. Finally the output edit quantities are discussed along with their definitions.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2021

Printing thermoelectric inks toward next-generation energy and thermal devices

The ability of thermoelectric (TE) materials to convert thermal energy to electricity and vice versa highlights them as a promising candidate for sustainable energy applications. Despite considerable increases in the figure of merit zT of thermoelectric materials in the past two decades, there is still a prominent need to develop scalable synthesis and flexible manufacturing processes to convert high-efficiency materials into high-performance devices. Scalable printing techniques provide a versatile solution to not only fabricate both inorganic and organic TE materials with fine control over the compositions and microstructures, but also manufacture thermoelectric devices with optimized geometric and structural designs that lead to improved efficiency and system-level performances. In this review, we aim to provide a comprehensive framework of printing thermoelectric materials and devices by including recent breakthroughs and relevant discussions on TE materials chemistry, ink formulation, flexible or conformable device design, and processing strategies, with an emphasis on additive manufacturing techniques. Additionally, we review recent innovations in the flexible, conformal, and stretchable device architectures and highlight state-of-the-art applications of these TE devices in energy harvesting and thermal management. Perspectives of emerging research opportunities and future directions are also discussed. While this review centers on thermoelectrics, the fundamental ink chemistry and printing processes possess the potential for applications to a broad range of energy, thermal and electronic devices.

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

Influences of δB contribution and parallel inertial term of energetic particles on MHD-kinetic hybrid simulations: a case study of the 1/1 internal kink mode

The magnetohydrodynamic-kinetic (MHD-kinetic) hybrid model (Park et al 1992 Phys. Fluids B 4 2033–7) has been widely applied in studying energetic particles (EPs) problems in fusion plasmas for past decades. The pressure-coupling scheme or the current-coupling scheme is adopted in this model. However, two noteworthy issues arise in the model application: firstly, the coupled term introduced in the pressure-coupling scheme, (∇•P h ) ⟂ , is often simplified by ∇•P h , which is equivalent to neglecting the parallel inertial term of EPs; secondly, besides the $δf$ contribution caused by changing in the EP distribution function, the magnetic field perturbation (the $δB$ contribution) generated during development of the instabilities should also be considered, but it is often ignored in existing hybrid simulations. In this paper, we derive the analytical formulations under these two coupling schemes and then numerically study the representative case of the linear stability of the $m/n$ = $1/1$ internal kink mode (IKM) (Fu et al 2006 Phys. Plasmas 13 052517) by using the CLT-K code. Further, it is found that the approximated models can still yield reasonable results when EPs are isotopically distributed. But it fails completely in cases with anisotropic EP distributions. In addition, we further investigate the influence of EP's orbit width on the stability of IKM and verify the equivalence between pressure-coupling scheme and the current-coupling scheme.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Exploring the Neutron Substructure with Advanced Polarized Helium-3 Targets (Or: How I Learned to Stop Worrying and Love Spectroscopy)

As we seek to understand the smallest, physical aspects of our universe, we cannot simply rely on our senses to probe the world around us as we did in the past. The smallest physical elements of our universe behave in strange, probabilistic ways and are completely invisible to the naked eye/ear/etc. So, we design clever experiments (such as scattering experiments) to probe these minute realms. Then, just as with the larger, observable world, we devise models and equations to describe what we think is happening. Due to the nature of the physical universe at the quantum scale and with the aid of symmetries such as Lorentz invariance, we can write down equations that describe the scattering, but the expressions contain functions, which we call ?form factors? and ?structure functions?, that we cannot compute from first principles. We can, however, formulate models that make predictions for these functions. By comparing our predictions with the observed data, we can gain insight into the validity of our models and thus a better physical understanding of what is happening at these minuscule scales. Studying the constituents inside of the nucleus of an atom adds another layer of difficulty if we can?t remove those components from the nucleus. This is the case with the neutron. When not bound in the nucleus with protons and other neutrons, the neutron will decay into a proton after about 15 minutes. So, we?re forced to study the neutron while it is still bound in the nucleus of an atom such as helium-3 (3He). For the last 1,000 years (rounding up), our group has developed high quality, polarized 3He targets made of an aluminosilicate glass. These targets are made in order to perform experiments at Jefferson Lab (JLab), experiments which let us determine the form factors and structure functions of the neutron by scattering polarized electrons from polarized neutrons (or rather polarized 3He). The specific experiments reported on in this thesis push the bounds of our understanding of the internal structure of the neutron. Good science is often about pushing experimental techniques to a new level. Toward that goal we study our polarized 3He targets both to advance the technology and to choose the best ones for our experiments. We do this using a process called nuclear magnetic resonance (NMR) to gauge the maximum polarization of a target and how fast the polarization decays with time. While these tests primarily provide us information that make analysis of our experimental scattering data possible, they also let us determine whether or not a target-cell is useful or even, dare I say, of spectacular quality. Our latest targets utilize a novel convection design allowing 3He to be polarized and quickly moved in front of the electron-beam, making it possible to use larger targets with higher electron-beam currents than ever before. This means more electrons scatter and we get more data. And by studying our targets in detail prior to using them in our experiments, we have found techniques to take effects which could have been detrimental to target quality and turn them to our advantage! It?s a real case of making lemonade out of lemons. We also use laser spectroscopy to study the absorption lines of alkali-metals in the target (potassium and rubidium, specifically). We add these alkali-metals to our target to facilitate polarizing the 3He. We can use the measurement of these pressure broadened absorption lines to determine the 3He density inside of the target with great precision. Historically, we understood the width of these lines would be dependent on the temperature of the target. Specifically, if I raise the temperature, the width should get bigger. I found that was not the case, which was very confusing at first, though very exciting now that I realize the data are self-consistent and suggestive of unexpected behavior. This thesis details the development of high quality, glass, polarized 3He targets for the 2020 An 1 /dn 2 and 2023 Gn E experiments, which utilized the first 3He convection targets and broke records in target quality. This thesis also covers the initial development of metal windows for the next-generation of 3He target-cells. Finally, this thesis documents the temperature dependence of the width of potassium (K) and rubidium (Rb) absorption lines as measured with laser spectroscopy.

Jantzi, Christopher↗
DOE OSTI2025

Establishing nationwide power system vulnerability index across US counties using interpretable machine learning

Power outages have become increasingly frequent, intense, and prolonged in the US due to climate change, aging electrical grids, and rising energy demand. However, largely due to the absence of granular spatiotemporal outage data, we lack data-driven evidence and analytics-based metrics to quantify power system vulnerability. This limitation has hindered the ability to effectively evaluate and address vulnerability to power outages in US communities. Here, in this work, we collected ∼179 million power outage records at 15-min intervals across 3022 US contiguous counties (96.15 % of the area) from 2014 to 2023. We developed a power system vulnerability assessment framework based on three dimensions (intensity, frequency, and duration) and applied interpretable machine learning models (XGBoost and SHAP) to compute Power System Vulnerability Index (PSVI) at the county level. Our analysis reveals a consistent increase in power system vulnerability across the US counties over the past decade. We identified 318 counties across 45 states as hotspots for high power system vulnerability, particularly in the West Coast (California and Washington), the East Coast (Florida and the Northeast area), the Great Lakes megalopolis (Chicago-Detroit metropolitan areas), and the Gulf of Mexico (Texas). Our heterogeneity analysis indicates that urban counties and those located along regional transmission boundaries tend to exhibit significantly higher vulnerability. Our results highlight the significance of the proposed PSVI for evaluating the vulnerability of communities to power outages. The findings underscore the widespread and pervasive impact of power outages across the country and offer crucial insights to support infrastructure operators, policymakers, and emergency managers in formulating policies and programs aimed at enhancing the resilience of the US power infrastructure.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2026

Avoidance of disruptions on KSTAR due to vertical displacement events via novel real-time stability assessment

Disruption avoidance via the DECAF approach has been achieved on KSTAR using a novel real-time vertical stability assessment and a multiactuator feedback control strategy. The development of disruption avoidance strategies with reactor-relevant reliability is an urgent activity, enabling future fusion power plants. The stability metric employed is based on a new formulation of a vertical force gradient balance metric evaluated across the poloidal cross section of the plasma, with parameters tuned using historical data. Evaluation of this metric on a validation set of 400 recent KSTAR shots indicates >82% of Vertical displacement events can be avoided via feedback control. Essential to its calculation is the two-dimensional toroidal current density distribution in the plasma. Measurement of this profile faster than fully-converged equilibrium reconstructions can deliver is found to improve forecaster performance and is achieved with a surrogate model that takes as input magnetic diagnostic measurements and outputs the current profile on a basis comprising the top principal components of historical current profiles (from past equilibrium reconstructions). This method solves the non-uniqueness problem typically faced when reconstructing current profiles directly from diagnostics, while improving computational time and accuracy. On average, profiles produced by this model reach coefficients of determination of >0.99 with respect to those from equilibrium reconstructions. The avoidance actuators employed include poloidal field coils and an electron cyclotron current drive system. The multiactuator approach, as shown in this first demonstration, allows disruption avoidance while minimizing impact to operational performance. This ability, along with its flexibility and speed, makes this new approach an attractive option for avoiding these types of disruptions in reactors.

Tobin, Matthew [Columbia Univ., New York, NY (Unit↗
DOE OSTI2023

NIMROD Development AND Applications for Advanced Simulations of Tokamak Plasmas

Magnetic fusion experiments have the largest gradients in the world that are in steady state. These gradients lead to instabilities, and thus, most of plasma theory over the past four decades have been devoted to the development of instability theory and their use in interpreting plasma results. Instabilities in tokamaks lead to 3 different phenomenological outcomes: coherent saturation, turbulent saturation, and sudden relaxation (which may include disruptions). For long-wavelength instabilities, the NIMROD extended MHD code has emerged as an important tool for understanding tokamak instabilities. Because of the long history of NIMROD, a new version of it has been started at Tech-X to be able to address multi-species capabilities as well as exploit modern GPU systems. The new version builds on the previous version and improves the workflow by enabling new equilibria from experiments to be resolved. The fundamental formulation of the equations underlying NIMROD has also been explored in multiples ways. First, understanding instabilities in tokamaks through numerical methods is aided by being able to understand the rich history of analytic studies. A paper was published which aid in understanding the literature by simplifying some of the analytic machinery inherent in these studies. Extended MHD, like gyrokinetics, are quasineutral models. From a theoretical point of view, plasma quasineutrality approximation is best understood as the same as the magneto-quasistatic approximation of the Maxwell equations. Finally, the standard model of tokamak theory is that of instability theory. A simple dynamical systems model has been developed to better illustrate the strengths and weaknesses of this model.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

A survey of nonlinear robust optimization

Robust optimization (RO) has attracted great attention from the optimization community over the past decade. RO is dedicated to solving optimization problems subject to uncertainty: design constraints must be satisfied for all the values of the uncertain parameters within a given uncertainty set. Uncertainty sets may be modeled as deterministic sets (boxes, polyhedra, ellipsoids), in which case the RO problem may be reformulated via worst-case analysis, or as families of distributions. The challenge of RO is to reformulate or approximate robust constraints so that the uncertain optimization problem is transformed into a tractable deterministic optimization problem. Most reformulation methods assume linearity of the robust constraints or uncertainty sets of favorable shape, which represents only a fraction of real-world applications. This survey addresses nonlinear RO and includes problem formulations and applications, solution approaches, and available software with code samples.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Site demonstration and performance evaluation of MPC for a large chiller plant with TES for renewable energy integration and grid decarbonization

Thermal energy storage (TES) for a cooling plant is a crucial resource for load flexibility. Traditionally, simple, heuristic control approaches, such as the storage priority control which charges TES during the nighttime and discharges during the daytime, have been widely used in practice, and shown reasonable performance in the past benefiting both the grid and the end-users such as buildings and district energy systems. However, the increasing penetration of renewables changes the situation, exposing the grid to a growing duck curve, which encourages the consumption of more energy in the daytime, and volatile renewable generation which requires dynamic planning. The growing pressure of diminishing greenhouse gas emissions also increases the complexity of cooling TES plant operations as different control strategies may apply to optimize operations for energy cost or carbon emissions. This paper presents a model predictive control (MPC), site demonstration and evaluation results of optimal operation of a chiller plant, TES and behind-meter photovoltaics for a campus-level district cooling system. The MPC was formulated as a mixed-integer linear program for better numerical and control properties. Compared with baseline rule-based controls, the MPC results show reductions of the excess PV power by around 25%, of the greenhouse gas emission by 10%, and of peak electricity demand by 10%.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗
DOE OSTI2022

Printable Fiber Reinforced Cement Composites – Feasibility Study

Additive manufacturing is enabling the manufacturability of structures with previously unattainable complexity or functionality, and there is growing interest in additive manufacturing of “printed” concrete structures. The focus of this Phase 1 Technical Collaboration (TC) project was to evaluate feasibility of printing hybrid cement composite structures reinforced with textile carbon fibers (tCF). This project leverages other (non-IACMI) projects on cement formulation and printing process development, as well as on the production process for tCF. This project’s primary focus was to explore cement composite mix design with textile carbon fibers to be manufactured by MonteFibre (TC partner) and evaluate suitable fiber-matrix interface or sizing for cement composites working with Michelman (TC Partner). This project supports IACMI’s goal of reducing the cost and embodied energy of carbon fiber composites. Cost is one of the fundamental challenges to carbon fiber reinforced cement composites. Cement is an extremely inexpensive material (approximately $\$$0.05/lb). Adding 1 wt% of conventional carbon fiber to cement quadruples its cost. Therefore, the need to use low-cost carbon fiber and ensure that the additional cost of the carbon fiber has a greater cost benefit to the final product. This was the first preliminary evaluation to integrate tCF reinforcement in cement composites, and such potential tCF utilization should significantly reduce materials cost. Cement composite production is energy and emissions intensive, thus by strengthening it less material will be required. Hence, the embodied energy and production time of the resulting structures will be reduced. Additionally, integrating these new materials into additive processes can enable selective use of the material in high load or stress areas. It is noteworthy that past work in this field of fiber reinforced cement composites did not consider the optimization of fiber-matrix interface using suitable sizing. Carbon fiber reinforcement offers potential added benefits of thermal conductivity (which affects cure rate) and flow behavior that could provide opportunities for site specific utilization of carbon fiber on hybrid cement structures (e.g. use the fiber reinforcement on outer surfaces to enhance strength and modulus and then infiltrating the internal structures with conventional concrete). MonteFibre was the industry lead and planned on supplying the tCF for this project. However, during the short Phase-1 duration of this project, MonteFibre was unable to produce tCF for this project due to manufacturing plant being off-line throughout the course of the project. The project team decided to pursue an alternate option which involved demonstrating printable concrete with steel fibers by the ORNL lead, Dr. Brian Post. The University of Tennessee collaboration team focused on evaluating the suitable chemical sizing for carbon fibers working with Michelman and also developed methods for material characterization of cement-based composites to evaluate the material response for compression, shear, flexure, and tension. The two milestones for University of Tennessee, Knoxville were realized related to identification of one sizing suitable for carbon fiber reinforced cement composite and developing data associated with mechanical behavior of unreinforced (neat) and carbon fiber reinforced cement composites. ORNL could not complete the task of carbon fiber reinforced printed cement composites due to the reasons mentioned earlier, but was able to replace tCF with steel fibers to demonstrate the feasibility of printing with fiber reinforced cement composites. The Project team reviewed possible sizing chemistry available in collaboration with Michelman for use on carbon fiber reinforcement in cement composites and concrete applications. Our initial goal was to identify a sizing most promising for formulation with textile carbon fibers (tCF) to deliver excellent mechanical properties in composite material state. Since tCF was not available for this project as originally envisioned, the team continued this task to identify a suitable sizing for carbon fiber applications by applying such sizing to lower cost carbon fibers currently available commercially from Zoltek called Panex fibers. At a future time this can be optimized for textile carbon fibers from Montefibre. The bulk of previous work on carbon fiber reinforced cement has neglected the importance of fiber-matrix adhesion on mechanical properties of the cement composite and identifying this missing link was an important accomplishment for future research. Tensile behavior of fiber reinforced concrete is important to evaluate in order to realize the dream of concrete products that do not need reinforcing steel. Important sample preparation and testing procedures were addressed in this study and it was concluded that substantial improvements in tensile behavior, without compromising compressive strength, and improved ductility can result from the use of carbon fiber reinforcement.

36 MATERIALS SCIENCE↗
DOE OSTI2026

ACT University Collaboration Proposal: Role of manufacturing defects on material failure under dynamic loading for developing enhanced failure models and theories (Final Report)

Research on the tensile behavior of additively manufactured 316L stainless steel coupon specimens at increasing strain rates was conducted over the past year at Penn State. Dynamic loading rates into the 1000 1/s loading rates were performed on specimens with nearly 1.32” gage lengths. Stress-strain plots show ductile and plastic behavior well beyond 12% strain with necked specimens, having smaller effective gage lengths showing up to 56% ultimate strain. Additional tests performed on compact tension specimens helped with simulation work to understand the deformation behavior. Using finite elements, it was possible to determine the feasibility of a comprehensive experimental validation study towards an improved model for AM material failure – specifically, the Bai-Wierzbicki approach which accounts for lode angle and triaxiality. A non-significant number of tests are projected for eight specific failure nodes, with additional replicates to provide strain-rate capabilities to the existing formulation. One such correction factor is explored for quasi-static, notched specimens. Finally, stress intensity factor was explored using a set of compact-tension specimens which would also provide useful validation data for any simulation work.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Mechanochemical formation of highly stable amorphous calcium carbonate in calcium silicates: Potential long-term storage of CO 2 in cement

Here, we have investigated mechanochemical reactions with the calcium silicate wollastonite to probe potential mechanisms of sequestration and long-term storage of CO 2 as mineral carbonates in cement pastes. Wollastonite, CaSiO 3 , was milled under ambient and 13 C enriched CO 2 atmospheres. Milling induced a structure change from monoclinic wollastonite-2M to triclinic-1T, consistent with high pressure treatment. The results from solid-state 13 C nuclear magnetic resonance (NMR) spectroscopy shows this process forms an extraordinarily stable amorphous calcium carbonate (ACC) phase that remains despite a high temperature dehydration treatment (130 °C) and persists in the same form after 3 yrs. of storage at ambient conditions. Thermogravimetric analysis (TGA) of these samples indicates that ACC is produced at the similar solid concentrations in both CO 2 -enriched and ambient atmospheres. The observation of meta-stable ACC due to wollastonite carbonation may play an important role in further exploration and optimization of mineral carbonation reactions for CO 2 capturing cement formulations.

13C solid-state NMR↗
DOE OSTI2021

Data assimilation empowered neural network parametrizations for subgrid processes in geophysical flows

In the past couple of years, there has been a proliferation in the use of machine learning approaches to represent subgrid-scale processes in geophysical flows with an aim to improve the forecasting capability and to accelerate numerical simulations of these flows. Despite its success for different types of flow, the online deployment of a data-driven closure model can cause instabilities and biases in modeling the overall effect of subgrid-scale processes, which in turn leads to inaccurate prediction. To tackle this issue, we exploit the data assimilation technique to correct the physics-based model coupled with the neural network as a surrogate for unresolved flow dynamics in multiscale systems. In particular, we use a set of neural network architectures to learn the correlation between resolved flow variables and the parametrizations of unresolved flow dynamics and formulate a data assimilation approach to correct the hybrid model during their online deployment. We illustrate our framework in a set of applications of the multiscale Lorenz 96 system for which the parametrization model for unresolved scales is exactly known, and the two-dimensional Kraichnan turbulence system for which the parametrization model for unresolved scales is not known a priori. Our analysis, therefore, comprises a predictive dynamical core empowered by (i) a data-driven closure model for subgrid-scale processes, (ii) a data assimilation approach for forecast error correction, and (iii) both data-driven closure and data assimilation procedures. We show significant improvement in the long-term prediction of the underlying chaotic dynamics with our framework compared to using only neural network parametrizations for future prediction. Moreover, we demonstrate that these data-driven parametrization models can handle the non-Gaussian statistics of subgrid-scale processes, and effectively improve the accuracy of outer data assimilation workflow loops in a modular nonintrusive way.

42 ENGINEERING↗
DOE OSTI2021

Scalable Pattern Matching in Metadata Graphs via Constraint Checking

Pattern matching is a fundamental tool for answering complex graph queries. Unfortunately, existing solutions have limited capabilities: They do not scale to process large graphs and/or support only a restricted set of search templates or usage scenarios. Moreover, the algorithms at the core of the existing techniques are not suitable for today’s graph processing infrastructures relying on horizontal scalability and shared-nothing clusters, as most of these algorithms are inherently sequential and difficult to parallelize. In this article we present an algorithmic pipeline that bases pattern matching on constraint checking. The key intuition is that each vertex and edge participating in a match has to meet a set of constraints implicitly specified by the search template. These constraints can be verified independently and typically are less expensive to compute than searching the full template. The pipeline we propose generates these constraints and iterates over them to eliminate all the vertices and edges that do not participate in any match, thus reducing the background graph to a subgraph that is the union of all template matches—the complete set of all vertices and edges that participate in at least one match. Additional analysis can be performed on this annotated, reduced graph, such as full match enumeration, match counting, or computing vertex/edge centrality. Furthermore, a vertex-centric formulation for constraint checking algorithms exists, and this makes it possible to harness existing high-performance, vertex-centric graph processing frameworks. This technique (i) enables highly scalable pattern matching in metadata (labeled) graphs; (ii) supports arbitrary patterns with 100% precision; (iii) enables tradeoffs between precision and time-to-solution, while always selects all vertices and edges that participate in matches, thus offering 100% recall; and (iv) supports a set of popular data analytics scenarios. We implement our approach on top of HavoqGT, an open-source asynchronous graph processing framework, and demonstrate its advantages through strong and weak scaling experiments on massive scale real-world (up to 257 billion edges) and synthetic (up to 4.4 trillion edges) labeled graphs, respectively, and at scales (1,024 nodes / 36,864 cores), orders of magnitude larger than used in the past for similar problems. This article serves two purposes: First, it synthesises the knowledge accumulated during a long-term project. Second, it presents new system features, usage scenarios, optimizations, and comparisons with related work that strengthen the confidence that pattern matching based on iterative pruning via constraint checking is an effective and scalable approach in practice. The new contributions include the following: (i) We demonstrate the ability of the constraint checking approach to efficiently support two additional search scenarios that often emerge in practice, interactive incremental search and exploratory search. (ii) We empirically compare our solution with two additional state-of-the-art systems, Arabsque and TriAD. (iii) We show the ability of our solution to accommodate a more diverse range of datasets with varying properties, e.g., scale, skewness, label distribution, and match frequency. (iv) We introduce or extend a number of system features (e.g., work aggregation, load balancing, and the ability to cap the generated traffic) and design optimizations and demonstrate their advantages with respect to improving performance and scalability. (v) We present bottleneck analysis and insights into artifacts that influence performance. (vi) We present a theoretical complexity argument that motivates the performance gains we observe.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

How does humidity data impact land surface modeling of hydrothermal regimes at a permafrost site in Utqia$\mathrm{\dot{g}}$vik, Alaska?

Humidity is a basic and crucial meteorological indicator commonly measured in several forms, including specific humidity, relative humidity, and absolute humidity. These different forms can be inter-derived based on the saturation vapor pressure (SVP). In past decades, dozens of formulae have been developed to calculate the SVP with respect to, and in equilibrium with, liquid water and solid ice surfaces, but many prior studies use a single function for all temperature ranges, without considering the distinction between over the liquid water and ice surfaces. These different approaches can result in humidity estimates that may impact our understanding of surface-subsurface thermal-hydrological dynamics in cold regions. In this study, we compared the relative humidity (RH) downloaded and calculated from four data sources in Alaska based on five commonly used SVP formulas. These RHs, along with other meteorological indicators, were then used to drive physics-rich land surface models at a permafrost-affected site. Here, we found that higher values of RH (up to 40%) were obtained if the SVP was calculated with the over-ice formulation when air temperatures were below freezing, which could lead to a 30% maximum difference in snow depths. The choice of whether to separately calculate the SVP over an ice surface in winter also produced a significant range (up to 0.2 m) in simulated annual maximum thaw depths. The sensitivity of seasonal thaw depth to the formulation of SVP increases with the rainfall rate and the height of above-ground ponded water, while it diminishes with warmer air temperatures. These results show that RH variations based on the calculation of SVP with or without over-ice calculation meaningfully impact physically-based predictions of snow depth, sublimation, soil temperature, and active layer thickness. Under particular conditions, when severe flooding (inundation) and cool air temperatures are present, care should be taken to evaluate how humidity data is estimated for land surface and earth system modeling.

54 ENVIRONMENTAL SCIENCES↗