Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “spectral 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 541 records · Page 30

Photosynthetic accessory pigments - Evidence for the influence of phycoerythrin on the submarine light field

Upwelled water-leaving spectral radiances and the submarine light-field source are theorized to have the physical scattering and absorption effects of photosynthetic accessory pigments. A laser is employed to induce phycoerythrin pigment fluorescence, and the presence of the photosynthetic accessory pigment is measured in oceanic waters with a spectral fluorescence method. Experimental data from active-passive correlation spectroscopy demonstrate that water-leaving spectral radiances of about 600 nm are highly correlated with phycoerythrin pigment fluorescence. Because the data on the phycoerythrin and chlorophyll fluorescences are entered randomly in the correlation the results are considered unambiguous. An estimate of the photopigment is derived that agrees with measurements of its laser-induced fluorescence.

Hoge, Frank E.↗

Uncertainty Propagation for Turbulent, Compressible Flow in a Quasi-1D Nozzle Using Stochastic Methods

This paper describes a fully spectral, Polynomial Chaos method for the propagation of uncertainty in numerical simulations of compressible, turbulent flow, as well as a novel stochastic collocation algorithm for the same application. The stochastic collocation method is key to the efficient use of stochastic methods on problems with complex nonlinearities, such as those associated with the turbulence model equations in compressible flow and for CFD schemes requiring solution of a Riemann problem. Both methods are applied to compressible flow in a quasi-one-dimensional nozzle. The stochastic collocation method is roughly an order of magnitude faster than the fully Galerkin Polynomial Chaos method on the inviscid problem.

Zang, Thomas A.↗

Preserving Superconvergence of Spectral Elements for Curved Domains via h and p-Geometric Refinement [Slides]

Spectral element methods (SEM) are extensions of finite element methods (FEM) that employ Gauss-Lobatto or similar nodes instead of equidistant nodes for high-order elements. SEM can deliver superior accuracy compared to equidistant FEM due to potential superconvergence. However, significant challenges remain for domains with curved boundaries, which have limited the advantages of SEM for real-world applications. In this work, we propose a novel approach to bolster the overall accuracy and preserve the superconvergence of SEM over curved domains.

97 MATHEMATICS AND COMPUTING↗

COMPAS: Compositional mineralogy with a photoacoustic spectrometer

There is an important need for an in situ method of mineral and rock identification and quantification that provides true absorption spectra for a wide spectral range for lunar lander/rover missions. Many common minerals such as feldspars, magnetite, ilmenite, and amorphous fine solids or glasses, can exhibit flat spectral reflectances in the 400-2500 nm spectral region that render inaccurate or difficult their spectral detection and quantitative analysis. Ideal rock and mineral spectra are, of course, pure absorption spectra that are independent of the spectral effects of scattering, particle size, and distribution that can result in a suppression or distortion of their spectral features. This ideal seldom pertains to real samples. Since sample preparation is difficult and may fundamentally alter the observed diffuse spectral reflectance, an in situ spectral measurement method for rocks and minerals on the Moon, insensitive to the sample morphology, would be invaluable. Photoacoustic spectroscopy is a well-established technique appropriate for this task that has been widely applied in condensed-phase spectral studies of complex, highly light scattering, unprepared samples of everything from coal to whole blood, including rock and mineral characterization. A Compositional Mineralogy Photoacoustic Spectrometer, or COMPAS, can enable in situ spectral measurement of rocks and minerals, bypassing the major limitations of diffuse reflectance spectroscopy. COMPAS spectral capabilities for rock and mineral samples will be incorporated into an instrument prototype specifically for lunar measurements, compatible with rover capabilities.

Smith, W. Hayden↗

NekRS, a GPU-accelerated spectral element Navier–Stokes solver

The development of NekRS, a GPU-oriented thermal-fluids simulation code based on the spectral element method (SEM) is described. For performance portability, the code is based on the open concurrent compute abstraction and leverages scalable developments in the SEM code Nek5000 and in libParanumal, which is a library of high-performance kernels for high-order discretizations and PDE-based miniapps. Critical performance sections of the Navier–Stokes time advancement are addressed. Performance results on several platforms are presented here, including scaling to 27,648 V100s on OLCF Summit, for calculations of up to 60B gridpoints.

97 MATHEMATICS AND COMPUTING↗

Two-Dimensional Fourier Transform Analysis of Helicopter Flyover Noise

