Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “solution multiplicity”

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 217 records · Page 12
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

The Technical and Economic Potential of the H2@Scale Hydrogen Concept within the United States

The U.S. energy system is evolving as society and technologies change. Renewable electricity generation - especially from wind and solar - is growing rapidly, and alternative energy sources are being developed and implemented across the residential, commercial, transportation, and industrial sectors to take advantage of their cost, security, and health benefits. Systemic changes present numerous challenges to grid resiliency and energy affordability, creating a need for synergistic solutions that satisfy multiple applications while yielding system-wide cost and emissions benefits. One such solution is an integrated hydrogen energy system. This is the focus of H2@Scale - a U.S. Department of Energy (DOE) initiative led by the Office of Energy Efficiency and Renewable Energy’s Hydrogen and Fuel Technologies Office. H2@Scale brings together stakeholders to advance affordable hydrogen production, transport, storage, and utilization in multiple energy sectors. The H2@Scale concept involves hydrogen as an energy intermediate. Hydrogen can be produced from various conventional and renewable energy sources including as a responsive load on the electric grid. Hydrogen has many current applications and many more potential applications, such as energy for transportation - used directly in fuel cell electric vehicles (FCEVs), as a feedstock for synthetic fuels, and to upgrade oil and biomass - feedstock for industry (e.g., for ammonia production, metals refining, and other end uses), heat for industry and buildings, and electricity storage. Owing to its flexibility and fungibility, a hydrogen intermediate could link energy sources that have surplus availability to markets that require energy or chemical feedstocks, benefiting both. This document builds upon a growing body of analyses of hydrogen as an energy intermediate by reporting the results from our initial analysis of the potential impacts of the H2@Scale vision by the mid-21 st century for the 48 contiguous U.S. states. Previous estimates have been based on expert elicitation and focused on hydrogen demands. We build upon them, first, by estimating hydrogen’s serviceable consumption potential for possible hydrogen applications and the technical potential for producing hydrogen from various resources. We define the serviceable consumption potential as the quantity of hydrogen that would be consumed to serve the portion of the market that could be captured without considering economics (i.e., if the price of hydrogen were $0/kg over an extended period); thus, it can be considered an upper bound for the size of the market. We define the technical potential as the resource potential constrained by real-world geography and system performance, but not by economics. We then compare the cumulative serviceable consumption potential with the technical potential of a number of possible sources. Second, we estimate economic potential: the quantity of hydrogen at an equilibrium price at which suppliers are willing to sell and consumers are willing to buy the same quantity of hydrogen. We believe this method provides a deeper understanding than was available in the previous analyses. We develop economic potentials for multiple scenarios across various market and technology-advancement assumptions.

08 HYDROGEN↗
NASA NTRS1987

Uniformly-valid asymptotic solutions to the Orr-Sommerfeld equation using multiple scales

The Orr-Sommerfeld equation is studied using the method of multiple scales. The leading-order uniformly-valid asymptotic solutions are obtained without recourse to heuristic methods for two different types of velocity profiles. Solutions were obtained for arbitrary velocity profiles without inflection or critical points and for profiles with simultaneous inflection and critical points. These are shown to reproduce previously obtained solutions in the inner and outer regions.

Long, Lyle N.↗
NASA NTRS1971

Probing atmospheric water vapor profiles via multiple scattering of electromagnetic waves

A theoretical analysis on the multiple scattering of electromagnetic waves propagating in a finite inhomogeneous medium is presented and applied to the study of wave propagation in a clear atmosphere (fine weather conditions). It is shown that the analysis offers a method of synthesizing the water vapor density profile in a clear atmosphere by measuring the resultant reflections from the density profile at several different frequencies. It is also shown that the resultant reflection emerges as the consequence of multiple scattering of partial reflections from various parts of the inhomogeneous medium. The solutions of the multiple scattering approach are shown to be more accurate than those of the WKB approach, which neglects the multiple scattering effects.

Tang, C. C. H.↗
NASA NTRS1993

Multiple scattering of light in a spherical cometary atmosphere with an axisymmetric dust jet. II - Image simulation

A numerical solution for the multiple light scattering in spherical axisymmetric geometry is applied to the simulation of images of a coma as it would appear to a near-flying satellite such as Giotto. The appearance of symmetric comas and dust jets is examined in detail; the nucleus visibility is studied; the effect of forward scattering is considered; and single and multiple scattering effects are quantified. Attention is given to simulated images of a coma with a hollow cone of dust, as predicted by dust-gas hydrodynamic modeling. The cone's appearance is very similar to the northern area of activity on Comet Halley, observed by the Giotto HMC.

