Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “MODELING”

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,351 records · Page 75

Surrogate model for third-integer resonance extraction at the Fermilab Delivery Ring

We present an ongoing work in which a surrogate model is being developed to reproduce the response dynamics of the third-integer resonant extraction process in the Delivery Ring (DR) at Fermilab. This is in pursuit of smoothly extracting circulating beam to the Mu2e Experiment’s production target, whereby the goal is to extract a uniform slice of the circulating 1e12 protons in the DR over 25,000 turns (43 ms). The DR contains 3 harmonic sextupoles that excite a third-integer resonance and three fast, tune-ramping quadrupole magnets that drive the horizontal tune towards the 29/3 resonance. In our initial work, the surrogate model trains on a semi-analytical simulation provided in the same format as live data. Using Reinforcement Learning (and other potential ML methods), the trained surrogate acts as the “environment” in which a simple ML control agent could learn to dynamically adjust the quadrupole ramp at 430 break points within the 43 microsecond spill window. The controller will be hosted on a dedicated Arria 10 FPGA. In this work, we report the accuracy and fidelity of the surrogate model in comparison to the response dynamics of the physics simulator.

Narayanan, Aakaash [Fermilab]↗

TriBITS Lifecycle Model: A Lean/Agile Software Lifecycle Model for Research-based Computational Science and Engineering and Applied Mathematical Software (V.1.0)

Software lifecycles are becoming an increasingly important issue for computational science and engineering (CSE) software. The process by which a piece of CSE software begins life as a set of research requirements and then matures into a trusted high-quality capability is both commonplace and extremely challenging. Although an implicit lifecycle is obviously being used in any effort, the challenges of this process - respecting the competing needs of research vs. production - cannot be overstated. Here we describe a proposal for a well-defined software lifecycle process based on modern Lean/Agile software engineering principles. What we propose is appropriate for many CSE software projects that are initially heavily focused on research but also are expected to eventually produce usable high-quality capabilities. The model is related to TriBITS, a build, integration and testing system, which serves as a strong foundation for this lifecycle model, and aspects of this lifecycle model are ingrained in the TriBITS system. Here, we advocate three to four phases or maturity levels that address the appropriate handling of many issues associated with the transition from research to production software. The goals of this lifecycle model are to better communicate maturity levels with customers and to help to identify and promote Software Engineering (SE) practices that will help to improve productivity and produce better software. An important collection of software in this domain is Trilinos, which is used as the motivation and the initial target for this lifecycle model. However, many other related and similar CSE (and non-CSE) software projects can also make good use of this lifecycle model, especially those that use the TriBITS system. Indeed this lifecycle process, if followed, will enable large-scale sustainable integration of many complex CSE software efforts across several institutions.

99 GENERAL AND MISCELLANEOUS↗

Reduced-order model based feedback control of the Modified Hasegawa-Wakatani model

In this work, the development of model-based feedback control that stabilizes an unstable equilibrium is obtained for the Modified Hasegawa-Wakatani (MHW) equations, a classic model in plasma turbulence. First, a balanced truncation (a model reduction technique that has proven successful in flow control design problems) is applied to obtain a low dimensional model of the linearized MHW equation. Then a model based feedback controller is designed for the reduced order model using linear quadratic regulators (LQR). Finally, a linear quadratic gaussian (LQG) controller, which is more resistant to disturbances is deduced. The controller is applied on the non-reduced, nonlinear MHW equations to stabilize the equilibrium and suppress the transition to drift-wave induced turbulence.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Final Report: Process Models of the Equilibrium Size & State of Organic/Inorganic Aerosols for the Development of Large Scale Atmospheric Models & the Analysis of Field Data

