Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “flow conservation”

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

Liquid Salt Combined-Cycle Pilot Plant Design

The work described in this report is responsive to the Office of Fossil Energy program ‘Energy Storage for Fossil Power Generation.’ This Phase I report has been prepared by Pintail Power LLC, with support from Nexant ECA, Electric Power Research Institute (EPRI) and Southern Company Services as a deliverable for the U.S. Department of Energy for NETL Award DE-FE-00320016. The Liquid Salt Combined Cycle™ (LSCC™) technology provides large-scale energy storage integrated with Fossil Electric Generating Units (FEGUs) to meet critical needs in the energy transition by providing: • the lowest cost large-scale storage for time-shifting of renewable energy, • superior fuel efficiency to reduce GHGs from dispatchable resources, • flexible capacity and ramping to balance variability of wind and solar resources, • essential grid stability services to assure reliability of a low-carbon grid. The LSCC approach: • employs equipment that has already been proven in utility service, • uses safe, non-toxic, non-degrading, perpetual-life storage medium, • leverages and repurposes existing FEGU assets, • expands the value stack of energy storage to reduce market, financing, and commodity risks. Pintail Power has developed the LSCC technology to meet the need for reliable, efficient, and cost-effective integration of Variable Renewable Energy (VRE) into a low-carbon electric grid by coupling proven thermal energy storage with proven gas turbines, steam turbines, and heat transfer equipment. This novel approach is intended to address the key issues facing the grid and operators of renewable and fossil generating units including: • Overgeneration and curtailment of renewables, • Need for fast ramping dispatchable resources, • Improved efficiency and flexibility of fossil units, • Additional peaking capacity to support electrification of transportation and heating, • Provision of reliability services to support high penetration of VRE, especially synchronous inertia and fast frequency response. A Technology Readiness assessment by EPRI confirmed that LSCC technology consists of commercially proven hardware used in industrial and utility applications. Although the novel LSCC approach has not yet been demonstrated as a complete system, interfaces between major components have been conservatively specified. A Phase III pilot is planned to demonstrate equipment integration and operation. The patented innovation is removal of the evaporator section from the exhaust heat recovery system, with the evaporation performed by stored energy in a separate steam generator. This arrangement couples renewable and fossil power generation via long-duration energy storage to deliver cost, performance, and operational synergies, including superior charging and discharging flexibility, reduced fuel consumption and lower CO 2 emissions compared to conventional Combined Cycle Power Plants, and low-cost, large-scale energy storage. The LSCC technology is composed of proven equipment integrated with gas turbine exhaust heat in a novel system. During charging, electric heaters raise the salt temperature as it flows from the Cold Salt Tank to the Hot Salt Tank. During discharging, hot salt produces steam from feedwater that is heated with gas turbine exhaust, which also superheats steam to drive a steam turbine. LSCC technology can be added to any combustion-turbine to integrate renewable energy, provide needed grid services, and increase the value of fossil electric generating units based on the technology’s following attributes: • Long-duration storage enables time-shifting of VRE to avoid curtailment and impairment of renewable assets. • Long storage duration combined with fast-charging capability increases arbitrage opportunities by storing more energy when the price is low and discharging more hours when the price is high. • Long storage duration allows resource adequacy to be supplied across multiple days to increase reliability and reduce risk. • The stored energy reduces fuel heat rate and GHG emissions, and increases merit, so the LSCC dispatches earlier and longer to increase the plant’s capacity factor and asset value. • The stored energy enables pre-heating and startup of the steam cycle, without operating the gas turbine, to enable fast startup and ramping when dispatched for discharge. • The steam turbine can operate without the gas turbine so it can provide valuable synchronous inertia during charging without consuming fuel. • Fast frequency response and regulation services can be provided during charging using solid-state heater and pump controls to vary the charge power input in response to grid signals. • The LSCC system can be configured for resilience including black start, islanded/micro-grid operation, and even self-recharging of storage using either gas turbine power or gas turbine exhaust heat. The commercialization plan is to add LSCC technology to existing simple cycle gas turbine power plants with the 50MW GE LM6000 aero-derivative gas turbine as the reference design basis. A Techno-economic assessment of the reference design evaluated the benefits (Levelized Avoided Cost of Energy) and costs (Levelized Cost of Energy). The plant definition included all major systems and budgetary vendor quotes. Pintail Power and NexantECA developed the overall cost estimate for the LSCC plant up to the total plant cost level, following the DOE-NETL cost estimate guidelines at AACE Class 3 (-20%/+30%). This includes the equipment cost, bulk material, direct and indirect labor costs to arrive at the bare erected cost. Engineering costs are factored from the BEC and added to it to arrive at the EPC cost. Process and project contingencies were then factored from the EPC cost and rolled-up to yield the total plant cost of $\$$184 million for 1746 MWh of discharge electricity. • At $\$$105/kWh, the reference plant costs less than any of the Energy Storage Systems evaluated by PNNL in 2020 for the Energy Storage Grand Challenge. Operations and Maintenance cost estimates were scaled from combined cycle practice, assuming that the LSCC unit was co-located with and sharing some labor expense with other units, to arrive at $\$$2.2 million per year. Plant economics were evaluated using prices from the ERCOT Day-Ahead Market for calendar year 2019 (excluding the market disruptions from the COVID pandemic and the February 2020 deep freeze event). Assuming economic dispatch in the ERCOT Day-Ahead market, the reference plant capacity factor would have discharged for 2777 hours at 91.9 MW, a 31.66% capacity factor, with a marginal cost of $\$$25.59/MWh, and a LACE of $\$$82.41/MWh. Fixed charges were calculated according to EIA guidelines to arrive at an LCOE of $\$$83.48. The benefit-to-cost ratio of 0.99 suggests that the reference plant would have been cost-effective and competitive in the market. EPRI interviewed selected utilities to gauge the need for, applicability of and interest in the LSCC system. Several utilities are currently managing increased load growth along with the inclusion of increasing levels of renewable generation, putting pressure on conventional generation by requiring increased turndown requirements and ultimately lower capacity factors. All of the utilities interviewed have CO 2 reduction targets in the 2030-2050 timeframe that will severely limit the participation of fossil generation and require better utilization of carbon free generation. While there is limited opportunity for storage in the current markets, the utilities interviewed stated that there will be a substantial need for long duration energy storage in the future given the expected trends. Utilizing an energy storage system will generally be preferred over new gas capacity in some cases, with the capabilities of the LSCC system being a potential option for retrofit to existing simple cycle gas turbine units, allowing them to deliver greater participation in the market with lower carbon intensity. A technology gap assessment and technology maturation plan identified a pilot-scale demonstration as the final step before commercialization. Key gaps to be addressed during the Phase II FEED (Front-End Engineering Design) are component selection and design, commissioning procedures, and operational procedures and the control system for LSCC charging and discharging. The project team has been expanded to include Wood Group PLC as EPC. The proposed Phase II work leads to a pilot-scale engineering demonstration (TRL 6) to be conducted at Southern Company’s Plant Rowan, where the prototype system will perform “all the functions that will be required of the operational system.” The proposed pilot will facilitate commercialization (TRL-9) by scale-up to utility-scale systems integrated with peaking GTs or directly to facility scale systems using industrial GTs. The conceptual design for the pilot plant focuses on the novel integration aspects of LSCC technology. A slipstream of gas turbine exhaust will feed a waste heat recovery unit coupled to a molten salt steam generator heated by stored energy. The pilot is intended to demonstrate all key operating modes of the LSCC technology during charging, discharging and standby. The pilot equipment will be approximately one-seventh scale of the LM6000 commercial target and is expected to have commercial off-ramp potential for facility-scale applications.

