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.

SourceResearch connectionFrom yearThrough year

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

At least 37 records · Page 2

All sources NASA NTRS DOE OSTI

NASA NTRS1982

Magnetopause modeling - Flux transfer events and magnetosheath quasi-trapped distributions

Three-dimensional distribution functions for energetic ions are studied numerically in the magnetosphere, through the magnetopause, and in the magnetosheath using a simple one-dimensional quasi-static model and ISEE 1 magnetopause crossing data for November 10, 1977. Quasi-trapped populations in the magnetosheath observed near flux transfer events (FTEs) are investigated, and it is shown that the population in the sheath appears to sandwich the FTE distributions. These quasi-trapped distributions are due to slow, large pitch angle, outward moving particles left behind by the outward rush of the ions more field-aligned at the time the flux was opened. It is found that sheath convective flows can map along the connected flux tube without drastically changing the distribution function, and results suggest that localized tangential fields above the upper limit may exist.

Speiser, T. W.↗

NASA NTRS1986

The equilibrium structure of thin magnetic flux tubes. II

The thermal structure of the medium inside thin, vertical magnetic flux tubes embedded in a given external atmosphere is investigated, assuming cylindrical symmetry and a depth-independent plasma beta. The variation with tube radius of the temperature on the tube axis is computed and the temperature on the tube wall is estimated. The temperature variation across the flux tube is found to be due to the depth variation of the intensity and to the density stratification of the atmosphere. Since the temperature difference between the axis and the wall is small in thin flux tubes (of the order of 10 percent), the horizontal temperature gradient may often be neglected and the temperature in a tube of given radius may be described by a single function of depth. Thus, a more detailed numerical treatment of the radiative transfer within thin flux tubes can be substantially simplified by neglecting horizontal temperature differences within the flux tube proper.

Kalkofen, W.↗

NASA NTRS1985

Precipitating electron interaction with the atmosphere. II - The dayside cusp region

Interaction of precipitating low-energy magnetosheath electrons with the atmosphere in the dayside cusp region has been studied. Both pitch angle and energy distributions of the fluxes as well as excitation functions for selected N2 and O UV emissions were obtained by numerically solving the multiangle equations of electron transport. There is some possibility that atmospheric emissions may be used for remote measurements of incident soft energy flux, because the ratios of molecular to atomic emission line intensities in the low-energy region are quite different from those in the high-energy region.

Prasad, S. S.↗

DOE OSTI2026

Neural entropy-stable conservative flux form neural networks for learning hyperbolic conservation laws

We propose a neural entropy-stable conservative flux form neural network (NESCFN) for learning hyperbolic conservation laws and their associated entropy functions directly from solution trajectories, without requiring any predefined numerical discretization. While recent neural network architectures have successfully integrated classical numerical principles into learned models, most rely on prior knowledge of the governing equations or assume a fixed discretization. Our approach removes this dependency by embedding entropy-stable design principles into the learning process itself, enabling the discovery of physically consistent dynamics in a fully data-driven setting. By jointly learning both the flux function and a corresponding entropy, NESCFN promotes conservation and entropy dissipation, which is critical for long-term stability and fidelity in the system of hyperbolic conservation laws. Furthermore, numerical results demonstrate that the method achieves stability and conservation over extended time horizons and accurately captures shock propagation speeds, even without oracle access to future-time solution profiles in the training data.

Conservative flux form↗

NASA NTRS1996

Toward a Turbulence Constitutive Relation for Rotating Flows

