Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “parameter sampling”

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 451 records · Page 25

Cosmological constraints from H ii starburst galaxy, quasar angular size, and other measurements

ABSTRACT We compare the constraints from two (2019 and 2021) compilations of H ii starburst galaxy (H iiG) data and test the model independence of quasar (QSO) angular size data using six spatially flat and non-flat cosmological models. We find that the new 2021 compilation of H iiG data generally provides tighter constraints and prefers lower values of cosmological parameters than those from the 2019 H iiG data. QSO data by themselves give relatively model-independent constraints on the characteristic linear size, lm, of the QSOs within the sample. We also use Hubble parameter [H(z)], baryon acoustic oscillation (BAO), Pantheon Type Ia supernova (SN Ia) apparent magnitude (SN-Pantheon), and DES-3 yr binned SN Ia apparent magnitude (SN-DES) measurements to perform joint analyses with H iiG and QSO angular size data, since their constraints are not mutually inconsistent within the six cosmological models we study. A joint analysis of H(z), BAO, SN-Pantheon, SN-DES, QSO, and the newest compilation of H iiG data provides almost model-independent summary estimates of the Hubble constant, $H_0=69.7\pm 1.2\ \rm {km\,s^{-1}\,Mpc^{-1}}$, the non-relativistic matter density parameter, $\Omega _{\rm m_0}=0.293\pm 0.021$, and lm = 10.93 ± 0.25 pc.

79 ASTRONOMY AND ASTROPHYSICS↗

Application of LANDSAT to the surveillance and control of lake eutrophication in the Great Lakes basin

The author has identified the following significant results. Computer techniques were developed for mapping water quality parameters from LANDSAT data, using surface samples collected in an ongoing survey of water quality in Saginaw Bay. Chemical and biological parameters were measured on 31 July 1975 at 16 bay stations in concert with the LANDSAT overflight. Application of stepwise linear regression bands to nine of these parameters and corresponding LANDSAT measurements for bands 4 and 5 only resulted in regression correlation coefficients that varied from 0.94 for temperature to 0.73 for Secchi depth. Regression equations expressed with the pair of bands 4 and 5, rather than the ratio band 4/band 5, provided higher correlation coefficients for all the water quality parameters studied (temperature, Secchi depth, chloride, conductivity, total kjeldahl nitrogen, total phosphorus, chlorophyll a, total solids, and suspended solids).

Rogers, R. H.↗

Linear system identification via backward-time observer models

Presented here is an algorithm to compute the Markov parameters of a backward-time observer for a backward-time model from experimental input and output data. The backward-time observer Markov parameters are decomposed to obtain the backward-time system Markov parameters (backward-time pulse response samples) for the backward-time system identification. The identified backward-time system Markov parameters are used in the Eigensystem Realization Algorithm to identify a backward-time state-space model, which can be easily converted to the usual forward-time representation. If one reverses time in the model to be identified, what were damped true system modes become modes with negative damping, growing as the reversed time increases. On the other hand, the noise modes in the identification still maintain the property that they are stable. The shift from positive damping to negative damping of the true system modes allows one to distinguish these modes from noise modes. Experimental results are given to illustrate when and to what extent this concept works.

Juang, Jer-Nan↗

Lessons from Climate Modeling on the Design and Use of Ensembles for Crop Modeling

Working with ensembles of crop models is a recent but important development in crop modeling which promises to lead to better uncertainty estimates for model projections and predictions, better predictions using the ensemble mean or median, and closer collaboration within the modeling community. There are numerous open questions about the best way to create and analyze such ensembles. Much can be learned from the field of climate modeling, given its much longer experience with ensembles. We draw on that experience to identify questions and make propositions that should help make ensemble modeling with crop models more rigorous and informative. The propositions include defining criteria for acceptance of models in a crop MME, exploring criteria for evaluating the degree of relatedness of models in a MME, studying the effect of number of models in the ensemble, development of a statistical model of model sampling, creation of a repository for MME results, studies of possible differential weighting of models in an ensemble, creation of single model ensembles based on sampling from the uncertainty distribution of parameter values or inputs specifically oriented toward uncertainty estimation, the creation of super ensembles that sample more than one source of uncertainty, the analysis of super ensemble results to obtain information on total uncertainty and the separate contributions of different sources of uncertainty and finally further investigation of the use of the multi-model mean or median as a predictor.

