Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “PERIODICITY”

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 631 records · Page 35
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Delay-locked-loop timing error mitigation

Systems, methods, and circuits provide delay-locked loop (DLL) timing error mitigation. A DLL false-lock detection system can include DLL circuitry configured to receive a reference clock signal having a time period. The system can include shift register circuitry and latched comparison circuitry which can determine a time period of a locked condition of the DLL delay line with respect to the reference clock signal time period. The system can determine whether the system is correctly locked to the base time period or incorrectly locked to a multiple of the base time period. A further system can operate to cause a phase detector circuitry in a DLL to ignore the first edge of a reference clock signal presented to the phase detector circuitry and thereby avoid stuck-lock conditions.

Filippini, Leo↗
DOE OSTI2021

Intact archeological human bones and age at death studied with transmission x–ray diffraction and small angle x–ray scattering

In this work, high-energy, wide-angle x-ray scattering (WAXS, x-ray diffraction) and small-angle x-ray scattering (SAXS) were used to study intact human second metacarpal bones (mc2) from two UK archeological sites. A novel method correcting for irregular mass distribution was applied in these transmission geometry experiments done at beamline 1-ID of the Advanced Photon Source. The authors asked whether there were age-at-death-related changes in carbonated apatite (cAp) lattice parameters and whether SAXS could detect collagen D-period peaks in the archeological mc2. For each of the two sites, Ancaster and Wharram Percy in England, six female mc2s were studied; for each site, two were from each of three age-at-death cohorts (young, 18-29 years; middle, 30-49 years; old >= 50 years) along with a modern control mc2. The Rietveld method was applied to the WAXS patterns to provide precise lattice parameter values. The cAp lattice parameters did not correlate with age-at-death estimated from dental wear. From WAXS and the 00.2 diffraction peak widths, four archeological mc2s possessed coherently scattering domain lengths (crystallite c-axis sizes) that matched that of the modern mc2; SAXS revealed the same four archeological mc2 had D-period peak intensities equivalent to that of the modern mc2. The other eight archeological mc2s had significantly larger crystallite sizes (than the modern mc2) and weak or absent D-period peaks, differences attributed to diagenetic changes. Based on these data, the authors suggest that WAXS 00.2 peak width and SAXS D-period peak intensity can be used with intact bones to select those likely to retain largely unaltered tissue nanostructure, which might be required for other analyses. Taken as a whole, the results suggest detecting age-related deterioration in nanostructural features may be difficult in bone showing significant bioerosion.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2025

Understanding how defects and dopant atoms in copper surface oxides affect reactivity

Copper and its oxides are key catalytic materials, on which reactions often occur at the metal/oxide interface. Here, in this work, we directly connect the induction period observed during methanol-driven reduction of thin-film copper oxides to their atomic-scale structural order. Using temperature-programmed desorption (TPD) methanol titrations combined with scanning tunneling microscopy, we show that highly ordered oxide phases – particularly the “29” structure with its low defect density – exhibit long induction periods and initially low reactivity. The induction period, defined as the number of methanol TPD cycles required to reach half of the maximum formaldehyde yield, scales with oxide order and oxygen coverage. Enhanced reactivity of well-ordered oxides emerges only after repeated methanol adsorption/desorption cycles generate oxygen vacancies and new Cu(111)/Cu x O interfacial sites. In contrast, disordered or sub-stoichiometric oxides, which contain more intrinsic defects and interfaces, are active from the first TPD cycle. We further examine how dilute Pt and Rh dopants influence oxide order and reactivity: 1% Pt increases defect density and catalytic activity, while 1% Rh promotes oxide ordering and longer induction periods. These findings demonstrate that dilute alloying provides a potential method for tuning the structure and reactivity of Cu(111)/Cu x O interfaces.

Cu(111)Methanol oxidation↗
DOE OSTI2022

Impurity screening and Friedel oscillations in Floquet-driven two-dimensional metals

We develop a theory for the non-equilibrium screening of a charged impurity in a two-dimensional electron system under a strong time-periodic drive. Our analysis of the time-averaged polarization function and dielectric function reveals that Floquet driving modifies the screened impurity potential in two main regimes. In the weak drive regime, the time-averaged screened potential exhibits unconventional Friedel oscillations with multiple spatial periods contributed by a principal period modulated by higher-order periods, which are due to the emergence of additional Kohn anomalies in the polarization function. In the strong drive regime, the time-averaged impurity potential becomes almost unscreened and does not exhibit Friedel oscillations. Furthermore, this tunable Friedel oscillations is a result of the dynamic gating effect of the time-dependent driving field on the two-dimensional electron system.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

