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.

At least 109 records · Page 6

Development and implementation of an in situ high-resolution isotopic microscope for measuring metabolic interactions in soil mesocosms (Final Report)

Many of our planet’s ecosystems rely on the activities of soil microbial communities. These microbes have significant and wide-ranging effects: they metabolize carbon and other nutrients, interact with plants and fungi, and perform other processes important to soil health. However, our ability to directly observe the enzymatic and metabolic activities of microbes within soil is currently limited, not only by the complexity of soil microbial communities themselves, but also by the lack of experimental tools to study them and their molecular interactions in situ. These challenges hinder our understanding of the life-sustaining processes of biomass decomposition and the manner in which it contributes to the movement of freed carbon within soil ecosystems. Here we aimed to develop a novel ultrahigh-resolution isotopic microscope that combines complementary imaging modalities to gain insights into metabolic cycling in soil. This final report covers both the portion of this work that was initially completed at the University of North Carolina at Chapel Hill from 2018 to 2019 (as award DE-SC0019012) and then from 2020 to 2022 (2023 in NCE) at the University of Massachusetts Chan Medical School. This goal was to build an integrated platform consisting of fluorescence microscopy, Raman microspectroscopy, and nanospray desorption electrospray ionization Fourier transform ion cyclotron resonance mass spectrometry (nanoDESI-FTICR-MS) to directly investigate microbial activities and molecular transformations occurring in soil by exploiting the use of both fluorescent labels and stable isotope probing. Specifically, our instrument was conceived as integrating: (A) fluorescence detection to localize soil microbes, identify bacteria taking up polysaccharides, and monitor gene expression of enzymes involved in decomposition; (B) Raman microspectroscopy to determine which microbes incorporate decomposition products into their biomass, as well as which decomposers are metabolically active; and (C) nanoDESI-FTICR-MS imaging to spatially probe, in real-time, the metabolites in the surrounding area, which we expect will reveal the distributions of the products of the enzymatic breakdown of polysaccharides, as well as specialized metabolites acting as cell-cell signals between decomposers. To enable these measurements, we exploited fluorescence-based probes to map the microbes that are enzymatically active (per A) and utilized stable isotope-labeled substrates to visualize both the microbial and molecular fate of decomposed biomass (B and C, respectively). To accomplish our goal, we combined the expertise of a multidisciplinary group of scientists to pursue the construction of this microscope and to investigate scientific questions that would be facilitated by such a capability. We aimed to employ this spatially informative, high-resolution isotopic microscope to visualize the critical steps of biomass degradation and the molecular fate of other environmentally relevant substrates within soil mesocosms. Accordingly, this technology will enhance our understanding of the microbial and metabolic interactions occurring within soil communities that are relevant to carbon degradation and other soil processes. The instrument has been developed and housed at EMSL, where it will be available to the entire EMSL User Base, enabling a variety of related DOE-relevant systems to be interrogated in the future by diverse scientific research groups.

59 BASIC BIOLOGICAL SCIENCES↗

Mapping microhabitats of lignocellulose decomposition by a microbial consortium

The leaf-cutter ant fungal garden ecosystem is a naturally evolved model system for efficient plant biomass degradation. Degradation processes mediated by the symbiotic fungus Leucoagaricus gongylophorus are difficult to characterize due to dynamic metabolisms and spatial complexity of the system. Herein, we performed microscale imaging across 12-µm-thick adjacent sections of Atta cephalotes fungal gardens and applied a metabolome-informed proteome imaging approach to map lignin degradation. This approach combines two spatial multiomics mass spectrometry modalities that enabled us to visualize colocalized metabolites and proteins across and through the fungal garden. Spatially profiled metabolites revealed an accumulation of lignin-related products, outlining morphologically unique lignin microhabitats. Metaproteomic analyses of these microhabitats revealed carbohydrate-degrading enzymes, indicating a prominent fungal role in lignocellulose decomposition. Integration of metabolome-informed proteome imaging data provides a comprehensive view of underlying biological pathways to inform our understanding of metabolic fungal pathways in plant matter degradation within the micrometer-scale environment.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Determination of latent dimensionality in international trade flow

