Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “fixed grid”

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 523 records · Page 29
All sourcesNASA NTRSDOE OSTI
NASA NTRS1989

Euler/Navier-Stokes calculations of transonic flow past fixed- and rotary-wing aircraft configurations

Computational fluid dynamics has an increasingly important role in the design and analysis of aircraft as computer hardware becomes faster and algorithms become more efficient. Progress is being made in two directions: more complex and realistic configurations are being treated and algorithms based on higher approximations to the complete Navier-Stokes equations are being developed. The literature indicates that linear panel methods can model detailed, realistic aircraft geometries in flow regimes where this approximation is valid. As algorithms including higher approximations to the Navier-Stokes equations are developed, computer resource requirements increase rapidly. Generation of suitable grids become more difficult and the number of grid points required to resolve flow features of interest increases. Recently, the development of large vector computers has enabled researchers to attempt more complex geometries with Euler and Navier-Stokes algorithms. The results of calculations for transonic flow about a typical transport and fighter wing-body configuration using thin layer Navier-Stokes equations are described along with flow about helicopter rotor blades using both Euler/Navier-Stokes equations.

Deese, J. E.↗
DOE OSTI2017

Reference Computational Meshing Strategy for Computational Fluid Dynamics Simulation of Departure from Nucleate Boiling

The objective of this effort is to establish a strategy and process for generation of suitable computational mesh for computational fluid dynamics simulations of departure from nucleate boiling in a 5 by 5 fuel rod assembly held in place by PWR mixing vane spacer grids. This mesh generation process will support ongoing efforts to develop, demonstrate and validate advanced multi-phase computational fluid dynamics methods that enable more robust identification of dryout conditions and DNB occurrence.Building upon prior efforts and experience, multiple computational meshes were developed using the native mesh generation capabilities of the commercial CFD code STAR-CCM+. These meshes were used to simulate two test cases from the Westinghouse 5 by 5 rod bundle facility. The sensitivity of predicted quantities of interest to the mesh resolution was then established using two evaluation methods, the Grid Convergence Index method and the Least Squares method. This evaluation suggests that the Least Squares method can reliably establish the uncertainty associated with local parameters such as vector velocity components at a point in the domain or surface averaged quantities such as outlet velocity magnitude. However, neither method is suitable for characterization of uncertainty in global extrema such as peak fuel surface temperature, primarily because such parameters are not necessarily associated with a fixed point in space. This shortcoming is significant because the current generation algorithm for identification of DNB event conditions relies on identification of such global extrema. Ongoing efforts to identify DNB based on local surface conditions will address this challenge.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2015

A critical assessment of stellar mass measurement methods

This is the second paper in a series aimed at investigating the main sources of uncertainty in measuring the observable parameters in galaxies from their spectral energy distributions (SEDs). In the first paper we presented a detailed account of the photometric redshift measurements and an error analysis of this process. In this paper we perform a comprehensive study of the main sources of random and systematic error in stellar mass estimates for galaxies, and their relative contributions to the associated error budget. Since there is no prior knowledge of the stellar mass of galaxies (unlike their photometric redshifts), we use mock galaxy catalogs with simulated multi-waveband photometry and known redshift, stellar mass, age and extinction for individual galaxies. The multi-waveband photometry for the simulated galaxies were generated in 13 filters spanning from U-band to mid-infrared wavelengths. Given different parameters affecting stellar mass measurement (photometric signal-to-noise ratios (S/N), SED fitting errors and systematic effects), the inherent degeneracies and correlated errors, we formulated different simulated galaxy catalogs to quantify these effects individually. For comparison, we also generated catalogs based on observed photometric data of real galaxies in the Great Observatories Origins Deep Survey-South field, spanning the same passbands. The simulated and observed catalogs were provided to a number of teams within the Cosmic Assembly Near-infrared Deep Extragalactic Legacy Survey collaboration to estimate the stellar masses for individual galaxies. A total of 11 teams participated, with different combinations of stellar mass measurement codes/methods, population synthesis models, star formation histories, extinction and age. For each simulated galaxy, the differences between the input stellar masses, M{sub input}, and those estimated by each team, M{sub est}, is defined as Δlog(M)≡log(M{sub estimated})−log(M{sub input}), and used to identify the most fundamental parameters affecting stellar mass estimate in galaxies, with the following results. (1) No significant bias in Δ log(M) was found among different codes, with all having comparable scatter (σ(Δlog(M))=0.136 dex). The estimated stellar mass values are seriously affected by low photometric S/N, with the rms scatter increasing for galaxies with H{sub AB}>26 mag; (2) A source of error contributing to the scatter in Δ log(M) is found to be due to photometric uncertainties (0.136 dex) and low resolution in age and extinction grids when generating the SED templates; (3) The median of stellar masses among different methods provides a stable measure of the mass associated with any given galaxy (σ(Δlog(M))=0.142 dex); (4) The Δ log(M) values are strongly correlated with deviations in age (defined as the difference between the estimated and expected values), with a weaker correlation with extinction; (5) The rms scatter in the estimated stellar masses due to free parameters (after fixing redshifts and initial mass function) are quantified and found to be σ(Δlog(M))=0.110 dex; (6) Using the observed data, we studied the sensitivity of stellar masses to both the population synthesis codes and inclusion of nebular emission lines and found them to affect the stellar mass by 0.2 and 0.3 dex respectively.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Large-Scale Simulation of Regional Demand Flexibility Implementation and Customer Economic Impact

