Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Scaling Studies”

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 775 records · Page 43
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Collisionless conduction in a high-beta plasma: a collision operator for whistler turbulence

The regulation of electron heat transport in high-β, weakly collisional, magnetized plasma is investigated. A temperature gradient oriented along a mean magnetic field can induce a kinetic heat-flux-driven whistler instability (HWI), which back-reacts on the transport by scattering electrons and impeding their flow. Previous analytical and numerical studies have shown that the heat flux for the saturated HWI scales as β$^{-1}_{e}$. These numerical studies, however, had limited scale separation and consequently large fluctuation amplitudes, which calls into question their relevance at astrophysical scales. To this end, we perform a series of particle-in-cell simulations of the HWI across a range of β$_e$ and temperature-gradient length scales under two different physical set-ups. The saturated heat flux in all of our simulations follows the expected β$^{-1}_{e}$ scaling, supporting the robustness of the result. We also use our simulation results to develop and implement several methods to construct an effective collision operator for whistler turbulence. The results point to an issue with the standard quasi-linear explanation of HWI saturation, which is analogous to the well-known 90° scattering problem in the cosmic-ray community. Despite this limitation, the methods developed here can serve as a blueprint for future work seeking to characterize the effective collisionality caused by kinetic instabilities.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Effects of Horizontal Resolution, Domain Size, Boundary Conditions, and Surface Heterogeneity on Coarse LES of a Convective Boundary Layer

Atmospheric properties in a convective boundary layer vary over a wide range of spatial scales and are commonly studied using large-eddy simulations (LES) in various configurations. We examine how the boundary layer depth and distribution of variability across scales are affected by LES grid spacing, domain size, inhomogeneity of surface properties, and external forcing. Two different setups of the Weather Research and Forecasting (WRF) model are analyzed. A semi-idealized configuration uses a periodic domain, flat surface, prescribed homogeneous surface heat fluxes, and horizontally uniform profiles of large-scale advective tendencies. A nested LES setup employs a larger domain and realistic initial and boundary conditions, including an interactive land surface model with representative topography and vegetation and soil types. Subdomains of identical size are analyzed for all simulations. Characteristic structure sizes are quantified using the variability scales L 50 and L 95 , defined such that features smaller than that contain 50% and 95% of the total variance, respectively. Progressive increase in L 50 from vertical velocity to temperature and moisture structures is systematically reproduced in all simulation configurations. This dependence of L 50 on the considered variable complicates the development of scale-aware parameterizations for models with grid spacing in the “terra incognita”. In simulations using a larger domain with heterogeneous surface properties, the development of internal mesoscale patterns significantly affects variance distributions inside analyzed subdomains. Sizes of boundary layer structures also strongly depend on the LES grid spacing and, in case of heterogeneous surface and topography, on location of the subdomain inside a larger computational domain.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

Robotics for Energy Infrastructure Construction on Earth and Beyond: Pilot Study Exploring Level 3-4 Semi-Autonomous Installation of a Reduced-Scale Photovoltaic Power Plant Prototype: Preprint

There is a growing need to build energy infrastructure at massive scale on both Earth and in space. On Earth, increasing electricity demand from applications such as AI data centers is driving rapid expansion of power systems, while in space, long-duration missions such as lunar operations require reliable and resilient energy infrastructure that can be deployed quickly. These needs outpace the capabilities of conventional construction methods. Robotic construction offers a path to faster, lower-cost, and more consistent deployment, yet most existing approaches focus on automating individual tasks rather than end-to-end assembly. The Industrialized Construction Innovation (ICI) team at the National Laboratory of the Rockies (NLR), in collaboration with the Colorado School of Mines, has been investigating robotics and automation strategies for rapid end-to-end energy infrastructure construction. This paper presents a pilot study demonstrating level 3-4 semi-autonomous, end-to-end robotic assembly of a reduced-scale photovoltaic power plant prototype under terrestrial conditions. The approach combines hardware redesign of photovoltaic components for robotic compatibility, discrete-event and physics-based robotic simulations in synthetic environments, and sim-to-real deployment to execute end-to-end robotic assembly of[NM1.1][NM1.2] tracking components, structural members, photovoltaic modules, and embedded electrical connections. The pilot identifies key limitations like tolerance management, perception and localization, end effector capability, fastening methods, and multi-physics simulation fidelity, that must be addressed to enable end-to-end robotic assembly of full-scale energy infrastructure. The study highlights how these insights inform scaling strategies for terrestrial energy infrastructure and adaptation to unstructured and extreme environments, including off-world scenarios such as robotic assembly of vertical tall lunar towers.