Currently, high-dimensional data is ubiquitous in data science, which necessitates the development of techniques to decompose and interpret such multidimensional (aka tensor) datasets. Finding a low dimensional representation of the data, that is, its inherent structure, is one of the approaches that can serve to understand the dynamics of low dimensional latent features hidden in the data. Moreover, decomposition methods with non-negative constraints are shown to extract more insightful factors. Nonnegative RESCAL is one such technique, particularly well suited to analyze self-relational data, such as dynamic networks found in international trade flows. Particularly, non-negative RESCAL computes a low dimensional tensor representation by finding the latent space containing multiple modalities. Furthermore, estimating the dimensionality of this latent space is crucial for extracting meaningful latent features. Here, to determine the dimensionality of the latent space with non-negative RESCAL, we propose a latent dimension determination method which is based on clustering of the solutions of multiple realizations of non-negative RESCAL decompositions. We demonstrate the performance of our model selection method on synthetic data. We then apply our method to decompose a network of international trade flows data from International Monetary Fund and shows that with a correct latent dimension determination, the resulting features are able to capture relevant empirical facts from economic literature.

97 MATHEMATICS AND COMPUTING↗

Spectral proper orthogonal decomposition of active wake mixing dynamics in a stable atmospheric boundary layer

Recent advancements in the use of active wake mixing (AWM) to reduce wake effects on downstream turbines open new avenues for increasing power generation in wind farms. However, a better understanding of the fluid dynamics underlying AWM is still needed to make wake mixing a reliable strategy for wind farm flow control. In this work, a spectral proper orthogonal decomposition (SPOD) is used to analyze the dynamics of coherent flow structures that are induced in the wake through blade pitch actuation. The data are generated using the ExaWind software suite to perform large eddy simulations of an NREL 2.8 MW turbine operating in a stable atmospheric boundary layer. SPOD tracks the modal behavior of flow structures from their generation in the turbine induction field through their growth in the near-wake region and to their subsequent evolution and energy transfers in the far wake. SPOD is shown to be a useful tool in the context of AWM because it translates the wavenumber and frequency inputs to the turbine controller to structures in the wake. A decomposition of the radial shear stress flux in the wake is also developed using SPOD to measure the contribution of coherent flow structures to mean flow turbulent entrainment and wake recovery. The effectiveness of AWM is connected to its ability to excite inherent structures in the wake of the turbine that arise using baseline controls. The effects of AWM on blade loading are also analyzed by connecting the axial force along the blade to the SPOD analysis of the turbine induction field. Lastly, the performance of different AWM strategies is demonstrated in a two-turbine array.

17 WIND ENERGY↗

Modal dynamics of wind turbine wake meandering from lidar observations

Horizontal scans from nacelle-mounted lidars provide time series measurements of wind turbine wakes across diverse atmospheric conditions, enabling analysis of coherent turbulent structures that influence wake meandering through proper orthogonal decomposition (POD). While low-order modes capture the most energetic turbulent structures, our analysis reveals that they do not necessarily dominate wake meandering dynamics. We evaluate more than 16,000 combinatorial reconstructions of the flow field for each inflow case, demonstrating that mode relevance depends on mode symmetry, turbulent kinetic energy content, and inflow characteristics. Cases with low turbulence intensity and large integral timescales show stronger correlations between POD modes and wake meandering, whereas higher turbulence conditions (turbulence intensity > 7%) are less effectively described by reduced-order models. However, the qualitative similarity of POD modes across varied atmospheric conditions suggests the potential existence of a semi-universal basis for representing wind turbine wakes, with implications for improving engineering wake models.

17 WIND ENERGY↗

Interface instabilities in hafnium hydride entrained iron metal matrix composites

