Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Building performance simulation”

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 1,189 records · Page 66
All sourcesNASA NTRSDOE OSTI
NASA NTRS1996

Pointing and Scanning Control of Optical Instruments using Rotating Unbalanced Masses

Correct pointing direction and scanning motions are essential in the operation of many flight payloads, such as balloon-borne telescopes and space-based X- ray and gamma-ray telescopes. Rotating unbalanced mass (RUM) devices have been recently proposed, implemented and successfully tested to produce a variety of scanning motions. Linear scans, raster scans, and circular scans have been successfully generated on a gimbaled payload using pairs of RUM devices. Theoretical analysis, computer simulations, and experiments have also been used to study the feasibility of using RUM devices to control instrument pointing direction, in addition to generating scanning motion. Dynamic modeling of a gimbaled payload equipped with a pair of RUM devices has been studied, and preliminary testing indicates that the pointing control is indeed feasible. However, there is also great potential for significant performance improvements through more advanced control system analysis, modeling and design. In this paper, modeling and control methods are described to achieve simultaneous scanning and pointing control of a gimbaled payload using rotating unbalance mass (RUM) devices. The model development work builds upon the results of Polites et al. and also some modeling approaches from robotics research. Results of some preliminary experiments are discussed and some nonlinear control methods will be proposed.

Bishop, Carlee A.↗
NASA NTRS2000

High Energy Flywheel Containment Evaluation

A flywheel testing facility is being constructed at the NASA Glenn Research Center. This facility is to be used for life cycle testing of various flywheel rotors. The lifecycle testing consists of spinning a rotor from a low rpm (approx. 20,000 ) to a high rpm (approx. 60,000) and then back to the low rpm. This spin cycle will model that which the rotor will see during use. To simulate the lifetime of the rotor, the spin cycle will be performed tens of thousands of times. A typical life cycle spin test is expected to last six months. During this time the rotor will be spun through a cycle every five minutes. The test will run continuously for the six month period barring a flywheel failure. Since it is not reasonable to have the surrounding area evacuated of personnel for the duration of the testing, the flywheel facility has to be designed to withstand a flywheel rotor failure and insure that there is no danger to any personnel in the adjacent buildings or surrounding areas. In order to determine if the facility can safely contain a flywheel rotor failure an analysis of the facility in conjunction with possible flywheel failure modes was performed. This analysis is intended as a worst case evaluation of the burst liner and vacuum tank's ability to contain a failure. The test chamber consists of a cylindrical stainless steel vacuum tank, two outer steel containment rings, and a stainless steel burst liner. The stainless steel used is annealed 302, which has an ultimate strength of 620 MPa (90,000 psi). A diagram of the vacuum tank configuration is shown. The vacuum tank and air turbine will be located below ground in a pit. The tank is secured in the pit with 0.3 m (12 in.) of cement along the base and the remaining portion of the tank is surrounded by gravel up to the access ports. A 590 kg (1300 lb.) bulkhead is placed on top of the pit during operation and the complete facility is housed within a concrete structure which has 7.5 cm (3 in.) thick walls. A cutaway of the facility is shown.

Colozza, Anthony J.↗
NASA NTRS2018

Risk-Based Explosion Assessment Internal Explosion of Exploration Upper Stage in Vehicle Assembly Building

This paper presents risk-based explosion analyses used to determine if an explosion of NASA's Space Launch System (SLS) Exploration Upper Stage (EUS) during fueling in Vehicle Assembly Building (VAB) High Bay 4 can hazard the Space Launch System (SLS) being processed in Bay 3. This paper focuses on blast fragmentation as it was the controlling hazard. An estimate of the amount of Net Explosion Weight (NEW) was used to develop a fragment list (mass, shape, takeoff velocity/angle) given an explosion occurs at ground level in Bay 4. Intervening superstructure that exists between Bay 4 and Bay 3 (that could potentially block EUS fragments) was modeled based on detailed visual images (VISSIM) and AUTOCAD/structural drawings. The fragment list was input to a 3D trajectory application that performs multiple random explosion simulations; each simulation varying fragment ballistic coefficient and takeoff angle/velocity. Individual fragments are tracked until stopped by the intervening superstructure or reach Bay 3, potentially impacting the SLS. The probability of impacting the SLS is then determined from the random simulations; the results are used to help decision-makers determine the efficacy of co-processing within the VAB and decide what mitigations may be appropriate.

