SEARCH · Engineering Papers

Results for “numerical stability”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 721 records · Page 40

All sources NASA NTRS DOE OSTI

DOE OSTI2021

A phase-field method for boiling heat transfer

Here we present a phase field method for heat transfer in two-phase flow with boiling. The vapor/liquid interface evolution is modeled by the Cahn-Hilliard equation. The phase change rate is determined by accounting for the heat conduction balance on either vapor or liquid side of the interfacial area, depending on which side the temperature is assumed to be maintained at the saturation temperature during boiling. The velocity correction scheme proposed by Dong & Shen [27] is extended to solve the Navier-Stokes equations for a non-solenoidal velocity field, and the entropy viscosity method is employed for stabilization. The phase change model is verified by two-dimensional simulations of a vapor bubble growing in super-heated liquid and in film boiling. In both cases, mesh independence of the results is systematically performed. Subsequently, the method is applied to predict the growth of three-dimensional vapor bubble in a rectangular microchannel with boiling flow, achieving good agreement with experimental measurements and available simulation results using the level-set method. The numerical experiments demonstrate that the required mesh resolution for the phase field method is comparable with that of volume of fluid (VoF) and level-set methods.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2021

Critical Role of Sc Substitution in Modulating Ferroelectricity in Multiferroic LuFeO 3

Understanding how individual dopants or substitutional atoms interact with host lattices enables us to manipulate, control, and improve the functionality of materials. However, because of the intimate coupling among various degrees of freedom in multiferroics, the atomic-scale influence of individual foreign atoms has remained elusive. Here, we unravel the critical roles of individual Sc substitutional atoms in modulating ferroelectricity at the atomic scale of typical multiferroics, Lu 1–x Sc x FeO 3 , by combining advanced microscopy and theoretical studies. Atomic variations in polar displacement of intriguing topological vortex domains stabilized by Sc substitution are directly correlated with Sc atom-mediated local chemical and electronic fluctuations. The local FeO 5 trimerization magnitude and Lu/Sc–O hybridization strength are found to be significantly reinforced by Sc, clarifying the origin of the strong dependence of improper ferroelectricity on Sc content. We find that this study could pave the way for correlating dopant-regulated atomic-scale local structures with global properties to engineer emergent functionalities of numerous chemically doped functional materials.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2025

An Open Benchmark of One Million High-Fidelity Cislunar Trajectories

Cislunar space spans from geosynchronous altitudes to beyond the Moon and will underpin future exploration, science, and security operations. We describe and release an open dataset of one million numerically propagated cislunar trajectories generated with the open-source Space Situational Awareness Python package (SSAPy). The model includes high-degree Earth/Moon gravity, solar gravity, and Earth/Sun radiation pressure; other planetary gravities are omitted by design for computational efficiency. Initial conditions uniformly sample commonly used osculating-element ranges, and each trajectory is propagated for up to six years under a single, fixed start epoch. The dataset is intended as a reusable benchmark for method development (e.g., space domain awareness, navigation, and machine-learning pipelines), a reference library for statistical studies of orbit families, and a starting point for community-driven extensions (e.g., alternative epochs). We report empirically observed stability trends (e.g., a band near ~5 GEO and persistence of some co-orbital classes including L4/L5 librators) as dataset descriptors rather than new dynamical results. The chief contribution is the scale, fidelity, organization (CSV/HDF5 with full state time series and metadata), and open availability, which together lower the barrier to comparative and data-driven studies in the cislunar regime.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2020

Impact of boundary layer simulation on predicting radioactive pollutant dispersion: a case study for HANARO Research Reactor using the WRF-MMIF-CALPUFF modeling system

