Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “correspondence”

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 793 records · Page 44
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

A Big Problem for Small Earthquakes: Benchmarking Routine Magnitudes and Conversion Relationships with Coda Envelope-Derived Mw in Southern Kansas and Northern Oklahoma

Earthquake magnitudes are widely relied upon measures of earthquake size. Although moment magnitude (M w ) has become the established standard for moderate and large earthquakes, difficulty in reliably measuring seismic moments for small (generally M w <4) earthquakes has meant that magnitudes for these events remain plagued by a patchwork of inconsistent measurement scales. Because of this, magnitudes of small earthquakes and statistics derived from them can be biased. Furthermore, because small earthquakes are much more numerous than large ones, many applications, such as seismic hazard modeling, depend critically on analysis of events characterized by magnitudes other than M w . Therefore, to assess this problem, we apply coda envelope analysis to reliably determine moment magnitudes for a case study of small earthquakes from northern Oklahoma and southern Kansas. Not surprisingly, we find significant differences among M L , m bLg , and M w for M ~2–4 earthquakes examined here. More troublingly, we find that relations designed to convert other magnitudes to M w , which are relied upon for important applications such as seismic hazard analysis, often increase rather than decrease this bias for our dataset. In our case study, we find that converted magnitudes can result in a systematic bias sometimes exceeding 0.5 magnitude units, a difference that typically corresponds to a factor of ~3 in seismicity rate. Moreover, we find a correspondingly large bias in Gutenberg–Richter b-values, controlled primarily by inaccurate magnitude scaling in the conversion relationships. Although this study focuses on a relatively small geographic area, we can expect that similar issues exist with varying severity in other regions. Therefore, magnitudes of small earthquakes and their associated statistics, including seismicity rates and b-values, should be treated with caution.

58 GEOSCIENCES↗
DOE OSTI2024

Investigation of hot-spots due to trapped flux in niobium superconducting radiofrequency cavities

One of the significant sources of residual losses in superconducting radio-frequency cavities is magnetic flux trapped during the cool-down due to the incomplete Messier effect. If the trapped vortices are non-uniformly distributed on the cavity surface, the temperature mapping revealed the “hotspots” at the location of high density of pinned vortices. Here, we performed a rf test on 1.3 GHz single cell cavity with the combination of the temperature mapping system. The temperature mapping reveled the development of the hot spots with the increase in rf field inside the cavity. When magnetic field is trapped locally on the surface of cavity, the hot-spots strength increase rapidly, showing the direct correlation of vortex induced hot spot and corresponding rf loss.

Accelerator Physics↗
DOE OSTI2025

Towards Improving luminosity using optics tuning and data-driven methods

The results of Run 24 experiments at Relativistic Heavy Ion Collider (RHIC) for improving luminosity using optics tuning are presented in this study. In the first experiment, MADx matching was used to output magnet strengths corresponding to specific s star movements around Interaction Region 8 (IR8). The corresponding Zero Degree Calorimeter (ZDC) signal was measured in place of luminosity, and Bayesian Optimization aids search of optimal movements. It was found that values retrieved from matching were inaccurate, resulting in negative feedback loops. The second experiment focused on calculating accurate s star movements. The matching method was replaced with a linear sensitivity matrix, directly relating optics to power supply, and its null space was used to fit constraints such as hysteresis effects. At the experiment, beam losses were observed at collimators around boundary of IR8, which were fixed for the third experiment. Dynamic mode decomposition was also introduced to improve quality of turn-by-turn (TBT) data as well as accuracy and consistency of optics measurements at IR8. These improvements will be tested in the experiment of next RHIC run for luminosity optimization.

Accelerator Physics↗
DOE OSTI2021

Quantum field theory and the Bieberbach conjecture

