Search NASA⌕ Search

SEARCH · Search NASA

Results for “time evolution operators”

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 631 records · Page 35
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Simulation of Continental Shallow Cumulus Populations Using an Observation-Constrained Cloud-System Resolving Model

Continental shallow cumulus (ShCu) clouds observed on 30 August 2016 during the Holistic Interactions of Shallow Clouds, Aerosols, and Land-Ecosystems (HI-SCALE) field campaign are simulated by using an observation-constrained cloud-system resolving model. On this day, ShCu forms over Oklahoma and southern Kansas and some of these clouds transition to deeper, precipitating convection during the afternoon. We apply a four-dimensional ensemble-variational (4DEnVar) hybrid technique in the Community Gridpoint Statistical Interpolation (GSI) system to assimilate operational data sets and unique boundary layer measurements including a Raman lidar, radar wind profilers, radiosondes, and surface stations collected by the U.S. Department of Energy's (DOE) Atmospheric Radiation Measurement (ARM) Southern Great Plains (SGP) atmospheric observatory into the Weather Research and Forecasting (WRF) model to ascertain how improved environmental conditions can influence forecasts of ShCu populations and the transition to deeper convection. Independent observations from aircraft, satellite, as well as ARM's remote sensors are used to evaluate model performance in different aspects. Several model experiments are conducted to identify the impact of data assimilation (DA) on the prediction of clouds evolution. The analyses indicate that ShCu populations are more accurately reproduced after DA in terms of cloud initiation time and cloud base height, which can be attributed to an improved representation of the ambient meteorological conditions and the convective boundary layer. Extending the assimilation to 18 UTC (local noon) also improved the simulation of shallow-to-deep transitions of convective clouds.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2015

Optimisation of sewage sludge anaerobic digestion through co-digestion with OFMSW: Effect of collection system and particle size

Highlights: • Methane production rate increased between 56% and 208% during OFMSW–SS codigestion. • The OFMSW particle size reduction from 20 to 8 mm did not affect the methane yield. • OFMSW–SS codigestion promoted β-oxidation and acetoclastic methanogenic activity. • The evolution of specific activity was a feasible tool to control the process. - Abstract: The effect of organic fraction of municipal solid waste (OFMSW) loading rate and particulate size on the sewage sludge (SS) mesophilic anaerobic co-digestion was assessed in continuous stirred tank reactor at hydraulic retention time of 20 days. The SS–OFMSW mixture composed by 54% of the volatile solids fed (inlet-VS), at OLR of 3.1 kg{sub COD} m{sup −3} d{sup −1} (1.9 kg{sub VS} m{sup −3} d{sup −1}), showed the highest increment on the volumetric methane production and yield of +200% and +59% respectively, under stable conditions. The effect of particulate size was assessed with the same mixture and same operational conditions but reducing the OFMSW particulate size from 20 mm to 8 mm with the aim to improve the hydrolysis step, but the results showed any influence in the OFMSW particulate size range analysed. In addition, specific biomass activity was assessed at the end of each co-digestion period. Results showed that OFMSW promoted β-oxidation syntrophic acetogens and the acetoclastic methanogens activity; although the last increase of the OFMSW percentage (from 47% to 54% inlet-VS) affected negatively the specific substrate activity, but not inhibitory effect was observed. Therefore, the results obtained in the continuous experiment could be related with some inhibitory or toxic effect and not due to hydrolysis limitation. The specific biomass activity test was demonstrated to be an interesting tool to evaluate and control the co-digestion process, especially when conventional parameters did not explain the behaviour of the biological system.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2022

Planar laser induced fluorescence mapping of a carbon laser produced plasma

We present measurements of ion velocity distribution profiles obtained by laser induced fluorescence (LIF) on an explosive laser produced plasma. The spatiotemporal evolution of the resulting carbon ion velocity distribution was mapped by scanning through the Doppler-shifted absorption wavelengths using a tunable, diode-pumped laser. The acquisition of these data was facilitated by the high repetition rate capability of the ablation laser (1 Hz), which allowed for the accumulation of thousands of laser shots in short experimental times. By varying the intensity of the LIF beam, we were able to explore the effects of fluorescence power against the laser irradiance in the context of evaluating the saturation vs the non-saturation regime. The small size of the LIF beam led to high spatial resolution of the measurement compared to other ion velocity distribution measurement techniques, while the fast-gate operation mode of the camera detector enabled the measurement of the relevant electron transitions.

