SEARCH · Engineering Papers
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.
Boundary element method (BEM) and WEC-Sim analysis of UMass Dartmouth's maximal asymmetric drag wave energy converter (MADWEC), including its tethered ballast system and PTO (power take-off).
A split simulation whole-building energy/three-dimensional earth contact model (termed the BUFETS/EnergyPlus Model or BEM) can model the full range of foundation systems in the target retrofit housing stock. This model was used to extensively evaluate the results produced by the NREL BEoptE+ tool. These foundation systems that include above-grade foundation walls, diabatic floors or slabs, and lookout or walkout walls could not be modeled in BEopt version 1.1 that was used in this project.
Infrared (IR) imagery can provide valuable data to building energy auditors and the developers of Building Energy Models (BEM). IR by itself is not a panacea to the challenges facing energy auditors and modelers, but when combined with other techniques and Open Studio software, IR can play an important role in the minimization of the burden of creating building energy models and the diagnosing of building energy problems. While IR can be used to detect and quantify many aspects of a building’s thermal and energy performance, in Phase I we focused on a single capability: determining the R value of exterior walls. Two key attributes of using IR methods for R value determination is that IR methods measure the actual, installed performance in-situ – as opposed to relying on manufacturer or handbook values – and they can integrate performance over a wide area and thus can accurately account for thermal bridging and spatial anomalies such as “spotty” insulation installation. There is a trade between user burden and accuracy when using IR to determine a wall’s R value. The least burdensome approach is computing R values from a single IR snapshot. The most burdensome method we considered was to combine IR images with thermocouple and heat flux sensor (HFS) measurements made over the course of several days. The more effort users expend in making the measurements, the more accurate and robust the method becomes. Through a series of field tests at the ThermoAnalytics Inc. (TAI) HQ at Calumet, MI and a simulation of a 1980s era office building in Flagstaff, AZ we investigated methods to compute the R value of exterior walls using data from IR cameras, thermocouples and heat flux sensors. To compute R values the heat flux through the wall, along with temperatures on either side of the wall, must be known, modeled, or estimated – and the temperature difference through the wall should be at least 10 C. If only outdoor IR images are used in the method, the inside ambient air temperature must be estimated, and a simulation model created to predict heat fluxes and the indoor wall temperature. The accuracy of the resulting R value depends on the accuracy of the simulation and the assumptions it employs. Applying thermocouple sensors, either as iButtons or as part of a wireless sensor network, to the walls allows for the acquisition of temperatures over time. The longer the duration of the measurements, the more accurate and robust the method will be. Even with this level of measurement effort, modeling is still required to calculate the heat flux through the wall. The addition of a heat flux sensor alleviates the need for any modeling, and that results in improved accuracy. Driving the accuracy of the methods are several factors that can be quantified. How well are the wall surface temperatures known? Modeling the heat flux through walls requires knowing the details of the surrounding radiant environment and the heat transfer due to convection. How well are these boundary conditions known? Even when the heat flux is measured, error can arise when some of the measured heat is not conducted through the wall but rather goes into the thermal mass of the wall. Thermal mass effects occur when conditions and the heat flux values are fluctuating. A technique that is not robust will suffer in accuracy when these fluctuations occur. The R value determination method must be incorporated into the Open Studio process. In our proposed process, IR images and other measured sensor data are input into a module that would send R value estimates to a building library developed by team member Skidmore, Owings, & Merrill (SOM). This library of parameterized building models would use these R values along with information about the building’s age, type, and location to query the Open Studio BCL to identify the most likely constructions for walls, roof, and windows. Being based on the IMF feature of Energyplus, the SOM building library can quickly produce Energyplus models of buildings based on a limited number of user inputs.
It is well known that airfoils under unsteady flow conditions with a periodically varying angle of attack exhibit aerodynamic characteristics different from those under steady flow conditions, a phenomenon commonly known as dynamic stall. It is also well known that the steady aerodynamic characteristics of airfoils in the inboard region of a rotating blade differ from those under steady two-dimensional (2D) flow conditions, a phenomenon commonly known as rotational augmentation. This paper presents an investigation of these two phenomena together in the inboard parts of wind turbine blades. This analysis is carried out using data from three sources: (1) the National Renewable Energy Laboratory’s Unsteady Aerodynamics Experiment Phase VI experimental data, including constant as well as continuously pitching blade conditions during axial operation, (2) data from unsteady Delayed Detached Eddy Simulations (DDES) carried out using the Technical University of Denmark’s in-house flow solver Ellipsys3D, and (3) data from a simplified model based on the blade element momentum method with a dynamic stall subroutine that uses rotationally augmented steady-state polars obtained from steady Phase VI experimental sequences, instead of the traditional 2D nonrotating data. The aim of this work is twofold. First, the blade loads estimated by the DDES simulations are compared to three select cases of the N sequence experimental data, which serves as a validation of the DDES method. Results show reasonable agreement between the two data in two out of three cases studied. Second, the dynamic time series of the lift and the moment polars obtained from the experiments are compared to those from the dynamic stall subroutine that uses the rotationally augmented steady polars. This allowed the differences between the stall phenomenon on the inboard parts of harmonically pitching blades on a rotating wind turbine and the classic dynamic stall representation in 2D flow to be investigated. Results from the dynamic stall subroutine indicated a good qualitative agreement between the model and the experimental data in many cases, which suggests that the current 2D dynamic stall model as used in BEM-based aeroelastic codes may provide a reasonably accurate representation of three-dimensional rotor aerodynamics when used in combination with a robust rotational augmentation model.
A split simulation whole building energy / 3-dimensional earth contact model (termed the BUFETS/EnergyPlus Model or BEM) capable of modeling the full range of foundation systems found in the target retrofit housing stock has been extensively tested. These foundation systems that include abovegrade foundation walls, diabatic floors or slabs as well as lookout or walkout walls, currently cannot be modeled within BEopt.
CCBlade predicts aerodynamic loading of wind turbine blades using blade element momentum (BEM) theory. CC stands for continuity and convergence. CCBlade was developed primarily for use in gradient-based optimization applications where C1 and robust convergence are essential.
The U.S. Department of Energy (DOE) supports the development of commercial building energy codes and standards by participating in industry review, update processes, and providing technical analyses to support both published model codes and potential changes. DOE publishes its findings to ensure transparency in its support, and to make its analysis available for both public review and use. In conjunction with this effort, DOE’s flagship building energy model (BEM) tools consist of the open-source EnergyPlus simulation engine and OpenStudio. These tools allow a user to modify a building, estimate energy use, and new releases are typically downloaded by over 20,000 users.
In the recent years, there has been a tremendous interest to numerical methods that approximate partial differential equations (PDEs) on computational meshes with arbitrarily-shaped polygonal/polyhedral (polytopal, for short) elements. A nonexhaustive list of such methods include the Mimetic Finite Difference method [4, 9, 21, 24, 28, 33, 35, 62, 64–67], the Polygonal Finite Element Method [68, 73, 74], the polygonal Discontinuous Galerkin Finite Element Methods [5, 7, 8, 10, 13–15, 18, 36, 38, 39, 41, 50], the Hybridizable Discontinuous Galerkin and Hybrid High-Order Methods [46–48, 51], the Gradient Discretization method [49, 53, 54], the Finite Volume Method [52], and the BEM-based FEM [71]. An alternative approach that also proved to be very successful is provided by the Virtual Element method (VEM), which was originally proposed in [19] for the numerical treatment of second-order elliptic problems [43, 44], and readily extended to linear and nonlinear elasticity [22, 25, 56], plate bending problems [37], Cahn-Hilliard equation [3], Stokes equations [2], Laplace-Beltrami equation [55], Darcy-Brinkam equation [76], discrete topology optimization problems [6], fracture networks problems [29], eigenvalue problems [59, 75]. The mixed virtual element formulation was proposed in [20, 34]; the nonconforming Virtual element formulations was proposed for second-order elliptic problems in [16], and later extended to general advection-reaction-diffusion problems, Stokes equation, the biharmonic problems, the iegenvalue problem, and the Schrodinger equation in [11, 30, 42, 45, 58, 80]. The p- and hp-version of the VEM were developed in [26, 27] and efficient multigrid methods for the resulting linear system of equations in [12]. A posteriori error estimates can be found in [40].
The INTEGRATE (Inverse Network Transformations for Efficient Generation of Robust Airfoil and Turbine Enhancements) project is developing a new inverse-design capability for the aerodynamic design of wind turbine rotors using invertible neural networks. This AI-based design technology can capture complex non-linear aerodynamic effects while being 100 times faster than design approaches based on computational fluid dynamics. This project enables innovation in wind turbine design by accelerating time to market through higher-accuracy early design iterations to reduce the levelized cost of energy. INVERTIBLE NEURAL NETWORKS Researchers are leveraging a specialized invertible neural network (INN) architecture along with the novel dimension-reduction methods and airfoil/blade shape representations developed by collaborators at the National Institute of Standards and Technology (NIST) learns complex relationships between airfoil or blade shapes and their associated aerodynamic and structural properties. This INN architecture will accelerate designs by providing a cost-effective alternative to current industrial aerodynamic design processes, including: - Blade element momentum (BEM) theory models: limited effectiveness for design of offshore rotors with large, flexible blades where nonlinear aerodynamic effects dominate - Direct design using computational fluid dynamics (CFD): cost-prohibitive - Inverse-design models based on deep neural networks (DNNs): attractive alternative to CFD for 2D design problems, but quickly overwhelmed by the increased number of design variables in 3D problems AUTOMATED COMPUTATIONAL FLUID DYNAMICS FOR TRAINING DATA GENERATION - MERCURY FRAMEWORK The INN is trained on data obtained using the University of Marylands (UMD) Mercury Framework, which has with robust automated mesh generation capabilities and advanced turbulence and transition models validated for wind energy applications. Mercury is a multi-mesh paradigm, heterogeneous CPU-GPU framework. The framework incorporates three flow solvers at UMD, 1) OverTURNS, a structured solver on CPUs, 2) HAMSTR, a line based unstructured solver on CPUs, and 3) GARFIELD, a structured solver on GPUs. The framework is based on Python, that is often used to wrap C or Fortran codes for interoperability with other solvers. Communication between multiple solvers is accomplished with a Topology Independent Overset Grid Assembler (TIOGA). NOVEL AIRFOIL SHAPE REPRESENTATIONS USING GRASSMAN SPACES We developed a novel representation of shapes which decouples affine-style deformations from a rich set of data-driven deformations over a submanifold of the Grassmannian. The Grassmannian representation as an analytic generative model, informed by a database of physically relevant airfoils, offers (i) a rich set of novel 2D airfoil deformations not previously captured in the data , (ii) improved low-dimensional parameter domain for inferential statistics informing design/manufacturing, and (iii) consistent 3D blade representation and perturbation over a sequence of nominal shapes. TECHNOLOGY TRANSFER DEMONSTRATION - COUPLING WITH NREL WISDEM Researchers have integrated the inverse-design tool for 2D airfoils (INN-Airfoil) into WISDEM (Wind Plant Integrated Systems Design and Engineering Model), a multidisciplinary design and optimization framework for assessing the cost of energy, as part of tech-transfer demonstration. The integration of INN-Airfoil into WISDEM allows for the design of airfoils along with the blades that meet the dynamic design constraints on cost of energy, annual energy production, and the capital costs. Through preliminary studies, researchers have shown that the coupled INN-Airfoil + WISDEM approach reduces the cost of energy by around 1% compared to the conventional design approach. This page will serve as a place to easily access all the publications from this work and the repositories for the software developed and released through this pr...
Electric water heaters represent 14% of the electricity consumption in residential buildings. An average household in the United States (U.S.) spends about USD 400–600 (0.45 ¢/L–0.68 ¢/L) on water heating every year. In this context, water heaters are often considered as a valuable asset for Demand Response (DR) and building energy management system (BEMS) applications. To this end, this study proposes a model-free deep reinforcement learning (RL) approach that aims to minimize the electricity cost of a water heater under a time-of-use (TOU) electricity pricing policy by only using standard DR commands. In this approach, a set of RL agents, with different look ahead periods, were trained using the deep Q-networks (DQN) algorithm and their performance was tested on an unseen pair of price and hot water usage profiles. The testing results showed that the RL agents can help save electricity cost in the range of 19% to 35% compared to the baseline operation without causing any discomfort to end users. Additionally, the RL agents outperformed rule-based and model predictive control (MPC)-based controllers and achieved comparable performance to optimization-based control.
WEC-Sim (Wave Energy Converter SIMulator) is an opensource software for modeling the motions, loads and power generation of wave energy converters. WEC-Sim performs simulations in the time domain using hydrodynamic coefficients calculated by boundary element method (BEM) frequency-domain potential flow solvers such as WAMIT, NEMOH, Capytaine, or Ansys AQWA. WEC-Sim development is ongoing, including various features, applications, and example cases used to demonstrate potential use-cases to meet the needs of the growing marine energy industry. Through input from a broad user base and an extensive team of developers and collaborators, new features of WEC-Sim are developed to expand the software’s use cases and improve overall functionality. Three new WEC-Sim features highlighted in this paper include updating WEC-Sim to be compatible with MoorDyn Version 2, incorporation of second order excitation loads (quadratic transfer functions) and allowing for dynamically changing hydrodynamics.
Abstract. Changes in land cover and aerosols resulting from urbanization may impact convective clouds and precipitation. Here we investigate how Houston urbanization can modify sea-breeze-induced convective cloud and precipitation through the urban land effect and anthropogenic aerosol effect. The simulations are carried out with the Chemistry version of the Weather Research and Forecasting model (WRF-Chem), which is coupled with spectral-bin microphysics (SBM) and the multilayer urban model with a building energy model (BEM-BEP). We find that Houston urbanization (the joint effect of both urban land and anthropogenic aerosols) notably enhances storm intensity (by ∼ 75 % in maximum vertical velocity) and precipitation intensity (up to 45 %), with the anthropogenic aerosol effect more significant than the urban land effect. Urban land effect modifies convective evolution: speed up the transition from the warm cloud to mixed-phase cloud, thus initiating surface rain earlier but slowing down the convective cell dissipation, all of which result from urban heating-induced stronger sea-breeze circulation. The anthropogenic aerosol effect becomes evident after the cloud evolves into the mixed-phase cloud, accelerating the development of storm from the mixed-phase cloud to deep cloud by ∼ 40 min. Through aerosol–cloud interaction (ACI), aerosols boost convective intensity and precipitation mainly by activating numerous ultrafine particles at the mixed-phase and deep cloud stages. This work shows the importance of considering both the urban land and anthropogenic aerosol effects for understanding urbanization effects on convective clouds and precipitation.
Abstract. Urban overheating and its ongoing exacerbation due to global warming and urban development lead to increased exposure to urban heat and increased thermal discomfort and heat stress. To quantify thermal stress, specific indices have been proposed that depend on air temperature, mean radiant temperature (MRT), wind speed, and relative humidity. While temperature and humidity vary on scales of hundreds of meters, MRT and wind speed are strongly affected by individual buildings and trees and vary on the meter scale. Therefore, most numerical thermal comfort studies apply microscale models to limited spatial domains (commonly representing urban neighborhoods with building blocks) with resolutions on the order of 1 m and a few hours of simulation. This prevents the analysis of the impact of city-scale adaptation and/or mitigation strategies on thermal stress and comfort. To solve this problem, we develop a methodology to estimate thermal stress indicators and their subgrid variability in mesoscale models – here applied to the multilayer urban canopy parameterization BEP-BEM within the Weather Research and Forecasting (WRF) model. The new scheme (consisting of three main steps) can readily assess intra-neighborhood-scale heat stress distributions across whole cities and for timescales of minutes to years. The first key component of the approach is the estimation of MRT in several locations within streets for different street orientations. Second, mean wind speed and its subgrid variability are downscaled as a function of the local urban morphology based on relations derived from a set of microscale LES and RANS simulations across a wide range of realistic and idealized urban morphologies. Lastly, we compute the distributions of two thermal stress indices for each grid square, combining all the subgrid values of MRT, wind speed, air temperature, and absolute humidity. From these distributions, we quantify the high and low tails of the heat stress distribution in each grid square across the city, representing the thermal diversity experienced in street canyons. In this contribution, we present the core methodology as well as simulation results for Madrid (Spain), which illustrate strong differences between heat stress indices and common heat metrics like air or surface temperature both across the city and over the diurnal cycle.
Building Energy Modeling (BEM) is a versatile tool for designing, retrofitting, ensuring code compliance, obtaining certifications, qualifying for incentives, and enabling real-time building control. However, capturing all the details of building geometry for thermal zoning can be time-consuming, costly, and sometimes computationally challenging. As a result, modelers have been applying zoning simplification based on factors such as space functions and internal loads, as well as relying on their experience and judgment while adhering to zoning rules outlined in industry standards. Despite the prevalence of this common practice, a notable gap exists in the literature regarding studies quantifying the influence of simplified thermal zoning on the evaluation of Energy Efficiency Measures (EEMs). Recognizing this gap, this paper seeks to contribute to the field by enhancing the understanding of how the simplification of thermal zoning influences the evaluation of EEMs against a baseline design. The study utilized a medium office prototype model with a detailed floor plan featuring over 20 zones per floor covering diverse functional spaces with varying internal loads and occupancy schedules. A standard thermal zoning strategy outlined in ASHRAE Standard 90.1 Appendix G was employed as the simplified zoning method. This strategy condenses the zoning into a core zone and four perimeter zones per floor. It was compared with the detailed zoning approach, which involves one zone per space. Common Energy EEMs, such as enhanced envelope, high-efficiency appliances and equipment, and HVAC controls, were individually implemented and evaluated. The results indicate that the performance comparison between the two zoning methods varies depending on the type of measures considered. Basic measures, such as adding wall insulation, demonstrate similar energy impacts, while advanced HVAC control measures, such as static pressure reset, exhibit a more substantial difference that cannot be overlooked.
The inclusion of body forces has received a good deal of attention in boundary element research. The consideration of such forces is essential in the desgin of high performance components such as fan and turbine disks in a gas turbine engine. Due to their critical performance requirements, optimal shapes are often desired for these components. The boundary element method (BEM) offers the possibility of being an efficient method for such iterative analysis as shape optimization. The implicit-differentiation of the boundary integral equations is performed to obtain the sensitivity equations. The body forces are accounted for by either the particular integrals for uniform body forces or by a surface integration for non-uniform body forces. The corresponding sensitivity equations for both these cases are presented. The validity of present formulations is established through a close agreement with exact analytical results.
The intent here is to discuss the status of the boundary element method (BEM) for structural analysis, both in terms of the present and anticipated capabilities of the method and in terms of the incorporation of the method in the design/analysis process, particularly for gas turbine engine components.
Recent advances in FEM techniques and applications are discussed in reviews and reports presented at the ICASE/LaRC workshop held in Hampton, VA in July 1986. Topics addressed include FEM approaches for partial differential equations, mixed FEMs, singular FEMs, FEMs for hyperbolic systems, iterative methods for elliptic finite-element equations on general meshes, mathematical aspects of FEMS for incompressible viscous flows, and gradient weighted moving finite elements in two dimensions. Consideration is given to adaptive flux-corrected FEM transport techniques for CFD, mixed and singular finite elements and the field BEM, p and h-p versions of the FEM, transient analysis methods in computational dynamics, and FEMs for integrated flow/thermal/structural analysis.
This paper presents a generalized numerical technique for evaluating the optimal performance of active noise controllers. In this technique, the indirect BEM numerical procedures are used to derive the active noise controllers for optimal control of enclosed harmonic sound fields where the strength of the noise sources or the description of the enclosure boundary may not be known. The performance prediction for a single-input single-output system is presented, together with the analysis of the stability and observability of an active noise-control system employing detectors. The numerical procedures presented can be used for the design of both the physical configuration and the electronic components of the optimal active noise controller.