Search NASA⌕ Search

SEARCH · Search NASA

Results for “variable order”

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,099 records · Page 61
All sourcesNASA NTRSDOE OSTI
DOE OSTI2017

An Analysis of Terrestrial and Aquatic Environmental Controls of Riverine Dissolved Organic Carbon in the Conterminous United States

Riverine carbon cycling is an important, but insufficiently investigated component of the global carbon cycle. Analyses of environmental controls on riverine carbon cycling are critical for improved understanding of mechanisms regulating carbon processing and storage along the terrestrial-aquatic continuum. Here, we compile and analyze riverine dissolved organic carbon (DOC) concentration data from 1402 United States Geological Survey (USGS) gauge stations to examine the spatial variability and environmental controls of DOC concentrations in the United States (U.S.) surface waters. DOC concentrations exhibit high spatial variability, with an average of 6.42 ± 6.47 mg C/ L (Mean ± Standard Deviation). In general, high DOC concentrations occur in the Upper Mississippi River basin and the Southeastern U.S., while low concentrations are mainly distributed in the Western U.S. Single-factor analysis indicates that slope of drainage areas, wetlands, forests, percentage of first-order streams, and instream nutrients (such as nitrogen and phosphorus) pronouncedly influence DOC concentrations, but the explanatory power of each bivariate model is lower than 35%. Analyses based on the general multi-linear regression models suggest DOC concentrations are jointly impacted by multiple factors. Soil properties mainly show positive correlations with DOC concentrations; forest and shrub lands have positive correlations with DOC concentrations, but urban area and croplands demonstrate negative impacts; total instream phosphorus and dam density correlate positively with DOC concentrations. Notably, the relative importance of these environmental controls varies substantially across major U.S. water resource regions. In addition, DOC concentrations and environmental controls also show significant variability from small streams to large rivers, which may be caused by changing carbon sources and removal rates by river orders. In sum, our results reveal that general multi-linear regression analysis of twenty one terrestrial and aquatic environmental factors can partially explain (56%) the DOC concentration variation. In conclusion, this study highlights the complexity of the interactions among these environmental factors in determining DOC concentrations, thus calls for processes-based, non-linear methodologies to constrain uncertainties in riverine DOC cycling.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2019

Modeling of the Kinetics of Carbonitride Precipitation Process in High-Strength Low-Alloy Steels Using Cellular Automata Method

Proposition of a cellular automata (CA) model of carbonitride precipitation to simulate image of microstructure in HSLA steels and to predict the mechanical properties is presented in the paper. CA model proved to be very efficient in modeling various phenomena in material science. Numerical modeling played an important role in development of new processing technologies taking advantage of precipitation. By modeling, we mean a mathematical description of the relation between the main process variables and the resulting material properties, based on sound physical principles. In microalloyed steels, the microalloying elements Ti, Nb and V are added in order to control their microstructure and mechanical properties. High chemical affinity of these elements for interstitials (C, N) results in precipitation of binary compound, nitrides and carbides and products of their mutual solubility—carbonitrides. The model accounts for an increase in dislocation density due to plastic deformation and predicts kinetics of precipitation as well as shape of precipitates.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Assessing Global and Local Radiative Feedbacks Based on AGCM Simulations for 1980–2014/2017

In order to avoid contamination of unobservable and uncertain effective radiative forcing (ERF) on the diagnosis of radiative feedbacks based on short-term climate variability, we verify the Kernel-Gregory feedback calculation method using atmospheric model experiments with prescribed ERF. We show that both clear-sky radiative fluxes and all-sky radiative feedbacks have a closure between model simulations and kernel derivations. A near-zero global mean net cloud feedback is found in the simulation with prescribed ERF, which results in a more negative global net climate feedback, -2 W m -2 K -1 . Consistent with AMIP6 ensemble mean results, the lapse rate feedback is the largest contributor among all feedbacks to temperature amplification over the three poles (Arctic, Antarctic and Tibetan Plateau), followed by surface albedo feedback and Planck feedback deviation from its global mean. Except for the higher surface albedo feedback in the Antarctic, other feedbacks are almost same between the Arctic and Antarctic.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Sparse thermal data for cellular automata modeling of grain structure in additive manufacturing

