SEARCH · Engineering Papers

Results for “power modeling and estimation”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,189 records · Page 66

All sources NASA NTRS DOE OSTI

NASA NTRS2012

Preliminary Assessment of Variable Speed Power Turbine Technology on Civil Tiltrotor Size and Performance

A Large Civil Tiltrotor (LCTR) conceptual design was developed as part of the NASA Heavy Lift Rotorcraft Systems Investigation in order to establish a consistent basis for evaluating the benefits of advanced technology for large tiltrotors. The concept has since evolved into the second-generation LCTR2, designed to carry 90 passengers for 1,000 nm at 300 knots, with vertical takeoff and landing capability. This paper performs a preliminary assessment of variable-speed power turbine technology on LCTR2 sizing, while maintaining the same, advanced technology engine core. Six concepts were studied; an advanced, single-speed engine with a conventional power turbine layout (Advanced Conventional Engine, or ACE) using a multi-speed (shifting) gearbox. There were five variable-speed power turbine (VSPT) engine concepts, comprising a matrix of either three or four turbine stages, and fixed or variable guide vanes; plus a minimum weight, twostage, fixed-geometry VSPT. The ACE is the lightest engine, but requires a multi-speed (shifting) gearbox to maximize its fuel efficiency, whereas the VSPT concepts use a lighter, fixed-ratio gearbox. The NASA Design and Analysis of Rotorcraft (NDARC) design code was used to study the trades between rotor and engine efficiency and weight. Rotor performance was determined by Comprehensive Analytical Model of Rotorcraft Aerodynamics and Dynamics (CAMRAD II), and engine performance was estimated with the Numerical Propulsion System Simulation (NPSS). Design trades for the ACE vs. VSPT are presented in terms of vehicle gross and empty weight, propulsion system weight and mission fuel burn for the civil mission. Because of its strong effect on gearbox weight and on both rotor and engine efficiency, rotor speed was chosen as the reference design variable for comparing design trades. Major study assumptions are presented and discussed. Impressive engine power-to-weight and fuel efficiency reduced vehicle sensitivity to propulsion system choice. The 10% weight penalty for multi-speed gearbox was more significant than most engine technology weight penalties to the vehicle design because drive system weight is more than two times engine weight. Based on study assumptions, fixed-geometry VSPT concept options performed better than their variable-geometry counterparts. Optimum design gross weights varied 1% or less and empty weights less than 2% among the concepts studied, while optimum fuel burns varied up to 5%. The outcome for some optimum configurations was so unexpected as to recommend a deeper look at the underlying technology assumptions.

Snyder, Christopher A.↗

DOE OSTI2024

Characterizing and improving the performance of molten-salt-steam heat exchangers in concentrating solar power plants

Shell-and-tube heat exchangers (HXs) for steam generation from molten salts in concentrating solar power (CSP) plants experience thermal fatigue due to significant temperature gradients and inherent transient operation. Molten salt-steam HX design lifespans exceed actual lifespans, and, as a consequence, designers overpredict plant profitability and operators neglect appropriate prescriptions to optimize these lifetimes. Here, this study refines HX lifespan estimates with data benchmarked against thermal-fluid mechanical modeling of stress and accumulated fatigue. Reduced-order thermal models of the molten salt-steam, shell-and-tube evaporator and superheater predict transient temperature profiles along the two HXs salt-steam flow paths. The modeled evaporator and superheater temperature profiles enable assessment of cyclic stresses within the HX tubesheets, where molten-salt HX failures are most common. Evaporator and superheater performance data from a current 110 MW elec commercial CSP plant provide a basis for validating the reduced-order HX models. HX life predictions derived from stochastic failure distributions serve as inputs for simulating and optimizing existing plant operations. The impact of the updated lifespans on overall plant revenue depends on operating scenarios. This study suggests that typical ramping rates for a CSP plant with a high-temperature Rankine cycle result in an evaporator and superheater life of approximately 10 and 25 years, respectively, compared to the design target of 30 years. Reduced HX lifespans decrease operational plant revenue on average by 4.6-5.1%. Furthermore, there may be as many as four HX replacements over the 30-year lifetime of the plant; and, purchase agreement loss due to failure to meet contractual production requirements can have ramifications that include the risk of bankruptcy.

