SEARCH · Engineering Papers

Results for “spectral element method”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 217 records · Page 12

All sources NASA NTRS DOE OSTI

DOE OSTI2016

Toward Linear Stability Analysis of Upper Plenum Flows

Linear stability analysis is a tool used to help identify the sensitivity of a particular domain to perturbations and evaluate the system dynamics as a transition, if any, occurs. The development of a time invariant base flow for future use in a linear stability analysis is the focus of this investigation. A test case, with a Reynolds number of 100, was brought to a steady state using the spectral element code Nek5000 and the method of selective frequency damping. The convergence criteria for this method is the volume averaged value of the forcing function which, after approximately 80,000 time-steps, decreased to 3.6x10{sup -4}. Analysis shows that the case has yet to have reached the tolerance desired, 1x10{sup -6}. With the simulation in the process of converging a preliminary perturbation analysis was conducted. Due to the low Reynolds number of this case, 100, it was hypothesized that the flow would be linearly stable. This preliminary stability analysis only focused on the largest fundamental mode and the energy growth of that perturbation over time. Upon evaluation of the perturbation energy, which was shown to be monotonically decreasing in time, it was concluded that, despite not yet being converged, the system is in fact linearly stable. Future work includes evaluating the domain at higher Reynolds numbers, solving the eigenvalue problem created when the linearized Navier-Stokes are written in operator form and calculating the energy evolution of the optimal perturbation. (authors)

42 ENGINEERING↗

NASA NTRS2016

Tensor-Product Preconditioners for Higher-Order Space-Time Discontinuous Galerkin Methods

space-time discontinuous-Galerkin spectral-element discretization is presented for direct numerical simulation of the compressible Navier-Stokes equat ions. An efficient solution technique based on a matrix-free Newton-Krylov method is developed in order to overcome the stiffness associated with high solution order. The use of tensor-product basis functions is key to maintaining efficiency at high order. Efficient preconditioning methods are presented which can take advantage of the tensor-product formulation. A diagonalized Alternating-Direction-Implicit (ADI) scheme is extended to the space-time discontinuous Galerkin discretization. A new preconditioner for the compressible Euler/Navier-Stokes equations based on the fast-diagonalization method is also presented. Numerical results demonstrate the effectiveness of these preconditioners for the direct numerical simulation of subsonic turbulent flows.

Diosady, Laslo T.↗

NASA NTRS2013

Design of a Variational Multiscale Method for Turbulent Compressible Flows

A spectral-element framework is presented for the simulation of subsonic compressible high-Reynolds-number flows. The focus of the work is maximizing the efficiency of the computational schemes to enable unsteady simulations with a large number of spatial and temporal degrees of freedom. A collocation scheme is combined with optimized computational kernels to provide a residual evaluation with computational cost independent of order of accuracy up to 16th order. The optimized residual routines are used to develop a low-memory implicit scheme based on a matrix-free Newton-Krylov method. A preconditioner based on the finite-difference diagonalized ADI scheme is developed which maintains the low memory of the matrix-free implicit solver, while providing improved convergence properties. Emphasis on low memory usage throughout the solver development is leveraged to implement a coupled space-time DG solver which may offer further efficiency gains through adaptivity in both space and time.

Design↗

DOE OSTI2024

A robust spectral element implementation of the $k - τ$ RANS model in Nek5000/NekRS

The $k - ω$ Reynolds Averaged Navier Stokes (RANS) model is one of the industry standard approaches for modeling of turbulent flows. It performs better than the $k - ϵ$ model for low Reynolds number flows and is also more suitable for boundary layers with adverse pressure gradients. Major drawback of the model, however, is that the asymptotic value of $ω$ at the walls is singular, necessitating the use of a contrived “sufficiently” large value for $ω$ as the boundary condition for its transport equation. Here, this invariably leads to the solution being sensitive to near wall grid spacing. While an acceptable solution for low order (finite volume) methods, the excessive near wall gradients lead to persistent numerical stability issues in high order codes. To alleviate the problem, specifically in the context of the high order spectral element code Nek5000, a regularized $k - ω$ approach was formulated in our prior work (Tomboulides et al., 2018). The formulation, however, relies on the use of wall distance and its gradients for modeling the closure terms and can pose problems for simulations in complex geometries. This work presents a novel implementation of the $k - τ$ RANS model in Nek5000, where $τ = 1/ω$, eliminating the need for regularization, owing to the asymptotically bounded behavior of the source terms in the $τ$ transport equation, and also eliminating dependence on wall distance. Robustness and stability of the $k - τ$ model is ensured through implicit treatment of the source terms and their careful numerical implementation and demonstrated through several cases aimed at verification and validation. Studies include both canonical and engineering relevant problems, viz., turbulent channel flow, pipe flow, backward facing step, flow over NACA 0012 airfoil and flow in a T-junction. Results from the $k - τ$ model are shown to be consistent with regularized $k - ω$ model and also with the $k - ω$ SST model in OpenFOAM (for select studies). Comparison with experimental data is also shown, where available, to bolster validation efforts for the $k - τ$ model implementation through prediction of key turbulent quantities of interest.

