Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “modal decomposition”

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.

122 records · Page 7

Sparsity of Radiating Characteristic Modes on Infinite Periodic Structures

Characteristic modes on infinite periodic structures are studied using spectral dyadic Green’s functions. This formulation demonstrates that, in contrast to the modal analysis of finite structures, the number of radiating characteristic modes is limited by unit cell size and incident wave vector (i.e., scan angle or phase shift per unit cell). Here, the reflection tensor is decomposed into modal contributions from radiating modes, indicating that characteristic modes are a predictably sparse basis in which to study reflection phenomena.

42 ENGINEERING↗

Model-Agnostic Algorithm for Real-Time Attack Identification in Power Grid using Koopman Modes

Malicious activities on measurements from sensors like Phasor Measurement Units (PMUs) can mislead the control center operator into taking wrong control actions resulting in disruption of operation, financial losses, and equipment damage. In particular, false data attacks initiated during power systems transients caused due to abrupt changes in load and generation can fool the conventional model-based detection methods relying on thresholds comparison to trigger an anomaly. In this paper, we propose a Koopman mode decomposition (KMD) based algorithm to detect and identify false data attacks in real-time. The Koopman modes (KMs) are capable of capturing the nonlinear modes of oscillation in the transient dynamics of the power networks and reveal the spatial embedding of both natural and anomalous modes of oscillations in the sensor measurements. The Koopman-based spatio-temporal nonlinear modal analysis is used to filter out the false data injected by an attacker. The performance of the algorithm is illustrated on the IEEE 68-bus test system using synthetic attack scenarios generated on GridSTAGE, a recently developed multivariate spatio-temporal data generation framework for simulation of adversarial scenarios in cyber-physical power systems.

Nandanoori, Sai Pushpak↗

Control of elastic robotic systems by nonlinear inversion and modal damping

Energy efficient, lightweight robot arms for space applications have considerable structural flexibility. For large and fast motions, both the nonlinear coupled dynamics and the elastic behavior of the robots must be considered in control system designs. This paper presents an approach to the control of a class of flexible robotic systems. A control law is derived which decouples the joint-angle motion from the flexible motion and asymptotically decomposes the elastic dynamics into two subsystems, representing the transverse vibrations of the elastic link in two orthogonal planes. This decomposition allows the design of an elastic mode stabilizer independently based on lower order models representing structural flexibility. The closed-loop system is shown to be globally asymptotically stable and robust to uncertainty in system parameters. Simulation results are presented to show that large, fast control of joint angles can be performed in spite of space vehicle motion and uncertainty in the payload.

Singh, S. N.↗

Non-Linear Structural Dynamics Characterization using a Scanning Laser Vibrometer

This paper presents the use of a scanning laser vibrometer and a signal decomposition method to characterize non-linear dynamics of highly flexible structures. A Polytec PI PSV-200 scanning laser vibrometer is used to measure transverse velocities of points on a structure subjected to a harmonic excitation. Velocity profiles at different times are constructed using the measured velocities, and then each velocity profile is decomposed using the first four linear mode shapes and a least-squares curve-fitting method. From the variations of the obtained modal \ielocities with time we search for possible non-linear phenomena. A cantilevered titanium alloy beam subjected to harmonic base-excitations around the second. third, and fourth natural frequencies are examined in detail. Influences of the fixture mass. gravity. mass centers of mode shapes. and non-linearities are evaluated. Geometrically exact equations governing the planar, harmonic large-amplitude vibrations of beams are solved for operational deflection shapes using the multiple shooting method. Experimental results show the existence of 1:3 and 1:2:3 external and internal resonances. energy transfer from high-frequency modes to the first mode. and amplitude- and phase- modulation among several modes. Moreover, the existence of non-linear normal modes is found to be questionable.

Pai, P. F.↗

Acoustic Mode Decomposition in Rectangular Ducts with Sheared Flow

The performance of new acoustic liner concepts are, in general, characterized and assessed in grazing flow rigs early in the development cycle. These test rigs expose an acoustic liner sample, installed on the side wall of the duct, to a grazing flow and incident acoustic field. The process to characterize these liners involves educing the impedance on the wall where the sample is installed and examining the acoustic power attenuation. Standard approaches to computing impedance or power attenuation generally consider only the effects of a 2D shear flow or uniform flow, on the acoustic field. In this study, the objective is to incorporate 3D shear flow effects in the analysis of acoustic mode attenuation in a rectangular duct flow rig. A modal analysis of microphone measurements obtained on the side walls of the duct upstream and downstream of the test section of the rig is developed. A Galerkin projection of the Pridmore-Brown equation is performed with Chebyshev basis functions in order to incorporate the effects of the Mach number profile on the computation of the axial wavenumber of each mode. Measurements of the Mach number profile are obtained in the test rig and used as input to compute the modes. Comparisons made between the sound field computed with traditional convective Helmholtz modes and the new procedure using Pridmore-Brown modes indicate that the computed acoustic field using Pridmore-Brown modes more accurately reconstructs the acoustic signal at each microphone in the array. The mode structure of the lowest-order mode is shown to be significantly impacted by shear flow refraction effects, and higher-order mode structures are also affected at higher frequencies and centerline Mach number. An assessment of the acoustic mode attenuation for two acoustic liner samples demonstrates that the computed mode amplitudes for both the traditional and new approach are in agreement for the lowest-order mode, but discrepancies arise when higher-order modes are the dominant component of the acoustic field.

