Search NASA⌕ Search

SEARCH · Search NASA

Results for “gradient”

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 361 records · Page 20
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Functional but not taxonomic diversity increases productivity of Populus in the southeastern United States

Plant interactions like competition and facilitation impact ecosystem function and resilience. Improving our understanding of the relationships between these interactions and community productivity has important implications for managers of production systems in forestry and agriculture as well as conservation science. Populus spp. are an excellent model system for exploring how inter‐ and intraspecific interactions impact ecosystem functions, such as productivity, in forest plantations. In this study, we compared aboveground productivity of six Populus clones from three different taxa grown in monoclonal and mixed‐clonal plots. The different mixture treatments were intended to experimentally test aboveground biomass response to contrasting levels of taxonomic diversity and functional diversity based on nitrogen use characteristics of Populus clones. We hypothesize that functional diversity would be more important than taxonomic diversity in increasing aboveground productivity of mixed‐clonal plantings compared to monocultures. In addition, a subset of treatments was carried out on additional sites representing a productivity gradient in order to determine if the relationship between biodiversity and productivity in these systems diminished at more productive sites as suggested by the stress‐gradient hypothesis. We found that functionally diverse mixtures of clones had greater yield of aboveground biomass than the average of their constituent monocultures, while more taxonomically diverse mixes of clones did not differ from the average of their constituent monocultures. However, when reestablished on sites with extremely high or low productivity, the best performing clone mixture also did not differ from the average of its constituent monocultures. Our results suggest that intimate clone mixtures of Populus have the potential to significantly increase productivity, but results vary by mixture and by site. To capitalize on positive biodiversity effects on yield in production systems, targeted mixtures based on divergent functional traits linked to different use and acquisition strategies for site‐specific limiting resources are most likely to be successful.

BEF↗
DOE OSTI2024

Adaptive immersed isogeometric level-set topology optimization

Here, this paper presents for the first time an adaptive immersed approach for level-set topology optimization using higher-order truncated hierarchical B-spline discretizations for design and state variable fields. Boundaries and interfaces are represented implicitly by the iso-contour of one or multiple level-set functions. An immersed finite element method, the eXtended IsoGeometric Analysis, is used to predict the physical response. The proposed optimization framework affords different adaptively refined higher-order B-spline discretizations for individual design and state variable fields. The increased continuity of higher-order B-spline discretizations together with local refinement enables direct control over the accuracy of the representation of each field while simultaneously reducing computational cost compared to uniformly refined discretizations. A flexible mesh adaptation strategy enables local refinement based on geometric measures or physics-based error indicators. These adaptive discretization and analysis approaches are integrated into gradient-based optimization schemes, evaluating the design sensitivities using the adjoint method. Numerical studies illustrate the features of the proposed framework with static, linear elastic, multi-material, two- and three-dimensional problems. The examples provide insight into the effect of refining the design variable field on the optimization result and the convergence rate of the optimization process. Using coarse higher-order B-spline discretizations for level-set fields promotes the development of smooth designs and suppresses the emergence of small features. Moreover, adaptive mesh refinement for state variable fields results in a reduction of overall computational cost. Higher-order B-spline discretizations are especially interesting when evaluating gradients of state variable fields due to their higher inter-element continuity.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Physicochemical and Performance Characterization of Six Commercial Organic Solvent Nanofiltration Membranes

This work introduces a novel, gradient-free metamaterial design method based on Gaussian process regression to represent the density field of a unit cell. The dimension of the design space is determined by the covariance matrix dimension in the Gaussian process regression. We propose compressing this matrix using an autoencoder, enabling the decoder to generate the density field and effectively reduce the originally large design space to a lower-dimensional subspace. In this compressed space, we employ an active learning method, Bayesian Adaptive Direct Search (BADS), for efficient exploration of the design space. We demonstrate that for simple 2D designs aimed at maximizing unit cell stiffness, our method yields results comparable to those of standard topology optimization. Furthermore, we extend our approach to various mechanical problems, from linear elasticity to hyperelastic large deformation and elasto-plasticity under finite deformation, to 3D metamaterial design. This illustrates the method’s versatility and effectiveness across a range of applications.