The Distribution System Operator with Transactive (DSO+T) study conducted a large-scale simulation of over 60,000 customers in a region the size of Texas to demonstrate the effective coordination of distributed energy resources (DERs) in commercial and residential buildings. The integrated simulation included both the bulk (wholesale generation and transmission) and distribution systems. The DERs (including batteries, electric vehicles, air conditioning, and water heaters) participated in a transactive energy retail market that was integrated into an existing competitive wholesale market. The engineering and economic performance of the resulting demand flexibility was evaluated over annual simulations for both moderate and high renewable generation scenarios. A detailed parametric cost model was developed to enable detailed economic analysis of key stakeholders. In addition, fixed and dynamic customer tariffs were designed and applied to the customer population. This allowed the impact on annual customer bills to be analyzed for various building types (residential versus commercial; single- versus multi-family). This paper presents results showing the relative flexibility of batteries, electric vehicles, and building loads throughout the year and under different renewable scenarios. This feeds a detailed breakdown of the impact this flexibility has on the operating costs of the grid and the resulting net economic benefit. Finally, the study showed that practically all customer classes (including non-participating customers) save money under the proposed demand flexibility scheme. The study found overall net annual economic savings of $3.3-5.0B for a region the size of Texas equating to average customer bill savings of 10-16%.

Reeve, Hayden M.↗
DOE OSTI2018

Description of the code ANVIL (ANisotropic Vhf Impulse Location)

In this paper, we have developed a general, three-dimensional method to locate sources of earth-directed radiation which takes into account the flight paths of reflected signals. While time-of-arrival algorithms exist for locating radio sources using line-of-sight propagation paths, radio sources originating above the surface of the earth will not necessarily emit strong power along direct paths to satellites. A combination of direct and ground-reflected pulses or only ground-reflected signals from such sources could be received by satellite-borne sensors. The work presented here applies to satellite detection of sub-ionospheric radio sources in a vacuum environment with ideal reflection off the earth's surface. Because satellites are not static receivers, their configuration with respect to a radio source is not always optimal. Therefore, a statistical study is performed using 1000 randomized configurations of 24 satellites in middle earth orbit. For each configuration, the radio source latitude and longitude is fixed, and its altitude is varied from 1 to 97 km. An analytic direct-path algorithm using the five satellites nearest the radio source provides an initial guess for the radio source latitude and longitude. We find that, using this approximation, the mean error in the calculated nadir location has a maximum value of 134 m for the 97-km source altitude. The coordinates of the initial guess are used to define the origin of a grid within which an all-points search refinement is performed. Finally, using this procedure, the overall maximum mean error in radio source position is found to be on the order of the computational grid size.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2015

Accurate determination of segmented X-ray detector geometry