In rapidly rotating turbulent flows the largest scales of the motion are in approximate geostrophic balance. Single-point turbulence closures, in general, cannot attain a geostrophic balance. This article addresses and resolves the possibility of constitutive relation procedures for single-point second order closures for a specific class of rotating or stratified flows. Physical situations in which the geostrophic balance is attained are described. Closely related issues of frame-indifference, horizontal nondivergence, Taylor-Proudman theorem and two-dimensionality are, in the context of both the instantaneous and averaged equations, discussed. It is shown, in the absence of vortex stretching along the axis of rotation, that turbulence is frame-indifferent. A derivation and discussion of a geostrophic constraint which the prognostic equations for second-order statistics must satisfy for turbulence approaching a frame-indifferent limit is given. These flow situations, which include rotating and nonrotating stratified flows, are slowly evolving flows in which the constitutive relation procedures are useful. A nonlinear non-constant coefficient representation for the rapid-pressure strain covariance appearing in the Reynolds stress and heat flux equations consistent with the geostrophic balance is described. The rapid-pressure strain model coefficients are not constants determined by numerical optimization but are functions of the state of the turbulence as parameterized by the Reynolds stresses and the turbulent heat fluxes. The functions are valid for all states of the turbulence attaining their limiting values only when a limit state is achieved. These issues are relevant to strongly vortical flows as well as flows such as the planetary boundary layers, in which there is a transition from a three-dimensional shear driven turbulence to a geostrophic or horizontal turbulence.

Ristorcelli, J. R.↗

NASA NTRS1989

Heat flux in a non-Maxwellian plasma

A hybrid numerical scheme is applied to solve the Landau equation for the electron distribution function over all velocity space. Evidence is presented for the first time of the degree and character of the failure of the classical Spitzer-Haerm heat flux approximation in a realistic solar coronal loop structure. In the loop model used, the failure is so severe at some points that the role of the heat flux in the plasma's energy balance is completely misinterpreted. In the lower corona the Spitzer-Haerm approximation predicts that the heat flux should act as an energy source, whereas the more accurate distribution functions calculated here show this to be an energy sink.

Ljepojevic, N. N.↗

NASA NTRSDate not supplied

CFD Predictions of Boiling Regime Transitions during Line Chilldown validated against a 1G LN2 Experiment

Introduction Before filling a propellant tank on the ground or in Space, the transfer line between the donor and receiver tanks must be cooled down preferably by sacrificing a minimum amount of the cryogenic fluid. The cryogenic line chill-down process involves a transition between different flow boiling regimes, namely, film boiling, transition film boiling, and nucleate boiling which are complex and may be quite gravity-dependent. Capturing these boiling phenomena and predicting the transition between them in a CFD framework is new and challenging both for 1g and microgravity applications. Materials & Methods The present work addresses this challenge by employing a two-phase Eulerian approach in the context of a homogeneous fluid mixture together with the Lee phase change model to capture the film boiling regime of the chill-down process using ANSYS Fluent®. The nucleate boiling regime is predicted by incorporating an in-house developed sub-grid model that accounts for bubble nucleation, bubble growth, bubble departure diameter, and their shedding frequency. The sub-grid model is encoded and implemented into Fluent via a user-defined function for the wall-fluid heat flux calculations. The mathematical formulation and numerical implementation of the CFD model are described in detail. The coupled CFD-Subgrid model is validated against published experimental data for liquid nitrogen chill-down of a heated stainless-steel pipe in 1g. Results Numerical simulation results show good agreements between the CFD predictions of the wall temperature evolution, rewetting temperature, and transition between film and nucleate boiling, with the experimental measurements published by Darr et al [2] for several different LN2 flowrates in the vertical pipe orientation. The CFD predictions for the wall temperature distribution indicate a rapid quenching of the wall at two upstream and downstream temperature sensing locations as compared to the experimental measurement. The only tuning parameter in the CFD model is the Lee mass transfer coefficient. The CFD Model predicts the Liedenfrost rewetting temperature in close agreement with the experiment. This marks a transition between stable and transitionary flow boiling regimes. The CFD-predicted boiling curve for the downstream sensor location is also compared against its experimental counterpart and indicates that the model is able to predict all the key temperature and heat flux parameters during the transitions from stable to transitionary film boiling to nucleate boiling in close agreement with the experiment. A sequence of predicted volume fraction, and temperature contours depicting these transitions will be presented.

Evaporation Condensation↗

NASA NTRS1995

Eight-moment approximation solar wind models

Heat conduction from the corona is important in the solar wind energy budget. Until now all hydrodynamic solar wind models have been using the collisionally dominated gas approximation for the heat conductive flux. Observations of the solar wind show particle distribution functions which deviate significantly from a Maxwellian, and it is clear that the solar wind plasma is far from collisionally dominated. We have developed a numerical model for the solar wind which solves the full equation for the heat conductive flux together with the conservation equations for mass, momentum, and energy. The equations are obtained by taking moments of the Boltzmann equation, using an 8-moment approximation for the distribution function. For low-density solar winds the 8-moment approximation models give results which differ significantly from the results obtained in models assuming the gas to be collisionally dominated. The two models give more or less the same results in high density solar winds.