Explosion Assessment↗
DOE OSTI2022

QUIC-URB and QUIC-fire extension to complex terrain: Development of a terrain-following coordinate system

Ensemble-based approaches to prescribed fire planning cannot be supported by CFD-based models like FIRETEC and WFDS because they are too computationally expensive and cannot leverage LES approaches like CAWFE and WRF-SFIRE because too coarse of resolution. QUIC-Fire was developed to fill this gap but it cannot currently address complex terrain, typical for instance of the Western United States. In this paper, we describe the extension of the diagnostic wind model QUIC-URB, the wind engine of QUIC-Fire, to a terrain-following coordinate system. In particular, the paper presents the mathematical derivation of the wind solver leading to a linear system of equations that are solved through the successive over-relaxation method. The model is validated against a standard test used in previous works (the Askervein Hill) and against a new dataset from measurements in the Socorro Mountains, New Mexico. The terrain-following implementation captured the correct phenomenology for the isolated Askervein Hill, with a wind speed up at the top of the hill. We report the model agreed well with measurements on the upwind side of the peak, but overestimated speed-up on the downwind side of the hill. This is due to the inability of the model to generate flow separation and wake-eddy dynamics. On a common laptop, the divergence-free wind field was obtained in 6 s, making the solver appealing for coupled fire–atmosphere simulations. The Socorro Mountain was highly complex, with many cliff faces, peaks, and valleys. Although the model captures the magnitude and direction of inlet and outlet areas of the domain, it performs rather poorly in the valley region and in the regions near the steep cliffs. Hence, the model shows good agreement with data in areas of open sloped terrain but lacks in areas where flow separation and thermally driven effects may be present (neither effect was addressed in this work). Results highlight that future work should focus on the implementation of parameterizations of wake-eddies, similar to QUIC-URB’s building parameterizations, and on thermodynamic-driven flow.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

Phase change material-to-refrigerant heat exchangers: Experimental validation and uneven melting analysis

Thermal energy storage (TES) using phase change materials (PCMs) enables load shifting and reduces mismatches between the building's thermal demand and the heat pump (HP) system's thermal capacity. While most TES studies have focused on single-phase heat transfer fluids, PCM-to-refrigerant heat exchangers (PRHX), in which the refrigerant undergoes liquid-vapor phase change, remain largely underexplored. Here, to address this gap, this study systematically investigates a shell-and-tube PRHX operating as the condenser in an HP through experiments and simulations. A two-dimensional enthalpy-based finite-volume PRHX model was developed and validated against experimental data, including refrigerant and PCM temperature profiles, PRHX capacity, and system-level performance. Results show three characteristic stages: Stage I, the condenser temperature rapidly increases with PCM in the solid state absorbing sensible heat; Stage II, stable operation during PCM phase change; and Stage III, performance decline once part of the PCM becomes fully melted. Heat transfer analysis revealed that uneven PCM melting along the condenser length was driven primarily by variations in the refrigerant-side heat transfer coefficient, rather than local approach temperature differences. Further study demonstrated that minimizing refrigerant outlet subcooling led to more uniform PCM melting and extended Stages I and II by up to 89%. This paper provides a validated PRHX model, clarifies the mechanisms of uneven PCM melting, and highlights subcooling control as an effective strategy to improve TES performance in HP systems without requiring secondary loops.

Graphite matrix↗
DOE OSTI2025

