Search NASA⌕ Search

SEARCH · Search NASA

Results for “Closed form solution”

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

Mineral Scaling in 3D Interfacial Solar Evaporators–A Challenge for Brine Treatment and Lithium Recovery

In this work, we analyzed the effects of mineral scaling on the performance of a 3D interfacial solar evaporator, with a focus on the cations relevant to lithium recovery from brackish water. The field has been rapidly moving toward resource recovery applications from brines with higher cation concentrations. However, the potential complications caused by common minerals in these brines other than NaCl have been largely overlooked. Therefore, in this study, we thoroughly examined the effects of two common cations (calcium and magnesium) on the long-term solar evaporation performance of a 3D graphene oxide stalk. The 3D stalk can achieve an evaporation flux as high as 17.8 kg m –2 h –1 under one-sun illumination, and accumulation of NaCl on the stalk surface has no impact. However, the presence of CaCl 2 and MgCl 2 significantly reduces the evaporative flux even in solutions lacking scale-forming anions. A close examination of scale formation during long-term evaporation experiments revealed that CaCl 2 and MgCl 2 tend to precipitate out within the stalk, thus blocking water transport through the stalk and significantly dropping the evaporation rates. These findings imply that research attention is needed to modify and optimize the internal water transport channels for 3D evaporators. Additionally, we emphasize the importance of testing realistic mixtures–including prominent divalent cations– and testing long-term operations in interfacial solar evaporation research and investigating approaches to mitigate the negative impacts of divalent cations.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2021

Solving the sample size problem for resource selection functions

Abstract Sample size sufficiency is a critical consideration for estimating resource selection functions (RSFs) from GPS‐based animal telemetry. Cited thresholds for sufficiency include a number of captured animals and as many relocations per animal N as possible. These thresholds render many RSF‐based studies misleading if large sample sizes were truly insufficient, or unpublishable if small sample sizes were sufficient but failed to meet reviewer expectations. We provide the first comprehensive solution for RSF sample size by deriving closed‐form mathematical expressions for the number of animals M and the number of relocations per animal N required for model outputs to a given degree of precision. The sample sizes needed depend on just 3 biologically meaningful quantities: habitat selection strength, variation in individual selection and a novel measure of landscape complexity, which we define rigorously. The mathematical expressions are calculable for any environmental dataset at any spatial scale and are applicable to any study involving resource selection (including sessile organisms). We validate our analytical solutions using globally relevant empirical data including 5,678,623 GPS locations from 511 animals from 10 species (omnivores, carnivores and herbivores living in boreal, temperate and tropical forests, montane woodlands, swamps and Arctic tundra). Our analytic expressions show that the required M and N must decline with increasing selection strength and increasing landscape complexity, and this decline is insensitive to the definition of availability used in the analysis. Our results demonstrate that the most biologically relevant effects on the utilization distribution (i.e. those landscape conditions with the greatest absolute magnitude of resource selection) can often be estimated with much fewer than animals. We identify several critical steps in implementing these equations, including (a) a priori selection of expected model coefficients and (b) regular sampling of background (pseudoabsence) data within a given definition of availability. We discuss possible methods to identify a priori expectations for habitat selection coefficients, effects of scale on RSF estimation and caveats for rare species applications. We argue that these equations should be a mandatory component for all future RSF studies.

Street, Garrett M.↗
DOE OSTI2026

Bayesian Physics Informed Spatio-Temporal Network for Streamflow Data Imputation

Reliable reconstruction of incomplete streamflow records is critical for improving hydrological forecasting, flood preparedness, and water resource management. However, large observational gaps and uncertainties in governing physical parameters limit the accuracy of traditional statistical and machinelearning imputation frameworks. To address these challenges, we develop a Bayesian Physics-Informed Spatio-Temporal Network (BPI-STNet) that jointly captures spatial and temporal dependencies while enforcing hydrologic consistency through embedded physical constraints. The framework integrates a GraphSAGE-LSTM architecture to model spatial connectivity across gauges and temporal flow dynamics, coupled with a Bayesian update mechanism to estimate uncertain parameters in a simplified water-balance framework. Unlike conventional physics-informed networks that rely on sampling-based posterior estimation, BPI-STNet derives an analytic solution to the inverse problem, allowing closed-form Bayesian updates of uncertain parameters Λ={α,β,k} using Gaussian priors and likelihoods. Applied to daily observations from the Susquehanna River Basin (1980-2022), BPI-STNet achieves substantial improvements over a purely data-driven RGNN baseline, which reduced RMSE by 23 % and MAE by 9 %, and achieving an average NSE values up to 0.96. The results demonstrate that coupling Bayesian inference with physics-informed learning yields physically consistent, uncertainty-aware reconstructions that preserve the temporal persistence and statistical distribution of observed flows. The proposed framework establishes a generalizable paradigm for data-sparse hydrologic systems where both data fidelity and physical interpretability are essential.

