Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “dynamic load model”

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 487 records · Page 27
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Calculating adsorption isotherms using the two-phase thermodynamic method and molecular dynamics simulations

We describe the calculation of adsorption isotherms from molecular dynamics simulations based on the two-phase thermodynamic (2PT) model. The 2PT model developed for bulk fluid phases treats the gas-like components as hard spheres (HSs), which correctly recovers the limiting behaviors of unconfined fluids. We showed that this treatment, however, does not always lead to the correct zero-loading behavior in strongly confining systems. For methane adsorption into zeolite MFI, the HS reference state underestimates entropy by up to 20% at low loadings and leads to an order-of-magnitude increase in the adsorption onset pressure. To fix these issues, we propose the use of ideal adsorbed gas (IAG) as the gas reference model, the properties of which can be computed using the Widom insertion method on an empty adsorbent. We further describe three routes to compute adsorption isotherms from the Helmholtz free energy at different loadings. Comparing against established Monte Carlo (MC) methods, we found that the adsorption isotherms obtained using the IAG reference state agrees to within 40%, which corresponds to deviations of <5% in adsorption free energy. The isotherms calculated using the HS reference state underestimate the adsorption uptake at low to medium loadings in strongly confining systems, but its accuracy improves at higher loadings and as the pore size increases relative to the sorbate diameter. The methods described here provide an alternative approach for computing adsorption isotherms when MC simulations in an open ensemble are undesirable and enable a direct comparison of computed adsorption thermodynamics with experiments.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

Aeroelastic Analysis Using Deforming Cartesian Grids

Ongoing work in air-vehicle design illustrates the potential of advanced concepts to provide significant improvements in efficiency; but with their incorporation of lightweight flexible structures, such configurations may require active control systems to ensure reliability and safety. However, many contemporary analysis methods are inefficient for aeroelastic analysis and design of such configurations. This paper describes the development of a new approach that automates the geometry setup, mesh generation, and assembly of fluid–structural coupling interfaces to enable efficient aeroelastic and aeroservoelastic analysis of advanced concepts. The core elements for this approach are a cut-cell Cartesian grid-based computational fluid dynamics solver, a nonlinear beam element structural model, a conservative fluid–structural interface treatment, and the formulation and implementation of a new deforming grid capability within the cut-cell Cartesian grid solver. In this paper, emphasis is on this latter component with detailed description given of the mesh motion strategy, evaluation of fluxes and structural loads at the surface, and computation of geometrical properties such as cell volume, directed face areas, centroids, and motion-induced fluxes for deforming Cartesian grids required to advance the flow states. Aeroelastic simulations exercising the capability show favorable agreement with data and predictions in the literature for subsonic and supersonic applications.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Modeling Framework to Predict Melting Dynamics at Microstructural Defects in TNT-HMX High Explosive Composites

Many high explosive (HE) formulations are composite materials whose microstructure is understood to impact functional characteristics. Interfaces are known to mediate the formation of hot spots that control their safety and initiation. Here, to study such processes at molecular scales, we developed all-atom force fields (FFs) for Octol, a prototypical HE formulation comprised of TNT (2,4,6-trinitrotoluene) and HMX (octahydro-1,3,5,7-tetranitro-1,3,5,7-tetrazocine). We extended a FF for TNT and recasted it in a form that can be readily combined with a well-established FF for HMX. The resulting FF was extensively validated against experimental results and density functional theory calculations. We applied the new combined TNT-HMX FF to predict and rank surface and interface energies, which indicate that there is an energetic driver for coarsening of microstructural grains in TNT-HMX composites. Finally, we assess the impact of several microstructural environments on the dynamic melting of TNT crystal under ultrafast thermal loading. We find that both free surfaces and planar material interfaces are effective nucleation points for TNT melting. However, MD simulations show that TNT crystal is prone to superheating by at least 50 K on subnanosecond time scales and that the degree of superheating is inversely correlated with surface and interface energy. The modeling framework presented here will enable future studies on hot spot formation processes in accident scenarios that are governed by strong coupling between microstructural interfaces, material mechanics, momentum and energy transport, phase transitions, and chemistry.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Spatial description of dislocation nucleation in the shock response of single-crystal aluminum