Nek5000↗

NASA NTRS1983

Spectral methods for the Euler equations

Spectral methods for compressible flows are introduced in relation to finite difference and finite element techniques within the framework of the method of weighted residuals. Current spectral collocation methods are put in historical context. The basic concepts of both Fourier and Chebyshev spectral collocation methods are provided. Filtering strategies for both shock-fitting and shock-capturing approaches are also presented. Fourier shock capturing techniques are evaluated using a one-dimensional, periodic astrophysical 'nozzle' problem. Examples of shock-fitting approaches include a shock/acoustic wave interaction, shock/vortex interaction, and the classical blunt body problem. While the shock capturing spectral method does not yet show a clear advantage over second-order finite differences, equivalent accuracy can be obtained using shock fitting with far fewer grid points.

Hussaini, M. Y.↗

NASA NTRS1984

Spectral methods for the Euler equations: Fourier methods and shock-capturing

Hussaini, M. Y.↗

NASA NTRS1985

Spectral methods for the Euler equations. I - Fourier methods and shock capturing

Hussaini, M. Y.↗

NASA NTRS1988

The analysis of spectra of novae taken near maximum

A project to analyze ultraviolet spectra of novae obtained at or near maximum optical light is presented. These spectra are characterized by a relatively cool continuum with superimposed permitted emission lines from ions such as Fe II, Mg II, and Si II. Spectra obtained late in the outburst show only emission lines from highly ionized species and in many cases these are forbidden lines. The ultraviolet data will be used with calculations of spherical, expanding, stellar atmospheres for novae to determine elemental abundances by spectral line synthesis. This method is extremely sensitive to the abundances and completely independent of the nebular analyses usually used to obtain novae abundances.

Stryker, L. L.↗

NASA NTRS1984

Spectral methods for exterior elliptic problems

Spectral approximations for exterior elliptic problems in two dimensions are discussed. As in the conventional finite difference or finite element methods, the accuracy of the numerical solutions is limited by the order of the numerical farfield conditions. A spectral boundary treatment is introduced at infinity which is compatible with the infinite order interior spectral scheme. Computational results are presented to demonstrate the spectral accuracy attainable. Although a simple Laplace problem is examined, the analysis covers more complex and general cases.

Canuto, C.↗

DOE OSTI2021

Separating Physics and Dynamics Grids for Improved Computational Efficiency in Spectral Element Earth System Models

Previous studies have shown that atmospheric models with a spectral element grid can benefit from putting physics calculations on a relatively coarse finite volume grid. Here we demonstrate an alternative high-order, element-based mapping approach used to implement a quasi-equal-area, finite volume physics grid in E3SM. Unlike similar methods, the new method in E3SM requires topology data purely local to each spectral element, which trivially allows for regional mesh refinement. Simulations with physics grids defined by 2 × 2, 3 × 3, and 4 × 4 divisions of each element are shown to verify that the alternative physics grid does not qualitatively alter the model solution. The model performance is substantially affected by the reduction of physics columns when using the 2 × 2 grid, which can increase the throughput of physics calculations by roughly 60%–120% depending on whether the computational resources are configured to maximize throughput or efficiency. A pair of regionally refined cases are also shown to highlight the refinement capability.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2023

Implementation of a High-Fidelity Interface Resolving Method in Nek5000

