Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Load flow”

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 991 records · Page 55

Local measurements of the spatial magnetic field distribution in a z-pinch plasma during and near stagnation using polarization spectroscopy

We present here the detailed measurements of radial distribution of the magnetic field in a gas-puff z-pinch plasma at the final stages of the implosion phase and at stagnation. While the measurements are chordal, the radial distribution of different charge states was utilized to measure the magnetic field locally for certain radii, so that unlike chordal measurements in general, the magnetic field radial distribution was obtained with no need for the Abel inversion of the data. The distribution was measured using the Zeeman effect via a novel spectroscopic technique, at several axial locations, and demonstrates striking features such as the peak field remaining at a radius much larger than the stagnation radius at all times. Furthermore, while the distribution observed is sometimes monotonic with respect to the radius, it is often not, a behavior that can be linked to 2D features in the plasma column resulting from the Rayleigh–Taylor instability. The current flowing through the stagnating plasma was found to be a small fraction of the total current, resulting in clearly insufficient magnetic pressure to balance the plasma pressure at stagnation. The magnetic field data, taken over several axial positions, are used to obtain the true inductance in the imploding plasma for the first time; it is found that the data cannot explain the current turnover at stagnation. A simulation with the MACH2-Tabular Collisional-Radiative Equilibrium magnetohydrodynamics code in the r–z plane shows that the peak of the magnetic field pinches to a much smaller radius than is observed in the spectroscopic data. Furthermore, the depth of the computed current turnover at stagnation is smaller than the measured one. The two observed features of a radially extended magnetic field at stagnation together with a deep current turnover are a challenge to match in simulations. Various calculations and estimates of the inductive and resistive load voltages are examined to ascertain if they are responsible for the observed current notch. Lastly, the results demonstrate that the knowledge of the true inductance in the driven load requires such magnetic-field-distribution measurements and that imaging data or electrical measurements are insufficient.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

A Novel Approach for Bounding the Stress Experienced by the Core of Utility-Scale Printed Circuit Heat Exchangers Under Thermohydraulic Loads

Printed circuit heat exchangers (PCHEs) have rapidly gained popularity since being introduced nearly three decades ago, and they are currently widely deployed in the petrochemical and aviation industry. Their compactness, thermohydraulic efficiency, inherent suitability for high temperature/pressure fluids containment, and demonstrated durability are some of the reasons the nuclear industry is seeking to adopt this technology as well. However, the relatively strict nuclear-related regulatory design code, especially when classified as critical to the safety of the reactors, is posing challenges to adopting the technology. From stress analysis point of view, one undesirable feature of PCHEs is their geometrical complexity, which is implied by their multilength-scale features. As a result, a full-scale model of a utility-scale exchanger cannot simply be solved on a computer because meshing such components results in a vast number of degrees-of-freedom. Here, this work seeks to address the challenge of stress analyses to PCHEs by presenting a method to simplify the geometry of PCHE designs. The models proposed by this work can be practically analyzed on a standard computer and provide a path for implementing ASME design rules. The analyses presented herein are divided into five separate investigations. Each is carried out to incrementally simplify the analyzed model by addressing features such as the shapes of the flow passages, the complex distribution of stress in large components, the three-dimensionality of the stress and strain, the thermal stresses caused by thermohydraulic operation observed experimentally and more.

42 ENGINEERING↗

Characteristics of Fluid‐Solid Interaction Constitutive Models Within Poroelastodynamics at Higher Strain‐Rates and Large Deformations Implemented in 1D