Nonequilibrium molecular dynamics simulations of shock loaded single-crystal Al in the $\langle$100$\rangle$, $\langle$110$\rangle$, $\langle$111$\rangle$, and $\langle$123$\rangle$ orientations are conducted to study elastic and plastic shockwave formation and details associated with dislocation activity. A computer vision-based approach is implemented to capture the presence of dislocations and describe their spatial characteristics in the zone of nucleation behind the propagating shockwave. The methodology developed relies on the sequences of images extracted during shock loading that show dislocation activity within a cross section of the sample. Results reveal that the spacing between activated slip systems is orientation dependent and exhibits a modest reduction for the $\langle$100$\rangle$ and $\langle$111$\rangle$ orientations as shock pressure increases. Comparisons are made to existing theoretical models. Such relationships between shock pressure and dislocation activity, extracted from molecular dynamics simulations, can be used to inform higher length scale simulations or modeling of dislocation-based plasticity during shock.

36 MATERIALS SCIENCE↗
DOE OSTI2022

DC-Link Current Minimization Control for Current Source Converter-Based Solid-State Transformer

This article proposes a fast predictive control method and a small DC-link inductor to minimize the DC-link current in current-source converter (CSC)-based solid-state transformer. The DC-link current minimization can significantly reduce power loss and improve efficiency. The challenge of this problem is on improving both steady-state and dynamic performance. PI control methods and large DC-link inductors are conventionally used in the CSC but have limited dynamic performance. A model predictive control (MPC) method is proposed to achieve switching-cycle-level settling time, and the DC-link inductor is sized for 40% ripple to enable fast current change. Importantly, this article also proposes to minimize the DC-link current by varying the current even within a line cycle under single-phase load to improve the steady-state performance, in contrast with the reduction to a constant value in the literature. The proposed MPC features a constant switching frequency without weighting factors. The MPC does not have a high computational burden and is implemented in a regular digital controller for a prototype of soft-switching solid-state transformer (S4T) with reduced conduction loss. The effectiveness of the proposed method has been experimentally verified on the SiC S4T prototype during steady-state and dynamics under different multiport power flow conditions up to 2 kV peak. Here, the DC-link current in the experiments is close to the minimum current with a short zero-vector duration, which further verifies the performance of the proposed method.

14 SOLAR ENERGY↗
DOE OSTI2025

Spatial Impacts of Electric Vehicle Charging on Power Grid Stability: A Downtown Atlanta Case Study

The rapid increase in electric vehicle (EV) charging demand poses a potential risk to power grid stability, particularly as the spatial distribution of this demand remains underexplored. Existing research often focuses on technical optimization models while overlooking the geographic and human dynamics that affect energy consumption. This study addresses this gap by incorporating mobility data to estimate both building energy use and EV charging demand while also considering geographic factors for a better understanding of grid load. Using agent-based simulations and the Open-Source Distribution System Simulator, the study evaluates the effect of various EV penetration scenarios on grid voltage and unbalance. The results show that, although voltage remains within acceptable limits at lower EV penetration rates, significant voltage drop and unbalance occur as EV penetration exceeds 40%, particularly in residential areas with high charging demand. This study offers a framework for integrating spatial analysis and mobility data in power network simulations, providing insights for future EV infrastructure planning.

Pan, Melrose [ORNL] (ORCID:000000031627448X)↗
DOE OSTI2024

Assessing Nutrient and Carbon Responses to Agricultural Conservation Practices in Two Midwest Watersheds