Olsen, Espen Lyngdal↗

NASA NTRS1994

Computing unsteady shock waves for aeroacoustic applications

The computation of unsteady shock waves, which contribute significantly to noise generation in supersonic jet flows, is investigated. The paper focuses on the difficulties of computing slowly moving shock waves. Numerical error is found to manifest itself principally as a spurious entropy wave. Calculations presented are performed using a third-order essentially nonoscillatory scheme. The effect of stencil biasing parameters and of two versions of numerical flux formulas on the magnitude of spurious entropy are investigated. The level of numerical error introduced in the calculation is quantified as a function of shock pressure ratio, shock speed, Courant number, and mesh density. The spurious entropy relative to the entropy jump across a static shock decreases with increasing shock strength and shock velocity relative to the grid, but is insensitive to Courant number. The structure of the spurious entropy wave is affected by the choice of flux formulas and algorithm biasing parameters. The effect of the spurious numerical waves on the calculation of sound amplification by a shock wave is investigated. For this class of problem, the acoustic pressure waves are relatively unaffected by the spurious numerical phenomena.

Meadows, Kristine R.↗

NASA NTRS1994

Computing unsteady shock waves for aeroacoustic applications

The computation of unsteady shock waves, which contribute significantly to noise generation in supersonic jet flows, is investigated. This paper focuses on the difficulties of computing slowly moving shock waves. Numerical error is found to manifest itself principally as a spurious entropy wave. Calculations presented are performed using a third order essentially nonoscillatory scheme. The effect of stencil biasing parameters and of two versions of numerical flux formulas on the magnitude of spurious entropy are investigated. The level of numerical error introduced in the calculation in quantified as a function of shock pressure ratio, shock speed, Courant number, and mesh density. The spurious entropy relative to the entropy jump across a static shock decreases with increasing shock strength and shock velocity relative to the grid, but is insensitive to Courant number. The structure of the spurious entropy wave is affected by the choice of flux formulas and algorithm biasing parameters. The effect of the spurious numerical waves on the calculation of sound amplification by a shock wave is investigated. For this class of problem, the acoustic pressure waves are relatively unaffected by the spurious numerical phenomena.

Meadows,, Kristine r.↗

NASA NTRS1991

Studies of magnetopause structure

From the International Sun-Earth Explorers (ISEE) 1 magnetopause crossings on November 10, 1977, three-dimensional distribution functions for energetic ions were studied in the magnetosphere, through the magnetopause, and in the magnetosheath (Speiser and Williams, 1982). The particle distributions were particularly examined at and near the times that Russell and Elphic (1978) identified as flux transfer events (FTE). Using a simple, one-dimensional, quasi-static model, particle orbits were followed numerically, from the magnetosphere into the sheath. The inner, trapped, distribution initializes the distribution function. Liouville's theorem allows the inner distribution to be mapped into the sheath following the orbits. This mapping is shown for four mangetosheath ion flows (MIF's) corresponding to four flux transfer events. Results from the studies are discussed. A brief discussion of current sheet particle motion is presented.

Speiser, Theodore W.↗

NASA NTRS1995

Radiative transfer in a semiinfinite medium with a specularly reflecting boundary

A consistent solution of the radiative transfer equation characterizing photon transport in a semi-infinite medium of refractive index greater than or equal to one is obtained following the method of Sobolev. Fresnel specular reflection, Snell's law and isotropic scattering are assumed. An algorithm is developed and its accuracy is demonstrated. A numerical Laplace transform inversion leads to an efficient evaluation for the interior flux and source function distributions.

Ganapol, B. D.↗

NASA NTRS1979

Fermi liquid viscosity in a finite geometry

