Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Gradient information”

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
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Distance‐based reconstruction of protein quaternary structures from inter‐chain contacts

Abstract Predicting the quaternary structure of protein complex is an important problem. Inter‐chain residue‐residue contact prediction can provide useful information to guide the ab initio reconstruction of quaternary structures. However, few methods have been developed to build quaternary structures from predicted inter‐chain contacts. Here, we develop the first method based on gradient descent optimization (GD) to build quaternary structures of protein dimers utilizing inter‐chain contacts as distance restraints. We evaluate GD on several datasets of homodimers and heterodimers using true/predicted contacts and monomer structures as input. GD consistently performs better than both simulated annealing and Markov Chain Monte Carlo simulation. Starting from an arbitrarily quaternary structure randomly initialized from the tertiary structures of protein chains and using true inter‐chain contacts as input, GD can reconstruct high‐quality structural models for homodimers and heterodimers with average TM‐score ranging from 0.92 to 0.99 and average interface root mean square distance from 0.72 Å to 1.64 Å. On a dataset of 115 homodimers, using predicted inter‐chain contacts as restraints, the average TM‐score of the structural models built by GD is 0.76. For 46% of the homodimers, high‐quality structural models with TM‐score ≥ 0.9 are reconstructed from predicted contacts. There is a strong correlation between the quality of the reconstructed models and the precision and recall of predicted contacts. Only a moderate precision or recall of inter‐chain contact prediction is needed to build good structural models for most homodimers. Moreover, GD improves the quality of quaternary structures predicted by AlphaFold2 on a Critical Assessment of Techniques for Protein Structure Prediction–Critical Assessments of Predictions of Interactions dataset.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2021

Learning viscoelasticity models from indirect data using deep neural networks

In this study, we propose a novel approach to model viscoelasticity materials, where rate-dependent and non-linear constitutive relationships are approximated with deep neural networks. We assume that inputs and outputs of the neural networks are not directly observable, and therefore common training techniques with input–output pairs for the neural networks are inapplicable. To that end, we develop a novel computational approach to both calibrate parametric and learn neural-network-based constitutive relations of viscoelasticity materials from indirect displacement data in the context of multiple-physics systems. We show that limited displacement data holds sufficient information to quantify the viscoelasticity behavior. We formulate the inverse computation – modeling viscoelasticity properties from observed displacement data – as a PDE-constrained optimization problem and minimize the error functional using a gradient-based optimization method. The gradients are computed by a combination of automatic differentiation and implicit function differentiation rules. The effectiveness of our method is demonstrated through numerous benchmark problems in geomechanics and porous media transport.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Engine Combustion System Optimization Using Computational Fluid Dynamics and Machine Learning: A Methodological Approach

Gasoline compression ignition (GCI) engines are considered an attractive alternative to traditional spark-ignition and diesel engines. Here, a Machine Learning-Grid Gradient Ascent (ML-GGA) approach was developed to optimize the performance of internal combustion engines. ML offers a pathway to transform complex physical processes that occur in a combustion engine into compact informational processes. The developed ML-GGA model was compared with a recently developed Machine Learning-Genetic Algorithm (ML-GA). Detailed investigations of optimization solver parameters and variable limit extension were performed in the present ML-GGA model to improve the accuracy and robustness of the optimization process. Detailed descriptions of the different procedures, optimization tools, and criteria that must be followed for a successful output are provided here. The developed ML-GGA approach was used to optimize the operating conditions (case 1) and the piston bowl design (case 2) of a heavy-duty diesel engine running on a gasoline fuel with a research octane number (RON) of 80. The ML-GGA approach yielded >2% improvements in the merit function, compared with the optimum obtained from a thorough computational fluid dynamics (CFD) guided system optimization. The predictions from the ML-GGA approach were validated with engine CFD simulations. This study demonstrates the potential of ML-GGA to significantly reduce the time needed for optimization problems, without loss in accuracy compared with traditional approaches.

33 ADVANCED PROPULSION SYSTEMS↗
DOE OSTI2022