01 COAL, LIGNITE, AND PEAT↗
DOE OSTI2024

Measuring Economy-Wide Circularity of the United States: An Input-Output Model in Mass Units

The ideal of creating closed material cycles by transforming the way we make, use and repurpose goods has become known as the circular economy. Besides its visionary appeal, material efficiency strategies central to circular economy can allow us to meet decarbonization goals that are otherwise out of reach. Production of basic materials such as cement, iron and steel and petrochemicals is one of the largest drivers of greenhouse gas emissions. Measuring circularity, and understanding its relationship with other sustainability metrics, however, is still difficult due to lack of data in mass units covering the entire economic system. Input-output (I-O) tables were originally created as a means of tracking the monetary flows that represent exchanges of goods and services within an economy, but have found additional uses in life cycle assessment and material flow accounts. In the traditional approach of environmentally extended I-O tables, emissions and energy data augment monetary flows. This approach can lead to price effects that distort physical quantities. Also, circular economy strategies and their goals relate to mass, not monetary flows. Therefore, it is best to simulate them using physical quantities. With these issues in mind, we have developed I-O tables in mass units to measure the flow of goods in the U.S. economy. This allows us to better measure how circular the economy really is, and how policy changes to affect this circularity may also affect decarbonization goals. Improving knowledge of these linkages could allow manufacturers to make changes in their production to achieve their sustainability and decarbonization targets. Our tool also provides a standard approach and public repository for data in physical units, making such data more available, useful, and meaningful. Following an established set of material flow metrics used by the European Union, we have used our tool to calculate material footprints over time differentiated by oil and gas versus other extractive industries. This approach also shows the relative trade balances of the U.S. in materials. We have developed a case study for the iron and steel sector, showing how different decarbonization scenarios affect not only economy-wide greenhouse gas emissions, but also total material use.