Duct acoustics↗

Acoustic Mode Decomposition in Rectangular Ducts with Sheared Flow

The performance of new acoustic liner concepts are, in general, characterized and assessed in grazing flow rigs early in the development cycle. These test rigs expose an acoustic liner sample, installed on the side wall of the duct, to a grazing flow and incident acoustic field. The process to characterize these liners involves educing the impedance on the wall where the sample is installed and examining the acoustic power attenuation. Standard approaches to computing impedance or power attenuation generally consider only the effects of a 2D shear flow or uniform flow, on the acoustic field. In this study, the objective is to incorporate 3D shear flow effects in the analysis of acoustic mode attenuation in a rectangular duct flow rig. A modal analysis of microphone measurements obtained on the side walls of the duct upstream and downstream of the test section of the rig is developed. A Galerkin projection of the Pridmore-Brown equation is performed with Chebyshev basis functions in order to incorporate the effects of the Mach number profile on the computation of the axial wavenumber of each mode. Measurements of the Mach number profile are obtained in the test rig and used as input to compute the modes. Comparisons made between the sound field computed with traditional convective Helmholtz modes and the new procedure using Pridmore-Brown modes indicate that the computed acoustic field using Pridmore-Brown modes more accurately reconstruct the acoustic signals at each microphone in the array. The mode structure of the lowest order mode is shown to be significantly impacted by shear flow refraction effects, and higher order mode structures are also effected at higher frequencies and centerline Mach number. An assessment of the acoustic mode attenuation for two acoustic liner samples demonstrates that the computed mode amplitudes for both the traditional and new approach are in agreement for the lowest order mode, but discrepancies arise when higher order modes are the dominant component of the acoustic field.

Duct acoustics↗

Chasing the multi-modal plasma response in MAST-U

Abstract Achieving edge localized modes (ELMs) suppression in spherical tokamaks by applying resonant magnetic perturbations (RMPs) has proven challenging. The poloidal spectrum of the applied RMP is a key parameter that has an impact on the capability to mitigate and eventually suppress ELMs. In this work the resistive magnetohydrodynamic code MARS-F (Liu et al 2000 Phys. Plasmas 7 3681) is used to evaluate the possibility of directly measuring the plasma response in MAST-U, and particularly its variation as function of the applied poloidal spectrum, in order to guide the experimental validation of the predicted best RMP configuration for ELM suppression. Toroidal mode number n = 2 RMP is considered to minimize the presence of sidebands, and to avoid the deleterious core coupling of n = 1. Singular Value Decomposition is used to highlight linearly independent structures in the simulated magnetic 3D fields and how those structures can be measured at the wall where the magnetic sensors are located. Alternative ways to measure the multimodal plasma response and how they can be used to infer the best RMP configuration to achieve ELM suppression are also presented, including the plasma displacement and the 3D footprints at the divertor plates.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Trend and Variability of China Precipitation in Spring and Summer: Linkage to Sea Surface Temperatures

Observational records in the past 50 years show an upward trend of boreal-summer precipitation over central eastern China and a downward trend over northern China. During boreal spring, the trend is upward over southeastern China and downward over central eastern China. This study explores the forcing mechanism of these trends in association with the global sea-surface temperature (SST) variations on the interannual and inter-decadal timescales. Results based on Singular Value Decomposition analyses (SVD) show that the interannual variability of China precipitation in boreal spring and summer can be well defined by two centers of actions for each season, which are co-varying with two interannual modes of SSTs. The first SVD modes of precipitation in spring and summer, which are centered in southeastern China and northern China, respectively, are linked to an ENSO-like mode of SSTs. The second SVD modes of precipitation in both seasons are confined to central eastern China, and are primarily linked to SST variations over the warm pool and Indian Ocean. Features of the anomalous 850-hPa winds and 700-Wa geopotential height corresponding to these modes support a physical mechanism that explains the causal links between the modal variations of precipitation and SSTs. On the decadal and longer timescale, similar causal links are found between the same modes of precipitation and SSTs, except for the case of springtime precipitation over central eastern China. For this case, while the interannual mode of precipitation is positively correlated with the interannual variations of SSTs over the warm pool and Indian Ocean; the inter-decadal mode is negatively correlated with a different SST mode, the North Pacific mode. The later is responsible for the observed downward trend of springtime precipitation over central eastern China. For all other cases, both the interannual and inter-decadal variations of precipitation can be explained by the same mode of SSTs. The upward trend of springtime precipitation over southeastern China and downward trend of summertime precipitation over northern China are attributable to the warming trend of the ENSO-like mode. The recent frequent summertime floods over central eastern China are linked to the warming trend of SSTs over the warm pool and Indian Ocean.