14 SOLAR ENERGY↗
DOE OSTI2024

Data for Pilot-Scale Processing of Miscanthus x giganteus for Recovery of Anthocyanins Integrated with Production of Microbial Lipids and Lignin-Rich Residue

Chemical-free hydrothermal pretreatment of Miscanthus x giganteus (Mxg) at the lab scale using high liquid-to-solid ratios resulted in the recovery of anthocyanins and enhanced enzymatic digestibility of residual biomass. In this study, the process is scaled up by using a continuous hydrothermal pretreatment reactor operated at a low liquid-to-solid ratio (50 % w/w solids) as an important step towards commercialization. Anthocyanin yield was 70 % w/w at the pilot scale (50 kg of Mxg), compared to the 94 % w/w yield achieved at the lab scale (0.5 g of Mxg). The pretreated biomass was subsequently refined mechanically using a disc mill to increase the accessibility of cellulose by cellulases. Enzymatic saccharification of the pretreated and disc-milled residue yielded 238 g/L sugar concentration by operating in fed-batch mode at 50 % w/v solids content. Two strains of Rhodosporidium toruloides were evaluated for converting the hydrolysate sugars into microbial lipids, and strain Y-6987 had the highest lipid titer (11.0 g/L). Further, the residue left after enzymatic saccharification was determined to be enriched 1.7-fold in the lignin content. This lignin-rich residue has value as a feedstock for the production of sustainable aviation fuel precursors and other high-value lignin-based chemicals. Hence the proposed biorefinery based on Mxg creates an opportunity for generating revenue from multiple high-value products. As the demand for biofuels and biobased products is rising, the biorefinery products from Mxg would create a niche in the industrial sector.

Conversion↗
DOE OSTI2021

First Simultaneous Acceleration of Multiple Charge States of Heavy Ion Beams in a Large-Scale Superconducting Linear Accelerator

Experimental studies of the simultaneous acceleration of three-charge-state 129 Xe 49+,50+,51+ beam from 17 to 180 MeV/nucleon in a superconducting linear accelerator are presented. The beam parameters for each individual- and multiple-charge-state beam were measured and compared with the particle tracking simulations. Detailed measurements were performed to characterize the multiple-charge-state beam’s recombination after a second-order achromat and isopath 180° bending system. As a result of the recombination of three charge states in the six-dimensional phase space, the xenon beam intensity was increased by 2.5-fold compared to the single-charge-state beam. Lastly, the results presented in the Letter fully validate the possibility to produce and utilize high-quality multiple-charge-state heavy-ion beams in a large-scale superconducting linac to increase the available beam power on an isotope production target.

43 PARTICLE ACCELERATORS↗
DOE OSTI2023

Community Centered Solar Development (CCSD) Case Study Interviews [Slides]

Large-scale solar (LSS, defined here as ground-mounted photovoltaic projects ≥1 MWDC) has grown rapidly in the U.S., accounting for nearly half of new electric generating capacity added to the U.S. grid in 2022. All sources of electricity bring positive and negative impacts to hosting communities and the rapid growth of LSS has increased the urgency to understand those impacts. Yet, information about the potential positive and negative impacts of LSS on host communities, and the factors or drivers leading to support or opposition to a project, is lacking. This information gap limits how project developers, municipalities, and local siting authorities can address community concerns and appropriately align proposed projects to best suit and benefit local communities. As part of Berkeley Lab’s Community-Centered Solar Development (CCSD) project, this research set out to explore deep insights and perceptions from LSS stakeholders that only qualitative data can provide to identify key factors driving project success or threatened failure. Case studies, such as those utilized in this research, are uniquely adept at capturing the subjective experience of individuals and at identifying variables, structures, and interactions between stakeholders. Our case studies included 54 semi-structured interviews across 7 different LSS sites, representing a diversity of geographies, project sizes (MW), site types (i.e., greenfield, agrivoltaic, and brownfield / contaminated sites), zoning jurisdiction types, and more (Table 1). In addition to local residents living in close proximity to these LSS sites, we interviewed other key stakeholders involved in the projects such as developers, decision-makers, utility representatives, landowners, and individuals from community-based organizations. The overarching aim of this case study research was two-fold: (1) to inform subsequent tasks in the CCSD research project (including an upcoming national survey of LSS neighbors), and (2) to provide insights into the following set of research questions: -What are the key positive and negative drivers leading to support and opposition to LSS projects? -To what extent do LSS projects exacerbate or mitigate perceived inequities and marginalization within hosting communities and how can those inequities be mitigated going forward? -What strategies can communities employ to align LSS development with local land-use plans and community needs and values? The research findings and next steps are described in this slide deck report.

