Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “linear projection”

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 163 records · Page 9

Compact Coaxial HOM Damper for SRF Cavities

SRF technology enable particle accelerators operate with greater average beam currents and higher duty cycles. In these regimes parasitic excitation of the cavity HOM spectrum become the limiting factor due to extra RF losses and instabilities appearing in the beam. We discuss the limitations of HOMs for large accelerator projects as the Fermilab MI upgrade and the proposed future circular and linear colliders EIC and ILC. A new compact HOM damper concept is suggested providing required HOM suppression. The idea is based on an oversized coaxial line with radial sections, which is directly connected to the axial region of the SRF cavity. Such radially sectioned coaxial works as a filter for the operating mode, while being transparent for the HOMs. Implementations of this design for the accelerator and crab cavities developed for the MI and EIC machines, respectively, are discussed.

Lunin, A. [Fermilab]↗

Discriminative Dimensionality Reduction using Deep Neural Networks for Clustering of LIGO Data

In this paper, leveraging the capabilities of neural networks for modeling the non-linearities that exist in the data, we propose several models that can project data into a low dimensional, discriminative, and smooth manifold. The proposed models can transfer knowledge from the domain of known classes to a new domain where the classes are unknown. A clustering algorithm is further applied in the new domain to find potentially new classes from the pool of unlabeled data. The research problem and data for this paper originated from the Gravity Spy project which is a side project of Advanced Laser Interferometer Gravitational-wave Observatory (LIGO). The LIGO project aims at detecting cosmic gravitational waves using huge detectors. However non-cosmic, non-Gaussian disturbances known as "glitches", show up in gravitational-wave data of LIGO. This is undesirable as it creates problems for the gravitational wave detection process. Gravity Spy aids in glitch identification with the purpose of understanding their origin. Since new types of glitches appear over time, one of the objective of Gravity Spy is to create new glitch classes. Towards this task, we offer a methodology in this paper to accomplish this.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

MFIX-Exa: Performance prediction of multiphase energy conversion devices

MFIX-Exa targets the high-fidelity CFD-DEM model in which particles are unresolved by the fluid grid (typically using a fluid mesh approximately twice the particle diameter) but discrete particle dynamics including collisions are resolved with a simple linear-spring dashpot model. The fluid field is solved with a modern low Mach number formation of a cell-centered, nodal-pressure approximate projection method using a Godunov scheme . Physics capabilities include: complex geometries via embedded boundaries (EBs), open and closed system ideal gas equations of state, species transport and heterogeneous chemi- cal reactions. MFIX-Exa is built on the AMReX software framework, the ECP Block-Structured Adaptive Mesh Refinement (AMR) Co-Design Center, which allows the code to be portable and performant.Further integration into the ECP ecosystem includes HYPRE linear solvers and Ascent for in situ visualization. Current project efforts focus on scaling up realistic simulations to the KPP of the challenge problem: NETs 50kW chemical looping reactor, shakedown and scaling on JSLE and OLCF TDS systems and reaching for any remaining performance improvements.

Musser, Jordan↗

Upgrading Fermilab s Accelerator Control System with ACORN

The Fermilab Accelerator Complex is the largest national user facility in the Office of High Energy Physics (DOE/HEP) program and the only national user facility operating at Fermilab. Fermilab serves as the host to the Long Baseline Neutrino Facility/Deep Underground Neutrino Experiment (LBNF/DUNE), the laboratory’s flagship project for neutrino science that is under construction. LBNF/DUNE will be powered by megawatt beams from an upgraded accelerator, the Proton Improvement Plan II (PIP-II) that will replace the laboratory’s aging linear accelerator with a new one based on superconducting radio-frequency cavities. The Accelerator Controls Operations Research Network (ACORN) Project will support LBNF/DUNE and PIP-II by modernizing the accelerator control system. The project is at the conceptual design phase and looking to achieve Critical Decision 1 (CD-1) later this year. The scope and structure of the project will be presented, along with an overview of how that has changed in the past year. Current design and technology choices will be shared. Specific challenges facing the project will be addressed, along with current thinking on solutions.

Roehrig, Christian [Fermilab]↗

Evaluation of global teleconnections in CMIP6 climate projections using complex networks

