Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Implicit functions”

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 109 records · Page 6
All sourcesNASA NTRSDOE OSTI
NASA NTRS1983

Derivation of the fundamental equation of sound generated by moving aerodynamic surfaces

Ffowcs Williams and Hawkings (1969) based their derivation of the fundamental equation of the sound generated by arbitrarily moving aerodynamic surfaces on the study of mass and momentum balance of a control volume imbedding a mathematical surface(s) exactly corresponding to real surface(s). These investigators also sketched an alternative method, employing generalized functions, for its derivation. This latter method, which was later developed by Farassat (1975), is purely mathematical and formal. Goldstein (1976) used the free-space Green function to produce an implicit derivation of the Ffowcs Williams and Hawkings equation. In the study presented here, Lowson's (1965) concept of moving point singularities is generalized to moving surface singularities, and a new derivation is given of the fundamental equation. The derivation is based on topological considerations of the underlying space, the fluid medium, and the integral properties of the Dirac delta function.

Aggarwal, H. R.↗
DOE OSTI2019

Implicit renormalization approach to the problem of Cooper instability

In the vast majority of cases, superconducting transition takes place at exponentially low temperature T c out of the Fermi liquid regime. We discuss the problem of determining T c from known system properties at temperatures T >> T c , and stress that this cannot be done reliably by following the standard protocol of solving for the largest eigenvalue of the original gap-function equation. However, within the implicit renormalization approach, the gap-function equation can be used to formulate an alternative eigenvalue problem, where solving leads to an accurate prediction for both T c and the gap function immediately below T c . Furthermore, with diagrammatic Monte Carlo techniques, this eigenvalue problem can be solved without invoking the matrix inversion or even explicitly calculating the four-point vertex function.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2026

Comparison between explicit and implicit discretization strategies for a dissipative thermal environment

We investigate strategies for simulating open quantum systems coupled to dissipative baths by comparing explicit wave function-based discretization [via multi-layer multi-configuration time-dependent Hartree (ML-MCTDH)] and the implicit density matrix-based master equation method [via tree tensor network hierarchical equations of motion (TTN-HEOM)]. For dissipative baths characterized by exponentially decaying bath correlation functions, the implicit discretization approach of HEOM—rooted in bath correlation function decompositions—proves significantly more efficient than explicit discretization of the bath into discrete harmonic modes. Explicit methods, like ML-MCTDH, require extensive mode discretization to approximate continuum baths, leading to computational bottlenecks. Case studies for two-level systems and a Fenna–Matthews–Olson complex model highlight TTN-HEOM’s superiority in capturing dissipative dynamics with relaxations with a minimal number of auxiliary modes, while the explicit methods are as exact as the HEOM in pure dephasing regimes. This comparison is enabled by the TENSO package, which has both ML-MCTDH and TTN-HEOM implemented using the same computational structure and propagation strategy.

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

Reaction Energetics and 13 C Fractionation of Alanine Transamination in the Aqueous and Gas Phases

The alanine transaminase enzyme catalyzes the transfer of an amino group from alanine to α-ketoglutarate to produce pyruvate and glutamate. Isotope fractionation factors (IFFs) for the reaction +H 3 NCH(CH 3 )COO – + – OOCCH 2 CH 2 C(O)COO – ↔ CH 3 C(O)COO – + + H 3 NCH(CH 2 CH 2 COO – )COO – (zwitterionic neutral alanine + doubly deprotonated α-ketoglutarate ↔ pyruvate + zwitterionic glutamate anion) were calculated from the partition functions of explicitly and implicitly solvated molecules at 298 K. Calculations were done for alanine (noncharge separated, zwitterion, deprotonated), pyruvic acid (neutral, deprotonated), glutamic acid (noncharge separated, zwitterion, deprotonated, doubly deprotonated), and α-ketoglutaric acid (neutral, deprotonated, doubly deprotonated). The computational results, calculated from gas phase- and aqueous-optimized clusters with explicit H 2 O molecules at the MP2/aug-cc-pVDZ and MP2/aug-cc-pVDZ/COSMO levels, respectively, predict that substitution of 13 C at the C2 position of alanine and pyruvic acid and their various forms leads to the C2 position of pyruvic acid/pyruvate being enriched in 13 C/ 12 C ratio by 9‰. Simpler approaches that estimate the IFFs based solely on changes in the zero-point energies (ZPEs) are consistent with the higher-level model. ZPE-based IFFs calculated for simple analogues formaldehyde and methylamine (analogous to the C2 positions of pyruvate and alanine, respectively) predict a 13 C enrichment in formaldehyde of 7–8‰ at the MP2/aug-cc-pVDZ and aug-cc-pVTZ levels. A simple predictive model using canonical functional group frequencies and reduced masses for 13 C exchange between R 2 C=O and R 2 CH–NH 2 predicted enrichment in R 2 C=O that is too large by a factor of two but is qualitatively accurate compared with the more sophisticated models. Our models are all in agreement with the expectation that pyruvate and formaldehyde will be preferentially enriched in 13 C because of the strength of their >C=O bond relative to that of ≡C–NH 2 in alanine and methylamine. 13 C/ 12 C substitution is also modeled at the methyl and carboxylic acid sites of alanine and pyruvic acid, respectively.

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