The United States is undertaking efforts to transition to a low-carbon economy in response to the heightened impacts of climate change, which are largely attributed to decades of carbon-intensive development. A concerted effort by the energy sector to achieve net-zero carbon emissions is aimed at strategically reducing the carbon footprint of the energy supply chain, particularly in the production of feedstocks for biofuels. The production of biofuels relies heavily on land use and management practices in biomass cultivation. Among other factors, soil organic carbon (SOC) plays a crucial role in the biofuel carbon cycle, impacting land productivity, greenhouse gas emissions, and water quality. The process of carbon drawdown during plant growth and storage helps to mitigate the release of carbon dioxide (CO 2 ) into the atmosphere. Another key parameter is the release of nitrous oxide (N 2 O) — a potent greenhouse gas with a global warming potential 273 times that of CO 2 — from soil. While prioritizing low-carbon production, sustainable bioenergy also necessitates improved water quality and ecosystem services. Agricultural conservation practices can contribute to enhanced soil health and reduce nutrient and sediment loss into streams. However, studies that explore the broader impact of these practices on soil carbon storage and N 2 O emissions at a watershed scale are limited. Specifically, our understanding of how low-carbon feedstock production and management affect nutrient cycle dynamics is incomplete. This study assesses watershed responses to land management practices in two agriculturally dominant watersheds: the Raccoon River watershed and the Southfork of Iowa River watershed. The study focuses on carbon and nutrient dynamics, characterizing spatial and temporal variations in SOC, N 2 O, and nutrient loadings to examine the relationship between nutrient and carbon responses and agricultural conservation practices. The study employs the newly developed Soil Water Analysis Tool for Carbon (SWAT-C) model, calibrated using 20 years’ of climate and water monitoring data, to simulate and evaluate two agricultural conservation practices: no-till and crop residue harvest with cover crop planting. We compared the calibrated SWAT-C model with a historical baseline model, which allowed us to analyze various aspects of the watershed, including stream flow, suspended sediments, nitrogen, phosphorus, organic carbon, SOC at different depths, and N 2 O emissions from topsoil. Our SWAT-C modeling results indicate that, compared with results from the historical baseline model, no-till practices are positively correlated with SOC accumulation, reduced N 2 O emissions, and decreased soil loss. We also observed no change or slightly increased nitrogen and phosphorus loss to water bodies in the watersheds compared with the baseline. Conversely, crop residue harvest with cover crop planting improved water quality by reducing nutrient and soil losses but also increased N 2 O emissions. This SWAT-C modeling study establishes the groundwork for further smaller-scale and/or sub-basin-level assessments of nitrogen, phosphorus, and carbon cycles. It provides valuable science-based insights to inform policy decisions — particularly in the context of transitioning to biofuel feedstock production in these watersheds — by addressing concerns about nutrient exports to local water bodies and, ultimately, the Mississippi River.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Sensitivity Analysis of Modal Parameters of a Jacket Offshore Wind Turbine to Operational Conditions

Accurate estimation of offshore wind turbine (OWT) modal parameters has a prominent effect on the design loads, lifetime prediction, and dynamic response of the system. Modal parameters can vary during the operation of OWTs. This paper studies the variation and sensitivity analysis of an OWT’s modal parameters with respect to operational and environmental conditions. Three finite element models of a jacket-supported OWT at the Block Island Wind Farm are created within the OpenSees, SAP2000, and OpenFAST platforms and validated using experimental measurements. The OpenFAST model is used to simulate the modal parameters of the turbine under various wind speed, rotor speed, power, yaw angle, mean sea level, blade pitch angle, and soil spring values. The model-predicted modal parameters of the first fore–aft (FA) and side–side (SS) modes are compared to those identified from experimental measurements. Results from the simulations show that the first FA natural frequency and damping ratio mostly depend on the rotor speed and wind speed, respectively, while yaw angle and mean sea level do not have a visible effect. It is observed that there is about 8% stiffening in the first FA frequency and an aerodynamic damping of 7.5% during the operation of the OWT.

17 WIND ENERGY↗
DOE OSTI2022

Effects of Port Angle on Scavenging of an Opposed Piston Two-Stroke Engine

Opposed-piston 2-stroke (OP-2S) engines have the potential to achieve higher thermal efficiency than a typical diesel engine. However, the uniflow scavenging process is difficult to control over a wide range of speeds and loads. Scavenging performance is highly sensitive to pressure dynamics, port timings, and port design. This study proposes an analysis of the effects of port vane angle on the scavenging performance of an opposed-piston 2-stroke engine via simulation. A CFD model of a three-cylinder opposed-piston 2-stroke was developed and validated against experimental data collected by Achates Power Inc. One of the three cylinders was then isolated in a new model and simulated using cycle-averaged and cylinder-averaged initial/ boundary conditions. This isolated cylinder model was used to efficiently sweep port angles from 12 degrees to 29 degrees at different pressure ratios. Results indicate that scavenging performance is correlated with the bulk swirl ratio generated by these port angles. Scavenging performance is also sensitive to the pressure ratio across the engine. It was concluded that, for a given pressure ratio, the smallest port angle produces the best scavenging; however, previous work indicates lower bulk swirl ratio during combustion can have an impact on emissions formation.