In climatological research, the evaluation of climate models is one of the central research subjects. As an expression of large-scale dynamical processes, global teleconnections play a major role in interannual to decadal climate variability. Their realistic representation is an indispensable requirement for the simulation of climate change, both natural and anthropogenic. Therefore, the evaluation of global teleconnections is of utmost importance when assessing the physical plausibility of climate projections. We present an application of the graph-theoretical analysis tool δ-MAPS, which constructs complex networks on the basis of spatio-temporal gridded data sets, here sea surface temperature and geopotential height at 500 hPa. Complex networks complement more traditional methods in the analysis of climate variability, like the classification of circulation regimes or empirical orthogonal functions, assuming a new non-linear perspective. While doing so, a number of technical tools and metrics, borrowed from different fields of data science, are implemented into the δ-MAPS framework in order to overcome specific challenges posed by our target problem. Those are trend empirical orthogonal functions (EOFs), distance correlation and distance multicorrelation, and the structural similarity index. δ-MAPS is a two-stage algorithm. In the first place, it assembles grid cells with highly coherent temporal evolution into so-called domains. In a second step, the teleconnections between the domains are inferred by means of the non-linear distance correlation. We construct 2 unipartite and 1 bipartite network for 22 historical CMIP6 climate projections and 2 century-long coupled reanalyses (CERA-20C and 20CRv3). Potential non-stationarity is taken into account by the use of moving time windows. The networks derived from projection data are compared to those from reanalyses. Our results indicate that no single climate projection outperforms all others in every aspect of the evaluation. But there are indeed models which tend to perform better/worse in many aspects. Differences in model performance are generally low within the geopotential height unipartite networks but higher in sea surface temperature and most pronounced in the bipartite network representing the interaction between ocean and atmosphere.

58 GEOSCIENCES↗

Initial conditions for simulations of beam physics in linear induction accelerators

Flash radiography of hydrodynamic experiments driven by high explosives is a well-known diagnostic technique in use at many laboratories. At Los Alamos, the DualAxis Radiographic Hydrodynamic Test (DARHT) facility two linear induction electron accelerators (LIAs) make the bremsstrahlung radiographic source spots for point projection radiographs from orthogonal views. A new LIA, called Scorpius, is presently under development to advance this technology. To better understand electron-beam physics in these LIAs, numerical simulations are frequently performed with the objective of improving the radiography. At Los Alamos we frequently use the TRAK ray-trace and LSP particle-in-cell (PIC) codes to simulate the injector, and the XTR and LAMDA envelope/centroid codes along with LSP to simulate transport of the accelerated beam through the LIAs. The LIA simulations need the injected beam parameters as initial conditions for calculating beam transport and stability. The determination of these initial conditions is the topic of this note.

43 PARTICLE ACCELERATORS↗

Correct Initial Conditions for Simulations of Beam Physics in Linear Induction Accelerators

Flash radiography of hydrodynamic experiments driven by high explosives is a well-known diagnostic technique in use at many laboratories. At Los Alamos, the Dual Axis Radiographic Hydrodynamic Test (DARHT) facility two linear induction electron accelerators (LIAs) make the bremsstrahlung radiographic source spots for point projection radiographs from orthogonal views. A new LIA, called Scorpius, is presently under development to advance this technology. To better understand electron-beam physics in these LIAs, numerical simulations are frequently performed with the objective of improving the radiography. At Los Alamos we frequently use the TRAK ray-trace and LSP particle-in-cell (PIC) codes to simulate the injector, and the XTR and LAMDA envelope/centroid codes along with LSP to simulate transport of the accelerated beam through the LIAs. The LIA simulations need the injected beam parameters as initial conditions for calculating beam transport and stability. The determination of these initial conditions is the topic of this note.

43 PARTICLE ACCELERATORS↗

Physics Basis, Optimization, and Control for Integrated 3D Edge Long-pulse Tokamak Scenarios

The aim of this project with international collaborations (PI: J.-K. Park of Princeton Plasma Physics Laboratory) is to develop a unified physics basis and predictive capability for the control of edge-localized modes (ELMs) with optimized non-axisymmetric (3D) fields, by leveraging the unique research capabilities of international KSTAR tokamak in Korea. At the University of California, Irvine, we use a fusion simulation code GTC to study the effects of the 3D magnetic fields (including both magnetic islands and stochastic magnetic field lines) on microturbulence. We focus on the effects of 3D equilibrium on turbulent and neoclassical transport that leads to ELM suppression. During the project period, the UCI team of this collaboration has made significant progress in understanding turbulent transport in 3D geometry including both tokamaks with RMP and stellarators. Four UCI postdoctoral researchers, Drs. Javier Nicolau, Gyungjin Choi, Pengfei Liu, and Xishuo Wei were partially supported by this project and collaborating with the entire project team. We have published 4 journal papers and gave several invited talks. This project advance physics understand of 3D effects on microturbulence, specially effects of magnetic islands. It will contribute to the DOE SciDAC ISEP project in the area of GTC development for general 3D geometry capability, rigorous verification and validation, and readiness for the emerging exascale computing. Through collaborations with computational scientists, the project helps to advance high performance computing in the area of GPU acceleration, linear solver, and data management.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

