Search NASA⌕ Search

SEARCH · Search NASA

Results for “NUMERICAL FLUX FUNCTION”

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 145 records · Page 8
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Performance Improvements of the Griffin Solvers in FY24

The Griffin code is a MOOSE-based reactor physics application jointly developed by Idaho National Laboratory and Argonne National Laboratory under the Department of Energy Office of Nuclear Energy Nuclear Energy Advanced Modeling and Simulation Program. This fiscal year, we have made significant efforts to improve the performance of transport solver options and cross-section generation for the efficient use of Griffin in advanced reactor applications. For the HFEM-PN solver, the residual evaluations of HFEM kernels were optimized by utilizing the pre- computed averaged cross sections for individual elements. Numerical integration involving the evaluation of basis functions at quadrature points was bypassed by facilitating precomputed element mass matrices for response matrices. Red-black iterations were improved by introducing a new generalized minimum residual based solver. The memory usage of response matrix storage was significantly reduced by applying basis function rotations on interfaces and calculating volumetric odd-parity moments on the fly. Additionally, the adjoint flux and transient calculation capabilities of the HFEM-PN solver were successfully implemented and verified using the TWIGL benchmark problem. For the DFEM-SN solver, memory footprint and computation time were significantly reduced by not treating angular flux vectors as the MOOSE nonlinear system vectors. Specifically for IQS, scalar adjoint weighting was introduced to further eliminate angular adjoint flux storage in the MOOSE auxiliary system. It was demonstrated through the three-dimensional Advanced Burner Test Reactor core problem that the memory usage for transient calculations with the IQS method was reduced by over 7.5× compared to before the optimizations. For the self-shielding application programming interface, a new double-heterogeneity treatment method, named the Bell Function-Based Analytic Two-Region Slowing Down Method, was developed to efficiently flux-volume homogenize TRISO particles with the matrix. Additionally, optimizations were made to hyper- fine group (HFG) slowing down calculations by pretabulating collision probability coefficients and grouping isotopes, significantly reducing the computational time for calculating scattering sources per HFG. Lastly, the pin power reconstruction module was extended to account for temporal behavior in a microreactor analysis problem, specifically for a control drum transient. Verification tests for each of these improvements demonstrated significant performance enhancements and memory reduction.

22 - GENERAL STUDIES OF NUCLEAR REACTORS↗
NASA NTRS1985

A new convective parameterization scheme for global models

A convective adjustment scheme is presented for parameterization of the vertical temperature and moisture structure for large-scale numerical weather models. The temperature and moisture structure is adjusted toward observed quasi-equilibrium structures to maintain a realistic vertical temperature and moisture structure in grid-scale global models during convection. Two reference thermodynamic structures are modeled analytically to account for shallow and deep convection. The thermodynamic structure obtained is permitted to change as a function of large-scale advective terms, radiation and surface fluxes. The reference profiles are applied when an energy integral provides a description of the shallow or deep convective processes occurring. Comparisons with GATE, BOMEX and ATEX data and a 10 day forecast for weather over the tropical ocean illustrate the climate and precipitation forecast accuracy improvements available with the adjustment scheme relative to standard parameterization schemes.

Miller, M. J.↗
NASA NTRS1979

Stellar convection 2: A multi-mode numerical solution for convection in spheres

The convective flow of a self gravitating sphere of Boussinesq fluid for small Reynolds and Peclet numbers is numerically determined. The decomposition of the equations of motion into modes is reviewed and a relaxation method is developed and presented to compute the solutions to these equations. The stable equilibrium flow for a Rayleigh number of 10 to the 4th power and a Prandtl number of 10 is determined. The 2 and 3 dimensional spectra of the kinetic and thermal energies and the convective flux as a function of wavelengths are calculated in terms of modes. The anisotropy of the flow as a function of wavelength is defined.

Marcus, P. S.↗
NASA NTRS1987

A generalized procedure for constructing an upwind based TVD scheme