Wind plays an important role in cases of unexpected radioactive pollutant dispersion, deciding distribution and concentration of the leaked substance. The accurate prediction of wind has been challenging in numerical weather prediction models, especially near the surface because of the complex interaction between turbulent flow and topographic effect. As such, in this study, we investigated the characteristics of atmospheric dispersion of radioactive material (i.e. 137 Cs) according to the simulated boundary layer around the HANARO research nuclear reactor in Korea using the Weather Research and Forecasting (WRF)-Mesoscale Model Interface (MMIF)-California Puff (CALPUFF) model system. We examined the impacts of orographic drag on wind field, stability calculation methods, and planetary boundary layer parameterizations on the dispersion of radioactive material under a radioactive leaking scenario. We found that inclusion of the orographic drag effect in the WRF model improved the wind prediction most significantly over the complex terrain area, leading the model system to estimate the radioactive concentration near the reactor more conservatively. We also emphasized the importance of the stability calculation method and employing the skillful boundary layer parameterization to ensure more accurate low atmospheric conditions, in order to simulate more feasible spatial distribution of the radioactive dispersion in leaking scenarios.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2026

Slow Wake Recovery and Low Turbulence Behind Wind Farms Parameterized in Mesoscale Simulations

Numerical weather prediction (NWP) and climate models equipped with wind-farm parameterizations (WFPs) can simulate cluster wake effects affecting downstream wind farms in both onshore and offshore environments. This study evaluates wake recovery behind a wind farm represented by the NWP-WFP approach in the Weather Research and Forecasting (WRF) model using either the Fitch et al. (2012) or Ma et al. (2022a, b) WFPs. Results are benchmarked against large-eddy simulations (LES) of an idealized offshore wind farm with aligned and staggered layouts under neutral atmospheric stability. Near-farm wake recovery is underestimated in NWP-WFP simulations due to its representation on a coarse mesoscale grid. This limitation leads to slow wake recovery through two interconnected mechanisms: (i) spatial gradients in the wind velocity field are weaker compared to LES and (ii) turbulence kinetic energy (TKE) remains low not because of excessive dissipation but due to insufficient shear production caused by these weakened gradients. For the scenario considered here, a wind-speed bias develops in the near-farm wake and persists into the far wake. Differences between the NWP-WFP simulations and LES emerge within a short distance downstream of the farm exit, where the mesoscale simulations recover too slowly. This reduced recovery contributes approximately 0.15-0.50 m s-1 to the near-farm wind-speed bias. The bias established in this region is not subsequently compensated for downstream but instead propagates into the far wake, where wind-speed differences of approximately 0.4-0.6 m s-1 remain up to 50 km downstream. Higher-resolution mesoscale simulations partially reduce this bias. Increasing turbine-added TKE or including subgrid wake effects provides additional improvement, but neither fully addresses the underlying cause. The slow wake recovery is not caused by limitations of the WFPs themselves, as it also occurs outside their region of influence, and adding subgrid wake effects does not significantly impact recovery. Rather, the slow wake recovery is a consequence of mesoscale flow representation. This behavior is not limited to regions downstream of the wind farm but is less visible within the farm, where wake recovery occurs simultaneously with turbine-induced momentum extraction. These results highlight the need for improved representations of wake recovery both within and downstream of wind farms. While enhanced subgrid modeling, shear-driven TKE production, and refined WFP formulations may improve intra-farm dynamics, accurately capturing near-farm wake recovery downstream remains challenging, as WFPs do not act in this region.

17 WIND ENERGY↗

DOE OSTI2023

Recent Beam Stability Analysis for the EIC

This document summarizes work done in WBS 6.02.02 during FY23. Additional details can be found in the EIC overview paper presented at IPAC’23 [1], and the references therein. Progress has been made in the ESR and HSR lattice design, dynamic aperture optimization, and the RCS lattice design. The impedance and Collective effects are progressing and include the impedance optimization process of the vacuum systems for RCS, HSR, and ESR; collective effects studies; collective effects and beam-beam interaction; coupled bunch instabilities and the crab cavities; low-level RF feedback system design and beam-ion instability. The reversed phasing RF system has been numerically studied for the ESR to mitigate Robinson instability and demonstrate reliable stable beam operation. Various codes, including C++, SPACE [2], ELEGANT [3], and Mbtrack2 [4] have been employed to benchmark the results. The simulation of HSR bunch splitting with beam loading has been performed at 275 GeV energy. The evaluation of Beam Position Monitors aimed to validate their expected performance, with the primary objective being to verify their accuracy. Calculations of Electron Polarization in the ESR, RCS, and HSR are showing good progress. In the RCS, preliminary studies indicate excellent polarization transmission over intrinsic spin resonances, achieving over 90% transmission with improved performance compared to the previous lattice. Preliminary simulations in the ESR indicate encouraging results in minimizing depolarization and improving equilibrium polarization.