Using Strain in the Vadose Zone During Sinusoidal-Rate Well Tests to Characterize an Aquifer

Well tests with sinusoidal variation in pumping rate are appealing because of the ability to filter competing signals and they also can be set up with no net removal of water, which is an advantage at contaminated sites. Current methods of interpreting sinusoidal-rate well tests analyze the phase lag and amplitude of pressure signals in the aquifer. This approach is restricted by the availability of monitoring wells. An alternative approach could be to measure the stain in the vadose zone instead of pressure in the aquifer. The goal of this research is to evaluate the feasibility of using signals from in-situ sensors in the vadose zone to characterize sinusoidal-rate well tests. Well tests with a variety of periods and amplitudes were performed using a well completed in an unconfined saprolite aquifer with a water table at 9m depth at our test site near Clemson, South Carolina. Vertical strain was measured every meter between 1 to 7 m using a distributed optical fiber sensing technique called coherence length gated microwave photonics interferometry (CMPI). A novel packaging system was designed to couple the optical fiber to the saprolite. The results show that strain measured in the vadose zone during sinusoidal-rate pumping varies periodically with amplitudes from 0.1 to 3 microstrain, depending on the amplitude of the sinusoidal pumping rate. Hydraulic heads at similar radial distances were also periodic with changing amplitudes ranging up to ±0.2 meters. The strain in the vadose zone responds to pressure in the underlying aquifer—the strain is tensile when the pressure is falling, and compressive when the pressure is rising. The phase lag between the pumping rate and the strain was 3 min with a pumping period of 15 minutes, and this gives a hydraulic diffusivity of 0.06 m2/s. This is the same diffusivity obtained by analyzing a constant rate pumping test, and it is slightly larger than the diffusivity (0.044 m2/s) obtained by analyzing the phase lag of the pressure in a nearby monitoring well. Strain measurements in the vadose zone could be made with much less cost and effort than installing a deep monitoring well. These results suggest that strain measurements in the vadose zone could be used to supplement monitoring wells, thereby increasing the resolution, and reducing costs of analyzing well tests.

well test, sinusoidal, periodic pumping, strain↗
DOE OSTI2020

The Tidal Clock Effects of the Lunisolar Gravitational Field and the Earth’s Tidal Deformation

State-of-the-art atomic and optical clocks have the great potential to precisely test fundamental physical assumptions and enhance our understanding of nature. Their widespread applications require us to rigorously deduce the relativistic frequency shift in the framework of general relativity. One interesting question for clocks is tidal field effect which contains variously periodic variations. By introducing tidally deformed Earth, the general algorithms for clock and frequency comparisons are derived. We investigate the effects of external gravitational bodies and Earth’s tidal deformation on the ground- and space-based clocks. The orbital elements of satellite and locations of laboratory are introduced for corresponding clocks, in which the secular and long-period terms and short-period terms are subsequently presented. These effects called tidal clock effects produce the non-negligible contributions in the modern clock experiments and can be directly evaluated from our parameterized formulas. In addition, we also demonstrate for tidal clock effects the position and distance dependences, as well as periodicity. These tidal effects and dependences can provide valuable information for the clock comparison experiments.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

An Asymmetric Eclipse Seen toward the Pre-main-sequence Binary System V928 Tau

K2 observations of the weak-lined T Tauri binary V928 Tau A and B show the detection of a single, asymmetric eclipse, which may be due to a previously unknown substellar companion eclipsing one component of the binary with an orbital period >66 days. Over an interval of about 9 hr, one component of the binary dims by around 60%, returning to its normal brightness about 5 hr later. From modeling of the eclipse shape, we find evidence that the eclipsing companion may be surrounded by a disk or a vast ring system. The modeled disk has a radius of 0.9923 ± 0.0005 R∗, with an inclination of 56.°78 ± 0.°03, a tilt of 41.°22 ± 0.°05, an impact parameter of -0.2506 ± 0.0002 R∗, and an opacity of 1.00. The occulting disk must also move at a transverse velocity of 6.637 ± 0.002 R∗ day -1 , which, depending on whether it orbits V928 Tau A or B, corresponds to approximately 73.53 or 69.26 km s -1 . A search in ground-based archival data reveals additional dimming events, some of which suggest periodicity, but no unambiguous period associated with the eclipse observed by K2. We present a new epoch of astrometry that is used to further refine the orbit of the binary, presenting a new lower bound of 67 yr, and constraints on the possible orbital periods of the eclipsing companion. The binary is also separated by 18' (~2250 au) from the lower-mass CFHT-BD-Tau 7, which is likely associated with V928 Tau A and B. We also present new high-dispersion optical spectroscopy that we use to characterize the unresolved stellar binary.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Defect‐Free Sound Insulator Using Single Metal‐Based Friction Stir Process Array

