Search NASA⌕ Search

SEARCH · Search NASA

Results for “pseudospectral 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.

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

At least 37 records · Page 2

Information content in spectral calculations

Analytical procedures for extracting piecewise smooth solutions of hyperbolic systems from raw oscillatory data obtained by pseudospectral methods are developed. The validity of the approach is demonstrated for the case of linear problems with constant coefficients, and plausibility arguments are presented which indicate its applicability to nonlinear operators when the steady state has been achieved. Numerical results for the development of an oblique shock when a wedge is inserted at zero angle of attack into a uniform supersonic flow of an ideal gas (the time-dependent two-dimensional Euler equations discretized in space by the pseudospectral Chebyshev method) are presented in tables and graphs.

Abarbanel, S.↗

Effect of spatial resolution on apparent sensitivity to initial conditions of a decaying flow as it becomes turbulent

In order to check for spurious chaos and obtain superior solutions for decaying Navier-Stokes flows, an investigation is conducted of the effect of spatial resolution on numerical results. The fourth-order finite difference method results obtained with grids of 32-cubed and 64-cubed points, and those of a pseudospectral method for 128-cubed points, indicate that the sensitivity of initially neighboring solutions to small changes in initial conditions increases with improving spatial resolution.

Deissler, Robert G.↗

Effect of spatial resolution on apparent sensitivity to initial conditions of a decaying flow as it becomes turbulent

Grids with 32(exp 3), 64(exp 3), and 128(exp 3) points are used in numerical solutions for a decaying flow. The sensitivity of initially neighboring solutions to small changes in initial conditions increases as the spatial resolution improves. A fourth-order finite-difference method is used for the solutions with 32(exp 3) and 64(exp 3) grid points, and a pseudospectral method is used for 128(exp 3) grid points. The latter solutions appear to be rather well-resolved, in spite of the formation of steep velocity gradients in the flow.

Deissler, Robert G.↗

Error analysis for spectral approximation of the Korteweg-De Vries equation

The conservation and convergence properties of spectral Fourier methods for the numerical approximation of the Korteweg-de Vries equation are analyzed. It is proved that the (aliased) collocation pseudospectral method enjoys the same convergence properties as the spectral Galerkin method, which is less effective from the computational point of view. This result provides a precise mathematical answer to a question raised by several authors in recent years.

Maday, Y.↗

Analysis of the low gravity tolerance of Bridgman-Stockbarger crystal growth. I - Steady and impulse accelerations

The effects of steady and impulse-type residual accelerations on dopant distributions during directional solidification in 2D and 3D 'generic' models of the Bridgman-Stockbarger technique are investigated using numerical methods. The calculations are based on the thermophysical properties of molten germanium doped with a low concentration of gallium. A Chebyshev collocation pseudospectral method is used for the solution of the governing momentum-, mass-, species-, and heat-transfer equations. Only convection caused by temperature gradients is considered. It is found that lateral nonuniformity in composition is very sensitive to the orientation of the steady component of the residual gravity vector and to the particular operating conditions under consideration. It is also found that laterally or radially averaged composition profiles are alone insufficient to describe the extent of residual convection in a spacecraft environment. The effects of impulse-type disturbances can be severe and can extend for times on the order of 1000 sec after the termination of the impulse.

Alexander, J. Iwan D.↗

Recent progress in numerical techniques for flow simulation

Recent developments in the use of numerical methods for fluid flow simulation show an increasing tendency to use numerical operators that can, in various ways, be factored. Use of methods having this property often increases the accuracy and efficiency of computer codes. Tracing the factorization property provides a unity to the basic concepts involved in the development of cyclic reduction, predictor-corrector, splitting, fast Fourier transform, and pseudospectral methods.-

Lomax, H.↗

Numerical simulation of thermocapillary flow under zero and low gravity conditions