43 PARTICLE ACCELERATORS↗

DOE OSTI2022

Methyl Carbamate-Lithium Salt Deep Eutectic Electrolyte for Lithium-Ion Batteries

Deep eutectic electrolytes (DEEs) represent a burgeoning field in electrolytes for energy storage applications. Compared to the more conventionally studied Li-ion electrolyte systems, these electrolytes offer numerous advantages, such as high ion concentration, costs, lower temperature operation, and safety. In this work, the formation of a novel eutectic electrolyte based on lithium salts and methyl carbamate was reported. These DEEs were formed via mixing lithium salts and low-cost methyl carbamate at ratios varying from 1:2 to 1:5 (mol:mol). Here, the DEE formed from methyl carbamate and lithium hexafluorophosphate (LiPF 6 ) at a 1:5 molar ratio generated a 4V stability window and a 25°C conductivity of up to 3.16E –3 S cm –1 . The DEE formed from methyl carbamate and lithium bis(trifluoromethane)sulfonimide (LiTFSI) at a 1:5 molar ratio generated a 3.2 V stability window and a 25°C conductivity of up to 2.87E –3 S cm –1 . The LiPF 6 DEE also demonstrated a discharge capacity of 128 mAh g –1 after 50 cycles with 93% discharge capacity retention at moderate cycling rates of 0.5C in a lithium titanate (LTO)/lithium iron phosphate (LFP) full cell. These results demonstrate the potential for methyl carbamate to generate eutectics from lithium salts and offer alternative electrolytes for use in low-temperature Li-ion battery applications.

25 ENERGY STORAGE↗

DOE OSTI2024

Porous Organic Cage Membranes for Molecular Gas Separations (Final Technical Report for DE-SC0021357)

This proposal aims at demonstrating the development of a novel family of membranes, composed of porous organic cages (POC) which offer the possibility of displaying high separation performance for challenging molecular gas separations relevant to natural gas purification, and olefin/paraffin separation. The proposed POCs synthesized in membrane form will display the most desirable properties of polymers (facile processability and flexibility) and inorganic materials (hierarchically ordered pores with molecular sieving properties) leading to highly selective and permeable membranes. POCs should display distinctive structural, compositional, adsorption and transport properties than those of conventional porous materials, opening the doors for a new research direction in membrane science, and gas separations. Our preliminary results demonstrate the feasibility of preparing POC crystals with controlled size, and continuous POC membranes with remarkable high permeances, and separation ability for CO 2 /CH 4 , N 2 /CH 4 and C 3 H 6 /C 3 H 8 separations serving as a solid foundation for our proposed work. Fundamentally, this proposal aims at elucidating separation mechanisms of different gas mixtures related to natural gas composition, and olefin/paraffin separation over porous organic cage membranes. The proposed research will result in fundamental understanding of adsorption and transport properties of industrially relevant gas molecules through novel microporous membranes, and may lead to the development of a cost effective membrane technology for natural gas purification, and olefin/paraffin separation surpassing the conventional benchmark technology distillation. Furthermore, we aim at demonstrating selective water transport through POC membranes, which can be positively impactful in numerous industrial applications in which water is present. The ability to fabricate thin, chemically and mechanically stable POC membranes for societal relevant gas separations constitute a new and distinctive direction in membrane science. Our proposed work aims at addressing some of the challenges recognized in the Research Agenda for Transforming Separation Science . Specifically: (a) advancing understanding of complex mixtures on separation performance; (b) exploring thermodynamic and kinetic mechanisms through the elucidation of separation mechanisms, and (c) study potential stability issues of the membranes to be assessed by evaluating the long term membrane stability and performance at various temperatures and pressures. The team is uniquely qualified to execute the proposed work. The PI has solid expertise in the rational molecular engineering design of porous crystalline membranes for molecular gas separations. The PNNL collaborator has extensive experience in the synthesis, characterization, and functional applications of microporous crystals, with particular emphasis on gas adsorption.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2024