A reduced-order model for nonlinear radiative transfer problems based on moment equations and POD-Petrov-Galerkin projection of the normalized Boltzmann transport equation

A data-driven projection-based reduced-order model (ROM) for nonlinear thermal radiative transfer (TRT) problems is presented. The TRT ROM is formulated by (i) a hierarchy of low-order quasidiffusion (aka variable Eddington factor) equations for moments of the radiation intensity and (ii) the normalized Boltzmann transport equation (BTE). The multilevel system of moment equations is derived by projection of the BTE onto a sequence of subspaces which represent elements of the phase space of the problem. Exact closure for the moment equations is provided by the Eddington tensor. A Petrov-Galerkin (PG) projection of the normalized BTE is formulated using a proper orthogonal decomposition (POD) basis representing the normalized radiation intensity over the whole phase space and time. The Eddington tensor linearly depends on the solution of the normalized BTE. By linear superposition of the POD basis functions, a low-rank expansion of the Eddington tensor is constructed with coefficients defined by the PG projected normalized BTE. The material energy balance (MEB) equation is coupled with the effective gray low-order equations which exist on the same dimensional scale as the MEB equation. The resulting TRT ROM is structure and asymptotic preserving. A detailed analysis of the ROM is performed on the classical Fleck-Cummings (F-C) TRT multigroup test problem in 2D geometry. Numerical results are presented to demonstrate the ROM's effectiveness in the simulation of radiation wave phenomena. Importantly, the ROM is shown to produce solutions with sufficiently high accuracy while using low-rank approximation of the normalized BTE solution. Essential physical characteristics of supersonic radiation wave are preserved in the ROM solutions.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Preconditioned least‐squares Petrov–Galerkin reduced order models

Abstract In this article, we introduce a methodology for improving the accuracy and efficiency of reduced order models (ROMs) constructed using the least‐squares Petrov–Galerkin (LSPG) projection method through the introduction of preconditioning. Unlike prior related work, which focuses on preconditioning the linear systems arising within the ROM numerical solution procedure to improve linear solver performance, our approach leverages a preconditioning matrix directly within the minimization problem underlying the LSPG formulation. Applying preconditioning in this way has the potential to improve ROM accuracy for several reasons. First, preconditioning the LSPG formulation changes the norm defining the residual minimization, which can improve the residual‐based stability constant bounding the ROM solution's error. The incorporation of a preconditioner into the LSPG formulation can have the additional effect of scaling the components of the residual being minimized to make them roughly of the same magnitude, which can be beneficial when applying the LSPG method to problems with disparate scales (e.g., dimensional equations, multi‐physics problems). Importantly, we demonstrate that an “ideal preconditioned” LSPG ROM (a ROM in which the preconditioner is the inverse of the Jacobian of its corresponding full order model) emulates projection of the full order model solution increment onto the reduced basis. This quantity defines a lower bound on the error of a ROM solution for a given reduced basis. By designing preconditioners that approximate the Jacobian inverse—as is common in designing preconditioners for solving linear systems—it is possible to obtain a ROM whose error approaches this lower bound. The proposed approach is evaluated on several mechanical and thermo‐mechanical problems implemented within the Albany HPC code and run in the predictive regime, with prediction across material parameter space. We demonstrate numerically that the introduction of simple Jacobi, Gauss‐Seidel, and ILU preconditioners into the proper orthogonal decomposition/LSPG formulation reduces significantly the ROM solution error, the reduced Jacobian condition number, the number of nonlinear iterations required to reach convergence, and the wall time (thereby improving efficiency). Moreover, our numerical results reveal that the introduction of preconditioning can deliver a robust and accurate solution for test cases in which the unpreconditioned LSPG method fails to converge.

Lindsay, Payton↗

Observable-projected ensembles

Measurements in many-body quantum systems can generate non-trivial phenomena, such as preparation of long-range entangled states, dynamical phase transitions, or measurement-altered criticality. Here, we introduce a new measurement scheme that produces an ensemble of mixed states in a subsystem, obtained by measuring a local Hermitian observable on part of its complement. We refer to this as the observable-projected ensemble . Unlike standard projected ensembles-where pure states are generated by projective measurements on the complement-our approach involves projective partial measurements of specific observables. This setup has two main advantages: theoretically, it is amenable to analytical computations, especially within conformal field theories. Experimentally, it requires only a linear number of measurements, rather than an exponential one, to probe the properties of the ensemble. As a first step in exploring the observable-projected ensemble, we investigate its entanglement properties in conformal field theory and perform a detailed analysis of the free compact boson.