14 SOLAR ENERGY↗

NASA NTRS2004

Improved Nuclear Reactor and Shield Mass Model for Space Applications

New technologies are being developed to explore the distant reaches of the solar system. Beyond Mars, solar energy is inadequate to power advanced scientific instruments. One technology that can meet the energy requirements is the space nuclear reactor. The nuclear reactor is used as a heat source for which a heat-to-electricity conversion system is needed. Examples of such conversion systems are the Brayton, Rankine, and Stirling cycles. Since launch cost is proportional to the amount of mass to lift, mass is always a concern in designing spacecraft. Estimations of system masses are an important part in determining the feasibility of a design. I worked under Michael Barrett in the Thermal Energy Conversion Branch of the Power & Electric Propulsion Division. An in-house Closed Cycle Engine Program (CCEP) is used for the design and performance analysis of closed-Brayton-cycle energy conversion systems for space applications. This program also calculates the system mass including the heat source. CCEP uses the subroutine RSMASS, which has been updated to RSMASS-D, to estimate the mass of the reactor. RSMASS was developed in 1986 at Sandia National Laboratories to quickly estimate the mass of multi-megawatt nuclear reactors for space applications. In response to an emphasis for lower power reactors, RSMASS-D was developed in 1997 and is based off of the SP-100 liquid metal cooled reactor. The subroutine calculates the mass of reactor components such as the safety systems, instrumentation and control, radiation shield, structure, reflector, and core. The major improvements in RSMASS-D are that it uses higher fidelity calculations, is easier to use, and automatically optimizes the systems mass. RSMASS-D is accurate within 15% of actual data while RSMASS is only accurate within 50%. My goal this summer was to learn FORTRAN 77 programming language and update the CCEP program with the RSMASS-D model.

Robb, Kevin↗

DOE OSTI2020

Planck 2018 results. V. CMB power spectra and likelihoods

We describe the legacy Planck cosmic microwave background (CMB) likelihoods derived from the 2018 data release. The overall approach is similar in spirit to the one retained for the 2013 and 2015 data release, with a hybrid method using different approximations at low ( ℓ < 30) and high ( ℓ ≥ 30) multipoles, implementing several methodological and data-analysis refinements compared to previous releases. With more realistic simulations, and better correction and modelling of systematic effects, we can now make full use of the CMB polarization observed in the High Frequency Instrument (HFI) channels. The low-multipole EE cross-spectra from the 100 GHz and 143 GHz data give a constraint on the ΛCDM reionization optical-depth parameter τ to better than 15% (in combination with the TT low- ℓ data and the high- ℓ temperature and polarization data), tightening constraints on all parameters with posterior distributions correlated with τ . We also update the weaker constraint on τ from the joint TEB likelihood using the Low Frequency Instrument (LFI) channels, which was used in 2015 as part of our baseline analysis. At higher multipoles, the CMB temperature spectrum and likelihood are very similar to previous releases. A better model of the temperature-to-polarization leakage and corrections for the effective calibrations of the polarization channels (i.e., the polarization efficiencies) allow us to make full use of polarization spectra, improving the ΛCDM constraints on the parameters θ MC , ω c , ω b , and H 0 by more than 30%, and n s by more than 20% compared to TT-only constraints. Extensive tests on the robustness of the modelling of the polarization data demonstrate good consistency, with some residual modelling uncertainties. At high multipoles, we are now limited mainly by the accuracy of the polarization efficiency modelling. Using our various tests, simulations, and comparison between different high-multipole likelihood implementations, we estimate the consistency of the results to be better than the 0.5 σ level on the ΛCDM parameters, as well as classical single-parameter extensions for the joint likelihood (to be compared to the 0.3 σ levels we achieved in 2015 for the temperature data alone on ΛCDM only). Minor curiosities already present in the previous releases remain, such as the differences between the best-fit ΛCDM parameters for the ℓ < 800 and ℓ > 800 ranges of the power spectrum, or the preference for more smoothing of the power-spectrum peaks than predicted in ΛCDM fits. These are shown to be driven by the temperature power spectrum and are not significantly modified by the inclusion of the polarization data. Overall, the legacy Planck CMB likelihoods provide a robust tool for constraining the cosmological model and represent a reference for future CMB observations.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2021