33 ADVANCED PROPULSION SYSTEMS↗
DOE OSTI2024

Wind farm structural response and wake dynamics for an evolving stable boundary layer: computational and experimental comparisons

Abstract. The wind turbine design process requires performing thousands of simulations for a wide range of inflow and control conditions, which necessitates computationally efficient yet time-accurate models, especially when considering wind farm settings. To this end, FAST.Farm is a dynamic-wake-meandering-based mid-fidelity engineering tool developed by the National Renewable Energy Laboratory targeted at accurately and efficiently predicting wind turbine power production and structural loading in wind farm settings, including wake interactions between turbines. This work is an extension of a study that addressed constructing a diurnal cycle evolution based on experimental data (Quon, 2024). Here, this inflow is used to validate the turbine structural and wake-meandering response between experimental data, FAST.Farm simulation results, and high-fidelity large-eddy simulation results from the coupled Simulator fOr Wind Farm Applications (SOWFA)–OpenFAST tool. The validation occurs within the nocturnal stable boundary layer when corresponding meteorological and turbine data are available. To this end, we compared the load results from FAST.Farm and SOWFA–OpenFAST to multi-turbine measurements from a subset of a full-scale wind farm. Computational predictions of blade-root and tower-base bending loads are compared to 10 min statistics of strain gauge measurements during 3.5 h of the evolving stable boundary layer, generally with good agreement. This time period coincided with an active wake-steering campaign of an upstream turbine, resulting in time-varying yaw positions of all turbines. Wake meandering was also compared between the computational solutions, generally with excellent agreement. Simulations were based on a high-fidelity precursor constructed from inflow measurements and using state-of-the-art mesoscale-to-microscale coupling.

17 WIND ENERGY↗
DOE OSTI2026

Dynamic crushing of metal lattice metamaterials: Shock mode diagrams and transition to topology-independent compaction regime

Additively manufactured lattice metamaterials offer design versatility in strength and energy absorption and provide an additional degree of freedom through the selection of the lattice topology. Under quasistatic loading, the unit cell structure can strongly affect the stiffness, yield, and post-yield behavior, but whether and to what degree the effect of lattice topology persists into dynamic loading scenarios, up to the compaction shock regime, has not been established. LLNL’ s ALE3D hydrocode was used to perform a computational investigation of dynamic loading in multiple lattice types, including the gyroid, octet, Schwarz D, and rhombic dodecahedron, under impact velocities from 0.25 to 2.25 km/s. Shock Hugoniots for each lattice topology are generated and compared, suggesting that above a critical velocity, distinctions between architectures may not persevere and compacted lattices behave similarly. Here, to investigate the transition between topology-dependent quasistatic compression and the topology-independent regime above the critical velocity, a one-dimensional elastic-linear hardening plasticity-densified solid (E-LHP-DS) shock model for lattice materials was developed that relies upon confined compression to link the quasistatic and shock mechanics. Unlike similar works, the model does not assume rigid behavior prior to yield or locking behavior at densification, allowing a richer exploration of lattice mechanics. With only six parameters, the analytical model simultaneously fit quasistatic confined compression simulations for relative densities 0.1 $≤ \bar{ρ} ≤$ 0.9 and predicted dynamic compaction behavior to traverse several distinct shock modes, each defined by a critical impact speed (equivalently, critical stresses). Comparing the numerical results to the one-dimensional E-LHP-DS shock model predictions suggests that the topology-independence under strong shocks is linked to the onset of densification, which can be predicted based on quasistatic confined compression results.

Cellular material↗
DOE OSTI2021

First Annual Report on Development of Microwave Resonant Cavity Transducer for Fluid Flow Sensing: Development of Sensor Performance Model of Microwave Cavity Flow Meter for Advanced Reactor High Temperature Fluids