Forced flow of a Fermi liquid is studied for a cell geometry consisting of two planes with a separation on the order of mean free path. An approximate transport equation is used to derive an integral equation for the velocity profile, which is solved numerically. Results for the total flux through the cell, which determines the dissipation, are given as a function of the Knudsen number N (ratio of cell thickness to mean free path). Effects of specular reflection at the boundaries are considered. It is found that the dissipation has a minimum at N approximately equal to 1/2, and behaves linearly for N greater than or equal to 3. Implications for present experimentation are discussed.

Jaffe, J. E.↗

NASA NTRS1979

Fermi liquid viscosity in a finite geometry

Forced flow of a Fermi liquid is studied for a cell geometry consisting of two planes with a separation on the order of the mean free path. An approximate transport equation is used to derive an integral equation for the velocity profile, which is solved numerically. Results for the total flux through the cell, which determines the dissipation, are given as a function of the Knudsen number kappa (ratio of cell thickness to mean free path). Effects of specular reflection at the boundaries are considered. It is found that the dissipation has a minimum at kappa of 1/2, and behaves linearly for kappa not less than 3. Implications for present experimentation are discussed.

Jaffe, J. E.↗

NASA NTRS1979

Numerical solution of transonic full stream function equations in conservation form

The stream function equation in conservation form is solved iteratively based on the artificial compressibility method. The density is not a unique function of the mass flux. In order to avoid the ambiguity near the sonic line, the density is updated in terms of the velocity, which is obtained through a simple integration of a first order equation step by step in the flow field. Iteration algorithms and finite difference approximations are discussed and numerical results of both conservative and nonconservative calculations are presented.

Hafez, M. M.↗

DOE OSTI2025

Including the vacuum energy in stellarator coil design

Being three-dimensional, stellarators have the advantage that plasma currents are not essential for creating rotational-transform; however, the external current-carrying coils in stellarators can have strong geometrical shaping, which can complicate the construction. Reducing the inter-coil electromagnetic forces acting on strongly shaped 3D coils and the stress on the support structure while preserving the favorable properties of the magnetic field is a design challenge. In this work, we recognize that the inter-coil forces are the gradient of the vacuum magnetic energy. We introduce an objective functional built on the usual quadratic flux on a prescribed target surface together with a weighed penalty on the vacuum energy. The Euler–Lagrange equation for stationary states is derived, and numerical illustrations are computed using a modern stellarator optimization framework. A study of the effect of the energy functional on the inter-coil forces is conducted and the energy is shown to be a promising quantity in producing coils with low forces.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

NASA NTRS1993

The distribution of flare parameters and implications for coronal heating

Assuming that solar hard X-ray flares are caused by bremsstrahlung from a distribution of accelerated electrons impinging on a thick target, we determine the relationship between the total flare energy released in the corona and observable flare parameters. It is shown that for the more numerous smaller flares to be energetically more important than larger ones, there must be an anticorrelation between observed flux and one or more of the photon spectral index, the total flare duration, the ratio of the total energy that goes into heating the corona to the energy that goes into acceleration of electrons, and the reciprocal of the low-energy cutoff of the accelerated electron spectrum. Using a new statistical test which properly takes account of data truncation due to observational selection effects, we find that neither spectral index nor duration show an anticorrelation with flux in ISEE 3/ICE hard X-ray data. Using a nonparametric method, we obtain single distributions of peak flux, spectral index, duration, and fluence, corrected for selection biases. The functional forms of these distributions do not support coronal heating by numerous small flares without very specific behavior by yet unobserved parameters.

Lee, Theodore T.↗

NASA NTRS1973

An improved separability approximation for line radiative transport in nonhomogeneous media.

A simple modification to the constant half-width approximation of Wilson and Greif, which is an extension of the well-known Curtis-Godson method to the treatment of temperature variations along the integration path, is introduced which permits a more accurate evaluation of line radiative transport in nonhomogeneous gases. To demonstrate the method's accuracy, comparisons are made with Wilson and Greif and numerical frequency integrated results for the line equivalent width and radiative flux in a planar slab with prescribed Planck function and line half-width spatial variations. These variations are chosen to represent typical shock layer conditions. It is found that the modified procedure reduces the inaccuracies inherent in Wilson and Greif's approximation by factors ranging from 5 to 10, while retaining the latter method's ease of application.

Boughner, R. E.↗