14 SOLAR ENERGY↗
DOE OSTI2021

Multi-modal imaging of a single mouse brain over five orders of magnitude of resolution

Mammalian neurons operate at length scales spanning six orders of magnitude; they project millimeters to centimeters across brain regions, are composed of micrometer-scale-diameter myelinated axons, and ultimately form nanometer scale synapses. Capturing these anatomical features across that breadth of scale has required imaging samples with multiple independent imaging modalities. Translating between the different modalities, however, requires imaging the same brain with each. Here, we imaged the same postmortem mouse brain over five orders of spatial resolution using MRI, whole brain micrometer-scale synchrotron x-ray tomography ($\mu$CT), and large volume automated serial electron microscopy. Using this pipeline, we can track individual myelinated axons previously relegated to axon bundles in diffusion tensor MRI or arbitrarily trace neurons and their processes brain-wide and identify individual synapses on them. This pipeline provides both an unprecedented look across a single brain's multi-scaled organization as well as a vehicle for studying the brain's multi-scale pathologies.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2021

A Pore‐Scale Investigation of Mineral Precipitation Driven Diffusivity Change at the Column‐Scale

Abstract Mineral precipitation affects the pore structure and thus transport properties of porous media. In this study, we investigated the pore‐scale dynamics of precipitation in diffusion controlled systems and the resulting impacts on the effective diffusivity, using a micro‐continuum reactive transport model. Forty two‐dimensional pore structures representing both idealized and realistic geometries were simulated with consideration of different precipitation scenarios and rates. A homogeneous nucleation scenario reproduced patterns observed in previous experimental study showing mixing‐induced precipitation, and a surface growth scenario captured the pattern for mineral precipitation on a substrate with the same or similar mineral structures. In all cases, local precipitation resulted in the reduction in the average porosity of the domain ( Φ ) until the cessation of diffusive transport and the termination of precipitation. The minimum porosity reached was referred to as the critical porosity ( Φ c ). The effective diffusivity ( D eff ) decreased with Φ and dropped sharply to effectively zero, that is, the critical effective diffusivity ( D eff c ), as Φ c was reached. These pore‐scale dynamics can be captured by a revised relationship that explicitly considers the critical porosity and the corresponding effective diffusivity, and the pre‐exponential coefficient and the exponent of the relationship varied with initial pore structure and the precipitation kinetics. Overall, the homogeneous nucleation scenario results in systematically larger Φ c and coefficients that give rise to a sharper decrease in diffusivity as Φ c is approached, compared to the surface growth scenario. The revised relationship was also implemented at continuum scale and used to examine column scale diffusivity change and reactions.

58 GEOSCIENCES↗
DOE OSTI2025

Impact of melt viscosity on filler dispersion in elastomeric nanocomposites

Compounding of commercial nanocomposites usually involves the addition of viscosity enhancers such as binder resins in ink jet inks, and paints. Contrary to this, plasticizers such as process oils are added to reduce the melt viscosity and ease processability of reinforced elastomers. Nanofillers such as silica and carbon black are typically added to reinforce rubber and enhance performance of automotive tire treads. Filler dispersion has traditionally been qualitatively (indirectly) assessed by measuring the properties of reinforced elastomers. While dispersion can be quantified by examining filler agglomeration through surface roughness measurements and microscopy, the size-scale dependence for these hierarchical fillers has usually been ignored. Here, we have recently devised a method to quantify nano-scale dispersion of fillers using Ultra small-angle X-ray scattering (USAXS) techniques. This method is advantageous since it directly links the controllable processing/compounding parameters such as the mixing speed, mixer geometry, residence time (or mixing duration), melt density, flow gap distance, and melt viscosity to nano-scale dispersion. While our previous studies have explored the impact of different processing parameters, this study specifically investigates the impact of melt viscosity on nano-scale filler dispersion in elastomer compounds. Commercially available polybutadienes with different Mooney viscosities were used in conjunction with different grades and amounts of process oils to modify the melt viscosity.

