Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Generalized 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 1,279 records · Page 71
All sourcesNASA NTRSDOE OSTI
NASA NTRSDate not supplied

Fracture Prediction of Epoxy Resin using Morse Bond Potential Embedded in GAFF using Molecular Dynamics Simulations

A generic method is developed to investigate the fracture behavior of an epoxy resin using atomistic molecular (MD) dynamics simulations. The epoxy system consisted of the stoichiometric mixture of tetra-/tri-functional epoxies and di-functional hardener molecules and was cured with a step-growth crosslinking algorithm. A hybrid force field (hFF) in which Morse bond function were added in second generation-general Amber Force Field (GAFF2) was used to capture bond breaking of the crosslinked epoxy system under uniaxial deformation. The Morse bond parameters for covalent bonds in the crosslinked backbone were fit to dissociation curves computed by CASPT2/6-311+G**. After the systems in various crosslink densities were obtained and equilibrated, the full length of stress-strain (𝜎-𝜀) curve was produced including initial elastic regime, yielding, plastic flow, strain hardening, and progressive failure. We found that the hFF was effective in studying the dissociation of crosslinked polymers and understanding the failure mechanism. In the results, the hFF delivered the same elastic property obtained from the unmodified force field GAFF. The stress-strain curve, however, began to be distinct in the plastic regime due to broken covalent bonds capable of altering the tensile behavior. Once the stress reached maximum, it progressively decreased until the ultimate failure (𝜎=0). Since hFF employed the existing GAFF functions with including a minimal change of its bond potential, this approach is easier to perform and computationally more efficient for unveiling the fracture behavior of polymer materials at the molecular level than other approaches with reactive force fields. We expect this approach will be utilized to accelerate the material-by-design process for thermosets by incorporate data from molecular models.

Changwoon Jang↗
DOE OSTI2023

Investigation of local distortion effects on X-ray absorption of ferroelectric perovskites from first principles simulations

Understanding the role of ferroelectric polarization in modulating the electronic and structural properties of crystals is critical for advancing these materials for overcoming various technological and scientific challenges. However, due to difficulties in performing experimental methods with the required resolution, or in interpreting the results of methods therein, the nanoscale morphology and response of these surfaces to external electric fields has not been properly elaborated. Here, in this work, we investigate the effect of ferroelectric polarization and local distortions in a BaTiO 3 perovskite, using two widely used computational approaches which treat the many-body nature of X-ray excitations using different philosophies, namely the many-body, delta-self-consistent-field determinant (mb-ΔSCF) and the Bethe–Salpeter equation (BSE) approaches. We show that in agreement with our experiments, both approaches consistently predict higher excitations of the main peak in the O–K edge for the surface with upward polarization. However, the mb-ΔSCF approach mostly fails to capture the L 2,3 separations at the Ti–L edge, due to the absence of spin–orbit coupling in Kohn–Sham density functional theory (KS-DFT) at the generalized gradient approximation level. On the other hand, and most promising, we show that application of the GW/BSE approach successfully reproduces the experimental XAS, both the relative peak intensities as well as the L 2,3 separations at the Ti–L edges upon ferroelectric switching. Thus simulated XAS is shown to be a powerful method for capturing the nanoscale structure of complex materials, and we underscore the need for many-body perturbation approaches, with explicit consideration of core-hole and multiplet effects, for capturing the essential physics in these systems.

77 NANOSCIENCE AND NANOTECHNOLOGY↗
DOE OSTI2021

Single-Site Metal Organic Complexes on Oxide Supports for Selective Alkane Functionalization

High levels of reaction selectivity for heterogeneous catalysis are generally difficult to achieve with traditional metal nanoparticle catalysts, due to the variety of metal binding sites available. In this project, we have developed a metal-ligand coordination strategy to form transition metal single-sites on high surface area supports as heterogeneous single-site catalysts. Metal centers can be atomically dispersed in the binding pockets of organic ligands on oxide supports. The systems have to be carefully designed so that tendencies for metal cluster formation or binding to surface defect sites are overcome by the attractive coordination environment provided by the organic ligand. This was confirmed by a comprehensive set of characterization methods, including XAS, XPS, XRD, TEM, and CO adsorption. Design and tuning of metal-support and support-ligand interactions were crucial for high structural uniformity. We demonstrate that supported metal-ligand SSCs are effective, recyclable catalysts for alkene hydrosilylation reactions. Compared with commercial homogeneous catalysts, they exhibit improved yield and selectivity, less metal aggregation and side reactions, and stronger tolerance with functionalized substrates. These results give some preliminary indication that coordination with organic ligands can improve activity and selectivity in metal/oxide heterogeneous catalysts.

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