Model ensembles↗

Simulation of Thematic Mapper performance as a function of sensor scanning parameters

The investigation and results of the Thematic Mapper Instrument Performance Study are described. The Thematic Mapper is the advanced multispectral scanner initially planned for the Earth Observation Satellite and now planned for LANDSAT D. The use of existing digital airborne scanner data obtained with the Modular Multispectral Scanner (M2S) at Bendix provided an opportunity to simulate the effects of variation of design parameters of the Thematic Mapper. Analysis and processing of this data on the Bendix Multispectral Data Analysis System were used to empirically determine categorization performance on data generated with variations of the sampling period and scan overlap parameters of the Thematic Mapper. The Bendix M2S data, with a 2.5 milliradian instantaneous field of view and a spatial resolution (pixel size) of 10-m from 13,000 ft altitude, allowed a direct simulation of Thematic Mapper data with a 30-m resolution. The flight data chosen were obtained on 30 June 1973 over agricultural test sites in Indiana.

Johnson, R. H.↗

Synchrotron self-Compton in a radiative-adiabatic fireball scenario: modelling the multiwavelength observations in some Fermi /LAT bursts

Energetic GeV photons expected from the closest and the most energetic Gamma-ray bursts (GRBs) provide a unique opportunity to study the very-high-energy emission as well as the possible correlations with lower energy bands in realistic GRB afterglow models. In the standard GRB afterglow model, the relativistic homogeneous shock is usually considered to be fully adiabatic, however, it could be partially radiative. Based on the external forward-shock scenario in both stellar wind and constant-density medium, we present a radiative-adiabatic analytical model of the synchrotron self-Compton (SSC) and synchrotron processes considering an electron energy distribution with a power-law index of $1\lt p\lt 2$ and $2\le p$. We show that the SSC scenario plays a relevant role in the radiative parameter $\epsilon$, leading to a prolonged evolution during the slow cooling regime. In a particular case, we derive the Fermi/LAT light curves together with the photons with energies $\ge 100$ MeV in a sample of nine bursts from the second Fermi/LAT GRB catalogue that exhibited temporal and spectral indices with $\gtrsim 1.5$ and $\approx 2$, respectively. These events can hardly be described with closure relations of the standard synchrotron afterglow model, and also exhibit energetic photons above the synchrotron limit. We have modelled the multiwavelength observations of our sample to constrain the microphysical parameters, the circumburst density, the bulk Lorentz factor, and the mechanism responsible for explaining the energetic GeV photons.

acceleration of particles↗

Data-Driven Tailoring Optimization of Thermoset Polymers Using Ultrasonics and Machine Learning

Thermoset polymers are highly demanded for their structural robustness, thermal stability, and chemical resistance. Tailoring the properties of these polymers for high-performance applications is often preferred to designing brand-new polymers. However, the traditional destructive techniques used to characterize their properties as a function of manufacturing parameters are expensive and time-consuming. A novel non-destructive, data-driven method leveraging ultrasonics and machine learning techniques to tailor the properties of thermosets as a function of the manufacturing parameters is demonstrated. Thermoset epoxy samples with varying curing temperatures (15–40 °C) and curing agent amounts (±40%) were manufactured and tested. Their curing kinetics were monitored by determining the sound speed in the material in real time, while the longitudinal modulus of the samples was determined post-cure. Machine learning models were developed using a k-nearest neighbors algorithm. These models were implemented to predict the curing and final elastic properties using the manufacturing parameters, i.e., stoichiometry and curing temperature, and vice versa. Understanding and modeling how these parameters affect the cure kinetics and final properties will allow for efficient and reliable optimization of thermoset tailoring and manufacturing.

