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 19 records
All sourcesNASA NTRSDOE OSTI
NASA NTRS1989

Some notes on shock resolving flux functions. Part 1: Stationary characteristics

Numerical flux functions for solving the Euler equations using exact and/or approximate solutions of the Riemann problem of gasdynamics are discussed. Under certain restrictive conditions, schemes using these flux functions produce systems of equations which can exhibit a single degree of freedom. In some instances, the solutions represented by this degree of freedom are unstable to perturbations. This local instability can seriously degrade the temporal convergence of numerical schemes. This point is demonstrated by numerical example.

Barth, Timothy J.↗
NASA NTRS1981

On the relation between upwind-differencing schemes of Godunov, Engquist-Osher and Roe

The upwind differencing first order schemes of Godunov, Engquist-Osher and Roe are discussed on the basis of the inviscid Burgers equations. The differences between the schemes are interpreted as differences between the approximate Riemann solutions on which their numerical flux functions are based. Special attention is given to the proper formulation of these schemes when a source term is present. Second order two step schemes, based on the numerical flux functions of the first order schemes are also described. The schemes are compared in a numerical experiment, and recommendations on their use are included.

Vanleer, B.↗
NASA NTRS1993

Grid-independent upwind scheme for multidimensional flow

Recent advances in the development of a grid-independent finite volume scheme for the Euler equations of gas dynamics are described. In the proposed method, flowfield gradient data are reconstructed locally (on a triangle) using five elementary planar waves, and an upwind numerical flux function for grid-oblique waves is developed to model the effect of the passage of these waves on the data in a cell. Numerical examples for several two-dimensional test problems are included. These results show high wave resolution and nearly monotone strong-wave transitions.

Parpia, Ijaz H.↗
NASA NTRS1983

Numerical approximations to nonlinear conservation laws with locally varying time and space grids

Numerical approximations to the initial value problem for nonlinear systems of conservation laws are considered. The considered system is said to be hyperbolic when all eigenvalues of every real linear combination of the Jacobian matrices are real. Solutions may develop discontinuities in finite time, even when the initial data are smooth. In the investigation, explicit finite difference methods which use locally varying time grids are considered. The global CFL restriction is replaced by a local restriction. The numerical flux function is studied from a finite volume viewpoint, and a differencing technique is developed at interface points between regions of distinct time increments.

Osher, S.↗
NASA NTRS1985

On the implementation of a class of upwind schemes for system of hyperbolic conservation laws

The relative computational effort among the spatially five point numerical flux functions of Harten, van Leer, and Osher and Chakravarthy is explored. These three methods typify the design principles most often used in constructing higher than first order upwind total variation diminishing (TVD) schemes. For the scalar case the difference in operation count between any two algorithms may be very small and yet the operation count for their system counterparts might be vastly different. The situation occurs even though one starts with two different yet equivalent representations for the scalar case.

Yee, H. C.↗
NASA NTRS1989

Sonic-point capturing

A prototype scheme that produces perfectly smooth transonic solutions to nozzle-flow problems is derived and tested. The basic upwind scheme is described as well as satisfying the entropy condition, treatment of the source term, and numerical verification. The analysis yielded a numerical flux function for use near a sonic point, which is based on a full model of a transonic expansion wave, and a matched treatment for the source term.

Van Leer, Bram↗
NASA NTRS1991

Numerical flux formulas for the Euler and Navier-Stokes equations. 2: Progress in flux-vector splitting

The accuracy of various numerical flux functions for the inviscid fluxes when used for Navier-Stokes computations is studied. The flux functions are benchmarked for solutions of the viscous, hypersonic flow past a 10 degree cone at zero angle of attack using first order, upwind spatial differencing. The Harten-Lax/Roe flux is found to give a good boundary layer representation, although its robustness is an issue. Some hybrid flux formulas, where the concepts of flux-vector and flux-difference splitting are combined, are shown to give unsatisfactory pressure distributions; there is still room for improvement. Investigations of low diffusion, pure flux-vector splittings indicate that a pure flux-vector splitting can be developed that eliminates spurious diffusion across the boundary layer. The resulting first-order scheme is marginally stable and not monotone.

