Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Refining”

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 955 records · Page 53

Quantitative characterization of gradient microstructures: A study on friction stir spot processing of pure cobalt

Heterogeneous microstructures in polycrystalline metals can enhance the strength and ductility, outperforming homogeneous structures of similar composition. This study investigates deformed cobalt via friction stir spot processing with varying dwell times to uncover the effects of plastic deformation and heat generation on the formation of morphological, phase, and grain boundary character gradients. A new approach to quantify the morphological gradients in materials, which describes grain morphology in terms of density followed by parametric regression, enables direct quantification of processing depth and gradient sharpness. Results show that longer processing times increase the steepness of morphological gradients and reduce the deformation depth for friction stir spot processing with low plunge depths and high tool rotational speeds. The amount of retained FCC is increased in the shorter processing conditions, primarily due to refined grain size, increased defect content, and reduced heat generation. Crystallographic texture analysis of the HCP phase indicated a dominant B-fiber described by (0001) ∥ shear plane normal in the extreme processing conditions and the formation of a P-fiber, shear direction ∥ ⟨11$\bar2$0⟩ for intermediate dwell times. The texture of the FCC phase for low processing times was a C texture {100}⟨011⟩ where longer processing times were dominated by a [001] fiber texture with a main {110}⟨100⟩ orientation and emergence of a slight [111] fiber in the longest processing condition. The approaches outlined in this work give insight into quantifying gradients and improve the understanding of highly deformed cobalt.

36 MATERIALS SCIENCE↗

Gradient-informed Hamiltonian Monte Carlo for multicomponent CALPHAD model optimization and uncertainty quantification

CALPHAD model parameter optimization is inherently challenging due to non-smooth objective functions, high-dimensional parameter spaces, and the need for uncertainty quantification (UQ). Traditional weighted nonlinear least squares approaches are computationally efficient but local, whereas black-box global optimizers and ensemble Markov Chain Monte Carlo (MCMC) methods provide broader exploration at substantial computational cost. The objective of this work is to combine the global exploration capability of gradient-informed Hamiltonian Monte Carlo – specifically the No-U-Turn Sampler (NUTS) – with local deterministic refinement using BFGS to efficiently optimize multicomponent CALPHAD models with minimal manual intervention. Analytic gradients are computed via the Jansson derivative framework. The methodology is demonstrated on the Cr—Fe binary system and extended to the Cr—Fe—Ni ternary system with 32 degrees of freedom. For Cr—Fe, NUTS achieves comparable or superior optimality relative to ensemble MCMC while requiring over an order-of-magnitude fewer likelihood evaluations. Parameter uncertainties are quantified through NUTS sampling and propagated to thermodynamic observables using local expansion, demonstrating a novel modular approach that combines binary and ternary parameter subsets without requiring global relaxation. These results establish gradient-informed exploration as a scalable strategy for multicomponent CALPHAD optimization and provide a practical route towards efficient higher-order database development with quantified uncertainty.

36 MATERIALS SCIENCE↗

Process-Structure-Property Relations for As-Deposited Solid-State Additively Manufactured High-Strength Aluminum Alloy

Solid-state additive manufacturing methods provide innovative solutions to circumvent problems associated with materials susceptible to hot cracking by avoiding liquid solid phase transformations. In this work, the process parameter influence on microstructural evolution and mechanical response of a fully dense aluminum alloy 7050 (AA7050) component manufactured via a rapid, solid-state additive manufacturing process known as Additive Friction Stir Deposition (AFS-D) was quantified for the first time. Three sections (starting dwell, transient, crossover of roads) of the deposition that exhibit differing thermomechanical processing histories were evaluated for the resulting microstructure and mechanical response. The microstructural characterization was performed on the as-deposited AA7050 via Electron Backscatter Diffraction (EBSD), TEM, optical microscopy, and Scanning Electron Microscopy (SEM). The microstructural characterization revealed refined constituent particles and grains throughout the as-deposited AA7050 microstructure. Furthermore, quasi-static tensile experiments were conducted in both the build and transverse directions, in order to quantify the orientation influence on tensile properties of the as-deposited AA7050 build. Spatially dependent tensile properties were observed in the material due to heat input variation coarsening of secondary phases towards the initial layers of the AFS-D build. Post-mortem analysis revealed that voids nucleated and coalesced from the overgrowth of the strengthening precipitates present in the material, resulting in fracture.