An intriguing correspondence between ingredients in geometric function theory related to the famous Bieberbach conjecture (de Branges' theorem) and the non-perturbative crossing symmetric representation of 2-2 scattering amplitudes of identical scalars is pointed out. Using the dispersion relation and unitarity, we are able to derive several inequalities, analogous to those which arise in the discussions of the Bieberbach conjecture. We derive new and strong bounds on the ratio of certain Wilson coefficients and demonstrate that these are obeyed in one-loop Φ 4 theory, tree level string theory as well as in the S-matrix bootstrap. Further, we find two sided bounds on the magnitude of the scattering amplitude, which are shown to be respected in all the contexts mentioned above. Translated to the usual Mandelstam variables, for large |s|, fixed t, the upper bound reads |M(s,t)| ≲ |s 2 |. We discuss how Szegö's theorem corresponds to a check of univalence in an EFT expansion, while how the Grunsky inequalities translate into nontrivial, nonlinear inequalities on the Wilson coefficients.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

Twistor space origins of the Newman-Penrose map

Recently, we introduced the “Newman-Penrose map”, a novel correspondence between a certain class of solutions of Einstein’s equations and self-dual solutions of the vacuum Maxwell equations, which we showed was closely related to the classical double copy. Here, we give an alternative definition of this correspondence in terms of quantities that are defined naturally on twistor space, and a shear-free null geodesic congruence on Minkowski space whose twistorial character is articulated by the Kerr theorem. The advantage of this reformulation is that it is purely geometrical in nature, being manifestly invariant under both spacetime diffeomorphisms and projective transformations on twistor space. While the original formulation of the map may be more convenient for most explicit calculations, the twistorial formulation we present here may be of greater theoretical utility.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

Hamiltonian structure of 2D fluid dynamics with broken parity

Isotropic fluids in two spatial dimensions can break parity symmetry and sustain transverse stresses which do not lead to dissipation. Corresponding transport coefficients include odd viscosity, odd torque, and odd pressure. We consider an isotropic Galilean invariant fluid dynamics in the adiabatic regime with momentum and particle density conservation. We find conditions on transport coefficients that correspond to dissipationless and separately to Hamiltonian fluid dynamics. The restriction on the transport coefficients will help identify what kind of hydrodynamics can be obtained by coarse-graining a microscopic Hamiltonian system. Interestingly, not all parity-breaking transport coefficients lead to energy conservation and, generally, the fluid dynamics is energy conserving but not Hamiltonian. We show how this dynamics can be realized by imposing a nonholonomic constraint on the Hamiltonian system.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

Lieb-Schultz-Mattis anomalies as obstructions to gauging (non-on-site) symmetries

We study ’t Hooft anomalies of global symmetries in 1+1d lattice Hamiltonian systems. We consider anomalies in internal and lattice translation symmetries. We derive a microscopic formula for the “anomaly cocycle” using topological defects implementing twisted boundary conditions. The anomaly takes value in the cohomology group H^3(G,U(1)) × H^2(G,U(1)) H 3 ( G , U ( 1 ) ) × H 2 ( G , U ( 1 ) ) . The first factor captures the anomaly in the internal symmetry group G, and the second factor corresponds to a generalized Lieb-Schultz-Mattis anomaly involving G and lattice translation. We present a systematic procedure to gauge internal symmetries (that may not act on-site) on the lattice. We show that the anomaly cocycle is the obstruction to gauging the internal symmetry while preserving the lattice translation symmetry. As an application, we construct anomaly-free chiral lattice gauge theories. We demonstrate a one-to-one correspondence between (locality-preserving) symmetry operators and topological defects, which is essential for the results we prove. We also discuss the generalization to fermionic theories. Finally, we construct non-invertible lattice translation symmetries by gauging internal symmetries with a Lieb-Schultz-Mattis anomaly.

Physics↗
DOE OSTI2024

Lattice regularization of reduced Kähler-Dirac fermions and connections to chiral fermions

We show how a path integral for reduced Kähler-Dirac fermions suffers from a phase ambiguity associated with the fermion measure that is an analog of the measure problem seen for chiral fermions. However, unlike the case of chiral symmetry, a doubler free lattice action exists which is invariant under the corresponding onsite symmetry. This allows for a clear diagnosis and solution to the problem using mirror fermions resulting in a unique gauge invariant measure. By introducing an appropriate set of Yukawa interactions which are consistent with ’t Hooft anomaly cancellation we conjecture the mirrors can be decoupled from low energy physics. Moreover, the minimal such Kähler-Dirac mirror model yields a light sector which corresponds, in the flat space continuum limit, to the Pati-Salam GUT model.

Physics↗
DOE OSTI2020

Pre-Ionization Controlled Laser Plasma Formation for Ignition Applications (Final Scientific/Technical Report)

This report consists of three grand funded papers. Abstracts for each follow. 1) In this work, we present experimental and modeling studies of air pre-ionization using ultraviolet (UV) laser pulses and its effect on laser breakdown of an overlapped near-infrared (NIR) pulse. Experimental studies are conducted with a 266 nm beam (fourth harmonic of Nd:YAG) for UV pre-ionization and an overlapped 1064 nm NIR beam (fundamental of Nd:YAG), both having pulse duration of 10 ns. Results show that the UV beam produces a pre-ionized volume which assists in breakdown of the NIR beam, leading to reduction in NIR breakdown threshold by factor of >2. Numerical modeling is performed to examine the ionization and breakdown of both beams. The modeled breakdown threshold of the NIR, including assist by pre-ionization, is in reasonable agreement with the experimental results. 2) The present contribution compares the energy absorption, optical emission, temperature, and fluid dynamics of ultraviolet (UV) λ = 266 nm and near infrared (NIR) λ = 1064 nm nanosecond laser induced plasmas in ambient air. For UV pulses at the conditions studied, energy absorption by the plasmas increases relatively gradually with laser pulse energy starting at delivered energy of E~8 mJ. Corresponding measurements of plasma luminosity show that the absorption of UV radiation does not necessarily result in visible plasma emission. For the NIR induced plasmas, the energy absorption profile is far more abrupt and begins at ~55 mJ. In contrast with UV, the absorption of NIR radiation is always accompanied by intense optical emission. The temperatures of both types of plasma have been measured with Rayleigh scattering thermometry (at times after the Thomson signal sufficiently diminishes). The UV plasmas can attain a wider range of temperatures, including lower temperatures, depending on the pulse energy (e.g., T ~400–2000K for E ~7–35 mJ at Δt = 10 μs after the pulse) while the NIR plasmas show only hotter temperatures (e.g., T ~ 12 000 K for E = 75 mJ at Δt = 10 μs after the pulse) as is consistent with the literature. Differences in the fluid dynamics for UV versus NIR pulses are shown with Schlieren imaging. In this work, the contrast in the UV and NIR plasma threshold behavior is attributed to differing roles of avalanche ionization and multiphoton ionization as is also illustrated by a simple numerical model. 3) The present contribution examines the impact of plasma dynamics and plasma-driven fluid dynamics on the flame growth of laser ignited mixtures and shows that a new dual-pulse scheme can be used to control the kernel formation process in ways that extend the lean ignition limit. We perform a comparative study between (conventional) single-pulse laser ignition (λ = 1064 nm) and a novel dual-pulse method based on combining an ultraviolet (UV) pre-ionization pulse (λ = 266 nm) with an overlapped near-infrared (NIR) energy addition pulse (λ = 1064 nm). We employ OH* chemiluminescence to visualize the evolution of the early flame kernel. For single-pulse laser ignition at lean conditions, the flame kernel separates through third lobe detachment, corresponding to high strain rates that extinguish the flame. In this work, we investigate the capabilities of the dual-pulse to control the plasma-driven fluid dynamics by adjusting the axial offset of the two focal points. In particular, we find there exists a beam waist offset whereby the resulting vorticity suppresses formation of the third lobe, consequently reducing flame stretch. With this approach, we demonstrate that the dual-pulse method enables reduced flame speeds (at early times), an extended lean limit, increased combustion efficiency, and decreased laser energy requirements.