circular economy↗
DOE OSTI2023

Adjustments to the law of the wall above an Amazon forest explained by a spectral link

Modification to the law of the wall represented by a dimensionless correction function ϕ RSL (z/h) is derived using atmospheric turbulence measurements collected at two sites in the Amazon in near-neutral stratification, where z is the distance from the forest floor and h is the mean canopy height. The sites are the Amazon Tall Tower Observatory for z/h∈ [1,2.3] and the Green Ocean Amazon (GoAmazon) site for z/h∈ [1,1.4]. Here, a link between the vertical velocity spectrum E ww (k) (k is the longitudinal wavenumber) and ϕ RSL is then established using a co-spectral budget (CSB) model interpreted by the moving-equilibrium hypothesis. The key finding is that ϕ RSL is determined by the ratio of two turbulent viscosities and is given as ν t,BL /ν t,RSL, where ν t,RSL = (1/A)∫$^{∞}_{0}$ τ(k)E ww (k)dk, ν t,BL = k v (z−d)u * , τ(k) is a scale-dependent decorrelation time scale between velocity components, A = C R /(1−C I ) = 4.5 is predicted from the Rotta constant C R = 1.8, and the isotropization of production constant C I = 3/5 given by rapid distortion theory, k v is the von Kármán constant, u * is the friction velocity at the canopy top, and d is the zero-plane displacement. Because the transfer of energy across scales is conserved in E ww (k) and is determined by the turbulent kinetic energy dissipation rate (ε), the CSB model also predicts that ϕ RSL scales with L BL /L d , where L BL is the length scale of attached eddies to z = d, and L d = u$_{*}^{3}$/ε is a macro-scale dissipation length.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Improving AHU Performance by Minimizing Approach Temperature, Reducing Air Maldistribution, and Efficiently Handling Sensible and Latent Loads (Phase 1 Interim Final Technical Report)

Residential air handling units (AHUs) have stayed the same in form and efficiency for the past 30+ years, with incremental improvements made to address safety, functionality, and energy-efficiency. The purpose of this research in Phase I, Topic 9a: Next Generation Residential Air Handlers, was to improve AHU performance by minimizing heat exchanger (HX) approach temperature, reducing air maldistribution, and developing alternative system configurations which more efficiently handle sensible and latent loads. In this research Optimized Thermal Systems (OTS) developed, modeled, and evaluated multiple alternative system concepts. A dual vapor compression system separate sensible and latent cooling (SSLC) concept was studied to inform work on alternative concepts and to show best-case performance benefit. System concepts included ejector enhanced vapor compression cycles, desiccant assisted dehumidification, dual evaporator SSLC, and alternative AHU HX configurations. A dual vapor compression system showed COP improvement of 20%, however, required additional components, increased unit size, and increased cost. Two types of ejector enhanced vapor compression cycles with dual evaporators improved system COP by 9 to 11%, and reduced AHU losses by as much as 18%, with design changes limited to the AHU, no unit physical size increase, and a moderate increase to system first cost. Desiccant assisted air-conditioning required increased air flow rate resulting in higher fan power and the desiccant wheel increased sensible heat load leading to increased compressor power and reduced system COP. Dual evaporator cycles were found to degrade performance due to increased expansion losses. Optimized single slab HX designs used in place of the traditional A-coil HX led to 44–49% reduction in aluminum, 47–60% less refrigerant charge, and improved HX velocity distribution.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2024

Steam Generator Model Design Parameter Sensitivity Study Using Advanced Optimization Tools

This study focuses on design parameter sensitivity studies pertaining to a steam generator (SG) model, using both Python and machine-learning tools. The SG model is a mathematical representation (including fluid flow and heat transfer equations/models/correlations) of a steam-generating unit in a pressurized water reactor (PWR)-type small modular reactor (SMR) system. Design studies involve changing the model’s input design parameters (e.g., temperature, pressure, mass flow rate) to observe the resulting effects on the output of the system (e.g., heat transfer coefficient [HTC], Nusselt number, heat transfer performance). Sensitivity studies analyze the degree to which system output and/or desired parameters (e.g., HTC or heat transfer performance) are sensitive to changes in input parameters. By using machine-learning tools such as the Risk Analysis Virtual Environment (RAVEN) developed at Idaho National Laboratory (INL), detailed design parametric sensitivity studies and model optimization were performed. Six input parameters—namely, the pressure, temperature, and mass flow rate for the inlet of the primary-side (hot fluid) and secondary-side (cold fluid) of the SG—were randomly perturbed via RAVEN’s Monte Carlo Sampler module, using uniform distributions (±1% relative changes). The analysis results give valuable insights into SG system performance and optimization, and provide justification for researching optimized sensor placement to effectively monitor and obtain experimental data.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2024