Microbial extracellular enzyme activity with simulated climate change

It is critical to understand the consequences of environmental change for the microbial regulation of carbon and nutrient cycling. Specifically, understanding microbial community traits, such as extracellular enzyme activity, can help inform nutrient cycling models and address knowledge gaps. We analyzed data on extracellular enzyme activities and litter decomposition from an 18-month experiment in which microbial communities were reciprocally transplanted along a climate gradient in Southern California. Communities were from desert, scrubland, grassland, pine–oak, and subalpine ecosystems. We aimed to test how enzyme activities responded to climate change following transplantation and how those responses related to decomposition rates. We hypothesized that microbial communities would specialize on their native climate conditions, resulting in higher enzyme activities when transplanted back into their native climate. We investigated the relationship between extracellular enzyme Vmax values, substrate mass loss, and microbial biomass as well as variation in these variables across the climate gradient. We found little evidence for climate specialization, and there was rarely a reduction in enzyme functioning after microbial communities were transplanted into new climate conditions. Moreover, observed differences in decomposition were not related to changes in extracellular enzyme potential, although there were significant differences in enzyme activities and decomposition rates across sites. These results suggest that direct, physiological impacts of climate are likely to be important for enzyme-mediated decomposition, but climate specialization will not constrain the microbial response to climate change in our system.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS1973

An iterative approach to the feature selection problem

The B-average divergence for m-distinct classes, resulting from the linear transformation y = Bx, is proposed as a feature selection criterion, where B is a k by n matrix of rank k not greater than n. It is shown that if the B-average divergence resulting from B is large enough, then the probability of misclassification, considered as a function f the class of all k by n matrices, is essentially minimized by B. A computer program, utilizing a gradient procedure, is developed to numerically maximize the B-average divergence and results are presented for the Cl flight line. For this example, corresponding to 9-distinct classes, most of the discriminatory information is found to lie in a 3-dimensional subspace, defined by an appropriately chosen 3 by 12 matrix B.

Decell, H. P., Jr.↗
NASA NTRS1983

Studies of convection in a solidifying binary mixture at reduced gravity

A great deal of interest was generated recently in the possibility of producing new materials in the reduced gravity environment provided during the forthcoming missions of Spacelab. The range of possibilities extend from producing large crystals of uniform properties to manufacturing materials with unique properties. Most of these processes involve the solidification of materials from the liquid state. Convective motions within the liquid during solidification can influence the local material composite and the shape of the solid-liquid interface which may result in solids with non-uniform properties and crystal defects. The microgravity environment of Spacelab is being viewed as one in which the buoyancy forces are eliminated so that convection driven by thermal gradients does occur, resulting in an improved solidification process. However, convection may occur for other reasons and whether convection is negligible or not during solidification constitutes processing in low-gravity environment. Little information exists presently on convection during solidification under such circumstances. A continuation of an analytical investigation into the nature of convective motion in a binary liquid layer due to surface tension forces during its solidification is reported. The onset of convection will be determined through a stability analysis which is described.

Antar, B. N.↗
NASA NTRS1991

Comparison of four stable numerical methods for Abel's integral equation

The 3-D image reconstruction from cone-beam projections in computerized tomography leads naturally, in the case of radial symmetry, to the study of Abel-type integral equations. If the experimental information is obtained from measured data, on a discrete set of points, special methods are needed in order to restore continuity with respect to the data. A new combined Regularized-Adjoint-Conjugate Gradient algorithm, together with two different implementations of the Mollification Method (one based on a data filtering technique and the other on the mollification of the kernal function) and a regularization by truncation method (initially proposed for 2-D ray sample schemes and more recently extended to 3-D cone-beam image reconstruction) are extensively tested and compared for accuracy and numerical stability as functions of the level of noise in the data.

Murio, Diego A.↗
NASA NTRS1996

More About Brazing Or Welding NiAl Without Filler