Grain growth in the wake of the melt pool formed during alloy-based additive manufacturing (AM) is complex and multifaceted, depending on parameters governing heat transport, fluid flow, and solidification itself. Cellular automata (CA) models have proven effective in providing computationally efficient and physically sound predictions of grain structure for several AM problems, but their efficiency is tied to the performance of heat transport models. CA models use only a small portion of the problem's temperature data (near the moving melt pool boundary), and much of the CA calculations do not affect the final result due to re-melting of material. Coupling of and communication between heat transport and solidification models, and eliminating operations irrelevant towards final grain structure prediction, will be necessary for using these methods for efficient simulation of large parts. In this work, we introduce a procedure of decoupling the CA from temperature field simulation, using files of relevant temperature data written by the heat transport model. This approach is validated against the standard coupling approach using data obtained through the computational fluid dynamics software OpenFOAM. Negligible differences are seen in grain size, volume, and texture distributions for multilayer simulation of test problems, while the quantity of temperature data for these test problems was reduced by four orders of magnitude (from 100s of GB to 10s of MB) and the code performance sped up by a factor of around 50. Furthermore, variability in microstructure as a function of cell size, substrate, time step, and nucleation parameters is studied, and it is found that cell sizes less than or equal to 1.67 μm and sufficiently small time steps yield Αstatistically equivalent microstructures. Finally, a potential use case for this CA approach—the layer-wise convergence in grain structure starting from extremes in initial grain size—is examined. This approach's ability to simulate expected trends in nucleation and epitaxial grain growth for large regions of microstructure, simulated independently of heat transport models themselves, should prove useful for investigation of various microstructure uncertainties and prediction of part-scale experimental results.

36 MATERIALS SCIENCE↗
DOE OSTI2022

A mixed, unified forward/inverse framework for earthquake problems: fault implementation and coseismic slip estimate

SUMMARY We introduce a new finite-element (FE) based computational framework to solve forward and inverse elastic deformation problems for earthquake faulting via the adjoint method. Based on two advanced computational libraries, FEniCS and hIPPYlib for the forward and inverse problems, respectively, this framework is flexible, transparent and easily extensible. We represent a fault discontinuity through a mixed FE elasticity formulation, which approximates the stress with higher order accuracy and exposes the prescribed slip explicitly in the variational form without using conventional split node and decomposition discrete approaches. This also allows the first order optimality condition, that is the vanishing of the gradient, to be expressed in continuous form, which leads to consistent discretizations of all field variables, including the slip. We show comparisons with the standard, pure displacement formulation and a model containing an in-plane mode II crack, whose slip is prescribed via the split node technique. We demonstrate the potential of this new computational framework by performing a linear coseismic slip inversion through adjoint-based optimization methods, without requiring computation of elastic Green’s functions. Specifically, we consider a penalized least squares formulation, which in a Bayesian setting—under the assumption of Gaussian noise and prior—reflects the negative log of the posterior distribution. The comparison of the inversion results with a standard, linear inverse theory approach based on Okada’s solutions shows analogous results. Preliminary uncertainties are estimated via eigenvalue analysis of the Hessian of the penalized least squares objective function. Our implementation is fully open-source and Jupyter notebooks to reproduce our results are provided. The extension to a fully Bayesian framework for detailed uncertainty quantification and non-linear inversions, including for heterogeneous media earthquake problems, will be analysed in a forthcoming paper.

58 GEOSCIENCES↗
DOE OSTI2023

Efficiency and Output Power Sensitivity to Control Parameter Variations in an LCC-Series Wireless Power Transfer System

In wireless power transfer (WPT) systems, output power (Pout) can be controlled through input parameters such as, input voltage (Vin) , duty cycle (d) , and switching frequency (fsw) . The impact of these control parameters on the overall system efficiency has not been explored in previous studies. This paper focuses on the efficiency (η) and Pout sensitivity to these control parameter variations in an LCC-Series tuned WPT system. There are multiple combinations of Vin, d , and fsw that can achieve a specified target Pout . However, it is important to determine the optimal combination of these control variables that will yield the maximum η for a target Pout . A MATLAB/Simulink model is built to calculate and measure η and Pout . In order to achieve this result, several simulations are run where Vin,d , and fsw are swept across a wide range to cover all possible combinations. MATLAB's Sensitivity Analysis toolbox is used to evaluate how the parameters influence the model's output. Finally, an experimental prototype is implemented. The optimal values acquired using the sensitivity analysis is then plugged into the hardware prototype and compared with the simulation. The results obtained confirm the validity of the proposed method.

Harave, Sudarshan↗
DOE OSTI2014

Predictability of the terrestrial carbon cycle