Performance Comparison of Clipping Detection Techniques in AC Power Time Series

In this research, a variety of methods were developed to detect clipping periods in AC power time series. AC power data streams associated with 36 unique systems across the United States were collected, and data points representing clipping periods were manually labeled by experts. Using this data set for training and validation, novel logic-based and machine learning (ML) approaches were developed to classify time series values as clipping or non-clipping. These approaches were compared to the RdTools method for detecting clipping periods. The logic-based and ML XGBoost approaches achieved F-scores of 85.0 and 77.6, respectively, when cross-validated against the manually labeled data, as compared to the current RdTools approach (F-score of 56.4), indicating a significant improvement at detecting clipping periods. Additionally, the effects of each clipping filter when evaluating system degradation rates were assessed, using 31 unique systems across the United States. Results indicate that estimated system degradation rate can vary based on the type of clipping filter used, by up to 0.6% degradation rate for some cases.

clipping↗

DOE OSTI2024

Geothermal Power Systems Analysis: Outcome of Industry Stakeholders Workshop: Preprint

Geothermal cost and performance evaluation implemented via technoeconomic assessment (TEA) modeling is critical for the Department of Energy (DOE) and other geothermal industry stakeholders in assessing the current state of geothermal technologies and to identify existing hurdles to commercially viable geothermal development. The Geothermal Electricity Technology Evaluation Model (GETEM) is a major TEA tool used in estimating the economic feasibility and levelized cost of energy (LCOE) of conventional hydrothermal systems and enhanced geothermal systems (EGS). Since 2021, GETEM has been transitioning from an intricate spreadsheet model to a user-friendly tool within the System Advisor Model (SAM) developed by the National Renewable Energy Laboratory (NREL). Apart from enabling an expanded visibility of the geothermal model among other renewable resources, having GETEM in SAM has the advantage of simulation automation, better usability, updates tracking, active user inputs/feedback, and extended financial modeling. GETEM is used in developing supply curves for the Annual Technology Baseline (ATB). The ATB data are inputs to the Renewable Energy Potential (reV) and the Regional Energy Deployment System (ReEDS) models. The geothermal module in NREL’s reV model assesses the geothermal energy potential in the conterminous United States by defining the geospatial intersection of geothermal resources with existing grid infrastructure within the constraint of land use characteristics. The ReEDS model is a capacity expansion model used for simulating the long-term build-out and operation of the US generation and transmission system based on current energy costs and policies. To ensure enhanced representation of current industry trends in our model transitions and development, we organized a two-day virtual workshop to elicit geothermal industry stakeholder input and recommendations on our current approaches and assumptions on technoeconomic, resource assessment, and deployment scenarios modeling of geothermal technologies. Participants included developers, operators, investors, regulatory agencies, system modelers, national laboratory researchers, consultants, and other stakeholders. In this workshop, we gained stakeholder insights on current geothermal plant performance (i.e., capacity factors), updated drilling costs and learning curves, and next generation technologies such as closed loop and superhot rock geothermal. Other outcomes from this workshop and its impact on future geothermal development feasibility, resource availability, and capacity expansion studies are compiled and discussed.

Annual Technology Baseline↗

NASA NTRS1997

Statistical Modeling of Landscape Pattern and the Effects of Coarse Spatial Resolution

