Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Perturbative methods”

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,027 records · Page 57
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Crossing the Streams – Sampler and the TemplateEngine [Slides]

This presentation discusses Sampler, which is a versatile UQ and parametric study tool that can be applied to any SCALE Sequence. Sampler can perturb any quantity in any SCALE input. Recent work at ORNL has developed new types of covariance data that allow Sampler UQ to be applied to nearly all SCALE applications, including reactor depletion, UNF fuel characterization, source term analysis, and decay heat calculation. In SCALE 6.2 releases, CE data in transport cannot be perturbed. Sampler was originally designed for stochastic sampling with any sequence within SCALE and Parametric capability added in SCALE 6.2.2. Sampler can be used for uncertainty quantification, including sample data in static or depletion calculations and sample inputs for uncertainties in compositions and dimensions. The SCALE TemplateEngine allows for expanding templates to full inputs and the combination provides a powerful UQ tool.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

Novel method to control explosive shock sensitivity: A mesoscale study to understand the effect of thermally expandable microsphere (TEM) inclusions in high explosives (HE) microstructure

When the void content and/or void structure of a high explosive (HE) is altered by some means (i.e., bulk heating or mechanical damage), the shock initiation behavior of the material changes. The ability to precisely predict the change in shock sensitivity after an HE has undergone microstructural changes is a crucial capability in multi-scale reactive flow models. Here, we utilize thermally expandable microspheres (TEMs) as a dopant in a polymer bonded explosive (PBX) matrix to alter the shock initiation properties in a controlled fashion. Here, using a mesoscale modeling approach, we evaluated how a single TEM (before and after thermal expansion) behaves under shock compression, as well as how the matrix PBX in the direct vicinity of the TEM is affected. We first examined the effect of an unexpanded TEM in the explosive matrix and found that its presence does not significantly perturb the bulk flow and by extension will not affect bulk sensitivity. Next, we examined the effect of an expanded TEM and found that its presence significantly perturbs the flow via hydrodynamic jetting, which causes a secondary shock wave with a strength that exceeds that of the incident wave. Finally, we showed that this secondary shock interacts with the downstream porosity to ignite a larger fraction of the overall pore volume, commensurate with the secondary shock strength and the affected volume, increasing the global (bulk) shock sensitivity.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Skewing the CMB×LSS: a fast method for bispectrum analysis

Upcoming cosmic microwave background (CMB) lensing measurements and tomographic galaxy surveys are expected to provide us with high-precision data sets in the coming years, thus paving the way for fruitful cross-correlation analyses. Here, in this paper we study the information content of the weighted skew-spectrum, a nearly-optimal estimator of the angular bispectrum amplitude, as a means to extract non-Gaussian information on both bias and cosmological parameters from the bispectra of galaxies cross-correlated with CMB lensing, while gaining significantly on speed. Our results show that for the combination of the Planck satellite and the Dark Energy Spectroscopic Instrument (DESI), the difference in the constraints on bias and cosmological parameters from the skew-spectrum and the bispectrum is at most 17%. We further compare and find agreement between our theoretical skew-spectra and those estimated from N-body simulations, for which it is important to include gravitational non-linearities beyond perturbation theory and the post-Born effect for CMB lensing. We define an algorithm to apply the skew-spectrum estimator to the data and, as a preliminary step, we use the skew-spectra to constrain bias parameters and the amplitude of shot noise from the simulations through a Markov chain Monte Carlo likelihood analysis, finding that it may be possible to reach percent-level estimates for the linear bias parameter b 1 .

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Synchronization in electric power networks with inherent heterogeneity up to 100% inverter-based renewable generation