A method to separate main rotor and tail rotor noise from a helicopter in flight is explored. Being the sum of two periodic signals of disproportionate, or incommensurate frequencies, helicopter noise is neither periodic nor stationary. The single Fourier transform divides signal energy into frequency bins of equal size. Incommensurate frequencies are therefore not adequately represented by any one chosen data block size. A two-dimensional Fourier analysis method is used to separate main rotor and tail rotor noise. The two-dimensional spectral analysis method is first applied to simulated signals. This initial analysis gives an idea of the characteristics of the two-dimensional autocorrelations and spectra. Data from a helicopter flight test is analyzed in two dimensions. The test aircraft are a Boeing MD902 Explorer (no tail rotor) and a Sikorsky S-76 (4-bladed tail rotor). The results show that the main rotor and tail rotor signals can indeed be separated in the two-dimensional Fourier transform spectrum. The separation occurs along the diagonals associated with the frequencies of interest. These diagonals are individual spectra containing only information related to one particular frequency.

SantaMaria, Odilyn L.↗

Two-Dimensional Fourier Transform Analysis of Helicopter Flyover Noise

A method to separate main rotor and tail rotor noise from a helicopter in flight is explored. Being the sum of two periodic signals of disproportionate, or incommensurate frequencies, helicopter noise is neither periodic nor stationary. The single Fourier transform divides signal energy into frequency bins of equal size. Incommensurate frequencies are therefore not adequately represented by any one chosen data block size. A two-dimensional Fourier analysis method is used to separate main rotor and tail rotor noise. The two-dimensional spectral analysis method is first applied to simulated signals. This initial analysis gives an idea of the characteristics of the two-dimensional autocorrelations and spectra. Data from a helicopter flight test is analyzed in two dimensions. The test aircraft are a Boeing MD902 Explorer (no tail rotor) and a Sikorsky S-76 (4-bladed tail rotor). The results show that the main rotor and tail rotor signals can indeed be separated in the two-dimensional Fourier transform spectrum. The separation occurs along the diagonals associated with the frequencies of interest. These diagonals are individual spectra containing only information related to one particular frequency.

SantaMaria, Odilyn L.↗

Two-Dimensional Fourier Transform Analysis of Helicopter Flyover Noise

A method to separate main rotor and tail rotor noise from a helicopter in flight is explored. Being the sum of two periodic signals of disproportionate, or incommensurate frequencies, helicopter noise is neither periodic nor stationary. The single Fourier transform divides signal energy into frequency bins of equal size. Incommensurate frequencies are therefore not adequately represented by any one chosen data block size. A two-dimensional Fourier analysis method is used to separate main rotor and tail rotor noise. The two-dimensional spectral analysis method is first applied to simulated signals. This initial analysis gives an idea of the characteristics of the two-dimensional autocorrelations and spectra. Data from a helicopter flight test is analyzed in two dimensions. The test aircraft are a Boeing MD902 Explorer (no tail rotor) and a Sikorsky S-76 (4-bladed tail rotor). The results show that the main rotor and tail rotor signals can indeed be separated in the two-dimensional Fourier transform spectrum. The separation occurs along the diagonals associated with the frequencies of interest. These diagonals are individual spectra containing only information related to one particular frequency.

SantaMaria, Odilyn L.↗

Recent Applications of Higher-Order Spectral Analysis to Nonlinear Aeroelastic Phenomena

Recent applications of higher-order spectral (HOS) methods to nonlinear aeroelastic phenomena are presented. Applications include the analysis of data from a simulated nonlinear pitch and plunge apparatus and from F-18 flight flutter tests. A MATLAB model of the Texas A&MUniversity s Nonlinear Aeroelastic Testbed Apparatus (NATA) is used to generate aeroelastic transients at various conditions including limit cycle oscillations (LCO). The Gaussian or non-Gaussian nature of the transients is investigated, related to HOS methods, and used to identify levels of increasing nonlinear aeroelastic response. Royal Australian Air Force (RAAF) F/A-18 flight flutter test data is presented and analyzed. The data includes high-quality measurements of forced responses and LCO phenomena. Standard power spectral density (PSD) techniques and HOS methods are applied to the data and presented. The goal of this research is to develop methods that can identify the onset of nonlinear aeroelastic phenomena, such as LCO, during flutter testing.

Silva, Walter A.↗

Turbulence Modeling with Nek5000/RS, SOD2D and Alya

We present Validation and Verification (V & V) study of two high-order spectral element method (SEM) based Computational Fluid Dyanimcs (CFD) codes that we will utilize for work on turbulence modeling: Nek5000 and SOD2D. While the former solves the incompressible form of Navier-Stokes equation, the latter works with compressible set of equations and uses an entropy-viscosity formulation to account for the discontinuities for high Mach number flows. We demonstrate the accuracy of these codes for two benchmark sub-sonic turbulent flows: periodic channel and pipe flow, by carrying out first and second order statistical analysis, grid convergence and turbulent structure analysis using wall-resolved large eddy simulations (WRLES). Later, we report on the implementation and testing of various wall-modeling strategies for large eddy simulation of turbulent flows in Nek5000. These include both classical log-law based and decision-tree based machine-learning models. Accuracy of these closure strategies are analyzed and necessary future work is outlined.