Analysis of classified satellite imagery was conducted to characterize errors in estimates of area based on coarse resolution satellite imagery which are due to distortions in sizes of small fragments, and to explore the feasibility of correcting for these errors using a statistical modeling approach. Sizes of bodies of open water on European Remote-Sensing (ERS-1) synthetic aperture radar (SAR) and fire scars on Landsat Multispectral Scanner (MSS) imagery were measured using geographic information system (GIS) software. Statistical analysis of the smaller scars and ponds as observed with this imagery of relatively fine resolution demonstrated that the distribution of the sizes could be modeled by either of two types of statistical distributions-a power distribution related to fractal processes or a simple exponential distribution. A comparison of the size distribution of small burn scars as observed with Landsat with the distribution observed with the much coarser Advanced Very High Resolution Radiometer (AVHRR) imagery revealed effects of coarse resolution, each of which were associated with a particular size range, and whose net effect accounted for the bias in total burn area due to coarse spatial resolution. Size dependent area correction procedures are therefore feasible. A statistical modeling approach to area estimation, based models of the underlying and observed size distributions of fragments in the scene, may lead to a procedure which requires little or no ancillary data.

Hlavka, C. A.↗

DOE OSTI2024

Assessment of thermal and radiation induced creep in the dual cooled lead lithium blanket

The creep in recently designed blanket for the fusion nuclear science facility is described by decomposing it into thermal, irradiation and cavity swelling creep. Initial estimates of thermal creep using elastic high temperature rules showed that the recent heat transfer improvements at the first wall (FW) is feasible, and the FW can withstand heat flux of 0.35 MW/m 2 at 3000 hrs without rupture. A viscoplastic model that is based on Norton's power law and relies on Multiphysics coupling of solid mechanics and heat transfer modules is used to capture the inelastic deformation in the blanket. Results showed that the design peak heat flux of 0.25 MW/m 2 produced maximum thermal creep strain of 0.45 % and the relaxation of the thermal stress at the FW. Irradiation creep is prescribed to be proportional to the displacement damage dose and applied stress. The displacements from irradiation creep radially decrease from the FW to the back wall, but the maximum deformation is found at the back wall that is connected to the Helium manifold due to the high stresses at the region. In conclusion, cavity swelling creep is dominant at steady state and combined with radiation creep to produce displacement of about 8 mm at the FW.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2020

Simulated wind farm wake sensitivity to configuration choices in the Weather Research and Forecasting model version 3.8.1

Wakes from wind farms can extend over 50 km downwind in stably stratified conditions. These wakes can undermine power production at downwind turbines, adversely impacting revenue. As such, wind farm wake impacts must be considered in wind resource assessments, especially in regions of dense wind farm development. The open-source Weather Research and Forecasting (WRF) numerical weather prediction model includes a wind farm parameterization to estimate wind farm wake effects, but model configuration choices can influence the resulting predictions of wind farm wakes. These choices include vertical resolution, horizontal resolution, and whether or not to include the addition of turbulent kinetic energy generated by the rotating wind turbines. Despite the sensitivity to model configuration, no clear guidance currently exists for these options. Here we compare simulated wind farm wakes produced by varying model configurations with meteorological observations near a land-based wind farm in flat terrain over several diurnal cycles. A WRF configuration comprised of horizontal resolutions of 3 km or 1 km paired with a vertical resolution of 10 m provides the most accurate representation of wind farm wake effects, such as the correct surface warming and elevated wind speed deficit. The inclusion of turbine-generated turbulence is also critical to produce accurate surface warming and should not be omitted.

17 WIND ENERGY↗

NASA NTRS1975

Power spectra of geoid undulations

Data from spacecraft altimeters are expected to contribute to an improved determination of the marine geoid. To better define altimeter system design requirements for geoid recovery, amplitudes of geoid undulations at short wavelengths were examined. Models of detailed geoids in selected areas around the earth, developed from a combination of satellite derived spherical harmonics and 1 deg-by-1 deg area mean free-air gravity anomalies, were subjected to a spectral analysis. The resulting undulation power spectra were compared to existing estimates for the magnitude of geoid undulations at short wavelengths. The undulation spectra were found to be consistent with Kaula's rule of thumb, following an inverse third power relationship with spatial frequency for wavelengths at least as small as 300 km. The requirements imposed by this relationship on altimeter accuracy, data rate, and horizontal resolution to meet the goal of a detailed geoid description are discussed.

Brown, R. D.↗

NASA NTRS1990

