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 883 records · Page 49
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Combined Precursor Engineering and Grain Anchoring Leading to MA‐Free, Phase‐Pure, and Stable α‐Formamidinium Lead Iodide Perovskites for Efficient Solar Cells

Abstract α‐Formamidinium lead iodide (α‐FAPbI 3 ) is one of the most promising candidate materials for high‐efficiency and thermally stable perovskite solar cells (PSCs) owing to its outstanding optoelectrical properties and high thermal stability. However, achieving a stable form of α‐FAPbI 3 where both the composition and the phase are pure is very challenging. Herein, we report on a combined strategy of precursor engineering and grain anchoring to successfully prepare methylammonium (MA)‐free and phase‐pure stable α‐FAPbI 3 films. The incorporation of volatile FA‐based additives in the precursor solutions completely suppresses the formation of non‐perovskite δ‐FAPbI 3 during film crystallization. Grains of the desired α‐phase are anchored together and stabilized when 4‐ tert ‐butylbenzylammonium iodide is permeated into the α‐FAPbI 3 film interior via grain boundaries. This cooperative scheme leads to a significantly increased efficiency close to 21 % for FAPbI 3 perovskite solar cells. Moreover, the stabilized PSCs exhibit improved thermal stability and maintained ≈90 % of their initial efficiency after storage at 50 °C for over 1600 hours.

Ling, Xufeng↗
DOE OSTI2022

High entropy alloys with hexagonal close-packed structure derived from thin film combinatorial approach

A thin film combinatorial synthesis approach was utilized to develop a compositional and structural library for the OsRuWCo alloy system. The survey of this alloy space produced 24 unique compositions and their corresponding crystal structures, which represent a phase map where amorphous and single-phase hexagonal close-packed alloys were identified. Based on the identification of a promising candidate film region, a new high entropy alloy based on OsRuWCo, but with non-equiatomic composition, was synthesized in bulk form. The alloy exhibited a single-phase hexagonal close-packed structure in the as-cast state. Three derivatives from this system were also developed by considering the heat of mixing, atomic size, and binary solubility. These new alloys are based on OsRuWCoIr and OsRuWCoFe, both of which exhibit single-phase hexagonal close-packed as-cast structures, as determined by X-ray diffraction and electron microscopy. Additionally, this large compositional space was utilized to evaluate conventional parameters that describe high entropy alloys (ΔS mix , ΔH mix , Ω, δ, and VEC), and these were found to elucidate the criteria for solid solution formation in this system. Additionally, trends illustrating the evolution from amorphous to crystalline phases are demonstrated and discussed. It is concluded that in this particular case, the dominant factors promoting crystallinity are thermodynamic or electronic, rather than geometric.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Spatial Correlations of the Poisson Model for Radiation Transport

Characterizing the relationship between bulk physical properties and mixing in randomly heterogeneous media is a central challenge across many areas of science and engineering. A benchmark model for such studies is the Poisson model, a random tessellation of space by a Poisson process of hyperplanes. In radiation transport studies, the lack of exact expressions for the Poisson model’s spatial multipoint functions has led to approximate methods being used, introducing unquantified sources of error. Here, we recently introduced an exact solution for the Poisson model’s multipoint functions and closely related conditional probability functions (CPFs), providing a new opportunity to understand and reduce these sources of error. In this paper, we enable a more rigorous investigation of radiation transport in stochastic media by applying the recently introduced exact solution for the Poisson model’s CPFs. This paper consists of three main contributions. First, we introduce a unified framework for CPFs of the Poisson model, encompassing the recently introduced exact CPFs as well as the previously introduced atomic mix, nearest-neighbor, and combination CPFs. This framework also includes existing pruning techniques for the approximate CPFs, such as angular exclusion, as well as a novel form of angular exclusion suitable for the exact CPFs. Second, we use the exact CPFs to characterize the spatial regions where each approximate three-point CPF is most accurate, thereby explaining the observed hierarchy of accuracy among the approximate models. Finally, we evaluate material transmittance, reflectance, and flux in a three-dimensional test problem using conditional point sampling, demonstrating the relationship between CPF accuracy and transport simulation accuracy.

Poisson model↗
NASA NTRS1974

A finite-element analysis for steady and oscillatory subsonic flow around complex configurations