Metals are excellent conductors for phonon transportation such as vibration, sound, and heat. Generally, metal sound insulators require multimaterial structure or defects and unimetal sound insulators are challenging. Therefore, a design of a defect-free sound insulator made by single alloys with multiple friction stir processes (FSPs) is proposed. Periodic friction stir processing can induce superlattice-like local mechanical properties’ modifications. By experimental acoustic characterization, it is observed that FSP can introduce clear acoustic–elastic property contrast on an aluminum plate by the presence of stir zone and heat-affected zones. In numerical simulations, the signature FSP-induced property profile is periodically and parallelly arranged on a long aluminum plate. The transmission gap frequencies are present on the frequency spectrum with the sound propagation direction perpendicular to the FSP paths. Disorder offsets on FSP periodicity are further introduced. Anderson localization is found on a resonance frequency, which provides −11 dB sound reduction by an exponential decay. Due to the finite design length, the slight disorder can also enhance sound insulation in the periodic transmission gap frequency. In conclusion, with analysis and comparison with different configurations, the best performance in the models can achieve −30 dB sound insulation in the 350 mm-long aluminum alloy plate with 14 parallel FSPs.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Response of the Quasi–Biennial Oscillation to a warming climate in global climate models

We compare the response of the Quasi–Biennial Oscillation (QBO) to a warming climate in eleven atmosphere general circulation models that performed time–slice simulations for present–day, doubled, and quadrupled CO 2 climates. No consistency was found among the models for the QBO period response, with the period decreasing by 8 months in some models and lengthening by up to 13 months in others in the doubled CO 2 simulations. In the quadrupled CO 2 simulations, a reduction in QBO period of 14 months was found in some models, whereas in several others the tropical oscillation no longer resembled the present–day QBO, although it could still be identified in the deseasonalized zonal mean zonal wind timeseries. In contrast, all the models projected a decrease in the QBO amplitude in a warmer climate with the largest relative decrease near 60 hPa. In simulations with doubled and quadrupled CO 2 , the multi–model mean QBO amplitudes decreased by 36 and 51%, respectively. Across the models the differences in the QBO period response were most strongly related to how the gravity wave momentum flux entering the stratosphere and tropical vertical residual velocity responded to the increases in CO 2 amounts. Likewise it was found that the robust decrease in QBO amplitudes was correlated across the models to changes in vertical residual velocity, parametrized gravity wave momentum fluxes, and to some degree the resolved upward wave flux. Here, we argue that uncertainty in the representation of the parameterized gravity waves is the most likely cause of the spread among the eleven models in the QBO's response to climate change.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Evaluation of the Quasi‐Biennial Oscillation in global climate models for the SPARC QBO‐initiative

Abstract Quasi‐biennial oscillations (QBOs) in thirteen atmospheric general circulation models forced with both observed and annually repeating sea surface temperatures (SSTs) are evaluated. In most models the QBO period is close to, but shorter than, the observed period of 28 months. Amplitudes are within ±20 % of the observed QBO amplitude at 10 hPa, but typically about half of that observed at lower altitudes (50 and 70 hPa). For almost all models, the oscillation's amplitude profile shows an overall upward shift compared to reanalysis and its meridional extent is too narrow. Asymmetry in the duration of eastward and westward phases is reasonably well captured, though not all models replicate the observed slowing of the descending westward shear. Westward phases are generally too weak, and most models have an eastward time mean wind bias throughout the depth of the QBO. The intercycle period variability is realistic and in some models is enhanced in the experiment with observed SSTs compared to the experiment with repeated annual cycle SSTs. Mean periods are also sensitive to this difference between SSTs, but only when parametrized non‐orographic gravity wave (NOGW) sources are coupled to tropospheric parameters and not prescribed with a fixed value. Overall, however, modelled QBOs are very similar whether or not the prescribed SSTs vary interannually. A portrait of the overall ensemble performance is provided by a normalized grading of QBO metrics. To simulate a QBO, all but one model used parametrized NOGWs, which provided the majority of the total wave forcing at altitudes above 70 hPa in most models. Hence the representation of NOGWs either explicitly or through parametrization is still a major uncertainty underlying QBO simulation in these present‐day experiments.