A generalized formulation for constructing second- and higher-order accurate TVD (total variation diminishing) schemes is presented. A given scheme is made TVD by limiting antidiffusive flux differences with some linear functions, so-called limiters. The general idea of the formulation and its mathematical proof of Harten's TVD conditions is shown by applying the Lax-Wendroff method to scalar nonlinear equations and a constant-coefficient system of conservation laws. For the system of equations, several definitions are derived for the argument used in the limiter function and present their performance in numerical experiments. The formulation is extended to the nonlinear system. It is demonstrated that the present procedure can easily convert existing central or upwind, and second- or higher-order differencing schemes to preserve monotonicity and yield physically admissible solutions. The formulation is simple mathematically as well as numerically; both matrix-vector multiplication and Riemann solver are avoided. Although the notion of TVD is based on the initial value problem, application to the steady Euler equations of the formulation is also made.

Liou, Meng-Sing↗
NASA NTRS1987

A generalized procedure for constructing an upwind-based TVD scheme

A generalized formulation for constructing second- and higher-order accurate TVD (total variation diminishing) schemes is presented. A given scheme is made TVD by limiting antidiffusive flux differences with some nonlinear functions, so-called limiters. The general idea of the formulation and its mathematical proof of Harten's TVD conditions is shown by applying the Lax-Wendroff method to a scalar nonlinear equation and constant-coefficient system of conservation laws. For the system of equations, several definitions are derived for the argument used in the limiter function and present their performance to numerical experiments. Then the formulation is formally extended to the nonlinear system of equations. It is demonstrated that use of the present procedure allows easy conversion of existing central or upwind, and second- or higher-order differencing schemes so as to preserve monotonicity and to yield physically admissible solutions. The formulation is simple mathematically as well as numerically; neither matrix-vector multiplication nor Riemann solver is required. Roughly twice as much computational effort is needed as compared to conventional scheme. Although the notion of TVD is based on the initial value problem, application to the steady Euler equations of the formulation is also made. Numerical examples including various ranges of problems show both time- and spatial-accuracy in comparison with exact solutions.

Liou, Meng-Sing↗
NASA NTRS1993

Propagation of electromagnetic waves parallel to the magnetic field in the nightside Venus ionosphere

The propagation of electromagnetic waves parallel to the magnetic field in the nightside Venus ionosphere is presented in a theoretical and numerical analysis. The model assumes a source of electromagnetic radiation in the Venus atmosphere, such as that produced by lightning. Specifically addressed is wave propagation in the altitude range z = 130-160 km at the four frequencies detectable by the Pioneer Venus Orbiter Electric Field Detector: 100 Hz, 730 Hz, 5.4 kHz, and 30 kHz. Parameterizations of the wave intensities, peak electron density, and Poynting flux as a function of magnetic field are presented. The waves are found to propagate most easily in conditions of low electron density and high magnetic field. The results of the model are consistent with observational data.

Huba, J. D.↗
NASA NTRS1995

Theoretical overview and modeling of the sodium and potassium atmospheres of mercury