Physics-based adaptivity of a spectral method for the Vlasov–Poisson equations based on the asymmetrically-weighted Hermite expansion in velocity space

We propose a spectral method for the 1D-1V Vlasov–Poisson system where the discretization in velocity space is based on asymmetrically-weighted Hermite functions, dynamically adapted via a scaling α and shifting u of the velocity variable. Specifically, at each time instant an adaptivity criterion selects new values of α and u based on the numerical solution of the discrete Vlasov–Poisson system obtained at that time step. Once the new values of the Hermite parameters α and u are fixed, the Hermite expansion is updated and the discrete system is further evolved for the next time step. The procedure is applied iteratively over the desired temporal interval. The key aspects of the adaptive algorithm are: the map between approximation spaces associated with different values of the Hermite parameters that preserves total mass, momentum and energy; and the adaptivity criterion to update α and u based on physics considerations relating the Hermite parameters to the average velocity and temperature of each plasma species. For the discretization of the spatial coordinate, we rely on Fourier functions and use the implicit midpoint rule for time stepping. The resulting numerical method possesses intrinsically the property of fluid-kinetic coupling, where the low-order terms of the expansion are akin to the fluid moments of a macroscopic description of the plasma, while kinetic physics is retained by adding more spectral terms. Moreover, the scheme features conservation of total mass, momentum and energy associated in the discrete, for periodic boundary conditions. A set of numerical experiments confirms that the adaptive method outperforms the non-adaptive one in terms of accuracy and stability of the numerical solution.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Stability Analysis of Substituted Cobaltocenium [Bis(cyclopentadienyl)cobalt(III)] Employing Chemistry-Informed Neural Networks

Cobaltocenium derivatives are promising components of the anion exchange membranes due to their excellent thermal and alkaline stability under the operating conditions of a fuel cell. Here we present an efficient modeling approach of assessing the chemical stability of substituted cobaltocenium CoCp 2 + based on the computed electronic structure enhanced by machine learning techniques. Within the aqueous environment, the positive charge of the metal cation is balanced by the hydroxide anion through formation of the CoCp 2 + OH¯ complexes, whose dissociation is studied within the implicit solvent employing density functional theory. The data set of about 118 species based on 42 substituent groups characterized by a range of electron- donating (ED) and electron-withdrawing (EW) properties is constructed and analyzed. Given 12 carefully chosen chemistry-informed descriptors of the complexes and relevant fragments, the stability of the complexes is found to strongly correlate with the energies of the highest occupied and lowest unoccupied molecular orbitals, modulated by a switching function of the Hirshfeld charge. The latter is used as a measure of the electron withdrawing-donating character of the substituents. Based on this observation from the conventional regression analysis, two fully connected, feed-forward neural network (FNN) models with different unit structures, called the chemistry-informed (CINN) and the quadratic (QNN) neural networks, are developed. Both models predict the bond dissociation energies of the cobaltocenium complexes with mean relative errors less than 5.40% and average absolute errors less than 0.94 kcal/mol. The results show the potential of QNN to efficiently capture more complex relationships. Here, the concept of incorporating the domain (chemical) knowledge/insight into the neural network structure paves the way to applications of machine learning techniques with small data sets, ultimately leading to better predictive models compared to the conventional regression analysis.

36 MATERIALS SCIENCE↗
DOE OSTI2018

Understanding the MXene Pseudocapacitance