Bushell, A. C.↗
DOE OSTI2020

Performance of the missing transverse momentum triggers for the ATLAS detector during Run-2 data taking

The factor of four increase in the LHC luminosity, from 0.5 × 10 34 cm -2 s -1 to 2.0 × 10 34 cm -2 s -1 , and the corresponding increase in pile-up collisions during the 2015–2018 data-taking period, presented a challenge for the ATLAS trigger, particularly for those algorithms that select events with missing transverse momentum. The output data rate at fixed threshold typically increases exponentially with the number of pile-up collisions, so the legacy algorithms from previous LHC data-taking periods had to be tuned and new approaches developed to maintain the high trigger efficiency achieved in earlier operations. A study of the trigger performance and comparisons with simulations show that these changes resulted in event selection efficiencies of > 98% for this period, meeting and in some cases exceeding the performance of similar triggers in earlier run periods, while at the same time keeping the necessary bandwidth within acceptable limits.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2026

Benchmarking Cylindrical Blast Wave Theory Against the OSIRIS-REx Sample Return Capsule Reentry

Weak shock theory based on cylindrical blast waves has been used to interpret meteor infrasound, but it has not been systematically benchmarked against a non-ablating hypersonic source with independently known parameters. The objective of this study is not to propose a new theoretical framework, but to evaluate the operational validity of the existing suite of blast radius formulations against a high-fidelity ground truth dataset. The OSIRIS-REx Sample Return Capsule reentry on 24 September 2023 provides such a benchmark because the capsule geometry, trajectory, and infrasound emission points are constrained from mission data and ray tracing, reducing source-side uncertainty associated with ablation. Using observations from 39 infrasound stations, this benchmarking study evaluates six published blast radius (${R}_{0}$ ) formulations and three weak-shock transition coefficients ( C ) within a stratified atmospheric propagation model to predict signal period and peak overpressure. The benchmarking identifies the Sakurai formulation as the best-performing formulation for non-ablating bodies, with the Jones/Plooster formulation performing comparably when a physically appropriate C is adopted. Sakurai and Jones/Plooster yield linear-period median absolute percentage residuals of 9% and 11%, respectively. The period predictions show only weak sensitivity to C at these propagation distances. The Mach-diameter approximation commonly used in meteor studies overestimates ${R}_{0}$ by more than a factor of 3 in the absence of ablation. Finally, these results establish a performance baseline for applying cylindrical blast wave theory to effectively non-ablating hypersonic bodies and demonstrate that the signal period is a robust observable for constraining ${R}_{0}$.

atmospheric entry↗
DOE OSTI2020

Robust late twenty-first century shift in the regional monsoons in RegCM-CORDEX simulations

Abstract We use an unprecedented ensemble of regional climate model (RCM) projections over seven regional CORDEX domains to provide, for the first time, an RCM-based global view of monsoon changes at various levels of increased greenhouse gas (GHG) forcing. All regional simulations are conducted using RegCM4 at a 25 km horizontal grid spacing using lateral and lower boundary forcing from three General Circulation Models (GCMs), which are part of the fifth phase of the Coupled Model Inter-comparison Project (CMIP5). Each simulation covers the period from 1970 through 2100 under two Representative Concentration Pathways (RCP2.6 and RCP8.5). Regional climate simulations exhibit high fidelity in capturing key characteristics of precipitation and atmospheric dynamics across monsoon regions in the historical period. In the future period, regional monsoons exhibit a spatially robust delay in the monsoon onset, an increase in seasonality, and a reduction in the rainy season length at higher levels of radiative forcing. All regions with substantial delays in the monsoon onset exhibit a decrease in pre-monsoon precipitation, indicating a strong connection between pre-monsoon drying and a shift in the monsoon onset. The weakening of latent heat driven atmospheric warming during the pre-monsoon period delays the overturning of atmospheric subsidence in the monsoon regions, which defers their transitioning into deep convective states. Monsoon changes under the RCP2.6 scenario are mostly within the baseline variability.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Photosynthesis phenology, as defined by solar-induced chlorophyll fluorescence, is overestimated by vegetation indices in the extratropical Northern Hemisphere