The large deformation, mixed formulation, finite element (FE) modeling approach presented in Irwin et al. 2024 is extended herein to include improved constitutive models for representing dynamic solid-fluid interactions at higher strain rates (𝒪⁢(1⁢0 2 −1⁢0 3 )⁢s −1 ) and larger overpressure magnitudes (𝒪⁡(1⁢0 2 )⁢kPa) within a biphasic soft porous material using Theory of Porous Media (TPM) at finite strain. Specifically, these constitutive modeling improvements are the following: (i) a more physically robust constitutive model for pore fluid seepage velocity via inclusion of pore fluid viscous stress, and (ii) a modified deformation-dependent-permeability model and updated hyperelastic constitutive model better suited for handling larger volumetric compressions and extensions. The novelty of the present work is mainly the contribution (i): inclusion of pore fluid viscous stress at higher strain-rate and large deformations, which requires 𝐶 1 continuity in the weak formulation, accomplished by employing Hermite cubic interpolation functions within a mixed nonlinear poromechanical finite element formulation. In (ii), the model is updated to weakly enforce solid phase incompressibility, such that this assumption is not violated numerically, which provides improved numerical stability for achieving larger overpressure magnitudes on 𝒪⁡(1⁢0 2 ) kPa, which were not achievable with the previous Kozeny–Carman model in Irwin et al. 2024. Also in (ii), the volumetric part of the solid skeleton free energy function is modified to ensure proper bounds on the solid skeleton Jacobian of deformation 𝐽 s related to incompressibility of the solid phase. Uniaxial strain, unidirectional flow examples at higher strain rates (𝒪⁢(1⁢0 2 −1⁢0 3 )⁢s −1 ) and larger deformations (up to 0.2 (or 20%) nominal axial strain) demonstrate the improved physical representation—and numerical stability—of these constitutive model improvements.

42 ENGINEERING↗

Hot, cold, or just right? An infrared biometric sensor to improve occupant comfort and reduce overcooling in buildings via closed-loop control

To improve occupant comfort and save energy in buildings, we have developed a closed-loop air conditioning (AC) sensor-controller that predicts occupant thermal sensation from the thermographic measurement of skin temperature distribution, then uses this information to reduce overcooling (cooling-energy overuse that discomforts occupants) by regulating AC output. Taking measures to protect privacy, it combines thermal-infrared (TIR) and color (visible spectrum) cameras with machine vision to measure the skin-surface temperature profile. Since the human thermoregulation system uses skin blood flow to maintain thermoneutrality, the distribution of skin temperature can be used to predict warm, neutral, and cool thermal states. We conducted a series of human-subject thermal-sensation trials in cold-to-hot environments, measuring skin temperatures and recording thermal sensation votes. We then trained random-forest classification machine-learning models (classifiers) to estimate thermal sensation from skin temperatures or skin-temperature differences. The estimated thermal sensation was input to a proportional integral (PI) control algorithm for the AC, targeting a sensation level between neutral and warm. Our sensor-controller includes a sensor assembly, server software, and client software. The server software orients the cameras and transmits images to the client software, which in turn assesses occupant skin temperature distribution, estimates occupant thermal sensation, and controls AC operation. A demonstration conducted in a conference room in an office building near Houston, TX showed that our system reduced overcooling, decreasing AC load by 42% when the room was occupied while improving occupant comfort (fraction of “comfortable” votes) by 15 percentage points.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

Proton beam power limits for stationary water-cooled tungsten target with different cladding materials

The proton beam power limit for a solid-tungsten spallation target is largely determined by beam induced thermomechanical structural loads and decay heat power deposition, while its lifetime is limited by radiation damage and fatigue life of the target materials. In this paper, we studied the power limits of a stationary water-cooled solid tungsten target concept. Tantalum clad tungsten was considered as a reference case. Being a low activation material, zircaloy 2 cladding option was studied and its decay heat driven power limit was compared with the reference case. Zirconium alloys have proven operations records in spallation target and nuclear fission environments, supported by materials data obtained from post irradiation examinations. Recent study also demonstrated feasibility of diffusion bonding zirconium to tungsten using vanadium foil inter layer. Particle transport simulations code FLUKA was used to calculate energy deposition and decay heat power deposition in the target, based on the beam parameters technically feasible at the Second Target Station of the Spallation Neutron Source at Oak Ridge National Laboratory. The energy deposition data were used for flow, thermal, and structural analyses to determine the beam intensity limit on the target concept studied. The decay heat deposition data were used to calculate the transient temperature evolution in the tungsten volumes in a loss of coolant accident (LOCA) scenario to determine its beam power limit. For a 1.3 GeV proton beam, the power limit on a stationary target was 400 kW for a tantalum clad target model and 800 kW for a zircaloy 2 clad target model.

Lee, Yong Joong↗

Effects of canopy midstory management and fuel moisture on wildfire behavior

