Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “inversion methods”

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

BayesMT: A Probabilistic Bayesian Framework for the Seismic Moment Tensor

Moment tensors (MTs) have long been used in earthquake and explosion source analysis, and there has been a renewed interest in how they can inform us about the seismic source, particularly in the geophysical monitoring community due to its application in event identification and yield analysis. However, parameter uncertainties in seismic MT inversion are rarely available. The inverse procedure often does not quantify MT model errors such as event location, data noise and Earth model that are essential for estimating solution robustness. To address this need, we propose to adopt the Bayesian probabilistic framework to incorporate uncertainties in MT inversions. In this study, we present the theoretical background of a probabilistic Bayesian framework for MT inversion accounting for model and measurements errors and illustrate the implementation of the method using a synthetic example.

58 GEOSCIENCES↗
DOE OSTI2025

Algorithm to extract direction in 2D discrete distributions and a continuous Frobenius norm

In this study, we present a novel algorithm for determining directionality in 2D distributions of discrete data. We compare a reference dataset with a known direction to a measured dataset with an unknown direction by the Frobenius norm of the difference (FND) to find the unknown direction. To generalize this concept, we develop a continuous Frobenius norm of the difference (CFND) as a continuous analog of the FND and derive its analytical expression. By relating fitted and normalized 2D Gaussian distributions, we show that the CFND approximates the FND, and we validate this relationship with computer simulations. We find that a first-order approximation of the CFND between two similar Gaussian distributions takes the form of an absolute sine function, offering a simple analytical form with potential for specialized applications in segmented inverse beta decay (IBD) neutrino detectors, astronomy, machine learning, and more. Although this method may easily extend to 3D scalar fields, our focus here is on 2D real-valued fields as it directly applies to directionality. Our methodology consists of modeling a 2D Gaussian distribution, binning the data into a histogram, and encoding it as a square matrix. Rotating this matrix around its geometric center and comparing it to a measured dataset using the FND gives us rotational data that we fit with an absolute sine function. The location of the minimum of this fit is the angle closest to the true angle of the direction in the measured dataset. We present the derivation and discuss initial applications of the CFND in our novel algorithm, demonstrating its success in approximating directionality in 2D distributions.

