Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “frequency quantification”

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 73 records · Page 4

Automated ICRF heating surrogate modeling via machine learning

This work introduces automated machine learning workflows that address critical bottlenecks in surrogate model development for Ion Cyclotron Range of Frequencies (ICRF) heating applications. The automated framework includes data analysis tools that transform raw datasets into actionable insights in seconds, replacing weeks of manual exploratory effort and ensuring consistent, reproducible dataset characterization. By integrating advanced hyperparameter optimization (HPO) methods including Bayesian optimization via BoTorch and Tree-structured Parzen Estimators (TPE), the framework significantly reduces model development time from weeks to hours, decreasing computational cost and required expertise, while enabling high-accuracy surrogate models. Compared to traditional hyperparameter scanning (HPS) techniques such as methodical, randomized, and grid searches, HPO methods achieve superior convergence and predictive performance, even when compared to already well-tuned reference models. On NSTX High Harmonic Fast Wave (HHFW) heating datasets, both Random Forest Regressor (RFR) and neural network surrogates demonstrate improved accuracy, achieving R 2 values beyond 0.97 and 0.98, respectively. The results show that while HPO gains are modest for robust architectures like RFR, they become essential for more sensitive models such as neural networks, highlighting the trade-offs across optimization strategies. Through automated workflows that eliminate manual hyperparameter tuning and require minimal ML expertise, this work enables widespread adoption of high-fidelity surrogate models across the fusion community for real-time plasma control, uncertainty quantification, rapid experimental scenario development, and integrated system optimization.

