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.

SourceResearch connectionFrom yearThrough year

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

At least 55 records · Page 3

All sources NASA NTRS DOE OSTI

DOE OSTI2018

Disconnected Hadronic Vacuum Polarization Contribution to the Muon g-2 with HISQ

We describe a computation of the contribution to the anomalous magnetic moment of the muon from the disconnected part of the hadronic vacuum polarization. We use the highly-improved staggered quark (HISQ) formulation for the current density with gauge configurations generated with four flavors of HISQ sea quarks. The computation is performed by stochastic estimation of the current density using the truncated solver method combined with deflation of low-modes. The parameters are tuned to minimize the computational cost for a given target uncertainty in the current-current correlation function. The calculation presented here is carried out on a single gauge-field ensemble of size $32^3\times48$ with an approximate lattice spacing of $\sim0.15$ fm and with physical sea-quark masses. We describe the methodology and the analysis procedure

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

NASA NTRS2000

ICEG2D: An Integrated Software Package for Automated Prediction of Flow Fields for Single-Element Airfoils with Ice Accretion

An integrated software package, ICEG2D, was developed to automate computational fluid dynamics (CFD) simulations for single-element airfoils with ice accretion. ICEG2D is designed to automatically perform three primary functions: (1) generating a grid-ready, surface definition based on the geometrical characteristics of the iced airfoil surface, (2) generating a high-quality grid using the generated surface point distribution, and (3) generating the input and restart files needed to run the general purpose CFD solver NPARC. ICEG2D can be executed in batch mode using a script file or in an interactive mode by entering directives from a command line. This report summarizes activities completed in the first year of a three-year research and development program to address issues related to CFD simulations for aircraft components with ice accretion. Specifically, this document describes the technology employed in the software, the installation procedure, and a description of the operation of the software package. Validation of the geometry and grid generation modules of ICEG2D is also discussed.

Thompson, David S.↗

DOE OSTI2017

Photonic Crystal-Based High-Power Backward Wave Oscillator

An electron beam traversing a slow wave structure can be used to either generate or amplify electromagnetic radiation through the interaction of the slow space charge wave on the beam with the slow wave structure modes. Here, a cylindrical waveguide with a periodic array of conducting loops is used for the slow wave structure. This paper considers operation as a backward wave oscillator. The dispersion properties of the structure are determined using a frequency-domain eigenmode solver. The interaction of the electron beam with the structure modes is investigated using a 2-D particle-in-cell (PIC) code. In conclusion, the operating frequency and growth rate dependence on beam energy and beam current are investigated using the PIC code and compared with analytic and scaling estimates where possible.

42 ENGINEERING↗

DOE OSTI2026

Novel application of the modal strain energy technique for state-of-the-art damping predictions

Structural damping, which measures the energy dissipation of a vibrating structure, is a key modeling input for lightweight structures but is notoriously hard to predict. This work utilizes vibration-based measurements of centimeter-scale coupons and the modal strain energy approach to predict structural damping of a lightly damped structure. The approach was originally validated with panels shorter than a meter in length. This work extends the validation to a 2.75 m beam made of unidirectional and biaxial glass fiber laminates bonded by adhesive. The comparison between three-dimensional finite element model predictions and full-scale experimental measurements of damping show an average error of 5.2% for the first five modes. Additionally, the modal strain energy approach is newly applied with a one-dimensional geometrically exact beam theory model and a two-dimensional sectional analysis solver. This beam approach accurately predicts the damping behavior of the first bending modes but loses accuracy for higher order modes that are dominated by three-dimensional effects. This novel approach provides faster simulations while allowing arbitrary beam cross sections. The paper also investigates traditional and high-force dynamic mechanical analysis to measure structural damping of coupons. Both alternatives show significant errors in attempted validation against the theoretical thermoelastic damping of aluminum coupons.

17 WIND ENERGY↗

NASA NTRS2004

Analysis Tools for CFD Multigrid Solvers

Analysis tools are needed to guide the development and evaluate the performance of multigrid solvers for the fluid flow equations. Classical analysis tools, such as local mode analysis, often fail to accurately predict performance. Two-grid analysis tools, herein referred to as Idealized Coarse Grid and Idealized Relaxation iterations, have been developed and evaluated within a pilot multigrid solver. These new tools are applicable to general systems of equations and/or discretizations and point to problem areas within an existing multigrid solver. Idealized Relaxation and Idealized Coarse Grid are applied in developing textbook-efficient multigrid solvers for incompressible stagnation flow problems.

Mineck, Raymond E.↗

NASA NTRS2004

Analysis Tools for CFD Multigrid Solvers

Mineck, Raymond E.↗

DOE OSTI2017

