SEARCH · Engineering Papers

Results for “simulation analytics”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 595 records · Page 33

All sources NASA NTRS DOE OSTI

DOE OSTI2025

An analytically tractable marked power spectrum

The increasing precision of cosmology data in the modern era is calling for methods to allow the extraction of non-Gaussian information using tools beyond two-point statistics. The marked power spectrum has the potential to extract beyond two-point information in a computationally efficient way while using much of the infrastructure already available for the power spectrum. In this work we explore the marked power spectrum from an analytical perspective. In particular, we explore a low-order polynomial for the mark that allows us to better control the theoretical uncertainties and we show that with minimal new degrees of freedom the analytical results match measurements from N-body simulations for both the matter field and biased tracers in redshift space. Finally, we show that even within the limited forms of mark that we consider, there are degeneracies that can be broken by inclusion of the marked auto-spectrum or the cross-spectrum with the unmarked field. I n conclusion, we discuss future theoretical developments that would enable us to apply this approach to survey data.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2021

Thermodynamics of Light Management in Near-Field Thermophotovoltaics

We evaluate near-field thermophotovoltaic (TPV) energy-conversion systems focusing in particular on their open-circuit voltage (V OC ). Unlike previous analyses based largely on numerical simulations with fluctuational electrodynamics, here, we develop an analytic model that captures the physics of near-field TPV systems and can predict their performance metrics. Using our model, we identify two opportunities of TPV systems operating in the near field. First, we show analytically that enhancement of radiative recombination is a natural consequence of operating in the near field. Second, we note that, owing to photon recycling and minimal radiation leakage in near-field operation, the PV cell used in near-field TPV systems can be much thinner compared to those used in solar PV systems. Since nonradiative recombination is a volumetric effect, use of a thinner cell reduces nonradiative losses per unit area. The combination of these two opportunities leads to increasingly large values of V OC as the TPV vacuum gap decreases. Hence, although operation in the near-field was previously perceived to be beneficial for electrical power-density enhancement, here, we emphasize that thin-film near-field TPVs are also significantly advantageous in terms of V OC and consequently conversion efficiency as well as power density. We provide numerical results for an InAs-based thin-film TPV that exhibits efficiency >50% at an emitter temperature as low as 1100 K.

42 ENGINEERING↗

DOE OSTI2020

Infrared plasmonics: STEM-EELS characterization of Fabry-Pérot resonance damping in gold nanowires

Materials possessing strong midinfrared responses are of current interest because of their potential application to long-wavelength metamaterials, photonic devices, molecular detection, and catalysis. In this work, we utilize high-energy resolution (80 cm -1 , 10 meV) electron-energy-loss spectroscopy (EELS) in a monochromated and aberration-corrected scanning transmission electron microscope (STEM) to resolve multipolar surface plasmon resonances (SPRs), sometimes called Fabry-Pérot (FP) resonances, in gold nanowires with mode energies spanning from ~1000 to 8000 cm -1 . STEM-EELS provides access to these mid- to near-IR responses in a single acquisition, avoiding the difficulties inherent in obtaining the same data using near-field optical techniques. The experimentally measured FP resonance energies and linewidths, together with analytical modeling and full-wave numerical electrodynamics simulations, provide a comprehensive picture of the radiative and intrinsic contributions to the total damping rates. We find some FP modes with dephasing times >60fs, which is almost twice the longest previously reported plasmon dephasing time for individual Au nanoparticles in the infrared. The long dephasing times and the broad tunability of the FP resonance energies throughout the infrared region suggest additional opportunities for harnessing infrared plasmonic energy before dephasing occurs.

36 MATERIALS SCIENCE↗

DOE OSTI2024

Signatures of primordial energy injection from axion strings

Axion strings are horizon-size topological defects that may be produced in the early Universe. Ultralight axion-like particles may form strings that persist to temperatures below that of big bang nucleosynthesis. Such strings have been considered previously as sources of gravitational waves and cosmic microwave background (CMB) polarization rotation. In this work we show, through analytic arguments and dedicated adaptive mesh refinement cosmological simulations, that axion strings deposit a subdominant fraction of their energy into high-energy Standard Model (SM) final states, for example, by the direct production of heavy radial modes that subsequently decay to SM particles. This high-energy SM radiation is absorbed by the primordial plasma, leading to novel signatures in precision big bang nucleosynthesis, the CMB power spectrum, and gamma-ray surveys. In particular, we show that CMB power spectrum data constrains axion strings with decay constants f a ≲ 10 12 GeV , up to model dependence on the ultraviolet completion, for axion masses m a ≲ 10 − 29 eV ; future CMB surveys could find striking evidence of axion strings with lower decay constants. Published by the American Physical Society 2024

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2023