Chick, Kenneth M.↗
NASA NTRS1985

A FORTRAN computer code for calculating flows in multiple-blade-element cascades

A solution technique has been developed for solving the multiple-blade-element, surface-of-revolution, blade-to-blade flow problem in turbomachinery. The calculation solves approximate flow equations which include the effects of compressibility, radius change, blade-row rotation, and variable stream sheet thickness. An integral equation solution (i.e., panel method) is used to solve the equations. A description of the computer code and computer code input is given in this report.

Mcfarland, E. R.↗
NASA NTRS1980

The radiation of sound by the instability waves of a compressible plane turbulent shear layer

The problem of acoustic radiation generated by instability waves of a compressible plane turbulent shear layer is solved. The solution provided is valid up to the acoustic far-field region. It represents a significant improvement over the solution obtained by classical hydrodynamic-stability theory which is essentially a local solution with the acoustic radiation suppressed. The basic instability-wave solution which is valid in the shear layer and the near-field region is constructed in terms of an asymptotic expansion using the method of multiple scales. This solution accounts for the effects of the slightly divergent mean flow. It is shown that the multiple-scales asymptotic expansion is not uniformly valid far from the shear layer. Continuation of this solution into the entire upper half-plane is described. The extended solution enables the near- and far-field pressure fluctuations associated with the instability wave to be determined. Numerical results show that the directivity pattern of acoustic radiation into the stationary medium peaks at 20 degrees to the axis of the shear layer in the downstream direction for supersonic flows. This agrees qualitatively with the observed noise-directivity patterns of supersonic jets.

Tam, C. K. W.↗
NASA NTRS1995

A Green's function method for heavy ion beam transport

The use of Green's function has played a fundamental role in transport calculations for high-charge high-energy (HZE) ions. Two recent developments have greatly advanced the practical aspects of implementation of these methods. The first was the formulation of a closed-form solution as a multiple fragmentation perturbation series. The second was the effective summation of the closed-form solution through nonperturbative techniques. The nonperturbative methods have been recently extended to an inhomogeneous, two-layer transport media to simulate the lead scattering foil present in the Lawrence Berkeley Laboratories (LBL) biomedical beam line used for cancer therapy. Such inhomogeneous codes are necessary for astronaut shielding in space. The transport codes utilize the Langley Research Center atomic and nuclear database. Transport code and database evaluation are performed by comparison with experiments performed at the LBL Bevalac facility using 670 A MeV 20Ne and 600 A MeV 56Fe ion beams. The comparison with a time-of-flight and delta E detector measurement for the 20Ne beam and the plastic nuclear track detectors for 56Fe show agreement up to 35%-40% in water and aluminium targets, respectively.

NASA Discipline Radiation Health↗
NASA NTRS1974

Stability analysis of nonlinear systems by multiple time scaling

The asymptotic solution for the transient analysis of a general nonlinear system in the neighborhood of the stability boundary was obtained by using the multiple-time-scaling asymptotic-expansion method. The nonlinearities are assumed to be of algebraic nature. Terms of order epsilon to the 3rd power (where epsilon is the order of amplitude of the unknown) are included in the solution. The solution indicates that there is always a limit cycle which is stable (unstable) and exists above (below) the stability boundary if the nonlinear terms are stabilizing (destabilizing). Extension of the solution to include fifth order nonlinear terms is also presented. Comparisons with harmonic balance and with multiple-time-scaling solution of panel flutter equations are also included.

Morino, L.↗
DOE OSTI2021

Component-wise reduced order model lattice-type structure design

Lattice-type structures can provide a combination of stiffness with light weight that is desirable in a variety of applications. Design optimization of these structures must rely on approximations of the governing physics to render solution of a mathematical model feasible. In this paper, we propose a topology optimization (TO) formulation that approximates the governing physics using component-wise reduced order modeling as introduced in Huynh et al. (2013); Eftang and Patera (2013), which can reduce solution time by multiple orders of magnitude over a full-order finite element model while providing a relative error in the solution of 1%. In addition, the offline training data set from such component-wise models is reusable, allowing its application to many design problems for only the cost of a single offline training phase, and the component-wise method is nearly embarrassingly parallel. We also show how the parameterization chosen in our optimization allows a simplification of the component-wise reduced order model (CWROM) not noted in previous literature, for further speedup of the optimization process. Furthermore, the sensitivity of the compliance with respect to the particular parameterization is derived solely at the component level. In numerical examples, we demonstrate a 1000x speedup over a full-order FEM model with relative error of 1% and show minimum compliance designs for two different cantilever beam examples, one smaller and one larger. Finally, error bounds for the displacement field, compliance, and compliance sensitivity of the CWROM are derived.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2026