The synchronized operation of power generators is the foundation of electric power network stability and a key to the prevention of undesired power outages and blackouts. Here, we derive the conditions that guarantee synchronization in power networks with inherent generator heterogeneity when subjected to small perturbations, and perform a parametric sensitivity analysis to understand synchronization with varied types of generators. As inverter-based resources, which are the primary interfacing technology for many renewable sources of energy, have supplanted synchronous generators in ever growing numbers, the center of attention on associated integration challenges have resided primarily on the role of declining system inertia. Our results instead highlight the critical role of generator damping in achieving a stable state of synchronization. Additionally, we report the feasibility of operating interconnected electric grids with up to 100% power contribution from inverter-based renewable generation technologies. Our study has important implications as it sets the basis for the development of advanced control architectures and grid optimization methods that ensure synchronization and further pave the path towards the decarbonization of the electric power sector.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2022

Development of theory and experimental operational framework for Coherent Thomson Scattering (Final Technical Report)

The research carried out explored and proved the feasibility and operational framework of a new diagnostic technique termed Coherent Thomson Scattering (CTS) for electrons in a low temperature plasma. The work is performed in collaboration with the Princeton Collaborative Research Facility (PCRF) at Princeton Plasma Physics Laboratory. The novel technique builds on an established and demonstrated single shot diagnostic method, called Coherent Rayleigh-Brillouin scattering, which has successfully been applied in neutral flows. The proposed novel four wave mixing diagnostic technique of CTS will allow for higher spatial resolution and lower detectable number densities for the electrons than conventional Thomson scattering. In this project we developed the theoretical framework for Coherent Thomson Scattering as well as the specification of the appropriate operational experimental parameters for successful CTS implementation in e.g. a low temperature plasma. Additionally, the mode of operation and the detection limits for a practical CTS experimental demonstration were explored. Ultimately, successful experimental demonstration of CTS can be seen as transformative in a multitude of plasma physics areas, since it will allow for detailed, non-perturbative measurements of electron density and temperature, previously unattainable by other measurement techniques. This project was the first successful step towards this direction.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Efficient continuous Energy-Multigroup hybrid depletion scheme using the Shift Monte Carlo code. Part I: Energy condensation sensitivity analysis

Monte Carlo (MC) codes coupled to depletion solvers are increasingly used to provide high fidelity fuel cycle modeling capabilities. Here, these coupled depletion-MC tools produce accurate results in general but can experience nonphysical spatial oscillations when time steps are large or when a system’s dominance ratio approaches unity. Two substepping techniques have been developed previously to remedy and dampen these spatial oscillations without needing to reduce step sizes. The first approach relied on higher-order techniques to account for spectral changes within steps (extrapolation and interpolation techniques). The second approach used the first order perturbation (FOP) theory to account for the change in the one-group spatial flux distribution within steps. This paper develops a hybrid depletion methodology which, in a way, combines how the flux is handled in both substepping techniques. Specifically, the multigroup (MG) MC Shift code is used to update the flux distribution within steps rather than a one-group FOP solver. A fully reflected pincell is investigated, which is not spatially dependent in the MG representation. Thus, the analysis in this paper is an initial demonstration of hybrid depletion. An upcoming companion paper will focus on how the hybrid depletion dampens spatial oscillations. The hybrid depletion approach is verified to be consistent with previous constant extrapolation depletion (CED) methods. This paper finds that the hybrid CED exhibits some error in the eigenvalue and one group constants within macro steps. To address this discrepancy, a simple interpolation scheme (CELI) is investigated. This work found that CELI sufficiently addresses the discrepancy in spectrum for macro steps up to 100 days. Overall, this work demonstrates that the hybrid depletion method can significantly reduce the number of high fidelity MC executions in a MC-coupled depletion with an acceptable eigenvalue error.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗
DOE OSTI2020

N 5 -Scaling Excited-State-Specific Perturbation Theory