We are investigating a microwave cavity-based transducer for in-core high-temperature fluid flow sensing in molten salt cooled reactors (MSCR) and sodium fast reactors (SFR). This sensor is a hollow metallic cylindrical cavity, which can be fabricated from stainless steel, and as such is expected to be resilient to radiation, high temperature and corrosive environment of MSCR and SFR. The principle of sensing consists of making one wall of the cylindrical cavity flexible enough so that dynamic pressure, which is proportional to fluid velocity, will cause membrane deflection. Membrane deflection causes cavity volume change, which leads to a shift in the resonant frequency. Feasibility of the sensor was initially investigated with analytical derivations and with COMSOL RF Module computer simulations of resonant frequency spectral shift due to uniform load. We also investigated the mechanical integrity of the flowmeter’s membrane through analytical modelling and COMSOL Structural Mechanics Module computer simulations. Both the analytic model and COMSOL model showed that maximum stresses on the plate, which are at the radial boundary of the plate, are three orders of magnitude smaller than the material’s yield strength and ultimate tensile strength. This indicates that the sensor is at a low risk of mechanical failure. Using results from models, we have developed an initial design for a microwave K-band sensor. A cylindrical resonator prototype was fabricated from brass for the initial tests. The external dimensions of the cavity are matched to the flange of a standard WR-42 waveguide. Microwave field is coupled into the resonant cavity through a subwavelength-size aperture. A test article was developed consisting of a piping Tee with a bulkhead WR-42 microwave waveguide installed in a leak-proof assembly. A microwave waveguide circulator was installed in the setup to suppress the effect of reflections at the cavity entrance by increasing the isolation between the input and the output port. Preliminary spectral characterization of cavity spectral response was performed with a portable PXIe chassis microwave VNA with a custom GUI. Preliminary dry tests of the transducer response were conducted with a set of calibrated weights. Transducer frequency shift was shown to be monotonically increasing with increasing pressure. The next steps will involve investigation of the transducer performance for water flow sensing.

42 ENGINEERING↗
DOE OSTI2026

Electromagnetic Transient Modeling of Large Data Centers for Grid-Level Studies: Beta Release

The magnitude and complexity of electricity usage patterns from large data centers are having significant impacts on the operation and dynamics of the power grid; grid operators and planners require a range of specialized data center models to properly evaluate these impacts and specify technical solutions as needed. Towards addressing this need, Pacific Northwest National Laboratory (PNNL) has developed a library of electromagnetic transient (EMT) models for grid-level studies of data centers called the data center model library (DML). This report describes how the DML was created and how it may properly be used. This report details the DML’s beta release, completed in July 2026. This is a revision and expansion of the alpha release, which was made available in January 2026 The models present in the DML are generic models; subject matter expertise and additional technical data are needed to modify these models before they can represent any real data center. However, they will significantly reduce the level of effort required to develop site-specific models and can serve as a common starting point to guide industry towards a more refined consensus. Most of the models within DML are dedicated to representing the power electronics interfaces commonly used in modern data centers, such as double-conversion uninterruptible power supplies and single-phase power factor correction converters. These models are intended for use in grid-level studies and are a simplified aggregation of many small components. That said, background material on the physical and electrical design of large data centers is provided as companion material so that users can be aware of many of the details which have been omitted or streamlined as a matter of practical necessity. Additionally, guidance on the application of EMT analysis for data center interconnection studies is provided, which aids users in identifying when the DML is necessary and what sort of additional model development may be necessary for conducting real-world studies.

electromagnetic transients↗
DOE OSTI2026

Vanishing dynamic strength measured in a transient high-pressure phase of tin

This work reports two independent experimental estimations of the dynamic strength of the body-centered tetragonal (bct) 𝛾 phase of tin, which occurs at pressures above about 9 GPa and transforms promptly back to the ambient 𝛽 phase upon pressure release. Measuring strength in such a transient high-pressure phase is challenging. One measurement used free-surface Richtmyer-Meshkov instabilities generated with gas gun impact. Another set of measurements used ramp-release loading in a pulsed-power facility. Strength estimations came from comparison to simulations using a comprehensive multiphase modeling framework that includes equation of state, shear moduli, and strength separately for each phase, and treats mixed-phase regions. Both experimental methods found the 𝛾-phase strength to be very low, within experimental uncertainty of zero. The existing literature on other materials, by contrast, almost universally reports higher strength in transient high-pressure phases compared with ambient phases. Recent Molecular-Dynamics simulations on tin in the literature showed almost zero deviatoric strength during a deformation-induced bct → bct transformation in which one of the 𝛾-phase 𝑎 axes “flips” to the 𝑐 axis. This reorientation currently provides the most plausible explanation for the low observed strength in 𝛾-phase tin.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Modelling of edge plasma dynamics with active wall boundary conditions