36 MATERIALS SCIENCE↗

Deep Neural Network Informed Markov Chain Monte Carlo Methods

In subsurface flow modeling, quantifying the uncertainty of model parameters and the corresponding uncertainly on output quantities is a crucial task for groundwater management. Markov chain Monte Carlo (MCMC) methods can take advantage of observed data to estimate parameters in a Bayesian setting. However, MCMC can be slow to converge and produce highly correlated samples when the dimensions of the parameters is high. Using gradients for the posterior distribution can help samplers explore the parameter space more efficiently, but obtaining gradients can be computationally challenging.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

User's guide to a system of finite-element supersonic panel flutter programs

The utilization and operation of a set of six computer programs for the prediction of panel flutter at supersonic speeds by finite element methods are described. The programs run individually to determine the flutter behavior of a flat panel where the finite elements which model the panel each have four degrees of freedom (DOF), a curved panel where the finite elements each have four DOF, and a curved panel where the finite elements each have six DOF. The panels are assumed to be of infinite aspect ratio and are subjected to either simply-supported or clamped boundary conditions. The aerodynamics used by these programs are based on piston theory. Application of the program is illustrated by sample cases where the number of beam finite elements equals four, the in-plane tension parameter is 0.0, the maximum camber to panel length ratio for a curved panel case is 0.05, and the Mach number is 2.0. This memorandum provides a user's guide for these programs, describes the parameters that are used, and contains sample output from each of the programs.

Woolley, Christine L.↗

A Simulation Atlas of Tidal Features in Galaxies

Detailed simulations of tidally induced structure in disk galaxies have either concentrated on specific systems or consisted of a few encounters with relatively small numbers of particles and no self-gravity. Observers need a 'dictionary' of simulations that covers many encounter parameters with fine morphological resolution and includes effects of self-gravitation. Observers can then search the dictionary for the parameters that best match a particular observed morphology. Alternatively, the dictionary can be used with observational samples for statistical studies of system parameters. To fill this need, we present a survey of model tidal encounters using a self-gravitating, 180,000 particle, two-component ('stars' and 'gas') disk. A wide variety of fascinating morphologies results. There are 86 different encounters that vary orbit tilt, perigalacticon distance, galaxy to companion mass ratio, and the amount of halo dark matter relative to the disk. For morphological comparisons, over 1700 images of the entire survey are available in video form. While there is a rich variety of tidal structure covering much of this parameter space, some general patterns may be remarked. There is a strong orbital inclination dependence of the symmetry of tidal patterns, most symmetric for planar orbits and nearly one-sided for polar encounters. Retrograde encounters produce only broad fanlike global patterns, but rich small-scale internal structure. In both kinds of encounter, our numerical resolution allows us to track internal spiral structure driven by the outer material arms, especially in the lighter halo simulations. We note also that polar encounters generate series of expanding, essentially non-rotating loops resembling shell structures in some respects.

Howard, Sethanne↗

Towards a Robust Sampling Approach: A Computational Review and Design