36 MATERIALS SCIENCE↗

Design of heterogeneous structured Al alloys with wide processing window for laser-powder bed fusion additive manufacturing

Required microstructural attributes of an alloy vary with structural applications. The microstructural fine-tuning capability of laser-powder bed fusion (L-PBF) additive manufacturing (AM) enables application specific manufacture of the components. Such manufacture with L-PBF AM requires alloys that exhibit wide processing window and are amenable to multiple deformation mechanisms. However, high hot cracking susceptibility of Al alloys poses a barrier to such printability-performance synergy. In this work we show that an integration of, a) grain refinement through heterogeneous nucleation, and b) eutectic solidification, leads to crack-free parts at wide range of process parameters, microstructural heterogeneity, and hierarchy in the Al-Ni-Ti-Zr alloy. Such an integration targets hot cracking at multiple stages of solidification in L-PBF as opposed to the contemporary alloy design strategies that target hot-cracking at only specific stages of solidification. The Al-Ni-Ti-Zr alloy exhibits excellent printability and a high as-built tensile performance. Due to the wide processing window and amenability to multiple deformation mechanisms, the alloy microstructure and subsequently the performance, can be fine-tuned. Here, such strategy opens the gateway for application-specific manufacture of Al alloys with L-PBF AM and establishes a fundamental shift in current methodologies for design of these alloys for L-PBF AM.

36 MATERIALS SCIENCE↗

Fatigue performance of laser powder bed fusion hydride-dehydride Ti-6Al-4V powder

We report hydride-dehydride (HDH) Ti-6Al-4V alloy with particle size distribution of 50–120 µm is laser powder bed fusion (L-PBF) processed using optimum processing parameters and a near-fully dense structure with a density of 99.9 % is achieved. Microstructural observations and phase analyses indicate formation of columnar β grains with acicular α/α' phases in as-built condition. The roughness of the as-fabricated samples is significant with an average roughness of R a = 15.71 ± 3.96 µm and a root mean square roughness of R rms = 108.4 ± 24.9 µm, however, both values are reduced to R a = 0.19 ± 0.04 µm and R rms = 4.9 ± 0.6 µm after mechanical grinding. Mechanical tests are carried out on as-fabricated specimens followed by stress relief treatment. All samples are tested to failure in fatigue, under fully-reversed tension-compression conditions of R = –1. The as-built samples failed from the surface with crack initiation mainly at micro-notches, whereas after mechanically grinding, crack initiation changed to subsurface defects such as pores. Minimizing surface roughness by mechanically grinding eliminates surface micro-notches which improves fatigue strength in the high cycle fatigue region. Fatigue notch factor calculations showed that the effect of surface roughness was significantly lower when HDH powder is used compared to standard spherical powder. X-ray diffraction analysis revealed an in-plane compressive stress, micro-strain and grain refinement on the surface of the mechanically ground samples. Fractography observations (macroscale) revealed a fully brittle fracture in the first stage of crack growth with a transition to a dominantly ductile fracture in the third stage of crack growth. On the other hand, at the micro scale, even the brittle fracture regions showed evidence of ductile fracture within the α' martensite laths.

36 MATERIALS SCIENCE↗

The Trajectory of Additively Manufactured Titanium Alloys with Superior Mechanical Properties and Engineered Microstructures

Additive manufacturing (AM) plays a crucial role in the manufacturing and development of advanced materials including titanium alloys for engineering applications. During the last two decades, a great deal of effort has been devoted to print titanium parts with desired properties. However, due to presence of non-equilibrium/metastable phases, anisotropic microstructure, and various porosities in the printed titanium alloys, achieving an excellent combination of mechanical properties has been challenging. Here this review paper aims to provide a comprehensive overview of multiple post-process and in-process approaches to achieve enhanced strength-ductility combinations by microstructure engineering. Approaches such as addition of alloying elements/nucleants, multi-step post heat treatment, in-situ heat treatment, forced inter-pass cooling, in situ peening/rolling, hybridizing microstructure, high-intensity acoustic vibration, thermo-hydrogen refinement of microstructure, deliberately introducing lack of fusion defects with subsequent hot isostatic pressing are discussed in detail. Finally, we conclude by highlighting the challenges associated with each method and potential areas of improvement that can pave the way to develop reliable components with enhanced performance.

36 MATERIALS SCIENCE↗