Two reports present additional information about two processes for joining, brazing, or welding workpieces made of nickel aluminide alloys, without use of filler metal. Joining processes involve uniform heating in vacuum-controlled furnace. Eliminates internal thermal gradients in workpieces joined and greatly reduces tendency toward cracking.

Moore, Thomas J.↗
NASA NTRS2001

Constraining the Mean Crustal Thickness on Mercury

The topography of Mercury is poorly known, with only limited radar and stereo coverage available. However, radar profiles reveal topographic contrasts of several kilometers over wavelengths of approximately 1000 km. The bulk of Mercury's geologic activity took place within the first 1 Ga of the planet's history), and it is therefore likely that these topographic features derive from this period. On Earth, long wavelength topographic features are supported either convectively, or through some combination of isostasy and flexure. Photographic images show no evidence for plume-like features, nor for plate tectonics; I therefore assume that neither convective support nor Pratt isostasy are operating. The composition and structure of the crust of Mercury are almost unknown. The reflectance spectrum of the surface of Mercury is similar to that of the lunar highlands, which are predominantly plagioclase. Anderson et al. used the observed center-of-mass center-of-figure offset together with an assumption of Airy isostasy to infer a crustal thickness of 100-300 km. Based on tidal despinning arguments, the early elastic thickness (T(sub e)) of the (unfractured) lithosphere was approximately equal to or less than 100 km. Thrust faults with lengths of up to 500 km and ages of about 4 Ga B.P. are known to exist on Mercury. Assuming a semicircular slip distribution and a typical thrust fault angle of 10 degrees, the likely vertical depth to the base of these faults is about 45 km. More sophisticated modelling gives similar or slightly smaller answers. The depth to the base of faulting and the elastic layer are usually similar on Earth, and both are thought to be thermally controlled. Assuming that the characteristic temperature is about 750 K, the observed fault depth implies that the heat flux at 4 Ga B.P. is unlikely to be less than 20 mW m(exp -2) for a linear temperature gradient. For an elastic thickness of 45 km, topography at 1000 km wavelength is likely to be about 60% compensated. There are thus likely to be considerable lateral variations in crustal thickness. Additional information is contained in the original extended abstract.

Nimmo, F.↗
NASA NTRS2024

Evaluation of Vertical Patterns in Chlorophyll-A Derived From A Data Assimilating Model of Satellite-Based Ocean Color

Satellite-based sensors of ocean color have become the primary tool to infer changes in surface chlorophyll, while BGC-Argo floats are now filling the information gap at depth. Here we use BGC-Argo data to assess depth-resolved information on chlorophyll-a derived from an ocean biogeochemical model constrained by the assimilation of surface ocean color remote sensing. The data-assimilating model replicates well the general seasonality and meridional gradients in surface and depth-resolved chlorophyll-a inferred from the float array in the Southern Ocean. On average, the model tends to overestimate float-based chlorophyll, particularly at times and locations of high productivity such as the beginning of the spring bloom, subtropical deep chlorophyll maxima, and non-iron limited regions of the Southern Ocean. The highest model RMSE in the upper 50 m with respect to the float array is of 0.6 mg Chl m −3 , which should allow the detection of seasonal changes in float-based biomass (varying between 0.01 and >1 mg Chl m −3 ) but might hinder the identification of subtle changes in chlorophyll at narrow local scales. Both model and float profiling data show good agreement with in situ data from station ALOHA, with model estimates showing a slight accuracy edge in inferring depth-resolved observations. Uncertainties in float bio-optical estimates impede their use as a reliable benchmark for validation, but the general qualitative agreement between model and float data provides confidence in the ability of model to replicate biogeochemical features below the surface, where data is not directly constrained by the assimilation of satellite ocean color.

Lionel A Quintero↗
DOE OSTI2025

Integrating Characteristic Arctic Vegetation in a Land Surface Model Improves Representation of Carbon Dynamics Across a Tundra Landscape: Modeling Archive