Carbon black↗
DOE OSTI2023

An Investigation on the Pollen-Induced Soiling Losses in Utility-Scale PV Plants

Here in this study, the impact of pollen as a PV soiling agent is investigated. The performance data of five utility-scale PV plants in North Carolina, USA, was collected and analyzed using two soiling extraction methods. Satellite and environmental data, including pollen counts, cropland, and vegetation, was also collected and analyzed to identify impacts to soiling losses. During the spring peak pollen season, performance losses of >15% were observed at all five sites. Partial performance recoveries following the pollen season were slow, with lack of correlation with rainfall. This means that the statistical soiling estimation methods that assume abrupt performance recovery from rain are not appropriate for pollen-impacted solar sites. When manual cleanings were performed on site the performance recovery ranged from 5% to 11% indicating persistent soiling impacts are present in this region. The results of this work provide new insights into the phenomenon of pollen deposition on PV systems, demonstrating that 1) soiling can also affect systems located in rainy locations and 2) that its effects cannot be determined using the current estimation methodologies.

14 SOLAR ENERGY↗
DOE OSTI2021

Impacts of Pore Network-Scale Wettability Heterogeneity on Immiscible Fluid Displacement: A Micromodel Study

The mixed-wet nature of reservoir formations imposes a wide range of rock wettability from strong resident-fluid wetting to strong invading-fluid wetting. The characteristics of two-phase flow in porous media composed of mixed-wetting surfaces remain poorly understood. In this study, we investigated the displacement of resident ethylene glycol (EG) by hexane in two mixed-wet micromodels of identical 2.5-D geometry heterogeneity, with uniformly or heterogeneously distributed patches strongly wetting to hexane. These patches are mixed among pores with unaltered EG-wetting surfaces. Along with control tests in the originally EG-wet micromodel, here we show the classic fingering and transitions in flow regimes at $logCa$ (capillary number) from -7.2 to -3.9. Moreover, pore-scale distributions of wettability and their spatial correlation influence displacement efficiency. In the two mixed-wet micromodels, we found (a) an increase of steady-state hexane saturation at the end of experiments by up to 0.12 in the capillary fingering regime and a decrease of at most by 0.06 in the viscous fingering regime, compared to the EG-wet micromodel, and (b) dispersed and fragmented hexane distribution after displacement. Brine drainage during supercritical CO 2 (scCO 2 ) injections in these micromodels occurs with lower wettability contrasts, and under similar viscosity ratios and interfacial tensions resulted in higher displacement efficiency relative to displacement of EG by hexane. While mixed-wettability can enhance displacement efficiency compared to uniform wettability, the dynamics of immiscible fluids in strong mixed-wet reservoirs are expected to be less pronounced in contributing to the efficiency of geological CO 2 sequestration, oil recovery, and remediation of hydrocarbon-contaminated aquifers.

2.5-D micromodel↗
DOE OSTI2023

Mutually beneficial outcomes for hydropower expansion and environmental protection at a basin scale

Reshaping the scale of planning for hydropower development, from reaches to basin-scales, has been recommended as a more effective way to ameliorate the environmental impacts of hydropower. One approach is identifying mutually exclusive areas where development is precluded for conservation purposes and areas of low conservation value that present fewer barriers to development. This strategy, however, is less adoptable in developed countries where hydropower is already widespread and large-scale construction of new dams is unlikely. To broaden the adoption of basin-scale planning, alternative approaches and planning tools are needed for identifying mutually beneficial opportunities for simultaneous increases in hydropower capacity while improving environmental conditions. In this study, we present the Basin Scale Opportunity Assessment as a methodology to improve environmental conditions through either direct (on-site) or indirect (off-site) mitigation. We assess whether direct or indirect mitigation activities lead to optimal results in terms of added hydropower, environmental improvement, and monetary cost at a basin scale. We present two case studies for the Connecticut River and Roanoke River Basins, USA. Significant opportunities for expanding hydropower generating capacity are numerous in both basins. Results suggest that total hydropower capacity could be increased 4 to 7 % in the Roanoke and Connecticut Basins, respectively, without new dam construction and with net improvements in environmental conditions. We found that environmentally and economically optimal win-win strategies for increasing hydropower capacity and improving environmental conditions included improving environmental conditions in rivers downstream of existing dams. Off-site mitigation opportunities, such as dam removal and wetland mitigation, were identified as optimum solutions for achieving net environmental improvements only when they were associated with new hydropower construction. Our results demonstrate that opportunities to increase hydropower capacity and improve environmental conditions are expanded by viewing cumulative benefits at basin scales; however, increasing regulatory flexibility may be required to realize these opportunities.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Strong Lensing by Galaxies