Our work addressed the following elements of the Call for Proposals: (i) “to improve the theoretical representation of aerosol processes studied in ASP laboratory or field studies”, (ii) “to enhance the incorporation of aerosol process information into modules suitable for large-scale or global atmospheric models”, and (iii) “provide systematic experimental validation of process model predictions ... using data from targeted laboratory and field experiments”. Achievements to the end of 2012 are described in four previous reports, and include: new models of densities and surface tensions of pure (single solute) and mixed aqueous solutions of typical aerosol composition under all atmospheric conditions (0 to 100% RH and T > 150 K); inclusion of these models into the widely used Extended Aerosol Inorganics model (E-AIM, http://www.aim.env.uea.ac.uk/aim/aim.php); the addition of vapor pressure calculators for organic compounds to the E-AIM website; the ability of include user-defined organic compounds and/or lumped surrogates in gas/aerosol partitioning calculations; the development of new equations to represent the properties of soluble aerosols over the entire concentration range (using methods based upon adsorption isotherms, and derived using statistical mechanics), including systems at close to zero RH. These results are described in publications 1-6 at the end of this report, and on the “News” page of the E-AIM website (http://www.aim.env.uea.ac.uk/aim/info/news.html). During 2012 and 2013 we have collaborated in a combined observation and lab-based study of the water uptake of the organic component of atmospheric aerosols (PI Gannet Hallar, of the Desert Research Institute). The aerosol samples were analyzed using several complementary techniques (GC/MS, FT-ICR MS, and ion chromatography) to produce a very complete organic “speciation” including both polar and non-polar compounds. Hygroscopic growth factors of the samples were measured, and we have just completed comparisons of the data with our process model predictions based upon the inorganic and organic composition of the samples.

54 ENVIRONMENTAL SCIENCES↗

Heat Pump Water Heater Using Wrapped-Tank Microchannel Condenser Coil for Charge Reduction – FY19 1st Quarter (Regular) Milestone Report: Development of Wrapped-Tank Condenser Model (Round Tube and Microchannel), Coupled with Water Tank Model

We recently grew the DOE/HPDM model library to enhance its heat pump water heater (HPWH) design and simulation capability. The new additions include wrapped-tank round tube coil, microchannel coil and a one-dimensional stratified water tank. A wrapped-tank condenser coil has strong coupling with a stratified water tank, which leads HPWH simulation a transient process. The new one-dimensional stratified water tank model is an improvement to the open-source EnergyPlus water tank model, by introducing a calibration factor to account for bulk mixing effect due to water draws and circulations, etc. The wrapped-tank condenser coil model, using a segment-to-segment modeling approach to simulate round tubes or microchannel ports. The heat exchanger models can accept arbitrary refrigerant types, e.g. HFCs, HFOs, natural refrigerants like CO 2 and propane, etc.

42 ENGINEERING↗

Model Package Report: Geoframework Model of the Hanford Site 100 Area

The purpose of the 100 Area Geoframework Model (GFM) is to provide a reasonable, consistent, and defensible three dimensional representation of the hydrostratigraphic units below the River Corridor at the Hanford Site to support contaminant fate and transport models. The GFM is a three dimensional representation of the subsurface geologic structure. From this, three dimensional geologic model-exported results, in the form of points or surfaces, are used as inputs to populate and assemble the various numerical model architectures. The objective of this report is to define the process used to produce a hydrostratigraphic model for the hydrostratigraphic units beneath the Hanford Site 100 Area. The GFM may support several other CH2M HILL Plateau Remediation Company project needs and objectives, including providing geologic information to support remedial investigations and actions and as a tool to present River Corridor geology and contaminant extents. The GFM was constructed based on information through 2019 available in the Integrated Data Management System, Hanford Environmental Information System, and the Hanford Site geologic contacts (GeoContacts) datasets. Revisions to the 100 Area GFM may occur to incorporate new data and information. Each version will be maintained in configuration control using date stamps for identifying supporting databases, figures, and interpretations. Supporting data include the final three dimensional geoframework surfaces, two dimensional structure and isopach maps, and all the geologic contact inputs and interpreted geologic log data. These final products provide a supporting set of geologic information that together allow the creation of the GFM. This information is managed, updated, and maintained via CH2M HILL Plateau Remediation Company under the Environmental Model Management Archive.

58 GEOSCIENCES↗

FY2021 Status Report on the Computing Systems for the Yucca Mountain Project TSPA-LA Models and Testing of Selected Process Models

Sandia National Laboratories continued evaluation of the total system performance assessment (TSPA) for License Application (LA) computing systems for the previously considered Yucca Mountain Project (YMP). This was done to maintain the operational readiness of the computing infrastructure (computer hardware and software) and knowledge capability for total system performance assessment) type analysis, as directed by the National Nuclear Security Administration (NNSA), DOE 2010. The FY21 task included continued operation of the cluster; maintenance of the TSPA-LA models (with GoldSim 9.60.300); continued assessment of the status of the Infiltration Model; (a process model that feeds the TSP -LA) and preliminary assessments of the Unsaturated Zone Flow Model and the Saturated Zone Flow and Transport Model Abstraction (process models that feed the TSPA-LA). The 2014 cluster and supporting software systems are currently fully operational to support TSPA-LA type analyses.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Model-Form Epistemic Uncertainty Quantification for Modeling with Differential Equations: Application to Epidemiology