Abstract Terrestrial ecosystems sequester roughly 30% of anthropogenic carbon emission. However this estimate has not been directly deduced from studies of terrestrial ecosystems themselves, but inferred from atmospheric and oceanic data. This raises a question: to what extent is the terrestrial carbon cycle intrinsically predictable? In this paper, we investigated fundamental properties of the terrestrial carbon cycle, examined its intrinsic predictability, and proposed a suite of future research directions to improve empirical understanding and model predictive ability. Specifically, we isolated endogenous internal processes of the terrestrial carbon cycle from exogenous forcing variables. The internal processes share five fundamental properties (i.e., compartmentalization, carbon input through photosynthesis, partitioning among pools, donor pool‐dominant transfers, and the first‐order decay) among all types of ecosystems on the Earth. The five properties together result in an emergent constraint on predictability of various carbon cycle components in response to five classes of exogenous forcing. Future observational and experimental research should be focused on those less predictive components while modeling research needs to improve model predictive ability for those highly predictive components. We argue that an understanding of predictability should provide guidance on future observational, experimental and modeling research.

Luo, Yiqi↗
DOE OSTI2014

Algorithms for Mathematical Programming with Emphasis on Bi-level Models (Final Scientific/Technical Report)

The research supported by this grant was focused primarily on first-order methods for solving large scale and structured convex optimization problems and convex relaxations of nonconvex problems. These include optimal gradient methods, operator and variable splitting methods, alternating direction augmented Lagrangian methods, and block coordinate descent methods.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2016

UQ Toolkit v. 3.0

The Uncertainty Quantification (UQ) Toolkit is a software library for the characterization and propagation of uncertainties in computational models. This library provides Bayesian inference tools to infer uncertain model parameters, Bayesian compressive sensing methods for discovering sparse representations of high-dimensional input-output response surfaces, methods for constructing Karhunen-Loève representations of stochastic processes, and global sensitivity analysis tools used to compute Sobol indices in order to characterize the importance of uncertain inputs or the interactions between them. The basis for many of these methods relies on representing random variables with Polynomial Chaos expansions (PCEs). This library implements several spectral basis function types (e.g. Hermite basis functions in terms of Gaussian random variables or Legendre basis functions in terms of uniform random variables) that can be used to represent random variables with PCEs. For the propagation of uncertainty, tools are provided to propagate PCEs that describe the input uncertainty through the computational model using either intrusive methods (Galerkin projection of equations onto basis functions) or non-intrusive methods (perform deterministic operation at sampled values of the random values and project the obtained results onto basis functions).

Sargsyan, Khachik↗
DOE OSTI2025

Advanced Method Optimization with Categorical and Constrained Continuous Parameters

Traditional approaches to analytical method optimization (e.g., univariate and “guess-and-check”) can be time-consuming, costly, and often fail to identify true optima within the parameter space. Previous work defined and implemented a generalized technique for method optimization for continuous method parameters, but a knowledge gap remains for the incorporation of categorical variables into these advanced method optimization schemes. This work presents and validates a generalized optimization approach that incorporates both continuous and categorical variables while also utilizing a multivariate, multiobjective optimization scheme with Karush–Kuhn–Tucker conditions to bound the optimization space to solutions within the physical limitations of the parameter space. Method optimization from a case study using GC–MS for the analysis of 11 analytical standards with objectives to minimize peak width and maximize peak height resulted in a 3 orders of magnitude improvement in the average peak height and a 2 orders of magnitude improvement in the average peak width compared to the least optimal (but reasonable) instrumental parameters utilized in this study. This approach to optimization allows for a customizable method optimization in which users can include both continuous and categorical variables to achieve objectives specific to their analytical goals. This approach significantly reduces the labor and cost associated with traditional method development approaches and can be applied in a variety of scientific fields across a range of laboratory techniques (e.g., instrument method development, sample preparation, and extraction techniques).

Amorphous materials↗
DOE OSTI2024

Measurement of jet substructure in boosted $t\overline{t}$ events with the ATLAS detector using 140 fb -1 of 13 TeV $pp$ collisions

Measurements of the substructure of top-quark jets are presented, using 140 fb -1 of 13 TeV pp collision data recorded with the ATLAS detector at the LHC. Top-quark jets reconstructed with the anti-k t algorithm with a radius parameter R = 1.0 are selected in top-quark pair ($t\overline{t}$) events where one top quark decays semileptonically and the other hadronically, or where both top quarks decay hadronically. The top-quark jets are required to have transverse momentum p T > 350 GeV, yielding large samples of data events with jet p T values between 350 and 600 GeV. One- and two-dimensional differential cross sections for eight substructure variables, defined using only the charged components of the jets, are measured in a particle-level phase space by correcting for the smearing and acceptance effects induced by the detector. The differential cross sections are compared with the predictions of several Monte Carlo simulations in which top-quark pair-production quantum chromodynamic matrix-element calculations at next-to-leading-order precision in the strong coupling constant α S are passed to leading-order parton shower and hadronization generators. The Monte Carlo predictions for measures of the broadness, and also the two-body structure, of the top-quark jets are found to be in good agreement with the measurements, while variables sensitive to the three-body structure of the top-quark jets exhibit some tension with the measured distributions.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2024