Stability of the modulator in a plasma-modulated plasma accelerator

We explore the regime of operation of the modulator stage of a recently proposed laser-plasma accelerator scheme [Phys. Rev. Lett. 127, 184801 (2021)], dubbed the plasma-modulated plasma accelerator (P-MoPA). The P-MoPA scheme offers a potential route to high-repetition-rate, GeV-scale plasma accelerators driven by picosecond-duration laser pulses from, for example, kilohertz thin-disk lasers. The first stage of the P-MoPA scheme is a plasma modulator in which a long, high-energy “drive” pulse is spectrally modulated by copropagating in a plasma channel with the low-amplitude plasma wave driven by a short, low-energy “seed” pulse. The spectrally modulated drive pulse is converted to a train of short pulses, by introducing dispersion, which can resonantly drive a large wakefield in a subsequent accelerator stage with the same on-axis plasma density as the modulator. In this paper we derive the 3D analytic theory for the evolution of the drive pulse in the plasma modulator and show that the spectral modulation is independent of transverse coordinate, which is ideal for compression into a pulse train. We then identify a transverse mode instability (TMI), similar to the TMI observed in optical fiber lasers, which sets limits on the energy of the drive pulse for a given set of laser-plasma parameters. We compare this analytic theory with particle-in-cell (PIC) simulations and find that even higher energy drive pulses can be modulated than those demonstrated in the original proposal.

43 PARTICLE ACCELERATORS↗

DOE OSTI2025

Flux Hypothesis for Odd Transport Phenomena

Onsager's regression hypothesis makes a fundamental connection between macroscopic transport phenomena and the average relaxation of spontaneous microscopic fluctuations. This relaxation, however, is agnostic to odd transport phenomena, in which fluxes run orthogonal to the gradients driving them. To account for odd transport, we generalize the regression hypothesis, postulating that macroscopic linear constitutive laws are, on average, obeyed by microscopic fluctuations, whether they contribute to relaxation or not. From this "flux hypothesis," Green-Kubo and reciprocal relations follow, elucidating the separate roles of broken time-reversal and parity symmetries underlying various odd transport coefficients. As an application, we derive and verify the Green-Kubo relation for odd collective diffusion in chiral active matter, first in an analytically tractable model and subsequently through molecular dynamics simulations of concentrated active spinners.

Hargus, Cory↗

DOE OSTI2026

Startup regime of high-efficiency tapering-enhanced FEL oscillator

In this paper, we present a design of a high-efficiency high-gain free-electron laser oscillator based on the use of a strongly tapered undulator for extracting energy from high-brightness electron beams. We provide an analytical model of the setup followed by numerical simulations for lasing at the wavelength of 13.5 nm. We discuss the optimization of the system in steady state and the conditions necessary for the pass-per-pass buildup of the power from shot noise level. We propose the use of fast phase shifters as a way to accelerate the buildup. Finally, we present time-dependent simulations of the oscillator and discuss the role of spectral filtering. The optimized working point yields a total energy conversion efficiency from the electron beam to output radiation above 1% at the wavelength of 13.5 nm.

Beam dynamics↗

DOE OSTI2020

Analytical Voltage Sensitivity Analysis for Unbalanced Power Distribution System

Large scale integration of distributed energy resources and electric vehicles in a transactive energy environment present new challenges in terms of voltage stability and fluctuations in a power distribution system. The impact of different level of DER/EV penetration on the voltages across the network is typically quantified through voltage sensitivity analyses. Existing methods of voltage sensitivity analysis are computationally expensive and prior efforts to develop analytical approximation lacks generality and have not been effectively validated. The objective of this work is to provide a new analytical method of voltage sensitivity analysis that has low computational cost and also allows for stochastic analysis of voltage change. This paper first derives an analytical approximation of change in voltage at a particular bus due to change in power consumption at other bus in a radial three phase unbalanced power distribution system. Then, the proposed method is shown to be valid for different load configurations, which demonstrates its generality. The results from our analytical approach is validated via classical load flow simulation of the test system based on IEEE 37 bus network. The proposed method is shown to have good accuracy, and computation complexity is of order O(1), compared to O(n3) in classical sensitivity analysis approaches.

Munikoti, Sai↗

DOE OSTI2020

Structural Analysis Approach for the Defense Programs Package 3 (DPP-3)

