Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “perturbation analysis”

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 523 records · Page 29
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Energy-Storing Cryogenic Carbon Capture™ for Utility and Industrial-scale Processes (Final Report)

This project leverages work completed under funding from previous U.S. Department of Energy (DOE) projects and other sources of funding around the Cryogenic Carbon Capture Process™ (CCC). The purpose of this project is to further develop a perturbation of the CCC process that stores energy called the Cryogenic Carbon Capture™ Energy Storing (CCC-ES) process. The successful completion of the tasks in this project has prepared the CCC-ES process for scaling and developed initial plans and a feasibility analysis for an engineering-scale or small-commercial scale system that can store at least 10 MWh of energy. The specific areas of work in this project are: (1) Develop a conceptual, site-specific integrated process flow diagram: Sustainable Energy Solutions (SES) developed a conceptual, site-specific plan including a process flow diagram for integrating CCC-ES at a specific location. The proposed system was based on a typical engineering-scale CCC system designed to capture 30-50 metric tonnes of CO 2 per day (TPD), with modifications and additional changes and equipment put in place to add the energy storage element to the system. The 30-50 TPD scale corresponds to the required 10 MWh of stored energy. This task includes updating simulation software for transient analysis. SES also coordinated with the Jim Bridger Power Plant in Point of Rocks, Wyoming to get site-specific data. Using these data, SES completed site-specific preliminary process flow diagrams, estimated capital and operating costs for the CCC-ES system, and used the upgraded software package to estimate the primary figures of merit. (2) Complete a technoeconomic analysis (TEA): SES completed a full-scale TEA of the CCC-ES system to meet specific success criteria. This TEA used process information from the host site selected above. (3) Technology and data gap assessment for the CCC-ES process: As part of the technology gap assessment, SES performed an assessment of the best alternative energy storage solutions including figures of merit and operational strengths and weaknesses. SES then compared the proposed system with a discussion of how it addresses the weaknesses of alternative systems, showing for example a cost of stored energy at less than $50/MWh and with a roundtrip efficiency of greater than 95%. This assessment includes technical and other risk analyses and technology gaps with research and development needs to commercialize the process by 2030. Finally, work includes a commercialization roadmap/development pathway outlining the major milestones to move to a commercially operating CCC-ES system. (4) Phase II pre-FEED project plan: SES completed a project plan for the Phase II pre-FEED analysis using the same host site as in the Phase I conceptual study.

03 NATURAL GAS↗
DOE OSTI2025

Generalized parton distributions and gravitational form factors at large momentum transfer

Within the soft collinear effective theory (SCET), we derive a factorization theorem which resums Sudakov logarithms (α s ln 2 ( –t)) n to all orders in the quark-in-quark generalized parton distribution (GPD) at large momentum transfer t, and perform a consistency check to one-loop. We show that the same Sudakov factor appears in the ‘Feynman’ contribution to the GPDs of the nucleon. Our result enables the resummation of all the large logarithms ln Q 2 and ln 2 t in exclusive processes with two hard scales Λ$^{2}_{QCD}$ ≪ |t| ≪ Q 2 . We also present a SCET power counting analysis of the Feynman contributions to the GPDs and show that the x-dependence of GPDs factorizes at large-t with controlled corrections. This in particular implies that any ratio of GPD moments such as the electromagnetic and gravitational form factors (GFF) is perturbatively calculable in this approximation. Furthermore, we identify a novel order α s power-law t-dependence in the GPD and the D-type GFF that will dominate over the standard order α$^{2}_{s}$ ‘leading twist’ asymptotic contribution in the phenomenologically relevant region of t.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2023

CSPlib: A performance portable parallel software toolkit for analyzing complex kinetic mechanisms

Computational singular perturbation (CSP) is a method to analyze dynamical systems. It targets the decoupling of fast and slow dynamics using an alternate linear expansion of the right-hand side of the governing equations based on eigenanalysis of the associated Jacobian matrix. This representation facilitates diagnostic analysis, detection and control of stiffness, and the development of simplified models. For this work, we have implemented CSP in a C++ open-source library CSPlib using the Kokkos parallel programming model to address portability across diverse heterogeneous computing platforms, i.e., multi/many-core CPUs and GPUs. We describe the CSPlib implementation and present its computational performance across different computing platforms using several test problems. Specifically, we test the CSPlib performance for a constant pressure ignition reactor model on different architectures, including IBM Power 9, Intel Xeon Skylake, and NVIDIA V100 GPU. The size of the chemical kinetic mechanism is varied in these tests. As expected, the Jacobian matrix evaluation, the eigensolution of the Jacobian matrix, and matrix inversion are the most expensive computational tasks. When considering the higher throughput characteristic of GPUs, GPUs performs better for small matrices with higher occupancy rate. CPUs gain more advantages from the higher performance of well-tuned and optimized linear algebra libraries such as OpenBLAS.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2026