Synergy of tensile strength and high cycle fatigue properties in a novel additively manufactured Al-Ni-Ti-Zr alloy with a heterogeneous microstructure

Alloy design strategies in additive manufacturing (AM) to achieve grain refinement and terminal eutectic solidification have been introduced to engineer Al alloys having microstructural hierarchy and heterogeneity. Such alloy design strategies enable crack-free builds with an expanded AM processing window and pushed the strength limit in Al alloys. However, fatigue performance of Al alloys made by AM is restricted by the presence of process induced defects and its stochasticity. In this work, tensile and high cycle fatigue (HCF) behavior of a novel Al-Ni-Ti-Zr alloy with a heterogeneous microstructure is studied in the as-built condition, supplemented by detailed microstructural and mechanical characterization. Excellent strength-ductility synergy of 342 MPa and 16% failure strain achieved in the alloy was associated with the microstructural attributes that pertain to the novel alloy. Additionally, the alloy showed excellent HCF performance with a fatigue endurance limit to ultimate tensile strength ratio of 0.29 in flexural fatigue mode. The study revealed the existence of multiple crack retardation mechanisms and favorable crack propagation pathways through the fine-grained regions which enabled good fatigue performance to the alloy. Further, a probabilistic model has been used to estimate the fatigue life of the alloy as a function of the stochastic microstructure by utilizing the statistical distribution of pores, solid-state inclusions, and grains in the AM Al alloy. Finally, the model parametric trends are consistent with the experimental observations.

36 MATERIALS SCIENCE↗

Additively manufactured Al-Ce-Ni-Mn alloy with improved elevated-temperature fatigue resistance

Here, the high cycle fatigue behavior of an additively manufactured (AM) Al-10.5Ce-3.1Ni-1.2Mn wt% alloy was evaluated at 350 °C. The measured fatigue strength of 50 MPa at this temperature is comparable to the most fatigue resistant wrought Al alloys (2618-T6 and 7079-T6) at 315 °C. Refinement of pore populations through advanced AM processing in the examined alloy led to oxide inclusions becoming dominant fatigue crack initiation sites. The distribution of microstructural features that determine crack initiation was measured by X-ray computed tomography and served as input to make fatigue strength model predictions with various methodologies. A Monte-Carlo based approach that was previously applied to predict the fatigue strength of cast aluminum alloys through pore size distributions yielded accurate predictions for the fatigue strength of this AM alloy. Given the tunability of defect distributions in AM alloys, the sensitivity of the Al-Ce-Ni-Mn alloy fatigue behavior to defect distributions and outstanding elevated-temperature fatigue resistance of this alloy; the results together suggest the possibility of a new generation of fatigue-resistant alloys produced by the additive manufacturing route.

36 MATERIALS SCIENCE↗

Grain structure and texture selection regimes in metal powder bed fusion

Additive manufacturing (AM) offers opportunities to produce complex part geometries not possible with conventional processing and in some cases even improve part performance. However, adoption has been slowed by difficulties assessing microstructure variability and there is no straightforward approach to relate processing to grain structure characteristics. In this study, datasets from AdditiveFOAM heat transport simulations of laser powder bed fusion (LPBF) are used to drive ExaCA simulations of grain structure. The GPU utilization of ExaCA and an algorithmic update for modeling melt pool overlap region solidification enabled rapid and parallel simulation across a wider range of process conditions than previously explored with cellular automata-based solidification models. A texture selection angle $θ_s$ is defined based on melt pool overlap geometry, and the range of $θ_s$ over which a commonly observed texture transition occurs in characterized AM builds was well-reproduced by ExaCA simulations over a wide range of melt pool shape, hatch spacing, and layer height. ExaCA simulations with 90 degree rotation of the scan direction on every other layer reproduced a number of trends from the AM literature including grain refinement, the dominance of layers with larger melt pools on the final grain structure, and the weakening or strengthening of texture depending on odd and even layer melt pool overlap geometry. EBSD data from a benchmark AM part is used to validate the simulated mechanism of a layer rotation-induced texture strengthening effect. Importantly, these results expand the understanding of the mechanisms for texture selection in alloys with cubic crystal symmetry and offer an approach to easily evaluate processing conditions. With this new understanding, these modeling tools will enable anticipation of previously unexpected variations in grain structure and target specific microstructures and properties.

36 MATERIALS SCIENCE↗

Uncovering grain and subgrain microstructure at the scale of additive manufacturing melt tracks with a scalable cellular automaton solidification model