The problem of potential subsonic flow around complex configurations is considered. The solution is given of an integral equation relating the values of the potential on the surface of the body to the values of the normal derivative, which is known from the boundary conditions. The surface of the body is divided into small (hyperboloidal quadrilateral) surface elements, which are described in terms of the Cartesian components of the four corner points. The values of the potential (and its normal derivative) within each element is assumed to be constant and equal to its value at the centroid of the element. The coefficients of the equation are given by source and doublet integrals over the surface elements. Closed form evaluations of the integrals are presented. The results obtained with the above formulation are compared with existing analytical and experimental results.

Chen, L. T.↗
DOE OSTI2023

Fitting Matérn smoothness parameters using automatic differentiation

The Mat$\acute{e}$rn covariance function is ubiquitous in the application of Gaussian processes to spatial statistics and beyond. Perhaps the most important reason for this is that the smoothness parameter $\nu$ gives complete control over the mean-square differentiability of the process, which has significant implications for the behavior of estimated quantities such as interpolants and forecasts. Unfortunately, derivatives of the Mat$\acute{e}$rn covariance function with respect to $\nu$ require derivatives of the modified second-kind Bessel function $K$ $\nu$ with respect to $\nu$. While closed form expressions of these derivatives do exist, they are prohibitively difficult and expensive to compute. For this reason, many software packages require fixing $\nu$ as opposed to estimating it, and all existing software packages that attempt to offer the functionality of estimating $\nu$ use finite difference estimates for $\partial$ $\nu$ $K$ $\nu$ . In this work, we introduce a new implementation of $K$$\nu$ that has been designed to provide derivatives via automatic differentiation (AD), and whose resulting derivatives are significantly faster and more accurate than those computed using finite differences. Here, we provide comprehensive testing for both speed and accuracy and show that our AD solution can be used to build accurate Hessian matrices for second-order maximum likelihood estimation in settings where Hessians built with finite difference approximations completely fail.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Semi-transparent p-type barium copper sulfide as a back contact interface layer for cadmium telluride solar cells

Optically transparent p-type materials play a critical role in transparent electronics including photovoltaic (PV) devices. P-type sulfide materials offer an alternative to oxides for PV application due to improved hole transport properties. Here, we report the solution-based synthesis of earth-abundant p-type transparent conducting barium copper sulfide (α-BaCu 4 S 3 , BCS) thin films. These films were characterized using scanning electron microscopy, X-ray diffraction, UV–Vis–NIR spectrophotometry, Raman spectroscopy, and spectroscopic ellipsometry. BCS films of ~100 nm thickness transmit >70% of visible light. We report on tests of the hole transport properties of these BCS films for cadmium telluride (CdTe) photovoltaics, finding that the BCS deposition process forms a beneficial tellurium (Te) rich surface on CdTe by selectively removing Cd from the surface. Based on our study, the BCS interface layer plays dual functions for CdTe PV devices as a hole transport material and as an etchant, enhancing the resulting device performance. We observed a significant increase in open-circuit voltage of CdTe solar cells with the BCS buffer layer. Furthermore, we discuss semitransparent CdTe solar cells with BCS as a hole transport layer and indium tin oxide as a finishing electrode. Semitransparent CdTe solar cells shows 13.3% conversion efficiency for the front side illumination and 1.2% efficiency for back side illumination, indicating high recombination of charge carriers generated close to the rear CdTe/BCS/ITO contact.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Solvent Dynamics Are Critical to Understanding Carbon Dioxide Dissolution and Hydration in Water

Simulations of carbon dioxide (CO 2 ) in water may aid in understanding the impact of its accumulation in aquatic environments and help advance technologies for carbon capture and utilization (via, e.g., mineralization). Quantum mechanical (QM) simulations based on static molecular models with polarizable continuum solvation poorly reproduce the energetics of CO 2 hydration to form carbonic acid in water, independent of the level of QM theory employed. Only with density-functional-theory-based molecular dynamics and rare-event sampling, followed by energy corrections based on embedded correlated wavefunction theory (in conjunction with density functional embedding theory), can a close agreement between theory and experiment be achieved. Such multilevel simulations can serve as benchmarks for simpler, less costly models, giving insight into potential errors of the latter. The strong influence of sampling/averaging over dynamical solvent configurations on the energetics stems from the difference in polarity of both the transition state and product (both polar) versus the reactant (nonpolar). When a solute undergoes a change in polarity during reaction, affecting its interaction with the solvent, careful assessment of the energetic contribution of the solvent response to this change is critical. Here, we show that static models (without structural sampling) that incorporate three explicit water molecules can yield far superior results than models with more explicit water molecules because fewer water molecules yield less configurational artifacts. Static models intelligently incorporating both explicit (molecules directly participating in the reaction) and implicit solvation, along with a proper QM theory, e.g., CCSD(T) for closed-shell systems, can close the accuracy gap between static and dynamic models.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS1961