Vegetation phenology is highly sensitive to climate change, although the data and methods used to estimate key phenological states can influence this sensitivity. Because of its direct relation to leaf photosynthetic carbon uptake, remotely sensed solar-induced chlorophyll fluorescence (SIF) can provide new insight assessing changes in vegetation phenology. In this work, we investigated the potential of using a SIF time series product named contiguous SIF (CSIF) to estimate spring, summer, and autumn phenology in the extratropical Northern Hemisphere (>30°N) and compared the results with those based on Moderate Resolution Imaging Spectroradiometer (MODIS) Normalized Difference Vegetation Index (NDVI) for the period 2001–2018. Overall, we found similar spatial patterns in phenological states. However, specific dates of key phenological events differed when using CSIF vs. MODIS NDVI data. NDVI data indicated that the growing season started earlier (by 10.1 days on average) and ended later (11.5 days on average) relative to CSIF data. This implies that actual periods for photosynthetic activity are shorter (by 21.6 days on average) than those estimated from vegetation indices more directly related to changes in canopy structure. These large differences between results from NDVI and that from CSIF suggest that vegetation indices such as NDVI seem to overestimate the period for active photosynthesis over the extratropical Northern Hemisphere. Furthermore, while phenology of the early growing season is dominated by temperature for both NDVI and CSIF data, phenology of the late growing season is mainly controlled by temperature for NDVI but by precipitation for CSIF. Our findings were further confirmed by other SIF (GOME-2 SIF) and vegetation index (MODIS EVI) datasets. Phenology modes in Earth system modelling are often parameterized using leaf unfolding and senescence from either station or satellite observations. Our results imply that canopy structure-based parameterization schemes may have overestimated photosynthesis active period, and thus productivity responses. We conclude that SIF data offers a novel and unique approach for assessing phenological change - one that is more directly tied to the carbon cycle and how it is being influenced by climate change.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Intra-annual variations of regional total column ozone, aerosol optical depth, and water vapor from ground-based, satellite-based and model-based observations

The Microtops II ozonometer has been deployed to continuously measure the total column ozone (TCO), pre-cipitable water content (PWC) and aerosol optical thickness (AOT) at Atigre village (16.74° N, 74.37° E, 604 m above sea level, masl) located on the southeastern slope of Indian western Ghats. The Microtops II ozonometer measurements during the period from September 2017 to June 2018 were analyzed along with the retrieved products of Ozone Monitoring Instrument (OMI) and Moderate Resolution Imaging Spectroradiometer (MODIS) onboard NASA's Terra satellites and also the ERA-interim reanalysis model. We found that the TCO (or PWC) data from the OMI (or MODIS) and the ERA reanalysis model data products are in good agreement. The mea-surements showed strong seasonal variability such that TCO was lowest in the winter season but increased in the summer season, and both PWC and AOT were lowest in the post-monsoon season and increased in the summer season. In continuous ERA-interim reanalysis TCO observations, we found the three types of periodicities in the daily mean TCO as well as PWC (i) The weekly (7 days), (ii) quasi-biweekly (14 to 16 days), (iii) Madden Julian Oscillations (MJO) (30–60 days) oscillations. These oscillations depend on the periodic weather changes in the troposphere. We also observe that the average day to day variability in daily mean TCO observations in the observational period is about 1.4%, 1.1%, 2.62% in the Microtops II ozonometer, ERA-interim reanalysis, and OMI data respectively. We found the positive correlation between daily mean AOT and PWC in all seasons (except monsoon), which is maximum in the winter season, showing the hygroscopic nature of aerosols. The sources of water vapor and aerosol at our location are also studied using back-trajectory analysis.

Nade, Dada R.↗
DOE OSTI2020

Dislocation-graphene interactions in Cu/graphene composites and the effect of boundary conditions: a molecular dynamics study