Wu, Haoran↗
DOE OSTI2026

PDE-constrained high-order mesh optimization

Here, we present a novel framework for PDE-constrained r-adaptivity of high-order meshes. The proposed method formulates mesh movement as an optimization problem, with an objective function defined as a convex combination of a mesh quality metric and a measure of the accuracy of the PDE solution obtained via finite element discretization. The proposed formulation achieves optimized, well-defined high-order meshes by integrating mesh quality control, PDE solution accuracy, and robust gradient regularization. We adopt the Target-Matrix Optimization Paradigm to control geometric properties across the mesh, independent of the PDE of interest. To incorporate the accuracy of the PDE solution, we introduce error measures that control the finite element discretization error. The implicit dependence of these error measures on the mesh nodal positions is accurately captured by adjoint sensitivity analysis. Additionally, a convolution-based gradient regularization strategy is used to ensure stable and effective adaptation of high-order meshes. We demonstrate that the proposed framework can improve mesh quality and reduce the error by up to 10 times for the solution of Poisson and linear elasto-static problems. The approach is general with respect to the dimensionality, the order of the mesh, the types of mesh elements, and can be applied to any PDE that admits well-defined adjoint operators.

Computer science↗
DOE OSTI2025

On the Sampling-Based Computation of Nash Equilibria Under Uncertainty via the Nikaido–Isoda Function

We consider the computation of an equilibrium of a stochastic Nash equilibrium problem, where the player objectives are assumed to be L 0 -Lipschitz continuous and convex, given rival decisions with convex and closed player-specific feasibility sets. To address this problem, we consider minimizing a suitably defined value function defined using the Nikaido–Isoda function. Such an avenue does not necessitate either monotonicity properties of the concatenated gradient map or potentiality requirements on the game but does require a suitable regularity requirement under which a stationary point is a Nash equilibrium. We design and analyze a sampling-enabled projected-gradient-response method, reliant on inexact resolution of a player-level best-response subproblem. Here, by deriving suitable Lipschitzian guarantees on the value function, we derive both asymptotic guarantees for the sequence of generated iterates as well as rate and complexity guarantees for computing a stationary point by appropriate choices of the sampling rate and inexactness sequence.

Nikaido-Isoda function↗
DOE OSTI2024

On the Departure from Monin–Obukhov Surface Similarity and Transition to the Convective Mixed Layer

Large-eddy simulations are used to evaluate mean profile similarity in the convective boundary layer (CBL). Particular care is taken regarding the grid sensitivity of the profiles and the mitigation of inertial oscillations in the simulation spin-up. The nondimensional gradients Φ for wind speed and air temperature generally align with Monin–Obukhov similarity across cases but have a steeper slope than predicted within each profile. The same trend has been noted in several other recent studies. The Businger-Dyer relations are modified here with an exponential cutoff term to account for the decay in Φ to first-order approximation, yielding improved similarity from approximately 0.05z i to above 0.3z i , where z i is the CBL depth. The necessity for the exponential correction is attributed to an extended transition from surface scaling to zero gradient in the mixed layer, where the departure from Monin–Obukhov similarity may be negligible at the surface but becomes substantial well below the conventional surface layer height of 0.1 z i .

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Efficient proximal subproblem solvers for a nonsmooth trust-region method

In [R. J. Baraldi and D. P. Kouri, Mathematical Programming, (2022), pp. 1-40], we introduced an inexact trust-region algorithm for minimizing the sum of a smooth nonconvex and nonsmooth convex function. The principle expense of this method is in computing a trial iterate that satisfies the so-called fraction of Cauchy decrease condition—a bound that ensures the trial iterate produces sufficient decrease of the subproblem model. In this paper, we expound on various proximal trust-region subproblem solvers that generalize traditional trust-region methods for smooth unconstrained and convex-constrained problems. We introduce a simplified spectral proximal gradient solver, a truncated nonlinear conjugate gradient solver, and a dogleg method. Finally, we compare algorithm performance on examples from data science and PDE-constrained optimization.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

