Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “BEM”

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 127 records · Page 7

Design Space Exploration for Novel Reduced-Vortex Turbine Rotors Using Free Vortex Wake Methods

In this work, we explore the use of mid-fidelity free vortex wake (FVW) tools, namely NREL's cOnvecting LAgrangian Filaments (OLAF) tool, to design and analyze wind turbine rotor blades. The goal of the exercise is to determine whether designs exist that out-perform traditional blade designs, which are optimized using the Blade-Element Momentum (BEM) method. We find that using design of experiments to generate simple modifications of the blade tip, OLAF predicts higher performance in terms of power coefficient for redesigned blades that are non-optimal according to BEM. We then use optimization methods with OLAF in the loop to attempt to automatically find OLAF-optimal designs. We analyze the results of these optimizations, demonstrating that some roadblocks exist before OLAF and FVW tools can be used to automatically find FVW-optimal blade designs that might exceed the performance of BEM-optimal blade designs.

blade design↗

Heating Ventilation Air-Conditioning (HVAC) System Analysis and Optimization for Revit with the OpenStudio Analytical Framework (CRADA Final Report)

This project proposes to integrate the NREL-developed OpenStudio (OS) and its accompanying Analytical Framework (OSAF) into Autodesk’s (ADSK’s) market leading suite of building information modeling (BIM) tools. As a result, BIM users will gain access to advanced (but user-friendly) building energy modeling (BEM), calibration, and design optimization capabilities. Following integration, the BIM software will automatically generate a consistent OS-based BEM model, allowing users to develop projects in an integrated manner, moving back and forth between architectural and engineering perspectives without redundancy or data loss, saving time and effort.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

Calibrating Building Energy Models Using Supercomputer Trained Machine Learning Agents

Building Energy Modeling (BEM) is an approach to model the energy usage in buildings for design and retrofit purposes. EnergyPlus is the flagship Department of Energy software that performs BEM for different types of buildings. The input to EnergyPlus can often extend in the order of a few thousand parameters which have to be calibrated manually by an expert for realistic energy modeling. This makes it challenging and expensive thereby making building energy modeling unfeasible for smaller projects. In this paper, we describe the Autotune research which employs machine learning algorithms to generate agents for the different kinds of standard reference buildings in the U.S. building stock. The parametric space and the variety of building locations and types make this a challenging computational problem necessitating the use of supercomputers. Millions of EnergyPlus simulations are run on supercomputers which are subsequently used to train machine learning algorithms to generate agents. These agents, once created, can then run in a fraction of the time thereby allowing cost-effective calibration of building models.

Sanyal, Jibonananda↗

A scoping review of non-destructive testing (NDT) techniques in building performance diagnostic inspections

Understanding building envelope thermodynamics is an essential foundation of building sciences, mainly due to the envelope’s role as a boundary layer for exterior environments, as well as a container and regulator of internal microclimates. This paper presents a scoping literature review of select Non-destructive Testing (NDT) techniques for building envelope scanning and surveying for thermodynamic diagnostics. The investigation focuses specifically on reviewing six NDT techniques: Ground Penetrating Radar (GPR), Light Detection and Ranging (LiDAR)/Laser Scanning, Thermography, Ultrasound, Close-Range Photogrammetry and Through Wall Imaging Radar (TWIR). The aim is to identify knowledge gaps in terms of their use in accurately characterizing envelope compositions for further integration in Building Energy Modeling (BEM). Each technique was evaluated according to set categories imbibed from the American Society of Heating, Refrigerating, and Air Conditioning Engineering (ASHRAE) Standard 211P that showcase the technique’s ability to extract various relevant information. A framework is then developed to inform users on how to use hybrid NDT-based workflows applied in building envelope energy audits. The study concludes by discussing possibilities of utilizing NDT in large-scale audit automation, BEM integration, and developing built environment policies focusing on increasing existing building performance through retrofitting design.