Mechanical engineering↗
DOE OSTI2020

Low-Temperature Plasma Treatment for Enhanced Recovery of Highly Valued Critical REEs from Coal (Final Report)

The principle objective of this project was to develop a novel process flowsheet using low-temperature plasma treatment (LTP) integrated with conventional hydrometallurgical processes to recover rare earth elements (REEs), especially highly valued REEs (i.e. scandium and critical REEs), from coal and coal byproducts. LTP treatment was used to pretreat the solid feed material prior to leaching in an effort to improve the leaching characteristics of the REEs and potentially reduce acid consumption. Initially, a detailed experimental program was performed on a laboratory LTP furnace to identify a set of process parameter values to be used for assessing the enhanced leaching performance achieved on four specific gravity fractions of West Kentucky No. 13 and Fire Clay coarse coal refuse materials. High-temperature oxidation (HTO)) was also evaluated as a means for improving REE leachability on the same fractions of both coal sources. The lixiviants used to assess leachability improvements included: i) untreated water, ii) 0.1 M ammonium sulfate solution and iii) 1.2 M sulfuric acid solution at a temperature of 75°C. The laboratory data was used to develop multiple process flow sheets using HTO only to produce a rare earth oxide mix concentrate or utilizing both LTP and HTO treatment processes to generate a high purity heavy REO mix concentrate and a light REO mix concentrate. A techno-economic analysis was performed on each circuit. Low-temperature plasma oxidation was found to enhance REE leaching characteristics for all specific gravity (SG) fractions of the West Kentucky No. 13 coarse refuse material. The leaching recovery improvement was preferentially achieved for the heavy REEs with higher values obtained when treating the lower density fraction material. For example, the leaching recovery for the heavy REEs in the untreated 1.6 SG Float fraction material was 8% in a 0.1 M ammonium sulfate solution. After LTP treatment for 32 hours, heavy REE recovery increased to 33%. The corresponding recovery value for the light REEs was 10%. These findings indicate a strong affinity of the heavy REEs with the organic material, which is the component being reduced in content during the LTP treatment. Scanning electron micrographs showed the development of a honeycomb structure with increased LTP treatment with higher surface area as confirmed from BET measurements. Higher surface area promotes improved leaching efficiency due to enhance chemical interactions. When using 1.2 M sulfuric acid, the LTP treatment was found to provide a heavy REE recovery increase from 23% to 53% for the 1.6 SG Float fraction which corresponded to an ash content increase from around 20% to 35%. However, LTP treatment provided no improvement in leaching characteristic for the 2.2 SG sink material in the acid solution, which was likely due to the low organic content. The LTP treatment of the Fire Clay SG fractions provided no improvement in REE leaching characteristic, which may be due to the association with REE minerals such as monazite and REE-associated mineral matter such as clays. High temperature oxidation using temperatures in the range of 600-750OC was found to significantly enhance the leaching characteristics of both the light and heavy REEs for all density fractions in both coarse coal refuse materials. Generally, the increase in light REE recovery for the West Kentucky material was from 15% for the untreated material to 90% for the HTO treated material at 600°C. Recovery values dropped significantly when using temperatures below 600°C and above 750°C. The impact of the excessive temperatures was most likely a result of sintering. For the heavy REEs, HTO significantly improved recovery in the lighter SG fractions, which may be due to an affinity with the organic content. Only limited heavy REE recovery improvement was obtained for the denser SG fractions. Another benefit of HTO is the conversion of calcium-based minerals to a very soluble calcium oxide, which is removable using a mild acid solution, and the conversion of Fe minerals to relatively insoluble minerals such as hematite and magnetite. The pregnant leaching solution produced from the leaching of all SG fractions obtained from both coal sources was treated using two stages of solvent extraction procedurally modeling a rougher-cleaner circuit followed by oxalic acid precipitation. Rare earth oxide products having a purity greater than 97% was produced from nearly all of the SG fractions. A novel process flowsheet was developed based on the selective leaching of heavy REEs that was made possible through LTP treatment of West Kentucky No. 13 coarse refuse followed by salt leaching. In the flowsheet, coal refuse material is ground to a top size of 1mm and then treated by LTP. The treated material is mixed into a 0.1 M ammonium sulfate solution and then filtered to produce a heavy REE enriched pregnant leach solution. The heavy REEs are concentrated by two stages of selective precipitation to produce a heavy REE mix oxalate and then roasted to obtain a heavy REO mix concentrate. The solid residues remaining after salt leaching is dried and then subjected to HTO treatment followed by concentration using selective precipitation and roasting to produce a light REO mix concentrate. An alternative circuit was developed using HTO treatment only to treat 100% of the coarse refuse material followed by selective precipitation and roasting to generate a high purity REO mix product. A techno-economic analysis performed for the treatment of West Kentucky No. 13 and Fire Clay coarse refuse based on a hypothetical 500-tph plant revealed that the net present value for all scenarios was significantly negative. Primary reasons include: i) a substantial decrease in the basket value of the final product relative to the feed and ii) high chemical costs relative to the amount of product generated. The basket value reduction was primarily due to little to no scandium recovery due to poor leaching characteristics when using low acid concentrations as necessitated by the low feed grades. Operating costs were primarily influenced by the chemical costs relative to the kilograms of concentrate produced and feed grade. For the Fire Clay coarse refuse material, the total operating cost was estimated to be $354 per kg of REO produced which was higher than the basket value of the REO product ($276/kg). Approximately 37% of the total cost was the sulfuric acid used during the leaching process. A comparative analysis performed to assess the economic viability of using low-temperature plasma furnaces found that market value enhancement for the generation of both heavy REO and light REO concentrates will be needed in addition to significant technology advancements to reduce the cost of plasma furnaces.