Field Validation of a Grid-Interactive Efficient Building Software Solution

The U.S. General Services Administration's (GSA's) Green Proving Ground (GPG) program, in partnership with the National Laboratory of the Rockies (NLR), completed a field study of a Grid-Interactive Efficient Buildings (GEB) software solution. The study focused on a single testbed facility to test the GEB functionality of the software solution, along with other features. The testbed facility - a courthouse - is a common building type in GSA's vast building portfolio, offering potentially impactful findings on a scalable level. The study evaluated Prescriptive Data's technology, Nantum OS, a connected building operating system ("GEB Solution") which aggregates multiple sources of previously siloed building data and combines that data with external sources, such as weather information or utility signals, into a single integrated platform. A GEB Solution is a type of Energy Management Information System (EMIS). EMIS is defined as a system of devices, data services, and software applications that communicates with any building system or third-party data source to aggregate and transform data into new capabilities to aid in the optimization of energy use at the building, campus, or agency level. This specific GEB Solution is an EMIS with ASO, automated system optimization, offering supervisory control of certain aspects of the Building Automation System (BAS). Multiple features were evaluated including, but not limited to, Continuous Demand Management to avoid setting new monthly kilowatt (kW) peaks, energy efficiency for reduction of kilowatt hours (kWh) and natural gas consumption, and automated demand response (ADR) for purposes of lowering demand during a utility called Demand Response (DR) event. The testbed facility was the Foley Federal Building and US Courthouse ("Foley Federal Building") located in Las Vegas, NV. This is a 209,496 sq. ft. building constructed in the 1960s with major renovations in 2004. The facility was a good candidate due to the large prevalence of office and courthouse spaces in the GSA portfolio of buildings. It also has many features which allow integration into and control of the building and a strong facilities team to assist with the study. Quantitative and qualitative performance objectives were developed using GSA's GPG GEB project template along with input from the vendor and building facility staff; these are outlined in Table 1. The quantitative performance objectives focused on continuous demand management, energy efficiency, and automated demand response. The qualitative performance objectives focused on the ease of installation and commissioning as well as the operability of the GEB solution. Other performance metrics that are reported on include carbon reduction, cost effectiveness, and occupant acceptance.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
NASA NTRS1995

Structure of the magnetopause inferred from one-dimensional hybrid simulations

The structure of the magnetopause is investigated by studying the interaction between two plasmas with solar wind and magnetospheric properties. Both Hall-MHD and hybrid (particle ions and fluid electrons) simulations are performed to compare and contrast the nature of the solutions in the fluid and the kinetic limits. It is shown that, in accordance with previous studies, the fluid solutions consist of multiple discontinuities and waves such as slow shocks and waves, as well as contact and rotational discontinuities. In contrast, the kinetic solutions consist of fewer discontinuities and include non-MHD boundaries. The difference between the two types of solutions are attributed to the absence of contact surfaces in collisionless plasmas and the possible Landau damping of slow waves. The kinetic solutions are found to be in a much better agreement with magnetopause observations, which have shown little evidence for the presence of slow shocks or contact discontinuities. The results of kinetic calculations suggest that the presence of a small but finite normal component of the magnetic field allows for the mixing of magnetosheath and magnetospheric plasmas and may, in part, be responsible for the formation of the boundary layer even during the periods of northward interplanetary magnetic field. These results also show that the necessary changes in the fields and the plasma parameters across the magnetopause do not occur simultaneously (i.e., in the same spatial location). As such, the magnetopause cannot be viewed as a single boundary but instead is a region which includes both abrupt and gradual changes in the fields and plasma parameters.

Omidi, N.↗
DOE OSTI2022

Potential of water quality wetlands to mitigate habitat losses from agricultural drainage modernization