36 MATERIALS SCIENCE↗

Building envelope anomaly characterization and simulation using drone time-lapse thermography

Defects in building envelopes deteriorate over time without being visible to the human eye, while significantly impacting energy performance due to unaccounted heat transfer. Defects can be characterized in the infrared (IR) spectrum. However, IR readings are typically recorded at singular points in time, when in several cases anomalies can only be revealed at specific times of the day, possibly in different seasons of the year. This paper presents a novel workflow for 3D envelope defect characterization and modeling using aerial time-lapse IR data collection using drones. A comprehensive envelope thermal profile is developed for a case study building employing the photogrammetry software Agisoft Photoscan, which generates temporal IR inspections of building skins using multiple thermography orthomosaics. Point-cloud data is then translated into a CAD model and thermal zones for whole Building Energy Modeling (BEM) using Honeybee as a frontend to EnergyPlus to showcase the potential of inclusion of detailed 4D data. Envelope contributions in this case study’s anomalies showed heat losses of 6447.6 kWh, and Energy Use Intensity (EUI) differences of ~2 kWh/m 2 /year from the baseline. Finally, why there is currently little translation of this work in BEM software is discussed, while identifying limitations and future research in the employment of time-lapse thermography using drones for more accurate building envelope inspection and modeling.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

A 2.5D boundary element formulation for modeling damped waves in arbitrary cross-section waveguides and cavities

Highlights: •Dispersive properties of viscoelastic waveguides and cavities are computed using a regularized 2.5D BEM. •Linear viscoelasticity is introduced at the constitutive level by means of frequency dependent complex moduli. •A contour integral algorithm is used to solve the nonlinear eigenvalue problem. •The Sommerfeld radiation condition is used to select the permissible Riemann sheets. •Attenuation of surface waves in cavities approaches the attenuation of Rayleigh waves. -- Abstract: A regularized 2.5D boundary element method (BEM) is proposed to predict the dispersion properties of damped stress guided waves in waveguides and cavities of arbitrary cross-section. The wave attenuation, induced by material damping, is introduced using linear viscoelastic constitutive relations and described in a spatial manner by the imaginary component of the axial wavenumber. The discretized dispersive wave equation results in a nonlinear eigenvalue problem, which is solved obtaining complex axial wavenumbers for a fixed frequency using a contour integral algorithm. Due to the singular characteristics and the multivalued feature of the wave equation, the requirement of holomorphicity inside the contour region over the complex wavenumber plane is fulfilled by the introduction of the Sommerfeld branch cuts and by the choice of the permissible Riemann sheets. A post processing analysis is developed for the extraction of the energy velocity of propagative guided waves. The reliability of the method is demonstrated by comparing the results obtained for a rail and a bar with square cross-section with those obtained from a 2.5D Finite Element formulation also known in literature as Semi Analytical Finite Element (SAFE) method. Next, to show the potential of the proposed numerical framework, dispersion properties are predicted for surface waves propagating along cylindrical cavities of arbitrary cross-section. It is demonstrated that the attenuation of surface waves approaches asymptotically the attenuation of Rayleigh waves.

97 MATHEMATICS AND COMPUTING↗

A hybrid Boundary Element Unstructured Transmission-line (BEUT) method for accurate 2D electromagnetic simulation

Time domain electromagnetic simulation tools have the ability to model transient, wide-band applications, and non-linear problems. The Boundary Element Method (BEM) and the Transmission Line Modeling (TLM) method are both well established numerical techniques for simulating time-varying electromagnetic fields. The former surface based method can accurately describe outwardly radiating fields from piecewise uniform objects and efficiently deals with large domains filled with homogeneous media. The latter volume based method can describe inhomogeneous and non-linear media and has been proven to be unconditionally stable. Furthermore, the Unstructured TLM (UTLM) enables modelling of geometrically complex objects by using triangular meshes which removes staircasing and unnecessary extensions of the simulation domain. The hybridization of BEM and UTLM which is described in this paper is named the Boundary Element Unstructured Transmission-line (BEUT) method. It incorporates the advantages of both methods. The theory and derivation of the 2D BEUT method is described in this paper, along with any relevant implementation details. The method is corroborated by studying its correctness and efficiency compared to the traditional UTLM method when applied to complex problems such as the transmission through a system of Luneburg lenses and the modelling of antenna radomes for use in wireless communications. - Graphical abstract:.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