01 COAL, LIGNITE, AND PEAT↗
DOE OSTI2020

Investigation of Thermolytic Hydrogen Generation Rate in Tank 44 Dissolved Saltcake Samples

Saltcake core samples collected from Tank 44 in 2006 were dissolved to provide material for HGR measurements applicable to F-Area dissolved saltcake material. Additionally, characterization was performed on the Tank 44 saltcake material. The following are key results from the Tank 44 saltcake characterization. The Tank 44 Upper Saltcake Composite, corresponding to the 171 to 285 inch tank level, contained by mass approximately 69% sodium nitrate, 11% sodium carbonate, 8% sodium nitrite, smaller amounts of other salts and components, and 9% unquantified (which includes water, water of hydration, oxygen/hydrogen content of oxides and hydroxides, and uncertainty). The Tank 44 Lower Saltcake Composite, corresponding to the 76 to 114 inch tank level, contained by mass approximately 49% sodium carbonate, 18% sodium nitrate, smaller amounts of other salts, at least 8% sludge, and 9% unquantified (see above). The dissolved saltcake contained free hydroxide less than quantifiable (<0.01 M) due to the limited quantity of material that could be removed from the Shielded Cells based on the sample radioactivity. Measurement by pH paper provided an approximate pH of 12. The following are key results from the Tank 44 HGR testing. During boiling at 106.7 °C, HGR for Tank 44 dissolved saltcake without added glycolate was 7.2×10 -8 ft 3 h -1 ga l-1 . During boiling at 106.9 °C, HGR for Tank 44 dissolved saltcake with 1000 mg/L of added glycolate was 8.2×10 -8 ft 3 h -1 gal -1 . ∙ For the test without added glycolate, the first several HGR measurements at 70, 85, and 100 °C gave indication of the release of dissolved hydrogen and should not be used to represent the sustained thermolytic HGR for those temperatures. The measurements at boiling are the best representation of thermolysis in this testing. Carbon dioxide was observed at concentrations up to 6 vol% in the flow-system offgas for the test at boiling. ∙ Methane generation was observed at 100 °C and boiling. Methane concentration in the total gas generated during testing remained well below the lower flammability limit for methane in air. The addition of 1000 mg/L of glycolate did not have a significant impact on the hydrogen generation rates measured during this testing. The low hydroxide concentration in the Tank 44 dissolved saltcake likely influenced the relatively low thermolytic HGR and high carbon dioxide release observations in this testing. Based on the observation that methane was generated or released upon heating SRS radioactive Tank 44 waste samples to 100 °C and above, we recommend gaining a greater understanding of the cause and mechanism of its generation. First, the applicable literature should be reviewed to reveal the thermolytic methane generation mechanisms of possible methane generating species in the SRS CSTF. If warranted, a plan should be developed for simulant tests with methylated siloxanes and other applicable compounds in order to gain a better mechanistic understanding of methane generation in the SRS CSTF.

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