Reduced‐Order Modeling for Linearized Representations of Microphysical Process Rates

Abstract Representing cloud microphysical processes in large scale atmospheric models is challenging because many processes depend on the details of the droplet size distribution (DSD, the spectrum of droplets with different sizes in a cloud). While full or partial statistical moments of droplet size distributions are the typical variables used in bulk models, prognostic moments are limited in their ability to represent microphysical processes across the range of conditions experienced in the atmosphere. Microphysical parameterizations employing prognostic moments are known to suffer from structural uncertainty in their representations of inherently higher dimensional cloud processes, which limit model fidelity and lead to forecasting errors. Here we investigate how data‐driven reduced‐order modeling can be used to learn predictors for microphysical process rates in bulk microphysics schemes in an unsupervised manner from higher dimensional bin distributions. Using simulations characteristic of marine stratiform clouds, we simultaneously learn lower dimensional representations of droplet size distributions and predict the evolution of the microphysical state of the system. Droplet collision‐coalescence, the main process for generating warm rain, is estimated to have an intrinsic dimension of three. This intrinsic dimension provides a lower limit on the number of degrees of freedom needed to accurately represent collision‐coalescence in models. We demonstrate how deep learning based reduced‐order modeling can be used to discover intrinsic coordinates describing the microphysical state of the system, where process rates such as collision‐coalescence are globally linearized. These implicitly learned representations of the DSD retain more information about the DSD than typical moment‐based representations.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Generating An Advanced Cross-section Library For HTGR Pebble Bed Depletion Calculations Using Reduced-Order Model Generation Techniques

For code development, Advanced Reactor Technologies - Gas Cooled Reactors Program (ART-GCR) rely on a collaboration with the Nuclear Energy Advanced Modeling and Simulation (NEAMS) program, but the cross sections generation and the methodology definition is part of this program area goals. Based on previous studies in FY23, the size of microscopic cross section libraries increases rapidly with the number of tabulations, requiring significant amount of memory and drastically slowing down the Griffin calculations when evaluating cross sections via the multivariate linear interpolation approach. Rising to these challenges, this work investigates constructing Reduced-order Models (ROMs) for the multi-group microscopic cross sections to accelerate the cross section evaluation in Griffin. A database of multigroup cross sections is first collected considering all possible parameters that a designer could change for optimization. Down-selection of the ROM techniques afterward shows Deep Neural Network (DNN) as the best candidate when jointly consider memory efficiency, predictive accuracy, computational cost, scalability, flexibility and ease of implementation of the algorithms in comparison to the multidimensional interpolation. This work develops a specific interface that enables the cross section predictions using pre-trained DNN models into Griffin leveraging the existing ROM capabilities. DNNs have been trained for all isotopes for use in Griffin. Preliminary Griffin testing shows that DNNs exhibit exceptional predictive accuracy and the use of DNNs provides orders of magnitude improvement in memory efficiency compared to conventional interpolation techniques. With such ROM techniques, it holds great promise to further increase the fidelity of the Pebble Bed Reactor (PBR) simulation by increasing the number of tabulations/state variables during cross section evaluation, while maintaining the computational cost affordable in Griffin.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2023

Evaluating Shallow Convection Parameterization Assumptions With a q t – w Quadrant Analysis

Abstract Uncertainties associated with the response of shallow clouds to global warming remain a great challenge for climate projection. Due to the small size of these clouds, parameterizations are required to represent them in both current‐ and next‐generation climate models. We present a quantitative evaluation of several important assumptions used in both mass‐flux and Assumed‐PDF Higher‐Order Closure parameterizations of shallow convection. We use large‐eddy simulations (LESs) of four different shallow convection regimes as benchmarks and apply a q t (total moisture)– w (vertical velocity) quadrant analysis technique to identify the “coherent structures” of moist and dry up/down‐drafts. The statistics of these coherent structures and the environment are then used to evaluate commonly used assumptions in mass‐flux parameterizations related to (a) the convective downdrafts and (b) the subplume variability of convective updrafts. For quantitative evaluation of the PDF closure used in AP‐HOCs, we perform offline calculations with the PDF closure used in the Cloud Layers Unified By Bi‐normals and the Simplified Higher‐Order Closure schemes by supplying the closure with statistical quantities directly calculated from the LESs. Then, the impact of parameterization assumptions embedded in the PDF closure in the same two categories is analyzed and understood through comparisons of the output from the closure with the statistics of the moist and dry up/down‐drafts from the LESs.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2015