97 MATHEMATICS AND COMPUTING↗

NeuroSEM: A hybrid framework for simulating multiphysics problems by coupling PINNs and spectral elements

Multiphysics problems that are characterized by complex interactions among fluid dynamics, heat transfer, structural mechanics, and electromagnetics, are inherently challenging due to their coupled nature. While experimental data on certain state variables may be available, integrating these data with numerical solvers remains a significant challenge. Physics-informed neural networks (PINNs) have shown promising results in various engineering disciplines, particularly in handling noisy data and solving inverse problems in partial differential equations (PDEs). However, their effectiveness in forecasting nonlinear phenomena in multiphysics regimes, particularly involving turbulence, is yet to be fully established. Here, this study introduces NeuroSEM, a hybrid framework integrating PINNs with the highfidelity Spectral Element Method (SEM) solver, Nektar++. NeuroSEM leverages the strengths of both PINNs and SEM, providing robust solutions for multiphysics problems. PINNs are trained to assimilate data and model physical phenomena in specific subdomains, which are then integrated into the Nektar++ solver. We demonstrate the efficiency and accuracy of NeuroSEM for thermal convection in cavity flow and flow past a cylinder. The framework effectively handles data assimilation by addressing those subdomains and state variables where the data is available. We applied NeuroSEM to the Rayleigh-B´enard convection system, including cases with missing thermal boundary conditions and noisy datasets. Finally, we applied the proposed NeuroSEM framework to real particle image velocimetry (PIV) data to capture flow patterns characterized by horseshoe vortical structures. Our results indicate that NeuroSEM accurately models the physical phenomena and assimilates the data within the specified subdomains. The framework’s plug-and-play nature facilitates its extension to other multiphysics or multiscale problems. Furthermore, NeuroSEM is optimized for efficient execution on emerging integrated GPU-CPU architectures. This hybrid approach enhances the accuracy and efficiency of simulations, making it a powerful tool for tackling complex engineering challenges in various scientific domains.

42 ENGINEERING↗

Development of MOSCATO: A CFD-Level Electrochemistry and Corrosion Simulator for Molten Salt Systems

For both coolant and fueled variants of molten salt reactors (MSRs), the corrosion of structural materials is a significant challenge. The corrosion stems from chemical and electrochemical reactions initiated by fissile material, fission products, and impurities in the salt. Lower-fidelity models rely on empirical correlations for mass transfer, simplified lumped temperature profiles, and similar assumptions. They do not capture detailed spatial variations in complex geometries, creating the need for high-fidelity modeling to bridge this gap.As we approach the demonstration and possible deployment of MSRs in this decade, the development of a high-fidelity, high-performance simulator becomes imperative. To simulate the complex electrochemical environment and corrosion within molten salt systems, we have developed the Molten Salt Chemistry And TranspOrt (MOSCATO) code. This endeavor is comprised of three essential components. First, mass transfer equations are coupled with the Navier-Stokes equations in order to account for the transport of species in the salt. Second, the diffusion of alloy constituents, such as Cr, Fe, Ni, etc. is simulated within the structural metals. Third, the alloy and salt domains are coupled to account for the heterogeneous chemical and electrochemical reactions that occur at the salt-alloy interface.MOSCATO manages all three components within the framework of the highly scalable, open-source spectral element method computational fluid dynamics code Nek5000/NekRS. This integration enables MOSCATO to harness the immense computational power of modern high-performance computing resources, ensuring both high fidelity and computational speed.In addition to code development, we have initiated a comprehensive verification and validation campaign, utilizing data from diverse sources. First, MOSCATO's electrochemical solver was verified with reference numerical data. Then validation occurred against experiments: one of a thermal galvanic cell and the other for corrosion in flowing molten salt of FLiNaK (LiF-NaF-KF). This campaign verified and validated MOSCATO as a reliable tool for simulating electrochemical environments and corrosion in molten salt systems.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Data analysis techniques: Spectral processing

The individual steps in the data processing scheme applied to most radars used for wind sounding are analyzed. This processing method uses spectral analysis and assumes a pulse Doppler radar. Improvement in the signal to noise ratio of some radars is discussed.

Strauch, R. G.↗

Spectral Analysis of Regular Material Point Method and its Application to Study High Pressure Reverse Osmosis Membrane Compaction and Embossing