Data Analysis, Statistics and Probability (physics↗
DOE OSTI2022

The weathered bedrock vadose zone: A hidden control on water availability in the western United States (Final Technical Report)

This report summarizes the research completed during the project entitled "The weathered bedrock vadose zone: A hidden control on water availability in the western United States" (Award: DE-SC0018039). The report covers work completed between 9/15/17 and 9/14/20. The objective of the exploratory project work was to evaluate the overarching hypothesis that the subsurface underlying actively eroding landscapes is systematically structured as a result of a balance between weathering and erosion and that this structure leads to dynamic moisture storage within weathered rock. To accomplish project objectives, we characterized hillslope weathering profile structure and conducted subsurface moisture monitoring across four study sites with comparable bedrock types but differing erosion rate and hydroclimate: The Angelo Coast Range Reserve (ACRR), Sagehorn-Russell Ranch (SRR), Antelope Valley Ranch (AVR), and the LBNL Watershed Function Scientific Focus Area (SFA). New characterization and monitoring infrastructure at the AVR and SFA sites was established through this project and enables long-term monitoring of hydrologic and carbon fluxes at the sites for comparison with long-term monitoring in other sites. Deep drilling, downhole logging, and characterization of chemistry and porosity of rock recovered from drilling reveals two weathering fronts across the sites. A shallower front is marked by the extent of pervasive fracturing and oxidation of pyrite and organic carbon. A deeper weathering front is marked by the appearance of open fractures and discoloration. At the AVR site, multiple ridge tops of different hillslope length were investigated to establish the relationship between the weathering front and topography. Drilling reveals that the fraction of the hillslope relief that is weathered scales with hillslope length, consistent with the theoretical framework that guided the project relating groundwater drainage with landscape evolution. With the establishment of new boreholes, we have documented significant water storage in weathered bedrock using low-field borehole nuclear magnetic resonance (NMR) across all four sites. In the upper 4 m of bedrock, the average water contents at ACRR, AVR, SRR, and SFA are approximately 25%, 16%, 25%, 23% respectively. At the ACRR site, where long-term monitoring via neutron probe surveys and time-domain transmission sensors are available, we have demonstrated that NMR reliably records moisture content and importantly, moisture content dynamics in weathered bedrock. Across all sites, our NMR monitoring has also revealed that a considerable fraction of the water storage occurs in fractures or pores that are significantly larger than the fine-grained bedrock matrix of the parent rock. To identify the proportion of water storage occurring in fractures and large pores, we analyzed the NMR signal by evaluating both the sum of echoes and an inversion of the NMR signal for a distribution of T2 relaxation times. Both methods reveal that at least 20% of the water storage occurs in fractures and larger pores across all sites. To constrain the interpretation of the in-situ field-scale NMR measurements, we compare the NMR response in the field to independent pore-structure information obtained via helium pycnometry. Results support the inference that dynamic seasonal water storage is dominantly restricted to the fracture network, which can reach 7% of the total volume. This dynamic range in water content is comparable to that of some soils, underscoring the need to understand how this region of the Earth System is structured, and functions as a hydrologic reservoir.

58 GEOSCIENCES↗
DOE OSTI2020

Microcanonical rates from ring-polymer molecular dynamics: Direct-shooting, stationary-phase, and maximum-entropy approaches

We address the calculation of microcanonical reaction rates for processes involving significant nuclear quantum effects using ring-polymer molecular dynamics (RPMD), both with and without electronically non-adiabatic transitions. After illustrating the shortcomings of the naive free-particle direct-shooting method, in which the temperature of the internal ring-polymer modes is set to the translational energy scale, we investigate alternative strategies based on the expression for the microcanonical rate in terms of the inverse Laplace transform of the thermal reaction rate. It is shown that simple application of the stationary-phase approximation (SPA) dramatically improves the performance of the microcanonical rates using RPMD, particularly in the low-energy region where tunneling dominates. Using the SPA as a Bayesian prior, numerically exact RPMD microcanonical rates are then obtained using maximum entropy inversion of the thermal reaction rates for both electronically adiabatic and non-adiabatic model systems. Finally, the direct-shooting method is revisited using the SPA-determined temperature for the internal ring-polymer modes, leading to a simple, direct-simulation method with improved accuracy in the tunneling regime. This work suggests a general strategy for the extraction of microcanonical dynamical quantities from RPMD (or other approximate thermal) simulations.

Tao, Xuecheng↗
DOE OSTI2025

Applying Gaussian Process Machine Learning and Modern Probabilistic Programming to Satellite Data to Infer CO 2 Emissions

Satellite data provides essential insights into the spatiotemporal distribution of CO 2 concentrations. However, many atmospheric inverse models fail to adequately incorporate the spatial and temporal correlations inherent in satellite observations and often lack rigorous methods for estimating parameters like spatial length scales. We introduce an inference model that processes the spatiotemporal covariance in satellite data and estimates hyperparameters such as covariance length scales. Our approach uses the Gaussian process (GP) machine learning (ML) and modern probabilistic programming languages (PPLs) to perform atmospheric inversions of emissions from satellite data. We develop a GP ML inversion system based on modern PPLs and the GEOS-Chem chemical transport model, simulating atmospheric CO 2 concentrations corresponding to the Orbiting Carbon Observatory-2/3 (OCO-2/3) data for July 2020. In our supervised learning framework, we treat the GEOS-Chem simulated data set as the target, with predictors derived by scaling the target with sector-specific factors hidden from the GP machine. Our results show that the GP model, combined with GPU-enabled PPLs, effectively retrieves true emission scaling factors and infers noise levels concealed within the data. This suggests that our method could be applied over larger areas with more complex covariance structures, enabling comprehensive analysis of the spatiotemporal patterns observed in OCO-2/3 and similar satellite data sets.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Search for magnetically-induced signatures in the arrival directions of ultra-high-energy cosmic rays measured at the Pierre Auger Observatory

We search for signals of magnetically-induced effects in the arrival directions of ultra-high-energy cosmic rays detected at the Pierre Auger Observatory. We apply two different methods. One is a search for sets of events that show a correlation between their arrival direction and the inverse of their energy, which would be expected if they come from the same point-like source, they have the same electric charge and their deflection is relatively small and coherent. We refer to these sets of events as “multiplets”. The second method, called “thrust”, is a principal axis analysis aimed to detect the elongated patterns in a region of interest. We study the sensitivity of both methods using a benchmark simulation and we apply them to data in two different searches. The first search is done assuming as source candidates a list of nearby active galactic nuclei and starburst galaxies. The second is an all-sky blind search. We report the results and we find no statistically significant features. We discuss the compatibility of these results with the indications on the mass composition inferred from data of the Pierre Auger Observatory.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Extraction of mechanical properties of materials through deep learning from instrumented indentation

Instrumented indentation has been developed and widely utilized as one of the most versatile and practical means of extracting mechanical properties of materials. This method is particularly desirable for those applications where it is difficult to experimentally determine the mechanical properties using stress–strain data obtained from coupon specimens. Such applications include material processing and manufacturing of small and large engineering components and structures involving the following: three-dimensional (3D) printing, thin-film and multilayered structures, and integrated manufacturing of materials for coupled mechanical and functional properties. Here, we utilize the latest developments in neural networks, including a multifidelity approach whereby deep-learning algorithms are trained to extract elastoplastic properties of metals and alloys from instrumented indentation results using multiple datasets for desired levels of improved accuracy. We have established algorithms for solving inverse problems by recourse to single, dual, and multiple indentation and demonstrate that these algorithms significantly outperform traditional brute force computations and function-fitting methods. Moreover, we present several multifidelity approaches specifically for solving the inverse indentation problem which 1) significantly reduce the number of high-fidelity datasets required to achieve a given level of accuracy, 2) utilize known physical and scaling laws to improve training efficiency and accuracy, and 3) integrate simulation and experimental data for training disparate datasets to learn and minimize systematic errors. The predictive capabilities and advantages of these multifidelity methods have been assessed by direct comparisons with experimental results for indentation for different commercial alloys, including two wrought aluminum alloys and several 3D printed titanium alloys.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Extraction of mechanical properties of materials through deep learning from instrumented indentation