Effect of Pressure on Structure and Extinction of Counterflow Premixed Ammonia/Air Flames

With the growing interest in ammonia as a carbon-free fuel, numerous recent research efforts on fundamental combustion characteristics of ammonia have been reported. Despite this attention, there remain fundamental aspects that lack experimental investigation or validated chemical kinetic models. Specifically, the understanding of ammonia flame structure and the quantification of extinction limits at elevated pressures are incomplete. Furthermore, the chemistry of NOx in relation to ammonia flames has not been thoroughly examined, especially at pressure, and it is useful to explore the operating conditions that exhibit favorable flame stability with minimal NOx formation. In the current work, key characteristics of ammonia/air combustion, including strain-induced extinction limits and NO* chemiluminescence, are experimentally investigated over a wide range of equivalence ratios and pressures. This paper presents initial results in a study directed to improved comprehension of ammonia combustion chemistry, and acquisition of new datasets under elevated pressure and temperature conditions, to better refine rate expressions for important reactions in ammonia combustion.

08 HYDROGEN↗

DOE OSTI2020

Dynamic Behavior of Natural Seep Vents: Analysis of Field and Laboratory Observations and Modeling (Final Scientific/Technical Report)

In this project, we have analyzed data collected by the U.S. Department of Energy (DOE), National Energy Technology Laboratory (NETL) in a high pressure water tunnel (HPWT) and data from two research cruises to natural seeps in the Gulf of Mexico to adapt and validate a numerical model to predict the dynamics of natural seeps in the deep oceans. The HPWT data include video observations of the shrinkage rate of individual methane and natural gas bubbles under simulated deep-water conditions. Field data were collected during two cruises by the Gulf Integrated Spill Research (GISR) Consortium led by Texas A&M University and funded by the Gulf of Mexico Research Initiative (GoMRI). These data included in situ observations from a remotely operated vehicle (ROV) of gas bubbles at two natural seep sites in the Gulf and acoustic observations of the natural seep bubble flares in the ocean water column. The acoustic data were from multibeam echosounders, one mounted in a forward-looking orientation on the ROV and another mounted down-looking in the haul of the ship. All of these laboratory and field data were focused on the dynamics of natural gas bubbles at temperatures and pressures favorable for clathrate hydrate formation between the gas and water. Our analyses of this data focused on understanding the mechanisms responsible for gas bubble dissolution within the hydrate stability zone (HSZ) of the oceans. Ice-like hydrate shells may form on the bubble-water interface under these conditions, and it was unknown how this might affect the mass transfer of gas into the ocean. We were able to extract bubble shrinkage rates from the HPWT datasets. Using this data we determined that mass transfer coefficients with and without a hydrate shell match empirical values for bubbles in contaminated systems (so-called dirty bubbles contaminated by naturally occurring surfactants). We also showed that free gas, and not gas hydrate, is the dominant dissolving phase when the hydrate sub-cooling is below 11 degree Celsius (temperature difference between hydrate the hydrate formation temperature and ambient temperator) or the pressure is reducing as bubbles rise through the ocean water column. Using this mass transfer model, our numerical model of bubble dissolution matched the over 200 HPWT experiments with an average error of 10% for predicting the bubble size at the end of an experiment. From field data in the literature, we also observed that gas bubbles dissolve faster when they are initially released, following mass transfer coefficients for so-called clean-bubbles (those not yet contaminated by surfactants). Shortly after release within the HSZ, a hydrate shell forms on the bubble-water interface, and the mass transfer reduces to rates matching those of dirty bubbles. We correlated this transition time from clean to dirty bubble behavior with the initial bubble surface area and the hydrate sub-cooling. With this model for hydrate formation time and using the mass transfer coefficients deduced from the HPWT data, we validated our numerical model for predicting the rise heights of natural seep flares in the oceans. Flare heights are commonly observed in haul-mounted acoustic multibeam data. The numerical model predicts bubbles to rise high in the ocean water column owing to the slower mass transfer rates for dirty bubbles that accompany the majority of their rise time. We found that the numerical model predictions matched the observed flare heights within 5% to 10% accuracy when we compared the rise heights of the largest bubbles released from the seafloor with the bubbles acoustically visible in the multibeam data. Bubbles become acoustically transparent as they shrink to sizes of order 1 mm in diameter for the multibeam frequencies used in the field. The forward-looking multibeam on the ROV also provided data on the lateral spreading of bubbles in natural seep flares. Our analysis of this data showed that spreading follows a diffusion process, with the effective diffusivity correlating with the wobbling length scale of these ellipsoidal bubbles. When we apply this diffusivity in a random displacement model of bubble spreading, our numerical simulations match closely the lateral spread observed by the M3 in the ocean water column. Finally, we compared the seep model predictions for the acoustic properties of these natural seep plumes with that observed by the acoustic instruments in the field. The M3 and EM 302 observations were converted to relative values of target strength using a calibration we obtained in the laboratory for the M3 and using an algorithm from the manufacturer for the EM 302. Comparing the numerical seep model to these data, we obtain good agreement over the whole height of rise of these bubble flares. This further validates the numerical model. Overall, our validated seep model captures the key dynamics of gas bubbles released from natural seeps in the oceans and helps to predict the fate of methane in the water column.