Estimation of Rotary Stability Derivatives at Subsonic and Transonic Speeds

The first part of this paper pertains to the estimation of subsonic rotary stability derivatives of wings. The unsteady potential flow problem is solved by a superposition of steady flow solutions. Numerical results for the damping coefficients of triangular wings are presented as functions of aspect ratio and Mach number, and are compared with experimental results over the Mach number range 0 to 1. In the second part, experimental results are used. to point out a close correlation between the nonlinear variations with angle of attack of the static pitching-moment curve slope and the damping-in-pitch coefficient. The underlying basis for the correlation is found as a result of an analysis in which the indicial function concept and. the principle of super-position are adapted to apply to the nonlinear problem. The form of the result suggests a method of estimating nonlinear damping coefficients from results of static wind-tunnel measurements.

Tobak, Murray↗
DOE OSTI2025

Model-based analysis of solute transport and potential carbon mineralization in the active layer of a hillslope underlain by permafrost with seasonal variability and climate change

Permafrost carbon, stored in frozen organic matter across vast Arctic and sub-Arctic regions, represents a substantial and increasingly vulnerable carbon reservoir. As global temperatures rise, the accelerated thawing of permafrost releases greenhouse gases, exacerbating climate change. However, freshly thawed permafrost carbon may also experience lateral transport by groundwater flow to surface water recipients such as rivers and lakes, increasing the terrestrial-to-aquatic transfer of permafrost carbon. Mobilization and subsurface transport mechanisms are poorly understood and not accounted for in global climate models, leading to high uncertainties in the predictions of the permafrost carbon feedback. Here, we focus on a hillslope in Endalen Valley, Svalbard, as a representative example of a high-Arctic hillslope underlain by continuous permafrost. We analyze solute transport in the form of a non-reactive tracer representing dissolved organic carbon (DOC) using a physics-based numerical model with the objective to study governing cryotic and hydrodynamic transport mechanisms relevant for warming permafrost regions. We first analyze transport times for DOC pools at different locations within the active layer under present-day climatic conditions and proceed to study susceptibility for deeper ancient carbon release in the upper permafrost due to thaw under different warming scenarios. Results suggest that DOC in the active layer near the permafrost table experiences rapid lateral transport upon thaw due to saturated conditions and lateral flow, while DOC close to the ground surface experiences slower transport due to flow in unsaturated soil. Deeper permafrost carbon release exhibits vastly different transport behaviors depending on warming and thaw rate. Gradual warming leads to small fractions of DOC being mobilized every year, while the majority moves vertically through percolation and cryosuction. Abrupt thaw resulting from a single very warm year leads to faster lateral transport times, similar to active layer DOC released in saturated conditions. Lastly, we analyze the potential susceptibility of DOC to mineralization to CO 2 prior to export due to soil moisture and temperature conditions. We find that high liquid saturation during transport coincides with very low mineralization rates and potentially inhibits mineralization into CO 2 before export. Overall, the results highlight the importance of subsurface hydrologic and thermal conditions for the retention and lateral export of permafrost carbon by subsurface flow.