Given widespread biodiversity declines, a growing global human population, and demands to improve water quality, there is an immediate need to explore land management solutions that support multiple ecosystem services. Agricultural water quality wetlands designed to provide both water quality benefits and wetland and grassland habitat are an emerging restoration solution that may reverse habitat declines in intensive agricultural areas. Installation of water quality wetlands in the Upper Midwest, USA, when considered alongside the repair and modification of aging agricultural tile drainage infrastructure, is a likely scenario that may mitigate nutrient pollution exported from agricultural systems and improve crop yields. The capacity of water quality wetlands to provide habitat within the wetland pool and the surrounding grassland is not well-studied, particularly with respect to potential habitat changes resulting from drainage infrastructure upgrades. For the current study, we produced spatially explicit models of 37 catchments distributed throughout an important region for agriculture and biodiversity, the Des Moines Lobe of Iowa. Four scenarios were considered - with and without improved drainage and with and without water quality wetlands - to estimate the net potential habitat implications of these scenarios for amphibians, grassland birds, and wild bees. Model results indicate that drainage modification alone will likely result in moderate direct losses of suitable amphibian habitat and large declines in overall habitat quality. However, inclusion of water quality wetlands at the catchment scale may mitigate these amphibian habitat losses while also increasing grassland bird and pollinator habitat. In conclusion, the impacts of water quality wetlands and drainage modernization on waterfowl in the region require additional study.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

Field testing and validation of a low-cost MPC for demand flexibility for grid-interactive K-12 schools

K-12 school buildings account for the highest energy consumption within the public sector. Implementing advanced HVAC controls in grid-interactive K-12 schools could bring substantial economic advantages and grid flexibility. Our previous study demonstrated that a low-cost model predictive control (MPC) solution, which coordinates multiple packaged units, can enable demand flexibility without major hardware upgrades. However, a significant gap remains between academic pilots and market-ready scalable solutions. This paper extends the previous single-site pilot to a multi-site demonstration involving three school campuses (95 total units) through a commercial technology transfer process. Addressing the challenge of verifying performance with sparse field data, we present a new statistical approach using Bayesian methods to estimate the MPC’s effect on peak demand. Unlike traditional methods, this approach robustly quantifies uncertainty in non-normal, limited datasets. The results confirm the solution’s replicability, achieving a 21.6–38.9% reduction in HVAC peak demand (10.8–22.1% at the site-level) with > 98% probability across diverse locations. Finally, we document critical barriers to scaling software-as-a-service (SaaS) solutions–such as API instability and diverse legacy systems–and offer practical strategies to accelerate the commercial adoption of grid-interactive efficient buildings.

Ham, Sang Woo↗
DOE OSTI2020

Learning and meta-learning of stochastic advection–diffusion–reaction systems from sparse measurements

Physics-informed neural networks (PINNs) were recently proposed in [18] as an alternative way to solve partial differential equations (PDEs). A neural network (NN) represents the solution, while a PDE-induced NN is coupled to the solution NN, and all differential operators are treated using automatic differentiation. Here, we first employ the standard PINN and a stochastic version, sPINN, to solve forward and inverse problems governed by a non-linear advection–diffusion–reaction (ADR) equation, assuming we have some sparse measurements of the concentration field at random or pre-selected locations. Subsequently, we attempt to optimise the hyper-parameters of sPINN by using the Bayesian optimisation method (meta-learning) and compare the results with the empirically selected hyper-parameters of sPINN. In particular, for the first part in solving the inverse deterministic ADR, we assume that we only have a few high-fidelity measurements, whereas the rest of the data is of lower fidelity. Hence, the PINN is trained using a composite multi-fidelity network, first introduced in [12], that learns the correlations between the multi-fidelity data and predicts the unknown values of diffusivity, transport velocity and two reaction constants as well as the concentration field. For the stochastic ADR, we employ a Karhunen–Loève (KL) expansion to represent the stochastic diffusivity, and arbitrary polynomial chaos (aPC) to represent the stochastic solution. Correspondingly, we design multiple NNs to represent the mean of the solution and learn each aPC mode separately, whereas we employ a separate NN to represent the mean of diffusivity and another NN to learn all modes of the KL expansion. For the inverse problem, in addition to stochastic diffusivity and concentration fields, we also aim to obtain the (unknown) deterministic values of transport velocity and reaction constants. The available data correspond to 7spatial points for the diffusivity and 20 space–time points for the solution, both sampled 2000 times. We obtain good accuracy for the deterministic parameters of the order of 1–2% and excellent accuracy for the mean and variance of the stochastic fields, better than three digits of accuracy. In the second part, we consider the previous stochastic inverse problem, and we use Bayesian optimisation to find five hyper-parameters of sPINN, namely the width, depth and learning rate of two NNs for learning the modes. Here, we obtain much deeper and wider optimal NNs compared to the manual tuning, leading to even better accuracy, i.e., errors less than 1% for the deterministic values, and about an order of magnitude less for the stochastic fields.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS2005