Sanchez-Villar, Alvaro [Princeton Plasma Physics L↗

Guided Wave Delamination Detection and Quantification With Wavefield Data Analysis

Unexpected damage can occur in aerospace composites due to impact events or material stress during off-nominal loading events. In particular, laminated composites are susceptible to delamination damage due to weak transverse tensile and inter-laminar shear strengths. Developments of reliable and quantitative techniques to detect delamination damage in laminated composites are imperative for safe and functional optimally-designed next-generation composite structures. In this paper, we investigate guided wave interactions with delamination damage and develop quantification algorithms by using wavefield data analysis. The trapped guided waves in the delamination region are observed from the wavefield data and further quantitatively interpreted by using different wavenumber analysis methods. The frequency-wavenumber representation of the wavefield shows that new wavenumbers are present and correlate to trapped waves in the damage region. These new wavenumbers are used to detect and quantify the delamination damage through the wavenumber analysis, which can show how the wavenumber changes as a function of wave propagation distance. The location and spatial duration of the new wavenumbers can be identified, providing a useful means not only for detecting the presence of delamination damage but also allowing for estimation of the delamination size. Our method has been applied to detect and quantify real delamination damage with complex geometry (grown using a quasi-static indentation technique). The detection and quantification results show the location, size, and shape of the delamination damage.

Tian, Zhenhua↗

Estimation of Time-Varying Pilot Model Parameters

Human control behavior is rarely completely stationary over time due to fatigue or loss of attention. In addition, there are many control tasks for which human operators need to adapt their control strategy to vehicle dynamics that vary in time. In previous studies on the identification of time-varying pilot control behavior wavelets were used to estimate the time-varying frequency response functions. However, the estimation of time-varying pilot model parameters was not considered. Estimating these parameters can be a valuable tool for the quantification of different aspects of human time-varying manual control. This paper presents two methods for the estimation of time-varying pilot model parameters, a two-step method using wavelets and a windowed maximum likelihood estimation method. The methods are evaluated using simulations of a closed-loop control task with time-varying pilot equalization and vehicle dynamics. Simulations are performed with and without remnant. Both methods give accurate results when no pilot remnant is present. The wavelet transform is very sensitive to measurement noise, resulting in inaccurate parameter estimates when considerable pilot remnant is present. Maximum likelihood estimation is less sensitive to pilot remnant, but cannot detect fast changes in pilot control behavior.

Zaal, Peter M. T.↗

A Model-Based, Bayesian Solution for Characterization of Complex Damage Scenarios in Aerospace Composite Structures

Ultrasonic damage detection and characterization is commonly used in nondestructive evaluation (NDE) of aerospace composite components. In recent years there has been an increased development of guided wave based methods. In real materials and structures, these dispersive waves result in complicated behavior in the presence of complex damage scenarios. Model-based characterization methods utilize accurate three dimensional finite element models (FEMs) of guided wave interaction with realistic damage scenarios to aid in defect identification and classification. This work describes an inverse solution for realistic composite damage characterization by comparing the wavenumber-frequency spectra of experimental and simulated ultrasonic inspections. The composite laminate material properties are first verified through a Bayesian solution (Markov chain Monte Carlo), enabling uncertainty quantification surrounding the characterization. A study is undertaken to assess the efficacy of the proposed damage model and comparative metrics between the experimental and simulated output. The FEM is then parameterized with a damage model capable of describing the typical complex damage created by impact events in composites. The damage is characterized through a transdimensional Markov chain Monte Carlo solution, enabling a flexible damage model capable of adapting to the complex damage geometry investigated here. The posterior probability distributions of the individual delamination petals as well as the overall envelope of the damage site are determined.

H. Reed↗

Satellite remote sensing of drought conditions

Multitemporal satellite data have application in the detection and quantification of drought through the ability of these data to estimate the photosynthetic capacity of the terrestrial surface and record microwave surface brightness at the 37 GHz frequency. With proper calibration and registration, comparisons can be made between and among years for specific months using the photosynthetic capacity and the 37 GHz microwave surface brightness for selected time periods or growing seasons. This technology has application in identifying and quantifying areas experiencing drought.

Tucker, Compton J.↗

Robust trajectory-constrained frequency control for microgrids considering model linearization error

Grid supportive modes integrated within inverter-based resources can improve the frequency response of renewable-rich microgrids. The synthesis of grid supportive modes to guarantee frequency trajectory constraints under a predefined disturbance set is challenging but essential. To tackle this challenge, a numerical optimal control (NOC)-based control synthesis methodology is proposed. Without loss of generality, a wind-diesel fed microgrid is studied, where we aim to design grid supportive functions in the wind turbine. In the control design, linearized models are used, and the linearization-induced errors are quantitatively analyzed by reachability and interval arithmetics and represented in the form of interval uncertainties. Then, the NOC problem can be formulated into a robust mixed-integer linear program. The control structure is strategically configured into two levels to realize online deployment. The proposed control is verified on the modified 33-node microgrid with a full-order three-phase nonlinear model in Simulink. In conclusion, the simulation results show the effectiveness of the proposed control paradigm and the necessity of considering linearization-induced uncertainty.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Gravity Wave Forcing of the Mesosphere and Lower Thermosphere: Mountain and Convective Waves Ascending Vertically (MaCWAVE)

The specific objectives of this research effort included the following: 1) Quantification of gravity wave propagation throughout the lower and middle atmosphere in order to define the roles of topographic and convective sources and filtering by mean and low-frequency winds in defining the wave field and wave fluxes at greater altitudes; 2) The influences of wave instability processes in constraining wave amplitudes and fluxes and generating turbulence and transport; 3) Gravity wave forcing of the mean circulation and thermal structure in the presence of variable motion fields and wave-wave interactions, since the mean forcing may be a small residual when wave interactions, anisotropy, and momentum and heat fluxes are large; 4) The statistical forcing and variability imposed on the thermosphere at greater altitudes by the strong wave forcing and interactions occurring in the MLTI.

Fritts, David C.↗

Vertical column dual-comb spectroscopy to a TBS

Open-path dual-frequency-comb spectroscopy (DCS) is a broadband, high spectral resolution, and high precision method for measuring gas concentrations over kilometer-scale paths. It has been used for detection and quantification of emissions of pollutants, hazardous gasses, and greenhouse gasses (GHGs). DCS has been shown to measure trace-gas mixing ratios with 0.14%-0.4% agreement between instruments. To achieve high signal-to-noise ratios (SNRs) with DCS, comb light is targeted onto a retroreflector at the end of the measurement path, which returns the signal to a detector co-located with the launch signal. Installing the retroreflector on a mobile platform such as a balloon or unmanned aerial vehicle (UAV) extends the capabilities of DCS by enabling variable path lengths, greater mobility, and access to higher altitudes. Mobile-target DCS has many uses in plume and leak detection, emissions modeling, and planetary boundary layer (PBL) studies. It is a promising method for observing vertical distributions of GHGs and mixing processes in the PBL, which are difficult to measure but important for pollution and climate monitoring as well as for understanding transport of gasses through the atmosphere. Tethered balloons are an intriguing platform because they enable longer flight durations and higher altitudes than easily obtainable with a UAV and thus allow for column measurements up to and above the PBL. New measurements completed in October/November 2024 reach the highest altitudes above ground level yet achieved by mobile-target DCS. For these measurements, the retroreflector was mounted on a 7 m diameter tethered helium balloon. An actively tracking gimbal on the ground holds the DCS launch telescope and keeps the 5 cm beam pointed onto the retroreflector while the balloon is lifted, lowered, and moved by wind and turbulence.

