A method for the numerical integration of coupled first-order differential equations with greatly different time constants.
Numerical integration of coupled first order ODE OF greatly differing time constant
DIFFERENTIAL EQUATION↗
SEARCH · Engineering Papers
Results for “coupling method”
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.
Hydrodynamic Coupling to a Homogenized Radiation Transport Method based on Young Measures
Resolving radiation transport fields subject to opacity profiles with strong, oscillatory line structure while potentially falling under intermediate optical depth conditions presents a numerical challenge in radiation transport modeling. The Young measure-based homogenization technique formulated by Haut et al. (2017) was investigated as a candidate method for resolving radiation fields under such conditions more accurately. The method was compared against frequently-utilized mean opacity methods as the Rosseland and Planck formulations. In this work, all methods were tested through radiation slab calculations separately comprised of aluminum, copper, and krypton, each for different thermodynamic conditions. Following these offline radiation slab calculations, demonstrations shifted towards the SCEPTRE radiation transport code and, subsequently, the multiphysics ALEGRA code for approximately-coupled radiation-material simulations. Throughout all the simulations shown in this study, for a fixed computational cost, the homogenized method was observed to be more accurate than any of the solutions determined through traditional mean opacity approaches.
73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
Quantum Solvers for Plane-Wave Hamiltonians: Abridging Virtual Spaces Through the Optimization of Pairwise Correlations
For many-body methods such as MCSCF and CASSCF, in which the number of one-electron orbitals is optimized and independent of the basis set used, there are no problems with using plane-wave basis sets. However, for methods currently used in quantum computing such as select configuration interaction (CI) and coupled cluster (CC) methods, it is necessary to have a virtual space that is able to capture a significant amount of electron-electron correlation in the system. The virtual orbitals in a pseudopotential plane-wave Hartree–Fock calculation, because of Coulomb repulsion, are often scattering states that interact very weakly with the filled orbitals. As a result, very little correlation energy is captured from them. The use of virtual spaces derived from the one-electron operators has also been tried, and while some correlations are captured, the amount is quite low. To overcome these limitations, we have been developing new classes of algorithms to define virtual spaces by optimizing orbitals from small pairwise CI Hamiltonians, which we term as correlation optimized virtual orbitals with the abbreviation COVOs. With these procedures, we have been able to derive virtual spaces, containing only a few orbitals, which are able to capture a significant amount of correlation. The focus in this manuscript is on using these derived basis sets to target full CI (FCI) quality results for H2 on near-term quantum computers. However, the initial results for this approach were promising. We were able to obtain good agreement with FCI/cc-pVTZ results for this system with just 4 virtual orbitals, using both FCI and quantum simulations. The quality of the results using COVOs suggests that it may be possible to use them in other many-body approaches, including coupled cluster and Møller–Plesset perturbation theories, and open up the door to many-body calculations for pseudopotential plane-wave basis set methods.
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
Energetics and mechanisms for the acetonyl radical + O 2 reaction: An important system for atmospheric and combustion chemistry
The acetonyl radical (•CH 2 COCH 3 ) is relevant to atmospheric and combustion chemistry due to its prevalence in many important reaction mechanisms. One such reaction mechanism is the decomposition of Criegee intermediates in the atmosphere that can produce acetonyl radical and OH. In order to understand the fate of the acetonyl radical in these environments and to create more accurate kinetics models, we have examined the reaction system of the acetonyl radical with O 2 using highly reliable theoretical methods. Structures were optimized using coupled cluster theory with singles, doubles, and perturbative triples [CCSD(T)] with an atomic natural orbital (ANO0) basis set. Energetics were computed to chemical accuracy using the focal point approach involving perturbative treatment of quadruple excitations [CCSDT(Q)] and basis sets as large as cc-pV5Z. The addition of O 2 to the acetonyl radical produces the acetonylperoxy radical, and multireference computations on this reaction suggest it to be barrierless. Additionally, no submerged pathways were found for the unimolecular isomerization of the acetonylperoxy radical. Besides dissociation to reactants, the lowest energy pathway available for the acetonylperoxy radical is a 1-5 H shift from the methyl group to the peroxy group through a transition state that is 3.3 kcal mol -1 higher in energy than acetonyl radical + O 2 . The ultimate products from this pathway are the enol tautomer of the acetonyl radical along with O 2 . Multiple pathways that lead to OH formation are considered; however, all of these pathways are predicted to be energetically inaccessible, except at high temperatures.
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
Six degrees of freedom compliant flexure and autonomous assembly using same
A flexure having six degrees of freedom is disclosed. The constraint and compliance of the flexure in all six degrees of freedom may be independently tuned for an application. The flexure may include various combinations of sections including, for example, straight arm and curvilinear sections. The constraint and compliance of the flexure is determined by the geometry, cross-sectional area, cross-sectional shape, and material used to form the various sections. A coupling mechanism and a corresponding method for using the coupling mechanism for automated or autonomous assembly of elements using tuned flexures is disclosed. The automated or autonomous assembly of elements employs a coupling mechanism having a motor driven grip and a coupler, in which the coupler includes a six degrees of freedom flexure. The motor driven grip and coupler may optionally include electrical or optical interconnections and self-aligning features.
Gillund, Daniel Paul↗
Dynamic Approach to Dependency Analysis in Human Reliability Analysis: Application in a Stream Generator Tube Rupture Scenario
Dependency analysis in human reliability analysis (HRA) is a method of adjusting the failure probability of a given action by considering the impact of the action preceding it. It plays a role in reasonably accounting for human actions in the context of probabilistic safety assessments (PSAs), preventing PSA results from being estimated too optimistically based on the HRA results. Nevertheless, the existing dependency methods present a couple of challenges in that the quantification approaches rarely explain the adjustment of human error probabilities (HEPs). For this reason, the authors’ previous research has pointed out challenges of the existing dependency approaches and conceptually, theoretically proposed a performance shaping factor (PSF)-based dynamic dependency analysis method for HRA in order to complement the existing dependency methods. The current paper explores the latest version of the method and guidance for applying it to a steam generator tube rupture (SGTR) scenario.
99 GENERAL AND MISCELLANEOUS↗
The Heat of Formation of HNO
The HNO molecule is of interest in both combustion and atmospheric chemistry. For example, Guadagnini et al. have recently presented ab initio potential energy surfaces for the three lowest lying electronic states of HNO and then used these in examining several chemical reactions that take place in the combustion of nitrogen containing fuels and in the oxidation of atmospheric nitrogen. We have previously studied the ground state potential energy surface (i.e., stationary points along the HNO rev. reaction HON path), vibrational spectrum (using an accurate quartic force field), zero-point energy, and bonding of HNO using coupled-cluster ab initio methods. HNO is also very interesting because of the unique nature of its bonding characteristics. That is, the potential energy surface is very flat along the H-N bonding coordinate thereby giving unusual harmonic and fundamental vibrational frequencies, and the H-N bond energy is rather weak in comparison to other H-N bond energies. In fact, using experimental heats of formation for HO, H, and NO, the H- bond energy is computed to be only 49.9 kcal/ mol (298 K). However, ab initio calculations of isodesmic reaction energies involving HNO, FNO, ClNO, and several other molecules have shown that there is an inconsistency in the experimental heats of formation of the XNO (X=H, F, and Cl) species. Hence the motivation for this study was to determine a very accurate(DELTA)H(sup o)(sub f) value for HNO using state of-the-art ab initio methods. Based on many recent studies it is evident that the singles and doubles coupled-cluster method that includes a perturbational estimate of the effects of connected triple excitations, denoted CCSD(T), in conjunction with large one-particle basis sets should be reliable to better than +0.8 kcal/mol for this quantity. The computational methodology is described in the next section followed by our results and discussion. Conclusions are presented in the final section.
Lee, Timothy J.↗
The Heat of Formation of HNO
The HNO molecule is of interest in both combustion and atmospheric chemistry. For example, Guadagnini et al. have recently presented ab initio potential energy surfaces for the three lowest lying electronic states of HNO and then used these in examining several chemical reactions that take place in the combustion of nitrogen containing fuels and in the oxidation of atmospheric nitrogen. We have previously studied the ground state potential energy surface (i.e., stationary points along the HNO reversible reaction HON path), vibrational spectrum (using an accurate quartic force field), zero-point energy, and bonding of HNO using coupled-cluster ab initio methods. HNO is also very interesting because of the unique nature of its bonding characteristics. That is, the potential energy surface is very flat along the H-N bonding coordinate thereby giving unusual harmonic and fundamental vibrational frequencies, and the H-N bond energy is rather weak in comparison to other H-N bond energies. In fact, using experimental heats of formation for HNO, H, and NO, the H- bond energy is computed to be only 49.9 kcal/ mol (298 K). However, ab initio calculations of isodesmic reaction energies involving HNO, FNO, ClNO, and several other molecules have shown that there is an inconsistency in the experimental heats of formation of the XNO (X double bond H, F, and Cl) species. Hence the motivation for this study was to determine a very accurate (Delta)H(sub f, sup o) value for HNO using state-of-the-art ab initio methods. Based on many recent studies it is evident that the singles and doubles coupled-cluster method that includes a perturbational estimate of the effects of connected triple excitations, denoted CCSD(T), in conjunction with large one-particle basis sets should be reliable to better than +/- 0.8 kcal/mol for this quantity. The computational methodology is described in the next section followed by our results and discussion. Conclusions are presented in the final section.
Lee, Timothy J.↗
Insights into Methodologies and Stochastic Optimization of Thermal Energy Storage-Coupled Advanced Reactor Systems: A Comparison of Methods for Accessing Long-Term Sub-System Sizing Adequacy
This paper investigates the potential of coupling Thermal Energy Storage (TES) with Advanced Reactors (ARs) to address uncertainties posed by climate change in deep decarbonized power systems. The TES Use-case Team at Idaho National Laboratory (INL) has examined the potential of storing thermal energy from ARs during low demand periods and optimizing discharge during peak-priced hours, in both steady-state and transient conditions. Building on this groundwork, this study bridges the gaps in optimal sizing of the sub-system of TES-coupled AR systems using Risk Analysis Virtual Environment (RAVEN) and Holistic Energy Resource Optimization Network (HERON), INL?s framework for grid optimization. By applying this framework, we present statistically-robust optimal charge, discharge including balance of plant (BOP), and storage sizing for the High-Temperature Gas-Cooled Reactor (HTGR) with 203 MWth output. To this end, we generated synthetic price samples for 30 years using 2018 ? 2021 real-time market data from ERCOT, PJM and MISO. Our results reveals that the TES-coupled HTGR system is highly effective in maximizing revenue from electricity sales. We observed a substantial increase of 40 % in ERCOT and a noteworthy 15 % increase in PJM and MISO when compared to the conventional BOP without TES. This improvement is achieved through regionally-tailored sub-system sizing, which ranges from 398 to 416 MWth for the discharge system and 610 to 1029 MWth for the TES. We find that the average electricity price directly impacts the overall economics, while price volatility influences storage size. Additional sensitivity analyses were performed to access the impact of key assumptions on system economics and sizing, differentiating the optimization window (i.e., 24 ? 219 hours of chronological observations) and by imposing storage continuity condition in tracking TES cycles. We observed that at the 120-hour of the optimization window, a reasonable balance between computation time and accuracy was achieved. Our analysis also highlights the significance of conducting multi-day cycle analysis (> 120-hour) for TES to capture interaction between electricity prices and storage dynamics, providing a comprehensive understanding of TES behavior that AR developers should integrate into their plant designs.
25 ENERGY STORAGE↗
Nonadiabatic couplings from a variational excited state method based on constrained DFT
Excited Costrained Density Functional Theory (XCDFT) is a variational excited state method that extends ground state DFT to the computation of low-lying excited states. It borrows much of the machinery of Constrained DFT (CDFT) with a crucial difference: the constraint imposes a population of one electron in the Hilbert space spanned by the virtuals of a reference ground state. Here, we present theory and implementation for evaluating nonadiabatic coupling vectors (NACVs) between the first excited state computed with XCDFT and the ground state. Our NACVs are computed analytically using density functional perturbation theory with a formalism that is general enough that could be applied to CDFT diabatic states. We showcase the new method with pilot NACV calculations for the conical intersection in H 3 , the avoided crossing in selenoacrolein, and the NACV magnitudes in azobenzene. Despite complications from the nonorthogonality of the wavefunctions, XCDFT’s energy surfaces and NACVs reproduce benchmark values and respect known sum rules within a reasonable degree. This shows that XCDFT is a viable method for nonadiabatic dynamics simulations.
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
Computation of viscous flows over airfoils, including separation, with a coupling approach
Viscous incompressible flows over single or multiple airfoils, with or without separation, were computed using an inviscid flow calculation, with modified boundary conditions, and by a method providing calculation and coupling for boundary layers and wakes, within conditions of strong viscous interaction. The inviscid flow is calculated with a method of singularities, the numerics of which were improved by using both source and vortex distributions over profiles, associated with regularity conditions for the fictitious flows inside of the airfoils. The viscous calculation estimates the difference between viscous flow and inviscid interacting flow, with a direct or inverse integral method, laminar or turbulent, with or without reverse flow. The numerical method for coupling determines iteratively the boundary conditions for the inviscid flow. For attached viscous layers regions, an underrelaxation is locally calculated to insure stability. For separated or separating regions, a special semi-inverse algorithm is used. Comparisons with experiments are presented.
Leballeur, J. C.↗
Solving Nonlinear Coupled Differential Equations
Harmonic balance method developed to obtain approximate steady-state solutions for nonlinear coupled ordinary differential equations. Method usable with transfer matrices commonly used to analyze shaft systems. Solution to nonlinear equation, with periodic forcing function represented as sum of series similar to Fourier series but with form of terms suggested by equation itself.
Mitchell, L.↗
Norton-Thevenin Receptance Coupling (NTRC) as a Payload Design Tool
The NASA Engineering and Safety Center (NESC) is funding a study to develop an alternate method for performing coupled loads analysis called Norton-Thevenin Receptance Coupling (NTRC). NTRC combines Receptance Coupling (RC), a frequency-domain synthesis method and Norton-Thevenin (NT) theory, an impedance based approach for simulating the interaction between dynamic systems. The goal of developing the NTRC method is to provide a tool that payload developers can use to reduce the conservatism in defining preliminary design loads, assess the impact of design changes between formal load cycles, and to perform trade studies for design optimization with a minimum amount of data required from the launch vehicle (LV) provider. NTRC also has the ability to perform parametric loads analysis where many different design configurations can be evaluated. This will result in cost and schedule benefits to the payload developer that are currently not possible under the standard coupled loads analysis (CLA) flow where typically only 2-3 official load cycles are performed by the LV provider over the life of a payload program. NTRC is not envisioned as a replacement for the official load cycles performed by the LV provider but rather as a means to address the types of design issues faced by the payload developer before and between official load cycles.The presentation provides an overview of the NTRC methodology and discusses how NTRC can be used to replicate the results from a standard LV CLA. The presentation covers the benchmarking that has been performed as part of the NESC study to demonstrate the accuracy of the technique for both frequency and time domain dynamic analyses. Future plans for benchmarking the NTRC approach against CLA results for NASAs Space Launch System (SLS) and commercial launch vehicles are discussed and the role that NTRC is envisioned to play in the payload development cycle.
Structural Dynamics↗
Resource-Efficient Chemistry on Quantum Computers with the Variational Quantum Eigensolver and The Double Unitary Coupled-Cluster approach
Applications of quantum simulation algorithms to obtain electronic energies of molecules on noisy intermediate-scale quantum (NISQ) devices require careful consideration of resources describing the complex electron correlation effects. In modeling second-quantized problems, the biggest challenge confronted is that the number of qubits scales linearly with the size of molecular basis. This poses a significant limitation on the size of the basis sets and the number of correlated electrons included in quantum simulations of chemical processes. To address this issue and to enable more realistic simulations on NISQ computers, we employ the double unitary coupled-cluster (DUCC) method to effectively downfold correlation effects into the reduced-size orbital space, commonly referred to as the active space. Using downfolding techniques, we demonstrate that properly constructed effective Hamiltonians can capture the effect of the whole orbital space in small-size active spaces. Combining the downfolding pre-processing technique with the Variational Quantum Eigensolver, we solve for the ground-state energy of H2 and Li2 in the cc-pVTZ basis using the DUCC-reduced active spaces. We compare these results to full configuration-interaction and high-level coupled-cluster reference calculations.
quantum computing, variational quantum solver, cou↗
Methodology to determine the coupling of continental clouds with surface and boundary layer height under cloudy conditions from lidar and meteorological data
Abstract. The states of coupling between clouds and surface or boundary layer have been investigated much more extensively for marine stratocumulus clouds than for continental low clouds, partly due to more complex thermodynamic structures over land. A manifestation is a lack of robust remote sensing methods to identify coupled and decoupled clouds over land. Following the idea for determining cloud coupling over the ocean, we have generalized the concept of coupling and decoupling to low clouds over land, based on potential temperature profiles. Furthermore, by using ample measurements from lidar and a suite of surface meteorological instruments at the U.S. Department of Energy's Atmospheric Radiation Measurement Program's Southern Great Plains site from 1998 to 2019, we have developed a method to simultaneously retrieve the planetary boundary layer (PBL) height (PBLH) and coupled states under cloudy conditions during the daytime. The new lidar-based method relies on the PBLH, the lifted condensation level, and the cloud base to diagnose the cloud coupling. The coupled states derived from this method are highly consistent with those derived from radiosondes. Retrieving the PBLH under cloudy conditions, which has been a persistent problem in lidar remote sensing, is resolved in this study. Our method can lead to high-quality retrievals of the PBLH under cloudy conditions and the determination of cloud coupling states. With the new method, we find that coupled clouds are sensitive to changes in the PBL with a strong diurnal cycle, whereas decoupled clouds and the PBL are weakly related. Since coupled and decoupled clouds have distinct features, our new method offers an advanced tool to separately investigate them in climate systems.
54 ENVIRONMENTAL SCIENCES↗
State-of-the-art methods for computing the electromagnetic interaction of lightning with aircraft
Nuclear electromagnetic pulse (NEMP) coupling codes and methods are evaluated and summarized. The differences between NEMP and lightning interaction with aircraft are discussed and critical parameters peculiar to lightning are examined.
Eriksen, F. J.↗
Design optimization of natural laminar flow fuselages in compressible flow
An optimization method has been developed to design axisymmetric body shapes such as fuselages, nacelles, and external fuel tanks with increased transition Reynolds numbers in subsonic compressible flow. The new design method involves a constraint minimization procedure coupled with analysis of the inviscid and viscous flow regions, and linear stability analysis of the compressible boundary layer. Boundary-layer transition is predicted by a 'hybrid' transition criterion based on Granville's transition criterion and a criterion using linear stability theory coupled with the e(n)-method. A tiptank of a business-jet is used as an example to illustrate that the method can be utilized to design an axisymmetric body shape with extensive natural laminar flow. Boundary layer transition is predicted to occur at a transition Reynolds number of 6.04 x 10 to the 6th on the original tiptank. On the designed body shape, a transition Reynolds number of 7.22 x 10 to the 6th is predicted using compressible linear stability theory coupled with e(n)-method.
Dodbele, Simha S.↗
Radiation/convection coupling in rocket motor and plume analysis
A method for describing radiation/convection coupling to a flow field analysis was developed for rocket motors and plumes. The three commonly used propellant systems (H2/O2, RP-1/O2, and solid propellants) radiate primarily as: molecular emitters, non-scattering small particles (soot), and scattering larger particles (Al2O3), respectively. For the required solution, the divergence of the radiation heat flux was included in the energy equation, and the local, volume averaged intensity was determined by a solution to the radiative transfer equation. A rigorous solution to this problem is intractable, therefore, solution methods which use the ordinary and improved differential approximation were developed. This radiation model was being incorporated into the FDNS code, a Navier-Stokes flowfield solver for multiphase, turbulent combusting flows.
Saladino, A. J.↗