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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 289 records · Page 16

All sources NASA NTRS DOE OSTI

NASA NTRS1992

Stochastic modeling of turbulent reacting flows

Direct numerical simulations of a single-step irreversible chemical reaction with non-premixed reactants in forced isotropic turbulence at R(sub lambda) = 63, Da = 4.0, and Sc = 0.7 were made using 128 Fourier modes to obtain joint probability density functions (pdfs) and other statistical information to parameterize and test a Fokker-Planck turbulent mixing model. Preliminary results indicate that the modeled gradient stretching term for an inert scalar is independent of the initial conditions of the scalar field. The conditional pdf of scalar gradient magnitudes is found to be a function of the scalar until the reaction is largely completed. Alignment of concentration gradients with local strain rate and other features of the flow were also investigated.

Fox, R. O.↗

DOE OSTI2024

Advanced surrogate model for electron-scale turbulence in tokamak pedestals

We derive an advanced surrogate model for predicting turbulent transport at the edge of tokamaks driven by electron temperature gradient (ETG) modes. Our derivation is based on a recently developed sensitivity-driven sparse grid interpolation approach for uncertainty quantification and sensitivity analysis at scale, which informs the set of parameters that define the surrogate model as a scaling law. Our model reveals that ETG-driven electron heat flux is influenced by the safety factor q, electron beta β e and normalized electron Debye length λ D , in addition to well-established parameters such as the electron temperature and density gradients. To assess the trustworthiness of our model's predictions beyond training, we compute prediction intervals using bootstrapping. The surrogate model's predictive power is tested across a wide range of parameter values, including within-distribution testing parameters (to verify our model) as well as out-of-bounds and out-of-distribution testing (to validate the proposed model). Overall, validation efforts show that our model competes well with, or can even outperform, existing scaling laws in predicting ETG-driven transport.

fusion plasma↗

NASA NTRS2022

A Dynamic Landslide Hazard Monitoring Framework for the Lower Mekong Region

The Lower Mekong region is one of the most landslide-prone areas of the world. Despite the need for dynamic characterization of landslide hazard zones within the region, it is largely understudied for several reasons. Dynamic and integrated understanding of landslide processes requires landslide inventories across the region, which have not been available previously. Computational limitations also hamper regional landslide hazard assessment, including accessing and processing remotely sensed information. Finally, open-source software and modelling packages are required to address regional landslide hazard analysis. Leveraging an open-source data-driven global Landslide Hazard Assessment for Situational Awareness model framework, this study develops a region-specific dynamic landslide hazard system leveraging satellite-based Earth observation data to assess landslide hazards across the lower Mekong region. A set of landslide inventories were prepared from high-resolution optical imagery using advanced image-processing techniques. Several static and dynamic explanatory variables (i.e., rainfall, soil moisture, slope, relief, distance to roads, distance to faults, distance to rivers) were considered during the model development phase. An extreme gradient boosting decision tree model was trained for the monsoon period of 2015–2019 and the model was evaluated with independent inventory information for the 2020 monsoon period. The model performance demonstrated considerable skill using receiver operating characteristic curve statistics, with Area Under the Curve values exceeding 0.95. The model architecture was designed to use near-real-time data, and it can be implemented in a cloud computing environment (i.e., Google Cloud Platform) for the routine assessment of landslide hazards in the Lower Mekong region. This work was developed in collaboration with scientists at the Asian Disaster Preparedness Center as part of the NASA SERVIR Program’s Mekong hub. The goal of this work is to develop a suite of tools and services on accessible open-source platforms that support and enable stakeholder communities to better assess landslide hazard and exposure at local to regional scales for decision making and planning.

Nishan Kumar Biswas↗

DOE OSTI2023

Distributed Power Allocation for 6-GHz Unlicensed Spectrum Sharing via Multi-agent Deep Reinforcement Learning