Instruments & Instrumentation↗
DOE OSTI2025

Experimental demonstration of real-time electron temperature profile control in DIII-D

Future tokamak reactor operation will require the ability to maintain a given plasma scenario for extended periods of time. This will necessitate the capability to react to changes in the plasma state and return the plasma to the target scenario; the principal method to achieve this is through feedback control. Thus, it is necessary to develop and test feedback controllers for the plasma profiles that define a target scenario. In this work, a feedback controller for the electron temperature (Te) profile is tested experimentally in DIII-D. This experiment relied on the ability to ascertain the electron temperature profile in real time, which was achieved using an observer algorithm. The observer relies on both diagnostic data and a predictive model of the electron temperature profile evolution; this predictive model includes contributions from neural network surrogate models. Because of these dependencies, a number of capabilities needed to be added to the real-time PCS for DIII-D in order to support the Te profile control experiment. The neural network surrogates needed to be integrated into the PCS to be called in real time. An observer algorithm for the Te profile needed to be added and connected to the Thomson scattering system to allow access to the current state of the profile in real time. When tested, the observer was shown to produce Te profiles that are consistent with the shape of the Thomson scattering data while rejecting much of the noise in the diagnostic data. Finally, the controller itself was tested in real time. This experiment showed that the controller is capable of tracking the electron temperature target at locations across the spatial profile.