Recent advances in X-ray detector technology have resulted in the introduction of segmented detectors composed of many small detector modules tiled together to cover a large detection area. Due to mechanical tolerances and the desire to be able to change the module layout to suit the needs of different experiments, the pixels on each module might not align perfectly on a regular grid. Several detectors are designed to permit detector sub-regions (or modules) to be moved relative to each other for different experiments. Accurate determination of the location of detector elements relative to the beam-sample interaction point is critical for many types of experiment, including X-ray crystallography, coherent diffractive imaging (CDI), small angle X-ray scattering (SAXS) and spectroscopy. For detectors with moveable modules, the relative positions of pixels are no longer fixed, necessitating the development of a simple procedure to calibrate detector geometry after reconfiguration. We describe a simple and robust method for determining the geometry of segmented X-ray detectors using measurements obtained by serial crystallography. By comparing the location of observed Bragg peaks to the spot locations predicted from the crystal indexing procedure, the position, rotation and distance of each module relative to the interaction region can be refined. Furthermore, we show that the refined detector geometry greatly improves the results of experiments.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
NASA NTRS2021

Agricultural Breadbaskets Shift Poleward Given Adaptive Farmer Behavior Under Climate Change

Modern food production is spatially concentrated in global “breadbaskets”. A major unresolved question is whether these peak production regions will shift poleward as the climate warms, allowing some recovery of potential climaterelated losses. While agricultural impacts studies to date have focused on currently cultivated land, the Global Gridded Crop Model Intercomparison Project (GGCMI) Phase 2 experiment allows us to assess changes in both yields and the location of peak productivity regions under warming. We examine crop responses under projected end-of-century warming using 7 process-based models simulating 5 major crops (maize, rice, soybeans, and spring and winter wheat) with a variety of adaptation strategies. We find that in no-adaptation cases, when planting date and cultivar choices are held fixed, regions of peak production remain stationary and yield losses can be severe, since growing seasons contract strongly with warming. When adaptations in management practices are allowed (cultivars that retain growing season length under warming and modified planting dates), peak productivity zones shift poleward and yield losses are largely recovered. While most growing-zone shifts are ultimately limited by geography, breadbaskets studied here move poleward over 600 km on average by end of the century under RCP8.5. These results suggest that agricultural impacts assessments can be strongly biased if restricted in spatial area or in the scope of adaptive behavior considered. Accurate evaluation of food security under climate change requires global modeling and careful treatment of adaptation strategies.

adaptation↗
DOE OSTI2016

SU-F-T-44: A Comparison of the Pre-Plan, Intra-Operative Plan, and Post-Implant Dosimetry for a Prostate Implant Case Using Prefabricated Linear Polymer-Encapsulated Pd-103

Purpose: To investigate the reproducibility and limitations of Pd-103 prostate brachytherapy using fixed length linear sources (CivaString). Methods: An LDR prostate brachytherapy case which was preplanned on MR images with prefabricated linear polymer-encapsulated Pd-103 sources (CivaString) was studied and compared with ultrasound based intra-operative planning and CT based post-implant dosimetry. We evaluated the following parameters among the three studies: prostate geometry (volume and cross sectional area), needle position and alignment deviations, and dosimetry parameters (D90). Results: The prostate volumes and axial cross sectional areas at center of prostate were measured as 41.8, 39.3 and 36.8 cc, and 14.9, 14.3, and 11.3 respectively on pre-plan MR, inter-op US, and post-implant CT studies. The deviation of prostate volumes and axial cross sectional areas measured on pre-planning MR and intra-operative US were within 5%. 17 out of 19 pre-planned needles were positioned within 5mm (the template grid size). One needle location was adjusted intra-operatively and another needle was removed due to proximity to urethra. The needle pathways were not always parallel to the trans-rectal probe due to the flexibility of CivaString. The angle of deviation was up to 10 degrees. Two pairs of needles were exchanged to better fit the length of prostate at the time of implant. This resulted in a prostate D90 of 153.8 Gy (124%) and 131.4 Gy (106.7%) for intra-op and PID respectively. Conclusion: Preplanning is a necessary part of implants performed with prefabricated linear polymer sources. However, as is often the case, there were real-time deviations from the pre-plan. Intra-op planning provides the ability conform to anatomy at the time of implant. Therefore, we propose to develop a systematic way to order extra strings of different length to provide the flexibility to perform intra-operative planning with fixed length strands.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2025

Scalable fabrication of an array-type fixed-target device for automated room temperature X-ray protein crystallography