We consider the problem of power allocation over the 6 GHz Unlicensed National Information Infrastructure (UNII)- 5 spectrum. We propose a novel deep Reinforcement Learning (DRL)-based distributed power allocation scheme which utilizes the multi-agent Deep Deterministic Policy Gradient (MADDPG) algorithm. In particular, we model the base stations (BSs) as DRL agents that simultaneously determine the transmit powers to their scheduled user equipment (UE) in a synchronized manner. The power decision of each BS is based on its own observation of the radio environment, which consists of several local interference measurements and a limited amount of information obtained from other BSs. One advantage of the proposed scheme is that it addresses the single-agent non-stationarity problem of RL in the multi-agent scenario by incorporating the actions and observations of other BSs into each BS’s own critic which helps it to gain a more accurate perception of the overall radio environment. A centralized-training-distributed execution framework is used to train the policies where the critics are trained over the joint actions and observations of all BSs while the actor of each BS only takes the local observation as input in order to produce the transmit power. Simulation shows that the proposed power allocation scheme can achieve better throughput performance than several state-of-the-art approaches.

99 GENERAL AND MISCELLANEOUS↗

DOE OSTI2023

Smoothing Lexis diagrams using kernel functions: A contemporary approach

Lexis diagrams are rectangular arrays of event rates indexed by age and period. Analysis of Lexis diagrams is a cornerstone of cancer surveillance research. Typically, population-based descriptive studies analyze multiple Lexis diagrams defined by sex, tumor characteristics, race/ethnicity, geographic region, etc. Inevitably the amount of information per Lexis diminishes with increasing stratification. Several methods have been proposed to smooth observed Lexis diagrams up front to clarify salient patterns and improve summary estimates of averages, gradients, and trends. In this article, we develop a novel bivariate kernel-based smoother that incorporates two key innovations. First, for any given kernel, we calculate its singular values decomposition, and select an optimal truncation point—the number of leading singular vectors to retain—based on the bias-corrected Akaike information criterion. Second, we model-average over a panel of candidate kernels with diverse shapes and bandwidths. The truncated model averaging approach is fast, automatic, has excellent performance, and provides a variance-covariance matrix that takes model selection into account. We present an in-depth case study (invasive estrogen receptor-negative breast cancer incidence among non-Hispanic white women in the United States) and simulate operating characteristics for 20 representative cancers. The truncated model averaging approach consistently outperforms any fixed kernel. Our results support the routine use of the truncated model averaging approach in descriptive studies of cancer.

60 APPLIED LIFE SCIENCES↗

DOE OSTI2021

Deep Kronecker neural networks: A general framework for neural networks with adaptive activation functions

Here we propose a new type of neural networks, Kronecker neural networks (KNNs), that form a general framework for neural networks with adaptive activation functions. KNNs employ the Kronecker product, which provides an efficient way of constructing a very wide network while keeping the number of parameters low. Our theoretical analysis reveals that under suitable conditions, KNNs induce a faster decay of the loss than that by the feed-forward networks. This is also empirically verified through a set of computational examples. Furthermore, under certain technical assumptions, we establish global convergence of gradient descent for KNNs. As a specific case, we propose the Rowdy activation function that is designed to get rid of any saturation region by injecting sinusoidal fluctuations, which include trainable parameters. The proposed Rowdy activation function can be employed in any neural network architecture like feed-forward neural networks, Recurrent neural networks, Convolutional neural networks etc. The effectiveness of KNNs with Rowdy activation is demonstrated through various computational experiments including function approximation using feed-forward neural networks, solution inference of partial differential equations using the physics-informed neural networks, and standard deep learning benchmark problems using convolutional and fully-connected neural networks.

97 MATHEMATICS AND COMPUTING↗

NASA NTRS1995

Optical diagnostics of solution crystal growth