Morosohk, Shira [Oak Ridge Associated Universities↗
DOE OSTI2021

Adaptive machine learning for time-varying systems: low dimensional latent space tuning

Machine learning (ML) tools such as encoder-decoder convolutional neural networks (CNN) can represent incredibly complex nonlinear functions which map between combinations of images and scalars. For example, CNNs can be used to map combinations of accelerator parameters and images which are 2D projections of the 6D phase space distributions of charged particle beams as they are transported between various particle accelerator locations. Despite their strengths, applying ML to time-varying systems, or systems with shifting distributions, is an open problem, especially for large systems for which collecting new data for re-training is impractical or interrupts operations. Particle accelerators are one example of large time-varying systems for which collecting detailed training data requires lengthy dedicated beam measurements which may no longer be available during regular operations. We present a novel method of adaptive ML for time-varying systems. Our approach is to map very high (N ≈ 100k) dimensional inputs (a combination of scalar parameters and images) into the low dimensional (N ≈ 2) latent space at the output of the encoder section of an encoder-decoder CNN. We then actively tune the low dimensional latent space-based representation of complex system dynamics by the addition of an adaptively tuned feedback vector directly before the decoder sections builds back up to our image-based high-dimensional phase space density representations. This method allows us to learn correlations within and to quickly tune the characteristics of incredibly large parameter space systems and to track their evolution in real time based on feedback without massive new data sets for re-training. We demonstrate that our method can accurately predict and track the phase space of charged particle beams at various locations in a particle accelerator by adaptively adjusting in real-time while the unknown input beam distribution of the accelerator is changing in shape, charge, and offset and while the RF system of the accelerator itself is also changing in an unpredictable way. For FACET-II we demonstrate that such an approach has the potential to use transverse deflecting cavity and energy spread spectrum beam measurements to accurately predict 2D projections of the 6D phase space of the electron beam at the plasma wakefield acceleration interaction point where such diagnostics are unavailable.

47 OTHER INSTRUMENTATION↗
DOE OSTI2016

Views on Long-Term Safety of Deep Borehole Disposal - 16657

The Swedish Nuclear Fuel and Waste Management Co., SKB, has since the end of the 1980s evaluated Deep Borehole Disposal (DBD), of spent nuclear fuel (SNF). This has been done in parallel with the development of a mined repository concept based on the KBS-3V method. The most recent evaluation involved a broad comparison between the KBS-3V method and DBD. The comparison highlights differences between the two methods in the whole chain of the SNF handling including premises for siting, licensing, repository construction, safety during the operational phase, post-closure safety, physical protection, nuclear safeguards, retrieval of the deposited SNF, time planning, need for development, costs and nontechnical project risks. The ambition has been to make an as objective comparison between the two methods as possible, bearing in mind that there are big differences between the two concepts in terms of technical maturity as well as in terms of both quantity and quality of available data. While the primary safety function of the KBS-3V method is long-term containment in a corrosion resistant copper canister, the safety of DBD relies on stagnant groundwater at depth. It is argued in the aforementioned study that this groundwater stagnancy may be jeopardised because of heat and gas evolution in the disposal boreholes, and that deformation of the boreholes make efficient sealing of the deposition boreholes difficult. SKB has taken the position that the DBD method is associated with too many question marks to make it interesting as an alternative method in the Swedish programme. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2019

Real-time capable modeling of neutral beam injection on NSTX-U using neural networks

A new model of heating, current drive, torque and other effects of neutral beam injection on NSTX-U that uses neural networks has been developed. The model has been trained and tested on the results of the Monte Carlo code NUBEAM for the database of experimental discharges taken during the first operational campaign of NSTX-U. By projecting flux surface quantities onto empirically derived basis functions, the model is able to efficiently and accurately reproduce the behavior of both scalars, like the total neutron rate and shine through, and profiles, like beam current drive and heating. With this being said, the model has been tested on the NSTX-U real-time computer, demonstrating a rapid execution time orders of magnitude faster than the Monte Carlo code that is well suited for the iterative calculations needed to interpret experimental results, optimization during scenario development activities, and real-time plasma control applications. Simulation results of a proposed design for a nonlinear observer that embeds the neural network calculations to estimate the poloidal flux profile evolution, as well as Z eff and fast ion diffusivity, are presented.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2018

Three dimensional optimization for subsonic axial turbines operating at high unsteady inlet Mach number

In this paper, a three dimensional optimization strategy is presented to improve the performance of a subsonic axial turbine exposed to high Mach number unsteady inflows encountered in pressure gain combustors. The optimization aims at altering the stator hub and shroud end wells. A multi-objective differential evolution optimizer was used to enhance specific work and efficiency. The optimization routine consisted of a multi-step approach where in a first step an optimization was carried out where the turbine was subjected to a steady high Mach inflow. The turbine stage was parametrized via an in-house turbine design tool based on Bezier curves. Seven design parameters were selected to contour the stator end walls and rotor hub end wall. The parametrized turbine stage was meshed and the steady Reynolds-Averaged Navier-Stokes equations were solved via the commercial solver CFD++ and calculation time of each stage was around 10 hours. Results of the steady optimization identified three major design parameters of stator end walls, which were taken as inputs for the next phase unsteady optimization. For the unsteady simulations, a large inlet fluctuation was imposed and the turbine was modeled with the mixing plane approach. Total calculation time of each individual was 48 hours. An efficiency improvement of 10% from the steady optimization and 13% from the unsteady optimization was found compared to the baseline case. Finally, full unsteady simulations with two stator and three rotor passages were performed on selected steady and unsteady optimized configurations.

03 NATURAL GAS↗
DOE OSTI2026

Microchannel-based Membrane-less Extraction of Li from Unconventional Lithium Sources & the Separation of REE

This final report provides an overview of the Project's entire duration, covering July 1, 2021 to December 31, 2023. It primarily focuses on the achievements, technological developments, and unique challenges the team faced while working on separating and extracting Lithium from produced waters. The project's primary aim was to create an integrated, high-throughput, membrane-less, and modular microfluidic platform that could extract Lithium from unconventional sources. We have successfully met all goals and milestones envisioned in the SOPO document. The most critical primary milestones, including the Go-No-Go milestone (refer to the Gantt chart in the Appendices), were successfully accomplished. We demonstrated phase separation (>90%) and extraction (>85%) performance in the MPSE using synthetic, and representative produced water composition feed at 50 ml/min total flow through MPSE 36. We have also performed a parametric study of the MPSE operations, beyond the scope of SOPO, exploring operating conditions of current and broader interest. The extended investigation of operational parameters is concurrent with our efforts to seek further development of the MPSE technology beyond the scope of the Project. Along these lines of development, we have made efforts to be responsive to DOE calls for technological developments of other types of resources (beyond PW) for the recovery of Critical Materials and higher TRL development (beyond TRL 4). During the work on this Project, we developed and implemented three innovative technical approaches that emerged from our efforts to successfully meet the Project milestones. The innovative & original technical approaches developed and implemented in this Project are now the contributions to process engineering that could be clearly credited to the Project. First, Convergent Design Approach is a comprehensive feedforward & feedback loop of four design phases: i) design for functionality, ii) design for manufacturing, iii) design for sustainability, and iv) design for market. Next was Process Intensification. A major aim of this Project was to create an innovative phase separation & extraction microscale-based technology for Li separation – thus the words microchannel-based in the Project title. A microscale-based technology is intrinsically in the center of the Process Intensification domain as defined by its unique principles. Therefore, Process Intensification was implicitly envisioned in the Project’s SOPO. Lastly, Time Scale Analysis is a novel tool for discovering the needs and directions of Process Intensification implementations in any process technology. This Project is fully credited for developing and implementing the three novel technical approaches mentioned above. These are general contributions to process engineering that emerged from this Project. Beyond the original SOPO scope, the OSU-U.Pitt research group utilized a Convergent Design methodology, integrating first-principles mathematical modeling with experimental validation on the Minimum Development Vehicle. By creating these Digital Twins, the team rapidly assessed manufacturing iterations to support TEA analysis. This framework further enabled the development of advanced Surface Modification Techniques, where hydrophobic and oleophobic coating strategies were optimized via Digital Twin tools and validated through rigorous 100-hour longevity testing. TEA Analysis: The closing efforts of this Project were focused on the TEA analysis. TEA analysis had two primary functions: i) enabling critical assessments of design variations withing 10 the Concurrent Design Approach, thus enabling evolution of the MPSE design to reach faster- better-cheaper alternatives; and ii) to create a bridge between the accomplishments of this Project and future projects of higher TRL, beyond TRL 6 level. It is important to note that the TEA model created in the Project stirred the technological solutions for the recovery of critical materials toward a vision of a very profitable modular plant that has unique zero-waste water discharge signature. More importantly, thanks to our experimental performance data and conservative assumptions, the TEA model predicts minimal technological and investment risks. Low cost of a modular unit of a nominal capacity of [1000 tons of Li 2 CO 3 /year] positions the MPSE based technology within the reach of community investors, thus offering a paradigm shift in the development of critical technologies. The project successfully navigated two primary challenges: solvent selection and manufacturing adaptation. Restricted by the SOPO to existing literature for lithium recovery, the team identified a critical need for a "material excellence program" to develop next-generation solvents, eventually concluding with a preliminary investigation into promising Ionic Liquids (ILs). Simultaneously, COVID-19 supply chain disruptions forced a pivot from traditional manufacturing to advanced additive methods at ATAMI-OSU. By transitioning from stainless steel to 3D-printed polymer substrates, the team achieved a transformative three-order-of- magnitude reduction in manufacturing costs and compressed prototyping timelines from several months to just two days. The MPSE technology offers significant energy, environmental, and economic advantages by overcoming the traditional bottlenecks of phase-separation hardware and contactor size. Unlike conventional mixer-settlers or membrane-based systems, MPSE operates without moving parts or fouling-prone membranes, achieving robust performance even with challenging, viscous, or particulate-heavy feeds. Key performance metrics include an energy intensity reduction of 5–50x (3–40 kJ/m 3 ) compared to incumbent technologies and a dramatic reduction of processing time to under 60 seconds, which drastically reduces the physical plant footprint. These technical efficiencies translate into superior economic outcomes; for a 100 t/year Li 2 CO 3 facility, implementing MPSE is projected to nearly halve contactor CAPEX (from $\$$6.08M to $\$$3.01M) and significantly increase the project's Net Present Value (NPV), derisking new investment and enabling distributed critical-mineral processing configurations. The commercialization of MPSE technology is being spearheaded by Vigsur Dynamics Inc., which has adopted a structured, parallel approach to technical and business development since its formation in January 2026. Following extensive customer discovery and engagement with the Oregon State University accelerator, Vigsur Dynamics is working to establish a business model that transitions from pilot demonstrations to modular hardware sales, ultimately aiming for a "build-own-operate" service strategy. Current technical milestones—including 100 hours of continuous operation, superior energy efficiency, and successful 6-unit modular scale-up— provide a foundation for this transition. Backed by ongoing IP licensing and a growing network of industrial and venture advisors, the company is actively de-risking the platform to replace conventional mixer-settler systems in the critical minerals market.

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