X-ray CT characterization and fracture simulation of ASR damage of glass particles in alkaline solution and mortar

The alkali-silica reaction (ASR) damage in reactive aggregates affects the long-term durability of the concrete infrastructure. The generated ASR gel can expand by imbibing water from the pore solution and the resulting expansion pressure causes the aggregate fracture. This study aims to simulate the development of ASR damage in glass particles of two types of samples (glass in alkali solution and glass mortars). The dynamic micro X-ray CT technique was conducted to monitor the crack propagation in glass aggregates at different reaction stages (up to 64 hrs). The Boundary Element Method (BEM) and Displacement Discontinuity Method (DDM) were used to efficiently simulate the crack propagation within irregular glass particles under gel expansion pressures. The aggregate boundaries were built with the BEM elements and the initial cracks were meshed with the DDM elements. The estimated expansion pressure was applied to the initial crack surfaces. The discontinuous displacements and stresses were calculated along the crack path and crack tip. The mixed-mode Stress Intensity Factors (SIFs) were calculated based on plane stress conditions. The maximum circumferential stress criteria was used to simulate the propagation of cracks along a specific angle. The simulation results include the simulated crack path and the combined SIFs changing with the increments. With the estimate expansion pressure, the glass particle damages were simulated within conditions of alkali solution and confined motar samples. The predicted crack propagation path was compared with the X-ray CT imaging data. The comparison results demonstrate the DDM has the ability to predicting the ASR damage propagation inside aggregates.

36 MATERIALS SCIENCE↗

Multiple scattering of surface waves by cavities in a half-space

Scattering of surface waves from multiple two-dimensional cavities at the surface of a homogenous, isotropic, linearly elastic half-space is analyzed in this work. For the case of multiple cavities, the scattered field is shown to be equivalent to the total radiation from the distributions of tractions, calculated from the incident wave, over the surfaces of the cavities. The multiple-scattering model is obtained from known single-scattering calculation for a cavity by the use of the self-consistent method. The second order approximation to the multiple-scattering problem by a random distribution of cavities is then considered and solved analytically. The vertical displacement at some distance from the cavities is calculated and verified by the solution of the same problem obtained by the boundary element method (BEM). The analytical and BEM results are graphically displayed and show good agreement when the depths of the cavities are small compared to the wavelength.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Effects of head modeling errors on the spatial frequency representation of MEG

Abstract Objectives. We aim to investigate the effects of head model inaccuracies on signal and source reconstruction accuracies for various sensor array distances to the head. This allows for the assessment of the importance of head modeling for next-generation magnetoencephalography (MEG) sensors, optically-pumped magnetometers (OPM). Approach. A 1-shell boundary element method (BEM) spherical head model with 642 vertices of radius 9 cm and conductivity of 0.33 S m −1 was defined. The vertices were then randomly perturbed radially up to 2%, 4%, 6%, 8% and 10% of the radius. For each head perturbation case, the forward signal was calculated for dipolar sources located at 2 cm, 4 cm, 6 cm and 8 cm from the origin (center of the sphere), and for a 324 sensor array located at 10 cm to 15 cm from the origin. Equivalent current dipole (ECD) source localization was performed for each of these forward signals. The signal for each perturbed spherical head case was then analyzed in the spatial frequency domain, and the signal and ECD errors were quantified relative to the unperturbed case. Main results. In the noiseless and high signal-to-noise ratio (SNR) case of approximately ≥6 dB, inaccuracies in our spherical BEM head conductor models lead to increased signal and ECD inaccuracies when sensor arrays are placed closer to the head. This is true especially in the case of deep and superficial sources. In the noisy case however, the higher SNR for closer sensor arrays allows for an improved ECD fit and outweighs the effects of head geometry inaccuracies. Significance. OPMs may be placed directly on the head, as opposed to the more commonly used superconducting quantum interference device sensors which must be placed a few centimeters away from the head. OPMs thus allow for signals of higher spatial resolution to be captured, resulting in potentially more accurate source localizations. Our results suggest that an increased emphasis on accurate head modeling for OPMs may be necessary to fully realize its improved source localization potential.