Abstract Increasing trends in wildfire severity can partly be attributed to fire exclusion in the past century which led to higher fuel accumulation. Mechanical thinning and prescribed burns are effective techniques to manage fuel loads and to establish a higher degree of control over future fire risk, while restoring fire prone landscapes to their natural states of succession. However, given the complexity of interactions between fine scale fuel heterogeneity and wind, it is difficult to assess the success of thinning operations and prescribed burns. The present work addresses this issue systematically by simulating a simple fire line and propagating through a vegetative environment where the midstory has been cleared in different degrees, leading to a canopy with almost no midstory, another with a sparse midstory and another with a dense midstory. The simulations are conducted for these three canopies under two different conditions, where the fuel moisture is high and where it is low. These six sets of simulations show widely different fire behavior, in terms of fire intensity, spread rate and consumption. To understand the physical mechanisms that lead to these differences, detailed analyses are conducted to look at wind patterns, mean flow and turbulent fluxes of momentum and energy. The analyses also lead to improved understanding of processes leading to high intensity crowning behavior in presence of a dense midstory. Moreover, this work highlights the importance of considering fine scale fuel heterogeneity, seasonality, wind effects and the associated fire-canopy-atmosphere interactions while considering prescribed burns and forest management operations.

54 ENVIRONMENTAL SCIENCES↗

TRUST Sensors in Environments: Thermocouples (SE-TC), Release FY25

The Delivery Environments Testbeds to Reduce Uncertainty in Simulations and Tests (TRUST) project is a broad project intended to analyze simplified problems experimentally and with modeling and simulation. The purpose of analyzing these simplified problems is to extend solution methods to more complex problems, as well as understand deficiencies and gaps in knowledge of methods currently used in more complex analyses. The TRUST project encompasses several smaller testbeds intended to isolate individual phenomena. The testbed under consideration in this report is the Sensors in Environments: Thermocouples testbed. In previous years, the purpose of this testbed was to quantify uncertainty of thermocouple sensors. To accomplish this, an aluminum plate was placed in a thermal chamber and subject to various types of thermal loading. Thermocouples were placed in various locations on the aluminum plate in various configurations (e.g., embedded in the plate, placed under Kapton tape), and an effort was made to quantify uncertainty in these measurements. Finite element simulations were performed to investigate how sensitive these measurements were to parameters such as the boundary conditions on the plate and material properties. However, a fundamental source of uncertainty in this analysis was the convective heat transfer from the plate. Convective heat transfer is a complex physical phenomenon comprised of a number of interacting sub-processes and is difficult to predict accurately a priori. As such, the main purpose of this testbed in FY25 was to better understand, both experimentally and numerically, the convective heat transfer from the plate. This is a highly applicable problem to several more complex problems, as convective heat transfer occurs in nearly all problems where a body is moving through air. Numerically, this required a two-step approach. First, the air flow in the thermal chamber was in vestigated using computational fluid dynamics. The commercial solver Fluent was used to perform these simulations. From these simulations, a heat transfer coefficient over the surface of the plate was calculated. This heat transfer was then used as boundary conditions for finite element heat transfer simulations within the plate, which were performed using Abaqus. Significant effort was devoted to automating the handoff between these two solvers. Experimentally, previous thermocouple results in the plate were used to validate the time-dependent thermal profiles produced from Abaqus. Further experimental efforts were performed both to help validate the Fluent simulations and to inform its boundary conditions. For example, hot-wire anemometers were used to measure the velocity in the chamber, which would be particularly useful in understanding the chamber inlet velocity. Thermocouple measurements were also taken in the chamber, instead of only on the plate, to serve as validation evidence for the Fluent simulations. Numerical results showed that the Fluent to Abaqus workflow matched previous plate thermocouple measurements well. This type of handoff is useful for more complex experiments, or those that are not able to be examined in as great of detail as this testbed, as it was performed without any experimental input. Experimental results, however, were more mixed. The anemometers proved unreliable, with inconsistent measurements across all anemometers, even at locations that were nearly identical. On the other hand, the thermocouples provided a relatively rich view of the temperature field in the chamber.

42 ENGINEERING↗

Filament extension atomizers: A study of integration potential into spray dryers for high-viscosity and non-Newtonian fluids

