Search NASA⌕ Search

SEARCH · Search NASA

Results for “Reduced Order Modeling”

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

FORCE-DISPATCHES Integration - Initial Demonstration

Integrated energy systems (IES) combine, in mutually beneficial ways, power from variable renewable energy sources and nuclear power plants (NPP) to improve economic viability under uncertain market and weather conditions. The open-source Framework for Optimization of Resources and Economics (FORCE) tool suite, developed at Idaho National Laboratory (INL), has enabled comprehensive modeling and simulation of IES. The capabilities within FORCE include grid portfolio optimization through the Holistic Energy Resource Optimization Network (HERON) and the transient process model analysis library HYBRID, among others. Continuous efforts and investments from the IES programs have been made to expand and improve the versatility of the FORCE toolset in fiscal year 2022. Code-coupling and cross-tool communication have been important methods for improving this versatility. This report focuses on an additional workflow in the HERON tool for capacity and dispatch stochastic optimization through integration with the external tool Design Integration and Synthesis Platform to Advance Tightly Coupled Hybrid Energy Systems (DISPATCHES). DISPATCHES was primarily developed by the National Energy Technology Laboratory, in collaboration with other national laboratories, which included INL, universities, and industry partners. It is coupled to a library of algebraic models for specific plant components, and to a framework for stochastic optimization different from that provided in the current Risk Analysis Virtual Environment (RAVEN)-running-RAVEN algorithm in HERON. HERON currently conducts stochastic optimization via an outer-inner loop: it optimizes over variable capacity on the outer loop, and at each step within the capacity parameter space, conducts an inner optimization over scenarios (of market signals, demand, and/or weather patterns) and hourly dispatch throughout a user-specified number of years. On the other hand, DISPATCHES conducts stochastic optimization via an “all-at-once” strategy in which capacity variables are optimized at the same level as dispatch variables, as all scenarios are considered at once. The latter method works especially well for projects of limited size and project length, as the necessary computational power and memory increases with the number of variables and scenarios. The new capability to use the DISPATCHES workflow in HERON enhances standalone simulations by leveraging FORCE tools—namely, the economic metrics from the Tool for Economic Analysis (TEAL) and reduced-order model (ROM) sampling from RAVEN. The initial demonstration of the DISPATCHES workflow simulates an existing nuclear-case flowsheet within the DISPATCHES repository—this models a NPP with a secondary revenue stream for hydrogen production. Electrical output from the plant is converted to hydrogen via a proton-exchange membrane (PEM) electrolyzer, hydrogen tanks are used for storage, and an additional turbine is added for hydrogen combustion. Continued work regarding this FORCE-DISPATCHES integration will include automatic generation of DISPATCHES models from HERON inputs, offering analysts the option of using either the RAVEN-runsRAVEN or DISPATCHES workflow to solve technoeconomic optimization problems.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2025

Influence of simple terrain on the spatial variability of a low-level jet and wind farm performance in the AWAKEN field campaign

In wind energy research, scientific challenges are often associated with complex terrain sites, where orography, vegetation, and buildings disrupt flow uniformity. However, even sites characterized as simple terrain can exhibit significant spatial variability in wind speed, particularly during stable boundary layers (SBLs) and low-level jets (LLJs). This study investigates these terrain interactions using both simulations and observations from the American WAKe ExperimeNt (AWAKEN). We employ a multiscale Weather Research and Forecasting (WRF) model simulation, integrating mesoscale forcing in the coarse domains and representing three rows of turbines from the King Plains wind farm as generalized actuator disks (GAD) in the large-eddy simulation (LES) domains. During a nocturnal LLJ event on 3 April 2023, the downstream, wake-affected turbine rows outperformed the upstream, unwaked row by 25 %–51 %. This counterintuitive result arises from terrain-induced streamwise variations in hub-height wind speed of approximately 4 m s −1 over 5 km – equivalent to ∼50 % of the upstream reference speed. This enhancement outweighs the wake-induced reduction in mean wind speed (∼12 %) and global blockage effects reported in the literature (∼1 %–3.4 %). The multiscale simulations capture the intra-farm spatial variability in power performance observed in SCADA data. Terrain-induced vertical displacement of the LLJ, coupled with large wind shear below the jet maximum, drives the substantial streamwise acceleration within the wind farm. These findings underscore the importance of accounting for spatial variability related to terrain, even in simple landscapes, particularly during LLJ conditions. Incorporating such effects into reduced-order modeling frameworks for wind farm design and control could significantly enhance their effectiveness.