H2OI95: A Stand-Alone Fortran Code for Evaluating the IAPWS-95 Equation-of-State Model for Water (Rev. 1)

H2OI95 is a stand-alone Fortran code for evaluating the IAPWS-95 equation-of-state model for the thermodynamic properties of water. It also evaluates the corresponding thermochemical properties of water consistent with the CODATA recommendations. The IAPWS-95 model is based on a model equation for the dimensionless Helmholtz energy for which the primary variables are the inverse reduced temperature (τ = T cr /T ) and reduced density (δ = ρ/ρ cr ). Here T is the absolute temperature (K), ρ is density (kg/m³), and the subscript “ cr ” refers to the critical point of water (647.096 K and 22.064 MPa pressure in this model, for which ρ cr is 322 kg/m³). The code solves four basic types of problems, distinguished by the specified inputs: 1. Temperature (K) and density (ρ kg/m³) or reduced density (δ) 2. Temperature (K) and pressure (MPa). 3. Temperature (K) on the saturation (liquid-vapor equilibrium) curve 4. Pressure (MPa) on the saturation curve. Each type of problem is run using a corresponding input (text) file. All but the first type of problem require iteration. For example, to solve for desired temperature and pressure, the reduced density must be adjusted to give the desired pressure. Iteration is accomplished using the Newton-Raphson method, though the secant method is also used in solving the fourth type of problem. For the last three types of problems, H2OI95 has been used to conduct numerical studies of convergence and the problem of multiple numerical solutions, some of which are not physically valid. Obtaining valid results depends mainly on appropriate choice of starting values for the density. The default values used in H2OI95 appear to consistently lead to generally desired results. With modification (not addressed here), H2OI95 can be used to support SUPCRT92 and similar codes that compute chemical thermodynamic properties of species and reactions over a wide range of temperature and pressure (273.16-1273K and 0-1000 MPa). Version 1.1 of this software calculates an extended set of water properties, including some derivable solely from IAPWS-95, such as the isobaric thermal expansivity and the isothermal compressibility, and still other properties based on supplementary models sanctioned by IAPWS, including the dynamic viscosity, thermal conductivity, surface tension, and static dielectric constant, and various other properties (notably Debye-Hückel parameters) derivable from the dielectric constant and its partial derivatives with respect to temperature and pressure. A coding error in Version 1.0 was identified while testing the implementation of the viscosity model. This error is described in Appendix C. It mainly affects calculation of “higher order” properties such as compressibility and heat capacity near the critical point.

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