Here, in this paper, systematic molecular dynamics simulations were employed to study dislocation-graphene interactions in Cu/graphene systems using two carefully selected configurations: (i) edge dislocations in the pileup interacting with graphene during nanoindentation and (ii) complex dislocations interacting with graphene in compression of layered nanopillars. The intrinsic and extrinsic size effects were investigated with respect to varying the Cu lamella thickness and pillar sizes, revealing an anomalous extrinsic size effect: the smaller, the weaker. To understand which boundary conditions are physically consistent, both free and periodic graphene were considered. It was seen that edge dislocations in the pileup continuously transmitted across the free graphene through Cu/graphene interfacial sliding and graphene reorientation, whereas transmission was very difficult for periodic graphene as its in-plane deformation was required. In compression of Cu/graphene nanopillars, the strengthening effect of free graphene was found to be less obvious than that of the periodic case, which could be attributed to free graphene edges acting as dislocation sources after the compressed Cu overflowed the graphene sheet from all directions. We therefore conclude that the strengthening effect of periodic graphene inclusions appears to be overestimated, such that free graphene should be more appropriate.

36 MATERIALS SCIENCE↗
DOE OSTI2020

ELSI — An open infrastructure for electronic structure solvers

Routine applications of electronic structure theory to molecules and periodic systems need to compute the electron density from given Hamiltonian and, in case of non-orthogonal basis sets, overlap matrices. System sizes can range from few to thousands or, in some examples, millions of atoms. Different discretization schemes (basis sets) and different system geometries (finite non-periodic vs. infinite periodic boundary conditions) yield matrices with different structures. The ELectronic Structure Infrastructure (ELSI) project provides an open-source software interface to facilitate the implementation and optimal use of high-performance solver libraries covering cubic scaling eigensolvers, linear scaling density-matrix-based algorithms, and other reduced scaling methods in between. In this paper, we present recent improvements and developments inside ELSI, mainly covering (1) new solvers connected to the interface, (2) matrix layout and communication adapted for parallel calculations of periodic and/or spin-polarized systems, (3) routines for density matrix extrapolation in geometry optimization and molecular dynamics calculations, and (4) general utilities such as parallel matrix I/O and JSON output. The ELSI interface has been integrated into four electronic structure code projects (DFTB+, DGDFT, FHI-aims, SIESTA), allowing us to rigorously benchmark the performance of the solvers on an equal footing. Based on results of a systematic set of large-scale benchmarks performed with Kohn–Sham density-functional theory and density-functional tight-binding theory, we identify factors that strongly affect the efficiency of the solvers, and propose a decision layer that assists with the solver selection process. As a result, we describe a reverse communication interface encoding matrix-free iterative solver strategies that are amenable, e.g., for use with planewave basis sets.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Techno-economic feasibility of borehole thermal energy storage system connected to geothermal heat pumps for seasonal heating load of two buildings in Fairbanks, Alaska

Borehole thermal energy storage (BTES) is an effective solution for managing imbalanced heating and cooling loads in cold regions. This study evaluated the long-term feasibility of a BTES system in the Fairbanks area, Alaska, through building energy modeling, resource characterization, and numerical modeling. The system was designed to store waste heat from a nearby coal power plant during summer and provide thermal energy during winter to geothermal heat pumps supplying heating loads in two buildings. Heating load profiles were modeled for the buildings using EnergyPlus, and the results indicated the annual heating load was 5.6 times greater than the cooling load. 40 borehole heat exchangers were pre-designed approximately 100 m away from the two buildings in terms of land availability and regulatorily optimized depth. The 20-year performance of the designed BTES system under two operational scenarios—one with a 5-year preheating period and one without—was numerically modeled using subsurface temperature and properties characterized through the literature review and thermal response tests. Both scenarios demonstrated that the BTES has the capacity to fully cover the heating loads in the two buildings throughout the 20-year lifetime. Production temperatures at central wells were 33 % higher on average than at outer wells in both scenarios. The 5-year preheating period increased subsurface and extraction temperatures, and correspondingly annual average and total thermal energy production was higher for 8 years than in the scenario without the preheating period. These results highlight the long-term reliability and sustainability of the BTES system in meeting heating demands over its lifetime, with the preheating period offering potential performance improvements. Implementing the BTES system in cold regions with high heating demand, such as Fairbanks, Alaska, could provide a long-term, sustainable energy solution for managing imbalanced heating and cooling loads.

15 GEOTHERMAL ENERGY↗