Material Point Method (MPM) is gaining widespread interest in applied continuum mechanics. The fact that all the continuum properties are stored on the particles (or material points) and the governing equations are solved on these material points makes MPM extremely suited to problems involving severe material deformations, such as crack propagation, soil movement, and fluid flows. Despite its popularity, only a few studies have focused on the numerical properties of MPM. This presentation introduces a global spectral analysis of the regular material point method. Contrary to previous studies, the analysis focuses on the numerical properties of the method in the spectral space. The amplification factor is derived as a function of the non- dimensional wave numbers. It provides insights into the stability and dissipative properties of the method for various CFL and Fourier numbers. The effect of the grid shape functions, number of particles per cell and their locations inside the grid cell are also analyzed. The EXAGOOP MPM solver (https://github.com/NREL/Exagoop.git) is developed at the National Renewable Energy Laboratory as a part of the NAWI UHPRO project and is based on the AMReX framework. A single-level, uniform cartesian grid is used as the background mesh, while the particle class in AMReX is used to manage the material point operations. Linear hat and B-splines are used as grid shape functions, while the time integration is performed using explicit Euler time integration. EXAGOOP is both CPU and GPU compatible and has been demonstrated to work well on multiple compute architectures. The performance of EXAGOOP on various computing architectures is presented along with its application to study compaction and embossing of high-pressure reverse osmosis membranes. The MPM solution accurately reproduces the membrane deformation. The deformed pore size and structure simulated using MPM also agree well with experimental SEM images.

material point method↗

Design of Unstructured Adaptive (UA) NAS Parallel Benchmark Featuring Irregular, Dynamic Memory Accesses

We describe the design of a new method for the measurement of the performance of modern computer systems when solving scientific problems featuring irregular, dynamic memory accesses. The method involves the solution of a stylized heat transfer problem on an unstructured, adaptive grid. A Spectral Element Method (SEM) with an adaptive, nonconforming mesh is selected to discretize the transport equation. The relatively high order of the SEM lowers the fraction of wall clock time spent on inter-processor communication, which eases the load balancing task and allows us to concentrate on the memory accesses. The benchmark is designed to be three-dimensional. Parallelization and load balance issues of a reference implementation will be described in detail in future reports.

Feng, Hui-Yu↗

Archetype-based Redshift Estimation for the Dark Energy Spectroscopic Instrument Survey

We present a computationally efficient galaxy archetype-based redshift estimation and spectral classification method for the Dark Energy Survey Instrument (DESI) survey. The DESI survey currently relies on a redshift fitter and spectral classifier using a linear combination of principal component analysis–derived templates, which is very efficient in processing large volumes of DESI spectra within a short time frame. However, this method occasionally yields unphysical model fits for galaxies and fails to adequately absorb calibration errors that may still be occasionally visible in the reduced spectra. Our proposed approach improves upon this existing method by refitting the spectra with carefully generated physical galaxy archetypes combined with additional terms designed to absorb data reduction defects and provide more physical models to the DESI spectra. We test our method on an extensive data set derived from the survey validation (SV) and Year 1 (Y1) data of DESI. Our findings indicate that the new method delivers marginally better redshift success for SV tiles while reducing catastrophic redshift failure by 10%–30%. At the same time, results from millions of targets from the main survey show that our model has relatively higher redshift success and purity rates (0.5%–0.8% higher) for galaxy targets while having similar success for QSOs. These improvements also demonstrate that the main DESI redshift pipeline is generally robust. Additionally, it reduces the false-positive redshift estimation by 5%–40% for sky fibers. We also discuss the generic nature of our method and how it can be extended to other large spectroscopic surveys, along with possible future improvements.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

International Conference on Numerical Methods in Fluid Dynamics, 11th, Williamsburg, VA, June 27-July 1, 1988, Proceedings

Recent advances in computational fluid dynamics (CFD) are discussed in reviews and reports. Topics addressed include CFD models in plasma dynamics, parallel computation for simulation studies, CFD for hypersonic airbreathing aircraft, multigrid methods for the steady incompressible Navier-Stokes equations, upwind differencing techniques, TV stable schemes for shock-interacting flows, Euler models of hypersonic vortex flows, parallel multilevel adaptive methods, and vortex methods for slightly viscous three-dimensional flows. Consideration is given to the accuracy of node-based solutions on irregular meshes, multigrid calculations for cascades, a finite-volume-element method for planar cavity flow, parallel heterogeneous mesh refinement for advection-diffusion equations, the convergence of the spectral-viscosity method for nonlinear conservation laws, and numerical simulations of Taylor vortices in a spherical gap.

Dwoyer, D. L.↗