A novel multistage estimation of signal parameters

A multistage estimation scheme is presented for estimating the parameters of a received carrier signal possibly phase-modulated by unknown data and experiencing very high Doppler, Doppler rate, etc. Such a situation arises, for example, in the case of the Global Positioning Systems (GPS). In the proposed scheme, the first-stage estimator operates as a coarse estimator of the frequency and its derivatives, resulting in higher rms estimation errors but with a relatively small probability of the frequency estimation error exceeding one-half of the sampling frequency (an event termed cycle slip). The second stage of the estimator operates on the error signal available from the first stage, refining the overall estimates, and in the process also reduces the number of cycle slips. The first-stage algorithm is a modified least-squares algorithm operating on the differential signal model and referred to as differential least squares (DLS). The second-stage algorithm is an extended Kalman filter, which yields the estimate of the phase as well as refining the frequency estimate. A major advantage of the is a reduction in the threshold for the received carrier power-to-noise power spectral density ratio (CNR) as compared with the threshold achievable by either of the algorithms alone.

Kumar, Rajendra↗

NASA NTRS2002

Power System for Venus Surface Exploration

A radioisotope power and cooling system is designed to provide electrical power for a probe operating on the surface of Venus. Most foreseeable electronics devices and sensors cannot operate at the 450 C ambient surface temperature of Venus. Because the mission duration is substantially long and the use of thermal mass to maintain an operable temperature range is likely impractical, some type of active refrigeration may be required to keep electronic components at a temperature below ambient. The fundamental cooling parameters are the cold sink temperature, the hot sink temperature, and the amount of heat to be removed. In this instance, it is anticipated that electronics would have a nominal operating temperature of 300 C. Due to the highly thermal convective nature of the high-density (90 bar CO2) atmosphere, the hot sink temperature was assumed to be 50 C, which provided a 500 C temperature of the cooler's heat rejecter to the ambient atmosphere. The majority of the heat load on the cooler is from the high temperature ambient surface environment on Venus, with a small contribution of heat generation from electronics and sensors. Both thermoelectric (RTG) and dynamic power conversion systems were analyzed, based on use of a standard isotope (General-purpose heat source, or GPHS) brick. For the radioisotope Stirling power converter configuration designed, the Sage model predicts a thermodynamic power output capacity of 478.1 watts, which slightly exceeds the required 469.1 watts. The hot sink temperature is 1200 C, and the cold sink temperature is 500 C. The required heat input is 1740 watts. This gives a thermodynamic efficiency of 27.48 %. It is estimated that the mechanical efficiency of the power converter design is on the order of 85 %, based on experimental measurements taken from 500-watt power class, laboratory-tested Stirling engines. The overall efficiency is calculated to be 23.36 %. The mass of the power converter is estimated at approximately 21.6 kg. Additional information is included in the original extended abstract.

Landis, Geoffrey A.↗

DOE OSTI2025

Impact of Geothermal District Heating System on Flexibility of Microgrid in Tuttle, Oklahoma

Flexibility is the capability of the power grid to maintain a balance between electricity generation and variable demand. This study presents preliminary results evaluating the impact of geothermal district heating systems on the flexibility of microgrid in Tuttle, Oklahoma. Heating demand profiles were modeled using EnergyPlus for the district that includes two schools and 250 single-family houses. Then, geothermal energy production was modeled using GEOPHIRES to estimate how much heating demand in the district can be supplied by five different geothermal system scenarios. The results indicated that geothermal energy production varied depending on the resource temperature at different depths, system configurations, and flow rates. For the grid flexibility analysis, electricity consumptions in the five geothermal systems were estimated for pump operations to circulate water from the wells to radiators, while electricity consumption by air-source heat pump in the base case was estimated to supply the same heating load. Electricity consumption in the geothermal systems was significantly lower than those in base cases. The electricity saved by the geothermal system was then incorporated into the microgrid electrical load profiles where variable renewable electricity generation is significantly high. The results visually showed that geothermal district heating system can improve grid flexibility as a baseload during the winter season. The results also highlighted potential opportunities to save energy costs that will be further analyzed in future study.