Filament extension atomizers (FEAs) are an emerging class of spray nozzles designed to atomize high-viscosity and non-Newtonian fluids that are challenging for conventional pressure or two-fluid nozzles. In this study, we investigate the influence of roller geometry, surface velocity, roller material, and fluid rheology on the atomization performance of FEA systems using concentrated whey protein suspensions (50–70 wt. %). Extensional and shear rheology experiments, along with high-speed imaging and particle image velocimetry, reveal that filament breakup dynamics are governed by competition between inertial, capillary, and viscoelastic stresses. High roller rotational velocity leads to narrower spray cones, contradicting rheology experiments and suggests a significant inertial contribution to filament breakup. Smaller rollers operating at the same rotational velocity led to broader spray cones consistent with expectations. An FEA nozzle was integrated into a conventional dryer producing particles in the 100 μm range. Results suggest that FEA technology enables atomization of highly viscous fluids at industrially relevant spray cone angles similar to those generated by pressure nozzles, offering a pathway to improve energy efficiency in spray drying by enabling higher solids loading feedstocks. Furthermore, these insights provide critical guidance for optimizing FEA nozzle designs and process parameters across a range of applications.

Energy efficiency↗

Demand Response Optimization and Management System for Real-TIme (DROMS-RT)

To design and demonstrate DROMS-RT, a highly distributed Demand Response Optimization and Management System for Real-Time (DROMS-RT) power flow control to support large scale integration of distributed renewable generation into the grid. AutoGrid developed a novel control and communications platform to allow highly dispatchable demand response (DR) services in time frames suitable for providing ancillary services to the transmission grid. These services will be substantially less expensive and more efficient than other forms of ancillary services options currently available to manage the intermittency associated with large-scale renewable integration. DROMSRT successfully leveraged Automated Demand Response (ADR) by fundamentally re-thinking the architecture of the DR platform from the ground up and by developing innovative new technologies in a number of areas related to DR. DROMS-RT leveraged the low-cost, open, interoperable DR signaling technology, OpenADR, and low-cost, internet-protocol based telemetry solutions to reduce the cost of hardware. This allowed DROMS-RT to provide dynamic price signals to millions of OpenADR clients. Statistically rigorous signal processing techniques were developed to reliably detect even small load reductions in the presence of noisy baseline profiles. Novel forecasting engines based on modern online machine learning algorithms enabled accurate individualized forecasts for customer loads in the presence of dynamic pricing signals, and a real-time decision engines enabled continuous optimization and optimal dispatch of DR resources across a large portfolio of heterogeneous loads that respond at varying time-scales. Moreover, the real-time optimization conducted by the decision engine can utilize grid physics to maximize load reduction at the transmission system in addition to the distribution sites, for more efficient grid operation. Finally, the Software-as-a-Service (SaaS) availability of the DROMS-RT platform has reduced the cost of deployment and enable participation of small commercial and residential customers in DR who otherwise would not be able to do so.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Involute Working Group – FSI Analysis of Fuel Plates Using Finite Volume and Finite Element Methods

The three involute plate research reactors RHF, HFIR, and FRM II have expressed an interest in using computational software to carry their steady-state safety analysis. Since these tools represent a significant departure from the methods used currently (one-dimensional), the acceptability of the new approach by regulators requires thorough verification and validation of these tools. Therefore, Argonne National Laboratory and the three involute-plate reactors formed an informal group called the Involute Working Group aiming at qualifying computational tools to perform steady-state safety analysis. The present report focuses on a comparison of finite volume and finite element methods to model solids in fluid-structure interaction problems with the goal to estimate the coolant flow-induced fuel plate deflections obtained with the two methods. The finite volume method will be obsoleted in STARCCM+ by the end of 2021, nevertheless, this evaluation is important because the method was used by ANL researchers to model the response of the fuel plates, despite its drawbacks, which are discussed in the report. It was essential to check how those estimates compare to the results obtained with the finite element method that is considered superior for structural analysis. Various geometries, i.e., flat, cylindrical and circle-involute fuel plates, as well as coolant flow speed, were considered. The comparison shows that, independently of the plate geometry, the finite volume method significantly underestimates the deflection as compared to finite element method for coarser meshes. When the discretization is developed as a result of a mesh sensitivity study using finite element method, the result obtained using finite volume method can be a few times smaller than the corresponding finite element method solution. A code-to-code comparison , between STAR-CCM+ and LS-DYNA was included in the analysis. Within the LS-DYNA models, two types of finite element formulations were used: solid and shell finite elements. Mesh sensitivity study showed that both approaches converge to a similar value that was obtained with STAR-CCM+ finite element solver. The evaluation of the computational solvers was extended by adding two benchmark cases from the STAR-CCM+ Verification Suite and presented in the Appendix A. The selected cases are: (1) bending of a cantilever beam under external load, and (2) cylindrical shell deformation analysis, known in the literature as ‘Scordelis-Lo roof’. The problems were solved with finite volume, and finite element methods, and the results confirmed the previously discussed findings. The analysis shows that the finite element solver is superior to the finite volume solver in terms of representation of model geometry and estimating the structural behavior of fuel plates. Depending on the ratio of the load to the flexibility of the plate, the finite volume solver can greatly under- or overestimate the structural response if a very carefully selected mesh is not used.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

