Search NASA⌕ Search

SEARCH · Search NASA

Results for “numerical 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 631 records · Page 35

GRMHD Simulations of Jet Formation with a New Code

We have developed a new three-dimensional general relativistic magnetohydrodynamic (GRMHD) code by using a conservative, high-resolution shock-capturing scheme. The numerical fluxes are calculated using the HLL approximate Riemann solver scheme. The flux-interpolated, constrained transport scheme is used to maintain a divergence-free magnetic field. Various one-dimensional test problems in both special and general relativity show significant improvements over our previous model. We have performed simulations of jet formations from a geometrically thin accretion disk near both nonrotating and rotating black holes. The new simulation results show that the jet is formed in the same manner as in previous work and propagates outward. In the rotating black hole cases, jets form much closer to the black hole's ergosphere and the magnetic field is strongly twisted due the frame-dragging effect. As the magnetic field strength becomes weaker, a larger amount of matter is launched with the jet. On the other hand, when the magnetic field strength becomes stronger, the jet has less matter and becomes poynting-flux dominated. We will also discuss how the jet properties depend on the rotation of a black hole.

Mizuno, Y.↗

GRMHD Simulations of Jet Formation with a Newly-Developed GRMHD Code

We have developed a new three-dimensional general relativistic magnetohydrodynamic code by using a conservative, high-resolution shock-capturing scheme. The numerical fluxes are calculated using the HLL approximate Riemann solver scheme. The flux-CT scheme is used to maintain a divergence-free magnetic field. Various 1-dimensional test problems show significant improvements over our previous GRMHD code. We have performed simulations of jet formations from a geometrically thin accretion disk near a non-rotating and a rotating black hole. The new simulation results show that the jet is formed by the same manner as in previous works and propagates outward. As the magnetic field strength becomes weaker, larger amount of matter launches with the jet. On the other hand when the magnetic field strength becomes stronger, the jet has less-matter and becomes poynting flux dominated. We will also discuss how the jet properties depend on the rotation of a black hole.

Mizuno, Y.↗

Hybrid-PIC Computer Simulation of the Plasma and Erosion Processes in Hall Thrusters

HPHall software simulates and tracks the time-dependent evolution of the plasma and erosion processes in the discharge chamber and near-field plume of Hall thrusters. HPHall is an axisymmetric solver that employs a hybrid fluid/particle-in-cell (Hybrid-PIC) numerical approach. HPHall, originally developed by MIT in 1998, was upgraded to HPHall-2 by the Polytechnic University of Madrid in 2006. The Jet Propulsion Laboratory has continued the development of HPHall-2 through upgrades to the physical models employed in the code, and the addition of entirely new ones. Primary among these are the inclusion of a three-region electron mobility model that more accurately depicts the cross-field electron transport, and the development of an erosion sub-model that allows for the tracking of the erosion of the discharge chamber wall. The code is being developed to provide NASA science missions with a predictive tool of Hall thruster performance and lifetime that can be used to validate Hall thrusters for missions.

Hofer, Richard R.↗

Flow Simulation of Supersonic Inlet with Bypass Annular Duct

A relaxed isentropic compression supersonic inlet is a new concept that produces smaller cowl drag than a conventional inlet, but incurs lower total pressure recovery and increased flow distortion in the (radially) outer flowpath. A supersonic inlet comprising a bypass annulus to the relaxed isentropic compression inlet dumps out airflow of low quality through the bypass duct. A reliable computational fluid dynamics solution can provide considerable useful information to ascertain quantitatively relative merits of the concept, and further provide a basis for optimizing the design. For a fast and reliable performance evaluation of the inlet performance, an equivalent axisymmetric model whose area changes accounts for geometric and physical (blockage) effects resulting from the original complex three-dimensional configuration is proposed. In addition, full three-dimensional calculations are conducted for studying flow phenomena and verifying the validity of the equivalent model. The inlet-engine coupling is carried out by embedding numerical propulsion system simulation engine data into the flow solver for interactive boundary conditions at the engine fan face and exhaust plane. It was found that the blockage resulting from complex three-dimensional geometries in the bypass duct causes significant degradation of inlet performance by pushing the terminal normal shock upstream.

Kim, HyoungJin↗

Magnetic Shielding of the Channel Walls in a Hall Plasma Accelerator