Coirier, William J.↗
NASA NTRS1991

Numerical flux formulas for the Euler and Navier-Stokes equations. II - Progress in flux-vector splitting

The accuracy of various numerical flux functions for the inviscid fluxes when used for Navier-Stokes computations is studied. The flux functions are benchmarked for solutions of the viscous, hypersonic flow past a 10 degree cone at zero angle of attack using first order, upwind spatial differencing. The Harten-Lax/Roe flux is found to give a good boundary layer representation, although its robustness is an issue. Some hybrid flux formulas, where the concepts of flux-vector and flux-difference splitting are combined, are shown to give unsatisfactory pressure distributions; there is still room for improvement. Investigations of low diffusion, pure flux-vector splittings indicate that a pure flux-vector splitting can be developed that eliminates spurious diffusion across the boundary layer. The resulting first-order scheme is marginally stable and not monotone.

Coirier, William J.↗
NASA NTRS2001

Energy Stable Flux Formulas For The Discontinuous Galerkin Discretization Of First Order Nonlinear Conservation Laws

We consider the discontinuous Galerkin (DG) finite element discretization of first order systems of conservation laws derivable as moments of the kinetic Boltzmann equation. This includes well known conservation law systems such as the Euler For the class of first order nonlinear conservation laws equipped with an entropy extension, an energy analysis of the DG method for the Cauchy initial value problem is developed. Using this DG energy analysis, several new variants of existing numerical flux functions are derived and shown to be energy stable.

Barth, Timothy↗
NASA NTRS2015

Development of the US3D Code for Advanced Compressible and Reacting Flow Simulations

Aerothermodynamics and hypersonic flows involve complex multi-disciplinary physics, including finite-rate gas-phase kinetics, finite-rate internal energy relaxation, gas-surface interactions with finite-rate oxidation and sublimation, transition to turbulence, large-scale unsteadiness, shock-boundary layer interactions, fluid-structure interactions, and thermal protection system ablation and thermal response. Many of the flows have a large range of length and time scales, requiring large computational grids, implicit time integration, and large solution run times. The University of Minnesota NASA US3D code was designed for the simulation of these complex, highly-coupled flows. It has many of the features of the well-established DPLR code, but uses unstructured grids and has many advanced numerical capabilities and physical models for multi-physics problems. The main capabilities of the code are described, the physical modeling approaches are discussed, the different types of numerical flux functions and time integration approaches are outlined, and the parallelization strategy is overviewed. Comparisons between US3D and the NASA DPLR code are presented, and several advanced simulations are presented to illustrate some of novel features of the code.

CFD↗
NASA NTRS1994

Parallel implicit unstructured grid Euler solvers

A mesh-vertex finite volume scheme for solving the Euler equations on triangular unstructured meshes is implemented on a multiple-instruction/multiple-data stream parallel computer. An explicit four-stage Runge-Kutta scheme is used to solve two-dimensional flow problems. A family of implicit schemes is also developed to solve these problems, where the linear system that arises at each time step is solved by a preconditioned GMRES algorithm. Two partitioning strategies are employed: one that partitions triangles and the other that partitions vertices. The choice of the preconditioner in a distributed memory setting is discussed. All of the methods are compared both in terms of elapsed times and convergence rates. It is shown that the implicit schemes offer adequate parallelism at the expense of minimal sequential overhead. The use of a global coarse grid to further minimize this overhead is also investigated. The schemes are implemented on a distributed memory parallel computer, the Intel iPSC/860.

TRT-THEORETICAL↗
NASA NTRS1986

Thermal balance of the atmospheres of Jupiter and Uranus