MXenes have attracted great attention as next-generation capacitive energy-storage materials, but the mechanisms underlying their pseudocapacitive behavior are not well understood. Here we provide a theoretical description of the surface redox process of Ti3C2Tx (T = O, OH), a prototypical MXene, in 1 M H2SO4 electrolyte, based on joint density functional theory with an implicit solvation model and the analysis of Gibbs free energy under a constant-electrode potential. From the dependence of the O/OH ratio (or the surface H coverage) and the surface charge on the applied potential, we obtain a clear picture of the capacitive energy-storage mechanism of Ti3C2Tx that shows good agreement with previous experimental findings in terms of the integral capacitance and Ti valence change. We find a voltage-dependent redox/double-layer co-charging behavior: the capacitive mechanism is dominated by the redox process, but the electric double-layer charge works against the redox process. This new insight may be useful in improving the capacitance of MXenes.

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

Aqueous-phase effects on ethanol decomposition over Ru-based catalysts

The effects of an aqueous phase on ethanol decomposition for hydrogen production over a Ru(0001) catalyst surface model have been investigated from first principles. Solvent effects on the reaction mechanism and kinetic parameters have been quantified with the help of a microkinetic reactor model, density functional theory, and an implicit solvation scheme (iSMS). Our calculations indicate that in both vapor- and aqueous-phase reaction environments, the ethanol decomposition starts with acetaldehyde formation on the surface, some of which further dehydrogenates to ketenyl species (CHCO), where the C–C bond cleaves to form methylidyne (CH) and CO. In the vapor phase, adsorbed CH gets hydrogenated to methane, and CO desorbs or undergoes methanation reducing the amount of hydrogen produced. In contrast, under aqueous phase reaction conditions, the methanation is inhibited, and the water–gas shift (WGS) reaction is accelerated, leading to complete conversion of CO to CO 2 and H 2 . Furthermore, calculations indicate that the observed reaction behavior under aqueous phase reforming conditions originates primarily from the higher water chemical potential, and implicit solvent models predict only a small solvation effect.

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

Novel functional and distributed approaches to data analysis available in ROOT

The bright future of particle physics at the Energy and Intensity frontiers poses exciting challenges to the scientific software community. The traditional strategies for processing and analyzing data are evolving in order to (i) offer higher-level programming model approaches and (ii) exploit parallelism to cope with the ever increasing complexity and size of the datasets. This contribution describes how the ROOT framework, a cornerstone of software stacks dedicated to particle physics, is preparing to provide adequate solutions for the analysis of large amount of scientific data on parallel architectures. The functional approach to parallel data analysis provided with the ROOT TDataFrame interface is then characterized. The design choices behind this new interface are described also comparing with other widely adopted tools such as Pandas and Apache Spark. The programming model is illustrated highlighting the reduction of boilerplate code, composability of the actions and data transformations as well as the capabilities of dealing with different data sources such as ROOT, JSON, CSV or databases. Details are given about how the functional approach allows transparent implicit parallelization of the chain of operations specified by the user. The progress done in the field of distributed analysis is examined. In particular, the power of the integration of ROOT with Apache Spark via the PyROOT interface is shown. In addition, the building blocks for the expression of parallelism in ROOT are briefly characterized together with the structural changes applied in the building and testing infrastructure which were necessary to put them in production.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2023

Variance Reduction within Implicit Monte Carlo Thermal Radiation Transport using the Local Importance Function Transform [Slides]

System of equations derived by Fleck and Cummings to model time-dependent thermal radiative transfer (TRT) problems: Solution of radiation specific intensity (I), non-linearly coupled to material internal energy (U); Absorption and re-emission approximated by effective scattering events. Useful for high-energy density physics simulations: i.e., astrophysics, inertial confinement fusion (ICF). Whereas deterministic methods are fully discretized and free of statistical noise, Monte Carlo methods allow for dynamic sampling of the phase space at the price of statistical noise; Monte Carlo methods are also characterized by a slow $\frac{1}{√Ν}$ convergence rate.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2023

Semi-Implicit Computation of Fast Modes in a Scheme Integrating Slow Modes by a Leapfrog Method Based on a Selective Implicit Time Filter