A self-consistent 2D model is presented for transport in boundary plasma and plasma-facing material walls. Plasma dynamics in the domain is represented by a 2D collisional plasma fluid model in the edge-plasma code UEDGE, and transport of hydrogen and heat in the wall is represented by a system of reaction–diffusion equations in the 1D wall code FACE. To account for variation of parameters along the wall, in the coupled model multiple instances of the FACE code run in parallel. Here, the coupled model provides a tool for investigating a range of dynamic plasma–material interactions phenomena in 2D. For demonstration of its capability, one application of particular interest is the role of active wall in tokamak strike point sweeping proposed for mitigation of divertor heat loads. In the present study, the coupled calculations are applied to investigation of the impact of heat and hydrogen transport in the material wall on the divertor plasma and target heat load during sweeping of the target strike point for parameters of a high-power tokamak.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Interplay between liver and blood stages of Plasmodium infection dictates malaria severity via γδ T cells and IL-17-promoted stress erythropoiesis

Plasmodium replicates within the liver prior to reaching the bloodstream and infecting red blood cells. Because clinical manifestations of malaria only arise during the blood stage of infection, a perception exists that liver infection does not impact disease pathology. By developing a murine model where the liver and blood stages of infection are uncoupled, we showed that the integration of signals from both stages dictated mortality outcomes. This dichotomy relied on liver stage-dependent activation of Vγ4 + γδ T cells. Subsequent blood stage parasite loads dictated their cytokine profiles, where low parasite loads preferentially expanded IL-17-producing γδ T cells. IL-17 drove extra-medullary erythropoiesis and concomitant reticulocytosis, which protected mice from lethal experimental cerebral malaria (ECM). Adoptive transfer of erythroid precursors could rescue mice from ECM. Modeling of γδ T cell dynamics suggests that this protective mechanism may be key for the establishment of naturally acquired malaria immunity among frequently exposed individuals.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Interaction of Soil pH and Mineralogy Controls Soil Organic Matter Persistence through Changes in the Composition and Amount of Microbial Necromass

Microbial necromass–mineral associations are key to long-term soil organic matter (SOM) persistence. However, how soil pH and mineralogy interact to regulate SOM stability remains poorly understood. Here, we used artificial soils to test how three clay minerals (bentonite, kaolinite, and goethite), adjusted to four pH levels (5–8), affect microbial activity (respiration), microbial physiology (carbon use efficiency, CUE), microbial-derived residue material (necromass), and the formation and stability of mineral-associated organic matter (MAOM). Artificial soils were inoculated with a rhizosphere-derived microbial community cultured under the same pH conditions and on two representative simulated exudate types (organic acids and carbohydrates) and incubated for 6 weeks. In two complementary experiments, we added necromass from known microbial taxa to the same minerals across pH levels to isolate the role of necromass chemistry and loading. We found that soil pH shaped MAOM chemistry by altering microbial activity and necromass composition. In interaction with mineral type, pH also controlled MAOM thermal stability. Higher necromass loading weakened mineral-organic bonding, reducing MAOM stability, consistent with zonal mineral–organic interaction models. Our results demonstrate that microbial activity, rather than carbon use efficiency, better predicts MAOM formation and that pH-dependent necromass composition and loading govern MAOM persistence. These findings advance mechanistic understanding of SOM stabilization and have implications for predicting soil carbon dynamics under shifting environmental conditions.

carbon use efficiency↗
DOE OSTI2025

GridLAB-D Technical Support Document: Residential Equivalent Thermal Parameter Model

GridLAB-D is a power distribution systems simulation and analysis tool developed by Pacific Northwest National Laboratory. Although GridLAB-D and its underlying residential load model has been used in a variety of power systems analyses, comprehensive documentation of the derivation of the residential house model has not been formally published. The purpose of this technical support document is to serve as a comprehensive resource for the thermal dynamics modeling and house definition implemented in GridLAB-D. It covers model derivation, parameter definitions, and implementation steps in the codebase.

24 POWER TRANSMISSION AND DISTRIBUTION↗