Metal additive manufacturing, characterized by rapid solidification, yields refined grains with a distinctive cellular subgrain microstructure that plays a pivotal role in determining material properties. Due to the significant computational expense demanded to simulate the required physics with submicron spatial resolution, their numerical simulations have been limited to proof-of-concept studies to either 2D or small subregions of a melt pool. In this study, an open-source, scalable, solidification code, muMatScale, based on the cellular automaton method, has been developed to predict the grain and the underlying subgrain microstructure over an entire melt pool. The model incorporates flexible parallelization schemes, utilizing MPI and OpenMP GPU Offloading, in addition to appropriate multi-physics specific to non-equilibrium rapid solidification in AM. The impact of nucleation parameters on grain microstructures was investigated with a focus on grain size variations and morphology transitions. With selected nucleation parameters, the simulation predicted the grain size, subgrain morphology, crystallographic orientation, and microsegregation aligned with experimental measurements. The model demonstrates that epitaxial grain growth is a dominant factor at the melt pool boundary, influencing grain size variation under different grain sizes in the build plate while maintaining consistent primary dendrite arm spacing under identical thermal conditions. Here, the highly efficient numerical model enables large-scale simulations with a spatial resolution of 100 nm or less, unveiling unprecedented insights into thermal and solutal diffusion driven grain growth, and the subgrains with microsegregation within grains in 3D across scales. muMatScale will enable the linking of submicron length-scale microstructure to part-level material behavior by investigating fundamental solidification problems at the intercellular scale in many-track and many-layer builds.

36 MATERIALS SCIENCE↗

What explains the year-to-year variation in growing season timing of boreal black spruce forests?

Amplified climate warming in high latitudes is expected to affect growing season timing of the vast boreal biome. It is unclear whether the presence of permafrost (perennially frozen ground) might have an influence on changes in growing season timing. Here, this study examined how different environmental variables explained, either directly or indirectly, the variation in growing season timing of boreal forest stands with and without permafrost. We expected that environmental variables explaining the variation in growing season timing differed or had different explanatory power depending on permafrost presence or absence. The growing season was delineated from daily gross primary productivity (GPP) time series derived from 40 site-year data of net ecosystem carbon dioxide exchange measured with eddy covariance techniques over five black spruce (Picea mariana [Mill.])-dominated boreal forest stands in North America. In permafrost-free forest stands, a combination of start in canopy ‘green-up’ in spring and the timing of air and soil temperature increasing above freezing explained the start-of-season (SOS GPP ). Results from commonality analysis and structural equation modeling suggest that canopy ‘green-up’ and air temperature directly affected SOS GPP in permafrost-free forest stands. In addition, soil temperature acted as mediator for an indirect effect of air temperature on SOS GPP . In contrast, none of the environmental variables, or their combination, explained the variation in SOS GPP in forest stands with permafrost. The explanatory power of environmental variables was more consistent regarding the end-of-season (EOS GPP ). In both, forest stands with and without permafrost, EOS GPP was directly explained by mean soil water content in the fall and the first day of continuous snowpack formation. A better understanding how environmental variables control SOS GPP and EOS GPP in forest stands with and without permafrost will help to refine parameterizations of the boreal biome in Earth system models.

54 ENVIRONMENTAL SCIENCES↗

Remote sensing of seasonal variation of LAI and fAPAR in a deciduous broadleaf forest

Climate change is affecting the phenology of terrestrial ecosystems. In deciduous forests, phenology in leaf area index (LAI) is the primary driver of seasonal variation in the fraction of absorbed photosynthetically active radiation (fAPAR), which drives photosynthesis. Remote sensing has been widely used to estimate LAI and fAPAR. However, while many studies have examined both empirical and model-based relationships among LAI, fAPAR, and spectral vegetation indices (SVI) from remote sensing, few studies have systematically and empirically examined how relationships among these variables change over the growing season. In this study, we examine how and why seasonal-scale covariation differs among time series of remotely sensed SVIs and both LAI and fAPAR based on current understanding and theory. To do this we use newly available remote sensing data sets in combination with time series of in-situ measurements and a canopy radiative transfer model to analyze how seasonal variation in canopy and environmental conditions affect relationships among remotely sensed SVIs, LAI, and fAPAR at a temperate deciduous forest site in central Massachusetts. Our results show that accounting for seasonal variation in canopy shadowing, which is driven by variation in solar zenith angle, improved remote sensing-based estimates of LAI, fAPAR, and daily total APAR. Specifically, we show that the phenology of SVIs is strongly influenced by seasonal variation in near infrared (NIR) reflectance arising from systematic variation in the canopy shadow fraction that is independent of changes in LAI or fAPAR. Therefore, results of this work provide a refined basis for understanding how remote sensing can be used to monitor and model the phenology of LAI, fAPAR, APAR, and gross primary productivity in temperate deciduous forests.