A general theoretical overview for the sources, sinks, gas-surface interactions, and transport dynamics of sodium and potassium in the exospheric atmsophere of Mercury is given. Information for these four factors, which control the spatial distribution of these two alkali-group gases about the planet, is incorporated in numerical models. The spatial nature and relative importance of the initial source atom atmosphere and the ambient (ballistic hopping) atom atmosphere are then examined and are shown to be controlled and coupled to a great extent by the extremely large and variable solar radiation acceleration experienced by sodium and potassium as they resonantly scatter solar photons. The lateral (antisunward) transport rate of thermally accommodated sodium and potassium ambient atoms is shown to be driven by the solar radiation acceleration and, over a significant portion of Mercury's orbit about the Sun, is sufficiently rapid to be competitive with the short photoionization lifetimes for these atoms when they are located on the summit surface near or within about 30 deg of the terminator. The lateral transport rate is characterized by a migration time determined by model calculations for an ensemble of atoms initially starting at a point source on the surface (i.e., a numerical spacetime dependent Green's function). Four animations for the spacetime evolution of the sodium (or potassium) atmosphere produced by a point source on the surface are presented on a videotape format. For extended surface sources for sodium and potassium, the local column density is determined by competition between the photoionization lifetimes and the lateral transport times of atoms originating from different surface source locations. Sodium surface source fluxes (referenced to Mercury at perihelion) that are required on the sunlit hemisphere to reproduce the typically observed several megarayleighs of D2 emission-line brightness and the inferred column densities of 1-2 x 10(exp 11) atoms per sq cm range from approximately 2-5 x 10(exp 7) atoms/sq cm/sec. The sodium model is applied to study observational data that document an anticorrelation in the average sodium column density and solar radiation acceleration. Lateral transport driven by the solar radiation acceleration is shown to produce this behavior for combinations of different sources and surface accomodation coefficients. The best fit model fits to the observational data require a significant degree of thermal accommodation of the ambient sodium atoms to the surface and a source rate that decreases as an inverse power of 1.5 to 2 in heliocentric distance.

Smyth, William H.↗
DOE OSTI2026

An immersed interface method for microstructure-scale electrochemical battery models: numerical formulation and performance portable implementation

We present the numerical formulation, verification, and performance portable implementation of an immersed interface method for microstructure scale electrochemical modeling of batteries. The innovation in this approach is the resolution of chemical species and electrostatic potential discontinuities at active interfaces without the use of interface conforming unstructured grids. A unified formulation on Cartesian grids for all domains (electrodes and electrolyte) is used with interfacial flux conditions applied using volume fraction or “color” function gradients. We have developed one dimensional and two dimensional test cases with analytic solutions for electrochemical modeling using which we verified the consistency and accuracy of our scheme. Our solver is also validated against solutions from a macroscale model and an unstructured multi-subdomain solver for a full lithium ion cell. We then demonstrated the utility of our solver on an image-based complex battery electrode microstructure at high charging rate. Our technique also exhibits good scalability on distributed memory architectures using central processing units (CPU), with problem sizes up to 1.8 billion degrees of freedom and with 5400 ranks. Initial performance studies of our open-source performance portable solver showed about 70 times speed up using a graphics processing unit (GPU) compared to single compute core for a problem with 4 million cells.

25 ENERGY STORAGE↗
NASA NTRS1979

Utility of a finite element solution algorithm for initial-value problems

The Galerkin criterion within a finite element Weighted Residuals formulation is employed to establish an implicit solution algorithm for an initial-value partial differential equation. Numerical solutions of a transient parabolic and a hyperbolic equation, obtained using linear, quadratic and two cubic finite element basis functions, are employed to quantize accuracy and confirm and refine theoretical convergence rate estimates. The linear basis algorithm for the hyperbolic equation displays excellent accuracy on a coarse computational grid and a high-order convergence rate with discretization refinement. Good accuracy and a strong convergence rate in surface flux are determined for a nonhomogeneous Neumann boundary constraint applied to a parabolic equation. The results amply demonstrate the impact of the nondiagonal finite element initial-value matrix structure on solution accuracy and/or convergence rate.

Baker, A. J.↗
NASA NTRS1993

Multi-resolution analysis for ENO schemes

Given a function u(x) which is represented by its cell-averages in cells which are formed by some unstructured grid, we show how to decompose the function into various scales of variation. This is done by considering a set of nested grids in which the given grid is the finest, and identifying in each locality the coarsest grid in the set from which u(x) can be recovered to a prescribed accuracy. We apply this multi-resolution analysis to Essentially Non-oscillatory Schemes (ENO) schemes in order to reduce the number of numerical flux computations which is needed in order to advance the solution by one time-step. This is accomplished by decomposing the numerical solution at the beginning of each time-step into levels of resolution, and performing the computation in each locality at the appropriate coarser grid. We present an efficient algorithm for implementing this program in the one-dimensional case; this algorithm can be extended to the multi-dimensional case with cartesian grids.