Multi-physics melt pool modeling and process optimization for laser direct energy deposition of Nb-based refractory C103: Defect formation, geometric precision, and process mapping

Recent developments in additive manufacturing (AM) technology have reignited interest in the fabrication of the Nb-based refractory C103 alloy offering solutions to the challenges posed by traditional manufacturing methods. However, the limited numerical and experimental studies on laser direct energy deposition (DED) of C103 have hindered the understanding of the relationships between process parameters and build quality. This has made it challenging to consistently produce parts with the desired quality and microstructure suitable for critical applications. In this study, we focus on optimizing the laser DED process for C103 by employing a hybrid approach that combines experimental techniques and computational fluid dynamics (CFD). This approach facilitates the development of process maps for defect detection and geometric precision. To achieve this, multi-layer C103 samples were fabricated using laser DED under various process parameters, enabling the creation of a process map for defect detection. Additionally, a multi-physics, multiphase simulation framework was developed within a high-performance computing (HPC) environment to establish process maps for geometric precision. Using these process maps, printability windows were identified for achieving both the desired geometric accuracy and defect-free prints. It was observed that prints with a power-to-velocity (P/V) ratio close to unity resulted in defect-free outcomes. This study provides a foundation for reducing design lead time and rejected parts, ultimately optimizing the laser DED process for C103.

Defect formation and geometric precision↗
DOE OSTI2026

The CP-PAW Code Package for First-Principles Calculations from a User’s Perspective

CP-PAW is a combined electronic structure and ab initio molecular dynamics code to perform mixed quantum and classical simulations of atomistic condensed phase systems, such as solids, liquids, and molecular systems. As the name suggests, the CP-PAW code unifies the all-electron projector augmented-wave (PAW) method with the Car–Parrinello (CP) approach to determine not only the electronic and nuclear ground states of condensed matter but also to study their properties and dynamics. In addition to briefly outlining the underlying theory, the focus will be on the unique aspects of CP-PAW and how to correctly employ them as a user. How to install CP-PAW using the new build system will also be briefly mentioned.