Asymptotic analysis of discrete schemes for non-equilibrium radiation diffusion

Motivated by providing well-behaved fully discrete schemes in practice, this paper extends the asymptotic analysis on time integration methods for non-equilibrium radiation diffusion in [2] to space discretizations. Therein studies were carried out on a two-temperature model with Larsen's flux-limited diffusion operator, both the implicitly balanced (IB) and linearly implicit (LI) methods were shown asymptotic-preserving. In this paper, we focus on asymptotic analysis for space discrete schemes in dimensions one and two. First, in construction of the schemes, in contrast to traditional first-order approximations, asymmetric second-order accurate spatial approximations are devised for flux-limiters on boundary, and discrete schemes with second-order accuracy on global spatial domain are acquired consequently. Then by employing formal asymptotic analysis, the first-order asymptotic-preserving property for these schemes and furthermore for the fully discrete schemes is shown. Finally, with the help of manufactured solutions, numerical tests are performed, which demonstrate quantitatively the fully discrete schemes with IB time evolution indeed have the accuracy and asymptotic convergence as theory predicts, hence are well qualified for both non-equilibrium and equilibrium radiation diffusion. - Highlights: • Provide AP fully discrete schemes for non-equilibrium radiation diffusion. • Propose second order accurate schemes by asymmetric approach for boundary flux-limiter. • Show first order AP property of spatially and fully discrete schemes with IB evolution. • Devise subtle artificial solutions; verify accuracy and AP property quantitatively. • Ideas can be generalized to 3-dimensional problems and higher order implicit schemes.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