In a qualification life test of a Hall thruster it was found that the erosion of the acceleration channel practically stopped after approx 5,600 h. Numerical simulations using a two-dimensional axisymmetric plasma solver with a magnetic field-aligned mesh reveal that when the channel receded from its early-in-life to its steady-state configuration the following changes occurred near the wall: (1) reduction of the electric field parallel to the wall that prohibited ions from acquiring significant impact kinetic energy before entering the sheath, (2) reduction of the potential fall in the sheath that further diminished the total energy ions gained before striking the material, and (3) reduction of the ion number density that decreased the flux of ions to the wall. All these changes, found to have been induced by the magnetic field, constituted collectively an effective shielding of the walls from any significant ion bombardment. Thus, we term this process in Hall thrusters "magnetic shielding."

ion density↗

Fluid Transient Analysis During Priming of Evacuated Line

Pressure surges are critical in the design of spacecraft propellant feed lines. The pressure transients that occur during priming of feed lines are very important in the design and analysis of liquid propulsion systems. During the start-up of the propulsion system of a spacecraft, the process of filling of an evacuated pipeline is called priming. Priming can generate severe pressure peaks due to the slam (water hammer) of the propellant against a closed thruster valve. The downstream conditions strongly affect the pressure surge. In space systems, satellites, or interplanetary probes, the propellant lines are vacuum-pumped or filled with low pressure helium or nitrogen before the launch. Before operations, these lines are primed with a vaporizing liquid, sometimes in the presence of a non-condensable gas (NCG), which produces water hammer phenomena. The objective of the current study is to use a finite volume based network flow solver (Generalized Fluid System Simulation Program, GFSSP) for the numerical simulation of Priming in (a) a straight feedline and (b) a flow network. The geometrical configurations and dimensions for the pipe and other components used for the current study are identical to experimental study of Prickett et al.

Bandyopadhyay, Alak↗

Low level coupling scheme between neutronics and thermal-hydraulics based on Anderson acceleration

The simulation of nuclear reactors is a multiphysics problem mixing, amongst other fields, neutron transport and thermal-hydraulics. The simplest and most used approach in multiphysics simulation is based on the coupling of single-physics codes in a black-box fashion. However, in order to reduce the computational time needed for such simulations, case-dependent optimizations are often required. In this paper, we aim at reducing the computational time required to solve a coupled neutronic/thermal-hydraulic steady-state problem on a simplified Pressurized Water Reactor (PWR) core. The idea is to deal simultaneously with the coupling of the energy groups of the deterministic neutronic description of the core and its thermal-hydraulic description with the Anderson acceleration. By doing so, the fission source terms are directly accelerated instead of the power map as done in most cases. The power method used to solve the k-eigenvalue problem inside the neutronic solver is thus accelerated with the Anderson acceleration. The numerical experimentations conducted in this work are performed using APOLLO3 and THEDI, and indicate that such coupling strategy improves the convergence rates in terms of number of iterations required and the total computational time. (authors)

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Euler equation computations for the flow over a hovering helicopter rotor

A numerical solution technique is developed for computing the flow field around an isolated helicopter rotor in hover. The flow is governed by the compressible Euler equations which are integrated using a finite volume approach. The Euler equations are coupled to a free wake model of the rotary wing vortical wake. This wake model is incorporated into the finite volume solver using a prescribed flow, or perturbation, technique which eliminates the numerical diffusion of vorticity due to the artificial viscosity of the scheme. The work is divided into three major parts: (1) comparisons of Euler solutions to experimental data for the flow around isolated wings show good agreement with the surface pressures, but poor agreement with the vortical wake structure; (2) the perturbation method is developed and used to compute the interaction of a streamwise vortex with a semispan wing. The rapid diffusion of the vortex when only the basic Euler solver is used is illustrated, and excellent agreement with experimental section lift coefficients is demonstrated when using the perturbation approach; and (3) the free wake solution technique is described and the coupling of the wake to the Euler solver for an isolated rotor is presented. Comparisons with experimental blade load data for several cases show good agreement, with discrepancies largely attributable to the neglect of viscous effects. The computed wake geometries agree less well with experiment, the primary difference being that too rapid a wake contraction is predicted for all the cases.

Roberts, Thomas Wesley↗

Realizability-preserving discontinuous Galerkin method for spectral two-moment radiation transport in special relativity