Steam Generator Model Design Parameter Sensitivity Study Using Advanced Optimization Tools

This study focuses on design parameter sensitivity studies pertaining to a steam generator (SG) model, using both Python and machine-learning tools. The SG model is a mathematical representation (including fluid flow and heat transfer equations/models/correlations) of a steam-generating unit in a pressurized water reactor (PWR)-type small modular reactor (SMR) system. Design studies involve changing the model’s input design parameters (e.g., temperature, pressure, mass flow rate) to observe the resulting effects on the output of the system (e.g., heat transfer coefficient [HTC], Nusselt number, heat transfer performance). Sensitivity studies analyze the degree to which system output and/or desired parameters (e.g., HTC or heat transfer performance) are sensitive to changes in input parameters. By using machine-learning tools such as the Risk Analysis Virtual Environment (RAVEN) developed at Idaho National Laboratory (INL), detailed design parametric sensitivity studies and model optimization were performed. Six input parameters—namely, the pressure, temperature, and mass flow rate for the inlet of the primary-side (hot fluid) and secondary-side (cold fluid) of the SG—were randomly perturbed via RAVEN’s Monte Carlo Sampler module, using uniform distributions (±1% relative changes). The analysis results give valuable insights into SG system performance and optimization, and provide justification for researching optimized sensor placement to effectively monitor and obtain experimental data.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2022

Stochastic hydrodynamics and hydro-kinetics: Similarities and differences

The hydro-kinetic formalism has been used as a complementary approach to solving the Stochastic Differential Equations (SDE) corresponding to noisy hydrodynamics. The hydro-kinetic formalism consists of a deterministic set of relaxation type equations that tracks the evolution of 2-point correlation functions of stochastic hydrodynamic quantities. Hence they are comparatively easier to solve than the SDEs, which are computationally intensive and need to deal with arbitrarily large gradients. This work compares the two approaches for the propagation and diffusion of conserved charge fluctuations in the Bjorken hydrodynamic model. For white noise, the two approaches agree. For colored Catteneo noise, which is causal, the two approaches diverge. This is because white noise only induces two-point correlations, while Catteneo noise also induces higher-order correlations. Furthermore, this difference is quantified from the effects of causal evolution and influence from higher-order correlations induced by the Catteneo noise.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

An Innovative Circuit Breaker for Renewable Power Grids

A key step towards a carbon-neutral future could be reached through dispersed power grids, featuring networks of local-scale renewable energy and battery storage plants. To prevent these power grids from damaging themselves and their surroundings when electrical faults arise, they must be integrated with ‘circuit breakers’, which temporarily interrupt the current flowing through them. However, currently available circuit breakers cannot handle the medium-voltage direct currents best suited for these grids. Through an innovative new circuit breaker design, Steve Schmalz and his colleagues at Eaton Corporation, the Illinois Institute of Technology (IIT), Virginia Tech and the National Energy Technology Laboratory (NETL), hope that this challenge will soon be overcome.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2023

Terrain-Influenced Winds and Fire-Fire Interactions in Wildland Fire Simulations [Dissertation]