As pointed out in several other works, the estimation of the reliability of the electrical grid can not be conducted without the estimation of the stochastic phenomena of electicity demand and electricity production by variable renewable sources. Therefore, sampling procedures have become integral in the design of engineering structures and analysis. Commonly referred to as Monte Carlo uncertainty analysis or integration, the general objective of these procedures is to establish specifics about the uncertainty of an output characteristic of such an engineering system, given uncertainty about its input characteristics. These sampling procedures are applied in a context in which establishing such specifics cannot be performed through other means. Variance-reduction techniques are designed to lessen the variability among estimators to estimate statistics of those output uncertainties. Importance-sampling techniques, on the other hand, are designed to reduce the number of samples needed to estimate a particular statistic—e.g., a tail probability. The combination of these approaches can reduce the computational burden considerably for a particular estimator and statistic. Importance sampling—geared and designed as it is toward improving a particular estimator—suffers, unfortunately, from the unintended consequence of reducing the performance of other estimators in terms of their variance. The objective of this paper is to offer an alternative sampling procedure where this variance does not grow unacceptably large for a suite of estimators. Moreover, it is anticipated that, with additional knowledge of how an engineering output characteristic responds to its input characteristics, tuning parameters of the input’s sampling procedure can be set to improve the output characteristic’s estimation. This work proves the effectiveness of the suggested alternative approach. Such positive outcome will lead to a decrease of the computational burden of performing stochastic optimization of integrated energy systems (e.g., components dispatch, and portfolio composition). In particular, capturing the contribution to the overall system cost of rare and unlikely events and patterns of the electricity demand and production will become less computationally expensive. This is due to the fact that the approach demonstrated here will allow the sampling of those rare occurrences more frequently without misrepresenting their probabilistic impacts.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Uncertainty quantification and sensitivity analysis of a nuclear thermal propulsion reactor startup sequence

The research presented in this article describes progress in applying stochastic methods, uncertainty quantification, parametric studies, and variance-based sensitivity analysis (also known as Sobol sensitivity analysis) to a full-core model of a nuclear thermal propulsion (NTP) system simulated via the radiation transport code Griffin to simulate neutronics. Our goal is to develop a reduced-order (surrogate) model that can be rapidly sampled with perturbations to multiple input parameters. In this NTP system, reactivity and power feedback affect the rotation of control drums (CDs), which is itself controlled by a hybrid proportional-integral-derivative (PID) controller actuated by the power demand and reactivity feedback from the numerical model. This model uses reactor kinetic feedback (mean generation time [Λ] and effective delayed neutron fraction [ β eff ] from a transient Griffin simulation executed via Griffin’s improved quasi-static solver to provide the kinetic parameters) as inputs to functions that control the CD rotation angle. By investigating numerous stochastic approaches, we developed a dual-purpose surrogate model of the NTP system, using polynomial regression in the Multiphysics Object-Oriented Simulation Environment (MOOSE) Stochastic Tools Module (STM). The trained model can be rapidly sampled while simultaneously perturbing various input parameters, such as coefficients on the PID control or temperature (directly affecting the neutron cross section). The surrogate model delivers accurate (within 5%) results at speeds orders of magnitude faster (minutes, not days of computational time) than the base model. Once the surrogate model has been trained, distributions of the uncertain parameters can be changed at will to investigate the effects of perturbing multiple inputs as well as the effects of these inputs on the model output. For example, coefficients used in the PID control system may vary due to some type of physical interference, or uncertainty may exist in the temperature of the neutron cross sections in various regions of the reactor. A distribution can be placed on these parameters, and operational boundaries can be determined. The goal of this work is to support development of an advanced control system for operating CDs in a functioning NTP system. This work is a scoping study of the MOOSE STM.

21 - SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLAN↗

Advances in Spectro‐Microscopy Methods and their Applications in the Characterization of Perovskite Materials

Perovskite materials are promising contenders as the active layer in light-harvesting and light-emitting applications if their long-term stability can be sufficiently increased. Chemical and structural engineering are shown to enhance long-term stability, but the increased complexity of the material system also leads to inhomogeneous functional properties across various length scales. Thus, scanning probe and high-resolution microscopy characterization techniques are needed to reveal the role of local defects and the results promise to act as the foundation for future device improvements. Here, a look at the parameter space: technique-specific sample penetration depth versus probe size highlights a gap in current methods. High spatial resolution combined with a deep penetration depth is not yet achievable. However, multimodal measurement technique may be the key to covering this parameter space. In this perspective, current advanced spectro-microscopy methods which have been applied to perovskite materials are highlighted.

Luo, Yanqi↗

CatMass : software for calculating optimal sample masses for X-ray absorption spectroscopy experiments involving complex sample compositions