17 WIND ENERGY↗
NASA NTRS2020

Multiple Shaker Placement for Ground Vibration Test of X-59 Aircraft Using Topology Optimization

A multiple shaker placement methodology is developed and tested using a topology optimization technique. Current multiple shaker placement methodology requires optimum accelerometer placement and optimum single-shaker placement techniques. The proposed methodology is tested using a finite element model of the X-59 Low Boom Flight Demonstrator aircraft. The effective independence and the driving point acceleration transfer function (DPATF) methods are used for the accelerometer placement study. In this study, four shakers are used to excite each mode more effectively during the ground vibration test; all the modes of interest thus are separated into four groups. Each shaker takes care of a separate group of modes. Grouping the modes of interest is performed utilizing topology optimization. The number of modes for each group therefore will be automatically decided during grouping. For each group of modes, perform the following two steps to determine optimal location of four shakers: 1) At each accelerometer location, compare the magnitude of DPATF values at natural frequencies, select the minimum value, and make a vector with these minimum values of the DPATF magnitudes for each group; and 2) Select the degrees of freedom corresponding to the maximum value of this vector. The objective function value is the maximum value of the vector with minimum value of the magnitude of the superposed acceleration transfer function. This objective function value is maximized by changing the modes for each group. Forty accelerometers are enough to have good correlation between mode shapes obtained from the reduced order model and the simulated ground vibration test.

Optimization Analysis↗
DOE OSTI2016

Galerkin v. least-squares Petrov–Galerkin projection in nonlinear model reduction

Least-squares Petrov–Galerkin (LSPG) model-reduction techniques such as the Gauss–Newton with Approximated Tensors (GNAT) method have shown promise, as they have generated stable, accurate solutions for large-scale turbulent, compressible flow problems where standard Galerkin techniques have failed. Furthermore, there has been limited comparative analysis of the two approaches. This is due in part to difficulties arising from the fact that Galerkin techniques perform optimal projection associated with residual minimization at the time-continuous level, while LSPG techniques do so at the time-discrete level. This work provides a detailed theoretical and computational comparison of the two techniques for two common classes of time integrators: linear multistep schemes and Runge–Kutta schemes. We present a number of new findings, including conditions under which the LSPG ROM has a time-continuous representation, conditions under which the two techniques are equivalent, and time-discrete error bounds for the two approaches. Perhaps most surprisingly, we demonstrate both theoretically and computationally that decreasing the time step does not necessarily decrease the error for the LSPG ROM; instead, the time step should be ‘matched’ to the spectral content of the reduced basis. In numerical experiments carried out on a turbulent compressible-flow problem with over one million unknowns, we show that increasing the time step to an intermediate value decreases both the error and the simulation time of the LSPG reduced-order model by an order of magnitude.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

Parameter identification methods for low-order gray box building energy models: A critical review

The body of knowledge on gray box building energy modeling (GBBEM) has been developed over the past few decades and has undergone some important changes recently. Starting with simple methods and simple buildings, the science of GBBEM has grown to encompass complex and more computationally intensive techniques and complex commercial buildings. Numerous works including a recent review have considered model structure and inputs in a nearly systematic way, but no extant work systematically reviews the approaches for GBBEM parameter search initialization and final identification, despite this being arguably the most difficult and impactful part of the modeling process. To this end, we critically review 55 extant works describing advantages, limitations, and domain of applicability of several classes of parameter initialization and optimization techniques specifically for GBBEM. We categorize the classes of methods and analyze the evidence of their applicability for different applications within the field of GBBEM. We find an emerging consensus that initialization of parameter searches for anything other than the simplest building elements is often challenging and sometimes requires a stochastic approach to begin the parameter identification process. After this initial process, faster methods have been used in some cases but often the nature of the problem requires stochastic methods for this portion of the process as well. For less complex systems, more deterministic and efficient methods have been shown to be effective. Finally, we draw conclusions as to the domain of applicability of different classes of initialization and optimization techniques for GBBEEM and offer suggestions for research directions that are likely to prove fruitful.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2017

Process-Structure Linkages Using a Data Science Approach: Application to Simulated Additive Manufacturing Data