The chemical interactions in Fe–HfH 2 metal matrix composites (MMCs) are studied across multiple length scales to elucidate the decomposition of the parent phases and corresponding reaction zone physics during direct current sintering. Fe–HfH 2 composites were synthesized with increasing as-mixed hydride contents of Fe–25% HfH 2 , Fe–40% HfH 2 , Fe–55% HfH 2 , and Fe–70% HfH 2 (all in vol. %) to demonstrate the ability to achieve sintered MMCs with target hydride contents. Samples were probed across multiple length scales through a multi-modal workflow employing x-ray diffraction, scanning electron microscopy and segmentation analysis, and synchrotron techniques including hard x-ray fluorescence mapping and nanoprobe x-ray absorption near-edge structure measurements. Under the selected sintering temperature and pressure conditions, hydrogen evolution is seen to evolve through parallel paths: thermal decomposition from during the transformation of HfH 2 to HfH x<2 and through subsequent reaction with the Fe matrix leading to intermetallic phase formation. Specifically, HfFe and HfFe 2 intermetallic formation accelerates the release of hydrogen with a subsequent HfO 2 phase forming at grain boundaries. For this MMC, the consumption or loss of hydrogen can be considerable in compacts with initial hydride loading of 25%–40% HfH 2 approaching 83% hydrogen loss for the lower volume fraction composites. Increasing the volume fraction of HfH 2 to 70% enhanced the retained hydrogen content to 53% and attributed to the reduced interfacial area intrinsic to the increased HfH 2 loading in this MMC.

36 MATERIALS SCIENCE↗

Harnessing on-machine metrology data for prints with a surrogate model for laser powder directed energy deposition

In this study, we leverage the massive amount of multi-modal on-machine metrology data generated from Laser Powder Directed Energy Deposition (LP-DED) to construct a comprehensive surrogate model of the 3D printing process. By employing Dynamic Mode Decomposition with Control (DMDc), a data-driven technique, we capture the complex physics inherent in this extensive dataset. This physics-based surrogate model emphasizes thermodynamically significant quantities, enabling us to accurately predict key process outcomes. The model ingests 21 process parameters, including laser power, scan rate, and position, while providing outputs such as melt pool temperature, melt pool size, and other essential observables. Furthermore, it incorporates uncertainty quantification to provide bounds on these predictions, enhancing reliability and confidence in the results. We then deploy the surrogate model on a new, unseen part and monitor the printing process as validation of the method. Our experimental results demonstrate that the predictions align with actual measurements with high accuracy, confirming the effectiveness of our approach. Furthermore, this methodology not only facilitates real-time predictions but also operates at process-relevant speeds, establishing a basis for implementing feedback control in LP-DED.

Digital twins↗

[Fe 4 S 4 ] cubane in sulfite reductases: new insights into bonding properties and reactivity

The dissimilatory sulfite reductase enzyme has very characteristic active site where the substrate binds to an iron site, ligated by a siroheme macrocycle and a thiol directly connected to a [Fe 4 S 4 ] cluster. This arrangement gives the enzyme remarkable efficiency in reducing sulfite and nitrite all the way to hydrogen sulfide and ammonia. For the first time we present a theoretical study where substrate binding modalities and activation are elucidated using active site models containing proton supply side chains and the [Fe 4 S 4 ] cluster. Density functional theory (DFT) was deployed in conjunction with the energy decomposition scheme (as implemented in AMS), the quantum theory of atoms in molecules (QTAIM), and conceptual DFT (cDFT) descriptors. We quantified the role of the electrostatic interactions inside the active site created by the side chains as well as the influence of the [Fe4S4] cluster on the substrate binding. Furthermore, using conceptual DFT results we shed light of the activation process, thus, laying foundation for further mechanistic studies. We found that the bonding of the ligands to the iron complex is dominated by electrostatic interactions, but the presence of the [Fe 4 S 4 ] cubane leads to substantial changes in electronic interaction. The spin state of the cubane, however, affects the binding energy only marginally. The conceptual DFT results show that the presence of the [Fe 4 S 4 ] cubane affects the reactivity of the active site as it is involved in electron transfer. This is corroborated by an increase in the electrophilicity index, thus making the active site more prone to react with the ligands. The interaction energies between the ligand and the siroheme group are also increased upon the presence of the cubane group, thus, suggesting that the siroheme group is not an innocent spectator but plays an active role in the reactivity of the dSIR active site.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Aeroelastic System Development Using Proper Orthogonal Decomposition and Volterra Theory