Cross section measurements of deuteron electro-disintegration at very high recoil momenta and large 4-momentum transfers (Q^2)

The 2H(e,e'p)n cross sections have been measured at negative 4-momentum transfers of Q^2 = 4.5 +/- 0.5 (GeV/c)^2 and Q^2 = 3.5 +/- 0.5 (GeV/c)^2 reaching neutron recoil (missing) momenta up to p_r ?1.0 GeV/c. The data have been obtained at fixed neutron recoil angles 5 <= theta_nq <= 95 degrees with respect to the 3-momentum transfer q. The new data agree well with the previous data which reached p_r ?550 MeV/c. At theta_nq = 35 and 45 degrees, final state interactions (FSI), meson exchange currents (MEC) and isobar configurations (IC) are suppressed and the plane wave impulse approximation (PWIA) provides the dominant cross section contribution. The new data are compared to recent theoretical calculations, and a significant disagreement for recoil momenta p_r > 700 MeV/c is observed. The experiment was carried out in experimental Hall C at the Thomas Jefferson National Accelerator Facility (TJNAF) and formed part of a group of four experiments that were used to commission the new Super High Momentum Spectrometer (SHMS). The experiment consisted of a 10.6 GeV electron beam incident on a liquid deuterium target which resulted in the break-up of the deuteron into a proton and neutron. The scattered electrons were detected by the SHMS in coincidence with the knocked-out protons detected in the previously existing High Momentum Spectrometer (HMS) and the recoiling neutrons were reconstructed from energy-momentum conservation laws. To ensure that the 2H(e,e'p)n reaction channel was selected, we required the missing energy of the system to be the binding energy of the deuteron (?2.22 MeV). The spectrometers? central angles and momenta were set to measure three central missing momentum settings of the neutron corresponding to p_r = 80, 580 and 750 MeV/c, which required the SHMS central angle and momentum to be fixed and the HMS to be rotated from smaller to larger angles corresponding to the lower and higher missing momentum settings, respectively. The experiment was carried out in a time period of six days with typical electron beam currents of 45-60 uA at about 50% beam efficiency.

Yero, Carlos↗
DOE OSTI2020

Development of Control System Functional Capabilities within the IES Plug-and-Play Simulation Environment