Yeo, Wan-Jin (ORCID:0000000241417123)↗

New Developments and Capabilities Within WEC-Sim

WEC-Sim is an open-source software for simulating wave energy converters and has been actively developed and applied since its initial release in 2014 to simulate a wide variety of device archetypes. WEC-Sim is developed jointly by the National Renewable Energy Laboratory and Sandia National Laboratories within the MATLAB/SIMULINK environment. A general wave-to-wire model begins with a deployment site resource characterization, which is used to complete the hydrodynamic simulation of wave energy converters (WEC), with the power generation profile imported to a grid simulator to understand the influence on the local electrical network. While modeling the entire wave-to-wire is difficult and encompasses multiple time scales and physics, WEC-Sim is focused on the hydrodynamics simulation to predict, analyze, and optimize WEC dynamics and power performance. WEC-Sim simulations are performed in the time domain based on the radiation and diffraction method using hydrodynamics coefficients derived from boundary element method (BEM)-based frequency-domain potential flow solvers (e.g., WAMIT, NEMOH, Capytaine, or ANSYS-AQWA). With this level of modeling fidelity, WEC-Sim can handle floating body hydrodynamics, mechanical and electrical power generation methods, advanced control implementation, mooring systems, and other unique applications such as desalination. Additional WEC-Sim functionalities include pre-built Simulink blocks and MATLAB scripts that can simulate a wide range of floating systems and the corresponding auxiliary subsystems. The developers of WEC-Sim continue to release new versions of the software, at least annually, with the latest release in September 2022. These releases include bug fixes, updates to software documentation, as well as new features to expand WEC-Sim's capabilities to model a wide range of WEC concepts. This publication will highlight the new features added to WEC-Sim between versions 4.1.0 to 5.0.1, which spans a 2-year period from June 2020 to September 2022. New features described here include topics such as continuous integration checks, revised Morison Element and nonlinear hydro implementations, run directly from Simulink (required for hardware-in-the-loop execution), BEMIO updates to import Capytaine BEM hydrodynamics, addition of cable blocks, and new wave visualization features.

TIDAL AND WAVE POWER↗

pnnl/openstudio-building-energy-standard-measures-gem

The OpenStudio Building Energy Standards Measure Repository is a comprehensive collection of innovative OpenStudio measures designed to facilitate energy code analysis. OpenStudio is a cutting-edge platform brings together physics-based building energy modeling (BEM), BEM process automation, and large-scale computing capabilities. Contains OpenStudio Standards Performance Rating Method Measure (IPID 32936) and OpenStudio Measure to Generate IPLV-specific Chiller Performance Curves for Chillers (IPID 32940)

Xu, Weili [PNNL]↗

Deep Design Data Portal (D3P) v0.01

The Deep Design Data Portal (D3P) tool was developed to demonstrate how readily accessible data sources, such as building energy model reports for design and baseline energy performance data for projects, can provide the data required for reporting to an industry initiative (AIA 2030 commitment), as well as more detailed data that makes the industry dataset more valuable to all stakeholders, enabling project level analysis and analysis of BEM industry trends. D3P provides an easier and less time-consuming way for firms to auto-extract data from this data source, compared to the current reporting workflows of the firms. The BEM reports are the first of several data sources that D3P could integrate. D3P also provides the ability for firms to review, compare, and evaluate the performance of their projects to not only their portfolio, but also to the larger anonymized industry dataset created each time a project is added to D3P. The intent of D3P is to become part of a data-sharing ecosystem to assist creating large anonymized industry datasets that are accessible to industry.