A novel data science workflow is developed and demonstrated to extract process-structure linkages (i.e., reduced-order model) for microstructure evolution problems when the final microstructure depends on (simulation or experimental) processing parameters. Our workflow consists of four main steps: data pre-processing, microstructure quantification, dimensionality reduction, and extraction/validation of process-structure linkages. These methods that can be employed within each step vary based on the type and amount of available data. In this paper, this data-driven workflow is applied to a set of synthetic additive manufacturing microstructures obtained using the Potts-kinetic Monte Carlo (kMC) approach. Additive manufacturing techniques inherently produce complex microstructures that can vary significantly with processing conditions. Using the developed workflow, a low-dimensional data-driven model was established to correlate process parameters with the predicted final microstructure. In addition, the modular workflows developed and presented in this work facilitate easy dissemination and curation by the broader community.

36 MATERIALS SCIENCE↗
DOE OSTI2021

AlphaBuilding ResCommunity: A multi-agent virtual testbed for community-level load coordination

Training and validating algorithms in a simulation testbed can accelerate research and applications of optimal control of residential loads to improve energy flexibility and grid resilience. We developed an open-source simulation environment, AlphaBuilding ResCommunity, that can be used to train and validate algorithms to control a single thermostatically controlled load (TCL) or coordinate a group of TCLs. We used reduced-order models to simulate the thermodynamics of TCLs, and the parameter values were determined from the connected smart thermostat data of real households. The environment was built upon the standardized OpenAI Gym interface. Ancillary functions, such as retrieving the parameters and weather forecasts, are provided to facilitate control strategies that require predictive information. Compared with existing efforts, AlphaBuilding ResCommunity has three advantages: (1) more realistic model settings because the parameter values are identified from actual household operating data, and modelling and measurement uncertainty are considered; (2) passive thermal storage control; and (3) ease of use due to a simple software dependency and standardized interface. We demonstrated the applications of the environment by implementing a Kalman Filter and Model Predictive Control on a single TCL and a Priority-Stack-Based Control and Alternating Direction Method of Multipliers to coordinate multiple TCLs for load tracking.

Wang, Z↗
DOE OSTI2017

Cost comparison of printed circuit heat exchanger to low cost periodic flow regenerator for use as recuperator in a s-CO 2 Brayton cycle

Supercritical Carbon Dioxide (sCO 2 ) power cycles have the potential to deliver high efficiency at low cost. However, in order for an sCO 2 cycle to reach high efficiency, highly effective recuperators are needed. These recuperative heat exchangers must transfer heat at a rate that is substantially larger than the heat transfer to the cycle itself and can therefore represent a significant portion of the power block costs. Regenerators are proposed as a cost saving alternative to high cost printed circuit recuperators for this application. A regenerator is an indirect heat exchanger which periodically stores and releases heat to the working fluid. The simple design of a regenerator can be made more inexpensively compared to current options. The objective of this paper is a detailed evaluation of regenerators as a competing technology for recuperators within an sCO 2 Brayton cycle. The level of the analysis presented here is sufficient to identify issues with the regenerator system in order to direct future work and also to clarify the potential advantage of pursuing this technology. A reduced order model of a regenerator is implemented into a cycle model of an sCO 2 Brayton cycle. An economic analysis investigates the cost savings that is possible by switching from recuperative heat exchangers to switched-bed regenerators. The cost of the regenerators was estimated using the amount of material required if the pressure vessel is sized using ASME Boiler Pressure Vessel Code (BPVC) requirements. The cost of the associated valves is found to be substantial for the regenerator system and is estimated in collaboration with an industrial valve supplier. The result of this analysis suggests that a 21.2% reduction in the contribution to the Levelized Cost of Electricity (LCoE) from the power block can be realized by switching to a regenerator-based system.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2024

Demonstration of a multi-channel fluidized bed particle–supercritical carbon dioxide heat exchanger for concentrating solar applications