This modeling archive is in support of the Next-Generation Ecosystem Experiments in the Arctic (NGEE Arctic) publication "Integrating Characteristic Arctic Vegetation in a Land Surface Model Improves Representation of Carbon Dynamics Across a Tundra Landscape", by Murphy et al. (2025). This archive contains model input files and outputs from landscape-scale simulations conducted using ELM, the land model component of the Department of Energy’s Energy Exascale Earth System Model (E3SM), at the Council NGEE Arctic field site (Council Road mile marker 71) on Alaska’s Seward Peninsula. Input data and model output from two sets of ELM simulations are provided. The first set of simulations were conducted with the two default ELM Arctic plant functional types (PFTs; broadleaf deciduous boreal shrub and a C3 grass) and the second set of simulations were conducted with a set of nine Arctic-specific PFTs including nonvascular mosses and lichens, graminoids, forbs, evergreen dwarf shrubs, three height classes of deciduous shrubs (dwarf, low, and low to tall), and deciduous alder shrubs (Sulman et al., 2021). Parameter names and major parameter changes in the Arctic-specific PFT configuration are described in Sulman et al. (2021) and archived in the Sulman et al. (2021) dataset (see below). Simulations were spatially explicit, covering an approximately 6.4X3.3 km domain at the Council site with a spatial resolution of 100 m for a total of 2,112 simulated grid cells under each ELM PFT configuration. The modeling archive contains meteorological forcing (seven *.nc files and one *.txt file), a domain definition file (one *.nc files), land surface configuration files (two *.nc files), parameter files (two *.nc files), annual ELM output files spanning 1980-2014 (68 *.nc files), and a User’s Guide (*pdf file). Additional information on the provided files is in the “Modeling Archive Contents” section of the User’s Guide. Model outputs are aggregated to the column scale (i.e. PFT-specific outputs are not provided here).

Murphy, Bailey [ORNL] (ORCID:0000000203995221)↗
DOE OSTI2025

Benchmarking optimization methods for materials research: Gradient descent and Bayesian optimization for lithium-ion battery aging diagnostics

Accurate and efficient parameter estimation is essential for battery diagnostics and aging analysis. Here, in this study, we compare two optimization-based approaches—gradient descent and Bayesian optimization—for extracting parameters from differential voltage analysis in lithium-ion batteries. While these techniques are widely used, their relative strengths and limitations for this application are not well understood. The study evaluates the trade-offs between these methods in terms of result quality, computational cost, and reliability within this specific application. The diagnostic results from our battery data suggest adopting gradient descent as an initial method for rapid and efficient analysis, while employing more stable optimization techniques, such as Bayesian optimization, as a verification step to mitigate potential instability. Comparing the two methods provides information on algorithmic choice, while inspiring further discussions on selecting appropriate techniques for specific research tasks.