In this paper, we show that by working in a basis similar to that of the natural transition orbitals and using a modified zeroth-order Hamiltonian, the cost of a recently introduced perturbative correction to excited-state mean field theory can be reduced from seventh to fifth order in the system size. The (occupied) 2 (virtual) 3 asymptotic scaling matches that of ground-state second-order Møller-Plesset theory but with a significantly higher prefactor because the bottleneck is iterative: it appears in the Krylov-subspace-based solution of the linear equation that yields the first-order wave function. Here, we discuss the details of the modified zeroth-order Hamiltonian we use to reduce the cost and the automatic code generation process we used to derive and verify the cost scaling of the different terms. Overall, we find that our modifications have little impact on the method's accuracy, which remains competitive with singles and doubles equation-of-motion coupled cluster.

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

Active control of resistive wall mode via modification of external tearing index

Modification of the external tearing index, Δ ext ′ , by magnetic feedback is analytically investigated for the purpose of controlling the resistive plasma resistive wall mode (RP-RWM). The matching method is pursued by deriving expressions for the close-loop Δ ext ′ and by linking it to the counterpart from the inner layer. Various feedback coil configurations are found to generally reduce Δ ext ′ and stabilize the RWM, with either proportional or derivative control. Feedback modification of Δ ext ′ is found to be generally independent of the inner layer resistive interchange index D R , confirming that feedback action primarily modifies the solution in the outer ideal region for the RP-RWM. Exception occurs when either the inner layer favorable curvature effect becomes sufficiently large or the feedback action is sufficiently strong to introduce a rotating RP-RWM in the static plasma, leading to complex-valued close-loop Δ ext ′ . The perturbed magnetic energy dissipation in the outer region, associated with the eddy current in the resistive wall, is identified as the key physics reason for feedback induced complex Δ ext ′ . Similar results are also obtained for active control of the external kink instability, whose open-loop growth rate is significantly reduced by inclusion of the plasma resistivity. Within the single poloidal harmonic approximation, which is most suitable for the matching approach, external active coils combined with poloidal sensors are often found to be more efficient for feedback stabilization of the mode at large proportional gain values. This counter-intuitive result is explained as the lack of (non-resonant) poloidal harmonics for proper description of the feedback coil geometry.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

High-order ( N = 4–6) multiphoton absorption and mid-infrared Kerr nonlinearity in GaP, ZnSe, GaSe, and ZGP crystals

We report a study of high-order multiphoton absorption, nonlinear refraction, and their anisotropy in four notable mid-infrared χ (2) crystals: GaP, ZnSe, GaSe, and ZGP using femtosecond pulses at 2.35 μm with peak intensity in excess of 200 GW/cm 2 . We found that the nonlinear absorption obeys a perturbation model with multiphoton absorption (MPA) orders from N = 4–6, in agreement with the bandgaps of the crystals. A study of the role of free carrier absorption, performed by changing the pulse duration from 30 to 70 fs while maintaining a constant peak intensity, showed that, at our intensity levels, free carriers generated during the MPA process absorb much more strongly than would be expected from their known linear absorption cross section. Possible mechanisms for this anomalous behavior are high-field effects, such as intravalley scattering in the conduction band and absorption to higher lying bands. Nonlinear refractive indices were measured using (i) closed aperture Z-scan and (ii) spectral broadening due to self-phase modulation, both methods agreeing with each other.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

Assessment of the Impact of Realistic Sensor Physics and the Integration of Ex-Core Sensors on Reactor Power Synthesis