03 NATURAL GAS↗

DOE OSTI2021

Photodepositing CdS on the Active Cyano Groups Decorated g‐C 3 N 4 in Z‐Scheme Manner Promotes Visible‐Light‐Driven Hydrogen Evolution

Abstract g‐C 3 N 4 /CdS heterojunctions are potential photocatalysts for hydrogen production but their traditional type‐II configuration generally leads to weak oxidative and reductive activity. How to construct the novel Z‐scheme g‐C 3 N 4 /CdS counterparts to address this issue remains a great challenge in this field. In this work, a new direct Z‐scheme heterojunction of defective g‐C 3 N 4 /CdS is designed by introducing cyano groups (NC‐) as the active bridge sites. Experimental observations in combination with density functional theory (DFT) calculations reveal that the unique electron‐withdrawing feature of cyano groups in the defective g‐C 3 N 4 /CdS heterostructure can endow this photocatalyst with numerous advantageous properties including high light absorption ability, strong redox performance, satisfactory charge separation efficiency, and long lifetime of charge carriers. Consequently, the resultant photocatalytic system exhibits more active performance than CdS and g‐C 3 N 4 under visible light and reaches an excellent hydrogen evolution rate of 1809.07 µmol h −1 g −1 , which is 6.09 times higher than pristine g‐C 3 N 4 . Moreover, the defective g‐C 3 N 4 /CdS photocatalyst maintains good stability after 40 h continuous test. This work provides new insights into design and construction of Z‐scheme heterojunctions for regulating the visible‐light‐induced photocatalytic activity for H 2 evolution.

Wang, Zhipeng↗

DOE OSTI2021

Direct numerical simulation of supercritical oxy-methane mixing layers with CO 2 substituted counterparts

Direct numerical simulations (DNS) of temporally developing, three dimensional, CH 4 /CO 2 , CH 4 /O 2 and CO 2 /O 2 mixing layers, are conducted at a supercritical pressure of 300 atm. To effectively model the supercritical regime, the employed formulation includes the compressible form of the governing equations, the cubic Peng-Robinson equation of state and a generalized formulation for heat and mass flux vectors derived from non-equilibrium thermodynamics and fluctuation theory. A linear inviscid stability analysis is also performed for each case, to determine its most unstable wavelength. Flow visualizations reveal the presence of high density gradient magnitude regions for all three mixing layers, with conditional averages indicating increased presence of heavier fluid species within these regions. No significant departures are observed from perfect gas behavior, with compressibility factors very close to unity for all three mixing cases. Applicability of presumed probability density function methods (PDF) is examined for the three supercritical mixing layers. An a priori analysis is also conducted to investigate various simplifying assumptions employed in modeling various subgrid scale (SGS) flux models. Two additional terms are identified in the large eddy simulations (LES) equations, the gradient of SGS contribution of pressure in the momentum equation and the gradient of SGS contribution of heat flux in energy equation, whose magnitudes are similar and comparable with their respective resolved terms. The performance of the scale similarity model to represent these additional terms is investigated. Lastly, the performance of Smagorinsky, gradient and scale similarity models are also investigated.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2022