Here we present a realizability-preserving numerical method for solving a spectral two-moment model to simulate the transport of massless, neutral particles interacting with a steady background material moving with relativistic velocities. The model is obtained as the special relativistic limit of a four-momentum-conservative general relativistic two-moment model. Using a maximum-entropy closure, we solve for the Eulerian-frame energy and momentum. The proposed numerical method is designed to preserve moment realizability, which corresponds to moments defined by a nonnegative phase-space density. The realizability-preserving method is achieved with the following key components: (i) a discontinuous Galerkin phase-space discretization with specially constructed numerical fluxes in the spatial and energy dimensions; (ii) a strong stability-preserving implicit-explicit time-integration method; (iii) a realizability-preserving conserved to primitive moment solver; (iv) a realizability-preserving implicit collision solver; and (v) a realizability-enforcing limiter. Component (iii) is necessitated by the closure procedure, which closes higher order moments nonlinearly in terms of primitive moments. The nonlinear conserved to primitive and the implicit collision solves are formulated as fixed-point problems, which are solved with custom iterative solvers designed to preserve the realizability of each iterate. With a series of numerical tests, we demonstrate the accuracy and robustness of this discontinuous-Galerkin-implicit-explicit method.

79 ASTRONOMY AND ASTROPHYSICS↗

Improve Problem Solving Skills through Adapting Programming Tools

There are numerous ways for engineers and students to become better problem-solvers. The use of command line and visual programming tools can help to model a problem and formulate a solution through visualization. The analysis of problem attributes and constraints provide insight into the scope and complexity of the problem. The visualization aspect of the problem-solving approach tends to make students and engineers more systematic in their thought process and help them catch errors before proceeding too far in the wrong direction. The problem-solver identifies and defines important terms, variables, rules, and procedures required for solving a problem. Every step required to construct the problem solution can be defined in program commands that produce intermediate output. This paper advocates improved problem solving skills through using a programming tool. MatLab created by MathWorks, is an interactive numerical computing environment and programming language. It is a matrix-based system that easily lends itself to matrix manipulation, and plotting of functions and data. MatLab can be used as an interactive command line or a sequence of commands that can be saved in a file as a script or named functions. Prior programming experience is not required to use MatLab commands. The GNU Octave, part of the GNU project, a free computer program for performing numerical computations, is comparable to MatLab. MatLab visual and command programming are presented here.

Shaykhian, Linda H.↗

Spectral Galerkin solver for intense beam Vlasov equilibria in nonlinear constant focusing channels

Here a numerical method is described for producing stationary solutions of the Vlasov-Poisson system describing a relativistic charged-particle beam in a constant focusing accelerator channel, confined transversely by a general (linear or nonlinear) focusing potential. The method utilizes a variant of the spectral Galerkin algorithm to solve a nonlinear partial differential equation (PDE) in two degrees of freedom for the beam space charge potential in equilibrium. Numerical convergence with an increasing number of computed spectral modes is investigated for several benchmark problems. Preservation of the stationary phase space density is verified using a strongly nonlinear focusing channel based on the Integrable Optics Test Accelerator at Fermi National Accelerator Laboratory.

43 PARTICLE ACCELERATORS↗

Fluid-structure interaction solver for transient dynamics of fracturing media

Computer implemented methods include providing a numerical simulation solid mesh comprising a plurality of solid nodes and a numerical simulation fluid mesh comprising a plurality of fluid nodes, and performing a numerical simulation with the solid mesh and fluid mesh including: determining a fluid-solid interaction between the fluid mesh and the solid mesh at a time step by (i) immersing boundaries of the solid and fluid meshes, (ii) computing current corrected fluid and solid velocities using an interaction parameter applied to an intermediate fluid-solid relative velocity, and (iii) computing a fluid-solid interaction force based on the current corrected fluid and solid velocities.

Rougier, Esteban↗

A block iterative LU solver for weakly coupled linear systems

A hybrid technique, called the block iterative LU solver, is proposed for solving the linear equations resulting from a finite element numerical analysis of certain fluid dynamics problems where the equations are weakly coupled between distinct sets of variables. Either the block Jacobi iterative method or the block Gauss-Seidel iterative solver is combined with LU decomposition.

Cooke, C. H.↗

Hundness versus Mottness in a three-band Hubbard–Hund model: On the origin of strong correlations in Hund metals☆