Modeling real-world phenomena to any degree of accuracy is a challenge that the scientific research community has navigated since its foundation. Lack of information and limited computational and observational resources necessitate modeling assumptions which, when invalid, lead to model-form error (MFE). The work reported herein explored a novel method to represent model-form uncertainty (MFU) that combines Bayesian statistics with the emerging field of universal differential equations (UDEs). The fundamental principle behind UDEs is simple: use known equational forms that govern a dynamical system when you have them; then incorporate data-driven approaches – in this case neural networks (NNs) – embedded within the governing equations to learn the interacting terms that were underrepresented. Utilizing epidemiology as our motivating exemplar, this report will highlight the challenges of modeling novel infectious diseases while introducing ways to incorporate NN approximations to MFE. Prior to embarking on a Bayesian calibration, we first explored methods to augment the standard (non-Bayesian) UDE training procedure to account for uncertainty and increase robustness of training. In addition, it is often the case that uncertainty in observations is significant; this may be due to randomness or lack of precision in the measurement process. This uncertainty typically manifests as “noisy” observations which deviate from a true underlying signal. To account for such variability, the NN approximation to MFE is endowed with a probabilistic representation and is updated using available observational data in a Bayesian framework. By representing the MFU explicitly and deploying an embedded, data-driven model, this approach enables an agile, expressive, and interpretable method for representing MFU. In this report we will provide evidence that Bayesian UDEs show promise as a novel framework for any science-based, data-driven MFU representation; while emphasizing that significant advances must be made in the calibration of Bayesian NNs to ensure a robust calibration procedure.

97 MATHEMATICS AND COMPUTING↗

A Multi-Model, Multi-Scale Research Program in Stressors, Responses, and Coupled Systems Dynamics at the Energy-Water-Land Nexus and for Concentrated, Interdependent Infrastructures: Toward Next Generation Capabilities in Integrated Impacts, Adaptation, and Vulnerability (I-IAV) Modeling and a Community of Practice

The goal of this research program was to build a next generation integrated suite of science-driven modeling and analytic capabilities, and a more expanded and connected community of practice, for analyses of the stressors, impacts, adaptations and vulnerabilities of global and regional change. The emphasis was on understanding energy-water-land interactions and feedbacks and interdependent infrastructures at appropriate regional and temporal scales. Although the scope spans many complex facets of data, modeling, and analysis, as well as scales appropriate for integrated impacts and adaptation research, the focus of this effort was the development of multi-model, multi-scale capabilities spanning the domains of Multi-Sector Dynamics (MSD) models; Impact, Adaptation, and Vulnerability (IAV) models; and Earth System Models (ESMs).

54 ENVIRONMENTAL SCIENCES↗

Modeling 2p-2h Interactions in Neutrino-Nucleus Scattering: A Comparative Overview of the Valencia and SuSAv2-MEC Models

2p-2h interactions are crucial for describing ν-A scattering, especially in the "dip region" between the Quasi-Elastic (QE) peak and Δ-resonance production. They involve the ejection of two nucleons, leading to two holes in the nuclear ground state. The underlying nuclear dynamics involve short-range correlations (SRC), long-range correlations (Random Phase Approximation, RPA), and the interplay of one- and two-body currents. Several models attempt to describe 2p-2h interactions, including Microscopic Models (e.g., Valencia), Scaling-Based Models (e.g., SuSAv2-MEC), and Transport Models (e.g., GiBUU). Despite progress, discrepancies exist between models and with data. I Will give an overview of 2p2h from the perspective of two Models Valencia and SuSAv2-MEC.

Hassinin, Karim [Houston U.]↗

The Aerosol Model Benchmarking Repository: A toolkit for model intercomparison

The Aerosol Model Benchmarking Repository and Standards (AMBRS) project was initiated to provide tools and to establish community standards for benchmarking aerosol models. This report describes a set of open-source tools for building, running, and analyzing aerosol box model simulations in a standardized framework. The framework consists of three core components: AMBuilder, a CMake-based build system that compiles supported models consistently; AMBRS, a Python module that defines unified numerical experiments and executes them with aligned inputs; and PyParticle, an aerosol analysis package that standardizes output, computes diagnostics, and visualizes simulation results. Together, these tools enable reproducible intercomparison of aerosol schemes and support process-level evaluation of how model simplifications affect predictions of size distributions, cloud condensation nuclei activity, and other relevant properties relevant for the Earth-Energy system. Beyond its role in benchmarking, AMBRS provides a platform for studying aerosol processes across scales and can be used to generate training data for AI/ML applications in support of a broader hierarchical aerosol modeling strategy.