This research combines Volterra theory and proper orthogonal decomposition (POD) into a hybrid methodology for reduced-order modeling of aeroelastic systems. The out-come of the method is a set of linear ordinary differential equations (ODEs) describing the modal amplitudes associated with both the structural modes and the POD basis functions for the uid. For this research, the structural modes are sine waves of varying frequency, and the Volterra-POD approach is applied to the fluid dynamics equations. The structural modes are treated as forcing terms which are impulsed as part of the uid model realization. Using this approach, structural and uid operators are coupled into a single aeroelastic operator. This coupling converts a free boundary uid problem into an initial value problem, while preserving the parameter (or parameters) of interest for sensitivity analysis. The approach is applied to an elastic panel in supersonic cross ow. The hybrid Volterra-POD approach provides a low-order uid model in state-space form. The linear uid model is tightly coupled with a nonlinear panel model using an implicit integration scheme. The resulting aeroelastic model provides correct limit-cycle oscillation prediction over a wide range of panel dynamic pressure values. Time integration of the reduced-order aeroelastic model is four orders of magnitude faster than the high-order solution procedure developed for this research using traditional uid and structural solvers.

Lucia, David J.↗

Development and implementation of an in situ high-resolution isotopic microscope for measuring metabolic interactions in soil mesocosms (Final Report)

Many of our planet’s ecosystems rely on the activities of soil microbial communities. These microbes have significant and wide-ranging effects: they metabolize carbon and other nutrients, interact with plants and fungi, and perform other processes important to soil health. However, our ability to directly observe the enzymatic and metabolic activities of microbes within soil is currently limited, not only by the complexity of soil microbial communities themselves, but also by the lack of experimental tools to study them and their molecular interactions in situ. These challenges hinder our understanding of the lifesustaining processes of biomass decomposition and the manner in which it contributes to the movement of freed carbon within soil ecosystems. Here we aim to develop a novel ultrahigh-resolution isotopic microscope that combines complementary imaging modalities to gain insights into metabolic cycling in soil. This final report covers the portion of this work that was completed at the University of North Carolina at Chapel Hill from 2018 to 2019; the research is continuing through 2022 at the University of Massachusetts Chan Medical School.

59 BASIC BIOLOGICAL SCIENCES↗

Design of a linear projector for use with the normal modes of the GLAS 4th order GCM