Ensemble-based approaches to prescribed fire planning cannot be supported by computational fluid dynamics based models like FIRETEC and the Wildland-Urban Interface Fire Dynamics Simulator (WFDS) because they are too computationally expensive and cannot leverage large eddy simulation approaches like CAWFE and WRF-SFIRE because they have too coarse of resolution. QUIC-Fire was developed to fill this gap but it cannot currently address complex terrain, that is typical for instance in the Western United States. This dissertation describes a variety of improvements made to QUIC-Fire and its various incorporated algorithms in an effort to make it a viable tool in simulating wildland and prescribed fires on terrain. The modifications made to QUIC-Fire are described in three chapters. The first chapter describes the extension of the diagnostic wind model QUIC-URB, the wind engine of QUIC-Fire, to a terrain-following coordinate system. The terraininfluenced winds it generates are analyzed and compared. In particular, this chapter presents the mathematical derivation of the wind solver leading to a linear system of equations that are solved through the successive over-relaxation method. The model is validated against a standard test used in previous works (the Askervein Hill) and against a new dataset from measurements in the Socorro Mountains, New Mexico. The terrain-following implementation captures the correct phenomenology for the isolated Askervein Hill, with a wind speed up at the top of the hill. The model agrees well with measurements on the upwind side of the peak, but overestimates speed-up on the downwind side of the hill. This is due to the inability of the model to generate flow separation and wake-eddy dynamics. On a common laptop, the divergence-free wind field is obtained in 6 s, making the solver appealing for coupled fire-atmosphere simulations. The Socorro Mountain is highly complex, with many cliff faces, peaks, and valleys. Although the model captures the magnitude and direction of inlet and outlet areas of the domain, it performs rather poorly in the valley region and in the regions near the steep cliffs. Hence, the model shows good agreement with data in areas of open sloped terrain but lacks in areas where flow separation and thermally driven effects may be present (neither effect is addressed in this work). In the second chapter the implementation of the terrain-following version of QUIC-URB into QUIC-Fire, and the necessary changes needed to include terrain are described. No changes to the underlying fire spread algorithm are made other than what is required to correctly account for the inclusion of terrain. Previously published FIRETEC results that use five different topographies that share the same centerline profile are compared to simulation results from the modified QUIC-Fire that use the same topographies and fuels. QUIC-Fire results show overall similar behaviors in terms of how the topographies affect fire shapes and trends in spread rates. Due to the terrain-following version of QUIC-URB being unable to generate flow separations at the crest of hills, fire spread rates in these regions across all topographies are over-predicted when compared to FIRETEC. Lateral fire growth shows similar trends with FIRETEC between topographies but does not capture the increase in spread due to a diagonal interface between grassland and forested fuel region of the domain. These results suggest that there are three algorithms within QUIC-Fire that could use improvement: how flame tilt angle is accounted for, the incorporation of non-local drag effects, and the inclusion of the wake-eddy parameterizations that are used in QUIC-URB. Lastly, the third chapter describes a modification to the initial guess used for the calculation of the QUIC-URB mass-conserved wind solution during fire simulations. The modification is aimed at improving fire-fire interactions in QUIC-Fire simulations. The modification consists of using the solution from the previous timestep as the starting point for the calculation of the solution for the next timestep. Fire-fire interactions is greatly improved by the change but a new source of error is introduced. Due to how plumes are modelled in QUIC-Fire the new solution contains errors where gaps in the plume structure are present. However, these errors are mostly limited to the upper atmosphere, where they do not affect fire behavior at the surface, and their magnitude isn’t significant enough to discount the amount of new fire phenomenology now captured in QUIC-Fire with the change.

58 GEOSCIENCES↗
DOE OSTI2022

Vehicle-to-manufacturing (V2M) system: A novel approach to improve energy demand flexibility for demand response towards sustainable manufacturing

The U.S. manufacturing sector accounts for 77% of the industrial energy consumption, but its participation in demand response (DR) programs is largely lagged behind. The limited flexibility in production scheduling under high capacity utilization and the lack of DR schemes that incorporate energy demand flexibility measures are deemed as major barriers. In this study, a framework of the interactive vehicles-to-manufacturing (V2M) energy sharing system is proposed to improve the energy demand flexibility of the aggregated system and enhance the DR effectiveness for manufacturers. The V2M-based DR scheme aims to reduce the energy cost by load shifting through joint production and energy sharing control, which can eventually promote manufacturing DR implementation even under high capacity utilization requirements and dynamic real-time electricity prices. The V2M system is modeled based on a discrete-Markov chain considering the complex interconnections among various manufacturing resources and multi-directional energy flows among manufacturing facilities, electric vehicles, and the power grid. Based on the system model, a mixed-integer nonlinear programming (MINLP) problem is formulated to identify the optimal DR scheme. The effectiveness of the proposed approach is validated through comparisons with traditional manufacturing DR schemes. Finally, the results show that a 2.1 to 6.5 times energy demand flexibility improvement and an additional 4.7% to 6.9% energy cost reduction can be achieved by the proposed approach.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2020

Flow and pressure estimators in a waste heat recovery system

An apparatus includes a pump circuit structured to receive pump data indicative of an operating characteristic of a pump feeding a fluid to a waste heat recovery (WHR) system; a flow circuit structured to receive valve position data indicative of a position of a valve downstream of the pump, estimate a flow rate of the fluid exiting the pump, and estimate the flow rate of the fluid exiting the valve; and a pressure circuit structured to receive pressure data indicative of the pressure of the fluid exiting the valve, estimate a change in pressure of the fluid across the WHR system, and determine a pressure of the fluid in a hot section of the WHR system based on the pressure of the fluid exiting the valve and the change in the pressure of the fluid across the WHR system.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2023

CCN data from TAMU TRACER campaign in the Houston TX region from July to September 2022