A randomized sketching trust-region secant method for low-memory dynamic optimization

The numerical solution of dynamic optimization problems is often limited by the memory required to store the state trajectory, which is used to evaluate the objective function and its derivatives. Recently, [R. Muthukumar et al., SIAM Journal on Optimization 31(2), pp. 1242–1275 (2021)] introduced a trust-region method for dynamic optimization that employs randomized sketching to compress the state trajectory, resulting in inexact derivative computations. By adaptively learning the sketch rank, the trust-region algorithm achieves rigorous convergence guarantees. Here, we extend this approach to use secant Hessian approximations. Due to the randomness introduced by the sketch, the traditional secant update formulae can produce poor Hessian approximations. In particular, the difference of two gradients, computed from two different sketches, may be inconsistent. To overcome this, we employ a sketched approximation of the Hessian application, in lieu of computing the gradient difference. We numerically demonstrate the improved stability of this approach on an example from PDE-constrained optimization.

dynamic optimization↗
DOE OSTI2025

Homogenization of Dendritic Structures and the High-Temperature Strength of the Refractory High-Entropy Alloy MoNbTaVW

MoNbTaVW, a pioneering refractory high-entropy alloy (RHEA), is renowned for its exceptional strength at elevated temperatures. Like most RHEAs, it solidifies into a dendritic microstructure with steep concentration gradients, necessitating heat treatment for equilibration. Steep concentration gradients typically complicate mechanical property analysis and modeling, key challenges for high-throughput alloy discovery. This study examines the effects of high-temperature homogenization (1800 °C for 8 hours) on the microstructure and mechanical properties of MoNbTaVW across temperatures up to 1200 °C. Postmortem microstructural analysis, coupled with chemical, crystallographic, and mechanical assessments, revealed that yield strength was surprisingly insensitive to homogenization, with less than a 3 pct difference between as-cast and heat-treated materials. Nanoindentation confirmed minimal nanohardness changes across the dendritic structure at room temperature. However, homogenization significantly enhanced high-temperature work hardening, producing higher peak compressive strength. Both as-cast and homogenized MoNbTaVW exhibited room-temperature strengths of 1400 MPa, exceeding previously reported values. These findings demonstrate that high-temperature treatment enables microstructural homogenization without compromising strength, a unique behavior among RHEAs. In conclusion, while not universally applicable, this insight highlights the importance of understanding microstructural development in experimental alloys and offers a pathway to design alloys that achieve homogeneous-like properties without heat treatment.