Cambium↗

DOE OSTI2021

Estimating Spatial Distribution Impacts of Rooftops Solar PV on Dynamic Hosting Capacity Evaluation for a Real Distribution Feeder

This paper models the sensitivity of dynamic hosting capacity analysis to the spatial distribution of distributed PV scenarios, random, close and far from substation deployments, at various penetration levels on a real distribution feeder, in a quasi-static time series (QSTS) power flow simulation using high resolution time-aware metrics. These spatial distribution scenarios were chosen to capture a wide range of potential operation impacts of these deployments using a set of thermaland voltage-based metrics such as instantaneous violations and moving averages for both thermal and voltage constraints. This study uses actual load and PV data up-scaled from 1-hour time step to 1-minute resolution to fully characterize the interaction between the daily changes in load and PV output, and their impacts on distribution system operations.

distribution system↗

DOE OSTI2025

Modeling and Optimization of a Rotating Packed Bed Contactor with a Tetraamine-Appended Metal−Organic Framework for CO 2 Capture

A potential contactor technology for sorbent-based CO 2 capture is the rotating packed bed that contains separate sections for continuous adsorption and desorption. A heat exchanger can be embedded to remove heat in the adsorption section and add heat in the desorption section. In this work, we develop a two-dimensional (2D) model of a rotating packed bed for use in CO 2 capture applications. Mass and energy balances for the model are developed based on a Ljungström-type air preheater, which accounts for the counter-current axial flow of gas phases in separate sections of the bed and the rotation of a solid sorbent, which cycles between adsorption and desorption sections. The sorbent used for this analysis is the tetraamine-appended metal−organic framework Mg 2 (dobpdc)(3−4− 3), chosen for its stability and affinity for CO 2 at low partial pressures, such as those from a natural gas power plant source. An optimization problem is solved that considers the trade-off between maximizing the productivity of the bed and minimizing energy consumption. Maximum productivity and minimum energy are found to be 8.53 kg/h/m 3 and 3.84 MJ/kg, respectively, when these objectives are optimized independently. It is observed that the flue gas pressure and bed rotational speed are the desired operating variables to vary for model-based design of experiments to reduce uncertainty in parameter estimation, as these two variables yielded the most information content based on the Fisher information matrix.

20 FOSSIL-FUELED POWER PLANTS↗

DOE OSTI2025

Geothermal Power Systems Analysis: Outcome of Industry Stakeholders Workshop

Geothermal cost and performance evaluation implemented via techno-economic assessment (TEA) modeling is critical for the U.S. Department of Energy (DOE) and other geothermal industry stakeholders in assessing the current state of geothermal technologies and to identify existing hurdles to commercially viable geothermal development. The Geothermal Electricity Technology Evaluation Model (GETEM) is a major TEA tool used in estimating the economic feasibility and levelized cost of energy (LCOE) of conventional hydrothermal systems and enhanced geothermal systems (EGS). Since 2021, GETEM has been transitioning from an intricate spreadsheet model to a user-friendly tool within the System Advisor Model (SAM) developed by the National Renewable Energy Laboratory (NREL). Apart from enabling an expanded visibility of the geothermal model among other renewable resources, having GETEM in SAM has the advantage of simulation automation, better usability, updates tracking, active user inputs/feedback, and extended financial modeling. GETEM is used in developing supply curves for NREL's Annual Technology Baseline (ATB), which provides inputs to the Renewable Energy Potential (reV) and the Regional Energy Deployment System (ReEDS) models. The geothermal module in NREL's reV model assesses the geothermal energy potential in the conterminous United States by defining the geospatial intersection of geothermal resources with existing grid infrastructure within the constraint of land use characteristics. The ReEDS model is a capacity expansion model used for simulating the long-term build-out and operation of the U.S. generation and transmission system based on current energy costs and policies. To ensure enhanced representation of current industry trends in our model transitions and development, we organized a two-day virtual workshop to elicit geothermal industry stakeholder input and recommendations on our current approaches and assumptions on techno-economic, resource assessment, and deployment scenarios modeling of geothermal technologies. Participants included developers, operators, investors, regulatory agencies, system modelers, national laboratory researchers, consultants, and other stakeholders. In this workshop, we gained stakeholder insights on current geothermal plant performance (i.e., capacity factors), updated drilling costs and learning curves, and next-generation technologies such as closed-loop and superhot rock geothermal. Other outcomes from this workshop and its impact on future geothermal development feasibility, resource availability, and capacity expansion studies are compiled and discussed.