Yang, Fanglin↗

Aeroelastic Flight Data Analysis with the Hilbert-Huang Algorithm

This report investigates the utility of the Hilbert Huang transform for the analysis of aeroelastic flight data. It is well known that the classical Hilbert transform can be used for time-frequency analysis of functions or signals. Unfortunately, the Hilbert transform can only be effectively applied to an extremely small class of signals, namely those that are characterized by a single frequency component at any instant in time. The recently-developed Hilbert Huang algorithm addresses the limitations of the classical Hilbert transform through a process known as empirical mode decomposition. Using this approach, the data is filtered into a series of intrinsic mode functions, each of which admits a well-behaved Hilbert transform. In this manner, the Hilbert Huang algorithm affords time-frequency analysis of a large class of signals. This powerful tool has been applied in the analysis of scientific data, structural system identification, mechanical system fault detection, and even image processing. The purpose of this report is to demonstrate the potential applications of the Hilbert Huang algorithm for the analysis of aeroelastic systems, with improvements such as localized online processing. Applications for correlations between system input and output, and amongst output sensors, are discussed to characterize the time-varying amplitude and frequency correlations present in the various components of multiple data channels. Online stability analyses and modal identification are also presented. Examples are given using aeroelastic test data from the F-18 Active Aeroelastic Wing airplane, an Aerostructures Test Wing, and pitch plunge simulation.

Brenner, Martin J.↗

Aeroelastic Flight Data Analysis with the Hilbert-Huang Algorithm

This paper investigates the utility of the Hilbert-Huang transform for the analysis of aeroelastic flight data. It is well known that the classical Hilbert transform can be used for time-frequency analysis of functions or signals. Unfortunately, the Hilbert transform can only be effectively applied to an extremely small class of signals, namely those that are characterized by a single frequency component at any instant in time. The recently-developed Hilbert-Huang algorithm addresses the limitations of the classical Hilbert transform through a process known as empirical mode decomposition. Using this approach, the data is filtered into a series of intrinsic mode functions, each of which admits a well-behaved Hilbert transform. In this manner, the Hilbert-Huang algorithm affords time-frequency analysis of a large class of signals. This powerful tool has been applied in the analysis of scientific data, structural system identification, mechanical system fault detection, and even image processing. The purpose of this paper is to demonstrate the potential applications of the Hilbert-Huang algorithm for the analysis of aeroelastic systems, with improvements such as localized/online processing. Applications for correlations between system input and output, and amongst output sensors, are discussed to characterize the time-varying amplitude and frequency correlations present in the various components of multiple data channels. Online stability analyses and modal identification are also presented. Examples are given using aeroelastic test data from the F/A-18 Active Aeroelastic Wing aircraft, an Aerostructures Test Wing, and pitch-plunge simulation.

Brenner, Marty↗

Agile Multi-Scale Decompositions for Automatic Image Registration

In recent works, the first and third authors developed an automatic image registration algorithm based on a multiscale hybrid image decomposition with anisotropic shearlets and isotropic wavelets. This prototype showed strong performance, improving robustness over registration with wavelets alone. However, this method imposed a strict hierarchy on the order in which shearlet and wavelet features were used in the registration process, and also involved an unintegrated mixture of MATLAB and C code. In this paper, we introduce a more agile model for generating features, in which a flexible and user-guided mix of shearlet and wavelet features are computed. Compared to the previous prototype, this method introduces a flexibility to the order in which shearlet and wavelet features are used in the registration process. Moreover, the present algorithm is now fully coded in C, making it more efficient and portable than the MATLAB and C prototype. We demonstrate the versatility and computational efficiency of this approach by performing registration experiments with the fully-integrated C algorithm. In particular, meaningful timing studies can now be performed, to give a concrete analysis of the computational costs of the flexible feature extraction. Examples of synthetically warped and real multi-modal images are analyzed.

Science Data Processing↗

Agile Multi-Scale Decompositions for Automatic Image Registration