X-ray crystallography is one of the leading tools to analyze the 3-D structure, and therefore, function of proteins and other biological macromolecules. Traditional methods of mounting individual crystals for X-ray diffraction analysis can be tedious and result in damage to fragile protein crystals. Furthermore, the advent of multi-crystal and serial crystallography methods explicitly require the mounting of larger numbers of crystals. To address this need, we have developed a device that facilitates the straightforward mounting of protein crystals for diffraction analysis, and that can be easily manufactured at scale. Inspired by grid-style devices that have been reported in the literature, we have developed an X-ray compatible microfluidic device that can be used to trap protein crystals in an array configuration, while also providing excellent optical transparency, a low X-ray background, and compatibility with the robotic sample handling and environmental controls used at synchrotron macromolecular crystallography beamlines. At the Stanford Synchrotron Radiation Lightsource (SSRL), these capabilities allow for fully remote-access data collection at controlled humidity conditions. Furthermore, we have demonstrated continuous manufacturing of these devices via roll-to-roll fabrication to enable cost-effective and efficient large-scale production.

chemical engineering↗
NASA NTRS1999

Using Spaceborne Ku-Band Scatterometer for Global Snow Cover Monitoring

We demonstrate for the first time the utility of spaceborne Ku-band scatterometer for global snow cover monitoring. Satellite radar data were collected over the globe by the NASA Scatterometer (NSCAT) operated at 14 GHz on board the Japanese ADEOS spacecraft from September 1996 to June 1997, spanning the 1997 seasonal snow season. First, we present backscatter signature of dry and wet snow to facilitate the interpretation of NSCAT backscatter evolution over snow cover regions. Surface field experiments indicated that dry snow backscatter at Ku band is approximately 40 times stronger than that at C band. Thus, Ku-band scatterometer measurements are sensitive to snow cover, which is typically transparent to C-band scatterometer returns. Furthermore, Ku-band backscatter does not saturate for most of natural snow depths as compared to radar responses at 19 GHz and 37 GHz or higher frequencies which have more limited penetration depths into snow. Ku-band backscatter is also sensitive to wetness in snow, which is appropriate to detect early snow melt conditions. Using the snow backscatter characteristics, we investigate NSCAT backscatter evolution over global snow cover regions throughout the 1997 snow season. The results reveal detail delineations between different regional snow areas. We show the correlation of these delineations with the boundaries of different global snow classes defined by the U.S. Army Cold Regions Research and Engineering Laboratory snow classification system. Using in-situ snow depth data from the U.S. National Climatic Data Center, we show that Ku-band backscatter corresponds very well to the trend of snow melt while snow mapping products (U.S. Climate Prediction Center gridded snow charts) from visible sensors does not reflect the fast snow melt trend. To illustrate the practical application of global snow monitoring with spaceborne Ku-band scatterometer, we present NSCAT backscatter response corresponding to the snow event leading to the 1997 Flood of the Century over the U.S. Northern plains and the Canadian prarie region, which caused loss of lives and several billion dollars in flood related damages and cleanup costs. Finally, we show that the fixed incidence configuration of Seawinds scatterometers, to be launched on QuickSCAT and on ADEOS-2 in the near future, is better for snow monitoring and then discuss the use of future high-resolution scatterometers for global snow mapping.

Nghiem, S. V.↗
NASA NTRS1996

Computation and Validation of Fluid/Structure Twin Tail Buffet Response

The buffet response of the flexible twin-tail/delta wing configuration, a multidisciplinary problem, is solved using three sets of equations on a multi-block grid structure. The first set is the unsteady, compressible, full Navier-Stokes equations which are used for obtaining the flow-field vector and the aerodynamic loads on the twin tails. The second set is the coupled aeroelastic equations which are used for obtaining the bending and torsional deflections of the twin tails. The third set is the grid-displacement equations which are used for updating the grid coordinates due to the tail deflections. The computational model is similar to the one used by Washburn et. al. which consists of a delta wing of aspect ratio one and twin tails with taper ratios of 0.23. The vortex of the twin tails are located at the wing trailing edge. The configuration is pitched at 30 deg angle of attack, and the freestream Mach number and Reynolds number are 0.3 and 1.25 million, respectively. With the twin tails fixed as rigid surfaces, the problem is solved for the initial flow conditions. Next, the problem is solved for the twin tail response for uncoupled bending and torsional vibrations due to the unsteady loads produced by the vortex breakdown flow of the leading-edge vortex cores. The configuration is investigated for three spanwise positions of the twin tails; inboard, midspan and outboard locations. The computational results are validated and are in very good agreement with the experimental data of Washburn, et. al.