Instrumented indentation has been developed and widely utilized as one of the most versatile and practical means of extracting mechanical properties of materials. This method is particularly desirable for those applications where it is difficult to experimentally determine the mechanical properties using stress–strain data obtained from coupon specimens. Such applications include material processing and manufacturing of small and large engineering components and structures involving the following: three dimensional (3D) printing, thin-film and multilayered structures, and integrated manufacturing of materials for coupled mechanical and functional properties. Here, we utilize the latest developments in neural networks, including a multifidelity approach whereby deep-learning algorithms are trained to extract elastoplastic properties of metals an d alloys from instrumented indentation results using multiple datasets for desired levels of improved accuracy. We have established algorithms for solving inverse problems by recourse to single, dual, and multiple indentation and demonstrate that these algorithms significantly outperform traditional brute force computations and function-fitting methods. Moreover, we present several multifidelity approaches specifically for solving the inverse indentation problem which 1) significantly reduce the number of high-fidelity datasets required to achieve a given level of accuracy, 2) utilize known physical and scaling laws to improve training efficiency and accuracy, and 3) integrate simulation and experimental data for training disparate datasets to learn and minimize systematic errors. The predictive capabilities and advantages of these multifidelity methods have been assessed by direct comparisons with experimental results for indentation for different commercial alloys, including two wrought aluminum alloys and several 3D printed titanium alloys.