Abstract A scheme for integration of atmospheric equations containing terms with differing time scales is developed. The method employs a filtered leapfrog scheme utilizing a fourth-order implicit time filter with one function evaluation per time step to compute slow-propagating phenomena such as advection and rotation. The terms involving fast-propagating modes are handled implicitly with an unconditionally stable method that permits application of larger time steps and faster computations compared to fully explicit treatment. Implementation using explicit and recurrent formulation is provided. Stability analysis demonstrates that the method is conditionally stable for any combination of frequencies involved in the slow and fast terms as they approach the origin. The implicit filter used in the method damps the computational modes without noticeably sacrificing the accuracy of the physical mode. TheO[(Δt 4 )] accuracy for amplitude errors achieved by the implicitly filtered leapfrog is preserved in applications where terms responsible for fast propagation are integrated with a semi-implicit method. Detailed formulation of the method for soundproof nonhydrostatic anelastic equations is provided. Procedures for implementation in global spectral shallow-water models are also given. Examples comparing numerical and analytical solutions for linear gravity waves demonstrate the accuracy of the scheme. The performance is also shown in more practical nonlinear applications, where numerical solutions accomplished by the method are evaluated against those computed from a scheme where the slow terms are handled by the third-order Runge–Kutta scheme. It demonstrates that the method is able to accurately resolve fine-scale dynamics of Kelvin–Helmholtz shear instabilities, the evolution of density current, and nonlinear drifts of twin tropical cyclones.

Meteorology & Atmospheric Sciences↗
NASA NTRS1996

Explicit and Implicit Resources: A Simplified Approach to User Requirements Modeling

The International Space Station (ISS) program is dependent on the concept of distributed planning, which places a large responsibility in the hands of the science users. The users assume the responsibility in submitting the operations scheduling requirements in the format required by the scheduling software. The concept of explicit and implicit resources, developed to facilitate the operations scheduling of ISS by non-scheduling experts, is presented. The proposed concepts are defined and examples of scheduling requirements, implemented using explicit and implicit resources, are described together with the requirements modeling and scheduler formatting functions. The software developed to support the implicit and explicit resources concept is discussed.

Hagopian, Jeff↗
DOE OSTI2021

Confinement effects on the solvation structure of solvated alkaline metal cations in a single-digit 1T-MoS 2 nanochannel: A first-principles study

Confinement plays an important role in determining ion transport in porous materials, which, in turn, may influence the performance of many energy storage and desalination devices. In this work, we combined density functional theory (DFT) with an implicit solvation model and ab initio molecular dynamics (AIMD) to investigate the effects of nanoconfinement on several solvated alkaline metal cations in a single-digit 1T-MoS 2 nanochannel. Our DFT calculations with a solvation model indicated that cations with stronger hydration energy introduce a higher number of co-intercalated water molecules into the channel, consistent with early experimental observation obtained for MXene (2D transition metal carbide) channels. The predicted optimal water numbers for the cations were then used for AIMD simulations that explicitly include the effects of the solvent. When compared with the cations in bulk solution, our simulations showed that the hydration structure and coordination number (CN) of the solvated cations confined in the MoS 2 channel can be significantly altered. Further, we found that larger cations with weaker hydration energy (K + , Rb + , and Cs + ) exhibited a distinctive CN decrease under confinement, while smaller cations (Li + and Na + ) retained a similar hydration shell as in the bulk solution. More specifically, the hydration shell of large cations (K + , Rb + , and Cs + ) in MoS 2 showed similar features of the coordination angle to the bulk, which suggests the partially broken hydration shell with no geometry change under confinement. Our simulations provided insights into the change of the hydration structure of alkaline metal cations under confinement, which may have important implications on their transport in the 1T-MoS 2 channel.

36 MATERIALS SCIENCE↗
DOE OSTI2015

A regulatory perspective on the response function method for determining allowable content of a radioactive materials transportation package or storage cask - 15023

In order to simplify shielding design, there has been an increasingly use of a so-called dose rate response function method for determining the allowable contents in radioactive materials (including spent fuel) transportation package and spent fuel storage cask designs. In essence, this approach uses a pre-calculated dose rate at a given location outside the storage cask or transportation package for a given source inside the cask. It assumes that there is a one-to-one relationship between a source inside the cask and the dose rate outside the cask, i.e., a dose rate response function of a source. Thus, the allowable content is obtained by simply dividing the targeted dose rate by the response function. This method converts the complex particle transport problem into a simple arithmetic calculation. If used correctly and prudently, this approach can greatly simplify the shielding design of casks. However, the use of this method must demonstrate the validity of the assumption, that is, the nuclear physics property of the media through which the particle transverses must be similar to that used in computing the response function. Because of the assumptions and approximations, either implicitly or explicitly, made in the models for dose rate response functions, to assure safety the users should also consider the errors in the models and the bias and uncertainties associated with the code used to calculate the response function. This paper reviewed the advantages and some of the major potential pitfalls of this method from the regulatory perspective. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2015

New implicitly solvable potential produced by second order shape invariance