Solution crystal growth monitoring of LAP/TGS crystals by various optical diagnostics systems, such as conventional and Mach-Zehnder (M-Z) interferometers, optical heterodyne technique, and ellipsometry, is under development. The study of the dynamics of the crystal growth process requires a detailed knowledge of crystal growth rate and the concentration gradient near growing crystals in aqueous solution. Crystal growth rate can be measured using conventional interferometry. Laser beam reflections from the crystal front as well as the back surface interfere with each other, and the fringe shift due to the growing crystal yields information about the growth rate. Our preliminary results indicate a growth rate of 6 A/sec for LAP crystals grown from solution. Single wavelength M-Z interferometry is in use to calculate the concentration gradient near the crystal. Preliminary investigation is in progress using an M-Z interferometer with 2 cm beam diameter to cover the front region of the growing crystal. In the optical heterodyne technique, phase difference between two rf signals (250 KHZ) is measured of which one is a reference signal, and the other growth signal, whose phase changes due to a change in path length as the material grows. From the phase difference the growth rate can also be calculated. Our preliminary results indicate a growth rate of 1.5 A/sec. the seed and solution temperatures were 26.46 C and 27.92 C respectively, and the solution was saturated at 29.0 C. an ellipsometer to measure the growth rate and interface layer is on order from JOBIN YVON, France. All these systems are arranged in such a manner that measurements can be made either sequentially or simultaneously. These techniques will be adapted for flight experiment.

Kim, Yongkee↗

DOE OSTI2024

Investigation of W-SiC compositionally graded films as a divertor material

W-SiC composite material is a promising plasma-facing material candidate alternative to pure W due to the low neutron activation, low impurity radiation, and low tritium diffusivity of SiC while leveraging the high erosion resistance of the W armor. Additionally, W and SiC have high thermomechanical compatibility given their similar thermal expansion rates. The present study addresses the synthesis and performance of compositionally graded W-SiC films fabricated by pulsed-DC magnetron sputtering. Compositional gradients were characterized using transmission electron microscopy (TEM) and energy-dispersive X-ray spectroscopy (EDS), and crystallographic information was obtained using electron diffraction and X-ray diffraction (XRD). Samples were exposed to L-mode deuterium plasma discharges in the DIII-D tokamak using the Divertor Material Evaluation System (DiMES). Post-mortem characterizations were performed using scanning electron microscopy (SEM) and XRD. Electron diffraction and XRD showed that the compositionally graded W-SiC films were composed of polycrystalline W and amorphous SiC with amorphous W+SiC interlayers. No macroscopic delamination or microstructural changes were observed under mild exposure conditions. Here, this study serves as a preliminary examination of W-SiC compositionally graded composites as a potential candidate divertor material in future tokamak devices.

36 MATERIALS SCIENCE↗

DOE OSTI2026

Data-based filtered dissipation rate modelling for multi-modal turbulent combustion: evaluating a priori model generalizability

Manifold-based models offer a computationally efficient alternative to directly transporting the thermochemical state in computational simulations of turbulent reacting flows, projecting the high-dimensional thermochemical state-space onto a low-dimensional manifold. Recent efforts have yielded a manifold-based model applicable to multi-modal combustion, enabling reconstruction of the thermochemical state from solutions to two-dimensional manifold equations in mixture fraction and generalized progress variable that are parameterised by three scalar dissipation rates. In coarse-grained simulations such as Large Eddy Simulation (LES), closure of the multi-modal manifold equations and subfilter variances/covariance requires closure of three filtered scalar dissipation rates. Here, the present work adopts a data-based approach, providing closure for the three filtered scalar dissipation rates via deep neural networks (DNNs). High-fidelity datasets corresponding to an autoigniting n-dodecane jet flame and a bluff body swirl-stabilized confined lifted spray flame of two aviation fuels (Jet-A and C1) with different ignition propensities are leveraged to generate training data that spans a diverse range of thermodynamic conditions and combustion modes, including low- and high-temperature ignition regimes in addition to premixed and nonpremixed behaviour. A final DNN model is trained to enforce inherent physical constraints by learning nonlinear functional transformations of the three filtered scalar dissipation rates. The generalizability of this constrained DNN model is demonstrated a priori via conditional statistics evaluated on the lifted spray flame with C1–a configuration that had not been included in the training data. Excellent DNN agreement with conditional DNS statistics is observed, and integrated gradients are computed to identify the most sensitive input variables. The similarity of the marginal PDFs of the most informative input variables and outputs across configurations are quantified via the Wasserstein metric, demonstrating that data-based models may successfully generalize to unseen parametric conditions so long as the most informative input variables share similar distributions across training and testing datasets.