The concept of an Integrated Energy System (IES) is meant to combine different energy technologies in synergistic ways to achieve a more secure and economical energy supply. The RAVEN-based HYBRID framework is used to find the optimal installed capacity and the optimal economical dispatch of each component of the IES. A new RAVEN (Risk Analysis Virtual ENvironment) plugin for grid and capacity optimization (HERON) has been developed for optimizing the production variables of the IES given the demand profile. Currently, only the limits that affect the production variables and their corresponding time rates of change are considered (explicit constraints). However, other variables are additionally subject to constraints, but the associated limits are not accounted for (implicit constraints). In particular, for the power dispatch problem, the optimization algorithm takes into account the limits on the electrical power output and the corresponding hourly power variations but does not consider other constraints on process variables whose response effects the service life of the IES. This report describes a scheme that allows accounting for implicit constraints without increasing the size of the optimization problem. The Reference Governor (RG) algorithm is traditionally used for enforcing state and control constraints by modifying the set-point trajectories supplied to the feedback regulators. In our application, the RG is coupled within an iterative loop with the HERON-power dispatcher to generate optimal trajectories that ensure the operational constraints are met. A data-driven procedure to derive a representation of the dynamics of the controlled system was developed. First, the variables that represented the state of the system are selected (PCA-based approach), and then state-space representation matrices are derived from the collected measurements (DMDc algorithm). A preliminary version of the developed workflow based on Linear Time Invariant matrices was assessed by adopting a two-unit test case. More sophisticated versions of this workflow foreseeing the on-line derivation of system matrices will be deployed in FY 2021. Finally, a “plug-and-play” library of controllers and state observers was developed in Dymola. Some aspects of the current configuration of the IES unit components, e.g., the encapsulation of the control schemes into dedicated blocks, are consistent with the “plug-and-play” philosophy. Other features, e.g., the system buses collecting the input and the output variables, are not. For this reason, once listed and described the limits of the current configuration, necessary modifications to the plant model interface are presented. As a test-case, the interfaces of the SES model in the RAVEN-based HYBRID framework were reworked accordingly, and two different control schemes were applied to the same plant model.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2020

Oxy-Combustion System Process Optimization

The overall objective of this work is to develop a new chemical absorbent-based, high pressure, CO 2 purification system to remove the residual oxygen that currently contaminates the recovered CO 2 , and to optimize the Pressurized Oxy-Combustion (POxC) process to minimize the Cost of Electricity (COE) generated in this advanced combustion process. TDA developed and validated the performance of the oxygen removal system for CO 2 purification. In collaboration with the Advanced Power & Energy Program (APEP) of University of California, Irvine (UCI), we optimized the POxC process, including thermal management, heat integration, and power cycle optimization using process design and modeling supported with Aspen Plus® process simulations. The techno-economic analysis results indicate that the pressurized oxycombustion coal power plant with Ion-transport membrane (ITM) air separation unit (ASU) (Case 2 – 30.55%) does not show an advantage over a cryogenic ASU (Case 1 – 31.24%) while TDA’s sorption-based ASU (Case 4 - 32.61%) shows a significant advantage over the cryogenic ASU (Case 1 – 331.24%). The specific plant costs show a wide range with a low of $2544/kW for Case 11C (TDA ASU, co-sequestering the SO x , and ultra-supercritical steam cycle) to a high of $2975/kW for Case 12A (cryogenic ASU and sCO 2 cycle). In general, the ITM ASU based cases have lower specific plant costs than the corresponding cryogenic ASU based cases while the TDA ASU based cases show the lowest specific plant costs. The main reason for lowering these costs is the higher overall plant thermal efficiency which decreases the plant cost on a $ per kW basis. Next comparing the cases with different power cycle working fluid conditions in terms of temperature and pressure while all utilizing steam, similar trends as the plant costs may be observed. However, with the supercritical CO 2 (sCO 2 ) cycle, the increase in thermal efficiency of the sCO 2 cycle was not able to offset its increase in plant cost making the plant costs higher than those of the corresponding steam cycle cases. The Cost of Electricity (COE) again shows similar trends as the specific plant costs. The COE for Case 11C at $110.1/MWh is also the lowest, but among all cases that do not co-sequester the SO x , Case 12C (TDA ASU and sCO 2 cycle) has the lowest COE at $\$$117.5/MWh while the highest is for Case 8A (cryogenic ASU and supercritical steam cycle) at $130.4/MWh.

20 FOSSIL-FUELED POWER PLANTS↗
DOE OSTI2020

Search for resonance in inclusive and b-tagged dijet mass spectra in proton-proton collision at $\sqrt{s}$ = 13 TeV and constraints on dark matter and other models