This paper presents software for calculating the optimal mass of samples with complex compositions ( e.g. supported metal catalysts) for X-ray absorption spectroscopy (XAS) and scattering measurements. The ability to calculate the sample mass and other relevant parameters needed for an XAS measurement allows experimentalists to be better prepared in terms of detector selection, energy range of scan and overall time needed to complete the measurement, thus increasing efficiency. CatMass builds on existing sample mass calculators allowing users to determine the optimum sample preparation, collection geometry, usable energy range for a scan and approximate edge step of the absorption event. Visualization tools present the absorption calculation results in a format familiar to XAS experimentalists, with the added ability to save calculations and plots for future reference or recalculation. CatMass is a program broadly applicable in catalysis and is helpful for users with complex samples due to composition/stoichiometry or multiple competing elements.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Method and apparatus for a Fabry-Perot multiple beam fringe sensor

A method and, in one embodiment of the invention, the resulting apparatus for implementing a unique multiple beam fringe sensor that is adapted to be interfaced with a low cost, compact fiber optic transmission system in order to provide an accurate digital representation of a physical parameter (e.g. temperature) of a remote sample. The sensor is fabricated so as to include a Fabry-Perot gap formed between the ends of two mated optical fibers. By examining the optical characteristics of light that is transmitted through the Fabry-Perot sensor gap, an indication of gap width can be ascertained. Accordingly, a change in Fabry-Perot sensor gap width is related to a change in the particular physical parameter to be measured. In another embodiment of the invention, a second unique multiple beam fringe sensor having a Fabry-Perot gap is disclosed that is also adapted to provide an accurate digital representation of a physical parameter (e.g. temperature) of a remote sample. The sensor may be fabricated in two segments. A fiber containing segment includes each of a driving optical fiber for supplying incident light signals to the Fabry-Perot gap and a sensing optical fiber for receiving output light signals that have been transmitted twice through the Fabry-Perot gap, the optical characteristics of which output signals provide an indication of the parameter to be sensed. A transducer segment includes the Fabry-Perot gap formed therein and means responsive to the physical parameter for changing the width of the Fabry-Perot gap and, accordingly, the optical characteristics of the light signals passing therethrough.

James, Kenneth A.↗

Differentiable Multiphysics Codes: A Breakthrough Technology for Simulation and Computing

This document summarizes the findings of a strategic planning exercise commissioned by the Weapons Simulation and Computing, Computational Physics (WSC/CP) program at the Lawrence Livermore National Laboratory (LLNL) in FY24. During the year, the committee met with multiple stakeholder communities to gather input, opinions, suggestions and concerns which have been incorporated throughout this document. The key findings from this exercise are summarized: • The development of multiphysics modelling and simulation (mod/sim) codes and software technologies, their deployment on exascale compute platforms, and their broad adoption across the NNSA is a major success of the Advanced Simulation and Computing (ASC) program and the Exascale Computing Project (ECP). Sustained investment in these core technologies is essential. • Today’s state of the art involves running ensembles of O(100K) simulations to perform uncertainty quantification (UQ) and design studies using multiple statistical methods such as Bayesian optimization to understand sensitivities of our models and explore parameterized design spaces. Even with exascale computing, we are practically limited to O(10) parameters in these studies since the number of simulations required to sample the space scales exponentially with the number of design parameters. • The data from these simulation ensembles is increasingly being used to train machine learned (ML) surrogates (or reduced order models, ROMs) which can then be used for optimization or real time design exploration. However, the trained surrogates are still limited in the number of parameters they can represent due to the sampling limitations previously noted. • Augmenting our suite of integrated multiphysics simulation codes, both current and emerging, with the ability to compute gradients (solution derivatives) of arbitrary simulation outputs with respect to (some or all) simulation inputs would be a breakthrough technology, opening the door to a new era of efficient and automated inverse design based on verified and validated mod/sim capabilities. • This capability, which we refer to as differentiable multiphysics codes (DMCs), would revolutionize both UQ and optimization studies by breaking the curse of dimensionality that presently limits our “gradient-free” ensemble based computing approach. A similar breakthrough occurred in the AI/ML community once the ability to compute gradients of arbitrary loss functions using back-propagation became commonplace. Gradient information from the multiphysics codes can also be used to dramatically improve the efficiency and scale of training of ML/ROM surrogates for rapid assessments. • Achieving this in our suite of codes will be a grand challenge, similar to the amount of effort that was required to transition from CPU to GPU computing. It will require buy-in from the entire WSC/CP program and beyond, including all integrated codes, physics and engineering models, third-party library dependencies and performance portability abstractions. It will also require investment in research and development of numerical methods for computing adjoints of coupled physics across multiple adaptively refined moving meshes and of stochastic (Monte Carlo) and mesh free (SPH) methods. • New software and numerical techniques, largely pioneered by the AI/ML community, make this feasible. Chief among these is automatic differentiation (AD), the ability to employ AD at point-wise locations in a physics calculation (instead of traditional black-box approaches) and the ability to perform “back-propagation in time” (or reverse mode AD) for non-linear partial differential equations (PDEs). Fundamentally, the conclusion of this strategic planning exercise is that the time is right to undertake a large scale effort in WSC, centered on the existing integrated codes, to continue the natural evolution of mod/sim in the age of AI/ML. Instead of attempting to replace mod/sim with purely data driven AI/ML models, we believe the key to success is to integrate AI/ML by building on top of the decades of hard-won knowledge and the verified/validated multiphysics modelling capability that is the hallmark of the ASC program.