Harten, Ami↗
NASA NTRS2004

The Investigation Of Carbon Contamination And Sputtering Effects Of Xenon Ion Thrusters

The Electro-Physics Branch of the NASA Glenn Research Center investigates the effect of atomic oxygen, environmental durability of high performance power materials and surfaces, and low earth orbit. One of its current projects involves the analysis of ion thrusters. Ion thrusters are devices that initiate a beam of ions to a target area. The type of ion thruster that I have been working with this Summer of 2004 emits positively charged Xenon (Xe(+)) atoms through two grids, the screen grid and the accelerator grid, after it enters an ionization chamber. Insulators are used to mechanically hold and separate these two grids. A propellant isolator, an instrument that closely resembles insulators, is placed in front of the ionization chamber. Both the insulator and isolator are made with a ceramic compound and filled with insulating beads. The main difference between the two devices is that the propellant isolator allows gas to flow through, in this case, the gas is Xe(+) and the insulators do not. In order to avoid carbon deposits and other contaminating chemicals to settle on the insulators and propellant isolator, a metal shadow shield is placed around them. These shadow shields function as a protectant and can be shaped in numerous configurations. Part of my job responsibility this summer is to investigate the effectiveness of different shadow shields that are utilized on three different ion engines: the NSTAR (NASA Solar Electric Propulsion Technology Application Readiness), JIMO (Jupiter Icy Moons Orbiter), and NEXIS (Nuclear Electric Xenon Ion System). Using calculus and other mathematical tactics, I was asked to find the total flux of carbon contamination that was able to pass the protectant shadow shield. I familiarized myself with the software program, MathCad2004, to help perform some mathematical computations such as complex integration. Another method of studying the probability of contamination is by experimental simulation. After attaining the precise parameters of the actual shadow shields, I created replicas of three types of shadow shielding to be used to undergo testing. It will be placed in a machine that produces carbon atoms at a high temperature of 200 C. or beam is aimed at a targeted material. As a result of this collision, atoms and other particles are ejected out of the target surface. Another part of my internship consisted of research on sputter ejection, or the angle distribution of sputtered material. This research entailed finding the past results of sputter ejection investigation as well as creating another type of mock simulation. Other minor projects include calculating the path of Xe(+) gas through the insulating beads of the isolators and assisting my mentor in collecting data for his paper for the Joint Propulsion Conference & Exhibit to be held July 11-14,2004 in Fort Lauderdale, Florida.

Prak, Moline K.↗
NASA NTRS2006

Development of a Countermeasure to Mitigate Postflight Locomotor Dysfunction