From Oxo to Oxyl to Biradical: Systematic Multireference Calculations of Methane Activation at MOF Nodes

Methane C–H activation at transition-metal sites often involves electronic structures that challenge conventional single-reference electronic structure descriptions. Although Kohn–Sham density functional theory (DFT) is widely used to study catalytic trends, its reliability for reactions involving strongly correlated species remains uncertain. Here we present a systematic multireference investigation of methane activation at metal–organic framework (MOF) node catalysts across the 3d transition-metal series. We introduce an automated workflow for active space selection to enable consistent application of multireference methods, including multiconfiguration pair-density functional theory and n-electron valence state perturbation theory, to these catalytic systems. These calculations show substantial static correlation in the C–H activation reaction step and predict activation barriers that differ from DFT by 30–70 kJ mol–1, with DFT often qualitatively disagreeing in barrier height trends across transition metals. Analysis of multireference wave functions shows that reactivity is governed by the electronic structure of the M–O moiety along a continuum from metal–oxo to oxyl radical and O biradical character. Increased oxygen-centered spin density and weakened M–O bonding are identified as descriptors of catalytic activity which correlate with lower activation barriers.

Wardzala, Jacob↗
DOE OSTI2020

Long-range quadrupole electron-phonon interaction from first principles

Lattice vibrations in materials induce perturbations on the electron dynamics in the form of long-range (dipole and quadrupole) and short-range (octopole and higher) potentials. The dipole Fröhlich term can be included in current first-principles electron-phonon (e-ph) calculations and is present only in polar materials. The quadrupole e-ph interaction is present in both polar and nonpolar materials, but currently it cannot be computed from first principles. In this work, we show an approach to compute the quadrupole e-ph interaction and include it in ab initio calculations of e-ph matrix elements. The accuracy of the approach is demonstrated by comparing with direct density functional perturbation theory calculations. We apply our method to silicon as a case of a nonpolar semiconductor and tetragonal PbTiO 3 as a case of a polar piezoelectric material. In both materials we find that the quadrupole term strongly impacts the e-ph matrix elements. Analysis of e-ph interactions for different phonon modes reveals that the quadrupole term mainly affects optical modes in silicon and acoustic modes in PbTiO 3 , although the quadrupole term is needed for all modes to achieve quantitative accuracy. The effect of the quadrupole e-ph interaction on electron scattering processes and transport is shown to be important. Our approach enables accurate studies of e-ph interactions in broad classes of nonpolar, polar, and piezoelectric materials.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Gluon distributions and their applications to Ioffe-time distributions

We investigate unpolarized and polarized gluon distributions and their applications to the Ioffe-time distributions, which are related to lattice QCD calculations of parton distribution functions. Guided by the counting rules based on the perturbative QCD at large momentum fraction x and the color coherence of gluon couplings at small x, we parametrize gluon distributions in the helicity basis. By fitting the unpolarized gluon distribution, the inferred polarized gluon distribution from our parametrization agrees with the one from global analysis. A simultaneous fit to both unpolarized and polarized gluon distributions is also performed to explore the model uncertainty. The agreement with the global analysis supports the (1 - x) power suppression of the helicity-antialigned distribution relative to the helicity-aligned distribution. The corresponding Ioffe-time distributions and their asymptotic expansions are calculated from the gluon distributions. Our results of the Ioffe-time distributions can provide guidance to the extrapolation of lattice QCD data to the region lacking precise gluonic matrix elements. Therefore, they can help regulate the ill-posed inverse problem associated with extracting the gluon distributions from discrete data from first-principle calculations, which are available in a limited range of the nucleon momentum and the spatial separation between the gluonic currents. Given various limitations in obtaining lattice QCD data at large Ioffe time, phenomenological approaches can provide important complementary information to extract the gluon distributions in the entire momentum fraction region, especially at small x. The possibility of investigating higher-twist effects and other systematic uncertainties in the contemporary first-principle calculations of parton distributions from phenomenologically well-determined Ioffe-time distributions in the large Ioffe-time region is also discussed.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2025

Visolis Microbial Chemical Intermediate Library Screening (CRADA Final Report)