Two-dimensional, radiative-convective-dynamical models of the visible atmospheres of Jupiter and Uranus are presented. Zonally-averaged temperatures and heat fluxes are calculated numerically as functions of pressure and latitude. In addition to radiative heat fluxes, the dynamical heat flux due to large-scale baroclinic eddies is included and is parametrized using a mixing length theory which gives heat fluxes similar to those of Stone. The results for Jupiter indicate that the internal heat flow is non-uniform in latitude and nearly balances the net radiative flux leaving the atmosphere. The thermal emission is found to be uniform in latitude in agreement with Pioneer and Voyager observations. Baroclinic eddies are calculated to transport only a small amount of the meridional heat flow necessary to account for the uniformity of thermal emission with latitude. The bulk of the meridional heat transfer is found to occur very deep in the stable interior of Jupiter as originally proposed by Ingersoll and Porco. The relative importance of baroclinic eddies vs. internal heat flow in the thermal balance of Uranus depends on the ratio of emitted thermal power to absorbed solar power. The thermal balance of Uranus is compared to that of Jupiter for different values of this ratio.

Friedson, A. J.↗
NASA NTRS1984

A shock-capturing finite element method

Finite element methods are developed for a one-dimensional singular perturbation problem involving a nonlinear flux function. In the limit of vanishing dissipation, the numerical flux is found to be almost identical to that due to Engquist and Osher, the only difference being in the treatment of transonic compression. This facilitates sharper resolution of stationary shocks, but is accompanied by loss of smoothness and monotonicity.

Hughes, T. J. R.↗
NASA NTRS2020

Studies in Astronomical Time Series Analysis: VII. An Enquiry Concerning Non-Linearity, the RMS-Mean Flux Relation, and log-Normal Flux Distributions

A broad and widely used class of stationary, linear, additive time series models can have statistical properties which many authors have asserted imply that the underlying process must be non-linear, non-stationary, multiplicative, or inconsistent with shot noise. This result is demonstrated with exact and numerical evaluation of the model flux distribution function and dependence of flux standard deviation on mean flux (here and in the literature called the rms-flux relation). These models can: (1) exhibit normal, log-normal or other flux distributions; (2) show linear or slightly non-linear rms-mean flux dependencies; as well as (3) match arbitrary second order statistics of the time series data. Accordingly the above assertions cannot be made on the basis of statistical time series analysis alone. Also idscussed are ambiguities in the meaning of terms relevant to this study – linear, stationary and multiplicative – and functions that can transform observed fluxes to a normal distribution as well or better than the logarithm.

Jeffrey D Scargle↗
DOE OSTI2025

Design of Hopfield Networks Based on Superconducting Coupled Oscillators

The global energy shortage has driven the development of many energy-efficient computational platforms beyond Moore's law, among which brain-inspired neuromorphic computing is one of the promising solutions. Associative memory and pattern recognition are important computations solved by brain-inspired Hopfield networks. Classical Hopfield networks store memories via fixed point attractors of their dynamics. In oscillatory Hopfield networks, these attractors are replaced by periodic orbits. Here, we design an oscillatory Hopfield network based on coupled superconducting oscillators. We first employ a mathematical phase reduction approach to map networks of coupled superconducting rapid single flux quantum (RSFQ) ring oscillators to coupled Kuramoto phase-oscillator networks. We use this theory to numerically optimize the hardware's mutual inductances in order to directly match the phase-reduced superconducting oscillators to a model of phase-oscillator-based Hopfield networks. The resulting network can store multiple oscillatory phase-locked memory patterns and recover the patterns based on the initial phase conditions. As different pattern recognition tasks, or learning, require tunable connectivity strengths between the oscillatory nodes, we further employ a coupler circuit that enables tuning the coupling strength between two oscillators by applying an external flux. We demonstrate the functionality of our design through numerical simulations of a small example network with oscillators operating at 86 GHz and recognizing patterns within 10 ns. Our approach enables the learning and retrieval of dynamical memory patterns with a wide range of applications where rhythmic dynamic output is beneficial.

Cheng, Ran↗
NASA NTRS1996

Two-Flux Green's Function Analysis for Transient Spectral Radiation in a Composite