High-temperature thermal energy storage in oxide particles at temperatures above 600°C can couple concentrated solar energy with high-efficiency thermal power cycles to provide dispatchable solar-driven electricity. Challenges remain in developing cost-effective primary heat exchangers, which require expensive alloys, to extract the high-temperature thermal energy from the particles to power cycle fluids, such as supercritical CO 2 (sCO 2 ) in recuperated Brayton cycles. To explore one pathway for cost-effective, high-temperature particle heat exchangers, the current study demonstrates a shell-and-plate, particle–sCO 2 heat exchanger with narrow- channel fluidized beds coupled with micro-channel sCO 2 flows in the heat exchanger walls. This study evaluates the feasibility of multiple parallel, narrow-channel fluidized beds in shell-and-plate particle–sCO 2 HXs, to achieve high bed-wall heat fluxes at elevated temperatures. A reduced-order model simulates the narrow- channel, fluidized-bed particle–sCO 2 heat exchanger to design the fluidized bed geometry, in terms of depth, height, and number of channels,for a nominal 40-kWth heat exchanger at particle and sCO 2 inlet temperatures up to 600 °C and 400 °C respectively. The resulting shell-and-plate heat exchanger design operates with bubbling fluidization of the downward-flowing oxide particles to enhance bed-wall heat transfer. The heat exchanger core is fabricated with etched sCO 2 micro-channels in thin wall plates that are diffusion bonded to spacer frames to form the shell-and-plate structure with 12 parallel, fluidized bed channels, 10.4 mm deep. The heat exchanger is tested at the National Solar Thermal Test Facility at Sandia National Laboratories with CARBOBEAD HSP particles at design particle flow rates of 0.20 kg s –1 and inlet temperatures up to 525 °C. Results show that fluidization across multiple parallel channel beds can maintain uniform particle inventory with a common freeboard zone above the heat exchanger core. Bubbling fluidization improves particle–wall heat transfer coefficients but also increases axial dispersion of particle thermal energy, which lowers the log- mean temperature difference such that total heat transfer remains relatively constant to within ±10% over a broad range of fluidization gas velocities. The axial dispersion required particle and sCO 2 flow rates to be increased by 25% over model-designed conditions to achieve the targeted 40 kWth, which indicates the importance of incorporating axial dispersion into heat exchanger design models and of deploying bed structures to suppress it. Furthermore, this study demonstrates the feasibility and preferred fluidizing gas conditions for particle heat exchangers for releasing high-temperature thermal energy storage systems.

14 SOLAR ENERGY↗
DOE OSTI2023

Continuous recovery of phosphoric acid and Rare-Earths containing particles from phosphoric acid sludge using a decanter centrifuge

Recovery of rare earth elements (REEs) from various industrial and natural streams currently draws significant attention in efforts to meet the demands of the manufacturing industry. Among many industrial byproducts and waste streams, phosphoric acid sludge could be one of the most economically feasible resources for the recovery of REEs because solid particles in the sludge contain relatively concentrated REEs, up to 3,000 ppm, while the liquid component of the sludge is valuable phosphoric acid (P 2 O 5 ) that can be recovered and returned to the main product. Due to high viscosity and large solids content (e.g., 30–40 %), however, this byproduct stream requires multistep separation and purification processes. In this study, a single-step process involving a continuous-flow decanter centrifuge (CFDC) was employed to investigate its feasibility for continuous solid/liquid separation from real phosphoric acid sludge. High centrifugal forces generated from up to 1500 G gravity acceleration separate solid particles from the sludge, generating a liquid-rich stream and a solids-rich stream at the exit of the CFDC. A single pass of phosphoric-acid sludge through the CFDC yielded 95 % liquid recovery and 90 % recovery of REEs-containing solids from 20 to 34 wt% solids-containing sludge. Finally, a reduced order model developed for the CFDC operation showed good agreement with experimental data, and preliminary technoeconomic analysis revealed potential process feasibility.

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

Probabilistic Risk Assessment of CO2 Trapping Mechanisms in a Sandstone CO2-EOR Field in Northern Texas, USA