Massey, Steven J.↗
DOE OSTI2010

Plug-in Electric Vehicle Interactions with a Small Office Building: An Economic Analysis using DER-CAM

It is generally believed that plug-in electric vehicles (PEVs) offer environmental and energy security advantages compared to conventional vehicles. Policies are stimulating electric transportation deployment, and PEV adoption may grow significantly. New technology and business models are being developed to organize the PEV interface and their interaction with the wider grid. This paper analyzes the PEVs' integration into a building's Energy Management System (EMS), differentiating between vehicle to macrogrid (V2M) and vehicle to microgrid (V2m) applications. This relationship is modeled by the Distributed Energy Resources Customer Adoption Model (DER-CAM), which finds optimal equipment combinations to meet microgrid requirements at minimum cost, carbon footprint, or other criteria. Results derive battery value to the building and the possibility of a contractual affiliation sharing the benefit. Under simple annual fixed payments and energy exchange agreements, vehicles are primarily used to avoid peak demand charges supplying cheaper off-peak electricity to the building during workdays.

24 POWER TRANSMISSION AND DISTRIBUTION↗
NASA NTRSDate not supplied

Validation Experiments of Incipient Turbulent Separation over an Axisymmetric Afterbody

Historically, the flow physics involved with most turbulent separated flows have presented fundamental challenges in validation between experimental and numerical approaches. As recognized by the CFD Vision 2030 study commissioned by NASA, validation of Reynolds-averaged Navier-Stokes (RANS) models and other scale-resolving methods for turbulent separated flow requires the support of advanced, high-fidelity experiments designed specifically for CFD implementation. In accordance with this effort, a new test platform, referred to as the NASA Axisymmetric Afterbody, was designed to obtain detailed information of the flow field over a smooth body, undergoing adverse pressure gradient induced separation for a fixed Reynolds number, Re = 180,000. The parametric body offers a range of flow states from fully attached,incipient separation, to fully separated flow based on variable afterbody geometries. In an initial effort to evaluate RANS turbulence model capabilities, the present configuration of the axisymmetric model features a mild adverse pressure gradient over the contoured boattail section, inducing incipient turbulent separation. Experiments include steady pressure measurements to serve as preliminary comparisons to simulation studies, which examine the effect of variable grid domains and RANS turbulence models, in an effort to understand and evaluate the critical variability between solutions for the present configuration. Results indicate potential discrepancies may be due to the effect of the square tunnel test section walls, relatively large blockage ratio, and inconsistent reference parameters. Ongoing work will focus on the experimental campaigns to obtain surface flow visualizations and high-resolution, off-body flow field measurements using Stereoscopic Particle Image Velocimetry (SPIV) and Laser Doppler Velocimetry (LDV) techniques to provide a detailed benchmark dataset to aid turbulence modelers.

axisymmetric afterbody↗
DOE OSTI2017

Exploring new topography-based subgrid spatial structures for improving land surface modeling

Topography plays an important role in land surface processes through its influence on atmospheric forcing, soil and vegetation properties, and river network topology and drainage area. Land surface models with a spatial structure that captures spatial heterogeneity, which is directly affected by topography, may improve the representation of land surface processes. Previous studies found that land surface modeling, using subbasins instead of structured grids as computational units, improves the scalability of simulated runoff and streamflow processes. In this study, new land surface spatial structures are explored by further dividing subbasins into subgrid structures based on topographic properties, including surface elevation, slope and aspect. Two methods (local and global) of watershed discretization are applied to derive two types of subgrid structures (geo-located and non-geo-located) over the topographically diverse Columbia River basin in the northwestern United States. In the global method, a fixed elevation classification scheme is used to discretize subbasins. The local method utilizes concepts of hypsometric analysis to discretize each subbasin, using different elevation ranges that also naturally account for slope variations. The relative merits of the two methods and subgrid structures are investigated for their ability to capture topographic heterogeneity and the implications of this on representations of atmospheric forcing and land cover spatial patterns. Results showed that the local method reduces the standard deviation (SD) of subgrid surface elevation in the study domain by 17 to 19 % compared to the global method, highlighting the relative advantages of the local method for capturing subgrid topographic variations. The comparison between the two types of subgrid structures showed that the non-geo-located subgrid structures are more consistent across different area threshold values than the geo-located subgrid structures. Altogether the local method and non-geo-located subgrid structures effectively and robustly capture topographic, climatic and vegetation variability, which is important for land surface modeling.