Thermal ion kinetic effects and Landau damping in fishbone modes

We report the kinetic–magnetohydrodynamic (MHD) hybrid simulation approach for macroscopic instabilities in plasmas can be extended to include the kinetic effects of both thermal ions and energetic ions. The new coupling scheme includes synchronization of the density and parallel velocity between thermal ions and MHD, in addition to pressure coupling, to ensure the quasineutrality condition and avoid numerical errors. The new approach has been implemented in the kinetic-MHD code M3D-C1-K, and was used to study the thermal ion kinetic effects and Landau damping in fishbone modes in both DIII-D and NSTX. It is found that the thermal ion kinetic effects can cause an increase of the frequencies of the non-resonant n = 1 fishbone modes driven by energetic particles for q min >1, and Landau damping can provide additional stabilization effects. A nonlinear simulation for n = 1 fishbone mode in NSTX is also performed, and the perturbation on magnetic flux surfaces and the transport of energetic particles are calculated.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2021

Manipulation of amorphous precursors to enhance zeolite nucleation

Crystallization in media comprised of amorphous precursors is becoming a more common phenomenon for numerous synthetic, biological, and natural materials that grow by a combination of classical and nonclassical pathways. Amorphous phases can exhibit a wide range of physicochemical properties that may evolve during the course of nucleation and crystal growth. This creates challenges for establishing causal relationships between amorphous precursor properties and their effect(s) on the selection of mechanistic pathways of crystallization and ultimately the properties of the crystalline product. In this study, we examine ways to manipulate the composition and colloidal stability of amorphous (alumino)silicate precursors that are prevalent in nanoporous zeolite syntheses. Changes in the amorphous precursor properties are evaluated on the basis of their ability to enhance rates of crystal formation. Here, we use fumed silica as the primary silicon source and examine the effects of infusing the source or growth medium with additional alkali metal, which serves as an inorganic structure-directing agent to facilitate the formation of porous crystal structures. We also assess the impact of adding a polymer additive, which reduces the colloidal stability of precursors, wherein we posit that the confined pockets of solution within the interstitial spaces of the precursor aggregates play an important role in regulating the rate of zeolite crystallization. Three commercially relevant zeolites (mordenite, SSZ-13, and ZSM-5) were selected for this study based on their diverse frameworks and methods of preparation. Here our findings reveal that alkali infusion significantly reduces the crystallization times for mordenite and SSZ-13, but has little impact on ZSM-5 synthesis. Conversely, we find that polymer addition markedly enhanced the rates of crystallization among all three zeolites, suggesting that this method may be a general approach to reduce zeolite synthesis times. Given the relatively high costs associated with commercial zeolite production, identifying new methods to improve the efficiency of hydrothermal syntheses can have significant practical implications beyond the fundamental benefits of developing new routes to tailor nonclassical crystallization.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2020

Turbulence model reduction by deep learning

A defining problem of turbulence theory is to produce a predictive model for turbulent fluxes. These have profound implications for virtually all aspects of the turbulence dynamics. In magnetic confinement devices, drift-wave turbulence produces anomalous fluxes via cross-correlations between fluctuations. In this work, we introduce an alternative, data-driven method for parametrizing these fluxes. The method uses deep supervised learning to infer a reduced mean-field model from a set of numerical simulations. We apply the method to a simple drift-wave turbulence system and find a significant new effect which couples the particle flux to the local gradient of vorticity. Notably, here, this effect is much stronger than the oft-invoked shear suppression effect. We also recover the result via a simple calculation. The vorticity gradient effect tends to modulate the density profile. In addition, our method recovers a model for spontaneous zonal flow generation by negative viscosity, stabilized by nonlinear and hyperviscous terms. We highlight the important role of symmetry to implementation of the new method.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2025