Magnitude and Variability of Controllable Charge Capacity Provided by Grid Connected Plug-in Electric Vehicles

As market penetration of plug-in electric vehicles (PEV) increases over time, the number of PEVs charging on the electric grid will also increase. As the number of PEVs increases, their ability to collectively impact the grid increases. The idea of a large body of PEVs connected to the grid presents an intriguing possibility. If utilities can control PEV charging, it is possible that PEVs could act as a distributed resource to provide grid services. The technology required to control charging is available for modern PEVs. However, a system for wide-spread implementation of controllable charging, including robust communication between vehicles and utilities, is not currently present. Therefore, the value of controllable charging must be assessed and weighed against the cost of building and operating such as system. In order to grasp the value of PEV charge control to the utility, the following must be understood: 1. The amount of controllable energy and power capacity available to the utility 2. The variability of the controllable capacity from day to day and as the number of PEVs in the market increases.

33 ADVANCED PROPULSION SYSTEMS↗
DOE OSTI2019

Propagating uncertainties in large-scale hemodynamics models via network uncertainty quantification and reduced-order modeling

Numerical simulations of the cardiovascular system are affected by uncertainties arising from a substantial lack of data related to the boundary conditions and the physical parameters of the mathematical models. Quantifying the impact of this uncertainty on the numerical results along the circulatory network is challenged by the complexity of both the morphology of the domain and the local dynamics. Here, we propose to integrate (i) the Transverse Enriched Pipe Element Methods (TEPEM) as a reduced-order model for effectively computing the 3D local hemodynamics; and (ii) a combination of uncertainty quantification via Polynomial Chaos Expansion and classical relaxation methods – called network uncertainty quantification (NetUQ) – for effectively propagating random variables that encode uncertainties throughout the networks. The findings demonstrate the computational effectiveness of computing the propagation of uncertainties in networks with nontrivial topology, including portions of the cerebral and the coronary systems.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Structural, magneto transport and magnetic properties of Ruddlesden–Popper La 2-2x Sr 1+2x Mn 2 O 7 (0.42≤x≤0.52) layered manganites

The magneto transport of Ruddlesden–Popper, La 2-2x Sr 1+2x Mn 2 O 7 (0.42≤x≤0.52), 2-dimensional bilayered manganites have been investigated in a broad temperature range. The samples have been synthesized using the solid-state reaction method. Rietveld refinement of the X-ray diffraction data indicates the tetragonal structure formation with I4/mmm space group. The resistivity curves of the samples present a general characteristic of metal-insulator (MI) transition at certain temperature (T MI ). Besides, for samples with x = 0.48, 0.50, and 0.52, at a certain temperature (T CO ) well below T MI the charge ordering is also evident. Furthermore, the samples display a shallow upturn in resistivity below T min due to the Kondo like spin scattering effect, weak localization and electron-phonon interaction. The high temperature semiconducting (T>T MI ) region, the resistivity curve follows the 3D Mott’s variable hopping transport mechanism. The overall suppression of resistivity accounts for a substantial magnetoresistance by applying a magnetic field and the characteristic change in T MI , T CO , and T min are discussed. Temperature-dependent magnetization demonstrates the suppression of ferromagnetism and evident of antiferromagnetic nature with an apparent charge ordering with increasing concentration of Sr 2+ .

36 MATERIALS SCIENCE↗
DOE OSTI2019

Renormalized Reduced Order Models with Memory for Long Time Prediction

We examine the challenging problem of constructing reduced models for the long time prediction of systems where there is no timescale separation between the resolved and unresolved variables. In previous work we focused on the case where there was only transfer of activity (e.g. energy, mass) from the resolved to the unresolved variables. Here we investigate the much more difficult case where there is {\it two-way} transfer of activity between the resolved and unresolved variables. Like in the case of activity drain out of the resolved variables, even if one starts with an exact formalism, like the Mori-Zwanzig (MZ) formalism, the constructed reduced models can become unstable. We show how to remedy this situation by using dynamic information from the full system to {\it renormalize} the MZ reduced models. In addition to being stabilized, the renormalized models can be accurate for very long times. We use the Korteweg-de Vries equation to illustrate the approach. The coefficients of the renormalized models exhibit rich structure, including algebraic time dependence and incomplete similarity.

Price, Jacob↗