In recent works, the first and third authors developed an automatic image registration algorithm based on a multiscale hybrid image decomposition with anisotropic shearlets and isotropic wavelets. This prototype showed strong performance, improving robustness over registration with wavelets alone. However, this method imposed a strict hierarchy on the order in which shearlet and wavelet features were used in the registration process, and also involved an unintegrated mixture of MATLAB and C code. In this paper, we introduce a more agile model for generating features, in which a flexible and user-guided mix of shearlet and wavelet features are computed. Compared to the previous prototype, this method introduces a flexibility to the order in which shearlet and wavelet features are used in the registration process. Moreover, the present algorithm is now fully coded in C, making it more efficient and portable than the MATLAB and C prototype. We demonstrate the versatility and computational efficiency of this approach by performing registration experiments with the fully-integrated C algorithm. In particular, meaningful timing studies can now be performed, to give a concrete analysis of the computational costs of the flexible feature extraction. Examples of synthetically warped and real multi-modal images are analyzed.

Science Data Processing↗

Striving to translate shale physics across ten orders of magnitude: What have we learned?

Shales will play an important role in the successful transition of energy from fossil-based resources to renewables in the coming decades. Aside from being a significant source of low-carbon intensity fuels, like natural gas, they also serve as geologic seals of subsurface formations that may be used to isolate nuclear waste, sequester CO 2 , or store intermittent energy (e.g., solar hydrogen). Despite their importance, shales pose significant engineering and environmental challenges due to their nanoporous structure and extreme heterogeneity that spans at least ~10 orders of magnitude in spatial scale. Two challenges inhibit a system-level understanding: (1) the physics of fluid flow and phase behavior in shales are poorly understood due to the dominant molecular interactions between minerals and fluids under confinement, and (2) the apparent lack of scale separation that prevents a reliable (closed) description of the physics at any single scale of observation. In this review, we focus on the latter issue and discuss scale translation, which in its broadest sense is transforming data or simulations from one spatiotemporal scale to another. While effective scale translation is not exclusive to shales, but all geologic porous media, the need for it is especially acute in shales given their high degree of heterogeneity. Classical theories like homogenization, while indispensable, fail when scales are not separated. Other methods, like numerical upscaling, scale-translate in only one direction: small to large, but not the reverse, called downscaling. However, the confluence of advances in three areas are bringing challenging problems such as shales within reach: increased computational power and scalable algorithms; high-resolution imaging and multi-modal data acquisition; and machine learning to process massive amounts of data. While these advances equip geoscientists with a wide array of experimental and computational tools, no individual tool can probe the entire gamut of heterogeneity in shales. Their effective use, therefore, requires an ability to bridge between various data types obtained at different scales. The aim of this review is to present a coherent account of computational and experimental methods that may be used to achieve just that, i.e., to perform scale translation. We provide a broader definition of scale translation, one that transcends classical homogenization and upscaling methods, but is consistent with them and accommodates notions like downscaling and data translation. After a brief introduction to homogenization, we review hybrid methods, numerical upscaling and its recent extensions, multiscale computing, high-resolution imaging, and machine learning. We place particular emphasis on multiscale computing and propose an algorithmic framework to bridge between the pore (micro) and Darcy (macro) scales. Throughout the paper, we draw comparisons between the various methods and highlight their (often hidden) similarities, differences, benefits, and pitfalls. We finally conclude with two case studies on shales that exemplify some of the methods presented.

58 GEOSCIENCES↗

Shape-constrained Input Estimation for Efficient Multi-shaker Vibration Testing

Multi-shaker vibration testing is gaining interest from structural dynamics test engineers as it can provide a much more accurate match to complicated field vibration responses than traditional single-axis shaker tests. However, the force capabilities of the small modal shakers typically used in multi-shaker vibration tests has limited the achievable response levels. To date, most multi-shaker vibration tests have been performed using a variety of standard, commercially-available control systems. While these control systems are adequate for a wide range of multiple-input/multiple-output tests, their control algorithms have not been tailored for the specific problem of multi-shaker vibration tests: efficiently coordinating the various shakers to work together to achieve a desired response. Here, a new input estimation algorithm is developed and demonstrated using simulations and actual test data. This algorithm, dubbed shape-constrained input estimation, is shown to effectively coordinate multiple shakers using a set of constraint vectors based on the deflection shapes of the test structure. This is accomplished by using the singular vector shapes of the system frequency response matrix, which allows the constraint vectors to automatically change as a function of frequency. Simulation and test results indicate a significant reduction in the input forces required to achieve a desired response. Finally, the results indicate that shape-constrained input estimation is an effective method to achieve higher response levels from limited shaker forces which will enable higher level multi-shaker vibration tests to be performed.

42 ENGINEERING↗