Krishnan Kutty Ambika, Anukesh [ORNL] (ORCID:00000↗
DOE OSTI2025

Dynamics of metastable contact soliton dissipative exchange flows in one-dimensional ferromagnetic channels

Dissipative exchange flows (DEFs) are large-amplitude boundary value solutions of ferromagnetic channels. In their low-injection limit, DEFs reduce to spin superfluids. However, in the strong injection limit, nonlinearities dominate close to the injection site and a soliton is formed; this solution has been termed a contact soliton dissipative exchange flow (CS-DEF). Here, in this work, we numerically investigate CS-DEF solutions in a moderate injection regime and a finite injection width. We find a solution where two metastable solitons coexist in the injection region. This solution is metastable in the sense that any perturbation to the system will eject one of the solitons out of the injection region. Moreover, soliton dynamics can be excited when two injection regions are separated by a certain distance. We find that the ensuing DEF between the solitons induces a steady-state dynamics in which metastable solitons are continually ejected and nucleated. Furthermore, and depending on the relative signs of the spin injections, the soliton dynamics possess a particular handedness and frequency related to the spin transfer torque delivered by the DEF. Our results provide insights into the transport of spin current by DEFs - where the interaction between DEFs and solitons suggests a mechanism for detaching contact-solitons from the injection boundary. Although this study focuses on the "nonlocal" interaction between solitons, it may lead to the investigation of new mechanisms for inserting solitons in a DEF, e.g., for discrete motion and transport of information over long distances.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Electric and magnetic black holes in a new nonlinear electrodynamics model

Highlights: • We introduce a new NED model which is comparable with the BI model in the weak-field limit. • We find an electric black hole solution in the context of Einstein’s gravity minimally coupled with the new NED model. • We study the thermal stability of the electric black hole solution. • Modified Smarr’s formula consistent with the first law of black hole thermodynamics is obtained. A new nonlinear electrodynamics (NED) model is introduced in the form of a nonpolynomial Lagrangian which admits static spherical and also plane wave solutions in a flat space. Upon coupling with gravity the electric field is finite and comparable with the Born–Infeld counterpart. The electric and magnetic black hole solutions in the Einstein’s gravity coupled with this NED model are presented. The solutions give both asymptotically and in the weak field limit Reissner–Nordström (RN) black hole and unlike the other known models our electric solution is expressed in terms of elementary functions in a closed form. We study the first law and derive the modified Smarr’s formula for the electric type extension of our model. Considerable rich structure, especially thermodynamic ones, ranging from first to the second order phase transitions are added to the RN black hole of linear electrodynamics with this NED model. Having the exact solution for the metric function at our disposal we investigate the stability of the electric black hole from both the thermodynamical and causal points of view.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Online State Estimation for Time-Varying Systems

The paper investigates the problem of estimating the state of a time-varying system with a linear measurement model; in particular, the paper considers the case where the number of measurements available can be smaller than the number of states. In lieu of a batch linear least-squares (LS) approach well-suited for static networks, where a sufficient number of measurements could be collected to obtain a full-rank design matrix the paper proposes an online algorithm to estimate the possibly time-varying state by processing measurements as and when available. The design of the algorithm hinges on a generalized LS cost augmented with a proximal-point-type regularization. With the solution of the regularized LS problem available in closed-form, the online algorithm is written as a linear dynamical system where the state is updated based on the previous estimate and based on the new available measurements. Conditions under which the algorithmic steps are in fact a contractive mapping are shown, and bounds on the estimation error are derived for different noise models. Numerical simulations are provided to corroborate the analytical findings.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2021

Quasi-linear theory of forced magnetic reconnection for the transition from the linear to the Rutherford regime

Abstract Using the in-viscid two-field reduced magneto-hydrodynamic (MHD) model, a new analytical theory is developed to unify the Hahm–Kulsrud–Taylor (HKT) linear solution and the Rutherford quasi-linear regime. Adopting a quasi-linear approach, we obtain a closed system of equations for plasma response in the Taylor problem. An integral form of an analytical solution is obtained for the forced magnetic reconnection, uniformly valid throughout the entire regimes from the HKT linear solution to the Rutherford quasi-linear solution. In particular, the quasi-linear effect can be described by a single coefficient K s ∝ S 8 / 5 ψ c 2 , where S = τ R τ A and ψ c are the Lundquist number and amplitude of external magnetic perturbation, respectively. The HKT linear solution for the response can be recovered when the index K s → 0. On the other hand, the quasi-linear effect plays a key role in the island growth when K s ~ 1. Our new analytical solution has also been compared with reduced MHD simulations with agreement.

Physics↗
DOE OSTI2022

Enhanced crevice corrosion of stainless steel 316 by degradation of Cr-containing hollandite crevice former

This study explores the corrosion interactions between stainless steel 316 (SS316), a candidate nuclear waste canister material, and Cr-containing hollandite (Ba 1.15 Cr 2.3 Ti 5.7 O 16 ), a model ceramic nuclear waste form. The two materials were corroded in close proximity in 0.6 M NaCl solution at 90 °C for 28 days. Severe crevice corrosion was observed on the stainless steel surface after the corrosion experiment, which is likely associated with the rapid accumulation of Cr 3+ cations originating from both SS316 and the Cr-containing hollandite. Furthermore, the enhanced corrosion of SS316 by the corrosion products of hollandite is quantitatively explained using a classical crevice corrosion model.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Reactive particle-tracking solutions to a benchmark problem on heavy metal cycling in lake sediments

Geochemical systems are known to exhibit highly variable spatiotemporal behavior. This may be observed both in non-smooth concentration curves in space for a single sampling time and also in variability between samples taken from the same location at different times. However, most models that are designed to simulate these systems provide only single-solution smooth curves and fail to capture the noise and variability seen in the data. We apply a recently developed reactive particle-tracking method to a system that displays highly complex geochemical behavior. When the method is made to most closely resemble a corresponding Eulerian method, in its unperturbed form, we see near-exact match between solutions of the two models. More importantly, we consider two approaches for perturbing the model and find that the spatially-perturbed condition is able to capture a greater degree of the variability present in the data. This method of perturbation is a task to which particle methods are uniquely suited and Eulerian models are not well-suited. Additionally, because of the nature of the algorithm, noisy spatial gradients can be highly resolved by a large number of mobile particles, and this incurs negligible computational cost, as compared to expensive chemistry calculations.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

The mechanism for acetate formation in electrochemical CO (2) reduction on Cu: selectivity with potential, pH, and nanostructuring

Nanostructured Cu catalysts have increased the selectivities and geometric activities for high value C–C coupled (C 2 ) products (ethylene, ethanol, and acetate) in the electrochemical CO (2) reduction reaction (CO (2) RR). The selectivity among the high-value C 2 products is also altered, where for instance the yield of acetate increases with alkalinity and is dependent on the catalyst morphology. The reaction mechanisms behind the selectivity towards acetate vs. other C 2 products remain controversial. In this work, we elucidate the reaction mechanism for acetate formation by using ab initio simulations, a coupled kinetic-transport model, and loading dependent experiments. We find that trends in acetate selectivity can be rationalized from variations in electrolyte pH and the local mass transport properties of the catalyst and not from changes in Cu's intrinsic activity. The selectivity mechanism originates from the transport of ketene, a stable (closed shell) intermediate, away from the catalyst surface into solution where it reacts to form acetate. While this type of mechanism has not yet been discussed in the CO (2) RR, variants of it may explain similar selectivity fluctuations observed for other stable intermediates like CO and acetaldehyde. Our proposed mechanism suggests that acetate selectivity increases with increasing pH, decreasing catalyst roughness and significantly varies with the applied potential.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2025

Shapes of ideal stalagmites

Stalagmites are isolated columns of calcium carbonate growing on a cave floor; their growth is driven by the constant dripping of supersaturated solutions from the roof of the cave. In this paper, we derive a closed-form expression for the shape of a steadily growing stalagmite. Our analysis gives rise to three distinct shapes, all of them observable in nature, with the shape characterized by a single dimensionless parameter. Transitions between different shapes occur at a specific value of this parameter, with additional selection rules determining the shape and size of stalagmites evolving under specific cave conditions. Our theory shows that the stalagmite shape influences the 13 C isotope shifts, which are an important source of paleoclimatic information.

invariant growth↗
DOE OSTI2025

Multipoint Correlations in Poisson Media

Multipoint correlations in randomly heterogeneous composite media are determined by the probability that a set of points belong to specific phases. They determine a wide range of macroscopic transport properties such as conductivity, dielectric constant, diffusion coefficient, and transmittance. The Poisson model—a random tesselation of space by hyperplanes—provides realistic descriptions of heterogeneous media in, e.g., radiation-transport applications; yet, until now, it has lacked closed-form expressions for its multipoint correlations. We resolve this problem by presenting an exact solution for the multipoint correlations in the Poisson model. By comparing it to Monte Carlo simulations of four-point correlations in three dimensions, we demonstrate the accuracy of our solution. In conclusion, we visualize the multipoint correlations and discuss their features.

Amorphous materials↗
DOE OSTI2024

Finite domain solution of a KGD hydraulic fracture in the viscosity-dominated regime

This paper describes a numerical algorithm for solving the classic problem of a plane strain (KGD) fracture propagating in an impermeable elastic medium with zero toughness. The method, which takes advantage of the self-similar nature of the solution, combines a domain-based scheme to solve the elasticity equations and a finite volume method to solve the nonlinear lubrication equation. This work represents a first step towards developing a model able to account for pore pressure diffusion in the medium and corresponding poroelastic effects, noting that these processes are more efficiently solved using a domain-based rather than a boundary integral method. To enhance the efficiency and accuracy of the numerical scheme, the far-field crack asymptotics is embedded in the discretized elastic relationship between the fluid pressure and the crack opening, while the coupled fluid-solid tip asymptote is enforced in a weak form when solving the nonlinear lubrication equation. The proposed technique yields results that closely match the analytical solution, even with a coarse mesh. This approach offers potential for addressing more complex hydraulic fracturing problems in the future.

Domain-based method↗
DOE OSTI2020

Data-driven learning of nonlocal models: from high-fidelity simulations to constitutive laws

We show that machine learning can improve the accuracy of simulations of stress waves in one-dimensional composite materials. We propose a data-driven technique to learn nonlocal constitutive laws for stress wave propagation models. The method is an optimization-based technique in which the nonlocal kernel function is approximated via Bernstein polynomials. The kernel, including both its functional form and parameters, is derived so that when used in a nonlocal solver, it generates solutions that closely match high-fidelity data. The optimal kernel therefore acts as a homogenized nonlocal continuum model that accurately reproduces wave motion in a smaller-scale, more detailed model that can include multiple materials. We apply this technique to wave propagation within a heterogeneous bar with a periodic microstructure. Several one-dimensional numerical tests illustrate the accuracy of our algorithm. The optimal kernel is demonstrated to reproduce high-fidelity data for a composite material in applications that are substantially different from the problems used as training data.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Pairing a Global Optimization Algorithm with EXAFS to Characterize Lanthanide Structure in Solution

Ensemble-average sampling of structures from ab initio molecular dynamics (AIMD) simulations can be used to predict theoretical extended X-ray absorption fine structure (EXAFS) signals that closely match experimental spectra. However, AIMD simulations are time-consuming and resource-intensive, particularly for solvated lanthanide ions, which often form multiple nonrigid geometries with high coordination numbers. Here, to accelerate the characterization of lanthanide structures in solution, we employed the Northwest Potential Energy Surface Search Engine (NWPEsSe), an adaptive-learning global optimization algorithm, to efficiently screen first-shell structures. As case studies, we examine two systems: Eu(NO 3 ) 3 dissolved in acetonitrile with a terpyridine ligand (terpyNO 2 ), and Nd(NO 3 ) 3 dissolved in acetonitrile. The theoretical spectra for structures identified by NWPEsSe were compared to both experimental and AIMD-derived EXAFS spectra. The NWPEsSe algorithm successfully identified the proper solvation structure for both Eu(NO 3 ) 3 (terpyNO 2 ) and Nd(NO 3 )(acetonitrile) 3 , with the calculated EXAFS signals closely matching the experimental spectra for the Eu-ligand complex and showing good similarity for the Nd salt; the better agreement with the ligand-containing structure is attributed to a less dynamic coordination environment due to the rigid ligand. The key advantage of the global optimization algorithm lies in its ability to sample the coordination environment across the potential energy surface and reduce the time required to identify structures from generally a month to within a week. Additionally, this approach is versatile and can be adapted to characterize main-group metal complexes.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2020

Corrosion interactions between stainless steel and lead vanado-iodoapatite nuclear waste form part I

This study explores the corrosion interactions between a metallic canister material, stainless steel (SS) 316, and an I-bearing ceramic waste form, lead vanado-iodoapatite (I-APT, Pb 9.85 (VO 4 ) 6 I 1.7 ), in a chloride solution. Crevice corrosion of the SS in close proximity to the I-APT resulted in the development of an aggressive environment at the interface of the two materials, which was acidic and enriched in Cl – anions. I-APT also corroded in the crevice region, primarily through ion-exchange between the I – ions from the I-APT matrix and anions from the environment. The enrichment of Cl – anions within the occluded crevice space as the result of SS crevice corrosion enhanced the corrosion of I-APT. The release of iodine from this apatite waste form could be accelerated owing to this mechanism. This is evidenced by a depletion of iodine from the I-APT matrix and a large amount of Cl-bearing precipitates on the surfaces of both SS and I-APT. On the other hand, the corrosion of I-APT leads to the precipitation of a V- and Pb-rich layer, which inhibits the localized corrosion of SS to an extent. This study advances the understanding of the near-field corrosion interactions between metallic canisters and ceramic waste forms.

36 MATERIALS SCIENCE↗
DOE OSTI2021

The superatomic state beyond conventional magic numbers: Ligated metal chalcogenide superatoms

The field of cluster science is drawing increasing attention due to the strong size and composition-dependent properties of clusters and the exciting prospect of clusters serving as the building blocks for materials with tailored properties. However, identifying a unifying central paradigm that provides a framework for classifying and understanding the diverse behaviors is an outstanding challenge. One such central paradigm is the superatom concept that was developed for metallic and ligand-protected metallic clusters. The periodic electronic and geometric closed shells in clusters result in their properties being based on the stability they gain when they achieve closed shells. This stabilization results in the clusters having a well-defined valence allowing them to be classified as superatoms – thus, extending the periodic table to a third dimension. This perspective focuses on extending the superatomic concept to ligated metal-chalcogen clusters that have recently been synthesized in solutions and form assemblies with counterions that have wide-ranging applications. Here we illustrate that the periodic patterns emerge in the electronic structure of ligated metal-chalcogenide clusters. The stabilization gained by the closing of their electronic shells allows for the prediction of their redox properties. Further investigations reveal how the selection of ligands may control the redox properties of the superatoms. These ligated clusters may serve as chemical dopants for two-dimensional semiconductors to control their transport characteristics. Superatomic molecules of multiple metal-chalcogen superatoms allow for the formation of nano pn junctions ideal for directed transport and photon harvesting. As a result, the perspective outlines future developments, including the synthesis of magnetic superatoms.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Asymptotic Expansion of the Impedance Per Unit Length for Rectangular Conductors

An iteration method is introduced to obtain the asymptotic form of the impedance per unit length of a rectangular conductor when the half side lengths are large compared to the skin depth. The first terms of the asymptotic expansion are extracted in closed form. The manner in which the corner corrections fit into the expansion are illustrated. The asymptotic results are compared to a numerical solution in the square limit. The odd corner correction for a right angle edge is also discussed.

42 ENGINEERING↗