In the work documented in this report, a weighting function–based core power synthesis method was applied to multiple Monte Carlo N-Particle (MCNP) reactor models, which are informed based on simulated self-powered neutron detector (SPND) responses. The weighting function method used has been coined the point-based iterative (PBI) method. The goal of this application is to assess the impact of considering realistic sensor physics in the generation of the simulated SPND outputs as well as to consider how the synthesis is impacted based on the inclusion of ex-core detectors in the model. The NuScale small modular reactor (SMR) and Westinghouse AP1000 pressurized water reactor (PWR) are the models that served as the testbeds for the assessment of realistic sensor physics; this was achieved by using Geant4 SPND models in comparison with analytical models, such that the effect of electron transport in realistic SPND geometries in the Geant4 model can be understood in terms of synthesis error and convergence time. The comparison was considered for fuel burnup–induced perturbations, for a range of sensor string densities and synthesized power distribution axial fidelities. The Texas A&M Testing, Research, Isotopes, General Atomics Reactor (TAMU TRIGA) reactor MCNP model was used to assess the impact of ex-core sensors; this was done by performing synthesis with and without the ex-core detectors and by quantifying the synthesis error and number of iterations associated with Gaussian-type perturbations in many locations in the core. The TAMU TRIGA model was particularly pertinent for this study because of the interest in future experimental tests with SPNDs in this reactor, as well as the ease of modifying the MCNP model to include ex-core detectors with heterogeneously described response functions. Results from the comparison between the Geant4 and analytical SPND models indicate that similar average and maximum synthesis errors were obtained for burnup-induced perturbations in both the NuScale SMR and the AP1000. This was true for a range of sensor string densities and axial fidelities. However, there were marked differences between both the Geant4 and analytically informed models in terms of the iterations required to converge on the synthesized power distribution. Namely, the Geant4-informed models tended to lead to fewer iterations, except for a few sensor–core configurations that had particularly numerous iterations. Results from the ex-core sensor assessment with the TAMU TRIGA model indicate that the inclusion of ex-core sensors drastically reduces the synthesis error of Gaussian-type perturbations close to the edge of the core, and it slightly reduces synthesis errors for perturbations closer to the center of the core. This was achieved with a minimal increase in computational cost—that is, the number of iterations required for convergence. The errors were identified to be in the same location as the perturbation in the core, indicating that the methodology remains robust for unperturbed regions of the core. A secondary result from this study with the TAMU TRIGA was yielded by analysis of the neutron flux levels in the in-core and ex-core sensor locations of the core; these flux levels indicate that SPNDs could be used as both in-core and ex-core sensors, so long as the emitter material is sensitive to thermal neutrons. The results from these studies provide a quantitative understanding of the importance of considering realistic sensor physics and including ex-core sensors to perform accurate and timely power distribution synthesis of a reactor core.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2025

Analytic reconstruction with massive particles: one-loop amplitudes for $0\to \overline{q} qt\overline{t}H$

We present an analytic reconstruction of one-loop amplitudes for the process $0\to \overline{q} qt\overline{t}H$. Our calculation is a novel use of analytic reconstruction, retaining explicit covariance in the massive spin states through the massive spinor-helicity formalism. The analytic reconstruction relies on embedding the massive five-point kinematics in a fully massless eight-point phase space while still building a minimal ansatz directly in the five-point phase space. In order to obtain compact analytic expressions it is necessary to identify suitable partial fraction decompositions and extract common numerator factors, which we achieve through careful inspection of limits in which pairs of denominators vanish. We find that the resulting amplitudes are more numerically efficient than ones computed using automatic methods but that the gains are not as significant as in the massless case, at least at present. The method opens the door to applications at two-loop order, where numerical efficiency and improvements in the reconstruction methodology are more crucial, especially with regards to the number of free parameters in the ansatz.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2020

A Generalized Variational Principle with Applications to Excited State Mean Field Theory

We present a generalization of the variational principle that is compatible with any Hamiltonian eigenstate that can be specified uniquely by a list of properties. This variational principle appears to be compatible with a wide range of electronic structure methods, including mean field theory, density functional theory, multireference theory, and quantum Monte Carlo. Like the standard variational principle, this generalized variational principle amounts to the optimization of a nonlinear function that, in the limit of an arbitrarily flexible wave function, has the desired Hamiltonian eigenstate as its global minimum. Unlike the standard variational principle, it can target excited states and select individual states in cases of degeneracy or near-degeneracy. As an initial demonstration of how this approach can be useful in practice, we employ it to improve the optimization efficiency of excited state mean field theory by an order of magnitude. With this improved optimization, we are able to demonstrate that the accuracy of the corresponding second-order perturbation theory rivals that of singles-and-doubles equation-of-motion coupled cluster in a substantially broader set of molecules than could be explored by our previous optimization methodology.

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