Bending as a control knob for the electronic and optical properties of phosphorene nanoribbons

We have assessed mechanical bending as a powerful controlling tool for the electronic structure and optical properties of phosphorene nanoribbons. We use state-of-the-art density functional approximations in our work. The overall performance of the recently developed meta-generalized gradient approximation (meta-GGA) mTASK [Phys. Rev. Materials 5, 063803 (2021)] functional establishes the method as a useful alternative to the screened hybrid HSE06 for better band gaps of phosphorene nanoribbons. We present a detailed analysis to interpret the optical absorption of bent phosphorene nanoribbons using the GW-Bethe-Salpeter approximation. Finally, we demonstrate the important role of the unoccupied in-gap state and conclude that this in-gap state in armchair nanoribbons introduced by bending can significantly affect the properties of low-energy excitons and add some useful opportunities for applications in optoelectronic devices.

36 MATERIALS SCIENCE↗
NASA NTRS1991

Deformation and stress response of composite laminated shells under internal pressure

This paper presents a theoretical study of the response of filament wound composite shells under internal pressure. Each layer of the material is generally cylindrically anisotropic. By using cylindrically anisotropic elasticity field equations and Lekhnitskii's stress functions, a system of sixth-order ordinary differential equations is obtained. The general expressions for the stresses and displacements in the laminated composite shells under internal pressure are discussed. Two composite systems, graphite/epoxy and glass/epoxy, are selected to demonstrate the influence of degree of material anisotropy and fiber orientations on the axial and induced twisting deformation. Stress distributions of (45/-45)s symmetric angle-ply fiber-reinforced laminated shells are shown to illustrate the effect of radius-to-thickness ratio.

Yuan, F. G.↗
DOE OSTI2020

Differential credibility assessment for statistical downscaling

Climate science is increasingly using (i) ensembles of climate projections from multiple models derived using different assumptions and/or scenarios and (ii) process-oriented diagnostics of model fidelity. Efforts to assign differential credibility to projections and/or models are also rapidly advancing. A framework to quantify and depict the credibility of statistically downscaled model output is presented and demonstrated. Here, the approach employs transfer functions in the form of robust and resilient generalized linear models applied to downscale daily minimum and maximum temperature anomalies at 10 locations using predictors drawn from ERA-Interim reanalysis and two global climate models (GCM; GFDL-ESM2M and MPI-ESM-LR). The downscaled time series are used to derive several impact relevant CLIMDEX temperature indices that are assigned credibility based on (1) the reproduction of relevant large-scale predictors by the GCMs (i.e. fraction of regression beta-weights derived from predictors that are well-reproduced) and (2) the degree of variance in the observations reproduced in the downscaled series following application of a new variance inflation technique. Credibility of the downscaled predictands varies across locations, between the two GCM and is generally higher for minimum temperature than maximum temperature. The differential credibility assessment framework demonstrated here is easy to use and flexible. It can be applied as is to inform decision makers regarding projection confidence, and/or extended to include other components of the transfer functions, and/or used to weight members of a statistically downscaled ensemble.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Time Evolution of Correlation Functions in Quantum Many-Body Systems

We give rigorous analytical results on the temporal behavior of two-point correlation functions—also known as dynamical response functions or Green’s functions—in closed many-body quantum systems. We show that in a large class of translation-invariant models the correlation functions factorize at late times $\langle$A ( t ) B$\rangle$ β → $\langle$ A$\rangle$ β $\langle$ B$\rangle$ β , thus proving that dissipation emerges out of the unitary dynamics of the system. We also show that for systems with a generic spectrum the fluctuations around this late-time value are bounded by the purity of the thermal ensemble, which generally decays exponentially with system size. For autocorrelation functions we provide an upper bound on the timescale at which they reach the factorized late time value. Remarkably, this bound is only a function of local expectation values and does not increase with system size. We give numerical examples that show that this bound is a good estimate in nonintegrable models, and argue that the timescale that appears can be understood in terms of an emergent fluctuation-dissipation theorem. Our study extends to further classes of two point functions such as the symmetrized ones and the Kubo function that appears in linear response theory, for which we give analogous results.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS1981