During TRACER, the Texas A&M Rapid Onsite Atmospheric Measurements Van (ROAM-V) was deployed to capture airmasses behind (maritime) and ahead (continental) of the passage of the sea-breeze front through Houston. On select sampling days, ROAM-V sampled in the morning/mid-day on the coast and then transited to a second inland site for the afternoon/evening. The suite of instruments deployed on ROAM-V included a Condensation Particle Counter (CPC; GRIMM Model 5.403 CPC), Scanning Mobility Particle Sizer (SMPS; TSI 3750 detector, TSI 3082 classifier, TSI 3088 neutralizer, TSI 3081A Differential Mobility Analyzer), Cloud Condensation Nuclei counter (Droplet Measurement Technologies CCN Counter), micro pulse lidar (Droplet Measurement Technologies Micro Pulse LiDAR (miniMPL)), and a Davis Rotating Uniform size-cut Monitor (DRUM; DRUMAir 4-DRUM). Before sampling at each location, the latitude and longitude were recorded using the GPS on the phone application “My Altitude”. Onboard the ROAM-V, aerosol samples are drawn through a shared isokinetic inlet at a flow rate ranging from 3.5 to 7.0 LPM. A portion of this flow is directed through a cyclone impactor (Brechtel, Inc. Model SCC 0.732) and 0.5 LPM is directed to the CCN. To calculate particle losses, we used a two-step method. First, the measured SMPS size distributions were used to calculate particle loss through the inlet during sampling. Second, the corrected SMPS data was used to calculate the average of the total losses per scan down the CCN line. Then, the correction was applied to the CCN data. This calculation was done separately for each deployment location due to changes in the measured size distributions between locations. Particle loss from diffusion (based on Kesten, 1991 and Gormley, 1949), inertial impaction in 90-degree bends (based on Aerosol Measurement, 2011 and Crane, 1977), and cyclone impactor efficiency (based on Dirgo, 1985) were included in the loss calculation. When the SMPS was not sampling at a location (in the case of an instrument malfunction or operator error), the reported CPC data was corrected with an average of the total losses for the entire campaign at the specified deployment location (e.g., if we needed to correct Galveston data, then the average of all calculated losses at Galveston was taken). These flatline corrections were used for all data on 22/07/13, 22/07/20, 22/07/22, and the data from Galveston on 22/08/09. The supersaturation uncertainty is estimated conservatively at +/- 0.03%, where variation in the inlet temperature, pressure, and calibration technique prevents a more accurate measurement. Confidence in the reported supersaturation measurements is based on a pre-campaign calibration (following the methods from our previous work and Deng, 2014 based on Rose, 2008) in addition to inter-comparisons with the DOE for two days (22/08/18 and 22/09/01) where TAMU was co-located with AMF1. The inter-comparisons show good agreement between our instrument and the DOEs instrument on both days at all supersaturations. After the last inter-comparison on 22/09/01, there was no indication of a malfunction by our instrument through the rest of the campaign. Unfortunately, the instrument was dropped during demobilization. A post-campaign calibration was conducted, which showed a substantial departure from the pre-campaign calibration. The drop may have damaged the instrument’s ability to produce the desired supersaturations. Therefore, we do consider the data after 22/09/01 to be correct, but it should be used with caution. The CCN counter sampled for 3 minutes at each supersaturation setpoint (0.2, 0.4, 0.6, 0.8, 1.0, and 1.2%). At the end of a cycle, the instrument was set to 0.01% supersaturation for 5 minutes. The data is comprised of the last 60 seconds of each supersaturation set point (0.2, 0.4, 0.6, 0.8, 1.0, and 1.2%) to ensure the instrument stabilized and was able to reach thermal equilibrium. We removed the data during the periods where there were operational difficulties, setup, or maintenance. This data was collected for ARM Field Campaign AFC07055 and supported by DOE ASR grant DE-SC0021047. For any further questions, please feel free to contact the instrument PI, Sarah D. Brooks, sbrooks@tamu.edu. Rose et. al. Calibration and Measurement Uncertainties of a Continuous-Flow Cloud Condensation Nuclei Counter (DMT-CCNC): CCN Activation of Ammonium Sulfate and Sodium Chloride Aerosol Particles in Theory and Experiment. Atmos. Chem. Phys., 8, 1153-1179, 2008. Deng et. al. Using Raman Microspectroscopy to Determine Chemical Composition and Mixing State of Airborne Marine Aerosols over the Pacific Ocean. Aerosol Science and Technology, Vol 48, Issue 2, 2014. Kesten et. al. Calibration of a TSI Model 3025 Ultrafine Condensation Particle Counter. Aerosol Science and Technology, 15:2, 107-111, 1991. Gormley et. al. Diffusion from a Stream Flowing through a Cylindrical Tube. Proceedings of the Royal Irish Academy, Vol 52, 163-169, 1948. Aerosol Measurement: Principles, Techniques, and Applications, Third Edition. John Wiley & Sons, Inc, 2011. Crane et. al. Inertial Deposition of Particles in a Bent Pipe. Journal of Aerosol Science, Vol 8, 161-170, 1977. Dirgo et. al. Cyclone Collection Efficiency: Comparison of Experimental Results with Theoretical Predictions. Aerosol Science and Technology, 4:4, 401-415, 1985.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Theoretical study of various nonlinear phenomena in plasma systems and scaling of magneto-inertial-fusion targets