Visolis is developing a commercial process for the bioproduction of a chemical intermediate, towards the derivative production of a portfolio of bio-based chemicals with large application potential, from drop-in commodity chemicals, to innovative specialty chemical offerings, to materials for multiple end markets and industries, which will ultimately benefit consumers and the public. Visolis already has several variants of a microbe that produce the chemical intermediate at distinct levels (i.e. high, medium, and low). They are currently pursuing transcriptomics analysis for some of these strains, towards a better understanding of the biology behind what makes certain strains perform better than others. There exist genome-scale library approaches for the microbe that could be used to screen gene disruption, overexpression, or repression candidates for perturbations to the production of the chemical intermediate. These approaches often reveal opportunities for further production improvement that are not otherwise accessible using hypothesis-driven metabolic engineering approaches. However, Visolis, while it could generate or procure such genome-scale libraries, does not have the automated strain engineering workflows required to screen thousands of variants, obtaining a production phenotype and a genotype for each. The purpose of this collaboration is to use LBNL and SNL capabilities to enable Visolis to effectively screen thousands of genome-scale library strain variants for phenotype/genotype relationships that will complement Visolis’ transcriptomic investigations.

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

Replacing Pyridine with Pyrazine in Molecular Cobalt Catalysts: Effects on Electrochemical Properties and Aqueous H2 Generation

Four new molecular Co(II)tetrapyridyl complexes were synthesized and evaluated for their activity as catalysts for proton reduction in aqueous environments. The pyridine groups around the macrocycle were substituted for either one or two pyrazine groups. Single crystal X-ray analysis shows that the pyrazine groups have minimal impact on the Co(II)–N bond lengths and molecular geometry in general. X-band EPR spectroscopy confirms the Co(II) oxidation state and the electronic environment of the Co(II) center are only very slightly perturbed by the substitution of pyrazine groups around the macrocycle. The substitution of pyrazine groups has a substantial impact on the observed metal- and ligand-centered reduction potentials as well as the overall H2 catalytic activity in a multimolecular system using the [Ru(2,2′-bipyridine)3]Cl2 photosensitizer and ascorbic acid as a sacrificial electron donor. The results reveal interesting trends between the H2 catalytic activity for each catalyst and the driving force for electron transfer between either the reduced photosensitizer to catalyst step or the catalyst to proton reduction step. The work presented here showcases how even the difference of a single atom in a molecular catalyst can have an important impact on activity and suggests a pathway to optimize the photocatalytic activity and stability of molecular systems.

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

Analysis of SCALE Criticality and Sensitivity Calculations for Reflected HEU Cylinders [Slides]

The SCALE code package offers multiple nuclear data libraries supporting Monte Carlo transport, with sensitivity and uncertainty methods derived from MC transport solutions. Several libraries are multigroup, which introduce bias differing by system. Previous work has shown poor S/U results in several reflector materials: ICSBEP benchmark HMF-084 was selected to analyze biases and S/U method applicability to a variety of reflectors. Prior and ongoing work found inaccuracies in CSAS and TSUNAMI results, which were further investigated utilizing the HEU-MET-FAST-084 ICSBEP critical benchmark, chosen for its geometrical simplicity and variety of reflector materials. Perturbation of reflector thickness across various reflector materials allowed for an assortment of materials is to be tested swiftly for each sequence and method. Observation was an increasing bias of MG $k_{eff}$ relative to CE, in both direction and magnitude. IFP produced extremely reliable results. >85% of CLUTCH cases were found satisfactory.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2026

Jet fragmentation function and groomed substructure of bottom quark jets in proton-proton collisions at 5.02 TeV

A measurement of the substructure of bottom quark jets (b jets) in proton-proton (pp) collisions is presented. The measurement uses data collected in pp collisions at $\sqrt{s}=5.02$ TeV, with a low number of simultaneous interactions per bunch crossing, recorded by the CMS experiment in 2017, corresponding to an integrated luminosity of 301 pb −1 . An algorithm to identify and cluster the charged decay daughters of b hadrons is developed for this analysis, which facilitates the exposure of the gluon radiation pattern of b jets using iterative Cambridge-Aachen declustering. The soft-drop-groomed jet radius, R g , and momentum balance, z g , of b quark jets are presented. These observables can be used to test perturbative quantum chromodynamics predictions that account for mass effects. Because the b hadron is partially reconstructed from its charged decay daughters, only charged particles are used for the jet substructure studies. In addition, a jet fragmentation function, z b,ch , is measured, which is defined as the distribution of the ratio of the transverse momentum (p T ) of the partially reconstructed b hadron with respect to the charged-particle component of the jet p T . The substructure variable distributions are unfolded to the charged-particle level. The b jet substructure is compared to the substructure of jets in an inclusive jet sample that is dominated by light-quark and gluon jets in order to assess the role of the b quark mass. A strong suppression of emissions at small R g values is observed for b jets when compared to inclusive jets, consistent with the dead-cone effect. The measurement is also compared with theoretical predictions from Monte Carlo event generators. This is the first substructure measurement of b jets that clusters together the b hadron decay daughters independent of the b hadron species and decay channel.