Demonstration Advanced Avionics System (DAAS). Phase 1 report

An integrated avionics system which provides expanded functional capabilities that significantly enhance the utility and safety of general aviation at a cost commensurate with the general aviation market is discussed. Displays and control were designed so that the pilot can use the system after minimum training. Functional and hardware descriptions, operational evaluation and failure modes effects analysis are included.

Source record↗
DOE OSTI2024

Exploring interlayer coupling in the twisted bilayer PtTe 2

We have investigated interlayer interactions in the bilayer PtTe 2 system, which influence the electronic energy bands near the Fermi level. Our diffusion Monte Carlo (DMC) calculations for the high-symmetry bilayer stackings (AA, AB, AC) manifest distinct interlayer binding characteristics among the stacking modes by revealing significantly different interlayer separations depending on the stacking, which is critical to understanding the interlayer coupling of the twisted bilayers consisting of various local stacking arrangements. Furthermore, a comparison between the interlayer separations obtained from DMC and density functional theory (DFT) shows that meta-generalized gradient approximation (GGA)-based van der Waals–DFT results agree with DMC for different layer stackings, including twisted bilayers, but only the ground-state AA stacking matches well with GGA-based DFT predictions. This underscores the importance of accurate exchange-correlation potentials even for capturing the stacking-dependent interlayer binding properties. We further show that the variability in DFT-predicted interlayer separations is responsible for the large discrepancy of band structures in the 21 . 79 ∘ twisted bilayer PtTe 2 , affecting its classification as metallic or insulating. These results demonstrate the importance of obtaining a correct description of stacking-dependent interlayer coupling in modeling delicate bilayer systems at finite twists. Published by the American Physical Society 2024

36 MATERIALS SCIENCE↗
DOE OSTI2021

Short-range correlation physics at low renormalization group resolution

Recent experiments have succeeded in isolating processes for which short-range correlation (SRC) physics is dominant and well accounted for by SRC phenomenology. But an alternative and compelling picture emerges from renormalization group (RG) evolution to low RG resolution. At high RG resolution, SRCs are identified as components in the nuclear wave function with relative pair momenta greater than the Fermi momentum. Scale separation results in wave-function factorization that can be exploited with phenomenologies such as the generalized contact formalism or the low-order correlation operator approximation. Evolution to lower resolution shifts SRC physics from nuclear structure to the reaction operators without changing the measured observables. We show how the features of SRC phenomenology manifested at high RG resolution are cleanly identified at low RG resolution using simple two-body operators and local-density approximations with uncorrelated wave functions, all of which can be systematically generalized. We verify that the experimental consequences to date follow directly at low resolution from well-established properties of nucleon-nucleon interactions such as the tensor force. Thus the RG reconciles the contrasting pictures of the same experiment and shows how to get correct results using wave functions without SRC components. Finally, our demonstration has implications for the analysis of knockout reactions for which SRC physics is not cleanly isolated.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
NASA NTRS1976

Thermodynamics and statistical mechanics

The basic thermodynamic properties of gases are reviewed and the relations between them are derived from the first and second laws. The elements of statistical mechanics are then formulated and the partition function is derived. The classical form of the partition function is used to obtain the Maxwell-Boltzmann distribution of kinetic energies in the gas phase and the equipartition of energy theorem is given in its most general form. The thermodynamic properties are all derived as functions of the partition function. Quantum statistics are reviewed briefly and the differences between the Boltzmann distribution function for classical particles and the Fermi-Dirac and Bose-Einstein distributions for quantum particles are discussed.

Source record↗
NASA NTRS1920

The general theory of blade screws including propellers, fans, helicopter screws, helicoidal pumps, turbo-motors, and different kinds of helicoidal blades