3D printed materials, stress-strain behavior, mult↗
DOE OSTI2023

Real-space inversion and super-resolution of ultrafast scattering

Ultrafast scattering using x-rays or electrons is an emerging method to obtain structure dynamics at the atomic lengthscales and timescales. However, directly resolving in real-space atomic motions is inherently limited by the finite detector range and the probe energy. So, as a result, the time-resolved signal interpretation is mostly done in reciprocal space and relies on modeling and simulations of specific structures and processes. Here, we introduce a model-free approach to directly resolve scattering signals in real space, surpassing the diffraction limit, using scattering kernels and signal priors that naturally arise from the measurement constraints. We demonstrate the approach on simulated and experimental data, recover multiple atomic motions at sub-angstrom resolutions, and discuss the recovery accuracy and resolution limits versus signal fidelity. The approach offers a robust path to obtain high-resolution real-space information of atomic-scale structure dynamics using current time-resolved x-ray or electron scattering sources.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Extension of the flux fit method for estimating power deposition profiles

The flux fit method is used to self-consistently estimate the power deposition profile and heat transport profiles from temperature measurements originating from perturbative experiments with a modulated source. Furthermore, this letter improves on this method by addressing the limitations and assumptions. The most crucial improvement is the additional freedom in the source deposition profile. Allowing for a variable central deposition location and height and including a skewness parameter produces deposition profiles more consistent with the measurement data, but still wider than equilibrium ray tracing in two different DIII-D discharges. Moreover, we show that the quality of the estimated deposition profile is key to the accuracy of diffusivity and connectivity estimates, but inversely, the estimated transport parameters hardly affect the quality of the power deposition estimate. Using this method, we show that the power deposition profile estimate is broadened with respect to ray-tracing by about 1.7-1.8 times in two DIII-D discharges.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Nuclear cusps and singularities in the nonadditive kinetic potential bifunctional from analytical inversion

The non-additive kinetic potential $v^{\text{NAD}}$ is a key quantity in density-functional theory (DFT) embedding methods, such as frozen density embedding theory and partition DFT. $v^{\text{NAD}}$ is a bi-functional of electron densities $\rho_{\rm B}$ and $\rho_{\rm tot} = \rho_{\rm A} + \rho_{\rm B}$. It can be evaluated using approximate kinetic-energy functionals, but accurate approximations are challenging. The behavior of $v^{\text{NAD}}$ in the vicinity of the nuclei has long been questioned, and singularities were seen in some approximate calculations. In this article, the existence of singularities in $v^{\text{NAD}}$ is analyzed analytically for various choices of $\rho_{\rm B}$ and $\rho_{\rm tot}$, using the nuclear cusp conditions for the density and Kohn-Sham potential. It is shown that no singularities arise from smoothly partitioned ground-state Kohn-Sham densities. We confirm this result by numerical calculations on diatomic test systems HeHe, HeLi$^+$, and H$_2$, using analytical inversion to obtain a numerically exact $v^{\rm NAD}$ for the local density approximation. We examine features of $v^{\rm NAD}$ which can be used for development and testing of approximations to $v^{\rm NAD}[\rho_{\rm B},\rho_{\rm tot}]$ and kinetic-energy functionals.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2023

Joint Inversion of Surface Electrical Resistivity Tomography and Seismic Refraction Data between the 200 Areas