Plasma physics is an exciting field of study with a wide variety of nonlinear processes that come into play. Examples of such processes include the interaction of small-scale turbulence with large-scale plasma structures and the nonlinear saturation of plasma instabilities, for example those of magneto-hydrodynamical nature. During this Truman LDRD project, I studied a collection of nonlinear problems that are of interest to the field of plasma physics. This LDRD report summarizes four main research accomplishments. First, a new statistical model for describing inhomogeneous drift-wave turbulence inter- acting with zonal flows was developed. This new model includes the effects of nonlinear wave-wave collisions, which are expected to change the spectrum of the underlying DW turbulence and therefore the generation of zonal flows. Second, a new mathematical formalism was proposed to systematically apply the non- linear WKB approximation to general field theories, including those often used in fluid dynamics. This formalism represents an interesting tool for studying physical systems that show an explicit scale separation. Third, a weakly nonlinear model was developed to describe the magneto-Rayleigh-Taylor instability. This instability is of paramount importance to understand as it can reduce the performance of magnetic-inertial-fusion (MIF) platforms. The developed models captures the effects of harmonic generation and saturation of the linear growth of the instability. Finally, a framework was proposed for scaling magneto-inertial fusion (MIF) targets to larger pulsed-power drivers. From this framework, a set of scaling rules were derived that conserve the physical regimes of MIF systems when scaling up in peak current. By doing so, deleterious nonlinear processes that affect MIF performance may be kept at bay.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Spectrally Stabilized Interface Capturing Formulation and Implementation in Nek5000/NekRS

This report documents the formulation of a novel level-set method for incompressible two-phase flows in the continuous Galerkin (CG) high order spectral element framework. The overall method hinges on a novel implementation of the spectral vanishing viscosity (SVV) operator for the stabilization of linear/non-linear hyperbolic problems. The multidimensional SVV convolution kernels, which in essence, have a similar effect as a high pass filter applied to the derivatives, are formulated by exploiting the tensor product form, analogous to the construction of the usual stiffness matrix system. The resulting kernels are directionally decoupled and ensure a linear, symmetric positive definite, elliptic matrix operator. The SVV formulation is demonstrated to provide a robust stabilizing mechanism through challenging linear and non-linear hyperbolic problems, including problems pertinent to the level-set formulation. The two-phase framework conceptualized herein is based on the conservative level-set (CLS) method which represents the interface between the fluids by the 0.5 iso-contour of the smoothed Heaviside function. The CLS method is augmented with a preconditioning procedure for interface normals using the signed distance function which precludes the manifestation of spurious oscillations in the vicinty of the interface. Further, the existing mixed explicit-implicit approach for the solution of Navier-Stokes equations in Nek5000, as described in Tomboulides et al, is augmented with a pressure coefficient splitting approach for the Poisson equation, which greatly accelerated the convergence of pressure solver for two-phase systems with large density ratio. The robustness and accuracy of the overall two-phase method is demonstrated through canonical challenging problems involving high density and viscosity ratios, with and without surface tension. The two-phase formulation is wholly implemented in Nek5000 and the SVV stabilization method is implemented in NekRS, which is the essential precursor to the two-phase framework, undergoing active development.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Dynamic modelling and control strategy of a temperature-driven metal hydride cooling system for buildings