Astronauts returning from space flight experience locomotor dysfunction following their return to Earth. Our laboratory is currently developing a gait adaptability training program that is designed to facilitate recovery of locomotor function following a return to a gravitational environment. The training program exploits the ability of the sensorimotor system to generalize from exposure to multiple adaptive challenges during training so that the gait control system essentially learns to learn and therefore can reorganize more rapidly when faced with a novel adaptive challenge. Evidence for the potential efficacy of an adaptive generalization gait training program can be obtained from numerous studies in the motor learning literature which have demonstrated that systematically varying the conditions of training enhances the ability of the performer to learn and retain a novel motor task. These variable practice training approaches have been used in applied contexts to improve motor skills required in a number of different sports. The central nervous system (CNS) can produce voluntary movement in an almost infinite number of ways. For example, locomotion can be achieved with many different combinations of joint angles, muscle activation patterns and forces. The CNS can exploit these degrees of freedom to enhance motor response adaptability during periods of adaptive flux like that encountered during a change in gravitational environment. Ultimately, the functional goal of an adaptive generalization countermeasure is not necessarily to immediately return movement patterns back to normal. Rather the training regimen should facilitate the reorganization of available sensory and motor subsystems to achieve safe and effective locomotion as soon as possible after long duration space flight. Indeed, this approach has been proposed as a basic feature underlying effective neurological rehabilitation. We have previously confirmed that subjects participating in an adaptive generalization training program using a variety of visuomotor distortions and throwing as the dependent measure can learn to enhance their ability to adapt to a novel sensorimotor environment (Roller et al., 2001). Importantly, this increased adaptability was retained even one month after completion of the training period. Adaptive generalization has been observed in a variety of other tasks requiring sensorimotor transformations including manual control tasks and reaching (Bock et al., 2001, Seidler, 2003) and obstacle avoidance during walking (Lam and Dietz, 2004). Taken together, the evidence suggests that a training regimen exposing crewmembers to variation in locomotor conditions, with repeated transitions among states, may enhance their ability to learn how to reassemble appropriate locomotor patterns upon return from microgravity. We believe exposure to this type of training will extend crewmembers locomotor behavioral repertoires, facilitating the return of functional mobility after long duration space flight. In other words, our proposed training protocol will compel subjects to develop new behavioral solutions under varying sensorimotor demands. Over time subjects will learn to create appropriate locomotor solution more rapidly enabling acquisition of mobility sooner after long-duration space flight. A gait adaptability training program can be superimposed on nominal treadmill exercise activities thus ensuring that no additional crew time is required to perform this type of training regimen and that it can be implemented with current in-flight exercise systems available on the International Space Station.

Bloomberg, J. J.↗
NASA NTRS1993

Numerical simulation of rarefied flow through a slit. I - Direct simulation Monte Carlo results

The pressure-driven flow of a rarefied monatomic gas through a two-dimensional slit is simulated using the direct simulation Monte Carlo technique. Of particular interest is the change in flow field structure as pressure ratio and Knudsen number are varied. Comparisons are made to quantify the limits of validity of free-molecular theory and approximate, nearly free-molecular iterative methods. Also addressed is the sensitivity of the numerical solutions to grid structure and boundary conditions. The free-molecular theory is found to predict quantitative flow field properties (e.g., centerline velocities or downstream flux profiles) reasonably well for large finite Knudsen number with the error dependent on the pressure ratio. The nearly free-molecular corrections are shown to have limited range of applicability. A previously derived parameter is found to correlate total mass flux well as a function of pressure ratio and Knudsen number over a large portion of the transitional regime.

Wadsworth, D. C.↗
NASA NTRS2015

Badhwar-O'Neill 2014 Galactic Cosmic Ray Flux Model Description

The Badhwar-O'Neill (BON) Galactic Cosmic Ray (GCR) model is based on GCR measurements from particle detectors. The model has mainly been used by NASA to certify microelectronic systems and the analysis of radiation health risks to astronauts in space missions. The BON14 model numerically solves the Fokker-Planck differential equation to account for particle transport in the heliosphere due to diffusion, convection, and adiabatic deceleration under the assumption of a spherically symmetric heliosphere. The model also incorporates an empirical time delay function to account for the lag of the solar activity to reach the boundary of the heliosphere. This technical paper describes the most recent improvements in parameter fits to the BON model (BON14). Using a comprehensive measurement database, it is shown that BON14 is significantly improved over the previous version, BON11.

O'Neill, P. M.↗
NASA NTRS1989

Solitary waves in mantle plumes

Numerical calculations in two dimensions show that solitary wave disturbances can propagate along thermal plumes in a homogeneous, viscous fluid with a thermally activated rheology similar to the mantle. Comparison of the numerical results with analytic results from an idealized two-fluid model indicates that the two-fluid model can be applied to solitary waves in the thermally activated plume. A two-fluid model of solitary waves on cylindrical conduits is accordingly applied to mantle plumes to estimate the propagation speeds, time durations, and pulse lengths of solitary waves in mantle plumes as a function of background mantle viscosity, plume flux, plume density deficit, plume viscosity, and the volume of material transported by the solitary wave. It is inferred that mantle plumes may be eposodic via solitary waves and that these disturbances might account for observed 0(10 m.y.) variations in the durations and spacings of episodes of enhanced hotspot volcanic activity. Solitary waves in mantle plumes could be generated by interactions among plumes and between plumes and the large-scale, time-dependent mantle circulation.