54 ENVIRONMENTAL SCIENCES↗

Building a DFT+U machine learning interatomic potential for uranium dioxide

Despite uranium dioxide (UO 2 ) being a widely used nuclear fuel, fuel performance models rely extensively on empirical correlations of material behavior, leveraging the historical operating experience of UO 2 . Mechanistic models that consider an atomistic understanding of the processes governing fuel performance (such as fission gas release and creep) will enable a better description of fuel behavior under non-prototypical conditions such as in new reactor concepts or for modified UO 2 fuel compositions. To this end, molecular dynamics simulation is a powerful tool for rapidly predicting physical properties of proposed fuel candidates. However, the reliability of these simulations depends largely on the accuracy of the atomic forces. Traditionally, these forces are computed using either a classical force field (FF) or density functional theory (DFT). While DFT is relatively accurate, the computational cost is burdensome, especially for f-electron elements, such as actinides. By contrast, classical FFs are computationally efficient but are less accurate. For these reasons, we report a new accurate machine learning interatomic potential (MLIP) for UO 2 that provides high-fidelity reproduction of DFT forces at a similar low cost to classical FFs. We employ an active learning approach that autonomously augments the DFT training data set to iteratively refine the MLIP. To further improve the quality of our predictions, we utilize transfer learning to retrain our MLIP to higher-accuracy DFT+U data. We validate our MLIPs by comparing predicted physical properties (e.g., thermal expansion and elastic properties) with those from existing classical FFs and DFT/DFT+U calculations, as well as with experimental data when available.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Agentic traffic intelligence: Augmented human-in-the-loop scenario generation for microscopic traffic simulation

Traditional microscopic traffic simulation generation often relies on static datasets and manual design, limiting its ability to simulate complex conditions easily. This paper presents a novel framework, Agentic Traffic Intelligence, which combines human approval large language models (LLMs), the Real-Twin tool, and multi-agent systems to perform realistic microscopic traffic simulation scenario generation. The proposed framework incorporates human-in-the-loop (HIL) control, retrieval-augmented generation (RAG), and multi-agent control mechanisms. HIL mechanisms are used to guide multiple LLMs focused on attributes for microscopic simulation generation and to improve the interpretability and transparency of LLM execution for users. RAG enhances context extraction by dynamically integrating external knowledge sources for traffic scenario generation foundations. A multi-agent architecture with supervisory control coordinates the interaction of simulation components, including traffic simulators, control logic, and calibration tools. This enables the synthesis of simulation-ready scenarios that reflect dynamic demand profiles and behavior controls. Furthermore, the framework fuses multisource traffic data with unstructured context and supports iterative refinement through interactive user feedback. Validated through microscopic simulation using Simulation of Urban Mobility, the generated scenarios demonstrate high-fidelity network generation with inflow and turn movement and behavioral calibration, offering a robust and efficient tool for stress-testing and optimizing urban mobility systems.

Hierarchical multi-agent control↗

Digital-twin-based improvements to diagnosis, prognosis, strategy assessment, and discrepancy checking in a nearly autonomous management and control system

We report the Nearly Autonomous Management And Control System (NAMAC) is a comprehensive control system to assist plant operations by furnishing control recommendations to operators in a broad class of situations. This study refines a NAMAC system for making reasonable recommendations during complex loss of flow (LOF) scenarios with a validated EBR-II simulator, digital twins improved by machine-learning algorithms, a multi-attribute decision-making scheme, and a discrepancy checker for identifying unexpected recommendation effects. The performance of each NAMAC component is assessed, while the capability of NAMAC is demonstrated and evaluated in a class of LOF scenarios.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Uncertainty quantification of a physics-informed model based on sparse identification of a Thermal Energy Distribution System