This study employs a polynomial chaos expansion (PCE) technique to develop the reduced order models (ROMs) integrated with the Monte Carlo simulations to quantify uncertainties of CO 2 storage associated with enhanced oil recovery (EOR) in a commercial-scale CO 2 -EOR field. The Morrow B sandstone at the Farnsworth Unit (FWU) EOR field in northern Texas was selected as a case study. A 3-D reservoir model, including simulation of a five-spot injection/production well pattern, was constructed based on the data from well logs and cores, and an updated seismic interpretation in the FWU. Porosity and permeability were considered as sources of uncertainty. The 1000-realization heterogeneous random fields of porosity and permeability were generated based on the variogram analysis on the core measurements of the FWU. Three model outputs of interest include the amount of CO 2 trapped by three trapping mechanisms, i.e. hydrodynamic trapping, oil dissolution trapping, and aqueous dissolution trapping. The total of 25 runs were conducted with the simulations of a 20-year CO 2 -EOR period, and a 20-year monitoring period to develop the ROMs between uncertain model input parameters (i.e., porosity and permeability in this study) and model outputs at the end of simulation for each grid cell using PCE approach. The high coefficient of determination (R 2 ) and small normalized root mean square error (NRMSE) between the simulated (using reservoir modeling) and predicted (using ROMs) CO 2 storage indicate that the ROMs are acceptable and reliable to be used for the predictions. Given the 1000 Monte Carlo samplings of the model input parameters, cumulative distribution functions (CDFs) and uncertainty bounds (10 th and 90 th percentiles) of model outputs were estimated based on ROMs. At the end of simulation, the most injected CO 2 at the FWU was stored by hydrodynamic trapping (between 121,400 and 166,860 tonnes), followed by oil dissolution trapping (between 55,303 and 76,816 tonnes), and aqueous dissolution trapping (between 1979 and 2751 tonnes).

58 GEOSCIENCES↗
DOE OSTI2019

Development of a topology analysis tool for fifth-generation district heating and cooling networks

Fifth generation district heating and cooling networks are characterized by supply temperatures in the ambient range of 15-25 degrees C, which not only reduces heat loss but also allows for integrating various kinds of low-temperature waste heat sources. The ability of these networks to absorb waste heat that is normally unrecoverable makes them an attractive solution for the future energy supply of urban areas. To enhance the adoption of these networks, this paper describes the development of a software tool to analyze the feasibility of fifth-generation of district heating and cooling systems in both new and existing districts. The research attempts to answer the question: 'Which buildings should be connected to a low-temperature hydraulic network given both the incremental benefits and costs relative to the scenario of decentralized (dedicated) heating and cooling systems for each building?' Therefore, all possible network layouts of the buildings are considered where the heating and cooling demands of each building are met by the fifth-generation of district heating and cooling network. For the heating and cooling load characterization of the buildings, reduced-order models are used while renewable energy and waste heat sources are included in the network. The focus of the research lies in the development of a hydraulic model for flexible use in the urban energy modeling and the optimization of the fifth-generation of district heating and cooling network topology for a given urban district. This research displays simulation results that quantify the performance of the fifth-generation of district heating and cooling network based on various output metrics, including primary energy usage, carbon dioxide emissions, and network implementation cost.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2022

Demonstration of RELAP5-3D for transient analysis of a dual coolant lead lithium fusion blanket concept

Our research aims to build a multiphysics framework for transient analysis of the Dual- Coolant Lead-Lithium (DCLL) blanket design of the proposed Fusion Nuclear Science Facility (FNSF) for all Onticipated Operational Occurrences (AOOs) using Reduced Order Models (ROMs). This paper presents an assessment of using RELAP5-3D for transient thermal-hydraulic blanket analysis of the DCLL design. Preliminary validation studies for the RELAP5-3D properties were conducted using a representative vertical flow loop model and a comparison study using experimental data from the thermal convection corrosion loops at ORNL. These studies inspected the thermal hydraulic response of Lead-Lithium Eutectic (PbLi) systems within the RELAP5-3D code; mainly pressure drop and heat transfer, as compared to analytical calculations and experimental data. The Magnetohydrodynamic (MHD) pressure drop effect for liquid metal under the effects of a uniform magnetic field was implemented into our models using an equivalent forms loss coefficient based on the phasic momentum equations. We verified our implementation of MHD pressure drop within RELAP5-3D using a systematic approach consisting of two studies. Models were developed to perform comparison studies based on experimental data from the Argonne Liquid metal EXperiment (ALEX) facility and the Magnetohydrodynamic PbLi Experiment (MaPLE) facility. Using the developed validation basis, we built a simplified DCLL channel model within RELAP5-3D, inspected MHD pressure drop and performed transient temperature analysis. The DCLL channel model used radial heating data from the International Thermonuclear Experimental Reactor (ITER) which was normalized to the FNSF design. Here, a representative startup transient was performed using a relative power curve from the initial ramp up of a long-term power pulse of the Tore Supra reactor. This study demonstrated that RELAP5-3D ROMs are capable of accurate transient analysis of the DCLL blanket. We have since utilized MCNP to develop heating profiles based on the FNSF design and are working toward implementing them into our DCLL model.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Uncertainty analysis for techno-economic and life-cycle assessment of wet waste hydrothermal liquefaction with centralized upgrading to produce fuel blendstocks