58 GEOSCIENCES↗
NASA NTRS1975

Radiation pattern synthesis of planar antennas using the iterative sampling method

A synthesis method is presented for determining an excitation of an arbitrary (but fixed) planar source configuration. The desired radiation pattern is specified over all or part of the visible region. It may have multiple and/or shaped main beams with low sidelobes. The iterative sampling method is used to find an excitation of the source which yields a radiation pattern that approximates the desired pattern to within a specified tolerance. In this paper the method is used to calculate excitations for line sources, linear arrays (equally and unequally spaced), rectangular apertures, rectangular arrays (arbitrary spacing grid), and circular apertures. Examples using these sources to form patterns with shaped main beams, multiple main beams, shaped sidelobe levels, and combinations thereof are given.

Stutzman, W. L.↗
NASA NTRS1999

Parallel Implementation of the Discontinuous Galerkin Method

This paper describes a parallel implementation of the discontinuous Galerkin method. Discontinuous Galerkin is a spatially compact method that retains its accuracy and robustness on non-smooth unstructured grids and is well suited for time dependent simulations. Several parallelization approaches are studied and evaluated. The most natural and symmetric of the approaches has been implemented in all object-oriented code used to simulate aeroacoustic scattering. The parallel implementation is MPI-based and has been tested on various parallel platforms such as the SGI Origin, IBM SP2, and clusters of SGI and Sun workstations. The scalability results presented for the SGI Origin show slightly superlinear speedup on a fixed-size problem due to cache effects.

Baggag, Abdalkader↗
DOE OSTI2017

High-order-harmonic generation of vibrating H 2 + and D 2 +

In this work, we study the isotope effect on high-order-harmonic generation (HHG) of vibrating H 2 + and D 2 + molecular ions aligned parallel to the polarization of the laser field. The time-dependent Schrödinger equation is solved accurately and efficiently by the time-dependent generalized pseudospectral and Fourier grid methods in three spatial coordinates, one of them being the internuclear separation and the other two describing the electronic motion. The laser pulses have a carrier wavelength of 800 nm and duration of 10 or 16 optical cycles. The peak intensities used in the calculations are 2×10 14 and 3×10 14 W/cm 2 . The effect of nuclear vibration is visible in both H 2 + and D 2 + but more pronounced in the lighter H 2 + molecule. Striking differences from the fixed nuclei case are a total disappearance of the traditional plateau in the HHG spectrum at the higher intensity and significant redshift of the harmonic peaks in the central part of the spectrum. These phenomena are explained based on the analysis of the dynamics of the nuclear vibrational wave packet.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2023

Intelligent Control of Thermal Energy Storage in the Manufacturing Sector for Plant-Level Grid Response

Industrial facilities are seeking new strategies that help in providing savings mechanisms for demand charges. Demand charges are the charges incurred by industrial facilities as a result of power usage. Thermal energy storage has advanced significantly with lots of new applications, garnering the interest of many industrial facilities. These applications could be used to shave the industrial facilities’ peak electric demand and reduce their demand charges. This paper aims to demonstrate the efficacy of thermal energy storage in reducing demand charges and highlight new developments in the integration of smart control systems with thermal energy storage. The study compares energy consumption and peak demand for a facility equipped with and without thermal energy storage tanks using a fixed schedule for charging and discharging. Additionally, the paper examines the impact of incorporating a smart controller to determine when to charge and discharge the tank based on the facility’s real-time power usage and a given setpoint. The results indicate cost savings from the use of thermal energy storage tanks under two proposed scenarios, reflected in the reduced cost of power consumption for the studied facility. The incorporation of a smart controller with the thermal energy storage tank in the facility studied could provide estimated savings of 3.3% per year of power consumption charges, without considering the contribution of any incentives. The estimated savings provided by the fixed schedule scenario are 2.7% per year.

25 ENERGY STORAGE↗