Zhao, Ziqing [Boston Univ., MA (United States)] (O↗
DOE OSTI2021

Spatially Resolved Potential and Li-Ion Distributions Reveal Performance-Limiting Regions in Solid-State Batteries

The performance of solid-state electrochemical systems is intimately tied to the potential and lithium distributions across electrolyte–electrode junctions that give rise to interface impedance. Here, we combine two operando methods, Kelvin probe force microscopy (KPFM) and neutron depth profiling (NDP), to identify the rate-limiting interface in operating Si-LiPON-LiCoO 2 solid-state batteries by mapping the contact potential difference (CPD) and the corresponding Li distributions. The contributions from ions, electrons, and interfaces are deconvolved by correlating the CPD profiles with Li-concentration profiles and by comparisons with first-principles-informed modeling. Furthermore, we find that the largest potential drop and variation in the Li concentration occur at the anode–electrolyte interface, with a smaller drop at the cathode–electrolyte interface and a shallow gradient within the bulk electrolyte. Correlating these results with electrochemical impedance spectroscopy following battery cycling at low and high rates confirms a long-standing conjecture linking large potential drops with a rate-limiting interfacial process.

25 ENERGY STORAGE↗
DOE OSTI2023

BEYONDPLANCK II. CMB mapmaking through Gibbs sampling

We present a Gibbs sampling solution to the mapmaking problem for cosmic microwave background (CMB) measurements that builds on existing destriping methodology. Gibbs sampling breaks the computationally heavy destriping problem into two separate steps: noise filtering and map binning. Considered as two separate steps, both are computationally much cheaper than solving the combined problem. This provides a huge performance benefit as compared to traditional methods and it allows us, for the first time, to bring the destriping baseline length to a single sample. Here, we applied the Gibbs procedure to simulated Planck 30 GHz data. We find that gaps in the time-ordered data are handled efficiently by filling them in with simulated noise as part of the Gibbs process. The Gibbs procedure yields a chain of map samples, from which we are able to compute the posterior mean as a best-estimate map. The variation in the chain provides information on the correlated residual noise, without the need to construct a full noise covariance matrix. However, if only a single maximum-likelihood frequency map estimate is required, we find that traditional conjugate gradient solvers converge much faster than a Gibbs sampler in terms of the total number of iterations. The conceptual advantages of the Gibbs sampling approach lies in statistically well-defined error propagation and systematic error correction. This methodology thus forms the conceptual basis for the mapmaking algorithm employed in the BEYONDPLANCK framework, which implements the first end-to-end Bayesian analysis pipeline for CMB observations.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Selecting durable building envelope systems with machine learning assisted hygrothermal simulations database

Hygrothermal simulations provide insight into the energy performance and moisture durability of building envelope components under dynamic conditions. The inputs required for hygrothermal simulations are extensive, and carrying out simulations and analyses requires expert knowledge. An expert system, the Building Science Advisor (BSA), has been developed to predict the performance and select the energy-efficient and durable building envelope systems for different climates. The BSA consists of decision rules based on expert opinions and thousands of parametric simulation results for selected wall systems. The number of potential wall systems results in millions, too many to simulate all of them. We present how machine learning can help predict durability data, such as mold growth, while minimizing the number of simulations needed to run. The simulation results are used for training and validation of machine learning tools for predicting wall durability. We tested Artificial Neural Network (ANN) and Gradient Boosted Decision Trees (GBDT) for their applicability and model accuracy. Models developed with both methods showed adequate prediction performance (root mean square error of 0.195 and 0.209, respectively). Finally, we introduce how the information supports guidance for envelope design via an easy-to-use web-based tool that does not require the end-user to run hygrothermal simulations.

Salonvaara, Mikael↗
NASA NTRS1980

A description of the NSSL cases used for a simulated VAS retrieval study

A documentation of eight National Severe Storm Laboratory severe storm cases, which serve as a basis for a simulated VISSR Atmospheric Sounder retrieval study, is presented in this paper. Six of the selected cases provide a control data set to complete the statistical information needed for retrieval techniques based upon the use of regression matrices. The other two cases are to be used in the actual retrieval experiments. The selection was based upon the presence of moisture gradients in the analysis region, the availability of satellite images at the selected time periods, and the extent of cloud cover within the observing network.

Mostek, A.↗
NASA NTRS2002

Low-Cost Evaluation of EO-1 Hyperion and ALI for Detection and Biophysical Characterization of Forest Logging in Amazonia (NCC5-481)

Major uncertainties exist regarding the rate and intensity of logging in tropical forests worldwide: these uncertainties severely limit economic, ecological, and biogeochemical analyses of these regions. Recent sawmill surveys in the Amazon region of Brazil show that the area logged is nearly equal to total area deforested annually, but conversion of survey data to forest area, forest structural damage, and biomass estimates requires multiple assumptions about logging practices. Remote sensing could provide an independent means to monitor logging activity and to estimate the biophysical consequences of this land use. Previous studies have demonstrated that the detection of logging in Amazon forests is difficult and no studies have developed either the quantitative physical basis or remote sensing approaches needed to estimate the effects of various logging regimes on forest structure. A major reason for these limitations has been a lack of sufficient, well-calibrated optical satellite data, which in turn, has impeded the development and use of physically-based, quantitative approaches for detection and structural characterization of forest logging regimes. We propose to use data from the EO-1 Hyperion imaging spectrometer to greatly increase our ability to estimate the presence and structural attributes of selective logging in the Amazon Basin. Our approach is based on four "biogeophysical indicators" not yet derived simultaneously from any satellite sensor: 1) green canopy leaf area index; 2) degree of shadowing; 3) presence of exposed soil and; 4) non-photosynthetic vegetation material. Airborne, field and modeling studies have shown that the optical reflectance continuum (400-2500 nm) contains sufficient information to derive estimates of each of these indicators. Our ongoing studies in the eastern Amazon basin also suggest that these four indicators are sensitive to logging intensity. Satellite-based estimates of these indicators should provide a means to quantify both the presence and degree of structural disturbance caused by various logging regimes. Our quantitative assessment of Hyperion hyperspectral and ALI multi-spectral data for the detection and structural characterization of selective logging in Amazonia will benefit from data collected through an ongoing project run by the Tropical Forest Foundation, within which we have developed a study of the canopy and landscape biophysics of conventional and reduced-impact logging. We will add to our base of forest structural information in concert with an EO-1 overpass. Using a photon transport model inversion technique that accounts for non-linear mixing of the four biogeophysical indicators, we will estimate these parameters across a gradient of selective logging intensity provided by conventional and reduced impact logging sites. We will also compare our physical ly-based approach to both conventional (e.g., NDVI) and novel (e.g., SWIR-channel) vegetation indices as well as to linear mixture modeling methods. We will cross-compare these approaches using Hyperion and ALI imagers to determine the strengths and limitations of these two sensors for applications of forest biophysics. This effort will yield the first physical ly-based, quantitative analysis of the detection and intensity of selective logging in Amazonia, comparing hyperspectral and improved multi-spectral approaches as well as inverse modeling, linear mixture modeling, and vegetation index techniques.