Two analyses searching for new heavy resonances and interpreting the result with dark matter are presented in this dissertation. The first analysis searches for resonances decaying into dijet final states using the proton-proton collision data at center mass of $\sqrt{s}=13$ TeV with the CMS detector at Large Hadron Collider in2016, corresponding to a total integrated luminosity of 36 fb -1 . No significant evidence for the production of new particles is observed. A simplified model of interactions between quarks and dark matter is used to interpret the result. The second analysis searches for resonances decaying into dijet final states using the proton-proton collision data at center mass of $\sqrt{s}=13$ TeV with the CMS detector at Large Hadron Collider in 2016, 2017, and 2018, corresponding to a total integrated luminosity of 137.5 fb -1 . A b-tagger named DeepJet is used to examine the events with jets containing B hadrons. No evidence of new heavy resonances is observed, as well. The 95% confidence level cross-section upper limit is set on excited b-quark b* and Coloron models. Limits are also set on a simplified model of interactions between b-quarks and dark matter.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2021

Validation and Demonstration of Control System Functional Capabilities within the IES Plug-and-Play Simulation Environment

The concept of an integrated energy system (IES) is meant to combine different energy technologies in synergistic ways to achieve a more secure and economical energy supply. The RAVEN-based HYBRID framework is used to find the optimal installed capacity and the optimal economical dispatch of each component of the IES. The new RAVEN plugin for grid and capacity optimization (HERON) only addresses the limits that affect the production variables and the corresponding rates of variation (explicit constraints). However, other variables are subject to constraints, and the associated limits should be accounted for (implicit constraints). In particular, for the power dispatch problem, the optimization algorithm takes into account the limits on the electrical power output and the corresponding hourly power variations but does not consider other constraints on process variables whose response affects the service life of the IES. This report describes a scheme that allows accounting for implicit constraints without increasing the size of the optimization problem. To obtain a more accurate approximation of the nonlinear dynamic behavior, a parametric version of the dynamic mode decomposition with control (DMDc) algorithm was developed to derive the state-space representation matrices of the IES components at different scheduling parameter. Thanks to this approach, a more accurate approximation of the system response can be obtained, the limits imposed by thermal mechanical implicit constraints can be translated into power dispatch limits, and the feedbacks to HERON power dispatcher can be provided. To assess the developed methodology, a power dispatching test case composed of three power generating and storage units (Balance of Plant, Secondary Energy Source, Thermal Energy Storage) was developed. The power output of each one of the three units was optimized to meet the imposed time-dependent load demand trajectory and to maximize the IES profitability by meeting both the explicit and implicit constraints.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2021

Evaluation of Two-Dimensional to One-Dimensional Site Response for Idaho National Laboratory

We perform two-dimensional (2D) site response analyses accounting for spatial variability of soil properties and subsurface geometry of the Eastern Snake River Plane (ESRP), and quantify their effects on ground surface motion relative to one-dimensional (1D) site response analyses at the Idaho National Laboratory (INL). We first present the development of random field idealizations of the repeated basalt lava flows, heterogeneously inter-layered with sediments, from seismic velocity data collected over four decades in the ESRP. Using realizations of the stochastic fields mapped on 2D deterministic finite element models, we perform 2D viscoelastic and equivalent-linear wave propagation simulations, and quantify the mean and variance of site response aggravation factors, defined as the response spectral ratio of 2D to 1D analyses on the ground surface. Results are shown to be insensitive to the constitutive material behavior considered here, for strains induced by rock outcrop peak ground acceleration (PGA) as high as 0.7g: viscoelastic and equivalent-linear analyses predict peak mean 2D/1D aggravation factor 1.05 at period T=0.075 sec (i.e. the 2D response spectrum is 5% higher than the corresponding 1D at that period, on average), which corresponds to the wavelength of the horizontal correlation length of the random field (50m). For periods longer than the fundamental period of the site (here, T 1 =0.3125 sec), the propagating wavelengths are too long to be affected by the 1D site response and the aggravation factor becomes equal to 1. The standard deviation of the natural logarithms of the 2D/1D aggravation factors is ~0.15 for periods shorter than the fundamental period of the site, and decays thereafter at a steady rate.

58 GEOSCIENCES↗