Search NASA⌕ Search

SEARCH · Search NASA

Results for “Mode Solver”

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
All sourcesNASA NTRSDOE OSTI
NASA NTRSDate not supplied

Aeroelastic Trim Drag Optimization of Mach 0.8 Transonic Truss-Braced Wing Aircraft using High-Lift Devices and Control Surfaces

This paper presents an aeroelastic trim drag optimization study of the Mach 0.8 Transonic Truss-Braced Wing (TTBW) aircraft using the High-lift devices and control surfaces. An aero-structural analysis solver VSPAERO with transonic small disturbance, integral boundary-layer, and wing-strut interference corrections coupled to mode shapes computed by NASTRAN using the Galerkin method is developed to provide a rapid aircraft aeroelastic performance evaluation. Three different flight conditions corresponding to Mach 0.8 are selected for the aeroelastic trim drag optimization at the design and off-design cruise lift coefficients. The preliminary optimization results show that the TTBW aircraft with the optimized deflection of the high-lift devices and control surfaces achieves a drag reduction of about 8.2 counts, 6.5 counts, and 9.7 counts corresponding to the lift coefficients 0.661, 0.695, and 0.729, respectively. A high-fidelity CFD solver FUN3D is used to verify the aeroelastic trim drag optimization.

TTBW↗
NASA NTRSDate not supplied

Time-Accurate Computational Fluid Dynamics Analysis of a Supersonic Turbine to Assess Structural Mode Excitation

Gas turbines can experience significant unsteady loading during operation which can cause damage such as blade cracking or bearing wear. While computational fluid dynamics (CFD) modeling of full three-dimensional gas turbine flow to verify meanline design performance is standard practice within the industry, the use of time-accurate CFD makes it possible to assess the unsteady loading which can be especially problematic at operating conditions where unsteady fluid loading aligns with structural modes. The Fluid Dynamics branch (ER42) at NASA Marshall Space Flight Center has developed a time-accurate CFD analysis methodology using the unstructured, density-based, Loci-CHEM solver which is capable of simulating the relative motion required for turbomachinery applications through sliding interface meshes. The highly-parallel unstructured solver allows for full three-dimensional, time-accurate simulation of turbines including upstream and downstream components such as manifolds, nozzles, and stators, any of which could excite a turbine structural mode. Recently, Loci-CHEM was used to analyze the FASTRAC supersonic turbine to assess the loading on the rotor. Overall performance and time-averaged loading on the rotor was computed, including the axial force over the entire rotor disk which is required for force balance calculations on the turbopump. The FASTRAC turbine has supersonic nozzles upstream of the rotor that result in significant unsteady loading on the rotor blades at multiples of the nozzle pass frequency. The unsteady blade pressure loading was decomposed into the frequency domain for direct use in structural forced response analysis. Utilizing this type of CFD analysis coupled with structural analysis has proven to be a valuable tool in identifying structural modes that are being excited by unsteady fluid loading.

M.V. Fischels↗
NASA NTRSDate not supplied

Time-Accurate Computational Fluid Dynamics Analysis of a Supersonic Turbine to Assess Structural Mode Excitation

Gas turbines can experience significant unsteady loading during operation which can cause damage such as blade cracking or bearing wear. While computational fluid dynamics (CFD) modeling of full three-dimensional gas turbine flow to verify meanline design performance is standard practice within the industry, the use of time-accurate CFD makes it possible to assess the unsteady loading which can be especially problematic at operating conditions where unsteady fluid loading aligns with structural modes. The Fluid Dynamics branch (ER42) at NASA Marshall Space Flight Center has developed a time-accurate CFD analysis methodology using the unstructured, density-based, Loci-CHEM solver which is capable of simulating the relative motion required for turbomachinery applications through sliding interface meshes. The highly-parallel unstructured solver allows for full three-dimensional, time-accurate simulation of turbines including upstream and downstream components such as manifolds, nozzles, and stators, any of which could excite a turbine structural mode. Recently, Loci-CHEM was used to analyze the FASTRAC supersonic turbine to assess the loading on the rotor. Overall performance and time-averaged loading on the rotor was computed, including the axial force over the entire rotor disk which is required for force balance calculations on the turbopump. The FASTRAC turbine has supersonic nozzles upstream of the rotor that result in significant unsteady loading on the rotor blades at multiples of the nozzle pass frequency. The unsteady blade pressure loading was decomposed into the frequency domain for direct use in structural forced response analysis. Utilizing this type of CFD analysis coupled with structural analysis has proven to be a valuable tool in identifying structural modes that are being excited by unsteady fluid loading.