97 MATHEMATICS AND COMPUTING↗

The rotation of small asteroids

The Binzel and Mulholland (1983) sample of photoelectrically determined rotational parameters for 17 main belt asteroids with diameters less than 30 km are compared with previous observations of asteroids of that size range. Rigorous statistical tests are applied to investigate bias effects and quantify results on asteroid rotation. The samples are described and compared for rotational frequency, rotational amplitude, frequency distribution, and diameter and frequency dependence. It is concluded that the observed rotational frequency distribution can be acceptably fit by two Maxwellian distributions, which is consistent with the hypothesis that there are separate populations of slow and fast rotating asteroids. The frequency distributions of main belt asteroids less than 15 km in diameter and earth and Mars crossers do not differ significantly, but the larger mean lightcurve amplitude of these crossers is statistically significant. No significant diameter dependence on rotational frequency is found among the sample asteroids.

Binzel, R. P.↗

High-Throughput Microfluidics Platform for Intracellular Delivery and Sampling of Biomolecules from Live Cells

Nondestructive cell membrane permeabilization systems enable the intracellular delivery of exogenous biomolecules for cell engineering tasks as well as the temporal sampling of cytosolic contents from live cells for the analysis of dynamic processes. Here, in this paper, we report a microwell array format live-cell analysis device (LCAD) that can perform localized-electroporation induced membrane permeabilization, for cellular delivery or sampling, and directly interfaces with surface-based biosensors for analyzing the extracted contents. We demonstrate the capabilities of the LCAD via an automated high-throughput workflow for multimodal analysis of live-cell dynamics, consisting of quantitative measurements of enzyme activity using self-assembled monolayers for MALDI mass spectrometry (SAMDI) and deep-learning enhanced imaging and analysis. By combining a fabrication protocol that enables robust assembly and operation of multilayer devices with embedded gold electrodes and an automated imaging workflow, we successfully deliver functional molecules (plasmid and siRNA) into live cells at multiple time-points and track their effect on gene expression and cell morphology temporally. Furthermore, we report sampling performance enhancements, achieving saturation levels of protein tyrosine phosphatase activity measured from as few as 60 cells, and demonstrate control over the amount of sampled contents by optimization of electroporation parameters using a lumped model. Lastly, we investigate the implications of cell morphology on electroporation-induced sampling of fluorescent molecules using a deep-learning enhanced image analysis workflow.

59 BASIC BIOLOGICAL SCIENCES↗