54 ENVIRONMENTAL SCIENCES↗

Fracture Analysis of Cohesive Zone Models for Modeling Residual Stress Induced Delamination in Composite Structures

A fracture study of coupon-scale composite cylinders with embedded defects was conducted with an objective to assess and validate a modeling approach using two available cohesive material models. The study included experimental and simulation evaluations of initiation of crack growth and progression. Interrupted thermal experiments used acoustic emissions monitoring to identify the onset of crack progression during each cooling interval and ultrasonic scanning provided images of defect growth. Verification, validation, and uncertainty quantification (VVUQ) processes were performed in the assessment of the simulation predicted temperature at which crack propagation begins (quantity of interest). The Sobol sensitivity analysis identified the hoop direction elastic modulus in the carbon fiber reinforced polymer (CFRP) plies as the most influential parameter for simulations using both cohesive models, accounting for at least 70% of the variation in the temperature at crack propagation. The UQ temperature range for the Tvergaard-Hutchinson model was higher (more conservative) than the experimental acoustic measurement indicators of crack progression, while the temperature range for the Thouless-Parmigiani model enveloped the experimental data points for the primary defect size of 0.75 x 1 in. The simulations could not capture the stable crack growth indicated in the experiments. This is likely due to the models’ inability to represent anisotropic fracture toughness attributed to the structure of the orthotropic fiber weave in a woven composite laminate.

42 ENGINEERING↗

Global Sensitivity Analysis and Estimation of Model Error, Toward Uncertainty Quantification in Scramjet Computations

The development of scramjet engines is an important research area for advancing hypersonic and orbital flights. Progress toward optimal engine designs requires accurate flow simulations together with uncertainty quantification. However, performing uncertainty quantification for scramjet simulations is challenging due to the large number of uncertain parameters involved and the high computational cost of flow simulations. These difficulties are addressed in this paper by developing practical uncertainty quantification algorithms and computational methods, and deploying them in the current study to large-eddy simulations of a jet in crossflow inside a simplified HIFiRE Direct Connect Rig scramjet combustor. First, global sensitivity analysis is conducted to identify influential uncertain input parameters, which can help reduce the system’s stochastic dimension. Second, because models of different fidelity are used in the overall uncertainty quantification assessment, a framework for quantifying and propagating the uncertainty due to model error is presented. In conclusion, these methods are demonstrated on a nonreacting jet-in-crossflow test problem in a simplified scramjet geometry, with parameter space up to 24 dimensions, using static and dynamic treatments of the turbulence subgrid model, and with two-dimensional and three-dimensional geometries.

33 ADVANCED PROPULSION SYSTEMS↗

Modelling Hyperglycaemia in an Epithelial Membrane Model: Biophysical Characterisation

Biomimetic models are valuable platforms to improve our knowledge on the molecular mechanisms governing membrane-driven processes in (patho)physiological conditions, including membrane permeability, transport, and fusion. However, current membrane models are over simplistic and do not include the membrane’s lipid remodelling in response to extracellular stimuli. Our study describes the synthesis of glycated dimyristoyl-phosphatidylethanolamine (DMPE-glyc), which was structurally characterised by mass spectrometry (ESI-MS) and quantified by NMR spectroscopy to be further incorporated in a complex phospholipid (PL) membrane model enriched in cholesterol (Chol) and (glyco)sphingolipids (GSL) designed to mimic epithelial membranes (PL/Chol/GSL) under hyperglycaemia conditions. Characterisation of synthesised DMPE-glyc adducts by tandem mass spectrometry (ESI-MS/MS) show that synthetic DMPE-glyc adducts correspond to Amadori products and quantification by 1H NMR spectroscopy show that the yield of glycation reaction was 8%. The biophysical characterisation of the epithelial membrane model shows that excess glucose alters the thermotropic behaviour and fluidity of epithelial membrane models likely to impact permeability of solutes. The epithelial membrane models developed to mimic normo- and hyperglycaemic scenarios are the basis to investigate (poly)phenol-lipid and drug–membrane interactions crucial in nutrition, pharmaceutics, structural biochemistry, and medicinal chemistry.

Reis, Ana (ORCID:0000000245379902)↗

Modeling viscoelasticity through spring–dashpot models in intermittent-contact atomic force microscopy