Strong gravitational lensing at the galaxy scale is a valuable tool for various applications in astrophysics and cosmology. Some of the primary uses of galaxy-scale lensing are to study elliptical galaxies’ mass structure and evolution, constrain the stellar initial mass function, and measure cosmological parameters. Since the discovery of the first galaxy-scale lens in the 1980s, this field has made significant advancements in data quality and modeling techniques. In this review, we describe the most common methods for modeling lensing observables, especially imaging data, as they are the most accessible and informative source of lensing observables. We then summarize the primary findings from the literature on the astrophysical and cosmological applications of galaxy-scale lenses. We also discuss the current limitations of the data and methodologies and provide an outlook on the expected improvements in both areas in the near future.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2026

Fracture Analysis of Cohesive Zone Models for Modeling Residual Stress Induced Delamination in Composite Structures

A fracture study of coupon-scale composite cylinders with embedded defects was conducted with an objective to assess and validate a modeling approach using two available cohesive material models. The study included experimental and simulation evaluations of initiation of crack growth and progression. Interrupted thermal experiments used acoustic emissions monitoring to identify the onset of crack progression during each cooling interval and ultrasonic scanning provided images of defect growth. Verification, validation, and uncertainty quantification (VVUQ) processes were performed in the assessment of the simulation predicted temperature at which crack propagation begins (quantity of interest). The Sobol sensitivity analysis identified the hoop direction elastic modulus in the carbon fiber reinforced polymer (CFRP) plies as the most influential parameter for simulations using both cohesive models, accounting for at least 70% of the variation in the temperature at crack propagation. The UQ temperature range for the Tvergaard-Hutchinson model was higher (more conservative) than the experimental acoustic measurement indicators of crack progression, while the temperature range for the Thouless-Parmigiani model enveloped the experimental data points for the primary defect size of 0.75 x 1 in. The simulations could not capture the stable crack growth indicated in the experiments. This is likely due to the models’ inability to represent anisotropic fracture toughness attributed to the structure of the orthotropic fiber weave in a woven composite laminate.

42 ENGINEERING↗
DOE OSTI2021

Investigation of Small-Scale Electron Density Irregularities Observed by the Arase and Van Allen Probes Satellites Inside and Outside the Plasmasphere

In situ electron density profiles obtained from Arase in the night magnetic local time (MLT) sector and from RBSP-B covering all MLTs are used to study the small-scale density irregularities present in the plasmasphere and near the plasmapause. Electron density perturbations with amplitudes >10% from background density and with time-scales less than 30-min are investigated here as the small-scale density irregularities. The statistical survey of the density irregularities is carried out using nearly 2 years of density data obtained from RBSP-B and 4 months of data from Arase satellites. The results show that density irregularities are present globally at all MLT sectors and L-shells both inside and outside the plasmapause, with a higher occurrence at L > 4. The occurrence of density irregularities is found to be higher during disturbed geomagnetic and interplanetary conditions. The case studies presented here revealed: (1) The plasmaspheric density irregularities observed during both quiet and disturbed conditions are found to coexist with the hot plasma sheet population. (2) During quiet periods, the plasma waves in the whistler-mode frequency range are found to be modulated by the small-scale density irregularities, with density depletions coinciding well with the decrease in whistler intensity. Our observations suggest that different source mechanisms are responsible for the generation of density structures at different MLTs and geomagnetic conditions.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Scaling Arctic landscape and permafrost features improves active layer depth modeling