Geologic stratigraphy on the Hanford Site influences groundwater and contaminant migration through the aquifer system and the vadose zone. The current geologic framework model (GFM) relies heavily on a sparse distribution of borehole data in some locations to map geologic contacts and hydrologic properties in the subsurface. Non-invasive geophysical methods such as electrical resistivity tomography (ERT), transient electromagnetic surveying, and seismic imaging are being used at Hanford to map subsurface structure in areas with limited well observations. This is to develop and mature the capability of geophysical methods to aid in GFM refinement, to identify regions of subsurface complexity, and for optimal well siting. A joint inversion of co-located seismic refraction and ERT data was carried out for data collected on a ~2.3-km profile between the 200 Areas on the Hanford Site. While ERT and seismic refraction images have sensitivity to overlapping physical properties (porosity, moisture content, lithology), the resolution and physics used to acquire each of these datasets are different and therefore information can be different or mutually complementary. Performing a joint inversion provides a reasonable option for a coherent, coupled interpretation for mutually complementary datasets. Between the 200 Areas, there are few boreholes to interpret the geologic framework model, and these data sets were obtained to provide a first line of evidence toward identifying stratigraphic structure. The seismic refraction and ERT data were independently inverted during fiscal year 2022 and broadly showed a two-layer structure with a trough-like feature that is ~1 km wide and upwards of 150 m deep. The depth of the trough feature was greater in the ERT image compared to the seismic image, which indicated a maximum depth of approximately 110 m. The objective of the joint inversion described in this report was to invert the seismic refraction and ERT data together while constraining the ERT image to be structurally similar to the seismic refraction image. The approach was applied using the geophysical inverse modeling program E4D, which has the capability to invert first-arrival times from seismic refraction data and ERT resistances using a “cross-gradient” constraint. The application of cross-gradient constraints with different weights produces ERT models that show a high degree of similarity within the upper 100 m (above ~120 m elevation). None of the ERT models show an improved structural similarity to the seismic result; therefore, it is recommended that further attempts to jointly interpret these models focus on petrophysics and image resolution. Petrophysical measurements of core samples would improve knowledge of what drives the ERT response in this region and, along with downhole geophysical measurements, could be used to “ground truth” the surface-based geophysical results. Image resolution studies would provide insight into which regions of the inverted images are reliable and which regions are poorly constrained.

58 GEOSCIENCES↗
DOE OSTI2020

Three-dimensional fracture continuum characterization aided by surface time-domain electromagnetics and hydrogeophysical joint inversion—proof-of-concept

Efficient and safe production of hydraulically fractured reservoirs benefits from the prediction of their geometrical attributes. Geophysical methods have the potential to provide data that are sensitive to fracture geometries, alleviating the typically sparse nature of in situ reservoir observations. Moreover, surface-based methods can be logistically and economically attractive since they avoid operational interference with the injection well infrastructure. Here, we investigate the potential of the surface-based time-domain electromagnetic (EM) method. EM methods can play an important role owing to their sensitivity to injection-induced fluid property changes. Two other advantageous factors are the EM signal-enhancing effect of vertical steel-cased wells and the fact that injected proppants can be enhanced to produce a stronger electrical conductivity contrast with the reservoir’s connate fluid. Nevertheless, an optimal fracture characterization will no doubt require the integration of EM and reservoir injection and production data. We hence carry out our investigations within a hydrogeophysical parameter estimation framework where EM data and injection flow rates are combined in a fully coupled way. Given the interdisciplinary nature of coupled hydrogeophysical inverse modeling, we dedicate one section to laying out key aspects in a didactic manner.

58 GEOSCIENCES↗
DOE OSTI2026

Analytical desmearing of Bonse–Hart ultra-small-angle neutron scattering data via truncated Abel inversion

A non-iterative analytical framework based on the truncated Abel inversion is developed for desmearing Bonse–Hart ultra-small-angle neutron scattering (USANS) data. The method directly inverts the slit-averaged intensity without empirical extrapolation or iterative regularization, establishing a closed-form relationship between the measured and intrinsic scattering profiles. Numerical benchmarks on representative models, including a rigid-line form factor, a Lorentzian function and a fractal structural model, demonstrate quantitative recovery of the ground-truth intensity across the full Q range. Application to a deuterated polystyrene/poly(2-vinylpyridine) blend further confirms that the approach yields smooth continuous profiles consistent with companion small-angle neutron scattering data. The truncated Abel inversion thus provides a stable, model-independent and physically transparent route for accurate desmearing of Bonse–Hart USANS measurements.