Controlling the numerical Cerenkov instability in PIC simulations using a customized finite difference Maxwell solver and a local FFT based current correction

In this study we present a customized finite-difference-time-domain (FDTD) Maxwell solver for the particle-in-cell (PIC) algorithm. The solver is customized to effectively eliminate the numerical Cerenkov instability (NCI) which arises when a plasma (neutral or non-neutral) relativistically drifts on a grid when using the PIC algorithm. We control the EM dispersion curve in the direction of the plasma drift of a FDTD Maxwell solver by using a customized higher order finite difference operator for the spatial derivative along the direction of the drift (1ˆ direction). We show that this eliminates the main NCI modes with moderate |k 1 |, while keeps additional main NCI modes well outside the range of physical interest with higher |k 1 |. These main NCI modes can be easily filtered out along with first spatial aliasing NCI modes which are also at the edge of the fundamental Brillouin zone. The customized solver has the possible advantage of improved parallel scalability because it can be easily partitioned along 1ˆ which typically has many more cells than other directions for the problems of interest. We show that FFTs can be performed locally to current on each partition to filter out the main and first spatial aliasing NCI modes, and to correct the current so that it satisfies the continuity equation for the customized spatial derivative. This ensures that Gauss’ Law is satisfied. Lastly, we present simulation examples of one relativistically drifting plasma, of two colliding relativistically drifting plasmas, and of nonlinear laser wakefield acceleration (LWFA) in a Lorentz boosted frame that show no evidence of the NCI can be observed when using this customized Maxwell solver together with its NCI elimination scheme.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2021

Drift reduced Landau fluid model for magnetized plasma turbulence simulations in BOUT++ framework

Recently the drift-reduced Landau fluid six-field turbulence model within the BOUT++ framework has been upgraded. In particular, this new model employs a new normalization, adds a volumetric flux-driven source option, the Landau fluid closure for parallel heat flux and a Laplacian inversion solver which is able to capture $n$ = 0 axisymmetric mode evolution in realistic tokamak configurations. As we report here, these improvements substantially extended model's capability to study a wider range of tokamak edge phenomena, and are essential to build a fully self-consistent edge turbulence model capable of both transient (e.g., ELM, disruption) and transport time-scale simulations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

NASA NTRS2011

Progress Toward Efficient Laminar Flow Analysis and Design

A multi-fidelity system of computer codes for the analysis and design of vehicles having extensive areas of laminar flow is under development at the NASA Langley Research Center. The overall approach consists of the loose coupling of a flow solver, a transition prediction method and a design module using shell scripts, along with interface modules to prepare the input for each method. This approach allows the user to select the flow solver and transition prediction module, as well as run mode for each code, based on the fidelity most compatible with the problem and available resources. The design module can be any method that designs to a specified target pressure distribution. In addition to the interface modules, two new components have been developed: 1) an efficient, empirical transition prediction module (MATTC) that provides n-factor growth distributions without requiring boundary layer information; and 2) an automated target pressure generation code (ATPG) that develops a target pressure distribution that meets a variety of flow and geometry constraints. The ATPG code also includes empirical estimates of several drag components to allow the optimization of the target pressure distribution. The current system has been developed for the design of subsonic and transonic airfoils and wings, but may be extendable to other speed ranges and components. Several analysis and design examples are included to demonstrate the current capabilities of the system.

Campbell, Richard L.↗

NASA NTRS2004

MPT Prediction of Aircraft-Engine Fan Noise

A collection of computer programs has been developed that implements a procedure for predicting multiple-pure-tone (MPT) noise generated by fan blades of an aircraft engine (e.g., a turbofan engine). MPT noise arises when the fan is operating with supersonic relative tip Mach No. Under this flow condition, there is a strong upstream running shock. The strength and position of this shock are very sensitive to blade geometry variations. For a fan where all the blades are identical, the primary tone observed upstream of the fan will be the blade passing frequency. If there are small variations in geometry between blades, then tones below the blade passing frequency arise MPTs. Stagger angle differences as small as 0.1 can give rise to significant MPT. It is also noted that MPT noise is more pronounced when the fan is operating in an unstarted mode. Computational results using a three-dimensional flow solver to compute the complete annulus flow with non-uniform fans indicate that MPT noise can be estimated in a relatively simple way. Hence, once the effect of a typical geometry variation of one blade in an otherwise uniform blade row is known, the effect of all the blades being different can be quickly computed via superposition. Two computer programs that were developed as part of this work are used in conjunction with a user s computational fluid dynamics (CFD) code to predict MPT spectra for a fan with a specified set of geometric variations: (1) The first program ROTBLD reads the users CFD solution files for a single blade passage via an API (Application Program Interface). There are options to replicate and perturb the geometry with typical variations stagger, camber, thickness, and pitch. The multi-passage CFD solution files are then written in the user s file format using the API. (2) The second program SUPERPOSE requires two input files: the first is the circumferential upstream pressure distribution extracted from the CFD solution on the multi-passage mesh, the second file defines the geometry variations of each blade in a complete fan. Superposition is used to predict the spectra resulting from the geometric variations.