Tundra ecosystems in the Arctic store up to 40% of global below-ground organic carbon but are exposed to the fastest climate warming on Earth. However, accurately monitoring landscape changes in the Arctic is challenging due to the complex interactions among permafrost, micro-topography, climate, vegetation, and disturbance. This complexity results in high spatiotemporal variability in permafrost distribution and active layer depth (ALD). Moreover, these key tundra processes interact at different scales, and an observational mismatch can limit our understanding of intrinsic connections and dynamics between above and below-ground processes. Consequently, this could limit our ability to model and anticipate how ALD will respond to climate change and disturbances across tundra ecosystems. In this paper, we studied the fine-scale heterogeneity of ALD and its connections with land surface characteristics across spatial and spectral scales using a combination of ground, unoccupied aerial system, airborne, and satellite observations. We showed that airborne sensors such as AVIRIS-NG and medium-resolution satellite Earth observation systems like Sentinel-2 can capture the average ALD at the landscape scale. We found that the best observational scale for ALD modeling is heavily influenced by the vegetation and landform patterns occurring on the landscape. Landscapes characterized by small-scale permafrost features such as polygon tussock tundra require high-resolution observations to capture the intrinsic connections between permafrost and small-scale land surface and disturbance patterns. Conversely, in landscapes dominated by water tracks and shrubs, permafrost features manifest at a larger scale and our model results indicate the best performance at medium resolution (5 m), outperforming both higher (0.4 m) and lower resolution (10 m) models. This transcends our study to show that permafrost response to climate change may vary across dominant ecosystem types, driven by different above- and below-ground connections and the scales at which these connections are happening. We thus recommend tailoring observational scales based on landforms and characteristics for modeling permafrost distribution, thereby mitigating the influences of spatial-scale mismatches and improving the understanding of vegetation and permafrost changes for the Arctic region.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

A Case Study about Energy and Cost Impacts for Different Community Scenarios Using a Community-Scale Building Energy Modeling Tool

The United States building sector consumed approximately 75% of electricity in 2019. By implementing renewable energy technologies and control strategies into buildings, future buildings will serve as energy generators as well as consumers. To accommodate this transition, communications among buildings and between buildings and the grid could provide more possibilities to optimize the energy performance of buildings. This paper develops a community-scale building energy model tool and conducts a case study adopting behind-the-meter distributed energy resources, sharing energy in different buildings, and using different electricity tariff structures. Three scenarios are studied: (1) electricity only supplied by the grid, (2) photovoltaic (PV) panels installed on and available to some but not all buildings, and (3) a connected community. To consider the impacts of locations and energy tariffs, this paper selects four cities and three electricity tariffs to evaluate the energy and cost performances of these three scenarios. The results show that the PV panels in Scenario 2 reduce 25% to 33% of the community-level electricity consumption and 20% to 30% of the community-level electricity cost compared with Scenario 1 in all studied locations and energy tariffs. By considering power management in the connected community (Scenario 3), the electricity consumption and cost can be further reduced by 6% to 7% and 5% to 11%, respectively, compared with Scenario 2.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗
DOE OSTI2021

Modelica-based system modeling for studying control-related faults in chiller plants and boiler plants serving large office buildings

System modeling is critical when studying operation faults in chiller plants and boiler plants. However, current fault models have difficulties faithfully representing the operation of chiller plants and boiler plants under the effects of those faults, especially the control-related ones. In this study, we present a systematic method to develop high-fidelity models for approximating the behaviors of chiller plants and boiler plants under faulty conditions. Compared to existing ones, the resulting fault models have two advantages: first, they better characterize the dynamic patterns in the system operation. In those models, control architecture and control logic are faithfully implemented. Thus, they can be used to study control-related faults, such as incorrect staging control due to sensor bias and mistuned feedback control. Second, they are readily extensible and can support large-scale investigations to explore different faulty conditions/scenarios. Those models are established in a hierarchical structure while modules in each layer can be redeclared and parameterized at upper layers. In such a case, modifications to the models can be realized through model modifiers and the process can be easily streamlined with scripts. We applied the proposed models in a comprehensive fault impact evaluation of the 13 control-related faults of chiller and boiler plants. In this evaluation, the proposed plant model is coupled with the EnergyPlus thermal load model to study the impact of various faulty scenarios. Based on the evaluation results, we identified the faults that have the most significant impacts on the operation of the chiller and boiler plants, respectively. We also found that the relationship between the impacts of the studied faults and the severity level of the faults are highly non-linear

Huang, Sen↗