A temperature-driven coupled metal hydride (MH) based thermal energy storage (TES) system can allow to shave and shift the peak energy demand in buildings. The high energy density and long-term (seasonal) energy storage capability are its major advantages over other energy storage methods. The dynamic nature of the MH operation, however, requires controlled hydrogen transfer between the coupled MHs at a rate needed to meet the building's transient load. While temperature-driven MH systems are studied in the literature, their application in buildings and control are scarcely reported. Here, this paper presents a control-based dynamic modeling of the temperature-driven coupled MH-TES system for building cooling applications. The dynamic model is developed in MATLAB(R) Simulink environment, considering the thermodynamic and kinetic behaviors of the MH systems. Based on a preliminary analysis of a property database of over 337 hydrides, we select around 1600 MH pairs suitable for building cooling applications. Each of these MH pairs is studied for their performance using the dynamic model, and among all, Zr 0.76 Ti 0.24 Ni 1.16 Mn 0.63 V 0.14 Fe 0.18 -Ti 0.85 Zr 0.15 Cr 1.2 Mn 0.8 MH pair showed fast dynamics along with high coefficient of performance (COP) of 0.71. A parametric investigation is performed on this MH pair to understand the effect of operating temperatures. Finally, three proportional-integral (PI) feedback controllers are investigated to regulate the temperature, pressure and mass exchange between the coupled MH pairs. The developed PI controller is sufficiently capable of rejecting the signal noise from the hydrogen flow and internal heat exchange processes with root mean square error of 5.78 W between reference and actual cooling load.

08 HYDROGEN↗
DOE OSTI2026

A cell-centered AMR-ALE framework for 3D multi-material hydrodynamics. Part I: Lagrangian and indirect Euler AMR algorithms

Many applications of physics and engineering involve wide ranges of time and spatial scales. The numerical simulation of localized small scales such as shock waves and material interfaces requires a large number of computational cells in these regions. For these applications, Lagrangian and Arbitrary-Lagrangian-Eulerian (ALE) related methods are engaging since the moving mesh feature naturally brings mesh cells on shock discontinuities and material interfaces are carefully captured. In addition, Adaptive-Mesh-Refinement (AMR) strategies aim to optimize computational resources by concentrating finer mesh cells only in areas of interest while using coarser cells elsewhere. A key but challenging AMR requirement consists in efficiently distributing the computational effort to achieve high accuracy without the prohibitive computational costs associated with uniformly fine grids. Here, in this document, the coupling of the p4est AMR library with a cell-centered Lagrangian scheme is presented with the goal to perform reliable 3D Lagrangian-AMR and indirect Euler-AMR multi-material simulations. In particular, it is shown that starting from a 3D indirect ALE code, the memory management and load balancing requirements can be delegated to an external library (here the p4est library) to unlock ALE-AMR capabilities. First, we present a strategy to transcribe the octant-based connectivity of the 3D AMR framework with that of an unstructured mesh of polygonal cells used in Lagrangian hydrodynamics. Then, we show how refinement and coarsening operations must be adapted to the particular Lagrangian framework to ensure the conservation of volume during those steps. Finally, several numerical test cases are presented that demonstrate the capabilities of the Lagrangian-AMR and indirect Euler-AMR algorithms.

3D cell-centered Lagrangian numerical scheme↗
DOE OSTI2024

Trust Model Utilization for Energy Grid Communication

The internet information that is used by the Energy Grid of Things requires both preventative security measures as well as surveillance measures. The preventative security measures include certificates, encryption, and all of the basic security protocols as defined by published standards. The surveillance measures include monitoring information flow activities and evaluating these messages for indications of potential security attacks. We describe in this paper the utilization of a Distributed Trust Model that was developed specifically for monitoring communication within an Energy Grid of Things. The goal for the Distributed Trust Models is to provide a level of aggregate trust that a Distributed Energy Resource Management System can meet its grid service obligations, as opposed to a detailed individual Distributed Energy Resources assessment.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2023

Physical Properties of Moist, Fermented Corn Kernels

A novel approach to producing corn stover biomass feedstock has been investigated. In this approach, corn grain and stover are co-harvested at moisture contents much less than typical corn silage. The grain and stover are conserved together by anaerobic storage and fermentation and then separated before end use. When separated from the stover, the moist, fermented grain had physical characteristics that differ from typical low-moisture, unfermented grain. A comprehensive study was conducted to quantify the physical properties of this moist, fermented grain. Six corn kernel treatments, either fermented or unfermented, having different moisture contents, were used. Moist, fermented kernels (26 and 36% w.b. moisture content) increased in size during storage. The fermented kernels’ widths and thicknesses were 10% and 15% greater, respectively, and their volume was 28% greater than the dry kernels (15% w.b.). Dry basis particle density was 9% less for moist, fermented kernels. Additionally, the dry basis bulk density was 29% less, and the dry basis hopper-discharged mass flow rate was 36% less. Moist, fermented grain had significantly greater kernel-to-kernel coefficients of friction and angles of repose compared to relatively dry grain. The friction coefficient on four different surfaces was also significantly greater for fermented kernels. Fermented corn kernels had lower individual kernel rupture strengths than unfermented kernels. These physical differences must be considered when designing material handling and processing systems for moist, fermented corn grain.

09 BIOMASS FUELS↗