Model-Free Control of Grid-Interactive Efficient Buildings Under Communication Time Delays

Grid-interactive efficient buildings (GEBs) have recently been used to enhance the reliability and stability of the electric grid through demand response (DR) programs. However, most existing DR control strategies require accurate modeling of the various building thermostatically controlled loads (TCLs) and are computationally expensive. To address these challenges, a model-free control (MFC)-based strategy has recently been introduced for coordinating and controlling GEBs. MFC is a data-enabled control strategy that is computationally efficient and does not require the analytical models of the various building equipment. In this paper, we numerically investigate the impact of communication time delays on the performance of MFC in maintaining the TCLs' temperatures within the desired comfort levels while meeting the assigned power allocation constraint.

Telsang, Bhagyashri [University of Tennessee, Knox↗

DOE OSTI2024

Identification of Climatological Representative Days in the Mid-Atlantic for High-Fidelity Offshore Wind Energy Modeling

The goal of reaching 30 GW of offshore wind energy by 2030 becomes more realistic with the continued approval of offshore wind energy areas by the Biden Administration. In the Mid-Atlantic, where wind energy projects are in the most advanced stages of development, there is increased research focus on the eventual interaction of these wind farms. These interactions, in the form of wakes and cluster wakes, or wakes from multiple wind farms, could have detrimental effects on power production and forecastability for downwind wind farms (Pryor et al. 2022, Golbazi et al. 2022, Rosencrans et al. 2023). To help alleviate these issues, numerical simulations in the form of numerical weather prediction (NWP) and large eddy simulations (LES) can provide insight into when cluster wake situations may occur, but running such simulations can be expensive and difficult to run for multiple years. In this study, we leverage and build upon existing techniques in the literature (Fischereit et al. 2022) to identify climatologically representative days for wind energy areas in the Mid-Atlantic where conditions would promote cluster wake situations. We select meteorological variables (wind speed, wind direction, atmospheric stability, boundary-layer height, TKE) critical to understanding wind energy production and wake propagation. We then consider two different NWP datasets of varying spatial and temporal resolution: ERA5 provides data at hourly intervals from 1940 to present at 0.25 deg (31 km) spatial resolution (Hersbach et al. 2020), and the NOW-23 dataset provides data at 5-minute resolution for 21 years at 2-km spatial resolution (Bodini et al. 2020). Our first step is to compare these two datasets for an overlapping 21-year time period. Initial results show that the required number of days to represent the long-term climate increases with each additional variable considered. In their study of the German Bight, Fischereit et al. (2022) found that they could represent the long-term wind and wave climate in a "near-perfect" way with -180 days, by reaching a Perkins Skill Score (PSS) of 0.9; our investigation of the mid-Atlantic wind resource region with ERA5 and NOW-23 data suggests that we will need -100 days to reach a PSS of 0.9. As we expand our parameter space to include multiple variables, the number of required days will likely grow. These results will ultimately be used to select case studies to best represent cluster wake conditions that apply to this region for the lifetime of likely wind farms in this mid-Atlantic region.

clusterwakes↗

DOE OSTI2020

Super-time-stepping schemes for parabolic equations with boundary conditions

In this work, we present a super-time-stepping scheme for numerically solving parabolic partial differential equations with Dirichlet boundary conditions (BC). Using the general Forward Euler scheme, one can show that by taking varying step sizes there is the potential of propagating the solution forward in time by a greater amount than with uniform step sizes, while maintaining the same order of accuracy. As shown in previous works, if one further requires that the scheme have the Convex Monotone Property (CMP), then there exists a scheme which results in linear, monotone stability of the solution. This monotone stability is highly desirable in many physical situations, such as thermal diffusion, where the physical system will not oscillate, but will behave monotonically. However, the schemes devised in previous works do not include situations that have a boundary condition, and the inclusion of boundary conditions will henceforth be our focus. It is shown that a particular Runge-Kutta-Gegenbauer class of schemes [5] will maintain the CMP even in the presence of Dirichlet BC.

97 MATHEMATICS AND COMPUTING↗