Milekhin, Alexey [California Institute of Technolo↗

Bench-Scale Testing of Monolithic PPI Structured Contactors for Direct Air Capture of CO 2

The overall project objective was to develop, optimize and bench-scale test the integrated embodiment of a leading direct air capture (DAC) sorbent composition, linear-poly(propylenimine) (l-PPI), in a structured material system, specifically a monolithic contactor for achieving low pressure drop, to justify its further scale-up in a subsequent program. The project team was led by CORMETECH, a leading monolithic gas/solid contactor manufacturer for environmental applications, and included Global Thermostat, a well-known DAC start-up, and Georgia Tech, a prominent US academic institution in carbon capture and DAC. The monolithic contactor for l-PPI in this project employed an advanced, cost-effective fabrication technique from CORMETECH, namely a porous monolith substrate directly impregnated with the amine, an approach contrasting to the traditional wash-coating methods used for activating monoliths. The l-PPI monolithic contractor is targeted for use in DAC systems such as the Global Thermostat DAC process, which involves the cyclic operation of ambient air flow over a monolithic amine contactor followed by steam-mediated thermal desorption and CO 2 collection, maximizing volumetric productivity while reducing the auxiliary power required to capture CO 2 from air. The Global Thermostat DAC process currently utilizes poly(ethyleneimine) (PEI) materials as baseline sorbents, which are subject to oxidative degradation and slow capacity fade at the elevated temperature required for efficient CO 2 removal. Georgia Tech, with Global Thermostat, had discovered in lab-scale efforts that l-PPI had superior resistance to oxidative degradation compared to PEI, with similar CO 2 adsorption capacity. While the novel l-PPI sorbent offers significant advantages over PEI for DAC, including potential design simplification and increased process efficiency, it is not commercially available, and prior to this project, had not yet been evaluated on the bench-scale in a structured contactor.

36 MATERIALS SCIENCE↗

The future of the El Niño–Southern Oscillation: using large ensembles to illuminate time-varying responses and inter-model differences

Future changes in the El Niño–Southern Oscillation (ENSO) are uncertain, both because future projections differ between climate models and because the large internal variability of ENSO clouds the diagnosis of forced changes in observations and individual climate model simulations. By leveraging 14 single model initial-condition large ensembles (SMILEs), we robustly isolate the time-evolving response of ENSO sea surface temperature (SST) variability to anthropogenic forcing from internal variability in each SMILE. We find nonlinear changes in time in many models and considerable inter-model differences in projected changes in ENSO and the mean-state tropical Pacific zonal SST gradient. We demonstrate a linear relationship between the change in ENSO SST variability and the tropical Pacific zonal SST gradient, although forced changes in the tropical Pacific SST gradient often occur later in the 21st century than changes in ENSO SST variability, which can lead to departures from the linear relationship. Single-forcing SMILEs show a potential contribution of anthropogenic forcing (aerosols and greenhouse gases) to historical changes in ENSO SST variability, while the observed historical strengthening of the tropical Pacific SST gradient sits on the edge of the model spread for those models for which single-forcing SMILEs are available. Our results highlight the value of SMILEs for investigating time-dependent forced responses and inter-model differences in ENSO projections. The nonlinear changes in ENSO SST variability found in many models demonstrate the importance of characterizing this time-dependent behavior, as it implies that ENSO impacts may vary dramatically throughout the 21st century.

54 ENVIRONMENTAL SCIENCES↗

In-Process Melt Separation of Depolymerized PET in PET/PE Blends for Upcycling via Twin-Screw Extrusion

Our previous work focused on depolymerizing polyethylene terephthalate (PET) in twin-screw extrusion, as part of a broader project to continuously separate PET from polyolefins in the melt. This study focused on linear low-density polyethylene (LLDPE) and PET films and the use of ethylene glycol (EG), diethylene glycol (DEG), triethylene glycol (TEG), and bis(2-hydroxyethyl) terephthalate (BHET) to depolymerize the PET in the extruder to levels above 90% Mw. In this work, the focus will shift to achieving separation of the two polymers in the twin-screw extruder, which is made possible due to a 90% reduction in the Mw of the PET, which caused a decrease of its viscosity by several orders of magnitude. Owing to the viscosity difference and pressure buildup in the die, the low-viscosity PET preferentially exited a degassing vent instead of going through the die. This is because the flow of the PET would travel through a non-pressure vent rather than through a high-pressure die. However, owing to the higher viscosity of the LLDPE, the pressure was too high to pass through such a small diameter vent hole. Supercritical CO 2 (SCCO 2 ) was used to assist in this extraction, but SCCO 2 negatively impacted the overall degree of separation. Through analysis of the separated materials, it was concluded that a high separation of the two materials was achieved. TGA and FTIR confirmed that the material separated from the vent was 100% PET. The material removed from the die was composed of 95% LLDPE and 5% PET.

36 MATERIALS SCIENCE↗

In‐Process Melt Separation of PE/PET Blends for Upcycling via Twin‐Screw Extrusion. Impact of Catalyst Reagents on PET Depolymerization

Part 1 of this work focused on depolymerizing polyethylene terephthalate (PET) in twin-screw extrusion, as part of a broader project to continuously separate PET from polyolefins in the melt. This study focused on linear low-density polyethylene (LLDPE) and PET films and the use of ethylene glycol (EG), diethylene glycol (DEG), triethylene glycol (TEG), and bis(2-hydroxyethyl) terephthalate (BHET) to depolymerize the PET in the extruder to levels above 90% Mw. In part 2, the focus will shift to achieving separation of the two polymers in the twin-screw extruder, which is made possible due to a 90% reduction in the Mw of the PET, which caused a decrease of its viscosity by several orders of magnitude. Owing to the viscosity difference and pressure buildup in the die, the low-viscosity PET preferentially exited a degassing vent instead of going through the die. This is because the flow of the PET would travel through a non-pressure vent rather than through a high-pressure die. However, owing to the higher viscosity of the LLDPE, the pressure was too high to pass through such a small diameter vent hole. Supercritical CO₂ (SCCO₂) was used to assist in this extraction, but SCCO₂ negatively impacted the overall degree of separation. Through analysis of the separated materials, it was concluded that a very high separation of the two materials was achieved. TGA and FTIR confirmed that the material separated from the vent was 100% PET. The material removed from the die was composed of 90% LLDPE and 10% PET.

36 MATERIALS SCIENCE↗

AMR-Wind [SWR-20-85]

AMR-Wind is a massively parallel, block-structured adaptive-mesh, incompressible flow solver for wind turbine and wind farm simulations. The solver is built on top of the AMReX library. AMReX is developed at LBNL , NREL , and ANL as part of the Block-Structured AMR Co-Design Center in DOE's Exascale Computing Project. AMReX library provides the mesh data structures, mesh adaptivity, as well as the linear solvers used for solving the governing equations. The primary applications for AMR-Wind are: performing large-eddy simulations (LES) of atmospheric boundary layer (ABL) flows, simulating wind farm turbine-wake interactions using actuator disk or actuator line models for turbines, and as a background solver when coupled with a near-body solver with overset methodology to perform blade-resolved simulations of multiple wind turbines within a wind farm.

Ananthan, Shreyas↗

Materials Data on BaTiO3 by Materials Project

BaTiO3 crystallizes in the tetragonal P4/mmm space group. The structure is three-dimensional. Ba2+ is bonded in a linear geometry to two equivalent O2- atoms. Both Ba–O bond lengths are 2.52 Å. Ti4+ is bonded in a square co-planar geometry to four equivalent O2- atoms. All Ti–O bond lengths are 2.06 Å. There are two inequivalent O2- sites. In the first O2- site, O2- is bonded in a distorted linear geometry to two equivalent Ti4+ and two equivalent O2- atoms. Both O–O bond lengths are 2.06 Å. In the second O2- site, O2- is bonded to two equivalent Ba2+ and four equivalent O2- atoms to form distorted corner-sharing OBa2O4 octahedra. The corner-sharing octahedral tilt angles are 0°.

36 MATERIALS SCIENCE↗

Materials Data on AgAuS2 by Materials Project

AuAgS2 crystallizes in the orthorhombic Pccm space group. The structure is three-dimensional. Au3+ is bonded in a linear geometry to two equivalent S2- atoms. Both Au–S bond lengths are 2.31 Å. Ag1+ is bonded in a 4-coordinate geometry to four equivalent S2- atoms. All Ag–S bond lengths are 2.62 Å. S2- is bonded to one Au3+, two equivalent Ag1+, and one S2- atom to form a mixture of edge and corner-sharing SAg2AuS tetrahedra. The S–S bond length is 2.13 Å.

36 MATERIALS SCIENCE↗