Hund metals are multi-orbital systems with moderate Coulomb interaction, U, among charges and sizeable Hund’s rule coupling, J(<U), that aligns the spins in different orbitals. They show strong correlation effects, like very low Fermi-liquid coherence scales and intriguing incoherent transport regimes, resulting in bad metallic behavior. But to what extent are these strong correlations governed by Mottness, i.e. the blocking of charge fluctuations close to a Mott insulator transition (MIT) induced by U, or by Hundness, a new route towards strong correlations induced by J? To answer this question, we study the full phase diagram of a degenerate three-band Hubbard–Hund model on a Bethe lattice at zero temperature using single-site dynamical mean-field theory and the numerical renormalization group as efficient real-frequency multi-band impurity solver. Hund metal behavior occurs in this minimal model for a filling close to n{sub d}=2, moderate U and sizeable J, the “Hund-metal regime”. In particular, strong correlations manifest themselves there by an unusually low quasiparticle weight. Generalizing previous results on this model, we show that “spin–orbital separation” (SOS) is a generic Hund’s-coupling-induced feature in the whole metallic regime of the phase diagram for 1<n{sub d}<3 and sizeable J. There orbital screening always occurs at much higher energies than spin screening below which Fermi-liquid behavior sets in. The low quasiparticle weight can then be directly explained in terms of the Hund’s-coupling-reduced Fermi-liquid scale. We carefully analyze the effect of J (Hundness), and the effect of the MIT at n{sub d}=2 and n{sub d}=3 (Mottness) on the energy scales and the nature of SOS. In the Hund-metal regime, far from any MIT, Hundness is shown to be the key player to induce strong correlations by constraining the spin rather than the charge dynamics. There, physical properties are governed by a broad incoherent energy regime of SOS where intriguing Hund metal physics occurs: large, almost unscreened spins are coupled to screened orbital degrees of freedom. With increasing proximity to an MIT correlations are further enhanced and the Fermi-liquid scale is further reduced. However, in the Hund-metal regime, this effect of Mottness is minor. In contrast, very close to the MIT at n{sub d}=2, the incoherent spin–orbital separation regime is strongly downscaled and becomes negligibly small, whereas Mottness – the localization of charges – becomes dominant in inducing strong correlations. Close to the MIT at n{sub d}=3, the SOS regime widens up because the orbital degrees of freedom get blocked by the formation of an S=(3/2) impurity spin, but its nature changes: the orbital and spin dynamics get decoupled. Our results confirm Hundness as a distinct mechanism towards strong correlations in the normal state of Hund metals, leading to various interesting implications for the nature of electronic transport.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Merging Intelligent Systems Technologies with CFD Analysis Strategies: Prototype Development at NASA Ames

A collaborative team of researchers from fields of Computational Fluid Dynamics (CFD), fluid physics, computer architectures, and computer science and knowledge engineering have begun work on a prototype system that addresses several of industry's concerns in using NASA-developed CFD codes as part of the design cycle. A major problem exists in the application of CFD technologies within the aeronautics design cycle due primarily to misunderstandings in the ranges of applicability of the various solver codes or turbulence models. Features that arise during the CFD solution process need to be discriminated and recognized as actual flow features with physical support in the geometry and flow conditions of the problem being solved, or as numerical or non-physical errors arising from mis-application of solver code and its parameters, gridding strategies, or discretization. interpolations. The fundamental concept is to develop an intelligent computational system that can accept the engineer's definition of the problem and construct an optimal CFD solution. To do this requires capturing both the knowledge of how to apply the various CFD tools and how to adapt the application of those tools to flow structures as they evolve during the flow simulation. Embedded within this adaptive system approach is the additional desire to automatically identify and quantify the quality of resolution of the pertinent flow structures, be they genuine or error-induced, and then to adjust the solution strategy accordingly. This paper discusses the status of that prototyping effort.

Thompson, David E.↗

Acoustics and Aeroperformance of Nozzles With Screwdriver Shaped and Axisymmetric Plugs