In-Field Testing of Components for Feedback and Control of the ITER Disruption Mitigation System

Here, the shattered pellet injection (SPI) method has been chosen as the disruption mitigation system (DMS) for ITER. To protect the device from plasma disruptions that cause damaging heat and electromagnetic loads, SPI is used to inject high-Z material into the plasma. The process of SPI utilizes cryogenic cooling to form solid pellets. Pellets are accelerated down a barrel and into an angled surface, causing the pellet to shatter prior to entering the tokamak chamber. For the DMS to function reliably, the 27 separate shattered pellet injectors planned for ITER must rely on many components to provide accurate feedback data and for control functions. Each component in the DMS is exposed to an elevated background magnetic field depending on its placement and proximity to the plasma chamber. A Helmholtz coil test stand that is operated at Oak Ridge National Laboratory was utilized to test the components in relevant background field levels to assess component performance. This paper details the test design and results for in-field component operation for a variety of components. This list includes the following components: two different network switches for camera connectivity, a VAT fast shutter valve intended to reduce the flow of SPI propellant gas into the torus, a solenoid control valve intended for use in the pellet formation process, pressure/vacuum switches to be used for feedback and control, a printed circuit board piezo pressure sensor to be used to measure breech pressure, and various relays for the high-voltage pulsed power supply used to drive the SPI propellant valve.

Disruption mitigation↗

Aerodynamic characterization of two tandem wind turbines under yaw misalignment control using actuator line model

Yaw control has proven to be promising in alleviating the wake effects that plague the efficiency of wind farms. In this work, the actuator line modeling (ALM) method is adopted to simulate the flows over two tandem turbines distanced by 3–7 rotor diameters, with the yaw angle of the upstream rotor varying from γ 1 = 0° to 50°. The aim is to provide a comprehensive aerodynamic characterization of this simple wind farm under yaw misalignment control. With increasing yaw angle, the power generated by the downstream rotor increases, compensating the power loss in the upstream rotor, and resulting in significantly higher total power of the two turbines than that without yaw control. The maximum power output is achieved as the upstream wake of the yawed rotor is redirected away from the downstream rotor plane. Behind the downstream rotor, the secondary steering phenomenon is observed, where the wake is also redirected from the centerline. The use of the actuator line model also reveal unsteady aerodynamic characteristics that cannot be captured by lower-fidelity models. For the upstream rotor, the yaw misalignment results in time-varying change in the local angle of attack on the blade, giving rise to unsteady loading. The downstream rotor is partially submerged in the deflected wake incurred by the yawed upstream rotor. As the blade revolves into and out of the wake deficit, the blade experiences cyclic loading, leading to even stronger fluctuations in the aerodynamic loads than the upstream rotor. These analysis provides a comprehensive understanding of the yaw control effects on the two tandem rotors from the perspectives of aerodynamic performance, wake profiles, and unsteady characteristics. In conclusion, the insights gained from the present study can aid the design of collective yaw control strategies of wind farms, and lay the foundation for assessing the fatigue damage associated with yaw misalignment.

17 WIND ENERGY↗

Investigations of correlation and coherence in turbulence from a large-eddy simulation

Abstract. Microscale flow descriptions are often given in terms of mean quantities, turbulent kinetic energy, and/or stresses. Those metrics, while valuable, give limited information about turbulent eddies and coherent turbulent structures. This work investigates the structure of an atmospheric boundary layer using coherence and correlation in space and time with a range of separation distances. We calculate spatial correlations over entire planes of velocity fluctuations, from which we can evaluate the correlation along different directions at different spacings. Similarly, coherence of the three velocity components over separations in the three directions is also investigated. We apply these analyses to a mesoscale–microscale coupled scenario with time-varying conditions and examine nuances in spatial correlations that are often overlooked. Through these analyses and results, this work highlights important differences observed in terms of coherence when comparing large-eddy simulation data to simpler models and suggests ways to improve these simpler models. We note that such differences are important for disciplines like wind energy structural dynamic analysis, in which blade loading and fatigue depend strongly on the structure of the turbulence. We emphasize the additional wealth of data that can be provided by typical atmospheric boundary layer large-eddy simulation when correlation and coherence analysis is included, and we also state the limitations of large-eddy simulation data, which inherently truncate the smaller scales of turbulence.