Scalable Implicit Solvers with Dynamic Mesh Adaptation for a Relativistic Drift-Kinetic Fokker–Planck–Boltzmann Model

In this work we consider a relativistic drift-kinetic model for runaway electrons along with a Fokker–Planck operator for small-angle Coulomb collisions, a radiation damping operator, and a secondary knock-on (Boltzmann) collision source. Here, we develop a new scalable fully implicit solver utilizing finite volume and conservative finite difference schemes and dynamic mesh adaptivity. A new data management framework in the PETSc library based on the p4est library is developed to enable simulations with dynamic adaptive mesh refinement (AMR), distributed memory parallelization, and dynamic load balancing of computational work. This framework and the runaway electron solver building on the framework are able to dynamically capture both bulk Maxwellian at the low-energy region and a runaway tail at the high-energy region. To effectively capture features via the AMR algorithm, a new AMR indicator prediction strategy is proposed that is performed alongside the implicit time evolution of the solution. This strategy is complemented by the introduction of computationally cheap feature-based AMR indicators that are analyzed theoretically. Numerical results quantify the advantages of the prediction strategy in better capturing features compared with nonpredictive strategies; and we demonstrate trade-offs regarding computational costs. The robustness with respect to model parameters, algorithmic scalability, and parallel scalability are demonstrated through several benchmark problems including manufactured solutions and solutions of different physics models. We focus on demonstrating the advantages of using implicit time stepping and AMR for runaway electron simulations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2019

A Strategy for Increasing the Efficiency of the Oxygen Reduction Reaction in Mn-Doped Cobalt Ferrites