Data-based modelling↗

DOE OSTI2024

Reducing measurement costs by recycling the Hessian in adaptive variational quantum algorithms

Abstract Adaptive protocols enable the construction of more efficient state preparation circuits in variational quantum algorithms (VQAs) by utilizing data obtained from the quantum processor during the execution of the algorithm. This idea originated with Adaptive Derivative-Assembled Problem-Tailored variational quantum eigensolver (ADAPT-VQE), an algorithm that iteratively grows the state preparation circuit operator by operator, with each new operator accompanied by a new variational parameter, and where all parameters acquired thus far are optimized in each iteration. In ADAPT-VQE and other adaptive VQAs that followed it, it has been shown that initializing parameters to their optimal values from the previous iteration speeds up convergence and avoids shallow local traps in the parameter landscape. However, no other data from the optimization performed at one iteration is carried over to the next. In this work, we propose an improved quasi-Newton optimization protocol specifically tailored to adaptive VQAs. The distinctive feature in our proposal is that approximate second derivatives of the cost function are recycled across iterations in addition to optimal parameter values. We implement a quasi-Newton optimizer where an approximation to the inverse Hessian matrix is continuously built and grown across the iterations of an adaptive VQA. The resulting algorithm has the flavor of a continuous optimization where the dimension of the search space is augmented when the gradient norm falls below a given threshold. We show that this inter-optimization exchange of second-order information leads the approximate Hessian in the state of the optimizer to be consistently closer to the exact Hessian. As a result, our method achieves a superlinear convergence rate even in situations where the typical implementation of a quasi-Newton optimizer converges only linearly. Our protocol decreases the measurement costs in implementing adaptive VQAs on quantum hardware as well as the runtime of their classical simulation.

Ramôa, Mafalda (ORCID:0000000302187801)↗

DOE OSTI2023

A review on experimentally observed mechanical and microstructural characteristics of interfaces in multi-material laser powder bed fusion

Additive manufacturing (AM) is a revolutionary technology. One of the key AM categories, metal powder-based fusion processes, has many advantages compared to conventional methods for fabricating structural materials, such as permitting increased geometric complexity. While single material metal powder AM has advanced significantly in the past decade, multi-material AM is gradually attracting more attention owing to the recent breakthrough in multi-material feedstock delivery and the growing interest of fabricating functionally graded components. Multi-material AM offers an alternative route for applications that require location dependent material properties and high geometrical complexity. The AM community has invented several ways to achieve compositional gradients and discrete boundaries in two and three dimensions using mechanical spreading, nozzle-based, electrophotographic, and hybrid techniques. This article reviews the current state of laser powder bed fusion based multi-material AM of metals with focuses on the characteristics of the material interface as well as the properties and performance of the AM built functionally graded materials. We show the common challenges and issues related to material transitions, such as defects, segregation, phase separation, and the efficacy of some potential solutions including material and process optimizations. Additionally, this study evaluates the applicability and limitations of the existing testing standards and methods for measuring mechanical performance of functionally graded materials. Finally, we discuss mechanical testing development opportunities, which can help multi-material AM move towards higher technological maturity. In general, we find that the link between gradient microstructure and mechanical properties is not well understood or studied and suggest several mechanical tests that may better inform this knowledge gap.

42 ENGINEERING↗

DOE OSTI2023

The Chemodynamics of the Stellar Populations in M31 from APOGEE Integrated-light Spectroscopy

