SEARCH · Engineering Papers
Results for “spectral expansion”
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.
Inversion of the Solar Rotation Rate Versus Depth and Latitude
Three differential inversion techniques were used to compute the internal solar rotation rate from several sets of n-averaged frequency splittings. An iterative variation of the spectral expansion method; the optimal averaging kernel method; and a piecewise constant constrained least square method were used to invert the data. Each computation was carried out independently. While they present similar trends, each of the solutions differs in detail. A consistent feature in all the inversions is the disappearance of differential rotation below the base of the convection zone. Also, a strong differential signature in the deeper part of the convection zone is present in most of the solutions. A slow decrease of the rotation rate with depth for the equatorial and midlatitude curves is significant in the spectral expansion and the least square results but only marginally apparent in the averaging kernel results.
Effect of Swirl on Turbulent Structures in Supersonic Jets
Direct Numerical Simulation (DNS) is used to study the mechanism of generation and evolution of turbulence structures in a temporally evolving supersonic swirling round jet and also to examine the resulting acoustic radiations. Fourier spectral expansions are used in the streamwise and azimuthal directions and a 1-D b-spline Galerkin representation is used in the radial direction. Spectral-like accuracy is achieved using this numerical scheme. Direct numerical simulations, using the b-spline spectral method, are carried out starting from mean flow initial conditions which are perturbed by the most unstable linear stability eigenfunctions. It is observed that the initial helical instability waves evolve into helical vortices which eventually breakdown into smaller scales of turbulence. 'Rib' structures similar to those seen in incompressible mixing layer flow of Rogers and Moserl are observed. The jet core breakdown stage exhibits increased acoustic radiations.
Effect of Swirl on Turbulent Structures in Supersonic Jets
Direct Numerical Simulation (DNS) is used to study the mechanism of generation and evolution of turbulence structures in a temporally evolving supersonic swirling round jet and also to examine the resulting acoustic radiations. Fourier spectral expansions are used in the streamwise and azimuthal directions and a 1-D b-spline Galerkin representation is used in the radial direction. Spectral-like accuracy is achieved using this numerical scheme. Direct numerical simulations, using the b-spline spectral method, are carried out starting from mean flow initial conditions which are perturbed by the most unstable linear stability eigenfunctions. It is observed that the initial helical instability waves evolve into helical vortices which eventually breakdown into smaller scales of turbulence. 'Rib' structures similar to those seen in incompressible mixing layer flow of Rogers and Moserl are observed. The jet core breakdown stage exhibits increased acoustic radiations.
Effect of Swirl on Turbulent Structures in Supersonic Jets
Direct numerical simulation (DNS) is used to study the mechanism of generation and evolution of turbulence structures in a temporally evolving supersonic swirling round jet and also to examine the resulting acoustic radiations. Fourier spectral expansions are used in the streamwise and azimuthal directions and a 1-D b-spline Galerkin representation is used in the radial direction. Spectral-like accuracy is achieved using this numerical scheme. Direct numerical simulations, using the b-spline spectral method, are carried out starting from mean flow initial conditions which are perturbed by the most unstable linear stability eigenfunctions. It is observed that the initial.helical instability waves evolve into helical vortices which eventually breakdown into smaller scales of turbulence. 'Rib' structures similar to those seen in incompressible mixing layer flow of Rogers and Moser are observed. The jet core breakdown stage exhibits increased acoustic radiations.
Spectral element simulations of laminar and turbulent flows in complex geometries
Spectral element methods are high-order weighted residual techniques based on spectral expansions of variables and geometry for the Navier-Stokes (NS) and transport equations. Here, practical aspects of these methods and their efficient implementation are examined, and several examples of flows in truly complex geometries are presented. The spectral element discretization for NS equations is introduced, and the convergence of the method is addressed. An efficient data management scheme is discussed in the context of parallel processing computations. The method is validated by comparing the spectral element solutions with the exact eigensolutions for the Orr-Sommerfeld equations in two and three dimensions. Computer-aided flow visualizations are presented for an impulsive flow past a sharp edge wedge. Three-dimensional states of channel flow disrupted by an array of cylindrical eddy promoters are studied, and the results of a direct simulation of the turbulent flow in a plane channel are presented.
A new numerical method for the simulation of three dimensional flow in a pipe
A new numerical technique for simulating three dimensional, unsteady, incompressible pipe flows is presented and its utility and accuracy is shown. Each vector function in the expansion of the velocity field is divergence free and satisfies the boundary conditions for viscous flow. Some of the benefits of the expansion technique are that pressure is eliminated from the dynamics, only two unknowns per mesh point are required, implicit treatment of the viscous terms is provided at no extra computational cost, and no fractional time steps are required. The method uses spectral expansions: Fourier series in the azimuthal and streamwise directions, and Jacobi polynominals in the radial direction.
A new numerical method for the simulation of three-dimensional flow in a pipe
A new numerical technique for simulating three dimensional, unsteady, incompressible pipe flows is presented and its utility and accuracy is shown. Each vector function in the expansion of the velocity field is divergence free and satisfies the boundary conditions for viscous flow. Some of the benefits of the expansion technique are that pressure is eliminated from the dynamics, only two unknowns per mesh point are required, implicit treatment of the viscous terms is provided at no extra computational cost, and no fractional time steps are required. The method uses spectral expansions: Fourier series in the azimuthal and streamwise directions, and Jacobi polynominals in the radial direction. Previously announced in STAR as N82-31644
Gravitational waves from binary neutron star mergers with a spectral equation of state
In numerical simulations of binary neutron star systems, the equation of state of the dense neutron star matter is an important factor in determining both the physical realism and the numerical accuracy of the simulations. Some equations of state used in simulations are C 2 or smoother in the pressure/density relationship function, such as a polytropic equation of state, but may not have the flexibility to model stars or remnants of different masses while keeping their radii within known astrophysical constraints. Other equations of state, such as tabular or piece-wise polytropic, may be flexible enough to model additional physics and multiple stars' masses and radii within known constraints, but are not as smooth, resulting in additional numerical error. We will study in this paper a recently developed family of equation of state, using a spectral expansion with sufficient free parameters to allow for a larger flexibility than current polytropic equations of state, and with sufficient smoothness to reduce numerical errors compared to tabulated or piece-wise polytropic equations of state. We perform simulations at three mass ratios with a common chirp mass, using two distinct spectral equations of state, and at multiple numerical resolutions. We evaluate the gravitational waves produced from these simulations, comparing the phase error between resolutions and equations of state, as well as with respect to analytical models. From our simulations we estimate that the phase difference at merger for binaries with a dimensionless weighted tidal deformability difference greater than Δ$\tilde{Λ}$=55 can be captured by the SpEC code for these equations of state.
Anti-symmetric and positivity preserving formulation of a spectral method for Vlasov-Poisson equations
We analyze the anti-symmetric properties of a spectral discretization for the one-dimensional Vlasov-Poisson equations. The discretization is based on a spectral expansion in velocity with the symmetrically weighted Hermite basis functions, central finite differencing in space, and an implicit Runge Kutta integrator in time. The proposed discretization preserves the anti-symmetric structure of the advection operator in the Vlasov equation, resulting in a stable numerical method. We apply such discretization to two formulations: the canonical Vlasov-Poisson equations and their continuously transformed square-root representation. The latter preserves the positivity of the particle distribution function. We derive analytically the conservation properties of both formulations, including particle number, momentum, and energy, which are verified numerically on the following benchmark problems: manufactured solution, linear and nonlinear Landau damping, two-stream instability, bump-on-tail instability, and ion-acoustic wave.
Gaussian process regression constrained by boundary value problems
We develop a framework for Gaussian processes regression constrained by boundary value problems. The framework may be applied to infer the solution of a well-posed boundary value problem with a known second-order differential operator and boundary conditions, but for which only scattered observations of the source term are available. Scattered observations of the solution may also be used in the regression. The framework combines co-kriging with the linear transformation of a Gaussian process together with the use of kernels given by spectral expansions in eigenfunctions of the boundary value problem. Furthermore, it benefits from a reduced-rank property of covariance matrices. We demonstrate that the resulting framework yields more accurate and stable solution inference as compared to physics-informed Gaussian process regression without boundary condition constraints.
A parallel-kinetic-perpendicular-moment model for magnetised plasmas
We describe a new model for the study of weakly collisional, magnetised plasmas derived from exploiting the separation of the dynamics parallel and perpendicular to the magnetic field. This unique system of equations retains the particle dynamics parallel to the magnetic field while approximating the perpendicular dynamics through a spectral expansion in the perpendicular degrees of freedom, analogous to moment-based fluid approaches. In so doing, a hybrid approach is obtained that is computationally efficient enough to allow for larger-scale modelling of plasma systems while eliminating a source of difficulty in deriving fluid equations applicable to magnetised plasmas. We connect this system of equations to historical asymptotic models and discuss advantages and disadvantages of this approach, including the extension of this parallel-kinetic-perpendicular moment beyond the typical region of validity of these more traditional asymptotic models. This paper forms the first of a multi-part series on this new model, covering the theory and derivation, alongside demonstration benchmarks of this approach that include shocks and magnetic reconnection.
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.
(2 + δ)-dimensional theory of the electromechanics of lipid membranes: Electrostatics
The coupling of electric fields to the mechanics of lipid membranes gives rise to intriguing electromechanical behavior, as, for example, evidenced by the deformation of lipid vesicles in external electric fields. Electromechanical effects are relevant for many biological processes, such as the propagation of action potentials in axons and the activation of mechanically gated ion channels. Currently, a theoretical framework describing the electromechanical behavior of arbitrarily curved and deforming lipid membranes does not exist. Purely mechanical models commonly treat lipid membranes as two-dimensional surfaces, ignoring their finite thickness. While holding analytical and numerical merit, this approach cannot describe the coupling of lipid membranes to electric fields and is thus unsuitable for electromechanical models. In a sequence of articles, we derive an effective surface theory of the electromechanics of lipid membranes, called the (2 + δ)-dimensional theory, which has the advantages of surface descriptions while accounting for finite thickness effects. The present article proposes a generic dimension reduction procedure relying on low-order spectral expansions. This procedure is applied to the electrostatics of lipid membranes to obtain the (2 + δ)-dimensional theory that captures potential differences across and electric fields within lipid membranes. Finally, this model is tested on different geometries relevant for lipid membranes, showing good agreement with the corresponding three-dimensional electrostatics theory.
A new approach for resonance treatment of doubly heterogeneous fuel using the RSE method
A new resonance calculation method for the doubly-heterogeneous (DH) fuels such as high- temperature gas-cooled reactor fuel is proposed based on the resonance calculation based on Spectral Expansion (RSE) method. The concept of pointwise disadvantage factor for fuel grain is taken into account to treat the DH fuels. The verification calculation is carried out for simplified single fuel cell and fuel compact consisting of five fuel cells and graphite moderator. The calculation results indicate that the present method can appropriately handle the space- dependent self-shielding effect for DH fuels. (authors)
The effect of nonlinear diffusive heat transport in a simple climate model
The effects of nonlinear diffusive heat transport in North's (1975) model is reexamined. By inserting the nonlinearity into the original energy balance equation rather than the spectrally transformed linear algebraic equations, and including one more mode in the spectral expansion, it is found that nonlinear diffusion does give a significantly different sensitivity for states near the current climate. For linear diffusion, an increase of 0.1% of the solar constant from the present value is sufficient to melt all the ice on the earth, whereas an increase of about five times as large is required for nonlinear diffusion. This difference in sensitivity shows the importance of dynamical modeling in simple climate models.
A spectral numerical method for the Navier-Stokes equations with applications to Taylor-Couette flow
A new spectral method for solving the incompressible Navier-Stokes equations in a plane channel and between concentric cylinders is presented. The method uses spectral expansions which inherently satisfy the boundary conditions and the continuity equation and yield banded matrices which are efficiently solved at each time step. In addition, the number of dependent variables is reduced, resulting in a reduction in computer memory requirements. Several test problems have been computed for the channel flow and for flow between concentric cylinders, including Taylor-Couette flow with axisymmetric Taylor vortices and wavy vortices. In all cases, agreement with available experimental and theoretical results is very good.
Parallel pseudospectral domain decomposition techniques
The influence of interface boundary conditions on the ability to parallelize pseudospectral multidomain algorithms is investigated. Using the properties of spectral expansions, a novel parallel two domain procedure is generalized to an arbitrary number of domains each of which can be solved on a separate processor. This interface boundary condition considerably simplifies influence matrix techniques.