atmosphere↗

Crack Imaging and Quantification in Aluminum Plates with Guided Wave Wavenumber Analysis Methods

Guided wavefield analysis methods for detection and quantification of crack damage in an aluminum plate are presented in this paper. New wavenumber components created by abrupt wave changes at the structural discontinuity are identified in the frequency-wavenumber spectra. It is shown that the new wavenumbers can be used to detect and characterize the crack dimensions. Two imaging based approaches, filter reconstructed imaging and spatial wavenumber imaging, are used to demonstrate how the cracks can be evaluated with wavenumber analysis. The filter reconstructed imaging is shown to be a rapid method to map the plate and any existing damage, but with less precision in estimating crack dimensions; while the spatial wavenumber imaging provides an intensity image of spatial wavenumber values with enhanced resolution of crack dimensions. These techniques are applied to simulated wavefield data, and the simulation based studies show that spatial wavenumber imaging method is able to distinguish cracks of different severities. Laboratory experimental validation is performed for a single crack case to confirm the methods' capabilities for imaging cracks in plates.

Yu, Lingyu↗

Piezoresponse amplitude and phase quantified for electromechanical characterization

Piezoresponse force microscopy (PFM) is a powerful characterization technique to readily image and manipulate the ferroelectric domains. PFM gives an insight into the strength of local piezoelectric coupling and polarization direction through PFM amplitude and phase, respectively. Converting measured arbitrary units into units of effective piezoelectric constant remains a challenge, and insufficient methods are often used. While most quantification efforts have been spent on quantifying the PFM amplitude signal, little attention has been given to the PFM phase, which is often arbitrarily adjusted to fit expectations. This is problematic when investigating materials with unknown or negative sign of the probed effective electrostrictive coefficient or strong frequency dispersion of electromechanical responses, because assumptions about the PFM phase cannot be reliably made. The PFM phase can, however, provide important information on the polarization orientation and the sign of the effective electrostrictive coefficient probed by PFM. Most notably, the orientation of the PFM hysteresis loop is determined by the PFM phase. Moreover, when presenting PFM data as a combined signal, the resulting response can be artificially lowered or asymmetric if the phase data have not been correctly processed. Here, we explain the PFM amplitude quantification process and demonstrate a path to identify the phase offset required to extract correct meaning from the PFM phase data. We explore different sources of phase offsets including the experimental setup, instrumental contributions, and data analysis. We discuss the physical working principles of PFM and develop a strategy to extract physical meaning from the PFM amplitude and phase.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Quantification of Dissolved O 2 in Bulk Aqueous Solutions and Porous Media Using NMR Relaxometry

Effects of dissolved paramagnetic oxygen (O 2 ) in water on 1 H nuclear magnetic resonance (NMR) Carr-Purcell-Meiboom-Gill (CPMG) experiments is evaluated at a 1 H Larmor frequency of 2 MHz. Dissolution of O 2 into water significantly reduces the 1 H transverse relaxation coefficient (T 2 ). For deoxygenated water, T 2 is 3388 ms, water at ambient atmospheric conditions (7.4 mg/L O2) exhibits a T2 of 2465 ms, and dissolution of 2710 mg/L O 2 further reduces T 2 to 36 ms. The results were fit with an empirical model to facilitate prediction of T 2 times for bulk water as a function of paramagnetic oxygen concentrations in solution. Dissolved O 2 also greatly influences 1 H NMR CPMG experiments of confined water in a model system composed of Berea sandstone. For this system, 90 mg/L O 2 in H 2 O enhances T 2 relaxation of bulk water such that the relaxation time is comparable to physically confined water in the sandstone pores. Given the sensitivity of NMR T 2 coefficients to paramagnetic oxygen, low-field NMR-based characterization of fluid and porous media structure requires control of dissolved oxygen, as geospatial variation in the partial pressure of O 2 alone is expected to perturb fluid and pore relaxation times by up to 60 and 36%, respectively.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Dispersion curve regression – assisted wideband local wavenumber analysis for characterizing three-dimensional (3D) profile of hidden corrosion damage