The design of a linear projector for use with the normal modes of a model of atmospheric circulation is discussed. A central element in any normal mode initialization scheme is the process by which a set of data fields - winds, temperatures or geopotentials, and surface pressures - are expressed ("projected') in terms of the coefficients of a model's normal modes. This process is completely analogous to the Fourier decomposition of a single field (indeed a FFT applied in the zonal direction is a part of the process). Complete separability in all three spatial dimensions is assumed. The basis functions for the modal expansion are given. An important feature of the normal modes is their coupling of the structures of different fields, thus a coefficient in a normal mode expansion would contain both mass and momentum information.

Bloom, S. C.↗

International Space Station Future Correlation Analysis Improvements

Ongoing modal analyses and model correlation are performed on different configurations of the International Space Station (ISS). These analyses utilize on-orbit dynamic measurements collected using four main ISS instrumentation systems: External Wireless Instrumentation System (EWIS), Internal Wireless Instrumentation System (IWIS), Space Acceleration Measurement System (SAMS), and Structural Dynamic Measurement System (SDMS). Remote Sensor Units (RSUs) are network relay stations that acquire flight data from sensors. Measured data is stored in the Remote Sensor Unit (RSU) until it receives a command to download data via RF to the Network Control Unit (NCU). Since each RSU has its own clock, it is necessary to synchronize measurements before analysis. Imprecise synchronization impacts analysis results. A study was performed to evaluate three different synchronization techniques: (i) measurements visually aligned to analytical time-response data using model comparison, (ii) Frequency Domain Decomposition (FDD), and (iii) lag from cross-correlation to align measurements. This paper presents the results of this study.

time synchronization↗

Modal Behavior of Hemispheric Sea Ice Covers

Recent papers have described 18-year trends and annual oscillations in the Arctic and Antarctic sea ice extents, areas, and enclosed open water areas based on a newly-formulated 18.2-year ice concentration time series. This time series includes data for the entire Arctic and Antarctic ice covers, as well as for previously defined subregions consisting of 5 sectors in the Antarctic and 9 regions in the Arctic. It was obtained by fine-tuning the sea ice algorithm tie points individually for each of the four sensors used to acquire the data. In this paper, we extend these analyses to an examination of the intrinsic modes of these time series, obtained by means of Empirical Mode Decomposition, with emphasis on periodicities greater than the annual cycle. Quasibiennial and quasiquadrennial oscillations observed with a different technique and reported earlier for the first 8.8 years of this time series were also observed in the present series. However, the intrinsic modes were not monochromatic; they feature frequency as well as amplitude modulation within their respective frequency bands. Modal periods of up to 18 years are observed, with important implications for the trend analyses published earlier. These results are compared with the oscillations in the Length-of-Day and North Atlantic Oscillation parameters similarly determined for the same 18.2-year period.

Gloersen, Per↗

A decentralized linear quadratic control design method for flexible structures

A decentralized suboptimal linear quadratic control design procedure which combines substructural synthesis, model reduction, decentralized control design, subcontroller synthesis, and controller reduction is proposed for the design of reduced-order controllers for flexible structures. The procedure starts with a definition of the continuum structure to be controlled. An evaluation model of finite dimension is obtained by the finite element method. Then, the finite element model is decomposed into several substructures by using a natural decomposition called substructuring decomposition. Each substructure, at this point, still has too large a dimension and must be reduced to a size that is Riccati-solvable. Model reduction of each substructure can be performed by using any existing model reduction method, e.g., modal truncation, balanced reduction, Krylov model reduction, or mixed-mode method. Then, based on the reduced substructure model, a subcontroller is designed by an LQ optimal control method for each substructure independently. After all subcontrollers are designed, a controller synthesis method called substructural controller synthesis is employed to synthesize all subcontrollers into a global controller. The assembling scheme used is the same as that employed for the structure matrices. Finally, a controller reduction scheme, called the equivalent impulse response energy controller (EIREC) reduction algorithm, is used to reduce the global controller to a reasonable size for implementation. The EIREC reduced controller preserves the impulse response energy of the full-order controller and has the property of matching low-frequency moments and low-frequency power moments. An advantage of the substructural controller synthesis method is that it relieves the computational burden associated with dimensionality. Besides that, the SCS design scheme is also a highly adaptable controller synthesis method for structures with varying configuration, or varying mass and stiffness properties.

Su, Tzu-Jeng↗

Data-driven recovery of hidden physics in reduced order modeling of fluid flows

In this article, we introduce a modular hybrid analysis and modeling (HAM) approach to account for hidden physics in reduced order modeling (ROM) of parameterized systems relevant to fluid dynamics. The hybrid ROM framework is based on using first principles to model the known physics in conjunction with utilizing the data-driven machine learning tools to model the remaining residual that is hidden in data. This framework employs proper orthogonal decomposition as a compression tool to construct orthonormal bases and Galerkin projection (GP) as a model to build the dynamical core of the system. Our proposed methodology hence compensates structural or epistemic uncertainties in models and utilizes the observed data snapshots to compute true modal coefficients spanned by these bases. The GP model is then corrected at every time step with a data-driven rectification using a long short-term memory (LSTM) neural network architecture to incorporate hidden physics. A Grassmannian manifold approach is also adopted for interpolating basis functions to unseen parametric conditions. The control parameter governing the system's behavior is thus implicitly considered through true modal coefficients as input features to the LSTM network. The effectiveness of the HAM approach is then discussed through illustrative examples that are generated synthetically to take hidden physics into account. Furthermore, our approach thus provides insights addressing a fundamental limitation of the physics-based models when the governing equations are incomplete to represent underlying physical processes.

42 ENGINEERING↗

Data-Driven Analysis of Multipactor Dynamics via Dynamic Mode Decomposition

Multipactor effect is a performance-limiting kinetic plasma effect that can occur in high-power microwave and radio frequency (RF) devices. Multipactor effect is of special concern in vacuum or near-vacuum conditions such as those in particle accelerators and spaceborne devices. In this work, we present a data-driven reduced-order model (ROM) based on dynamic mode decomposition (DMD) for modeling of multipactor effects. We study multipactor effects and the resulting nonlinear harmonic generation by processing high-fidelity data generated from electromagnetic particle-in-cell (EMPIC) simulations using the DMD algorithm. We also investigate time-delay embedding extensions of DMD with improved generalizability and accuracy for modeling the electron plasma current density behavior. Here, the results show that DMD provides valuable insights into multipactor phenomena by extracting relevant modal spatiotemporal patterns and frequencies. In addition, DMD offers the potential to time extrapolate EMPIC simulations at a minimal cost, thereby reducing overall simulation time.

43 PARTICLE ACCELERATORS↗

Nonstationary Dynamics Data Analysis with Wavelet-SVD Filtering

Nonstationary time-frequency analysis is used for identification and classification of aeroelastic and aeroservoelastic dynamics. Time-frequency multiscale wavelet processing generates discrete energy density distributions. The distributions are processed using the singular value decomposition (SVD). Discrete density functions derived from the SVD generate moments that detect the principal features in the data. The SVD standard basis vectors are applied and then compared with a transformed-SVD, or TSVD, which reduces the number of features into more compact energy density concentrations. Finally, from the feature extraction, wavelet-based modal parameter estimation is applied.

Brenner, Marty↗

Modal dynamics of high-frequency transverse combustion instabilities

This paper investigates the modal dynamics of self-excited high frequency transverse instabilities in a multi-nozzle, swirl-stabilized can combustor. It utilizes results from multiple pressure sensors to decompose the disturbance into clockwise (CW) and counter-clockwise (CCW) waves, enabling a reconstruction of the pressure amplitude, phase, anti-nodal line location, and spin ratio. While these types of decompositions have been demonstrated on small scale combustors, particularly annular combustors, this is the first study to analyze these results for a high power, multi-nozzle can combustor. Its particular focus is in characterizing the system evolution in phase space and the nature of the system attractors under linearly unstable conditions; results are presented in spin ratio-phase difference space, so that limit cycle oscillations appear as fixed points. Results are shown for conditions where three different attractors dominate (standing, CW, and CCW) but are taken for very long time periods (~50,000 cycles) so that the system's interaction with these attractors can be observed. For the standing mode-dominated case, the phase portrait shows spiral trajectories converging to a single fixed point. For the spinning mode dominated cases, however, multiple fixed points and saddle points were observed, and the trajectories showed a more complicated structure. It was also shown that even though the CW mode is dominant, this mode intermittently switches to standing or CCW modes for short period. These results suggest that the strength of each attractor depends on test conditions and that random noise stochastically perturbs the system so that large regions of the phase space are explored. Lastly, one interesting observation is that the angular velocity of the anti-nodal line linearly varies with the spin ratio under standing-mode dominant conditions.

42 ENGINEERING↗