Finite-temperature vibrational full configuration interaction

Thermodynamic functions of an ideal molecular gas due to its anharmonic vibrations are evaluated in a wide range of temperature (T) by the vibrational full-configuration-interaction (FCI) method using a quartic force field and a finite number (N) of harmonic-oscillator basis functions along each normal mode. The thermodynamic functions considered are the grand potential (Ω), internal energy (U), and entropy (S). They are compared with those obtained from the Bose–Einstein theory with or without truncation of the harmonic-oscillator basis functions after quantum number N–1. The comparison reveals that the finite-basis-set errors in Ω and U are, respectively, k B Tln(k B T/Nℏω) and k B T per mode in the high-T limit, obscuring anharmonic effects when k B T > ℏω, where ω is the lowest mode frequency. Here, the benchmark data for several low-order perturbation corrections to Ω, U, and S are also obtained as the numerical derivatives of their FCI values with respect to dimensionless perturbation strength, and the domain of T and N in which these data are reliable (for the N → ∞ limits) is discussed.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2022

Asymptotic preserving methods for fluid electron-fluid models in the large magnetic field limit with mathematically guaranteed properties (Final Report)

The current manuscript is a final report on the activities carried out under the Project LDRD-CIS #226834. In scientific terms, the work reported in this manuscript is a continuation of the efforts started with Project LDRD-express #223796 with final report of activities SAND2021-11481, see [83]. In this section we briefly explain what pre-existing developments motivated the current body of work and provide an overview of the activities developed with the funds provided. The overarching goal of the current project LDRD-CIS #226834 and the previous project LDRD-express #223796 is the development of numerical methods with mathematically guaranteed properties in order to solve the Euler-Maxwell system of plasma physics and generalizations thereof. Even though Project #223796 laid out general foundations of space and time discretization of Euler-Maxwell system, overall, it was focused on the development of numerical schemes for purely electrostatic fluid-plasma models. In particular, the project developed a family of schemes with mathematically guaranteed robustness in order to solve the Euler-Poisson model. This model is an asymptotic limit where only electrostatic response of the plasma is considered. Its primary feature is the presence of a non-local force, the electrostatic force, which introduces effects with infinite speed propagation into the problem. Even though instantaneous propagation of perturbations may be considered nonphysical, there are plenty of physical regimes of technical interest where such an approximation is perfectly valid.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

The many-body expansion for metals. I. The alkaline earth metals Be, Mg, and Ca

We examine the many-body expansion (MBE) for alkaline earth metal clusters, Be n , Mg n , Ca n ( n = 4, 5, 6), at the Møller–Plesset second order perturbation theory, coupled-cluster singles and doubles with perturbative triples, multi-reference perturbation theory, and multi-reference configuration interaction levels of theory. The magnitude of each term in the MBE is evaluated for several geometrical configurations. We find that the behavior of the MBE for these clusters depends strongly on the geometrical arrangement and, to a lesser extent, on the level of theory used. Another factor that affects the MBE is the in situ (ground or excited) electronic state of the individual atoms in the cluster. For most geometries, the three-body term is the largest, followed by a steady decrease in absolute energy for subsequent terms. Though these systems exhibit non-negligible multi-reference effects, there was little qualitative difference in the MBE when employing single vs multi-reference methods. Useful insights into the connectivity and stability of these clusters have been drawn from the respective potential energy surfaces and quasi-atomic orbitals for the various dimers, trimers, and tetramers. Through these analyses, we investigate the similarities and differences in the binding energies of different-sized clusters for these metals.

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