M V Fischels↗
NASA NTRS1988

A simplified analysis of the multigrid V-cycle as a fast elliptic solver

For special model problems, Fourier analysis gives exact convergence rates for the two-grid multigrid cycle and, for more general problems, provides estimates of the two-grid convergence rates via local mode analysis. A method is presented for obtaining mutigrid convergence rate estimates for cycles involving more than two grids (using essentially the same analysis as for the two-grid cycle). For the simple cast of the V-cycle used as a fast Laplace solver on the unit square, the k-grid convergence rate bounds obtained by this method are sharper than the bounds predicted by the variational theory. Both theoretical justification and experimental evidence are presented.

Decker, Naomi H.↗
DOE OSTI2022

Numerical algorithms for water waves with background flow over obstacles and topography

Abstract We present two accurate and efficient algorithms for solving the incompressible, irrotational Euler equations with a free surface in two dimensions with background flow over a periodic, multiply connected fluid domain that includes stationary obstacles and variable bottom topography. One approach is formulated in terms of the surface velocity potential while the other evolves the vortex sheet strength. Both methods employ layer potentials in the form of periodized Cauchy integrals to compute the normal velocity of the free surface, are compatible with arbitrary parameterizations of the free surface and boundaries, and allow for circulation around each obstacle, which leads to multiple-valued velocity potentials but single-valued stream functions. We prove that the resulting second-kind Fredholm integral equations are invertible, possibly after a physically motivated finite-rank correction. In an angle-arclength setting, we show how to avoid curve reconstruction errors that are incompatible with spatial periodicity. We use the proposed methods to study gravity-capillary waves generated by flow around several elliptical obstacles above a flat or variable bottom boundary. In each case, the free surface eventually self-intersects in a splash singularity or collides with a boundary. We also show how to evaluate the velocity and pressure with spectral accuracy throughout the fluid, including near the free surface and solid boundaries. To assess the accuracy of the time evolution, we monitor energy conservation and the decay of Fourier modes and compare the numerical results of the two methods to each other. We implement several solvers for the discretized linear systems and compare their performance. The fastest approach employs a graphics processing unit (GPU) to construct the matrices and carry out iterations of the generalized minimal residual method (GMRES).

Ambrose, David M.↗
DOE OSTI2024

Constrained Local Approximate Ideal Restriction for Advection-Diffusion Problems

Herein this paper focuses on developing a reduction-based algebraic multigrid (AMG) method that is suitable for solving general (non)symmetric linear systems and is naturally robust from pure advection to pure diffusion. Initial motivation comes from a new reduction-based AMG approach, $\ell \text{AIR}$ (local approximate ideal restriction), that was developed for solving advection-dominated problems. Though this new solver is very effective in the advection-dominated regime, its performance degrades in cases where diffusion becomes dominant. This is consistent with the fact that in general, reduction-based AMG methods tend to suffer from growth in complexity and/or convergence rates as the problem size is increased, especially for diffusion-dominated problems in two or three dimensions. Motivated by the success of $\ell \text{AIR}$ in the advective regime, our aim in this paper is to generalize the AIR framework with the goal of improving the performance of the solver in diffusion-dominated regimes. To do so, we propose a novel way to combine mode constraints as used commonly in energy-minimization AMG methods with the local approximation of ideal operators used in $\ell \text{AIR}$. The resulting constrained $\ell \text{AIR}$ algorithm is able to achieve fast scalable convergence on advective and diffusive problems. In addition, it is able to achieve standard low complexity hierarchies in the diffusive regime through aggressive coarsening, something that was previously difficult for reduction-based methods.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