The Pacific Northwest National Laboratory (PNNL) is the design authority for a new Type B hazardous materials transportation package designated as the Defense Programs Package 3 (DPP-3) for the U.S. Department of Energy (DOE) National Nuclear Security Administration (NNSA). The DPP-3 has been developed using similar materials and fabrication methods employed in previous U.S. Nuclear Regulatory Commission (NRC), DOE, and NNSA certified packages. The DPP-3 design criteria are derived from the American Society of Mechanical Engineers (ASME) Boiler and Pressure Vessel Code (BPVC), NNSA guidance and NRC regulatory guides in order to safely and securely transport a variety of payloads. Final regulatory approval by the NNSA will require physical testing to demonstrate that the containment vessel (CV) remains leaktight after enduring the entire regulatory testing sequence prescribed in Title 10 of the Code of Federal Regulations Part 71 (10 CFR 71). In order to gain confidence that the DPP-3 will remain leaktight after testing, the DPP-3 has been structurally analyzed using the Finite Element Analysis (FEA) software LS-DYNA. The FEA analyses serve two general purposes: first, they aid in design and development of the package, and second, they advise as to which drop orientations are expected to cause the most damage during regulatory testing. This paper will discuss how the design criteria are incorporated into analytical techniques needed to evaluate the FEA structural simulation results for 10 CFR 71 conditions to give confidence the DPP-3 testing campaign will be successful.

Sakalaukus Jr., Peter J.↗

DOE OSTI2026

Feynman diagrams for matter wave interferometry

We introduce a new theoretical framework based on Feynman diagrams to compute phase shifts in matter wave interferometry. The method allows for analytic computation of higher order quantum corrections, beyond the traditional semi-classical approximation. These additional terms depend on the finite size of the initial matter wavefunction and/or have higher order dependence on ℏ. We apply the method to compute the response of matter wave interferometers to power law potentials and potentials with an arbitrary spatial dependence. The analytic expressions are validated by comparing to numerical simulations, and estimates are provided for the scale of the quantum corrections to the phase shift response to the gravitational field of the earth, anharmonic trapping potentials, and gravitational fields from local proof masses. We also find that for certain experimentally feasible parameters, these corrections are large enough to be measured and could lead to systematic errors if they are not mitigated. We find that to first order in a spatially dependent potential, quantum corrections vanish when the initial matter wavepacket has spherical symmetry and the potential satisfies Laplace's equation. We anticipate these quantum corrections will be especially important for trapped matter wave interferometers and for free-space matter wave interferometers in the presence of proof masses. These interferometers are becoming increasingly sensitive tools for mobile inertial sensing, gravity surveying, tests of gravity and its interplay with quantum mechanics, and searches for dark energy.