Hamm, Alexandra [Stockholm Univ. (Sweden)] (ORCID:↗
NASA NTRS1989

Recursive forward dynamics for multiple robot arms moving a common task object

Recursive forward dynamics algorithms are developed and presented for an arbitrary number of robot arms moving a commonly held object. The multiarm forward dynamics problem is to find the angular accelerations at the joints and the contact forces that the arms impart to the task object. The problem also involves finding the acceleration of this object. The multiarm forward dynamics solutions provide a thorough physical and mathematical understanding of the way several arms behave in response to a set of applied joint moments. Such an understanding simplifies and guides the subsequent control design and experimentation process. The forward dynamics algorithms also provide the necessary analytical foundation for conducting analysis and simulation studies. The multiarm algorithms are based on the filtering and smoothing approach recently advanced for single-arm dynamics, and they can be built up modularly from the single-arm algorithms. The algorithms compute recursively the joint angle accelerations, the contact forces, and the task-object accelerations. Algorithms are also developed to evaluate in closed form the linear transformations from the active joint moments to the joint angle accelerations, to the task object accelerations, and to the task-object contact forces. A possible computing architecture is presented as a precursor to a more complete investigation of the computational performance of the dynamics algorithms.

Rodriguez, Guillermo↗
DOE OSTI2025

Streptococcus pneumoniae HtrA is a dynamic and monomeric virulence factor capable of forming larger oligomeric complexes

Abstract High‐temperature requirement A (HtrA) proteases are a conserved family of serine proteases central to protein quality control and bacterial virulence. While Gram‐negative and human HtrAs are structurally well studied, Gram‐positive homologs remain essentially uncharacterized. Here, we present the first integrated structural and mechanistic analysis of a Gram‐positive HtrA, from Streptococcus pneumoniae , a virulence factor essential for adhesion and infection in vivo. Proteomic profiling of an htrA knockout and cleavage assays demonstrate that S. pneumoniae HtrA is required for protein quality control, with the PDZ domain mediating substrate recognition. Biochemically, S. pneumoniae HtrA exists exclusively as a monomer in solution, a striking divergence from canonical trimeric HtrAs that we show is shared with other Gram‐positive homologs. NMR analyses reveal that the monomer dynamically samples open and closed conformations, while cryo‐EM of a catalytic mutant identifies a hexamer stabilized by a unique LoopA–PDZ interaction. Together, these findings define S. pneumoniae HtrA as a dynamic monomer with interdomain coupling between its protease and PDZ domains, establishing Gram‐positive HtrAs as a mechanistically divergent subgroup within the HtrA family.

Lee, Eunjeong [Department of Biochemistry and Mole↗
DOE OSTI2020

On the use of charged particles to characterize precipitation in irradiated reactor pressure vessel steels with a wide range of compositions

Nuclear reactor lifetimes may be limited by nano-scale Cu-Mn-Ni-Si precipitates (CRPs and MNSPs) that form under neutron irradiation (NI) of pressure vessel (RPV) steels, resulting in hardening and ductile to brittle transition temperature increases (embrittlement). Physical models of embrittlement must be based on characterization of precipitation as a function of the combination of metallurgical and irradiation variables. Here we focus on rapid and convenient charged particle irradiations (CPI) to both: a) compare to precipitates formed in NI; and, b) use CPI to efficiently explore precipitation in steels with a very wide range of compositions. Atom probe tomography (APT) comparisons show NI and CPI for similar bulk steel solute contents yield nearly the same precipitate compositions, albeit with some differences in their number density, size and volume fraction (f) dose (dpa) dependence. Yet, the overall precipitate evolutions are very similar. Advanced high Ni (>3 wt%) RPV steels, with superior unirradiated properties, were also investigated at high CPI dpa. For typical Mn contents, MNSPs have Ni 16 Mn 6 Si 7 or Ni 3 Mn 2 Si phase type compositions, with f values that are close to the equilibrium phase separated values. In contrast, in steels with very low Mn and high Ni, Ni 2-3 Si silicide phase type precipitate compositions are observed; and when Ni is low, the precipitate compositions are close to the MnSi phase field. Low Mn significantly reduces, but does not eliminate, precipitation in high Ni steels. A comparison of dispersed barrier model predictions with measured hardening data suggests that the Ni-Si dominated precipitates are weaker dislocation obstacles than the G phase type MNSPs.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
NASA NTRS1988

Recursive forward dynamics for multiple robot arms moving a common task object

Recursive forward dynamics algorithms are developed for an arbitrary number of robot arms moving a commonly held object. The multiarm forward dynamics problem is to find the angular accelerations at the joints and the contact forces that the arms impart to the task object. The problem also involves finding the acceleration of this object. The multiarm forward dynamics solutions provide a thorough physical and mathematical understanding of the way several arms behave in response to a set of applied joint moments. Such an understanding simplifies and guides the subsequent control design and experimentation process. The forward dynamics algorithms also provide the necessary analytical foundation for conducting analysis and simulation studies. The multiarm algorithms are based on the filtering and smoothing approach recently advanced for single-arm dynamics, and they can be built up modularly from the single-arm algorithms. The algorithms compute recursively the joint-angle accelerations, the contact forces, and the task-object accelerations. Algorithms are also developed to evaluate in closed form the linear transformations from the active joint moments to the joint-angle accelerations, to the task-object accelerations., and to the task-object contact forces. A possible computing architecture is presented as a precursor to a more complete investigation of the computational performance of the dynamics algorithms.

Rodriguez, G.↗
NASA NTRS1974

A finite element formulation for supersonic flows around complex configurations