Connell, Stuart D.↗

NASA NTRS2002

A Method for Optimizing Non-Axisymmetric Liners for Multimodal Sound Sources

Central processor unit times and memory requirements for a commonly used solver are compared to that of a state-of-the-art, parallel, sparse solver. The sparse solver is then used in conjunction with three constrained optimization methodologies to assess the relative merits of non-axisymmetric versus axisymmetric liner concepts for improving liner acoustic suppression. This assessment is performed with a multimodal noise source (with equal mode amplitudes and phases) in a finite-length rectangular duct without flow. The sparse solver is found to reduce memory requirements by a factor of five and central processing time by a factor of eleven when compared with the commonly used solver. Results show that the optimum impedance of the uniform liner is dominated by the least attenuated mode, whose attenuation is maximized by the Cremer optimum impedance. An optimized, four-segmented liner with impedance segments in a checkerboard arrangement is found to be inferior to an optimized spanwise segmented liner. This optimized spanwise segmented liner is shown to attenuate substantially more sound than the optimized uniform liner and tends to be more effective at the higher frequencies. The most important result of this study is the discovery that when optimized, a spanwise segmented liner with two segments gives attenuations equal to or substantially greater than an optimized axially segmented liner with the same number of segments.

Watson, W. R.↗

NASA NTRS2003

Computational Evaluation of the Steady and Pulsed Jet Effects on the Performance of a Circulation Control Wing Section

Circulation Control Wing (CCW) technology is a very effective way of achieving very high lift coefficients needed by aircraft during take-off and landing. This technology can also be used to directly control the flow field over the wing. Compared to a conventional high-lift system, a Circulation Control Wing (CCW) can generate the required values of lift coefficient C(sub L,max) during take-off/landing with fewer or no moving parts and much less complexity. Earlier designs of CCW configurations used airfoils with a large radius rounded trailing edge to maximize the lift benefit. However, these designs also produced very high drag. These high drag levels associated with the blunt, large radius trailing edge can be prohibitive under cruise conditions when Circulation Control is no longer necessary. To overcome this difficulty, an advanced CCW section, i.e., a circulation hinged flap was developed to replace the original rounded trailing edge CC airfoil. This concept developed by Englar is shown. The upper surface of the CCW flap is a large-radius arc surface, but the lower surface of the flap is flat. The flap could be deflected from 0 degrees to 90 degrees. When an aircraft takes-off or lands, the flap is deflected as in a conventional high lift system. Then this large radius on the upper surface produces a large jet turning angle, leading to high lift. When the aircraft is in cruise, the flap is retracted and a conventional sharp trailing edge shape results, greatly reducing the drag. This kind of flap does have some moving elements that increase the weight and complexity over an earlier CCW design. But overall, the hinged flap design still maintains most of the Circulation Control high lift advantages, while greatly reducing the drag in cruising condition associated with the rounded trailing edge CCW design. In the present work, an unsteady three-dimensional Navier-Stokes analysis procedure has been developed and applied to this advanced CCW configuration. The solver can be used in both a 2-D and a 3-D mode, and can thus model airfoils as well as finite wings. The jet slot location, slot height, and the flap angle can all be varied easily and individually in the grid generator and the flow solver. Steady jets, pulsed jets, the leading edge and trailing edge blowing can all be studied with this solver.

Liu, Yi↗

DOE OSTI2022

Hybrid gyrokinetic ion/fluid electron simulation of toroidal tearing modes

The effects of toroidicity and kinetic ions on the resistive tearing mode are systematically studied with the gyrokinetic particle-in-cell simulation code GEM and compared with analytic theory. A new field solver in toroidal geometry has been developed for the simulation of low-n (n = 1, 2) modes in tokamaks. It is found that the toroidal effect significantly reduces the growth rate of the tearing mode. The toroidal effect can also increase the radial width of the tearing mode and change the scaling between the radial mode width and resistivity due to the toroidal pressure term in the electron continuity equation. The kinetic effects of ions can decrease the growth rate of the tearing mode. The plasma flux-surface shaping is found to have significant effect on the tearing mode.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

NASA NTRSDate not supplied

Aeroelastic Analysis of Highly Flexible Wings with Linearized Frequency-Domain Aerodynamics