Alkaline fuel cells have drawn increasing attention as next-generation energy-conversion devices for electrical vehicles, since high pH enables the use of non-precious-metal catalysts. Herein, we report on a family of rationally designed Mn-doped cobalt ferrite (MCF) spinel nanocrystals, with an optimal composition Mn 0.8 (CoFe 2 ) 0.73 O 4 (MCF-0.8), that are effective electrocatalysts for the oxygen reduction reaction. MCF-0.8 exhibits a half-wave potential (E 1/2 ) of 0.89 V vs RHE in 1 M NaOH, only 0.02 V less than that of commercial Pt/C under identical testing conditions and, to the best of our knowledge, one of the highest recorded values in the literature. Moreover, MCF-0.8 exhibits remarkable durability (ΔE 1/2 = 0.014 V) after 10 000 electrochemical cycles. In situ X-ray absorption spectroscopy (XAS) reveals that the superior performance of the trimetallic MCF-0.8 originates from the synergistic catalytic effect of Mn and Co, while Fe helps preserve the spinel structure during cycling. We employed in situ XAS to track the evolution of the oxidation states and the metal–oxygen distances not only under constant applied potentials (steady state) but also during dynamic cyclic voltammetry (CV) (nonsteady state). The periodic conversion between Mn(III, IV)/Co(III) and Mn(II, III)/Co(II) as well as the essentially constant oxidation state of Fe during the CV suggests collaboration efforts among Mn, Co, and Fe. Mn and Co serve as the synergistic coactive sites to catalyze the oxygen reduction, apparently resulting in the observed high activity, while Fe works to maintain the integrity of the spinel structure, likely contributing to the remarkable durability of the catalyst. Furthermore, these findings provide a mechanistic understanding of the electrocatalytic processes of trimetallic oxides under real-time fuel cell operating conditions. This approach provides a new strategy to design high-performance non-precious-metal electrocatalysts for alkaline fuel cells.

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

Snow ALbedo eVOlution (SALVO) Campaign Spectral Albedo and Related Measurements from April - June, 2024 in Utqiagivk, AK

A field-portable spectroradiometer, referred to herein as an ‘ASD’, was used to make spatially-distributed spectral albedo (350 – 2500 nm) measurements on tundra and sea ice surfaces. The ASD detector is carried in a backpack and controlled via a computer mounted on the front of the operator (see Figure 1). The ASD measures the spectral irradiance from a fiber optic cable that is routed from the backpack to a custom, gooseneck cosine collector mounted on the end of a 1-m long boom (Grenfell and Perovich, 2008). The boom was held at hip height (approximately 1 m) and had an integrated bubble level for levelling. To make an albedo measurement, first the operator collect an incident (down-welling) irradiance, followed by a reflected (up-welling) measurement. The time between incident and reflected measurements was typically between 11 and 26 seconds (interquartile range). For each measurement, 10 spectra are averaged together. Albedo is calculated as the ratio of the reflected to incident measurement, which obviates the need for absolute radiometric calibration. Albedo measurements were taken parallel to the 200-m albedo lines at 5-m increments (41 measurements) ~1 m south of the line. While the ASD operator was making measurements, an assistant kept notes on the scan number associated with each measurement, the surface type (see below), and collected photos of each measurement (see companion oblique photos data archive). Measurements were made within 3 hours of solar noon.

ASD Spectroradiometer↗
DOE OSTI2020

Soft fragmentation on the celestial sphere

We develop two approaches to the problem of soft fragmentation of hadrons in a gauge theory for high energy processes. The first approach directly adapts the standard resummation of the parton distribution function’s anomalous dimension (that of twist-two local operators) in the forward scattering regime, using k T -factorization and BFKL theory, to the case of the fragmentation function by exploiting the mapping between the dynamics of eikonal lines on transverse-plane to the celestial-sphere. Critically, to correctly resum the anomalous dimension of the fragmentation function under this mapping, one must pay careful attention to the role of regularization, despite the manifest collinear or infra- red finiteness of the BFKL equation. The anomalous dependence on energy in the celestial case, arising due to the mismatch of dimensionality between positions and angles, drives the differences between the space-like and time-like anomalous dimension of parton densities, even in a conformal theory. The second approach adapts an angular-ordered evolution equation, but working in 4 – 2ϵ dimensions at all angles. The two approaches are united by demanding that the anomalous dimension in 4 – 2ϵ dimensions for the parton distribution function determines the kernel for the angular-ordered evolution to all orders.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2016

Eulerian formulation of the interacting particle representation model of homogeneous turbulence