The problem of small perturbation potential supersonic flow around complex configurations is considered. This problem requires the solution of an integral equation relating the values of the potential on the surface of the body to the values of the normal derivative, which is known from the small perturbation boundary conditions. The surface of the body is divided into small (hyperboloidal quadrilateral) surface elements which are described in terms of the Cartesian components of the four corner points. The values of the potential (and its normal derivative) within each element are assumed to be constant and equal to its value at the centroid of the element. This yields a set of linear algebraic equations whose coefficients are given by source and doublet integrals over the surface elements. Closed form evaluations of the integrals are presented.

Morino, L.↗
NASA NTRS1992

Characteristics of obliquely propagating, nonlinear Alfven waves

Consideration is given to properties of obliquely propagating (with respect to the large-scale magnetic field) nonlinear Alfven and fast magnetosonic waves as described by the derivative nonlinear Schroedinger equation. The focus is on those waves which in the parallel-propagating limit are circularly polarized, exact solutions to the MHD equations. Results are presented of a study of the changes in the structure of these waves as the propagation angle is increased. Attention is also given to the evolution of initial conditions which do not correspond to stationary wave propagations. It is found that, as the propagation angle with respect to the mean field increases, the functional forms of the wave magnetic field components in the two directions orthogonal to the propagation direction become increasingly dissimilar. The angle about the mean field within which stationary waves resemble parallel-propagating, circularly polarized waves becomes smaller as the wavelength increases. If the 'initial condition' wave train is close to parallel-propagating, then there is little subsequent evolution in its form.

Spangler, Steven R.↗
NASA NTRS1973

A Finite-Element Formulation for Subsonic Flows Around Complex Configurations

The problem of potential steady subsonic flow around complex configurations is considered. This problem requires the solution of an integral equation relating the values of the potential on the surface of the body to the values of the normal derivative, which is known from the boundary conditions. The surface of the body is divided into small (hyperboloidal quadrilateral) surface elements, which are described in terms of the Cartesian components of the four corner points. The values of the potential (and its normal derivative) within each element is assumed to be constant and equal to its value at the centroid of the element. This yields a set of linear algebraic equations. The coefficients of the equation are given by source and doublet integrals over the surface elements. Closed form evaluations of the integrals are presented.

Morino, L.↗
NASA NTRS2017

Unconventional Imaging with Contained Granular Media

Typically, the cost of a space-borne imaging system is driven by the size and mass of the primary aperture. The solution that we propose uses a method to construct an imaging system in space in which the nonlinear optical properties of a cloud of micron-sized particles, shaped into a specific surface by electromagnetic means, and allows one to form a very large and lightweight aperture of an optical system, hence reducing overall mass and cost. Recent work at JPL has investigated the feasibility of a granular imaging system, concluding that such a system could be built and controlled in orbit. We conducted experiments and simulation of the optical response of a granular lens. In all cases, the optical response, measured by the Modulation Transfer Function, of hexagonal reflectors was closely comparable to that of a conventional spherical mirror. We conducted some further analyses by evaluating the sensitivity to fill factor and grain shape, and found a marked sensitivity to fill factor but no sensitivity to grain shape. We have also found that at fill factors as low as 30%, the reflection from a granular lens is still excellent. Furthermore, we replaced the monolithic primary mirror in an existing integrated model of an optical system (WFIRST Coronagraph) with a granular lens, and found that the granular lens that can be useful for exoplanet detection provides excellent contrast levels. We will present our testbed and simulation results in this paper.

Sidick, Erkin↗
NASA NTRS1987

An evaluation of higher-order model methods for calculating transient structural response

A higher-order modal method proposed by Leung for transient structural analysis entitled the force-derivative method is evaluated. This method repeatedly integrates by parts with respect to time the convolution-integral form of the structural response to produce successively better approximations to the contribution of the higher modes which are neglected in the modal summation. Comparisons are made of the force-derivative, the mode-displacement, and the mode-acceleration methods for several numerical example problems for various times, levels of damping, and forcing functions. The example problems include a tip-loaded cantilevered beam and a simply-supported multispan beam. The force-derivative method is shown to converge to an accurate solution in fewer modes than either the mode-displacement or the mode-acceleration methods. In addition, for problems in which there are a large number of closely-spaced frequencies whose mode shapes have a negligible contribution to the response, the force-derivative method is very effective in representing the effect of the important, but otherwise neglected, higher modes.

Camarda, Charles J.↗