Flutter analysis of configurations with geometric structural nonlinearities typically is done with time-domain analysis. The results from this process are computationally expensive and can yield cumbersome results that may be difficult to manage/interpret. Compared to time-domain methods, frequency-domain flutter analysis can provide additional insight into the characteristics of a flutter stability problem. By linearizing the aeroelastic problem about the nonlinear equilibrium state, this work applies frequency-domain aeroelastic analysis to the Pazy wing, the subject of the Large Deformation Working Group in the Aeroelastic Prediction Workshop. Generalized aerodynamic forces (GAFs) are computed with both a doublet-lattice method and a computational fluid dynamics solver at a range of reduced frequencies as well as a range of dynamic pressures to account for the dependence of the mode shapes on the nonlinear equilibrium state. These GAFs are used in a p-k flutter solver, which is modified to handle the dependence of the stiffness matrix and GAFs on the dynamic pressure.

Bret K Stanford↗

DOE OSTI2023

Data-driven time-scale separation of ODE right-hand sides using dynamic mode decomposition and time delay embedding

Multi-physics simulation often involves multiple different scales. The ARKODE ODE solver package in the SUNDIALS library addresses multi-scale problems with a multi-rate time-integrator that can work with a right-hand side that has fast scale and slow scale components. In this report, we use dynamic mode decomposition and time delay embedding to extract the fast and slow components of the right-hand sides of a simple ODE from data. We then use the extracted components to solve the ODE with ARKODE. Finally, to move towards a real-world use case, we attempt to extract fast and slow scale dynamics from synthetic seismic modeling data.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2023

Accurate simulation of direct laser acceleration in a laser wakefield accelerator

In a laser wakefield accelerator (LWFA), an intense laser pulse excites a plasma wave that traps and accelerates electrons to relativistic energies. When the pulse overlaps the accelerated electrons, it can enhance the energy gain through direct laser acceleration (DLA) by resonantly driving the betatron oscillations of the electrons in the plasma wave. The traditional particle-in-cell (PIC) algorithm, although often the tool of choice to study DLA, contains inherent errors due to numerical dispersion and the time staggering of the electric and magnetic fields. Furthermore, conventional PIC implementations cannot reliably disentangle the fields of the plasma wave and laser pulse, which obscures interpretation of the dominant acceleration mechanism. Here, a customized field solver that reduces errors from both numerical dispersion and time staggering is used in conjunction with a field decomposition into azimuthal modes to perform PIC simulations of DLA in an LWFA. Comparisons with traditional PIC methods, model equations, and experimental data show improved accuracy with the customized solver and convergence with an order-of-magnitude fewer cells. Furthermore, the azimuthal-mode decomposition reveals that the most energetic electrons receive comparable energy from DLA and LWFA.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2026

Modeling Temperature Profiles in the Pedestal of NSTX with Reduced Models

This paper describes new modeling capabilities for predicting H-mode pedestal profiles in spherical tokamaks. Temperature profiles for NSTX discharges 132543 and 132588 are modeled by coupling the \textsc{astra} transport solver with neoclassical transport and gyrokinetic-based reduced models for electron temperature gradient (ETG) and kinetic ballooning mode (KBM) instabilities. A quasi-linear surrogate model for ion-scale transport is developed using linear \textsc{gene} simulations, requiring only a single free parameter calibrated to one discharge. Time-evolving the temperatures with fixed density yields good agreement with experiments for both discharges. Systematic analysis of the transport mechanisms reveals that neoclassical transport is huge across the entire pedestal region for the ion channel. ETG turbulence is large in the plasma edge and low density gradient region, contributing substantially to the electron channel. However, KBM/MHD-like modes also drive significant transport in both the ion and electron thermal channels, making them essential for accurate pedestal modeling. Further refinements, including explicit E×B shear suppression and scaled ETG transport, produce quantitative but not qualitative improvements. This work lays the foundation for predictive modeling of future devices. This paper is on arxiv and has been submitted to Nuclear Fusion.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

NASA NTRS2009

Development of Unsteady Aerodynamic and Aeroelastic Reduced-Order Models Using the FUN3D Code

Recent significant improvements to the development of CFD-based unsteady aerodynamic reduced-order models (ROMs) are implemented into the FUN3D unstructured flow solver. These improvements include the simultaneous excitation of the structural modes of the CFD-based unsteady aerodynamic system via a single CFD solution, minimization of the error between the full CFD and the ROM unsteady aero- dynamic solution, and computation of a root locus plot of the aeroelastic ROM. Results are presented for a viscous version of the two-dimensional Benchmark Active Controls Technology (BACT) model and an inviscid version of the AGARD 445.6 aeroelastic wing using the FUN3D code.

Silva, Walter A.↗