The Interacting Particle Representation Model (IPRM) of homogeneous turbulence incorporates information about the morphology of turbulent structures within the con nes of a one-point model. In the original formulation [Kassinos & Reynolds, Center for Turbulence Research: Annual Research Briefs, 31{51, (1996)], the IPRM was developed in a Lagrangian setting by evolving second moments of velocity conditional on a given gradient vector. In the present work, the IPRM is re-formulated in an Eulerian framework and evolution equations are developed for the marginal PDFs. Eulerian methods avoid the issues associated with statistical estimators used by Lagrangian approaches, such as slow convergence. A specific emphasis of this work is to use the IPRM to examine the long time evolution of homogeneous turbulence. We first describe the derivation of the marginal PDF in spherical coordinates, which reduces the number of independent variables and the cost associated with Eulerian simulations of PDF models. Next, a numerical method based on radial basis functions over a spherical domain is adapted to the IPRM. Finally, results obtained with the new Eulerian solution method are thoroughly analyzed. The sensitivity of the Eulerian simulations to parameters of the numerical scheme, such as the size of the time step and the shape parameter of the radial basis functions, is examined. A comparison between Eulerian and Lagrangian simulations is performed to discern the capabilities of each of the methods. Finally, a linear stability analysis based on the eigenvalues of the discrete differential operators is carried out for both the new Eulerian solution method and the original Lagrangian approach.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

Parallel Implicit Hydrodynamics with Material Strength for High Explosive Burn Calculations

High explosives are almost always evolving in some form of metal containment vessel. This fact requires that a materials model for the evolution of the metal containment vessel be part of any simulation of the HE. Since the actual form of the strength model to be used is an open question, we incorporate the material strength properties in a relatively agnostic fashion, which will accommodate many strength models. The only restriction is that the stress components be a function of density, specific internal energy and velocities. In addition, the timescales and rates of the chemical reactions in the HE vary between thousands of seconds and nanoseconds. A hydrodynamics capable of operating in anticipation of the eventual violent release of energy via these chemical reactions is a necessity. The central issue in creating a hydrodynamics capable of spanning these timescales is to forgo the use of a Courant time step control necessary in an explicit hydrodynamics. The Courant time-step arises because of characteristic velocities associated with the material, such as the sound speed or, in materials with strength, characteristic longitudinal and transverse speeds associated with compression and shear. The method used to circumvent the need for a Courant time-step limit is to develop an implicit calculation of the advanced particle pressure and shear modulus derived from the conservation laws of mass, momentum and specific internal energy. This process creates the implicit hydrodynamics needed for the HE calculations of interest to this study.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Reversible solid oxide cells: Early performance and microstructural evolution during electrolysis and switched mode operation

Incorporating more renewables into our energy mix will require a diversity of technological solutions. Reversible solid oxide cell systems (RSOCs) operating reversibly between fuel cell and electrolysis modes with high round-trip efficiencies offer one such solution. However, questions remain about the stability of the oxygen and fuel electrodes under reversible operation. Here, in this study, single cells with the configuration Ni-YSZ|YSZ|GDC10|NNO-NDC50|NNO are operated by reversible cycling under potentiostatic conditions. Test conditions are chosen to be representative of average stack conditions that are likely to be encountered during operation of a stack constructed of such cells. The results of reversible cycling over 500 h (21 cycles) are compared with operating the cells under potentiostatic conditions in the electrolysis-only mode (i.e., without cycling). Voltage-current density (DC), AC complex impedance spectroscopy, distribution of relaxation times (DRT) analysis, and scanning electron microscopy (SEM) indicate excellent reversibility and stability of the neodymium nickelate-based oxygen electrodes in both electrolysis and reversible cycling modes. Degradation behavior is attributed to the loss and coarsening of connected Ni particles in the fuel electrode, particularly in the active layer. However, reversible operation mitigates this degradation compared to operation in the electrolysis-only (SOEC-only) mode.

25 ENERGY STORAGE↗
DOE OSTI2020

Size of bulk fermions in the SYK model

The study of quantum gravity in the form of the holographic duality has uncovered and motivated the detailed investigation of various diagnostics of quantum chaos. One such measure is the operator size distribution, which characterizes the size of the support region of an operator and its evolution under Heisenberg evolution. In this work, we examine the role of the operator size distribution in holographic duality for the Sachdev-Ye-Kitaev (SYK) model. Using an explicit construction of AdS 2 bulk fermion operators in a putative dual of the low temperature SYK model, we study the operator size distribution of the boundary and bulk fermions. Our result provides a direct derivation of the relationship between (effective) operator size of both the boundary and bulk fermions and bulk SL(2;$\mathbb{R}$) generators.

2D gravity↗