Huang, Guan-Rong [National Tsing Hua University, T↗
DOE OSTI2026

Importance of finite-size corrections for accurate ab initio modeling of carrier capture at semiconductor defects: A case study of substitutional C N in GaN

In ab initio studies of carrier-capture processes in defective semiconductor materials, the single-effective-mode formalism and the static-coupling approximation have become the predominant theoretical approaches for determining carrier-capture coefficients. The single-mode formalism relies on accurate nonequilibrium defect energies obtained from density-functional theory (DFT), where required inputs are a series of configurationally displaced, defect-containing supercells obtained using an interpolative ansatz, and where the DFT outputs are corresponding total energies that have traditionally been postprocessed using a long-established ground-state formulation of finite-size corrections and defect-formation energies. This formulation remains commonly used even though the defects that form a configuration-coordinate (CC) diagram typically exist as structures that are displaced from the ground state. To remedy this inconsistency, Kumagai has recently proposed novel methods for implementing finite-size corrections specifically intended for DFT calculations of the defect energies used to construct CC diagrams and implement the single-mode formalism [Y. Kumagai, Phys. Rev. B 107, L220101 (2023)]. Kumagai's approach builds on the latest finite-size-correction methods introduced to describe vertical charge-state transitions for charge-localizing point defects in semiconductors and insulators [T. Gake et al., Phys. Rev. B 101, 020102 (2020); S. Falletta et al., Phys. Rev. B 102, 041115 (2020)]. The newly identified finite-size artifact treated in these studies is the polarization charge induced on a configurationally frozen defect and its subsequent interaction with a vertical transition in charge state. In this work, we evaluate Kumagai's proposed methodology by applying it in a high-precision DFT study of carrier capture by substitutional C N in GaN, a well-characterized and technologically relevant defect and material. We have rigorously calculated C N defect energies across various supercell sizes for each defect configuration and charge state on the hole-capture CC diagram of C N (𝑞=−1), enabling a direct comparison of the slopes of the defect energies versus inverse cell size with those predicted by Kumagai. The most consequential prediction of Kumagai's method is that these slopes distinctly vary as the square of the linear-interpolation parameter used to construct the nonequilibrium defect configurations. Our results quantitatively support this prediction. Moreover, with these new finite-size corrections and multiple-cell-size DFT calculations in place, we find that the classical energy barrier for hole capture by C N (𝑞=−1) in GaN decreases to 0.092–0.127 eV. This finding confirms the recent ≈ 0.1 eV prediction of Reshchikov based on the weak temperature dependence for hole capture observed in photoluminescence experiments [M. A. Reshchikov, J. Appl. Phys. 129, 121101 (2021)]. These results stand in stark contrast to previously calculated barriers of 0.486 and 0.73 eV, which also used the single-mode formalism but were obtained by instead using ground-state-based finite-size corrections. Our reduced classical barrier for capture increases the temperature-dependent hole-capture coefficient of a C N (𝑞=−1) defect by more than two to four orders of magnitude for temperatures of 100–600 K, compared to the previous 0.486 eV results. While other defects may not be as dramatically affected as here, we suggest that incorporating proper finite-size corrections for the vertical-transition-like states embedded within CC diagrams is an essential, yet previously unrecognized, component of accurate modeling of carrier-capture when using the single-effective-mode formalism.

dielectric properties↗
DOE OSTI2022

Image-domain DAS 3D VSP elastic transmission tomography

SUMMARY Full-wavefield elastic imaging of active-source seismic data acquired by downhole receivers commonly offers higher-resolution subsurface images in the vicinity of a borehole compared to conventional surface seismic data sets, which can lack higher-frequency wavefield components due to longer travel paths and increased attenuation. An increasingly used approach for downhole acquisition is vertical seismic profiling (VSP), which has become more attractive when coupled with distributed acoustic sensing (DAS) using optical fibres installed in wells. The main difficulty for generating high-quality images with full-wavefield imaging tools for DAS VSP data, though, is the need for an accurate velocity model. To build plausible velocity models using active-source DAS VSP data, we adopt a 3-D image-domain elastic transmission tomography technique, originally developed for surface-recorded passive (microseismic) data, by exchanging the source and receiver positions (i.e. reciprocity) to mimic a passive-seismic surface monitoring scenario. The inversion approach exploits various images for each source constructed through time-reverse imaging (TRI) of downgoing P- and S-wave first-arrival waveforms. The TRI process uses the kinetic term of the (extended) PS energy imaging condition that exhibits sufficient sensitivity to velocity model errors. The method automatically updates the P- and S-wave velocity models to optimize image focusing via adjoint-state inversion. We illustrate the efficacy of the adopted elastic inversion technique using an active-source DAS 3-D VSP field data set acquired in the North Slope of Alaska. The numerical experiments demonstrate that the inverted elastic velocity models can be further used in full-wavefield acoustic/elastic imaging algorithms to obtain accurate subsurface images.

Oren, Can (ORCID:000000019534265X)↗
DOE OSTI2022

Tropopause Characteristics Based on Long-Term ARM Radiosonde Data: A Fine-Scale Comparison at the Extratropical SGP Site and Arctic NSA Site

The variations in the characteristics of the tropopause are sensitive indicators for the climate system and climate change. By using Atmospheric Radiation Measurement (ARM) radiosonde data that were recorded at the extratropical Southern Great Plains (SGP) and Arctic North Slope of Alaska (NSA) sites over an 18-year period (January 2003 to December 2020), this study performs a fine-scale comparison of the climatological tropopause features between these two sites that are characterized by different climates. The static stability increases rapidly above the tropopause at both sites, indicating the widespread existence of a tropopause inversion layer. The structures of both the tropopause inversion layer and the stability transition layer are more obvious at NSA than at SGP, and the seasonal variation trends of the tropopause inversion layer and stability transition layer are distinctly different between the two sites. A fitting method was used to derive the fitted tropopause height and tropopause sharpness (λ). Although this fitting method may determine a secondary tropopause rather than the primary tropopause when multiple tropopause heights are identified on one radiosonde profile, the fitted tropopause heights generally agree well with the observed tropopause heights. Broad tropopause sharpness values (λ > 2 km) occur more frequently at SGP than at NSA, resulting in a greater average tropopause sharpness at SGP (1.0 km) than at NSA (0.6 km). Significant positive trends are exhibited by the tropopause heights over the two sites, with rates of increase of 23.7 ± 6.5 m yr ⁻1 at SGP and 28.0 ± 4.0 m yr ⁻1 at NSA during the study period.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Modeling of transient conduction in building envelope assemblies: A review

As buildings age, retrofits are becoming an increasingly important topic for the ever-growing and aging existing building stock. To compare designs or evaluate in-service building envelopes, thermal modeling is utilized to evaluate the thermal performance of envelope assemblies; however, it can be difficult to model the thermal performance of as-built assemblies due to degradation or missing documentation. To address this issue, inverse modeling can be applied to infer the properties of as-built envelope assemblies. This paper presents a review of the literature and published methods to model and infer the transient conductive performance of building envelopes. This review serves as a survey of existing transient conduction algorithms to evaluate performance, computational speed, and relevance for inverse modeling applications. In addition to the literature review, this work also evaluates the computational performance of the most prevalent transient conduction algorithms against the ASHRAE 1052RP toolkit to assess inverse modeling potential. This methodology serves as a foundation for future research to characterize the transient thermal performance of as-built building envelope assemblies.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