Wet waste hydrothermal liquefaction is a promising technology for producing transportation fuels with much lower greenhouse gases emissions than petroleum-based fuels. However, its techno-economic and life cycle assessment are primarily based on laboratory scale testing data, subject to considerable uncertainties, and even bias, due to knowledge gaps. Here, a preliminary uncertainty analysis of key economic measures was conducted based on the 2019 state-of-technology model for biocrude production. Building on the preliminary analysis, this work presents a comprehensive uncertainty analysis in both economic and environmental measures of the entire supply chain of wet waste hydrothermal liquefaction to fuel blendstocks including biocrude upgrading based on the 2021 state-of-technology model. The analysis includes the most recent developments in hydrothermal liquefaction and biocrude upgrading technologies and Monte Carlo simulation based on an integrated model system including an improved reactor yield model, reduced-order process model, discounted cash flow economic model and simplified life-cycle assessment model. The estimated biocrude yield ranges from 42.2% to 52.4% with a median of 47.3%. The estimated fuel yield ranges from 34.7% to 42.7% with a median of 38.7%. The estimated minimum fuel selling price ranges from $\$ $2.28/gge to $\$ $3.45/gge with a median of $\$ $2.80/gge. Relative to petroleum-derived diesel, the estimated reduction in supply chain greenhouse gas emissions ranges from 73.4% to 81.8% with a median of 77.7%. Compared to the 2019 state-of-technology analysis, a significant improvement in biocrude selectivity and economic measures and reduction in uncertainties were achieved due to the incorporation of additional continuous experimental data sets, technology development and de-risking, and improvement in model accuracy.

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

Nuclear safety Enhanced: A Deep dive into current and future RAVEN applications

As the horizon of nuclear energy expands with the advent of small modular reactors, IV generation reactors, and fusion reactors, there is a growing perspective that the licensing process could benefit from a more comprehensive approach. Moving beyond traditional deterministic and PRA analysis might pave the way for a novel safety analysis paradigm propelled by the increasing computational power at our disposal. This paper explores different methodologies that can improve the outcomes of nuclear safety analysis. These range from uncertainty quantification techniques, aimed at enhancing the precision of safety margins, to deploying dynamic event trees by driving system code simulations, capturing the potential evolutions of severe accidents. These methodologies introduce innovative dimensions to safety analysis, considering the consequences of postulated events and the dynamics of accident sequences. However, they also bring forth challenges, especially in managing the complexity and sheer volume of potential scenarios. The paper touches upon some strategies to counter these challenges, emphasizing the importance of adaptability and continuous evolution in the face of emerging nuclear safety concerns. Additionally, the paper sheds light on the need for advanced tools to apply these methodologies. Among these tools is RAVEN, an open-source software designed for parametric and probabilistic analyses. Its core components, including distribution, sampler, and reduced order model, enable various applications, from risk assessment and mitigation to dynamic learning and plant control logic simulations.

97 - MATHEMATICS AND COMPUTING↗
DOE OSTI2017

Targeted energy transfer in laminar vortex-induced vibration of a sprung cylinder with a nonlinear dissipative rotator