Asner, Gregory P.↗
NASA NTRS1992

Modelling TOVS radiances of synoptic systems

Two operational analysis models, produced by ECMWF and NMC, are compared against observed satellite radiances by converting the analyzed variables into synthetic radiances using a radiative transfer model (RTM). Observed TOVS satellite radiances are used as a ground truth, since they are the only source of continuous, synoptic-scale resolution observations available over the data-sparse tropical Pacific Ocean. It is assumed that if the analysis models correctly simulate atmospheric conditions, the observed and modelled radiance fields will have similar patterns and magnitudes. In each of seven TOVS channels examined, the mean temperature difference between the models was small, but the differences between the models and the observations were large. The mean temperature difference within each channel was larger than the standard deviation in most cases. The moisture channels of the RTM were sensitive to changes in moisture content, but none of the channels were significantly affected by changes in temperature. Single-day modelled radiance analyses contain synoptic information similar to the satellite observations. ECMWF is often too noisy, producing features not verified by the satellite; NMC is often too smooth. NMC more closely resembles the TOVS contours. Modelled water vapor channel analyses resemble observed monthly and composited observations well, but do not accurately reproduce the highly-variable single-day fields. ECMWF generally has more accurate gradients. Neither model simulates the tropics well. The water vapor channels, especially channel eleven (weighted at 700 mb), most closely resembled the tropical plume composite described by McGuirk and Ulsh. The only plume evidence found in the thermal or microwave channels, in either the observations or the models, was a composite trough situated in the correct location, extending from the mid-latitudes into the subtropics. This feature was not readily discernable until the definition and mature stages. Additional information besides the modelled temperatures would be required to produce an accurate forecast for the tropics. Model patterns more closely resembled the observations in the quiescent atmosphere, though the magnitudes were not any closer. ECMWF simulates the TOVS quiescent state better than NMC. The moisture channels continue to present the most realistic patterns. The models contained more detail in the tropics than was verified by the satellite observations, apparently trying to force an active ITCZ. Lag correlations showed TOVS to be more persistent than the models. TOVS also showed a smoother lag trend, but all products appeared to have adequate ability to react to changes in the atmosphere over time.

Coe, Thomas Eddy↗