Blöchl, Peter E [Institute for Theoretical Physic↗
DOE OSTI2021

Potential design problems for ITER fusion device

The international thermonuclear experimental reactor (ITER) is a worldwide project currently being built in France for the demonstration of the feasibility of thermonuclear technologies for future realization of successful commercial fusion energy. ITER is of the tokamak based design using strong magnetic fields to confine the very hot plasma needed to induce the fusion reaction. Tokamak devices are currently the front leading designs. Building a successful magnetic fusion device for energy production is of great challenge. A key obstacle to such design is the performance during abnormal events including plasma disruptions and so-called edge-localized modes (ELMs). In these events, a massive and sudden release of energy occurs quickly, due to loss of full or partial plasma confinement, leading to very high transient power loads on the reactor surface boundaries. A successful reactor design should tolerate several of these transient events without serious damages such as melting and vaporization of the structure. This paper highlights, through comprehensive state-of-the-art computer simulation of the entire ITER interior design during such transient events, e.g., ELMs occurring at normal operation and disruptions during abnormal operation, potential serious problems with current plasma facing components (PFCs) design. The HEIGHTS computer package is used in these simulations. The ITER reactor design was simulated in full and exact 3D geometry including all known relevant physical processes involved during these transient events. The current ITER divertor design may not work properly and may requires significant modifications or new innovative design to prevent serious damage and to ensure successful operation.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

High frequency beam oscillation keyhole dynamics in laser melting revealed by in-situ x-ray imaging

The metal additive manufacturing industry is actively developing instruments and strategies to enable higher productivity, optimal build quality, and controllable as-built microstructure. A beam controlling technique, laser oscillation has shown potential in all these aspects in laser welding; however, few attempts have been made to understand the underlying physics of the oscillating keyholes/melt pools which are the prerequisites for these strategies to become a useful tool for laser-based additive manufacturing processes. Here, to address this gap, we utilized a synchrotron-based X-ray operando technique to image the dynamic keyhole oscillation in Ti-6Al-4V using a miniature powder bed fusion setup. We found good agreement between the experimental observations and simulations performed with a validated Lattice Boltzmann multiphysics model. The study revealed the continuous and periodic fluctuations in the characteristic keyhole parameters that are unique to the oscillating laser beam processing and responsible for the chevron pattern formation at solidification. In particular, despite the intrinsic longer-range fluctuation, the oscillating technique displayed potential for reducing keyhole instability, mitigating porosity formation, and altering surface topology. These insights on the oscillating keyhole dynamics can be useful for the future development and application of this technique.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Automated Control for Nuclear Thermal Propulsion Start-Up using MOOSE-based Applications

This report presents a Griffin/Bison/RELAP-7 numerical model of a prototypical NTP system that features fuel assemblies arranged in rings, and which was designed to simulate rapid startup transients. The physics modeled include full-core neutronics, assembly-wise heat conduction, and conjugate heat transfer, with the balance of plant mainly imposed through boundary conditions. In addition, various forms of automated reactivity control were deployed by using the MOOSE to autonomously drive the model and simulate the reactor transitioning from assumed initial conditions to nominal power in a fraction of a minute. To generate the cross-sections of the neutronics model, and in an effort to simultaneously account for the tremendous axial temperature gradients in the reactor and to limit the number of state points required for cross-section generation, the average component temperatures and hydrogen densities in the cooling channels were correlated to the average fuel and moderator temperatures, and fixed axial profiles were derived for nominal conditions and then used during the transient. With this approximation, a tractable cross-section library tabulated with fuel/moderator temperatures and CD angles was generated using Serpent. The full-core SPH correction procedure and the CD decusping technology in Griffin, respectively, ensure preservation of the multiplication factor and reaction rates at state points, along with a reasonably accurate reactivity worth between tabulated CD angles, despite using a coarse mesh. Feedback from other physics was calculated by modeling one representative fuel assembly per ring, along with the corresponding fuel and moderator cooling channels. To limit power overshoots during startup, another layer of multiphysics coupling was added to the model in order to automatically control the drums. Two different technologies presented herein showed outstanding performance in this regard: (1) a novel hybrid PID controller based on both power and reactivity signals, and (2) a PGC that relies on kinetics parameters and reactivity coefficients to predict future behavior and adjust the desired signal accordingly. A challenging benchmark was devised, featuring a power demand curve that exponentially increases by a factor of 500 within 30 seconds, then levels out after that. Both control approaches create a simulated power curve that closely follows the power demand curve and limits power overshoots to 1% or less. While the former approach requires more tuning of the internal parameters, the latter requires additional knowledge of the reactivity feedback coefficients and rates of change of the corresponding variables, including fuel and moderator temperature, which could be difficult to dynamically measure for a real NTP system. Fortunately, some inaccuracy in these quantities will not drastically degrade the PGC performance. Subsequently, a more realistic startup sequence was considered, in which the mass flow rate and outlet pressures are ramped up to model bootstrap and thrust build-up phases prior to reaching steady-state conditions, demonstrating the ability of the hybrid PID and PGCs to handle such transients, with both types of controllers exhibiting very similar behavior. Nevertheless, a significant chamber temperature overshoot was observed, caused by the demanded power signal and assumed mass flow rate. This issue could be mitigated by deploying a reactor controller that follows the chamber temperature signal and actuates both the control valves and drums (rather than using a power signal based solely on the drums to control reactivity). Enhancement of the hydrogen fluid properties available in MOOSE, as well as a better understanding of prototypical initial conditions, are also needed to further enhance this startup model. Finally, a study was performed to model decay heat post-shutdown, and to prepare for extending this model to predict shutdown behavior and post-shutdown pulsed cooling requirements.

33 ADVANCED PROPULSION SYSTEMS↗
NASA NTRSDate not supplied

New Methods to Determine the Electromagnetic Effects from Precipitation Static

The Electrostatics and Surface Physic Laboratory has devised a new method to determine the electromagnetic effects from precipitation static or p-static as its commonly called. P-static occurs as a result of triboelectric charging of granular materials or ice crystals bombarded onto the surfaces of air and/or space vehicles at high speeds. The electrostatic charge that builds up on the outer surface can have deleterious effects that arise during discharge events. These discharges can cause loss of communication to and from the vehicle from ground stations, EMI interference, radio frequency noise, etc… There are even cases of total loss of vehicles caused by p-static charging. We have developed a method to simulate the extensive tribocharging of insulating materials that occurs during p-static through the use of corona charging which allows testing to occur within a lab environment. These tests simulate the discharges that occur on rocket vehicles during ascent as a result of gas breakdown by performing the corona charging tests inside a vacuum chamber and subsequently removing the air. The simultaneous gas discharges emit electromagnetic interference (EMI) which can be captured and monitored.

Charles R Buhler↗
DOE OSTI2020

High-Performance Integrated Window and Façade Solutions for California

The researchers developed a new generation of high-performance façade systems and supporting design and management tools to support industry in meeting California’s greenhouse gas reduction targets, reduce energy consumption, and enable an adaptable response to minimize real-time demands on the electricity grid. The project resulted in five outcomes: (1) The research team developed an R-5, 1-inch thick, triplepane, insulating glass unit with a novel low-conductance aluminum frame. This technology can help significantly reduce residential cooling and heating loads, particularly during the evening. (2) The team developed a prototype of a windowintegrated local ventilation and energy recovery device that provides clean, dry fresh air through the façade with minimal energy requirements. (3) A daylight-redirecting louver system was prototyped to redirect sunlight 15–40 feet from the window. Simulations estimated that lighting energy use could be reduced by 35–54 percent without glare. (4) A control system incorporating physics-based equations and a mathematical solver was prototyped and field tested to demonstrate feasibility. Simulations estimated that total electricity costs could be reduced by 9-28 percent on sunny summer days through adaptive control of operable shading and daylighting components and the thermostat compared to state-of-the-art automatic façade controls in commercial building perimeter zones. (5) Supporting models and tools needed by industry for technology R&D and market transformation activities were validated. Attaining California’s clean energy goals require making a fundamental shift from today’s ad-hoc assemblages of static components to turnkey, intelligent, responsive, integrated building façade systems. These systems offered significant reductions in energy use, peak demand, and operating cost in California.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2022

Wall Upgrades for Energy Retrofits: A Techno-Economic Study

Homes built before 1992, when the U.S. Department of Energy’s (DOE) Building Energy Codes Program was established, represent approximately 68% of the residential building stock in the country. Up to 43% of these homes have little to no insulation in the walls and have very high air leakage rates of 10 or more air changes per hour at 50 pascals of pressure (ACH50). These issues can represent a substantial portion of unnecessary money spent on utility bills for homeowners, especially in the colder climates. There is a significant need for cost-effective, reliable retrofit methods for these homes that include air, moisture, and vapor controls which are considered best practices for high-performance new home construction. Well-tested and documented wall retrofit systems can help to achieve substantial energy savings and also improve durability, comfort, health, and resilience. In 2018, DOE’s Building Technologies Office awarded Pacific Northwest National Laboratory, Oak Ridge National Laboratory, and the University of Minnesota funding to complete a 3-year project to compare a range of residential wall retrofit systems that prioritized affordability, durability and energy savings potential. In addition to these core criteria, the ease-of-construction and the wide-scale applicability of the solutions also were considered. In this project, the research team identified, constructed, tested, simulated, and analyzed the feasibility and economics of 16 wall retrofit assemblies (14 test configurations and two baseline configurations).

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2022

Comparison of scattering ratio profiles retrieved from ALADIN/Aeolus and CALIOP/CALIPSO observations and preliminary estimates of cloud fraction profiles

The space-borne active sounders have been contributing invaluable vertically resolved information of atmospheric optical properties since the launch of Cloud-Aerosol Lidar and Infrared Pathfinder Satellite Observation (CALIPSO) in 2006. To build long-term records from space-borne lidars useful for climate studies, one has to understand the differences between successive space lidars operating at different wavelengths, flying on different orbits, and using different viewing geometries, receiving paths, and detectors. In this article, we compare the results of Atmospheric Laser Doppler INstrument (ALADIN) and Cloud-Aerosol Lidar with Orthogonal Polarization (CALIOP) lidars for the period from 28 June to 31 December 2019. First, we build a dataset of ALADIN–CALIOP collocated profiles (Δdist<1°; Δtime<6 h). Then we convert ALADIN's 355 nm particulate backscatter and extinction profiles into the scattering ratio vertical profiles SR(z) at 532 nm using molecular density profiles from Goddard Earth Observing System Data Assimilation System, version 5 (GEOS-5 DAS). And finally, we build the CALIOP and ALADIN globally gridded cloud fraction profiles CF(z) by applying the same cloud detection threshold to the SR(z) profiles of both lidars at the same spatial resolution. Before comparing the SR(z) and CF(z) profiles retrieved from the two analyzed lidar missions, we performed a numerical experiment to estimate the best achievable cloud detection agreement CDA norm (z) considering the differences between the instruments. We define CDA norm (z) in each latitude–altitude bin as the occurrence frequency of cloud layers detected by both lidars, divided by a cloud fraction value for the same latitude–altitude bin. We simulated the SR(z) and CF(z) profiles that would be observed by these two lidars if they were flying over the same atmosphere predicted by a global model. By analyzing these simulations, we show that the theoretical limit for CDA$_{norm}^{theor}$(z) for a combination of ALADIN and CALIOP instruments is equal to 0.81±0.07 at all altitudes. In other words, 19 % of the clouds cannot be detected simultaneously by two instruments due to said differences. The analyses of the actual observed CALIOP–ALADIN collocated dataset containing ~78 000 pairs of nighttime SR(z) profiles revealed the following points: (a) the values of SR(z) agree well up to ~3 km height. (b) The CF(z) profiles show agreement below ~3 km, where ~80 % of the clouds detected by CALIOP are detected by ALADIN as expected from the numerical experiment. (c) Above this height, the CDA$_{norm}^{obs}$(z) reduces to ~50 %. (d) On average, better sensitivity to lower clouds skews ALADIN's cloud peak height in pairs of ALADIN–CALIOP profiles by ~0.5±0.6 km downwards, but this effect does not alter the heights of polar stratospheric clouds and high tropical clouds thanks to their strong backscatter signals. (e) The temporal evolution of the observed CDA$_{norm}^{obs}$(z) does not reveal any statistically significant change during the considered period. This indicates that the instrument-related issues in ALADIN L0/L1 have been mitigated, at least down to the uncertainties of the following CDA$_{norm}^{obs}$(z) values: 68±12 %, 55±14 %, 34±14 %, 39±13 %, and 42±14 % estimated at 0.75, 2.25, 6.75, 8.75, and 10.25 km, respectively.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRSDate not supplied

Solar Cruiser TRAC Boom Development

The TRAC boom studied is a 30 meter, high-aspect ratio and highly nonlinear structural element planned for the NASA Solar Cruiser spacecraft where it would serve as the solar sail’s skeletal system designed to deploy the packaged reflective sail membrane as well as support operational sail tension loads for flight. Key design drivers for this boom are understanding and characterizing load-deformation behavior and quantifying buckling performance–under flight conditions (zero gravity). Knowing that the full-scale solar sail system cannot be feasibly tested for flight conditions under Earth’s gravity, a building block approach is methodically followed to verify finite element analysis (FEA) predictions and to evolve the finite element model (FEM) accordingly, beginning with the lowest level complexity to match global elastic behavior to the highest level complexity to bound theoretical responses under conservative expected flight conditions. Although the testing campaign for the TRAC boom has not developed yet to produce a full-scale test data point that loads the boom to failure, the fidelity of simulations have evolved with success and warrant further interest and work on this high-aspect ratio, but highly stable, deployable structural boom.

Lee Nguyen↗
DOE OSTI2024

Practical Hamiltonian learning with unitary dynamics and Gibbs states

We study the problem of learning the parameters for the Hamiltonian of a quantum many-body system, given limited access to the system. In this work, we build upon recent approaches to Hamiltonian learning via derivative estimation. We propose a protocol that improves the scaling dependence of prior works, particularly with respect to parameters relating to the structure of the Hamiltonian (e.g., its locality k). Furthermore, by deriving exact bounds on the performance of our protocol, we are able to provide a precise numerical prescription for theoretically optimal settings of hyperparameters in our learning protocol, such as the maximum evolution time (when learning with unitary dynamics) or minimum temperature (when learning with Gibbs states). Thanks to these improvements, our protocol has practical scaling for large problems: we demonstrate this with a numerical simulation of our protocol on an 80-qubit system.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

A synthetic building operation dataset

Abstract This paper presents a synthetic building operation dataset which includes HVAC, lighting, miscellaneous electric loads (MELs) system operating conditions, occupant counts, environmental parameters, end-use and whole-building energy consumptions at 10-minute intervals. The data is created with 1395 annual simulations using the U.S. DOE detailed medium-sized reference office building, and 30 years’ historical weather data in three typical climates including Miami, San Francisco, and Chicago. Three energy efficiency levels of the building and systems are considered. Assumptions regarding occupant movements, occupants’ diverse temperature preferences, lighting, and MELs are adopted to reflect realistic building operations. A semantic building metadata schema - BRICK, is used to store the building metadata. The dataset is saved in a 1.2 TB of compressed HDF5 file. This dataset can be used in various applications, including building energy and load shape benchmarking, energy model calibration, evaluation of occupant and weather variability and their influences on building performance, algorithm development and testing for thermal and energy load prediction, model predictive control, policy development for reinforcement learning based building controls.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2026

Algebraic Multigrid with Filtering: An Efficient Preconditioner for Interior Point Methods in Large-Scale Contact Mechanics Optimization

Large-scale contact mechanics simulations are crucial in many engineering fields such as structural design and manufacturing. In the frictionless case, contact can be modeled by minimizing an energy functional; however, these problems are often nonlinear, nonconvex, and increasingly difficult to solve as mesh resolution increases. In this work, we employ a Newton-based interior-point (IP) filter line-search method, an effective approach for large-scale constrained optimization. While this method converges rapidly, each iteration requires solving a large saddle-point linear system that becomes ill-conditioned as the optimization process converges, largely due to IP treatment of the contact constraints. Such ill-conditioning can hinder solver scalability and increase iteration counts with mesh refinement. Here, to address this, we introduce a novel preconditioner, algebraic multigrid with filtering (AMGF), tailored to the Schur complement of the saddle-point system. Building on the classical AMG solver, commonly used for elasticity, we augment it with a specialized subspace correction that filters near null space components introduced by contact interface constraints. Through theoretical analysis and numerical experiments on a range of linear and nonlinear contact problems, we demonstrate that the proposed solver achieves mesh independent convergence and maintains robustness against the ill-conditioning that notoriously plagues IP methods. These results indicate that AMGF makes contact mechanics simulations more tractable and broadens the applicability of Newton-based IP methods in challenging engineering scenarios. More broadly, AMGF is well suited for problems, optimization or otherwise, where solver performance is limited by a low-dimensional subspace, such as those arising from localized constraints, interface conditions, or model heterogeneities. This makes the method widely applicable beyond contact mechanics and constrained optimization.

Mathematics and Computing↗