We examine different approaches to model viscoelasticity within atomic force microscopy (AFM) simulation. Our study ranges from very simple linear spring–dashpot models to more sophisticated nonlinear systems that are able to reproduce fundamental properties of viscoelastic surfaces, including creep, stress relaxation and the presence of multiple relaxation times. Some of the models examined have been previously used in AFM simulation, but their applicability to different situations has not yet been examined in detail. The behavior of each model is analyzed here in terms of force–distance curves, dissipated energy and any inherent unphysical artifacts. We focus in this paper on single-eigenmode tip–sample impacts, but the models and results can also be useful in the context of multifrequency AFM, in which the tip trajectories are very complex and there is a wider range of sample deformation frequencies (descriptions of tip–sample model behaviors in the context of multifrequency AFM require detailed studies and are beyond the scope of this work).

36 MATERIALS SCIENCE↗

PROPERTIES OF CARBON–OXYGEN WHITE DWARFS FROM MONTE CARLO STELLAR MODELS

We investigate properties of carbon–oxygen white dwarfs with respect to the composite uncertainties in the reaction rates using the stellar evolution toolkit, Modules for Experiments in Stellar Astrophysics (MESA) and the probability density functions in the reaction rate library STARLIB. These are the first Monte Carlo stellar evolution studies that use complete stellar models. Focusing on 3 M{sub ⊙} models evolved from the pre main-sequence to the first thermal pulse, we survey the remnant core mass, composition, and structure properties as a function of 26 STARLIB reaction rates covering hydrogen and helium burning using a Principal Component Analysis and Spearman Rank-Order Correlation. Relative to the arithmetic mean value, we find the width of the 95% confidence interval to be ΔM{sub 1TP} ≈ 0.019 M{sub ⊙} for the core mass at the first thermal pulse, Δt{sub 1TP} ≈ 12.50 Myr for the age, Δlog(T{sub c}/K) ≈ 0.013 for the central temperature, Δlog(ρ{sub c}/g cm{sup −3}) ≈ 0.060 for the central density, ΔY{sub e,c} ≈ 2.6 × 10{sup −5} for the central electron fraction, ΔX{sub c}({sup 22}Ne) ≈ 5.8 × 10{sup −4}, ΔX{sub c}({sup 12}C) ≈ 0.392, and ΔX{sub c}({sup 16}O) ≈ 0.392. Uncertainties in the experimental {sup 12}C(α,γ){sup 16}O, triple-α, and {sup 14}N(p,γ){sup 15}O reaction rates dominate these variations. We also consider a grid of 1–6 M{sub ⊙} models evolved from the pre main-sequence to the final white dwarf to probe the sensitivity of the initial–final mass relation to experimental uncertainties in the hydrogen and helium reaction rates.

79 ASTRONOMY AND ASTROPHYSICS↗

UNCERTAINTIES IN GALACTIC CHEMICAL EVOLUTION MODELS

We use a simple one-zone galactic chemical evolution model to quantify the uncertainties generated by the input parameters in numerical predictions for a galaxy with properties similar to those of the Milky Way. We compiled several studies from the literature to gather the current constraints for our simulations regarding the typical value and uncertainty of the following seven basic parameters: the lower and upper mass limits of the stellar initial mass function (IMF), the slope of the high-mass end of the stellar IMF, the slope of the delay-time distribution function of Type Ia supernovae (SNe Ia), the number of SNe Ia per M {sub ⊙} formed, the total stellar mass formed, and the final mass of gas. We derived a probability distribution function to express the range of likely values for every parameter, which were then included in a Monte Carlo code to run several hundred simulations with randomly selected input parameters. This approach enables us to analyze the predicted chemical evolution of 16 elements in a statistical manner by identifying the most probable solutions, along with their 68% and 95% confidence levels. Our results show that the overall uncertainties are shaped by several input parameters that individually contribute at different metallicities, and thus at different galactic ages. The level of uncertainty then depends on the metallicity and is different from one element to another. Among the seven input parameters considered in this work, the slope of the IMF and the number of SNe Ia are currently the two main sources of uncertainty. The thicknesses of the uncertainty bands bounded by the 68% and 95% confidence levels are generally within 0.3 and 0.6 dex, respectively. When looking at the evolution of individual elements as a function of galactic age instead of metallicity, those same thicknesses range from 0.1 to 0.6 dex for the 68% confidence levels and from 0.3 to 1.0 dex for the 95% confidence levels. The uncertainty in our chemical evolution model does not include uncertainties relating to stellar yields, star formation and merger histories, and modeling assumptions.

79 ASTRONOMY AND ASTROPHYSICS↗