We computationally investigate the dynamics of a linearly-sprung circular cylinder immersed in an incompressible flow and undergoing transverse vortex-induced vibration (VIV), to which is attached a rotational nonlinear energy sink (NES) consisting of a mass that freely rotates at constant radius about the cylinder axis, and whose motion is restrained by a rotational linear viscous damper. The inertial coupling between the rotational motion of the attached mass and the rectilinear motion of the cylinder is “essentially nonlinear”, which, in conjunction with dissipation, allows for one-way, nearly irreversible targeted energy transfer (TET) from the oscillating cylinder to the nonlinear dissipative attachment. At the intermediate Reynolds number , the NES-equipped sprung cylinder undergoes repetitive cycles of slowly decaying oscillations punctuated by intervals of chaotic instabilities. During the slowly decaying portion of each cycle, the dynamics of the cylinder is regular and, for large enough values of the ratio of the NES mass to the total mass (i.e., NES mass plus cylinder mass), can lead to significant vortex street elongation with partial stabilization of the wake. As approaches zero, no such vortex elongation is observed and the wake patterns appear similar to that for a sprung cylinder with no NES. Here we apply proper orthogonal decomposition (POD) to the velocity flow field during a slowly decaying portion of the solution and show that, in situations where vortex elongation occurs, the NES, though not in direct contact with the surrounding fluid, has a drastic effect on the underlying flow structures, imparting significant and continuous passive redistribution of energy among POD modes. We construct a POD-based reduced-order model for the lift coefficient to characterize energy transactions between the fluid and the cylinder throughout the slowly decaying cycle. We introduce a quantitative signed measure of the work done by the fluid on the cylinder over one quasi-period of the slowly decaying response and find that vortex elongation is associated with a sign change of that measure, indicating that a reversal of the direction of energy transfer, with the cylinder “leaking energy back” to the flow, is responsible for partial stabilization and elongation of the wake. We interpret these findings in terms of the spatial structure and energy distribution of the POD modes, and relate them to the mechanism of transient resonance capture into a slow invariant manifold of the fluid–structure interaction dynamics.

42 ENGINEERING↗
DOE OSTI2025

Modal dynamics of wind turbine wake meandering from lidar observations

Horizontal scans from nacelle-mounted lidars provide time series measurements of wind turbine wakes across diverse atmospheric conditions, enabling analysis of coherent turbulent structures that influence wake meandering through proper orthogonal decomposition (POD). While low-order modes capture the most energetic turbulent structures, our analysis reveals that they do not necessarily dominate wake meandering dynamics. We evaluate more than 16,000 combinatorial reconstructions of the flow field for each inflow case, demonstrating that mode relevance depends on mode symmetry, turbulent kinetic energy content, and inflow characteristics. Cases with low turbulence intensity and large integral timescales show stronger correlations between POD modes and wake meandering, whereas higher turbulence conditions (turbulence intensity > 7%) are less effectively described by reduced-order models. However, the qualitative similarity of POD modes across varied atmospheric conditions suggests the potential existence of a semi-universal basis for representing wind turbine wakes, with implications for improving engineering wake models.

17 WIND ENERGY↗
DOE OSTI2018

Self similarities in desalination dynamics and performance using capacitive deionization

Charge transfer and mass transport are two underlying mechanisms which are coupled in desalination dynamics using capacitive deionization (CDI). We developed simple reduced-order models based on a mixed reactor volume principle which capture the coupled dynamics of CDI operation using closed-form semi-analytical and analytical solutions. We use the models to identify and explore self-similarities in the dynamics among flow rate, current, and voltage for CDI cell operation including both charging and discharging cycles. The similarity approach identifies the specific combination of cell (e.g. capacitance, resistance) and operational parameters (e.g. flow rate, current) which determine a unique effluent dynamic response. We here demonstrate self-similarity using a conventional flow between CDI (fbCDI) architecture, and we hypothesize that our similarity approach has potential application to a wide range of designs. We performed an experimental study of these dynamics and used well-controlled experiments of CDI cell operation to validate and explore limits of the model. For experiments, we used a CDI cell with five electrode pairs and a standard flow between (electrodes) architecture. Guided by the model, we performed a series of experiments that demonstrate natural response of the CDI system. We also identify cell parameters and operation conditions which lead to self-similar dynamics under a constant current forcing function and perform a series of experiments by varying flowrate, currents, and voltage thresholds to demonstrate self-similarity. Based on this study, we hypothesize that the average differential electric double layer (EDL) efficiency (a measure of ion adsorption rate to EDL charging rate) is mainly dependent on user-defined voltage thresholds, whereas flow efficiency (measure of how well desalinated water is recovered from inside the cell) depends on cell volumes flowed during charging, which is determined by flowrate, current and voltage thresholds. Results of experiments strongly support this hypothesis. Results show that cycle efficiency and salt removal for a given flowrate and current are maximum when average EDL and flow efficiencies are approximately equal. We further explored a range of CC operations with varying flowrates, currents, and voltage thresholds using our similarity variables to highlight trade-offs among salt removal, energy, and throughput performance.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