DNS of Flame Stabilization Dynamics of a Swirl-Stabilized Spray Burner Using Sustainable Aviation Fuels

Concern for emission reduction has motivated the development of new cost-effective alternative sustainable aviation fuels (SAFs). Drop-in SAFs with blends of the certified and the alternative fuels are beneficial because they do not require engine modifications for use in current aviation engines. Combustion characteristics are of particular concern when comparing a new fuel with Jet-A. Lean blow-off and flame stabilization have been identified to be governed by the fuel propensity to autoignition, i.e. derived cetane number. However, fundamental investigation is required to identify the impact of derived cetane number on the combustion mode of flame stabilization in a realistic combustor. The combustion mode has also a direct impact on turbulent combustion modeling closures. DNS simulations are performed in the low-Mach solver of the Pele Suite called PeleLMeX. Lagrangian multi-phase modeling is used to capture the liquid spray injection of Jet-A (reference fuel) and C1 as a representative of a low cetane number SAF. Adaptive Mesh Refinement (AMR) is used to enable a more efficient simulation of a more realistic domain size and embedded boundary treatment is used to model a bluff-body geometry. Local extinction and edge flame propagation were observed for both fuels. The edge flame propagation mode was quantified in terms of a Damkohler number defined as the ratio between progress variable reaction rate and its diffusive flux. Initial analysis suggests that a mixed mode combustion occurs for the edge flame propagation, with flame propagation assisted by ignition and autoignition co-existing for both Jet-A and C1 flames. The analysis of the extinction region shows a larger progress variable for Jet-A which can further increase the local displacement speed and present a faster reignition of the stoichiometric mixture fraction.

adaptive mesh refinement↗
DOE OSTI2017

Final Technical Report on “Center for Nonlinear Simulation of Energetic Particles in Burning Plasmas (CSEP)”

The gyrokinetic turbulence code GEM is used to study reverse shear Alfven eigenmodes in the DIII-D tokamak, cylindrical tearing mode, the nonlinear saturation mechanism of Alfven eigenmodes, and verify other gyrokinetic codes. It is used to study the stability of Alfven waves in the presence of fusion product alpha particles in ITER. New algorithm features and physics capabilities are added to GEM, including a novel discretization method for the Poisson equation, a vorticity approach to the fluid electron model and a low-n solver in cylindrical geometry.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2016

CLUMPY DISKS AS A TESTBED FOR FEEDBACK-REGULATED GALAXY FORMATION

We study the dependence of fragmentation in massive gas-rich galaxy disks at z >1 on stellar feedback schemes and hydrodynamical solvers, employing the GASOLINE2 SPH code and the lagrangian mesh-less code GIZMO in finite mass mode. Non-cosmological galaxy disk runs with the standard delayed-cooling blastwave feedback are compared with runs adopting a new superbubble feedback, which produces winds by modeling the detailed physics of supernova-driven bubbles and leads to efficient self-regulation of star formation. We find that, with blastwave feedback, massive star-forming clumps form in comparable number and with very similar masses in GASOLINE2 and GIZMO. Typical clump masses are in the range 10{sup 7}–10{sup 8} M {sub ⊙}, lower than in most previous works, while giant clumps with masses above 10{sup 9} M {sub ⊙} are exceedingly rare. By contrast, superbubble feedback does not produce massive star-forming bound clumps as galaxies never undergo a phase of violent disk instability. In this scheme, only sporadic, unbound star-forming overdensities lasting a few tens of Myr can arise, triggered by non-linear perturbations from massive satellite companions. We conclude that there is severe tension between explaining massive star-forming clumps observed at z >1 primarily as the result of disk fragmentation driven by gravitational instability and the prevailing view of feedback-regulated galaxy formation. The link between disk stability and star formation efficiency should thus be regarded as a key testing ground for galaxy formation theory.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS1985

Multigrid solutions to quasi-elliptic schemes