The recent experimental and numerical tests of corrugated nozzles have shown some acoustic and thrust benefits relative to traditional round nozzles. For example, a Bluebell nozzle which was obtained by 3D nozzle design incorporating a corrugated cross section nozzle shape with a sinusoidal lip line nozzle edge, can provide an acoustic benefit up to 4dB with about a 1% thrust augmentation. In references, this effect was explained as being the result of the corrugated design producing more efficient mixing of the exhausted jet with ambient air. Based on this argument, the authors have proposed the application of this concept for a centerbody (plug) which can form several vortices downstream from the centerbody. Several different corrugated designs are proposed and described in detail in this paper. The main design is a Screwdriver shaped centerbody or plug (SCR) which was tested experimentally and numerically. The acoustic tests were conducted in the anechoic chamber of the Central AeroHydrodynamics Institute (TsAGI, Moscow) under Civilian Research and Development Foundation (CRDF) grant. These experiments have shown an essential acoustic benefit of about 10-13% with the application of the co-annular nozzles by comparison with the reference round nozzle with the same mass flow rate. However, the expected acoustic benefits with the application of the 4-petal Screwdriver shaped centerbody were not obtained by comparison with the reference axisymmetric centerbody (CON) having the same length and the same cross section areas at the same distance from the nozzle throat. Moreover, for some angles (Theta = 60 deg and 90 deg) noise increase was observed (about 1-3%). These tests will be continued with the goal of obtaining better acoustic results. In particular, acoustic characteristics are hoped to be improved by moving t lie centerbody into the nozzle and using penetrable walls for the SCR and/or for the main nozzle. Preliminary results for such approach are very promising. Aeroperformance effects were analyzed numerically. The nozzle thrust, calculations were based on a full Navier--Stokes equations solver (NSE), and both full and marching Euler codes: CFL3D, CRAFT, and Krayko-Godunov . Grid preparation and its optimization were conducted using GRIDGEN and our own codes. The general conclusion of this numerical analysis is some thrust loss with the application of SCR design (about 1 - 1.5%). But, again, some constructive features of SCR design give some promising perspectives for its application in aviation and domestic industries.

Gilinsky, M.↗

Computational Analysis of a Boundary Layer Ingesting Tailcone Thruster Configuration Using the LAVA Curvilinear Solver

Boundary layer ingesting (BLI) propulsion systems are one of the many technologies currently under investigation within the aerospace community to meet NASA's Advanced Air Transport Technology project goal of a sustainable future in aviation. To that end, an experimental campaign was conducted in the National Transonic Facility (NTF) wind tunnel to study the propulsion-airframe integration effects of the Boundary Layer Ingesting Tailcone System installed in the aft portion of a 2.7% scale design of the NASA Common Research Model (CRM). This test was accompanied by a numerical simulation effort using the Launch, Ascent, and Vehicle Aerodynamics (LAVA) curvilinear solver to validate the applicability of current best-practice Reynolds-Averaged Navier-Stokes (RANS) models in accurately predicting the relevant flow physics in free-air. A total of 205 steady RANS simulations were performed using LAVA, covering the range of conditions present in the NTF test matrix. Flow conditions ranged from 5- to 15-million Reynolds number, Mach numbers of 0.75, 0.80 and 0.85, and angles of attack between -3 and 4 degrees. Four different mass-flow plugs were also tested to assess propulsor operating condition effects. Sensitivity to these flow conditions is analyzed in detail, especially in terms of the nacelle flow distortion which was the main subject of the study. Aerodynamic load coefficients, aft body boundary layer measurements and fuselage pressure tap results are also discussed, providing a wide range of validation results from this experimental campaign. This validation effort shows that current best-practices using RANS models within the LAVA curvilinear solver flow solver are well suited for predicting the inlet flow distortion characteristics of novel aircraft configurations employing BLI at the aft end of the fuselage. This talk will cover the combined experimental and numerical efforts that resulted from this BLI-focused investigation.

ARMD↗

PETSc/TAO developments for GPU-based early exascale systems

The Portable Extensible Toolkit for Scientific Computation (PETSc) library provides scalable solvers for nonlinear time-dependent differential and algebraic equations and for numerical optimization via the Toolkit for Advanced Optimization (TAO). PETSc is used in dozens of scientific fields and is an important building block for many simulation codes. During the U.S. Department of Energy’s Exascale Computing Project, the PETSc team has made substantial efforts to enable efficient utilization of the massive fine-grain parallelism present within exascale compute nodes and to enable performance portability across exascale architectures. We recap some of the challenges that designers of numerical libraries face in such an endeavor, and then discuss the many developments we have made, which include the addition of new GPU backends, features supporting efficient on-device matrix assembly, better support for asynchronicity and GPU kernel concurrency, and new communication infrastructure. In conclusion, we evaluate the performance of these developments on some pre-exascale systems as well as the early exascale systems Frontier and Aurora, using compute kernel, communication layer, solver, and mini-application benchmark studies, and then close with a few observations drawn from our experiences on the tension between portable performance and other goals of numerical libraries.

Exascale Computing Project (ECP)↗