Schubert, G.↗
NASA NTRS2008

Tungsten and Barium Transport in the Internal Plasma of Hollow Cathodes

The effect of tungsten erosion, transport and redeposition on the operation of dispenser hollow cathodes was investigated in detailed examinations of the discharge cathode inserts from an 8200 hour and a 30,352 hour ion engine wear test. Erosion and subsequent re-deposition of tungsten in the electron emission zone at the downstream end of the insert reduces the porosity of the tungsten matrix, preventing the flow of barium from the interior. This inhibits the interfacial reactions of the barium-calcium-aluminate impregnant with the tungsten in the pores. A numerical model of barium transport in the internal xenon discharge plasma shows that the barium required to reduce the work function in the emission zone can be supplied from upstream through the gas phase. Barium that flows out of the pores of the tungsten insert is rapidly ionized in the xenon discharge and pushedback to the emitter surface by the electric field and drag from the xenon ion flow. Thisbarium ion flux is sufficient to maintain a barium surface coverage at the downstream endgreater than 0.6, even if local barium production at that point is inhibited by tungsten deposits. The model also shows that the neutral barium pressure exceeds the equilibrium vapor pressure of the impregnant decomposition reaction over much of the insert length,so the reactions are suppressed. Only a small region upstream of the zone blocked by tungsten deposits is active and supplies the required barium. These results indicate that hollowcathode failure models based on barium depletion rates in vacuum dispenser cathodes are very conservative.

Hollow cathodes↗
NASA NTRS1993

A new flux conserving Newton's method scheme for the two-dimensional, steady Navier-Stokes equations

A new numerical method is developed for the solution of the two-dimensional, steady Navier-Stokes equations. The method that is presented differs in significant ways from the established numerical methods for solving the Navier-Stokes equations. The major differences are described. First, the focus of the present method is on satisfying flux conservation in an integral formulation, rather than on simulating conservation laws in their differential form. Second, the present approach provides a unified treatment of the dependent variables and their unknown derivatives. All are treated as unknowns together to be solved for through simulating local and global flux conservation. Third, fluxes are balanced at cell interfaces without the use of interpolation or flux limiters. Fourth, flux conservation is achieved through the use of discrete regions known as conservation elements and solution elements. These elements are not the same as the standard control volumes used in the finite volume method. Fifth, the discrete approximation obtained on each solution element is a functional solution of both the integral and differential form of the Navier-Stokes equations. Finally, the method that is presented is a highly localized approach in which the coupling to nearby cells is only in one direction for each spatial coordinate, and involves only the immediately adjacent cells. A general third-order formulation for the steady, compressible Navier-Stokes equations is presented, and then a Newton's method scheme is developed for the solution of incompressible, low Reynolds number channel flow. It is shown that the Jacobian matrix is nearly block diagonal if the nonlinear system of discrete equations is arranged approximately and a proper pivoting strategy is used. Numerical results are presented for Reynolds numbers of 100, 1000, and 2000. Finally, it is shown that the present scheme can resolve the developing channel flow boundary layer using as few as six to ten cells per channel width, depending on the Reynolds number.

Scott, James R.↗
NASA NTRS1973

Measurement of the lunar photoelectron layer in the geomagnetic tail

Stable photoelectron fluxes (with energies between 40 and 200 eV) were observed by the Apollo 15 Charged Particle Lunar Environment Experiment (CPLEE). These observations, made in the magnetotail under near vacuum conditions, are compared with numerically calculated photoemission spectra to determine the approximate potential difference between ground and CPLEE'S apertures (a distance of 26 cm). Numerically calculated density and potential distributions are compared with the measured values to provide an estimate of the photoelectron yield function of the dust layer covering the moon.

Reasoner, D. L.↗