The development and utilization of computational fluid dynamics (CFD) models for large, high-temperature electric melters in the Waste Treatment and Immobilization Plant (WTP) in eastern Washington State have proven to be valuable for various purposes. These models allow for a better understanding of the physio-chemical processes occurring within the melter vessels and can contribute to improving operational efficiency, throughput, and addressing operational issues related to vitrification. The CFD models employed for these melter vessels incorporate multiphase fluid flow and heat transfer simulations in different regions, including the plenum, cold cap, and molten glass regions. As the tank waste and glass formers are introduced into the melter, a reacting batch layer known as the cold cap forms on top of the molten glass. To enhance the melt rate, forced convection bubblers located at the bottom of the melters generate convection currents that help homogenize the molten glass and provide heat to the cold cap. As the bubbles rise through the highly viscous glass, they adopt a spherical-cap shape [1]. Meanwhile, the conversion of the batch to glass generates significant amounts of gases (such as water vapor, carbon dioxide, sulfur dioxide, and NOx) due to thermal decomposition [2]. These gases become trapped between the cold cap and molten glass, forming a foam layer [3]. For modeling multiphase flow in CFD and heat transfer simulations of waste glass melters at different scales, efforts are underway to augment the capabilities of the Nek5000 [4] and NekRS [5] open-source codes [6]. Nek5000/NekRS is a scalable and efficient spectral element code that has been successfully applied to a wide range of fluid dynamics problems. By leveraging the Nek5000/NekRS software, it becomes possible to model the melter systems more affordably and with lower computational requirements compared to currently utilized commercial CFD software. The specific objectives of this ongoing effort include: 1. Implementation of a level set method in Nek5000/NekRS: The level set method is a numerical technique commonly employed in CFD simulations to track and represent the interface between different phases or materials accurately. By incorporating this method into Nek5000/NekRS, the ability to simulate multiphase flows in waste glass melters at a high level of fidelity can be achieved. 2. Demonstration of capability for air bubbling through molten glass: As part of the development process, a specific case of air bubbling through molten glass will be simulated using the augmented Nek5000/NekRS code. This demonstration aims to showcase the ability of the software to accurately capture and analyze the complex phenomena involved in the multiphase flow within waste glass melters. By achieving these objectives, the improved Nek5000/NekRS code will offer a powerful computational tool for simulating and analyzing waste glass melter systems, enabling better understanding, optimization, and troubleshooting of these vitrification processes. The ability to accurately model and simulate multiphase flows has broad relevance across many industries and scientific domains, and the improved functionality can contribute to advancements in various fields beyond waste glass melter simulations.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

DOE OSTI2016

NO THERMAL INVERSION AND A SOLAR WATER ABUNDANCE FOR THE HOT JUPITER HD 209458B FROM HST /WFC3 SPECTROSCOPY

The nature of the thermal structure of hot Jupiter atmospheres is one of the key questions raised by the characterization of transiting exoplanets over the past decade. There have been claims that many hot Jupiters exhibit atmospheric thermal inversions. However, these claims have been based on broadband photometry rather than the unambiguous identification of emission features with spectroscopy, and the chemical species that could cause the thermal inversions by absorbing stellar irradiation at high altitudes have not been identified despite extensive theoretical and observational effort. Here we present high-precision Hubble Space Telescope WFC3 observations of the dayside thermal emission spectrum of the hot Jupiter HD 209458b, which was the first exoplanet suggested to have a thermal inversion. In contrast to previous results for this planet, our observations detect water in absorption at 6.2 σ confidence. When combined with Spitzer photometry, the data are indicative of a monotonically decreasing temperature with pressure over the range of 1–0.001 bars at 7.7 σ confidence. We test the robustness of our results by exploring a variety of model assumptions, including the temperature profile parameterization, presence of a cloud, and choice of Spitzer data reduction. We also introduce a new analysis method to determine the elemental abundances from the spectrally retrieved mixing ratios with thermochemical self-consistency and find plausible abundances consistent with solar metallicity (0.06–10 × solar) and carbon-to-oxygen ratios less than unity. This work suggests that high-precision spectrophotometric results are required to robustly infer thermal structures and compositions of extrasolar planet atmospheres and to perform comparative exoplanetology.

79 ASTRONOMY AND ASTROPHYSICS↗

NASA NTRS1988

Computational structural mechanics and fluid dynamics: Advances and trends; Proceedings of the Symposium, Washington, DC, Oct. 17-19, 1988

Recent advances in computational structural and fluid dynamics are discussed in reviews and reports. Topics addressed include fluid-structure interaction and aeroelasticity, CFD techniques for reacting flows, micromechanics, stability and eigenproblems, probabilistic methods and chaotic dynamics, and perturbation and spectral methods. Consideration is given to finite-element, finite-volume, and boundary-element methods; adaptive methods; parallel processing machines and applications; and visualization, mesh generation, and AI interfaces.

Noor, Ahmed K.↗

NASA NTRS1997

High precision computing with charge domain devices and a pseudo-spectral method therefor

The present invention enhances the bit resolution of a CCD/CID MVM processor by storing each bit of each matrix element as a separate CCD charge packet. The bits of each input vector are separately multiplied by each bit of each matrix element in massive parallelism and the resulting products are combined appropriately to synthesize the correct product. In another aspect of the invention, such arrays are employed in a pseudo-spectral method of the invention, in which partial differential equations are solved by expressing each derivative analytically as matrices, and the state function is updated at each computation cycle by multiplying it by the matrices. The matrices are treated as synaptic arrays of a neural network and the state function vector elements are treated as neurons. In a further aspect of the invention, moving target detection is performed by driving the soliton equation with a vector of detector outputs. The neural architecture consists of two synaptic arrays corresponding to the two differential terms of the soliton-equation and an adder connected to the output thereof and to the output of the detector array to drive the soliton equation.

Barhen, Jacob↗

NASA NTRS1992

A spectral overlay method for dissimilar materials