Rietema, C. J. [Lawrence Livermore National Labora↗
DOE OSTI2025

Effects of an External Magnetic Field on Keyhole Mode Laser Melting of 316 Stainless Steel

Keyhole-mode laser melting is an efficient method for joining or cutting large, thick components, but controlling keyhole depth and fluctuations has remained challenging. Applying an external magnetic field can control melt pool flows and indirectly influence keyhole morphology and dynamics. The induced Lorentz force, comprising Seebeck and damping components, plays a crucial role in the melt pool dynamics, depending on temperature gradient, flow rates, and magnetic field orientation and magnitude. This research investigates the effects of an external magnetic field on keyhole behavior during laser spot melting of 316 stainless steel using synchronized high-speed synchrotron X-ray and thermal imaging. Findings revealed that a longitudinal magnetic field (120 mT) increased keyhole depth but exacerbated lateral fluctuations, resulted in a 20% increase in the melt pool temperature gradient and a 27% decrease in cooling rate. Conversely, a transverse magnetic field (760 mT) reduced keyhole depth and improved porosity formation. The findings suggest that a decrease in keyhole depth correlates with a decrease in fluctuations, and vice versa. These insights enhance understanding of external magnetic fields’ impact on laser melting, with implications for improving part quality.

Alamdari, Aslan Bafahm [The Ohio State Univ., Colu↗
DOE OSTI2026

An Integrated Computational Materials Engineering (ICME) Approach to Design Nonlinear Transition Zones Between Dissimilar Metals

Current approaches to designing graded transition joints (GTJs) between dissimilar metals often rely on linear changes in both composition profiles and thickness of each sublayer. This increases fabrication cost and may not be optimal with respect to residual stress or the formation of undesirable phases. Here, in this study, GTJs between P91 ferritic/martensitic steel and 347H austenitic stainless steel were designed using Integrated Computational Materials Engineering (ICME) principles with nonlinear composition and length profiles. Guided by inputs from classical mechanics and CALPHAD predictions of carbon chemical potential, a novel transition zone consisting of five discrete compositions was proposed, with the thickness of each sublayer varying according to a brachistochrone-inspired distribution. In addition to carbon potential gradients, CALPHAD was used to predict coefficients of thermal expansion, which were incorporated into finite element models to evaluate stress evolution. The proposed nonlinear design resulted in a smoother carbon potential gradient, lower carbon depletion at the P91 interface, and a comparable residual stress under long-term thermal exposure, compared to a conventional linear design using ten sublayers with equal thickness. This work introduces a brachistochrone-inspired distribution for GTJ design, offering a general framework for optimizing graded interfaces between dissimilar metals.

Directed Energy Deposition↗
DOE OSTI2025

Effect of artificial viscosity on shocked particle-laden flows for staggered grid Lagrangian methods

Abstract Shocked particle-laden flows are important to many natural and industrial processes. When simulating these systems, artificial viscosity is often required to prevent numerical artifacts, such as ringing, from arising in the pressure and density fields. The linear and quadratic coefficients of the artificial viscosity determine the amount of smoothing that occurs in these fields. For particle-laden flows, however, many of the fluid–particle interaction forces, for example, the pressure gradient force and unsteady forces, depend on gradients in the fluid fields. Furthermore, while the shock passes over a particle, these forces can be more dominant than drag. This means that the artificial viscosity coefficients affect how a particle and fluid interact when simulating shocked particle systems. Here this effect is investigated for isolated particles and for a particle curtain using a staggered grid Lagrangian approach. The artificial viscosity coefficients have a significant impact on the maximum force that a fluid imparts to a particle, which is important for determining whether a particle will break up in response to the shock. Furthermore, it is found that the density ratio between the particle and the fluid is important in determining whether the artificial viscosity coefficients have a significant impact on the particle’s motion.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2026

Modeling the influence of the solid electrolyte interphase on the sand’s time and dendrite formation on lithium metal electrodes

Lithium metal is a sought after battery material for its high energy density due to the low electrochemical potential and density. However, lithium metal is also highly reactive, which results in a strong propensity for dendrite formation. The Sand’s time has previously been used to predict the time of dendrite initiation on metals that do not form a solid-electrolyte interphase (SEI), but it has been shown that the Sand’s time is not accurate for lithium electrodes when using transport parameters associated with the electrolyte. Thus, we built a numerical model to simulate lithium ion transport through a growing SEI to predict the Sand’s time. The numerical model is shown to be more accurate than previous analytical solutions, especially for low current densities. We then analyze the sensitivity of the Sand’s time to different SEI properties and the chemical potential gradients present in the SEI, driving lithium transport. The results showed that high lithium concentration has a greater impact at high current density, while fast diffusivity is more important at low current density. Lastly, we modeled the influence of surface roughness on the plating evolution and chemical potential gradients when an SEI is present in comparison to the electrolyte. As a result, we demonstrate that the SEI plays a critical role in lithium electrode stability, and that improved characterization techniques are needed to better understand transport through the SEI and increase lithium metal utilization in energy storage devices.

Chemistry↗
DOE OSTI2025

The rigorous upscaling of advection-dominated transport in heterogeneous porous media via the Method of Finite Averages

Systems involving advection-dominated transport through heterogeneous porous and fractured media are ubiquitous in subsurface engineering applications. However, upscaling such systems continues to challenge rigorous modeling efforts, particularly when advection is stronger than diffusion at fine spatial scales (i.e., when the Péclet number is greater than one at length scales that characterize a system’s unit-cells, representative elementary volumes, or averaging regions). Here, in this work, we propose and validate a strategy for extending the Method of Finite Averages (MoFA), a rigorous upscaling methodology for heterogeneous porous media, to upscale transport systems experiencing stronger advection than diffusion at fine scales (i.e., fine-scale Péclet numbers greater than one). We detail the strategy, the physical conditions under which it can be applied while retaining a priori modeling error guarantees, and implement the strategy to obtain a MoFA model for advective-diffusive transport that accommodates advective physics at fine spatial scales. We then perform two numerical experiments considering systems with system-scale Péclet numbers of 300 and 1000 — which correspond to fine-scale Péclet numbers of 30 and 100, respectively — to verify that the error guarantees are met under the strategy. After, we conduct a numerical study to demonstrate the strategy’s advantages over the original MoFA methodology. The results suggest that rigorously-upscaled transport models for heterogeneous porous media experiencing advective physics at finer spatial scales can be derived through MoFA and resolved orders of magnitude faster than their pore-scale counterparts. The results also suggest that the presented strategy is limited to modeling shallow concentration gradients when there are large differences between the time scales related to advection and a system’s temporally-varying boundary conditions. This limitation hinders the strategy’s practicality in modeling more advective systems, and as such, opportunity exists for developing additional strategies that accommodate rapidly-varying boundary conditions — and consequentially, steeper concentration gradients — while modeling advective systems with MoFA.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Reduced soil diazotroph diversity decreases nitrogen fixation rates, but depends on land management

Soil diazotrophs convert atmospheric nitrogen into plant-available ammonium through free-living nitrogen fixation (FLNF). This sustainable nitrogen source can reduce our dependence on synthetic fertilizer inputs in conventional agricultural systems. However, we know little about the effect of diazotroph diversity on FLNF, especially given that FLNF is intermediate within the broad-narrow functional spectrum. Here, we determined how management-mediated shifts in diazotroph diversity would impact their ecosystem function (FLNF) by quantifying diazotroph diversity across a long-term management gradient during and after the growing season. In addition to field observations, we leveraged the same management gradient to manipulate diversity in soil microcosms via chloroform fumigation exposure. In the field, diazotroph diversity was significantly higher after the growing season, and the biologically-based annual cropping system harbored the highest diazotroph diversity. However, perennial cropping systems maintained the highest FLNF despite lower diazotroph diversity, and both soil moisture and temperature were stronger predictors of FLNF. Based on these results, integrating diverse perennial crops into agricultural landscape could result in greater N from FLNF, particularly at the end of the growing season. When we reduced biodiversity in a manipulation experiment, the diversity-FLNF association was stronger than in the field experiment, suggesting that FLNF communities are not as functionally redundant as taxonomically ‘broad’ ecosystem functions. Here, the strength of diversity-FLNF correlation varied by previous land management. Diazotroph diversity better predicted FLNF in annual and forest soil microcosms, and microbial biomass carbon better predicted FLNF in perennial soil microcosms. Taken together, our results show that while diazotroph diversity influences FLNF, especially under extreme environmental disturbances, abiotic factors like soil moisture and temperature are stronger constraints on FLNF in the field.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2025

Mapping tree height in complex terrain of northern China using ultra-high-resolution images

Tree height is a key parameter for estimating forest biomass and carbon sequestration. In recent years, notable progress has been made in mapping tree height using satellite imagery. However, existing tree height products show low accuracy in mountainous and complex terrains, and few studies typically addressed tree height estimations in mountain areas. This study examines the Mentougou district of Beijing, China, characterized by complex terrain and mountainous landscapes. We analyzed two methods for estimating tree height: one using only spectral features and another combining spectral features with topographic factors (elevation, slope, aspect). We used 3-m resolution PlanetScope 8-band multispectral imagery, with 710 field-measured individual tree heights averaged to obtain 471 pixel-level tree height values as ground-truth, to develop tree height prediction models using eXtreme Gradient Boosting (XGBoost), Random Forest (RF), and Gradient Boosting Machine (GBM) models. The results show that the XGBoost model consistently presented the highest accuracy for both methods evaluated. Specifically, the XGBoost model that combined spectral data with elevation and slope variables with an R² of 0.75 and an RMSE of 2.69 m. Using the XGBoost model, we generated the tree height map for the Mentougou area at 3 m resolution, showing tree heights ranging from 0.5 to 30.4 m, and the model’s prediction error standard deviations ranged from 2.50 to 4.71 m, indicating reliable performance across varied terrain. Additionally, we compared and evaluated the global tree height products, identifying limitations in the accuracy within complex terrains. This study demonstrates the potential for accurately predicting tree heights by combining high-resolution multispectral satellites with a terrain factor modeling approach.

Complex terrain↗
DOE OSTI2025

LABQ3: Bayesian method for quantification of mineral compositions and nano-scale elemental mapping of 3D synchrotron XCT data

Quantitative analysis of mineral compositions is essential in understanding geochemical, mineralogical and environmental processes. Fine-resolution 3D imaging is widely done using synchrotron X-ray computed tomography (XCT), but existing analyses are limited to visualization and segmentation. This paper presents a new method, Linear Attenuation Bayesian Quantitative 3D-mapper (LABQ3), based on the linearity of X-ray attenuation with respect to elemental concentrations. To address the random variability in attenuation measurements, LABQ3 employs Bayesian decision theory to minimize classification error, using reference attenuation distributions from scans of pure mineral standards. To demonstrate LABQ3 and test its performance, we studied precipitated carbonate samples. XCT scans were done at multiple energies using the transmission X-ray microscope (TXM) at beamline 32-ID-C of the Advanced Photon Source at Argonne National Laboratory. The reconstructed 3D images have a voxel size of 20 nm. Analyses revealed rich nano-scale compositional heterogeneity within individual particles. A mixture of calcium and cadmium produced an overall stoichiometric composition of (Ca 0.78 ,Cd 0.22 )CO 3 , with some voxels containing nearly pure CdCO 3 . The addition of zinc led to an overall stoichiometric composition of 33% Ca, 28% Cd, 39% Zn, with a nearly pure CaCO 3 core and compositional zonation through the rim. These compositional gradients are related to temporal sequences of carbonate mineral formation where Cd precipitated at the beginning in (Ca,Cd)CO 3 , while Cd and Zn precipitated at the end in (Ca, Cd,Zn)CO 3 . Results differ from bulk analyses using Inductively Coupled Plasma-Mass Spectrometry (ICP-MS), showing that LABQ3 provides particle-specific insights. LABQ3 distinguishes itself by quantifying chemical compositions along a continuum, making it different from XCT analyses based on segmentation. LABQ3 allows simultaneous acquisition of morphology and chemical composition in 3D, facilitating the interpretation of chemical gradients of trace elements, quantification of solid solution compositions, inferences about temporal sequences of mineral precipitation, and addressing other concerns about solid-phase chemistry.

58 GEOSCIENCES↗
DOE OSTI2024

Analytically differentiable metrics for phase stability

Here, in this work, a long-established but sparsely documented method of obtaining semi-analytic derivatives of thermodynamic properties with respect to equilibrium conditions is briefly reviewed and rigorously derived. This procedure is then leveraged to construct general forms of derivatives of the residual driving force, a metric for measuring phase stability used in CALPHAD model optimization, with respect to overall system and individual phase compositions. Applied examples – calculating heat capacity in the Al-Fe system, thermodynamic factors in the Nb-V-W system, and residual driving force derivatives in the Ni-Ti system – demonstrate the versatility, accuracy, and extensibility of this method. Using the developed method, residual driving force gradients can be applied directly in CALPHAD model optimizers, as well as in materials design frameworks, to identify regions of phase stability with an efficient, gradient-based approach.

36 MATERIALS SCIENCE↗