annual technology baseline↗

DOE OSTI2026

Hierarchical semi-Markov models with duration-aware dynamics for activity sequences

Residential electricity demand at granular scales is driven by what people do and for how long. Accurately forecasting this demand for applications like microgrid management and demand response therefore requires generative models for activities that can produce realistic daily activity sequences, capturing both the timing and duration of human behavior. This paper develops a generative model of human activity sequences using nationally representative time-use diaries at a 10-min resolution. We use this model to quantify which demographic factors are most critical for improving predictive performance. We propose a hierarchical semi-Markov framework that addresses two key modeling challenges. First, a time-inhomogeneous Markov router learns the patterns of “which activity comes next.” Second, a semi-Markov hazard component explicitly models activity durations, capturing “how long” activities realistically last. To ensure statistical stability when data are sparse, the model pools information across related demographic groups and time blocks. The entire framework is trained and evaluated using survey design weights to ensure our findings are representative of the U.S. population. On a held-out test set, we demonstrate that explicitly modeling durations with the hazard component provides a substantial and statistically significant improvement over purely Markovian models. Furthermore, our analysis reveals a clear hierarchy of demographic factors: Sex, Day-Type, and Household Size provide the largest predictive gains, while Region and Season, though important for energy calculations, contribute little to predicting the activity sequence itself. The result is an interpretable and robust generator of synthetic activity traces, providing a high-fidelity foundation for downstream energy systems modeling.

24 POWER TRANSMISSION AND DISTRIBUTION↗

DOE OSTI2023

The influence of physical and algorithmic factors on simulated far-field waveforms and source–time functions of underground explosions using unsupervised machine learning

SUMMARY Characterizing explosion sources and differentiating between earthquake and underground explosions using distributed seismic networks becomes non-trivial when explosions are detonated in cavities or heterogeneous ground material. Moreover, there is little understanding of how changes in subsurface physical properties affect the far-field waveforms we record and use to infer information about the source. Simulations of underground explosions and the resultant ground motions can be a powerful tool to systematically explore how different subsurface properties affect far-field waveform features, but there are added variables that arise from how we choose to model the explosions that can confound interpretation. To assess how both subsurface properties and algorithmic choices affect the seismic wavefield and the estimated source functions, we ran a series of 2-D axisymmetric non-linear numerical explosion experiments and wave propagation simulations that explore a wide array of parameters. We then inverted the synthetic far-field waveform data using a linear inversion scheme to estimate source–time functions (STFs) for each simulation case. We applied principal component analysis (PCA), an unsupervised machine learning method, to both the far-field waveforms and STFs to identify the most important factors that control variance in the waveform data and differences between cases. For the far-field waveforms, the largest variance occurs in the shallower radial receiver channels in the 0–50 Hz frequency band. For the STFs, both peak amplitude and rise times across different frequencies contribute to the variance. We find that the ground equation of state (i.e. lithology and rheology) and the explosion emplacement conditions (i.e. tamped versus cavity) have the greatest effect on the variance of the far-field waveforms and STFs, with the ground yield strength and fracture pressure being secondary factors. Differences in the PCA results between the far-field waveforms and STFs could possibly be due to near-field non-linearities of the source that are not accounted for in the estimation of STFs and could be associated with yield strength, fracture pressure, cavity radius and cavity shape parameters. Other algorithmic parameters are found to be less important and cause less variance in both the far-field waveforms and STFs, meaning algorithmic choices in how we model explosions are less important, which is encouraging for the further use of explosion simulations to study how physical Earth properties affect seismic waveform features and estimated STFs.

58 GEOSCIENCES↗