This paper discusses the numerical solution methods and results of steady and unsteady thermocapillary (surface-tension) and buoyancy driven flows in 2D cavities and liquid columns. The 2D cavity was assumed to be square with one free surface with a zero Capillary number (i.e., the free surface was constrained to be flat). A pseudospectral method was used to solve steady and unsteady surface tension-driven and mixed buoyancy-surface tension flows in a square cavity. For the liquid column a finite-difference scheme based on a Picard iteration was used to solve for the flow, temperature and free surface shape. The surface of the liquid column was allowed to deform and, as for the 2D cavity, the surface tension was assumed to depend on temperature.

Alexander, J. I. D.↗

On the boundary treatment in spectral methods for hyperbolic systems

Spectral methods were successfully applied to the simulation of slow transients in gas transportation networks. Implicit time advancing techniques are naturally suggested by the nature of the problem. The correct treatment of the boundary conditions are clarified in order to avoid any stability restriction originated by the boundaries. The Beam and Warming and the Lerat schemes are unconditionally linearly stable when used with a Chebyshev pseudospectral method. Engineering accuracy for a gas transportation problem is achieved at Courant numbers up to 100.

Canuto, C.↗

On the boundary treatment in spectral methods for hyperbolic systems

Spectral methods were successfully applied to the simulation of slow transients in gas transportation networks. Implicit time advancing techniques are naturally suggested by the nature of the problem. The correct treatment of the boundary conditions is clarified in order to avoid any stability restriction originated by the boundaries. The Beam and Warming and the Lerat schemes are unconditionally linearly stable when used with a Chebyshev pseudospectral method. Engineering accuracy for a gas transportation problem is achieved at Courant numbers up to 100.

Canuto, Claudio↗

Accurate numerical simulations of open quantum systems using spectral tensor trains

Decoherence between qubits is a major bottleneck in quantum computations. Decoherence results from intrinsic quantum and thermal fluctuations as well as noise in the external fields that perform the measurement and preparation processes. With prescribed colored noise spectra for intrinsic and extrinsic noise, we present a numerical method, Quantum Accelerated Stochastic Propagator Evaluation (Q-ASPEN), to solve the time-dependent noise-averaged reduced density matrix in the presence of intrinsic and extrinsic noise. Q-ASPEN is arbitrarily accurate and can be applied to provide estimates for the resources needed to error-correct quantum computations. We employ spectral tensor trains, which combine the advantages of tensor networks and pseudospectral methods, as a variational ansatz to the quantum relaxation problem and optimize the ansatz using methods typically used to train neural networks. Here, the spectral tensor trains in Q-ASPEN make accurate calculations with tens of quantum levels feasible. We present benchmarks for Q-ASPEN on the spin-boson model in the presence of intrinsic noise and on a quantum chain of up to 32 sites in the presence of extrinsic noise. In our benchmark, the memory cost of Q-ASPEN scales as a low-order polynomial in the size of the system once the number of system states surpasses the number of basis functions used in the spectral expansion.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Large eddy simulation of incompressible turbulent channel flow

The three-dimensional, time-dependent primitive equations of motion were numerically integrated for the case of turbulent channel flow. A partially implicit numerical method was developed. An important feature of this scheme is that the equation of continuity is solved directly. The residual field motions were simulated through an eddy viscosity model, while the large-scale field was obtained directly from the solution of the governing equations. An important portion of the initial velocity field was obtained from the solution of the linearized Navier-Stokes equations. The pseudospectral method was used for numerical differentiation in the horizontal directions, and second-order finite-difference schemes were used in the direction normal to the walls. The large eddy simulation technique is capable of reproducing some of the important features of wall-bounded turbulent flows. The resolvable portions of the root-mean square wall pressure fluctuations, pressure velocity-gradient correlations, and velocity pressure-gradient correlations are documented.

Moin, P.↗

Classifying photonic topology using the spectral localizer and numerical K -theory