The spectral overlay method is applied to examine the stress distribution along the interface between the two different materials. The essential feature of the spectral overlay method is that the high resolution of localized steep gradients can be achieved by overlaying a spectral interpolant on a standard finite element mesh. In order to evaluate this method and compare it with conventional finite element method, a test problem with high gradients in Poisson equation has been devised. The results show that this method is very powerful in capturing the shape of the gradient field and its peak value. Solutions for two problems with interfaces between dissimilar materials are then given to illustrate the effectiveness of this method.

Belytschko, T.↗

DOE OSTI2017

Pillars of the Mantle: Imaging the Interior of the Earth with Adjoint Tomography. In: PEARC17 Proceedings of the Practice and Experience in Advanced Research Computing 2017 on Sustainability, Success and Impact, Article No. 75

In this work, we investigate global seismic tomographic models obtained by spectral-element simulations of seismic wave propagation and adjoint methods. Global crustal and mantle models are obtained based on an iterative conjugate-gradient type of optimization scheme. Forward and adjoint seismic wave propagation simulations, which result in synthetic seismic data to make measurements and data sensitivity kernels to compute gradient for model updates, respectively, are performed by the SPECFEM3D_GLOBE package [1] [2] at the Oak Ridge Leadership Computing Facility (OLCF) to study the structure of the Earth at unprecedented levels. Using advances in solver techniques that run on the GPUs on Titan at the OLCF, scientists are able to perform large-scale seismic inverse modeling and imaging. Using seismic data from global and regional networks from global CMT earthquakes, scientists are using SPECFEM3D_GLOBE to understand the structure of the mantle layer of the Earth. Visualization of the generated data sets provide an effective way to understand the computed wave perturbations which define the structure of mantle in the Earth.

Pugmire, David↗

NASA NTRS1988

Semi-discrete approximations to nonlinear systems of conservation laws; consistency and L(infinity)-stability imply convergence

A convergence theory for semi-discrete approximations to nonlinear systems of conservation laws is developed. It is shown, by a series of scalar counter-examples, that consistency with the conservation law alone does not guarantee convergence. Instead, a notion of consistency which takes into account both the conservation law and its augmenting entropy condition is introduced. In this context it is concluded that consistency and L(infinity)-stability guarantee for a relevant class of admissible entropy functions, that their entropy production rate belongs to a compact subset of H(loc)sup -1 (x,t). One can now use compensated compactness arguments in order to turn this conclusion into a convergence proof. The current state of the art for these arguments includes the scalar and a wide class of 2 x 2 systems of conservation laws. The general framework of the vanishing viscosity method is studied as an effective way to meet the consistency and L(infinity)-stability requirements. How this method is utilized to enforce consistency and stability for scalar conservation laws is shown. In this context we prove, under the appropriate assumptions, the convergence of finite difference approximations (e.g., the high resolution TVD and UNO methods), finite element approximations (e.g., the Streamline-Diffusion methods) and spectral and pseudospectral approximations (e.g., the Spectral Viscosity methods).

Tadmor, Eitan↗

NASA NTRS2004

Parallel 3D Mortar Element Method for Adaptive Nonconforming Meshes

High order methods are frequently used in computational simulation for their high accuracy. An efficient way to avoid unnecessary computation in smooth regions of the solution is to use adaptive meshes which employ fine grids only in areas where they are needed. Nonconforming spectral elements allow the grid to be flexibly adjusted to satisfy the computational accuracy requirements. The method is suitable for computational simulations of unsteady problems with very disparate length scales or unsteady moving features, such as heat transfer, fluid dynamics or flame combustion. In this work, we select the Mark Element Method (MEM) to handle the non-conforming interfaces between elements. A new technique is introduced to efficiently implement MEM in 3-D nonconforming meshes. By introducing an "intermediate mortar", the proposed method decomposes the projection between 3-D elements and mortars into two steps. In each step, projection matrices derived in 2-D are used. The two-step method avoids explicitly forming/deriving large projection matrices for 3-D meshes, and also helps to simplify the implementation. This new technique can be used for both h- and p-type adaptation. This method is applied to an unsteady 3-D moving heat source problem. With our new MEM implementation, mesh adaptation is able to efficiently refine the grid near the heat source and coarsen the grid once the heat source passes. The savings in computational work resulting from the dynamic mesh adaptation is demonstrated by the reduction of the the number of elements used and CPU time spent. MEM and mesh adaptation, respectively, bring irregularity and dynamics to the computer memory access pattern. Hence, they provide a good way to gauge the performance of computer systems when running scientific applications whose memory access patterns are irregular and unpredictable. We select a 3-D moving heat source problem as the Unstructured Adaptive (UA) grid benchmark, a new component of the NAS Parallel Benchmarks (NPB). In this paper, we present some interesting performance results of ow OpenMP parallel implementation on different architectures such as the SGI Origin2000, SGI Altix, and Cray MTA-2.

Feng, Huiyu↗