17 WIND ENERGY↗

Pressure Spike in The LBNF Absorber Core s Gun Drilled Cooling Channel from an Accident Beam Pulse

he LBNF absorber consists of thirteen 6061-T6 aluminum core blocks. The core blocks are water cooled with de-ionized (DI) water which becomes radioactive during beam operations. The cooling water flows through gun-drilled channels in the core blocks. The cooling water is supplied by the LBNF absorber radioactive water (RAW) cooling system which is designed as per ASME B31.3 normal fluid service. An uninhibited beam accident pulse striking the water channels was identified as a credible accident scenario. In this study, it is assumed that the beam pulse hits the absorber directly without interacting with any of the other upstream beamline components. The beam parameters used for the LBNF beam are 120 GeV, 2.4 MW with a 1.2 s cycle time. The accident pulse lasts for 10 µs. The maximum energy is deposited in the 3rd aluminum core block. For the sake of simplicity, it is assumed that the accident pulse strikes the 1 in. ID water channel directly. The analysis here simulates the pressure rise in the water during and after the beam pulse and its effects on the aluminum piping components that deliver water to the core blocks. The weld strengths as determined by the Load and Resistance Factor Design (LRDF) and the Allowable Strength Design (ASD) are compared to the forces generated in the weld owing to the pressure spike. A transient structural analysis was used to determine the equivalent membrane, peak, and bending stresses and they were compared to allowable limits.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Pressure Spike in The LBNF Absorber Core s Gun Drilled Cooling Channel from an Accident Beam Pulse

The LBNF Absorber consists of thirteen 6061-T6 aluminum core blocks. The core blocks are water cooled with de-ionized (DI) water which becomes radioactive during beam operations. The cooling water flows through gun-drilled channels in the core blocks. The cooling water is supplied by the LBNF Absorber Radioactive Water (RAW) cooling system which is designed as per ASME B31.3 Normal Fluid Service [1]. An uninhibited beam accident pulse striking the water channels was identified as a credible accident scenario. In this study, it is assumed that the beam pulse hits the Absorber directly without interacting with any of the other upstream beamline components. The beam parameters used for the LBNF beam are 120 GeV, 2.4 MW with a 1.2 s cycle time. The accident pulse lasts for 10 μs. The maximum energy is deposited in the 3rd aluminum core block. For the sake of simplicity, it is assumed that the accident pulse strikes the 1 in. ID water channel directly. The analysis here simulates the pressure rise in the water during and after the beam pulse and its effects on the aluminum piping components that deliver water to the core blocks. The weld strengths as determined by the Load and Resistance Factor Design (LRDF) and the Allowable Strength Design (ASD) are compared to the forces generated in the weld owing to the pressure spike. A transient structural analysis was used to determine the equivalent membrane, peak, and bending stresses and they were com-pared to allowable limits.

43 PARTICLE ACCELERATORS↗

Hydrogen Storage for Load-Following and Clean Power: Duct-firing of Hydrogen to Improve the Capacity Factor of NGCC Plants (Final Report, Phase II Pre-Front End Engineering Design Study)