Integrated energy systems (IES)s are crucial for enhancing the economy and efficiency of power generation sources (e.g., nuclear energy) necessary to unleash American energy dominance. These systems can be integrated with thermal energy storage (TES) and intermittent renewable energies to optimize overall energy use, peak-load regulation, and demand-side responses. However, the stabilization of energy generation, transport, and utilization introduces operational complexities that exceed the challenges of managing each sub-component individually. Currently, though IESs rely on human operators for efficiency and stability, reducing human error risk and enhancing performance through automation is highly desirable. Recent advances at Idaho National Laboratory have demonstrated successful control of the Thermal Energy Distributed System (TEDS). However, the automatic control system depends on a deterministic Sparse Identification of Nonlinear Dynamics with Control (SINDyC) model, which are trained based on simulation data from physics-based simulations. Because of uncertainties in physics-based simulation, SINDyC model results in large discrepancies against experimental data and cannot be reliably used in automatic control. In this paper, we present an innovative approach to address these discrepancies by quantifying uncertainties and developing a more robust model. We first generated trajectories by using first-principles physics codes to encapsulate the experiment. Next, we trained thousands of models by randomly sampling these trajectories. We then collapsed all those models into one probabilistic SINDyC by fitting a multivariate Gaussian distribution onto the resulting coefficient’s distribution. Despite its simplicity, our approach successfully produced 95% confidence intervals that captured the experimental trajectories. It even did so with a higher probability and better U-pooling score across six of the seven relevant quantities of interest (QoIs), as compared to other classical approaches. In conclusion, ongoing research is focusing on generating new experimental trajectories to validate this approach, and on employing Bayesian calibration to refine parametric uncertainties and guide future model development efforts.

22 - GENERAL STUDIES OF NUCLEAR REACTORS↗

Efficient Computation of Doppler-Broadened Elastic Scattering Kernel Moments Using Ladder-Operator Formulation

Anefficient routine for computing Legendre moments of the Doppler-broadened elastic scattering kernel, including resonance scattering effects, has been implemented in the ISOXML module of Griffin. Isotropic scattering in the center-of-mass system and the ideal gas model for target motion are assumed. A ladder-operator formulation is introduced to compute all Legendre moments from order 0 to N simultaneously, enabling near-linear scaling of computational cost with respect to the maximum Legendre order. A physics-based strategy for constructing outgoing energy grids has also been developed, in which a tailored base grid is combined with adaptive refinement to maintain accuracy while limiting the number of outgoing energy points. For energies between resonances, a constant cross-section model is employed to further reduce computational cost. In addition, a quantitative criterion is derived to determine isotope-wise cut-off incident energies based on a prescribed up-scattering probability coverage. For 238U, up to incident energies of approximately 75, 230, and 661 eV at 294, 900, and 2500 K (corresponding to a 2% up-scattering probability threshold), computation of P0 kernels requires 1–8 s and computation of P0–P5 kernels requires 0.4–4 min using a single thread, while maintaining 1–3% relative error in up-scattering probability. These results demonstrate that the proposed formulation enables accurate and computationally practical Doppler-broadened kernel generation for online multigroup cross-section production in Griffin.

Doppler-broadening↗

Modeling occupancy-driven building loads for large and diversified building stocks through the use of parametric schedules

Building energy modeling provides a fundamental tool to assess the potential for energy efficiency to contribute to reducing world energy consumption and global emissions. Occupancy-related operations are a key source of uncertainty for building energy analysis, particularly for aggregated building stocks. At a district or city level, it is critical to estimate aggregated power load profiles for sizing power grid infrastructure, power plant capacity allocation, and energy efficiency measures. The stochastic nature of behavior-related operations complicates the creation of models that accurately capture building load profiles for entire building stocks. This research introduces a new methodology called parametric schedules to model occupancy-driven schedules for large and diverse building stocks. In contrast to computationally expensive methodologies proposed in the literature, our work does not use a recursive time-consuming step. Occupancy is estimated by the extrapolation of operation times directly from metered electric consumption data; occupancy-related schedules are stochastically assigned to each building model, guaranteeing diversity of operation times in the stock. Our procedure has been tested on a large, diversified data-set of 25,000 commercial buildings in Los Angeles, California. It proved to be able to adequately represent the stochastic schedules diversity of the stock and to refine the stock calibration process by 1%. This innovative approach represents a useful asset for utility companies, grid operators, urban planners, and balancing authorities, seeking to improve building stock modeling and better estimate the impact of energy conservation measures. – This work is part of a larger stock modeling tool called ComStock, which is under development by NREL.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