Lattice-QCD Computable Quark Correlation Functions at Three-Loop Order and Extraction of Splitting Functions

We present the first complete next-to-next-to-next-to-leading-order calculation of the matching coefficients that link unpolarized flavor nonsinglet parton distribution functions with lattice QCD computable correlation functions. By using this high-order result, we notice a reduction in theoretical uncertainties compared to relying solely on previously known lower-order matching coefficients. Furthermore, based on this result we have extracted the three-loop unpolarized flavor nonsinglet splitting function, which is in agreement with the state-of-the-art result. Because of the simplicity of our method, it has the potential to advance the calculation of splitting functions to the desired four-loop order.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2023

Light-sheet autofluorescence lifetime imaging with a single-photon avalanche diode array

Significance: Fluorescence lifetime imaging microscopy (FLIM) of the metabolic co-enzyme nicotinamide adenine dinucleotide (phosphate) [NAD(P)H] is a popular method to monitor single-cell metabolism within unperturbed, living 3D systems. However, FLIM of NAD(P)H has not been performed in a light-sheet geometry, which is advantageous for rapid imaging of cells within live 3D samples. Aim: We aim to design, validate, and demonstrate a proof-of-concept light-sheet system for NAD(P)H FLIM. Approach: A single-photon avalanche diode camera was integrated into a light sheet microscope to achieve optical sectioning and limit out-of-focus contributions for NAD(P)H FLIM of single cells. Results: An NAD(P)H light-sheet FLIM system was built and validated with fluores cence lifetime standards and with time-course imaging of metabolic perturbations in pancreas cancer cells with 10 s integration times. NAD(P)H light-sheet FLIM in vivo was demonstrated with live neutrophil imaging in a larval zebrafish tail wound also with 10 s integration times. Finally, the theoretical and practical imaging speeds for NAD(P)H FLIM were compared across laser scanning and light-sheet geometries, indicating a 30× to 6× acquisition speed advantage for the light sheet compared to the laser scanning geometry. Conclusions: FLIM of NAD(P)H is feasible in a light-sheet geometry and is attrac tive for 3D live cell imaging applications, such as monitoring immune cell metabolism and migration within an organism.

47 OTHER INSTRUMENTATION↗
DOE OSTI2023

Synchrotron radiography of Richtmyer–Meshkov instability driven by exploding wire arrays

We present a new technique for the investigation of shock-driven hydrodynamic phenomena in gases, liquids, and solids in arbitrary geometries. The technique consists of a pulsed power-driven resistive wire array explosion in combination with multi-MHz synchrotron radiography. Compared to commonly used techniques, it offers multiple advantages: (1) the shockwave geometry can be shaped to the requirements of the experiment, (2) the pressure (P > 300 MPa) generated by the exploding wires enables the use of liquid and solid hydrodynamic targets with well-characterized initial conditions (ICs), (3) the multi-MHz radiography enables data acquisition to occur within a single experiment, eliminating uncertainties regarding repeatability of the ICs and subsequent dynamics, and (4) the radiographic measurements enable estimation of compression ratios from the x-ray attenuation. In addition, the use of a synchrotron x-ray source allows the hydrodynamic samples to be volumetrically characterized at a high spatial resolution with synchrotron-based microtomography. This experimental technique is demonstrated by performing a planar Richtmyer–Meshkov instability (RMI) experiment on an aerogel–water interface characterized by Atwood number A0∼−0.8 and Mach number M∼1.5. The qualitative and quantitative features of the experiment are discussed, including the energy deposition into the exploding wires, shockwave generation, compression of the interface, startup phase of the instability, and asymptotic growth consistent with Richtmyer's impulsive theory. Additional effects unique to liquids and solids—such as cavitation bubbles caused by rarefaction flows or initial jetting due to small perturbations—are observed. It is also demonstrated that the technique is not shape dependent by driving a cylindrically convergent RMI experiment.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