An analysis is developed for obtaining transient temperatures in a two-layer semitransparent composite with spectrally dependent properties. Each external boundary of the composite is subjected to radiation and convection. The two-flux radiative transfer equations are solved by deriving a Green's function. This yields the local radiative heat source needed to numerically solve the transient energy equation. An advantage of the two-flux method is that isotropic scattering is included without added complexity. The layer refractive indices are larger than one. This produces internal reflections at the boundaries and the internal interface; the reflections are assumed diffuse. Spectral results using the Green's function method are verified by comparing with numerical solutions using the exact radiative transfer equations. Transient temperature distributions are given to illustrate the effect of radiative heating on one side of a composite with external convective cooling. The protection of a material from incident radiation is illustrated by adding scattering to the layer adjacent to the radiative source.

Siegel, Robert↗
NASA NTRS2013

Orbital Debris Engineering Model (ORDEM) v.3

A model of the manmade orbital debris environment is required by spacecraft designers, mission planners, and others in order to understand and mitigate the effects of the environment on their spacecraft or systems. A manmade environment is dynamic, and can be altered significantly by intent (e.g., the Chinese anti-satellite weapon test of January 2007) or accident (e.g., the collision of Iridium 33 and Cosmos 2251 spacecraft in February 2009). Engineering models are used to portray the manmade debris environment in Earth orbit. The availability of new sensor and in situ data, the re-analysis of older data, and the development of new analytical and statistical techniques has enabled the construction of this more comprehensive and sophisticated model. The primary output of this model is the flux [#debris/area/time] as a function of debris size and year. ORDEM may be operated in spacecraft mode or telescope mode. In the former case, an analyst defines an orbit for a spacecraft and "flies" the spacecraft through the orbital debris environment. In the latter case, an analyst defines a ground-based sensor (telescope or radar) in terms of latitude, azimuth, and elevation, and the model provides the number of orbital debris traversing the sensor's field of view. An upgraded graphical user interface (GUI) is integrated with the software. This upgraded GUI uses project-oriented organization and provides the user with graphical representations of numerous output data products. These range from the conventional flux as a function of debris size for chosen analysis orbits (or views), for example, to the more complex color-contoured two-dimensional (2D) directional flux diagrams in local spacecraft elevation and azimuth.

Matney, Mark↗
DOE OSTI2025

Probabilistic flux limiters

The stable numerical integration of shocks in compressible flow simulations relies on the reduction or elimination of Gibbs phenomena (unstable, spurious oscillations). A popular method to virtually eliminate Gibbs oscillations caused by numerical discretization in under-resolved simulations is to use a flux limiter. A wide range of flux limiters have been studied in the literature, with recent interest in their optimization via machine learning methods trained on high-resolution datasets. The common use of flux limiters in numerical codes as plug-and-play blackbox components makes them key targets for design improvement. Even for deterministic dynamical models, numerical uncertainty is introduced via coarse-graining required by insufficient computational power to solve all scales of motion. Conventional flux limiters are deterministic and lack the capacity to address uncertainties, both aleatoric (inherent randomness) and epistemic (modeling uncertainty due to limited knowledge), which arise in coarse-grained numerical simulations. Here, we introduce a conceptually distinct type of flux limiter that is designed to handle the effects of randomness in the model and uncertainty in model parameters. Unlike traditional single-function flux limiters, these new probabilistic flux limiters incorporate multiple flux limiting functions, each applied with a learned probability drawn from high-resolution data to mitigate the effects of uncertainty in numerical simulations. This approach departs from traditional single-function limiters by explicitly modeling and incorporating uncertainty into the shock capturing process. Using the example of Burgers' equation as a testbed, we show that a machine learned, probabilistic flux limiter may be used in a shock capturing code to more accurately capture shock profiles. In particular, we show that our probabilistic flux limiter outperforms standard limiters and can be successively improved upon (up to a point) by expanding the set of probabilistically chosen flux limiting functions.

97 MATHEMATICS AND COMPUTING↗