Quasi-elliptic schemes arise from central differencing or finite element discretization of elliptic systems with odd order derivatives on non-staggered grids. They are somewhat unstable and less accurate then corresponding staggered-grid schemes. When usual multigrid solvers are applied to them, the asymptotic algebraic convergence is necessarily slow. Nevertheless, it is shown by mode analyses and numerical experiments that the usual FMG algorithm is very efficient in solving quasi-elliptic equations to the level of truncation errors. Also, a new type of multigrid algorithm is presented, mode analyzed and tested, for which even the asymptotic algebraic convergence is fast. The essence of that algorithm is applicable to other kinds of problems, including highly indefinite ones.

Brandt, A.↗
NASA NTRS1985

Multigrid solutions to quasi-elliptic schemes

Quasi-elliptic schemes arise from central differencing or finite element discretization of elliptic systems with odd order derivatives on non-staggered grids. They are somewhat unstable and less accurate than corresponding staggered-grid schemes. When usual multigrid solvers are applied to them, the asymptotic algebraic convergence is necessarily slow. Nevertheless, it is shown by mode analyses and numerical experiments that the usual FMG algorithm is very efficient in solving quasi-elliptic equations to the level of truncation errors. Also, a new type of multigrid algorithm is presented, mode analyzed and tested, for which even the asymptotic algebraic convergence is fast. The essence of that algorithm is applicable to other kinds of problems, including highly indefinite ones.

Brandt, A.↗
DOE OSTI2025

Dynamic mode decomposition for gyrokinetic eigenmode analysis

Dynamic mode decomposition (DMD) is a post-processing approach to decompose a complex time series into a set of modes via spectral analysis. DMD provides a new and powerful method to recover gyrokinetic drift-wave eigenfrequencies and eigenfunctions based only on the solution of the gyrokinetic-Maxwell initial value problem with almost no added cost to the initial value solver. In the present paper, DMD is applied to the CGYRO gyrokinetic code using a newly-developed CGYRO-DMD post-processor. CGYRO-DMD is numerically efficient, even on a single CPU. It does not set any restrictions on the plasma shape, beta (ratio of the plasma pressure to the magnetic field pressure), collisionality or number of species, and allows one to resolve numerous eigenmodes, even of comparable growth rates. In addition, DMD is not limited to unstable modes, but rather can capture stable and unstable branches simultaneously. In this work, we illustrate the accuracy of DMD through gyrokinetic analysis of mode transition for electromagnetic drift wave instabilities.

drift-wave eigenmodes↗
DOE OSTI2020

A fluid-structure interaction solver for investigating torsional galloping in solar-tracking photovoltaic panel arrays

Solar-tracking photovoltaic arrays are susceptible to aeroelastic fluttering during high-wind events. This dynamic fluttering behavior can grow in amplitude until the panels enter an unstable mode known as torsional galloping which can lead to panel failure or total array destruction. To better understand the physics of the torsional galloping phenomenon and to inform the discussion around panel design and recommended panel stow positions during high wind events, a fluid-structure interaction solver composed of a simulated atmospheric boundary layer with simplified panel structural responses was designed. The simulation choices and features of this solver were informed by the geometry and physical properties of an experimental panel array known to exhibit torsional galloping behavior during hind-wind events. Furthermore, these simulations revealed that the torsional galloping instability is driven by a combination of cyclic vortex shedding from the sun-facing side of the panel and the elastic properties of the torque tube linking the panel assemblies. Testing different stow angles across a range of wind speeds indicates that panels are generally more stable when stowed at negative angles where the leading edge is closer to the ground, hypothesized to be due to ground-blocking effects. These results are supplemented by a discussion of stability trends noted during testing and possible implications when considering multi-row array interactions.

14 SOLAR ENERGY↗
DOE OSTI2020

A Fluid Solver for Studying Torsional Galloping in Solar-Tracking PV Panel Arrays