GTI Energy (GTI) and team members Southern Company Services (SCS), Pacific Gas & Electric (PG&E) and the Electric Power Research Institute (EPRI) performed a Phase II Pre-Feed Study under contract DE-FE0032008 for Hydrogen Storage for Load-Following and Clean Power. The configuration of the proposed system consists of subsystems for on-site H 2 production, on-site H 2 storage (up to 54 MWth in commercial vessels) and H 2 combustion in a duct burner in a Heat Recovery Steam Generator (HRSG) integrated with an existing fossil asset. Here, the firing rate of the duct burner is varied to let the plant respond to fluctuations of electrical load and H 2 production is relatively constant by storing H 2 . The proposed system is an improvement over alternate low carbon dispatchable power options. Technoeconomic analyses show H 2 produced with GTI’s patented Compact Hydrogen Generator (CHG) with inherent carbon capture will be lower cost (Levelized Cost of Hydrogen, LCOH) relative to hydrogen produced with a Steam Methane Reformer (SMR) with an amine system for carbon capture. This lower cost hydrogen enables our integrated system to deliver electricity (Levelized Cost of Electricity, LCOE) at 17.4% lower cost relative to a SMR sourced H 2 -fired HRSG (with an amine system for carbon capture) – Steam Turbine Generator (STG). Our analysis using EPRI’s US REGEN macroeconomic model shows our system will have significant demand in the power market and therefore require significant capacity expansion (CHG plants built to deliver hydrogen) to deliver low cost, low carbon power. In Phase II, our Team has completed a detailed system definition including the development of process models, definition of battery limits, Process Flow Diagrams (PFDs), Piping and Instrumentation Drawings (P&IDs), and a plant layout. A preliminary design for the key components of the CHG was developed including component lists & specifications. An evaluation of environmental and permitting considerations was completed and included the development of an Environmental Information Volume (EIV). The duct burners, which are flexible and can burn hydrogen and/or natural gas, were defined and initial CFD analyses were completed.

03 NATURAL GAS↗

Geomagnetically Induced Current Field Test on Large Grid-Connected Power Transformers: Analysis, Model Development, and Simulations

Geomagnetic-induced current (GIC) flow in power grids can cause undesirable effects such as transformer overheating, harmonics, higher reactive power demand, etc. Many simulation models have been developed to study these effects, but real-world verification on modern transformer designs is rare. Here, this paper presents the first long-duration GIC field test in the U.S. performed on high-voltage, grid-connected transformers featuring winding clamps and tie rods instead of conventional tie bars. Field measurements were taken to evaluate GIC effects. These measurements also aided in developing and validating thermal and electromagnetic transient (EMT) models of the transformers. During the test, significant current and voltage distortions were observed along with considerable transformer reactive power losses. Analysis of the field measurements showed that the transformers’ hottest spot was at the inner windings, and their k-factors were close to factory test and software default values. Thermal simulations indicated that the transformers would not violate their thermal limits even for a GIC waveform that peaks at about 200 A/phase. EMT simulations revealed that increased transformer loading may reduce GIC-induced reactive power demand and harmonics in certain scenarios. The study also highlighted potential inaccuracies in using the k-factor method to calculate transformer reactive power losses.

EMTDC↗

Multi-X-point radiation and its dynamics in up/down asymmetry in plasma detachment regimes in Wendelstein 7-X

In Wendelstein 7-X (W7-X), stable plasma detachment in ECRH hydrogen plasmas is routinely achieved with the standard magnetic configuration, containing five magnetic islands in the scrape-off layer (SOL). Plasma detachment induced by intrinsic carbon impurities is characterized by a high radiation fraction (ƒ rad ) with significantly reduced divertor heat load and particle flux. In this work, we present the radiation dynamics during the detached plasma (DP) phase, such as the shift of radiation zones in radial and poloidal directions towards the X-points as ƒ rad increases (∼0.6–0.9), as well as the penetration of the radiation layer into the confinement region in the deep DP phase with ƒ rad approaching unity. In particular, a structure of multi-X-point radiation (multi-XPR) with an up/down asymmetry in the DP phase is highlighted, which is revealed by bolometer tomography and is further confirmed by video diagnostics. The multi-XPR structure forms helical 3D bands aligned with W7-X field periodicity. Field reversal experiments show that the brightest XPR displaces between the upper and lower SOL regions, suggesting the E× B drift effects. This paper presents the multi-XPR structure observed in the W7-X plasma in detail for the first time. A simplified model considering the influence of the poloidal E× B drift ( V d ) on the impurity flow in the SOL shows that the downstream drift toward the target or target shadow region ( V d < 0) decreases the impurity content, while the upstream drift toward the LCFS ( V d > 0) increases the impurity content. On this basis, the poloidal drift potentially leads to an up/down asymmetry of impurity density in the SOL despite the symmetry magnetic topology. The dynamics of the up/down asymmetry in the multi-XPR structure is also related to the magnitude V d /D (normalized to the impurity diffusivity), with an additional effect owing to the radial inward shift of the emission zone. These results provide new insights into impurity-induced detachment dynamics and provide a basis for improving 3D modeling of impurity transport, which typically does not consider drifts.

E × B drift↗