Recently, the spectral localizer framework has emerged as an efficient approach for classifying topology in photonic systems featuring local nonlinearities and radiative environments. In nonlinear systems, this framework provides rigorous definitions for concepts such as topological solitons and topological dynamics, where a system’s occupation induces a local change in its topology due to nonlinearity. For systems embedded in radiative environments that do not possess a shared bulk spectral gap, this framework enables the identification of local topology and shows that local topological protection is preserved despite the lack of a common gap. However, as the spectral localizer framework is rooted in the mathematics of C*-algebras, and not vector bundles, understanding and using this framework requires developing intuition for a somewhat different set of underlying concepts than those that appear in traditional approaches for classifying material topology. In this tutorial, we introduce the spectral localizer framework from a ground-up perspective and provide physically motivated arguments for understanding its local topological markers and associated local measure of topological protection. In doing so, we provide numerous examples of the framework’s application to a variety of topological classes, including crystalline and higher-order topology. We then show how Maxwell’s equations can be reformulated to be compatible with the spectral localizer framework, including the possibility of radiative boundary conditions. To aid in this introduction, we also provide a physics-oriented introduction to multi-operator pseudospectral methods and numerical K-theory, two mathematical concepts that form the foundation for the spectral localizer framework. Finally, we provide some mathematically oriented comments on the C*-algebraic origins of this framework, including a discussion of real C*-algebras and graded C*-algebras that are necessary for incorporating physical symmetries. Looking forward, we hope that this tutorial will serve as an approachable starting point for learning the foundations of the spectral localizer framework.

97 MATHEMATICS AND COMPUTING↗

Resolution requirements for numerical modeling of neutrino quantum kinetics

Neutrino quantum kinetics is a rapidly evolving field in computational astrophysics, with a primary focus on collective neutrino oscillations in core-collapse supernovae and postmerger phases of binary neutron star mergers. In recent years, there has been considerable debate concerning resolution dependence in numerical simulations. In this paper, we conduct a comprehensive resolution study in both angular- and spatial directions by using two independent schemes of quantum kinetic neutrino transport: finite volume and pseudospectral methods. We complement our discussion by linear stability analysis including inhomogeneous modes. Our result suggests that decreasing spatial resolutions underestimates the growth of flavor instability, and then leads to wrong asymptotic states of flavor conversions, which potentially has a critical impact on astrophysical consequences. We further delve into numerical results of low resolution simulations, that reveals the underlying mechanism responsible for numerical artifacts caused by insufficient resolutions. Finally, this study settles the debate on requirements of resolutions and serves as a guideline for numerical modeling of quantum kinetic neutrino transport.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Shock-fitted Euler solutions to shock vortex interactions

The interaction of a planar shock wave with one or more vortexes is computed using a pseudospectral method and a finite difference method. The development of the spectral method is emphasized. In both methods the shock wave is fitted as a boundary of the computational domain. The results show good agreement between both computational methods. The spectral method is, however, restricted to smaller time steps and requires use of filtering techniques.

Salas, M. D.↗

Shock-fitted Euler solutions to shock-vortex interactions

The interaction of a shock wave with a hot spot, a single vortex and a vortex street is studied within the framework of the two dimensional compressible Euler equations. The numerical results obtained by the pseudospectral method and the finite difference MacCormack method are compared. In both the methods the shock wave is fitted as a boundary of the computational domain.

Salas, M. E.↗

Pseudospectral solution of two-dimensional gas-dynamic problems

Chebyshev pseudospectral methods are used to compute two dimensional smooth compressible flows. Grid refinement tests show that spectral accuracy can be obtained. Filtering is not needed if resolution is sufficiently high and if boundary conditions are carefully prescribed.

Kopriva, D. A.↗

Modeling of unsteady small disturbance transonic flow using parametric differentiation, pseudospectral analysis and finite-differencing

A procedure for solving the nonlinear unsteady small disturbance transonic equation is formulated. This procedure is a synthesis of a pseudospectral method, parametric differentiation, and finite differencing. It is equally applicable to lifting and nonlifting airfoils. The procedure is particularly well suited to aeroelastic stability studies that may require a set of affine solutions.

Brown, M.↗

Shock-fitted Euler solutions to shock-vortex interactions

The interaction of a planar shock wave with one or more vortexes is computed using a pseudospectral method and a finite difference method. The development of the spectral method is emphasized. In both methods the shock wave is fitted as a boundary of the computational domain. The results show good agreement between both computational methods. The spectral method is, however, restricted to smaller time steps and requires use of filtering techniques. Previously announced in STAR as N82-28061

Salas, M. D.↗