Solar-tracking photovoltaic arrays are susceptible to aeroelastic fluttering during high-wind events. This dynamic fluttering behavior can grow in amplitude until the panels enter an unstable mode known as torsional galloping which can lead to panel failure or total array destruction. To better understand the physics of the torsional galloping phenomenon, and to inform the discussion around panel design and recommended panel stow positions during high wind events, a fluid-structure interaction solver composed of a simulated atmospheric boundary layer with simplified panel structural responses was designed. The simulation choices and features of this solver were informed by the geometry and physical properties of an experimental panel array known to exhibit torsional galloping behavior during hind-wind events. These simulations revealed that the torsional galloping instability is driven by a combination of cyclic vortex shedding from the sun-facing side of the panel and the elastic properties of the torque tube linking the panel assemblies. Testing different stow angles across a range of wind speeds indicates that panels are generally more stable when stowed at negative angles where the leading edge is closer to the ground, hypothesized to be due to ground-blocking effects.

41 EE - Solar Energy Technologies Office (EE-4S)↗
DOE OSTI2023

The Generalized Green’s function Cluster Expansion: A Python package for simulating polarons

We present an efficient implementation of the Generalized Green’s function Cluster Expansion (GGCE), which is a new method for computing the ground-state properties and dynamics of polarons (single electrons coupled to lattice vibrations) in model electron-phonon systems. The GGCE works at arbitrary temperature and is well suited for a variety of electron-phonon couplings, including, but not limited to, site and bond Holstein and Peierls (Su-SchriefferHeeger) couplings, and couplings to multiple phonon modes with different energy scales and coupling strengths. Quick calculations can be performed efficiently on a laptop using solvers from NumPy and SciPy, or in parallel at scale using the PETSc sparse linear solver engine.

36 MATERIALS SCIENCE↗
NASA NTRS1996

Development of a Driver Code for the WICS Project

Wall Interference Correction System (WICS) is a computational technique to evaluate the wind tunnel wall interference corrections for blockage, Mach number, dynamic pressure and angle of attack. The objective is to predict the effects of the bounding walls on wind tunnel measurements [1]. The procedure attempts to compute these corrections in real time so that the test engineer can implement these by adjusting wind tunnel controls for angle of attack and dynamic pressure while the test is in progress in the 12 ft. pressure wind tunnel at the Ames Research Center. The technique utilizes an ideal flow solver PMARC (Panel Method developed at the Ames Research Center) to compute influence coefficients in an internal flow mode by representing a paneled model of the tunnel for fullspan and semispan tests.

Ghosh, Amitabha↗
DOE OSTI2018

Multiscale two-stage solver for Biot’s poroelasticity equations in subsurface media

In this paper, we propose a two-stage preconditioner for accelerating the iterative solution by a Krylov subspace method of Biot’s poroelasticity equations based on a displacement-pressure formulation. The spatial discretization combines a finite element method for mechanics and a finite volume approach for flow. The fully implicit backward Euler scheme is used for time integration. The result is a 2 × 2 block linear system for each timestep. The preconditioning operator is obtained by applying a two-stage scheme. The first stage is a global preconditioner that employs multiscale basis functions to construct coarse-scale coupled systems using a Galerkin projection. This global stage is effective at damping low-frequency error modes associated with long-range coupling of the unknowns. The second stage is a local block-triangular smoothing preconditioner, which is aimed at high-frequency error modes associated with short-range coupling of the variables. Lastly, various numerical experiments are used to demonstrate the robustness of the proposed solver.

58 GEOSCIENCES↗
DOE OSTI2017

A fast parallel 3D Poisson solver with longitudinal periodic and transverse open boundary conditions for space-charge simulations

A three-dimensional (3D) Poisson solver with longitudinal periodic and transverse open boundary conditions can have important applications in beam physics of particle accelerators. Here, we present a fast efficient method to solve the Poisson equation using a spectral finite-difference method. This method uses a computational domain that contains the charged particle beam only and has a computational complexity of $O(N_u (logN_{mode}))$, where $N_u$ is the total number of unknowns and $N_{mode}$ is the maximum number of longitudinal or azimuthal modes. This saves both the computational time and the memory usage of using an artificial boundary condition in a large extended computational domain. The new 3D Poisson solver is parallelized using a message passing interface (MPI) on multi-processor computers and shows a reasonable parallel performance up to hundreds of processor cores.

97 MATHEMATICS AND COMPUTING↗