Regnier, Cynthia [Lawrence Berkeley National Labor↗

Towards Digital and Performance-Based Supervisory HVAC Control Delivery

Upgrading supervisory HVAC control in commercial buildings is one of the most attractive decarbonization tools at our disposal. Modern controls are software programs and can in theory be deployed at scale and with a low up-front carbon "pulse". In practice, however, control delivery is a disjointed and inefficient process, dominated by manual handoffs of imprecise English language documents. A particularly high barrier exists between control implementation and building energy modeling (BEM) which results in control sequences typically not being tested for correctness or performance before implementation. Together with industry partners, DOE and the national labs are developing an ecosystem of tools and standards that can support fully digital performance-based control delivery workflows. This paper describes this ecosystem, which consists of three mutually supportive efforts. Semantic models of buildings and their systems enable automatic configuration and installation of control software. Platform-neutral control descriptions separate control algorithms from control platforms and enable the creation of libraries of reference control implementations. Dynamic whole-building energy-control simulation that can execute physically realistic control sequences makes it possible to test and evaluate the performance of control sequences and then directly compile them for installation and execution in control systems. In addition to digitizing and streamlining project-level control delivery, these standards and related software support benchmarking of control algorithms, both rule-based and optimization-based, and help to both advance the state of the art and to implement ratings and programs that encourage the adoption of high-performance control.

building controls↗

Simplified Performance Rating Method - Review of Existing Tools, Rulesets, and Programs. Literature and technical review

The Performance Rating Method (PRM) in ASHRAE Standard 90.1 (ASHRAE 2016) is a simulation ruleset for establishing minimum code compliance and for rating a building’s beyond-code applications. Building energy modeling (BEM) to comply with the PRM is often expensive and time consuming due to the complexity of code requirements. The goal of this project is to develop and codify a simplified PRM (S-PRM) approach for commercial buildings that will expand the use of BEM for small or simple buildings by defining a low-cost, simplified approach for creating robust and detailed models. This report summarizes a review of (i) the more commonly used existing codes, (ii) incentive programs, and (iii) the existing simplified energy modeling tools to identify which requirements are good candidates for simplification.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

Improved Air-Conditioning Demand Response of Connected Communities over Individually Optimized Buildings

Connected communities potentially offer much greater demand response capabilities over singular building energy management systems (BEMS) through an increase of connectivity. The potential increase in benefits from this next step in connectivity is still under investigation, especially when applied to existing buildings. This work utilizes EnergyPlus simulation results on eight different commercial prototype buildings to estimate the potential savings on peak demand and energy costs using a mixed-integer linear programming model. This model is used in two cases: a fully connected community and eight separate buildings with BEMS. The connected community is optimized using all zones as variables, while the individual buildings are optimized separately and then aggregated. These optimization problems are run for a range of individual zone flexibility values. The results indicate that a connected community offered 60.0% and 24.8% more peak demand savings for low and high flexibility scenarios, relative to individually optimized buildings. Energy cost optimization results show only marginally better savings of 2.9% and 6.1% for low and high flexibility, respectively.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗

Progress in the validation of rotor aerodynamic codes using field data

Within the framework of the fourth phase of the International Energy Agency (IEA) Wind Task 29, a large comparison exercise between measurements and aeroelastic simulations has been carried out featuring three simulation cases in axial, sheared and yawed inflow conditions. Results were obtained from more than 19 simulation tools originating from 12 institutes, ranging in fidelity from blade element momentum (BEM) to computational fluid dynamics (CFDs) and compared to state-of-the-art field measurements from the 2 MW DanAero turbine. More than 15 different variable types ranging from lifting-line variables to blade surface pressures, loads and velocities have been compared for the different conditions, resulting in over 250 comparison plots. The result is a unique insight into the current status and accuracy of rotor aerodynamic modeling. For axial flow conditions, a good agreement was found between the various code types, where a dedicated grid sensitivity study was necessary for the CFD simulations. However, compared to wind tunnel experiments on rotors featuring controlled conditions, it remains a challenge to achieve good agreement of absolute levels between simulations and measurements in the field. For sheared inflow conditions, uncertainties due to rotational and unsteady effects on airfoil data result in the CFD predictions standing out above the codes that need input of sectional airfoil data. However, it was demonstrated that using CFD-synthesized airfoil data is an effective means to bypass this shortcoming. For yawed flow conditions, it was observed that modeling of the skewed wake effect is still problematic for BEM codes where CFD and free vortex wake codes inherently model the underlying physics correctly. The next step is a comparison in turbulent inflow conditions, which is featured in IEA Wind Task 47. Doing this analysis in cooperation under the auspices of the IEA Wind Technology Collaboration Program (TCP) has led to many mutual benefits for the participants. The large size of the consortium brought ample manpower for the analysis where the learning process by combining several complementary experiences and modeling techniques gave valuable insights that could not be found when the analysis is carried out individually.

17 WIND ENERGY↗

Numerical Modeling of a Two-Stage Ocean Current Turbine

The Equinox Ocean Turbines (EQOT) current energy converter has a unique design with power generation in two small-diameter turbines attached to the tips of a large-diameter passive rotor. This configuration offers some key advantages for capturing ocean currents. With no centrally placed generator, almost no reaction torque is required at the nacelle of the main large-diameter rotor, and the small-diameter tip turbine generators operate at a higher speed and lower torque. The physics that determine the performance and loads on the turbine are also unique. The interactions of the flow field between the two stages and the general architecture of the system cannot be captured with traditional mid-fidelity modeling tools. For design iterations and large sets of load cases, it is important to have mid-fidelity models that can capture the important phenomenon with enough accuracy to identify global trends. This work uses a limited set of high-fidelity computational fluid dynamics (CFD) simulations to help inform the selection of and construction of a custom mid-fidelity model. Mid-fidelity modeling approaches were verified by comparing key turbine performance quantities to those found with the CFD model. Hydrodynamic interactions of the two-stage rotor were identified through high-fidelity CFD modeling. This highlighted the impact of the main rotor tip vortex and wake on the secondary rotor apparent inflow. This results in a relative flow rotation and sharp deficit, that change the optimal secondary rotor rotation speed and adds unsteadiness to the blade loading respectively. Multiple mid-fidelity approaches were evaluated for their ability to capture these effects. A simple approximation of the combined-stage performance based on single-stage BEM provides a reasonable rough prediction, especially near the peak TSR values, with some larger discrepancy at higher TSRs. Predicting the combined-stage performance based on single-stage CFD data improves this prediction across the TSR range. Although the combined-stage modeling in OLAF was not successful in this stage of the project, it showed promise as a mid-fidelity method, assuming the parameters can be tuned to account for the significant differences in time and length scales between the main and secondary rotors. This may be addressed through code changes in future work. A significant finding from the OLAF work was the agreement between the vortex core radius values found independently via a parameter space search and via CFD. The technique of using single-stage secondary rotor BEM, with a custom inflow taken from single-stage main rotor CFD or OLAF, provides an efficient method to capture one-way coupled flow interactions. This method provided generally good predictions of the impact of the flow rotation on the secondary rotor but struggled to accurately predict the peaks of the unsteady load progression. Future work could include some superposition of a tuned main rotor trailing edge viscous wake into the custom inflow to better predict this interaction.

16 TIDAL AND WAVE POWER↗