We present an analysis of nearly 1000 near-infrared, integrated-light spectra from APOGEE in the inner ~7 kpc of M31. We utilize full-spectrum fitting with A-LIST simple stellar population spectral templates that represent a population of stars with the same age, [M/H], and [α/M]. With this, we determine the mean kinematics, metallicities, α abundances, and ages of the stellar populations of M31's bar, bulge, and inner disk (~4–7 kpc). We find a nonaxisymmetric velocity field in M31 resulting from the presence of a bar. The bulge of M31 is less metal-rich (mean [M/H] = $0.149$$^{+0.067}_{-0.081}$ dex) than the disk, features minima in metallicity on either side of the bar ([M/H] ~ -0.2), and is enhanced in α abundance (mean [α/M] = $0.281$$^{+0.035}_{-0.038}$). The disk of M31 within ~7 kpc is enhanced in both metallicity ([M/H] = $-0.023$$^{+0.050}_{-0.052}$) and α abundance ([α/M] = $0.274$$^{+0.020}_{-0.025}$). Both of these structural components are uniformly old at ≃12 Gyr. We find the mean metallicity increases with distance from the center of M31, with the steepest gradient along the disk major axis (0.043 ± 0.021 dex kpc -1 ). This gradient is the result of changing light contributions from the bulge and disk. The chemodynamics of stellar populations encodes information about a galaxy's chemical enrichment, star formation history, and merger history, allowing us to discuss new constraints on M31's formation. Our results provide a stepping stone between our understanding of the Milky Way and other external galaxies.

79 ASTRONOMY AND ASTROPHYSICS↗

NASA NTRS1972

A proposed lunar orbiting gravity gradiometer experiment.

Analysis of the gravity gradiometer developed by Forward and Bell (1970) suggest that an accuracy, in the range 0.1 to 0.5 EU can be expected in a lunar orbiter application. This accuracy will allow gradient anomalies associated with mascons to be mapped with 1% accuracy and should reveal a great deal of new information about the lunar gravity field. The proposed experiment calls for putting such a gradiometer into a closely circular polar orbit at an average height of about 30 km above the lunar surface. This orbit allows the entire lunar surface to be covered in fourteen days, the gradiometer to be checked twice per revolution and results in successive passes above the lunar surface being spaced at about the resolution limit of about 30 km set both by the satellite altitude and instrumental integration time.

Debra, D. B.↗

NASA NTRS2000

Sensing Clouds via Spacecraft Radio Occultation Observations

The presence of clouds dramatically alters the opacity and radiative transfer within Earth's atmosphere at both short and long wavelengths. Knowledge of cloud top and base is needed to estimate the Outgoing Long wave Radiation (OLR) to space and the net radiation at the surface from a given atmospheric profile. Satellite observations are required to provide the global perspective needed for climate. Cloud top pressure can be determined to some accuracy from spaceborne radiance measurements when cloud opacities are sufficiently large and cloud top temperatures are sufficiently low to readily them from the surface. Cloud base is particularly difficult to determine from space. The relatively small sensitivity of long wavelengths must be used to penetrate the clouds while at the same time providing sufficient sensitivity to detect the cloud base. We are studying the indirect sensitivity of radio occultation observations to clouds through their impact on the refractivity structure. The tradewind inversion is an excellent example where a very sharp refractivity structure coincides with the top of the tradewind cumulus and stratus clouds. In general, any cloud with sufficient IR opacity will have large transmissivity gradient at cloud top (base) which will result in large cooling (heating) which will create a thermal inversion at cloud top (base). Both liquid and ice clouds can reach this critical opacity. The thermal inversion and sharp change in specific humidity will cause a sharp change in the refractivity gradient which can be identified in the radio occultation results, particularly the signal amplitude. The occultations yield very accurate information on the height of such features. The heating at cloud base drives convection causing the air to follow a moist adiabat within the Cloud which further helps constrain the interpretation of the observations. The upward expansion of such opaque clouds due to convection may be an important mechanism supplying moisture into the upper troposphere. We will present simulations using high resolution radiosondes from field campaigns representing the expected signatures of such features. We will also present initial results of comparisons between GPS observations and GOES-derived cloud tops to assess the utility of this concept.

Kursinski, E. R.↗

DOE OSTI2021

Machine Learning Predicts the Timing and Shear Stress Evolution of Lab Earthquakes Using Active Seismic Monitoring of Fault Zone Processes