Glick, Jonah [Northwestern U.; Fermilab] (ORCID:00↗

DOE OSTI2020

Theory of the tertiary instability and the Dimits shift within a scalar model

The Dimits shift is the shift between the threshold of the drift-wave primary instability and the actual onset of turbulent transport in magnetized plasma. It is generally attributed to the suppression of turbulence by zonal flows, but developing a more detailed understanding calls for consideration of specific reduced models. The modified Terry--Horton system has been proposed by St-Onge [J. Plasma Phys. {\bf 83}, 905830504 (2017)] as a minimal model capturing the Dimits shift. Here, we use this model to develop an analytic theory of the Dimits shift and a related theory of the tertiary instability of zonal flows. We show that tertiary modes are localized near extrema of the zonal velocity $U(x)$, where $x$ is the radial coordinate. By approximating $U(x)$ with a parabola, we derive the tertiary-instability growth rate using two different methods and show that the tertiary instability is essentially the primary drift-wave instability modified by the local $U''$. Then, depending on $U''$, the tertiary instability can be suppressed or unleashed. The former corresponds to the case when zonal flows are strong enough to suppress turbulence (Dimits regime), while the latter corresponds to the case when zonal flows are unstable and turbulence develops. This understanding is different from the traditional paradigm that turbulence is controlled by the flow shear $U'$. Our analytic predictions are in agreement with direct numerical simulations of the modified Terry--Horton system.

fusion plasmas↗

DOE OSTI2024

Global Pathway Selection with Zero-RK

Global Pathway Selection (GPS) is an algorithm to effectively generates reduced (skeletal) chemistry mechanisms, which speeds up simulations and can be used as a systematic analytics tool to extract insights from complex reacting system. This release is an extension of the original code to run in parallel and to use LLNL's Zero-RK solver for fast solution of chemical problems.

Whitesides, RussellA↗

DOE OSTI2025

Global Pathway Selection with Zero-RK v0.5

Whitesides, RussellA [Lawrence Livermore National ↗

DOE OSTI2023

Attack Surface Analysis of the Digital Twins interface with Advanced Sensor and Instrumentation Interfaces: Cyber Threat Assessment and Attack Demonstration for Digital Twins in Advanced Reactor Architectures

A digital twin is a virtual representation of a physical system or object using real-time data that can predict and analyze how the system or object performs. This relatively new technology can be applied to the field of nuclear power generation, to aid in the design and development of new nuclear power plants and reduce operation costs using predictive maintenance and other data analytical methods. While there are already companies utilizing simulation software to train operators and technicians in the nuclear industry, some are now transitioning to utilizing their existing technology, software, and methods to develop digital twin solutions for the next generation of nuclear power plants, offering their services to utilities and government organizations around the world.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

DOE OSTI2025

Euler equations and the Sod shock tube problem

The Euler equations are a subset of the magnetohydrodynamic (MHD) equations in the infinitely collisional, unmagnetized limit. MHD modeling is central to many areas of plasma physics, ranging from low-temperature glow discharges to inertial confinement fusion. An important aspect of the Euler equations is their ability to describe states with discontinuities, such as shock waves. A standard benchmark test for numerical implementation of the Euler equations is the Sod shock tube. In this test, the system is initialized at rest with a pressure and density discontinuity, which results in a shock wave traveling into the low-pressure region and a rarefaction wave traveling into the high-pressure region. Starting with the presentation of the Euler equations, a numerical algorithm is presented here to solve these equations in one dimension. This is followed by an overview of the Sod shock tube problem that includes the precise initial setup and the analytic solution. Finally, the analytic solution is compared with results from numerical simulations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2026

Dynamic Multipoles in the RCS

The Rapid Cycling Synchrotron (RCS) of the Electron--Ion Collider (EIC) will accelerate polarized electrons up to 18 GeV at a repetition rate of 1 Hz. Dynamic multipoles driven by eddy currents in the vacuum chamber can affect beam dynamics and polarization preservation during the energy ramp. Building on the analytical treatment of S.Y. Lee for a ramped normal-conducting dipole magnet, we use conformal transformations to derive compact expressions for dynamic multipoles in higher-order ramped magnets with circular beam pipes. We further derive analytical frequency-domain transfer functions for both self-response and cross-response multipoles, providing compact dynamical models for their evolution during arbitrary ramp waveforms. The resulting analytical predictions are benchmarked against finite-element simulations, including the effects of practical pole geometry and beam-pipe misalignment. Finally, the expected dynamic multipoles in the RCS magnets are estimated and their implications for beam dynamics and polarization preservation are discussed.

43 PARTICLE ACCELERATORS↗

DOE OSTI2024

Hydrogen Bonding in Amorphous Indomethacin

Amorphous Indomethacin has enhanced bioavailability over its crystalline forms, yet amorphous forms can still possess a wide variety of structures. Here, Empirical Potential Structure Refinement (EPSR) has been used to provide accurate molecular models on the structure of five different amorphous Indomethacin samples, that are consistent with their high-energy X-ray diffraction patterns. It is found that the majority of molecules in amorphous Indomethacin are non-bonded or bonded to one neighboring molecule via a single hydrogen bond, in contrast to the doubly bonded dimers found in the crystalline state. The EPSR models further indicate a substantial variation in hydrogen bonding between different amorphous forms, leading to a diversity of chain structures not found in any known crystal structures. The majority of hydrogen bonds are associated with the carboxylic acid group, although a significant number of amide hydrogen bonding interactions are also found in the models. Evidence of some dipole–dipole interactions are also observed in the more structurally ordered models. The results are consistent with a distribution of Z-isomer intramolecular type conformations in the more disordered structures, that distort when stronger intermolecular hydrogen bonding occurs. The findings are supported by 1H and 2H NMR studies of the hydrogen bond dynamics in amorphous Indomethacin.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2022

Relativistic Alfvén Waves Entering Charge-starvation in the Magnetospheres of Neutron Stars

Instabilities in a neutron star can generate Alfvén waves in its magnetosphere. Propagation along the curved magnetic field lines strongly shears the wave, boosting its electric current j A . We derive an analytic expression for the evolution of the wavevector k and the growth of j A . In the strongly sheared regime, j A may exceed the maximum current j0 that can be supported by the background e ± plasma. We investigate these charge-starved waves, first using a simplified two-fluid analytic model, then with first-principles kinetic simulations. We find that the Alfvén wave is able to propagate successfully even when κ ≡ j A /j 0 $\ll$ 1. It sustains j A by compressing and advecting the plasma along the magnetic field lines with an increasing Lorentz factor, γ ≳ κ 1/2 . The simulations show how plasma instabilities lead to gradual dissipation of the wave energy. Our results suggest that an extremely high charge-starvation parameter κ ≳ 10 4 may be required in order for this mechanism to power the observed fast radio bursts (FRBs) from SGR 1935+2154. However, cosmological FRBs with much higher luminosities are unlikely to be a result of charge-starvation.

79 ASTRONOMY AND ASTROPHYSICS↗