Sub-pixel Spatial Resolution Micro-roughness Measurements with Interlaced Stitching

In this paper we describe a method to increase the spatial resolution of surface micro-roughness measurements. As the surface specifications for precision optics become more demanding, the metrology instruments must cover a broad spatial frequency range. Generally, multiple instruments are used to cover the full range of the specifications. For example, a Fizeau interferometer would be used to test low spatial frequency surface errors, a white light interferometer would be used for mid spatial frequency errors, and an AFM would be used for high spatial frequency errors. For some precision optics, three instruments would be necessary. However, in many applications, an increase in the spatial resolution of a metrology instrument could negate the need for multiple instruments. A solution to increase the spatial resolution of micro-roughness measurements obtained using white light interferometry with "interlaced stitching" will be presented. In "interlaced stitching", multiple interferometric measurements are made as the optic under test (or the CCD array) is laterally shifted at sub-pixel increments. The measurements are then combined to construct a measurement with higher spatial resolution. Initial results obtained while implementing a similar process used to increase the spatial resolution of measurements made with a commercially available Fizeau interferometer will be presented.

Mooney, James T.↗
NASA NTRS1986

Provably conservative approximations to complex reliability models

Complex models can be the bases for derivation of provably conservative and optimistic reliability models that incorporate a reduced state space and fewer transitions; they accordingly possess solutions that are more cost-effective than those of the original complex models. Design space can thereby be extensively explored without incurring the expense of multiple complex model solutions. A conservative-optimistic pair of derived models produces a band that includes the solution to the complex model. Sensitivity analysis can be performed on this pair of models to determine those parameters of the original model that are most sensitive to change and therefore require further expense in obtaining tighter specifications.

Smotherman, M.↗
NASA NTRS2020

Improving Navigation Analysis with OD-D: The Visually Interactive Orbit Determination Dashboard

Orbit determination requires iterative analysis with the goal of converging on as accurate of a solution as possible, a process that can be time and labor intensive. To increase the efficiency of this analysis, we provide a diagnostic tool capable of comprehensively displaying multiple orbit determination solutions by leveraging data visualization techniques. This work details the design and visualization decisions made in the creation of the Orbit Determination Dashboard, a tool aimed at giving users an interactive workspace for understanding how changes to input parameters of orbit determination models affect their solutions.

Jah, Moriba↗
DOE OSTI2021

An Analytic Benchmark for Neutron Boltzmann Transport with Downscattering—Part I: Flux and Eigenvalue Solutions

Computing in the energy dimension is one of the greatest challenges confronting present-day deterministic neutron transport solvers. Accurately resolving the neutron flux as neutrons downscatter across resonances in the nuclear cross sections currently requires considerable computing power and suffers from approximation errors. Flux uncertainty resulting from the uncertainty of the resonance structure is the single-largest cause of reactivity uncertainty. Any additional reference solution for the critical neutron downscattering problem with resonance phenomena would be a boon to verification and validation of neutronics codes. This paper establishes a benchmark to verify the accuracy of neutron transport criticality solvers along the energy dimension. For the first time, the analytic solution of the flux amplitude is derived in the particular case of an infinite homogeneous medium with isotropic scattering in the center of mass and an arbitrary number of no-threshold, neutral particle reaction resonances (e.g., radiative capture, fission, and resonance scattering). Furthermore, original analytic expressions are established to quantify the discrepancy between the ψ k (E) and ψ α (E) flux amplitudes, respective solutions of the multiplication factor k, or the exponential time-evolution frequency α eigenproblems. The physical study of these relations led to analysis of their first-order relative difference near the criticality condition α=0. Finally, numerical solutions are provided to a benchmark problem constituted of the first resonance of 239 Pu, the 6.67-eV resonance of 238 U, and a scattering isotope with a flat cross section, allowing for the computational verification of the energy resolution of current neutron transport criticality codes. Through these novel results, this analytic benchmark can serve as a reference to verify the energy resolution and sensitivity analysis of neutron transport criticality calculations.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