boosted jets↗
DOE OSTI2026

Clustering analysis of medium-band selected high-redshift galaxies

Next-generation large-scale structure spectroscopic surveys will probe cosmology at high redshifts (2.3 < z < 3.5), relying on abundant galaxy tracers such as Lyα emitters (LAEs) and Lyman break galaxies (LBGs). Medium-band photometry has emerged as a potential technique for efficiently selecting these high-redshift galaxies. In this work, we present clustering analysis of medium-band selected galaxies at high redshift, utilizing photometric data from the Intermediate Band Imaging Survey (IBIS) and spectroscopic data from the Dark Energy Spectroscopic Instrument (DESI). We interpret the clustering of such samples using both Halo Occupation Distribution (HOD) modeling and a perturbation theory description of large-scale structure. Our modeling indicates that the current target sample is composed from an overlapping mixture of LAEs and LBGs with emission lines. Despite differences in target selection, we find that the clustering properties are consistent with previous studies, with correlation lengths r 0 ≃ 3-4 h -1 Mpc and a linear bias of b ∼ 1.8-2.5. Finally, we discuss the simulation requirements implied by these measurements and demonstrate that the properties of the samples would make them excellent targets to enhance our understanding of the high-z universe.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Temperature Effects on Droplet Oscillation Decay with Application to Fuel Property Measurement

Observation of oscillation decay in droplets has been shown to be an effective approach in determining physical properties such as viscosity and surface tension for emerging biofuels using μl quantities. Herein this work extends the approach to higher temperatures relevant to fuel injection conditions for internal combustion engines. Experiments are conducted that use high-resolution strobed imaging of moving, heated, μm-sized, fuel droplets to capture shape oscillation decay through image processing and analysis. Two fuels are investigated, iso-butanol and a primary reference fuel (PRF 84), which is a mixture of iso-octane and n-heptane. A piezoelectric droplet generator is used to generate a continuous train of single droplets, which are given an initial perturbation and subsequently undergo damped oscillations, captured using strobed imaging. Surface tension and viscosity calculated at four different temperatures ranging from 30°C-55°C using the frequency and decay time associated with the fundamental mode are found to be within ~ 10% of reference values obtained from the literature. Complementary numerical simulations are performed that utilize a volume-of-fluid approach to track transient droplet oscillation phenomena along with heat and mass transfer in a 2D axisymmetric domain. Simulations, where droplet size and temperature can be independently varied, capture an expected decrease in oscillation frequency and increase in decay time, with increase in fuel temperature. The simulations are further used to study the relative contribution to deviations in surface tension and viscosity predictions due to heat and mass transfer from the droplet, as well as viscous effects violating the inviscid flow assumption in the droplet oscillation theory. Mass loss effects are found to be negligible. Temperature change due to heat transfer is found to have the next highest sensitivity, particularly for the more volatile PRF 84, which has a lower viscosity and associated Ohnesorge number. For isobutanol, which has a higher viscosity and Ohnesorge number, viscous effects contribute the most to deviation in fuel property predictions.

09 BIOMASS FUELS↗
DOE OSTI2022

Resistive wall tearing mode disruptions in DIII-D and ITER tokamaks

Disruptions are a serious problem in tokamaks, in which thermal and magnetic energy confinement is lost. This paper uses data from the DIII-D experiment, theory, and simulations to demonstrate that resistive wall tearing modes (RWTMs) produce the thermal quench (TQ) in a typical locked mode shot. Analysis of the linear RWTM dispersion relation shows the parameter dependence of the growth rate, particularly on the resistive wall time. Linear simulations of the locked mode equilibrium show that it is unstable with a resistive wall and stable with an ideally conducting wall. Nonlinear simulations demonstrate that the RWTM grows to sufficient amplitude to cause a complete thermal quench. The RWTM growth time is proportional to the thermal quench time. Here, the nonlinearly saturated RWTM magnetic perturbation amplitude agrees with experimental measurements. The onset condition is that the q = 2 rational surface is sufficiently close to the resistive wall. Collectively, this identifies the RWTM as the cause of the TQ. In ITER, RWTMs will produce long TQ times compared to present-day experiments. ITER disruptions may be significantly more benign than previously predicted.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Development and application of two-step uncertainty propagation and sensitivity analysis methodology for fast reactor safety analysis