Corrosion as common damage in civil, petrochemical, nuclear, and aerospace structures affects the integrity and safety of the structures and may lead to catastrophic failures. This paper presents dispersion curve regression-assisted local wavenumber analysis method, which can analyze the time–space wavefield containing wideband information of wave-damage interaction and further extract the structural information carried by such wavefield for characterizing hidden corrosion damage in an isotropic plate. To acquire the time–space wavefield, a noncontact scanning laser Doppler vibrometer is used. In our analysis method, the acquired time–space wavefield is firstly processed to generate a local wavenumber function. The derived wavenumber function is further analyzed with a dispersion curve regression step, which searches in a set of theoretical frequency-wavenumber dispersion curves for different plate thicknesses and identifies the optimal plate thickness whose theoretical curve best matches the frequency-wavenumber relation contained in the local wavenumber function. By this means, a 3D profile including both in plane and thickness dimensions of the structure can be constructed for corrosion visualization and quantification. The experimental study demonstrates that our method can quantify the profile of hidden uniform corrosion in metal plates in terms of its in-plane shape and size as well as its out-of-plane depth in the subwavelength scale. Moreover, compared to previous Lamb wave-based corrosion inspection methods, our technique allows for noncontact 3D characterization of hidden corrosion from a far distance to the structure. Finally, we believe this work will inspire new 3D damage quantification methods that are based on wavefield analysis as well as enable potential applications for the quantification of hidden corrosion in civil, petrochemical, nuclear, and aerospace structures.

42 ENGINEERING↗

Data-driven gradient optimization for field emission management in a superconducting radio-frequency linac

Field emission can cause significant problems in superconducting radio-frequency linear accelerators (linacs). When cavity gradients are pushed higher, radiation levels within the linacs may rise exponentially, causing degradation of many nearby systems. This research aims to utilize machine learning with uncertainty quantification to predict radiation levels at multiple locations throughout the linacs and ultimately optimize cavity gradients to reduce field emission-induced radiation while maintaining the total linac energy gain necessary for the experimental physics program. The optimized solutions show over 40% reductions for both neutron and gamma radiation from the standard operational settings. Published by the American Physical Society 2025

43 PARTICLE ACCELERATORS↗

Impedance Response Influenced by Variability in the Random Distribution of Physical Properties of Coated Materials in Two-Dimensional Space

Heterogeneous physical characteristics of a system featuring a single-layer film on a metallic surface have been explored via its impedance response. The Nyquist plot showed a distorted semicircle, indicative of the system's unique distribution characteristics. Utilizing a copula-based probability method, a two-dimensional deterministic impedance model was successfully integrated, accounting for spatial physical properties such as permittivity and electrical conductivity. This strategy enabled in-depth exploration and mechanistic quantification of a broad spectrum of properties. A quantitative understanding of impedance signal alterations, characterized by normally or log-normally correlated variables, was achieved through the variation in aspect ratio and characteristic frequency of the impedance spectra. Log-normally distributed electrical properties provided a superior representation of the distorted impedance spectra. As coefficient of variation (CV) values fluctuated, the aspect ratio and characteristic frequency showed heightened sensitivity to log-normal permittivity compared to log-normal electrical conductivity. Notably, a marked positive linear correlation between electrical properties resulted in an impedance response that approximated perfect semicircular spectra. The variability in the electrical properties' distribution was demonstrated by considering the correlation coefficient between electrical conductivity and the z-direction position. Furthermore, the highest aspect ratio of the impedance spectra was observed when the electrical conductivity was randomly distributed across the z-direction space.

36 MATERIALS SCIENCE↗

Improving the Quasi‐Biennial Oscillation via a Surrogate‐Accelerated Multi‐Objective Optimization