Abstract Machine learning (ML) techniques have become increasingly important in seismology and earthquake science. Lab‐based studies have used acoustic emission data to predict time‐to‐failure and stress state, and in a few cases, the same approach has been used for field data. However, the underlying physical mechanisms that allow lab earthquake prediction and seismic forecasting remain poorly resolved. Here, we address this knowledge gap by coupling active‐source seismic data, which probe asperity‐scale processes, with ML methods. We show that elastic waves passing through the lab fault zone contain information that can predict the full spectrum of labquakes from slow slip instabilities to highly aperiodic events. The ML methods utilize systematic changes in P‐wave amplitude and velocity to accurately predict the timing and shear stress during labquakes. The ML predictions improve in accuracy closer to fault failure, demonstrating that the predictive power of the ultrasonic signals improves as the fault approaches failure. Our results demonstrate that the relationship between the ultrasonic parameters and fault slip rate, and in turn, the systematically evolving real area of contact and asperity stiffness allow the gradient boosting algorithm to “learn” about the state of the fault and its proximity to failure. Broadly, our results demonstrate the utility of physics‐informed ML in forecasting the imminence of fault slip at the laboratory scale, which may have important implications for earthquake mechanics in nature.

58 GEOSCIENCES↗

DOE OSTI2021

Manifold Sampling for Optimizing Nonsmooth Nonconvex Compositions

Here we propose a manifold sampling algorithm for minimizing a nonsmooth composition $f= h\circ F$, where we assume $h$ is nonsmooth and may be inexpensively computed in closed form and $F$ is smooth but its Jacobian may not be available. We additionally assume that the composition $h\circ F$ defines a continuous selection. Manifold sampling algorithms can be classified as model-based derivative-free methods, in that models of $F$ are combined with particularly sampled information about $h$ to yield local models for use within a trust-region framework. We demonstrate that cluster points of the sequence of iterates generated by the manifold sampling algorithm are Clarke stationary. We consider the tractability of three particular subproblems generated by the manifold sampling algorithm and the extent to which inexact solutions to these subproblems may be tolerated. Numerical results demonstrate that manifold sampling as a derivative-free algorithm is competitive with state-of-the-art algorithms for nonsmooth optimization that utilize first-order information about $f$.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2020

Deformation accommodating periodic computational domain for a uniform velocity gradient

Many multiscale methods for granular materials use periodic computational domains to consider particle scale interactions and then calculate the stress to drive the continuum scale calculations. For problems involving large material deformations, the computation domain often needs to be reinitialized because of the distortion, causing the loss of the history information of the system. This work introduces an algorithm to accommodate a large deformation of the material while maintaining the computational domain cuboid to avoid the domain reinitialization during the computation. The algorithm uses a rotating frame of reference, in which the velocity gradient can be represented by an upper triangular matrix. The deformation caused by the upper triangular matrix is treated by the image system implied in the periodicity to maintain the computational domain cuboid. The effect from the rotation of the reference frame is considered using the inertial forces. Finally, simulations of simple and pure shear motions are carried out to illustrate the algorithm.

42 ENGINEERING↗

DOE OSTI2023

Bayesian sequential optimal experimental design for nonlinear models using policy gradient reinforcement learning

We present a mathematical framework and computational methods for optimally designing a finite sequence of experiments. This sequential optimal experimental design (sOED) problem is formulated as a finite-horizon partially observable Markov decision process (POMDP) under a Bayesian setting and with information-theoretic utilities. The formulation is general and may accommodate continuous random variables, non-Gaussian posteriors, and nonlinear forward models. The sOED design policy incorporates elements of feedback and lookahead simultaneously, and we show it to generalize the commonly-used batch and greedy design strategies. We solve for the sOED policy using the policy gradient (PG) method from reinforcement learning, and provide a derivation for the PG expression in the sOED context. Adopting an actor-critic approach, the policy and value functions are parameterized using deep neural networks and improved via PG estimates produced from simulated episodes of designs and observations. The new PG-sOED algorithm is first validated on a linear-Gaussian benchmark, and then compared against other design baselines on a sensor movement problem for contaminant source inversion in a convection-diffusion field. As a result, we provide explanation for the policy behaviors using knowledge of the underlying physical process.

97 MATHEMATICS AND COMPUTING↗