The procedure proposed recently by Bougie et al. (2010) to study the general form of shape invariant potentials in one-dimensional Supersymmetric Quantum Mechanics (SUSY QM) is generalized to the case of Higher Order SUSY QM with supercharges of second order in momentum. A new shape invariant potential is constructed by this method. It is singular at the origin, it grows at infinity, and its spectrum depends on the choice of connection conditions in the singular point. The corresponding Schrödinger equation is solved explicitly: the wave functions are constructed analytically, and the energy spectrum is defined implicitly via the transcendental equation which involves Confluent Hypergeometric functions. - Highlights: • New potential with 2nd order irreducible shape invariance was constructed. • The connection conditions at the singularity of potential were obtained. • The explicit expressions for all wave functions were derived. • The implicit equation for the energy spectrum was obtained.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

CoordNet: Data Generation and Visualization Generation for Time-Varying Volumes via a Coordinate-Based Neural Network

Although deep learning has demonstrated its capability in solving diverse scientific visualization problems, it still lacks generalization power across different tasks. To address this challenge, we propose CoordNet, a single coordinate-based framework that tackles various tasks relevant to time-varying volumetric data visualization without modifying the network architecture. The core idea of our approach is to decompose diverse task inputs and outputs into a unified representation (i.e., coordinates and values) and learn a function from coordinates to their corresponding values. We achieve this goal using a residual block-based implicit neural representation architecture with periodic activation functions. We evaluate CoordNet on data generation (i.e., temporal super-resolution and spatial super-resolution) and visualization generation (i.e., view synthesis and ambient occlusion prediction) tasks using time-varying volumetric data sets of various characteristics. Here, the experimental results indicate that CoordNet achieves better quantitative and qualitative results than the state-of-the-art approaches across all the evaluated tasks. Source code and pre-trained models are available at https://github.com/stevenhan1991/CoordNet.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2019

Analysis of the Spectral Stability of the Generalized Runge–Kutta Methods Applied to Initial-Boundary-Value Problems for Equations of the Parabolic Type. II. Implicit Methods

We consider specific realizations of different implicit generalized Runge–Kutta methods as applied to the numerical integration with respect to time of initial-boundary-value problems for the second-order parabolic equations and investigate their spectral stability. It is shown that all implicit generalized Runge–Kutta methods are unconditionally spectrally stable but some of them have the conditional property of monotonicity of the numerical solution with respect to time. The functions of spectral stability of the implicit generalized Runge–Kutta methods are rational. We compare the analytic solution of the nonstationary one-dimensional problem of heat conduction with the numerical solutions of this problem obtained by different implicit generalized Runge–Kutta methods. It is shown that, in this case, the application of the one-stage Radau methods with subsequent discretization of the problem with respect to the space variable leads to the classical forward finite difference scheme (Laasonen scheme), whereas the use of the one-stage Gauss–Legendre method leads to a six-point symmetric scheme (Crank–Nicolson scheme). It is shown that diagonally implicit generalized Nørsett and Burrage methods are realized in almost the same way as the one-stage Radau and Gauss–Legendre methods but their accuracy in the time step is 10–1000 times higher. On the basis of comparison of the numerical and analytic solutions, we conclude that, in order to get practically suitable numerical solutions without any restrictions on the time step, it is reasonable to use one- and three-stage generalized Radau methods or two- and four-stage Lobatto IIIC methods. All other explicit and implicit generalized Runge–Kutta methods require certain restrictions imposed on the time step.

42 ENGINEERING↗
DOE OSTI2015

Quantum Effects on the Capacitance of Graphene-Based Electrodes

We recently measured quantum capacitance for electric double layers (EDL) at electrolyte/graphene interfaces. However, the importance of quantum capacitance in realistic carbon electrodes is not clear. Toward understanding that from a theoretical perspective, here we studied the quantum capacitance and total capacitance of graphene electrodes as a function of the number of graphene layers. The quantum capacitance was obtained from electronic density functional theory based on fixed band approximation with an implicit solvation model, while the EDL capacitances were from classical density functional theory. We found that quantum capacitance plays a dominant role in total capacitance of the single-layer graphene both in aqueous and ionic-liquid electrolytes but the contribution decreases as the number of graphene layers increases. Moreover, the total integral capacitance roughly levels off and is dominated by the EDL capacitance beyond about four graphene layers. Finally, because many porous carbons have nanopores with stacked graphene layers at the surface, this research provides a good estimate of the effect of quantum capacitance on their electrochemical performance.

36 MATERIALS SCIENCE↗