Report presents a theory which gives a complete picture and an exact quantitative analysis of the whole phenomenon of the working of blade screws, but also unites in a continuous whole the entire scale of states of work conceivable for a blade screw. Chapter 1 is devoted to the establishment of the system of fundamental equations relating to the blade screw. Chapter 2 contains the general discussion of the 16 states of work which may establish themselves for a blade screw. The existence of the vortex ring state and the whirling phenomenon are established. All the fundamental functions which enter the blade-screw theory are submitted to a general analytical discussion. The general outline of the curve of the specific function is examined. Two limited cases of the work of the screw, the screw with a zero constructive pitch and the screw with an infinite constructive pitch, are pointed out. Chapter 3 is devoted to the study of the propulsive screw or propeller. (author)

De Bothezat, George↗
NASA NTRS1994

F2 peak electron density at Millstone Hill and Hobart: Comparsion of theory and measurement at solar maximum

This paper compares the observed behavior of the (F2) layer of the ionosphere at Millstone Hill and Hobart with calculations from the field line interhemispheric plasma (FLIP) model for solar maximum, solstice conditions in 1990. During the study period the daily F(sub 10.7) index varied by more than a factor of 2 (123 to 280), but the 81-day mean F(sub 10.7) (F(sub 10.7 A)) was almost constant near 190. Calculations were performed with and without the effects of vibrationally excited N2 (N(sup *)(sub 2) which affects the loss rate of atomic oxygen ions. In the case without N(sup *)(sub 2) there is generally good agreement between the model and measurement for the daytime, peak density of the F region (NmF2). Both the model and the measurement show a strong seasonal anomaly with the winter noon densities a factor of 3 to 4 greater than the summer noon densities at Millstone Hill and a factor of 2 greater at Hobart. The seasonal anomaly in the model is caused by changes in the neutral composition as given by the mass spectrometer and incoherent scatter (MSIS) 86 neutral density model. There is generally little or no increase in the observed noon NmF2 as a function of daily F(sub 10.7) except at Millstone Hill in winter. In contrast to the generally good agreement between model and data at noon, the model badly underestimates the density at night at Millstone Hill at all seasons. At Hobart the model reproduces the nighttime density variations well in both winter and summer. The international reference ionosphere (IRI) model generally provides a good representation of the average behavior of noon NmF2 and hmF2 but because the data show a lot of day-to-day variability, there are often large differences. The FLIP model is able to reproduce this variability when hmF2 is specified. The IRI model peak densities are better than the FLIP densities at night, but the IRI model does not represent the Millstone Hill summer data very well at night in 1990.

Richards, P. G.↗
NASA NTRS1996

Variational Methods in Sensitivity Analysis and Optimization for Aerodynamic Applications

Variational methods (VM) sensitivity analysis, which is the continuous alternative to the discrete sensitivity analysis, is employed to derive the costate (adjoint) equations, the transversality conditions, and the functional sensitivity derivatives. In the derivation of the sensitivity equations, the variational methods use the generalized calculus of variations, in which the variable boundary is considered as the design function. The converged solution of the state equations together with the converged solution of the costate equations are integrated along the domain boundary to uniquely determine the functional sensitivity derivatives with respect to the design function. The determination of the sensitivity derivatives of the performance index or functional entails the coupled solutions of the state and costate equations. As the stable and converged numerical solution of the costate equations with their boundary conditions are a priori unknown, numerical stability analysis is performed on both the state and costate equations. Thereafter, based on the amplification factors obtained by solving the generalized eigenvalue equations, the stability behavior of the costate equations is discussed and compared with the state (Euler) equations. The stability analysis of the costate equations suggests that the converged and stable solution of the costate equation is possible only if the computational domain of the costate equations is transformed to take into account the reverse flow nature of the costate equations. The application of the variational methods to aerodynamic shape optimization problems is demonstrated for internal flow problems at supersonic Mach number range. The study shows, that while maintaining the accuracy of the functional sensitivity derivatives within the reasonable range for engineering prediction purposes, the variational methods show a substantial gain in computational efficiency, i.e., computer time and memory, when compared with the finite difference sensitivity analysis.

Ibrahim, A. H.↗
DOE OSTI2020

Real-time coupled-cluster approach for the cumulant Green's function

Green’s function methods within many-body perturbation theory provide a general framework for treating electronic correlations in excited states. Here we investigate the cumulant form of the one-electron Green’s function based on the coupled-cluster equation of motion approach in an extension of our previous study. The approach yields a non-perturbative expression for the cumulant in terms of the solution to a set of coupled first order, non-linear differential equations. The method thereby adds non-linear corrections to traditional cumulant methods linear in the self energy. The approach is applied to the core-hole Green’s function and illustrated for a number of small molecular systems. For these systems we find that the non-linear contributions lead to significant improvements both for quasiparticle properties such as core-level binding energies, as well as the satellites corresponding to inelastic losses observed in photoemission spectra.

Vila, Fernando D.↗
DOE OSTI2024

High-Pressure Rate Rules for Ether Alkylperoxy Radical Isomerization

The first isomerization reaction of an alkylperoxy (RO 2 ) radical holds significant importance in low-temperature oxidation, as it governs the branching ratios of the hydroperoxyalkyl (QOOH) radicals, which influence the competition between the chain-propagation and chain-branching reactions. In this study, we systematically calculated high-pressure rate rules for the RO 2 isomerization reaction of monoethers, exploring 5-, 6-, 7-, and 8-membered ring transition states. Primary, secondary, and tertiary carbon sites, where both the abstracting peroxy group and the abstracted hydrogen are located, were considered, with particular emphasis on distinguishing between secondary carbons adjacent (alpha) and nonadjacent to the ether functional group. Using the G4//B3LYP/6-311++G(2df,2pd) level of theory and the transition state theory, we estimated the rate constants and the Arrhenius coefficient for over 120 possible isomerization reactions. We examined the effect of ring size and ring atoms, revealing that 6- and 7-membered ring isomerizations were generally the fastest. The impact of the ether functional group on transition states was investigated by comparing reactions with identical ring size, peroxy, and radical positions, but with the ether functional group positioned either outside (i.e., out) or inside (i.e., in) the transition state ring, leading to differences in the rate constants. When comparing to analogous alkane rate constants, differences of up to an order of magnitude were observed, underscoring the need for caution when assigning rate rules by analogy. We applied our rate constants in the di-iso-butyl ether kinetic model and evaluated their influence on low-temperature chemistry finding that they altered the branching ratios by up to a factor of 9, highlighting the significance of site-specific rate constants for more accurate low-temperature modeling.

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

Convergence of Newton's method for a single real equation

Newton's method for finding the zeroes of a single real function is investigated in some detail. Convergence is generally checked using the Contraction Mapping Theorem which yields sufficient but not necessary conditions for convergence of the general single point iteration method. The resulting convergence intervals are frequently considerably smaller than actual convergence zones. For a specific single point iteration method, such as Newton's method, better estimates of regions of convergence should be possible. A technique is described which, under certain conditions (frequently satisfied by well behaved functions) gives much larger zones where convergence is guaranteed.

Campbell, C. W.↗
DOE OSTI2024

HHL algorithm with mapping function and enhanced sampling for model predictive control in microgrids

Here, this paper presents a refined quantum Harrow Hassidim Lloyd (HHL) algorithm for microgrid control. The first novelty of the developed method is that a mapping shift function enables the original HHL algorithm to handle general linear equations with non-singular and indefinite matrix. Second, a method of Matrix Extension for Amplifying Sampling Probabilities of Intended Solution (ME-ASPI) is proposed to design the reformulated linear algebraic equations, allowing for improved sampling efficiency of the quantum tomography in the refined HHL algorithm. Then, we applied the method to solve the model predictive control (MPC) problem in nonlinear dynamical microgrids. Specifically, with the ME-ASPI method, the refined HHL algorithm can effectively obtain the intended partial optimal control inputs for MPC. The optimization of quadratic programming problem in each time step of MPC is transformed into a linear system problem, which is addressed by the proposed quantum solver through using only partial information, with the time complexity improved from $\mathscr{O}(\mathscr{N}^{2.37286})$ classically to $\mathscr{O}(\mathscr{N}^{2} log \mathscr{N}$ x $p$ log $p)$ in quantum. Numerical examples have validated the effectiveness of the refined HHL algorithm with the proposed mapping function and the ME-ASPI method. By leveraging quantum properties, the proposed method provides a hybrid quantum–classical framework for microgrid control. This generic method can also potentially tackle many other challenges in analyzing and controlling general complex engineered systems.

24 POWER TRANSMISSION AND DISTRIBUTION↗