Accurate simulation of the quasi-biennial oscillation (QBO) is challenging due to uncertainties in representing convectively generated gravity waves. We develop an end-to-end uncertainty quantification workflow that calibrates these gravity wave processes in E3SM for a realistic QBO. Central to our approach is a domain knowledge-informed, compressed representation of high-dimensional spatio-temporal wind fields. By employing a parsimonious statistical model that learns the fundamental frequency from complex observations, we extract interpretable and physically meaningful quantities capturing key attributes. Building on this, we train a probabilistic surrogate model that approximates the fundamental characteristics of the QBO as functions of critical physics parameters governing gravity wave generation. Leveraging the Karhunen–Loève decomposition, our surrogate efficiently represents these characteristics as a set of orthogonal features, capturing cross-correlations among multiple physics quantities evaluated at different pressure levels and enabling rapid surrogate-based inference at a fraction of the computational cost of full-scale simulations. Finally, we analyze the inverse problem using a multi-objective approach. Our study reveals a tension between amplitude and period that constrains the QBO representation, precluding a single optimal solution. To navigate this, we quantify the bi-criteria trade-off and generate a set of Pareto optimal parameter values that balance the conflicting objectives. This integrated workflow improves the fidelity of QBO simulations and offers a versatile template for uncertainty quantification in complex geophysical models.

54 ENVIRONMENTAL SCIENCES↗

Quantification and prediction of error in MIMO tests

Dynamic testing of systems is most realistic of real-world conditions when multiple input and multiple output (MIMO) techniques are used. To replicate measured environmental conditions, a series of desired outputs (responses) on the system must be realized by inverting the frequency response functions (FRF) matrix for input estimation. System identification is dependent on the number of inputs and outputs, which affect the type of solutions. In input estimation, the auto spectra of the determined inputs are affected by both the auto spectra and cross spectra of the desired outputs. This research evaluates errors and investigates two sources in an experimental setting. The first significant error source is test-to-test variability in FRFs from system identification tests, quantified by magnitude variability from long-duration tests. The second significant error source is the realization of time histories from frequency domain inputs. The spectral content of the realized random process deviates from the desired content and is also quantified. The structure is a linear three-story frame fixed at the base; two inputs are provided uniaxially with the use of suspended electrodynamic shakers. The proposed model for error sources was found to be effective in predicting errors in a SISO (single-input single-output) and a square (2-input 2-output) MIMO test.

Yousef, Odey↗

Dual-Comb Spectroscopy to Tethered Balloon System (DCS to TBS) Field Campaign Report

There is a critical need for the ability to detect, quantify, and attribute sources and sinks of greenhouse gases (GHGs) with high temporal and spatial resolution. Current methods are typically based on continuously operating near-surface point sensors or intermittent aerial/satellite observations. Dual optical frequency comb spectroscopy (“dual-comb spectroscopy”; DCS) is a laser-based technique that can bridge this sensing gap by providing high-precision, high-accuracy, and high-temporal-resolution detection and quantification of trace gases over kilometer-length scales, connecting local and regional observations. A single DCS system can also provide spatial resolution for source/sink attribution with the addition of a few minor optical elements, enabling cross-validation of other observations (e.g., Total Carbon Column Observing Network [TCCON], Orbiting Carbon Observatory [OCO-2/3], Greenhouse Gases Observing Satellite [GOSAT-1/2][).

54 ENVIRONMENTAL SCIENCES↗

A Probabilistic Asteroid Impact Risk Model

Asteroid threat assessment requires the quantification of both the impact likelihood and resulting consequence across the range of possible events. This paper presents a probabilistic asteroid impact risk (PAIR) assessment model developed for this purpose. The model incorporates published impact frequency rates with state-of-the-art consequence assessment tools, applied within a Monte Carlo framework that generates sets of impact scenarios from uncertain parameter distributions. Explicit treatment of atmospheric entry is included to produce energy deposition rates that account for the effects of thermal ablation and object fragmentation. These energy deposition rates are used to model the resulting ground damage, and affected populations are computed for the sampled impact locations. The results for each scenario are aggregated into a distribution of potential outcomes that reflect the range of uncertain impact parameters, population densities, and strike probabilities. As an illustration of the utility of the PAIR model, the results are used to address the question of what minimum size asteroid constitutes a threat to the population. To answer this question, complete distributions of results are combined with a hypothetical risk tolerance posture to provide the minimum size, given sets of initial assumptions. Model outputs demonstrate how such questions can be answered and provide a means for interpreting the effect that input assumptions and uncertainty can have on final risk-based decisions. Model results can be used to prioritize investments to gain knowledge in critical areas or, conversely, to identify areas where additional data has little effect on the metrics of interest.

Asteroid↗