Uncertainty quantification (UQ) in nuclear reactors for transients is directly linked with safety assessment through the cross-sections uncertainties, provided as a covariance matrix, which are propagated through the reactor system to output of interest pertaining to reactor safety, such as peak temperatures in fuel/clad/coolant. Using a two-step approach, uncertainties are first quantified and propagated from basic input variables (such as reaction cross-sections) to intermediate quantities (such as reactivity feedback coefficients) through lattice level calculations. Uncertainties of intermediate quantities (from the first step) are then propagated through the system transient calculations, in the second step, to obtain uncertainties on reactor safety output parameters of interest. The scope of this work consists of Uncertainty Quantification & Propagation of nuclear data uncertainties that are highly correlated through unprotected transient overpower and unprotected loss of flow to assess their impact on core safety parameters. This two-step approach in the presence of covariance renders the sensitivity analysis very challenging. In fact, usually the sensitivity analysis is restricted to each step, which limits its application since the sensitivities between the system output quantities and the basic input variables are difficult to obtain. Here, in this work, we address this issue by proposing a simple, general methodology to combine the sensitivity indices obtained in each step by assuming the model behavior being linear. For the first step Generalized Perturbation theory based indices are used while in the second step the recently studied Johnson indices. The uncertainty quantification and sensitivity methodologies discussed here are demonstrated on a generic LFR design which is based on the 500 MWth demonstration Lead-cooled fast reactor (DLFR) using oxide fuel, developed by Westinghouse Electric Company (WEC).

42 - ENGINEERING↗
DOE OSTI2024

Fermions at finite density in the path integral approach

Abstract We study relativistic fermionic systems in 3 + 1 spacetime dimensions at finite chemical potential and zero temperature, from a path-integral point of view. We show how to properly account for theiεterm that projects on the finite density ground state, and compute the path integral analytically for free fermions in homogeneous external backgrounds, using complex analysis techniques. As an application, we show that the U(1) symmetry is always linearly realized for free fermions at finite charge density, differently from scalars. We study various aspects of finite density QED in a homogeneous magnetic background. We compute the free energy density, non-perturbatively in the electromagnetic coupling and the external magnetic field, obtaining the finite density generalization of classic results of Euler-Heisenberg and Schwinger. We also obtain analytically the magnetic susceptibility of a relativistic Fermi gas at finite density, reproducing the de Haas-van Alphen effect. Finally, we consider a (generalized) Gross-Neveu model forNinteracting fermions at finite density. We compute its non-perturbative effective potential in the large-Nlimit, and discuss the fate of the U(1) vector and$$ {\mathbb{Z}}_2^A $$ ℤ 2 A axial symmetries.

Physics↗
DOE OSTI2025

Single-Determinant Ground State in Ce 4+ Imidophosphorane Complexes

X-ray spectroscopy techniques are critical in the electronic structure analysis of high-valent lanthanides. The interpretation of multipeaked features at the lanthanide L 3 -edge has remained a challenging question, as it is observed across a range of material classes. A series of structurally related Ce 4+ complexes were prepared to probe the potential ligand field perturbation of the ground state within the series. The tuning of relative 4f and ligand orbital energies in homoleptic and heteroleptic tetravalent Ce imidophosphorane complexes is achieved through ligand derivatization and is clearly demonstrated by UV−vis spectroscopy and electrochemically measured redox potentials. However, Ce L 3 - edge high-energy resolution fluorescence-detected (HERFD) X-ray absorption near-edge structure (XANES) spectra present features at consistent numbers and energies across the range of complexes. Resonant inelastic X-ray scattering (RIXS) is employed to visualize the observed features in the HERFD-XANES spectra. Large complete active space configuration interaction singles and doubles (CASCISD) calculations demonstrate that the ground-state wave function of all complexes can be described employing a single determinant. As a result, the multielectron feature at the L 3 -edge observed in this study for the Ce 4+ imidophosphorane complexes is described as excited-state multiconfigurational behavior that is independent of ligand variation, in systems where the ground state is described using the simple single determinant wave function.

Lanthanides↗