Search NASA⌕ Search

SEARCH · Search NASA

Results for “Numerical methods and modeling”

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 127 records · Page 7
All sourcesNASA NTRSDOE OSTI
NASA NTRSDate not supplied

Modeling and Impact of Solar Energetic Particles in the Heliosphere and Geospace

Understanding the radiation environment due to solar energetic particles in the heliosphere and the Earth’s magnetosphere is a challenging and practically important task. Exposure to energetic particles often leads to malfunctions and unexpected failures of electronics onboard spacecraft. The most vulnerable are exploratory missions when outside of the Earth’s magnetosphere. Geomagnetic field deflects Solar Energetic Particles (SEPs) moving through geospace though some of these particles propagate to LEO and have a high penetrating capability, thus producing significant radiation hazards for human spaceflight. Solar energetic particles also have an essential effect on the composition and dynamics of the Earth’s atmosphere. Precipitating SEPs enhance the atmospheric concentration of NOx and HOx, which play a crucial role in the ozone balance in the middle atmosphere by destroying odd oxygen through catalytic reactions. Numerical modeling of the radiation environment due to SEPs in the inner heliosphere and geospace is a multifold problem. That includes simulating 1) solar wind dynamics and the interplanetary magnetic field, 2) global modeling of the Earth’s magnetosphere, and 3) modeling transport and acceleration of SEPs in the inner heliosphere and geospace. The lecture will outline the key physical processes that control the behavior of SEPs in the heliosphere and geospace, along with the contemporary numerical methods used for their modeling. It will primarily focus on describing the SEP population in geospace across different altitudes, ranging from Low Earth Orbit (LEO) through Medium Earth Orbit (MEO) and Geostationary Orbit (GEO), up to the magnetopause, while considering the realistic geomagnetic field. Additionally, the lecture will cover how the SEP population in geospace varies in response to geomagnetic activity. This includes the temporal trapping of SEPs in geospace and the reduction of the rigidity cutoff during geomagnetic storms.

solar energetic particles↗
NASA NTRS1982

Interactive initialization of heat flux parameters for numerical models using satellite temperature measurements

A method for obtaining patterns of moisture availability (and net evaporation) from satellite infrared measurements employs Carlson's boundary layer model and a variety of image processing routines executed by a minicomputer. To test the method with regard to regional scale moisture analyses, two case studies were chosen because of the availability of HCMM data and because of the presence of a large horizontal gradient in antecedent precipitation and crp moisture index. Results show some correlation in both cases between antecedent precipitation and derived moisture availability. Apparently, regional-scale moisture availability patterns can be determined with some degree of fidelity but the values themselves may be useful only in the relative sense and significant to within plus or minus one category of dryness over a range of 4 or 5 categories between absolutely dry and field saturation. Preliminary results suggest that the derived moisture values correlate best with longer-term precipitation totals, suggesting that the infrared temperatures respond more sensitively to a relatively deep substrate layer.

Carlson, T. N.↗
NASA NTRS2007

Semianalytical Propagation of Satellite Orbits about an Arbitrary Central Body

Precision mean element (PME) satellite theories play a key role in orbit dynamics analyses. These theories employ: nonsingular orbital elements comprehensive force models Generalized Method of Averaging Numerical interpolation concepts The Draper Semianalytical Satellite Theory (DSST) (Refs. 1 - 6), whose development was led by the author, and the independently-developed Universal Semianalytical Method (USM) (Ref. 7) are examples of such theories. These theories provide the capability to tailor the force modeling to meet the desired computational speed vs. accuracy trade-off. The flexibility of such theories is demonstrated by their ability to include complicated atmosphere density models and spacecraft models in the perturbation theory context. The value of high speed satellite theories, in this era of computational plenty, is that they allow new ways of looking at astrodynamical problems such as orbit design (Refs. 8, 9) and atmosphere density updating (Refs. 10, 11). In the mid to late-1980 s, the geodynamics community led the development of very precise geopotential models such as GEM T2 and GEM T3 (Ref. 12), and with the subsequent analysis of the TOPEX flight data, JGM-2 and JGM-3 (Ref. 13). These were high degree and order geopotentials, at least 50 x 50. In 1993, the DSST implementation in the GTDS program was extended to include the 50 x 50 geopotential models (Ref. 14). The 50 x 50 geopotential, J2000 integration coordinate system, and solid Earth tide capabilities were integrated in GTDS by Scott Carter (Ref. 15). This capability demonstrated 1 m accuracy versus the TOPEX Precise Orbit Ephemerides. Subsequently the DSST Standalone program was also extended to include high degree and order geopotential models (Ref. 5). More recently GTDS has been hosted in the Linux PC environment. However, all of these efforts have been limited to modeling the motion of an artificial Earth satellite. They did not consider the additional complexities associated with lunar, planetary, or other natural satellite orbiters. Such complexities include: additional coordinate systems (associated with the direction of the north pole of rotation and the prime meridian of the new central bodies) (Ref. 16) normalized gravity model coefficients (desirable for high degree and order fields) (Ref. 17) indirect oblateness

Cefola, Paul J.↗
NASA NTRS1990

Bipolar supernova remnants and the obliquity dependence of shock acceleration

The diffusive shock acceleration mechanism proposed to explain the bipolarity observed in the synchrotron radio emission of young adiabatically expanding shell SNRs is investigated by means of numerical simulations. The theoretical basis of the SNR models and the numerical computation methods are explained, and the results are presented in graphs and synthetic radio maps and discussed in detail. It is found that the efficiency of the acceleration process depends on the obliquity angle theta(Bn) between the shock normal and the uniform magnetic field: models with theta(Bn) of about 90 deg can reproduce the observed azimuthal intensity ratios in most cases, but models with theta(Bn) near 0 deg cannot.

Fulbright, Michael S.↗
NASA NTRS1985

Spectral methods for modeling supersonic chemically reacting flow fields

A numerical algorithm was developed for solving the equations describing chemically reacting supersonic flows. The algorithm employs a two-stage Runge-Kutta method for integrating the equations in time and a Chebyshev spectral method for integrating the equations in space. The accuracy and efficiency of the technique were assessed by comparison with an existing implicit finite-difference procedure for modeling chemically reacting flows. The comparison showed that the procedure presented yields equivalent accuracy on much coarser grids as compared to the finite-difference procedure with resultant significant gains in computational efficiency.

Drummond, J. P.↗
NASA NTRS1984

A method for modeling finite-core vortices in wake-flow calculations

A numerical method for computing nonplanar vortex wakes represented by finite-core vortices is presented. The approach solves for the velocity on an Eulerian grid, using standard finite-difference techniques; the vortex wake is tracked by Lagrangian methods. In this method, the distribution of continuous vorticity in the wake is replaced by a group of discrete vortices. An axially symmetric distribution of vorticity about the center of each discrete vortex is used to represent the finite-core model. Two distributions of vorticity, or core models, are investigated: a finite distribution of vorticity represented by a third-order polynomial, and a continuous distribution of vorticity throughout the wake. The method provides for a vortex-core model that is insensitive to the mesh spacing. Results for a simplified case are presented. Computed results for the roll-up of a vortex wake generated by wings with different spanwise load distributions are presented; contour plots of the flow-field velocities are included; and comparisons are made of the computed flow-field velocities with experimentally measured velocities.

Stremel, P. M.↗
NASA NTRS1992

Recovery of the Navier-Stokes equations using a lattice-gas Boltzmann method

A lattice Boltzmann model is presented which gives the complete Navier-Stokes equation and may provide an efficient parallel numerical method for solving various fluid problems. The model uses the single-time relaxation approximation and a particular Maxwell-type distribution. The model eliminates exactly (1) the non-Galilean invariance caused by a density-dependent coefficient in the convection term and (2) a velocity-dependent equation of state.

Chen, Hudong↗
NASA NTRS1981

Advantages of multigrid methods for certifying the accuracy of PDE modeling

Numerical techniques for assessing and certifying the accuracy of the modeling of partial differential equations (PDE) to the user's specifications are analyzed. Examples of the certification process with conventional techniques are summarized for the three dimensional steady state full potential and the two dimensional steady Navier-Stokes equations using fixed grid methods (FG). The advantages of the Full Approximation Storage (FAS) scheme of the multigrid technique of A. Brandt compared with the conventional certification process of modeling PDE are illustrated in one dimension with the transformed potential equation. Inferences are drawn for how MG will improve the certification process of the numerical modeling of two and three dimensional PDE systems. Elements of the error assessment process that are common to FG and MG are analyzed.

Forester, C. K.↗
NASA NTRS1980

Mn-53 and Al-26 evidence for solar cosmic ray constancy - An improved model for interpretation

A general, numerical method for calculating activities of solar cosmic ray produced radionuclides at any point within an irregularly shaped lunar rock of known surface contour and lunar surface orientation is described. This method is then used to predict the activities of Mn-53 andAl-26 as a function of postion within lunar rock 68815 for various assumed values of solar cosmic ray flux (J), rigidity (R sub 0), and rock erosion rate (ER). The predicted activities agree with the measured activities of Kohl et al. (1978) when values of R sub 0 = 100 MV, J = 70 p/sq cm-sec (4 pi, E greater than 10 MeV), ER not greater than 1 mm/m.y. and a total exposure time of 2 m.y. are assumed. These values are in agreement with those found for rocks exposed for not less than 10 m.y. and provide no evidence for variation of the average solar cosmic ray parameters between the last 2 and 10 m.y. intervals. When interpreted with the improved model the activity vs. depth profiles for three faces of 68815 show no evidence of SCR anisotropy or differential erosion.

Russ, G. P., III↗
NASA NTRS2017

Integrated Model of the Eye/Optic Nerve Head Biomechanical Environment

Visual Impairment and Intracranial Pressure (VIIP) syndrome is a concern for long-duration space flight. Previously, it has been suggested that ocular changes observed in VIIP syndrome are related to the cephalad fluid shift that results in altered fluid pressures [1]. We are investigating the impact of changes in intracranial pressure (ICP) using a combination of numerical models, which simulate the effects of various environment conditions, including finite element (FE) models of the posterior eye. The specific interest is to understand how altered pressures due to gravitational changes affect the biomechanical environment of tissues of the posterior eye and optic nerve sheath. METHODS: Additional description of the numerical modeling is provided in the IWS abstract by Nelson et al. In brief, to simulate the effects of a cephalad fluid shift on the cardiovascular and ocular systems, we utilized a lumped-parameter compartment model of these systems. The outputs of this lumped-parameter model then inform boundary conditions (pressures) for a finite element model of the optic nerve head (Figure 1). As an example, we show here a simulation of postural change from supine to 15 degree head-down tilt (HDT), with primary outcomes being the predicted change in strains at the optic nerve head (ONH) region, specifically in the lamina cribrosa (LC), retrolaminar optic nerve, and prelaminar neural tissue (PLNT). The strain field can be decomposed into three orthogonal components, denoted as the first, second and third principal strains. We compare the peak tensile (first principal) and compressive (third principal) strains, since elevated strain alters cell phenotype and induces tissue remodeling. RESULTS AND CONCLUSIONS: Our lumped-parameter model predicted an IOP increase of c. 7 mmHg after 21 minutes of 15 degree HDT, which agreed with previous reports of IOP in HDT [1]. The corresponding FEM simulations predicted a relative increase in the magnitudes of the peak tensile and compressive strains in the lamina cribrosa of 42 and 43, respectively (Fig. 2). The corresponding changes in the optic nerve strains were 17 and 39, while in the PLNT they were 47 and 43. These magnitudes of relative elevations in peak strains may induce a phenotypic response in resident mechano-responsive resident cells [2]. This approach may be expanded to investigate other environmental changes (e.g. parabolic flight). Through our VIIP SCHOLAR project, we will validate and improve these integrated models by measuring patient-specific changes in optic nerve sheath geometry in patients with idiopathic intracranial hypertension before and after lumbar puncture and CSF removal.

biomechanics↗
NASA NTRS1976

Turbulence and numerical wind tunnels

The feasibility and possible significance of aerodynamic simulations are considered. Despite certain limitations, there is reason to believe that many useful and important flows containing turbulence can be computed with an accuracy acceptable for use in aircraft designs. Attention is given to the equations for unsteady fluid flow, the direct calculation of unsteady viscous flows, the direct calculation of nonlinear equations, the computation of flows about bodies using turbulence modeling, future developments in numerical methods, and future developments in computers.

Lomax, H.↗
NASA NTRS1978

A scheme for computing surface layer turbulent fluxes from mean flow surface observations

A physical model and computational scheme are developed for generating turbulent surface stress, sensible heat flux and humidity flux from mean velocity, temperature and humidity at some fixed height in the atmospheric surface layer, where conditions at this reference level are presumed known from observations or the evolving state of a numerical atmospheric circulation model. The method is based on coupling the Monin-Obukov surface layer similarity profiles which include buoyant stability effects on mean velocity, temperature and humidity to a force-restore formulation for the evolution of surface soil temperature to yield the local values of shear stress, heat flux and surface temperature. A self-contained formulation is presented including parameterizations for solar and infrared radiant fluxes at the surface. Additional parameters needed to implement the scheme are the thermal heat capacity of the soil per unit surface area, surface aerodynamic roughness, latitude, solar declination, surface albedo, surface emissivity and atmospheric transmissivity to solar radiation.

Hoffert, M. I.↗
NASA NTRS1978

Status and future prospects of using numerical methods to study complex flows at High Reynolds numbers

The calculation of flow fields past aircraft configuration at flight Reynolds numbers is considered. Progress in devising accurate and efficient numerical methods, in understanding and modeling the physics of turbulence, and in developing reliable and powerful computer hardware is discussed. Emphasis is placed on efficient solutions to the Navier-Stokes equations.

Maccormack, R. W.↗
NASA NTRS1981

Studies in astronomical time series analysis. I - Modeling random processes in the time domain

Several random process models in the time domain are defined and discussed. Attention is given to the moving average model, the autoregressive model, and relationships between and combinations of these models. Consideration is then given to methods for investigating pulse structure, procedures of model construction, computational methods, and numerical experiments. A FORTRAN algorithm of time series analysis has been developed which is relatively stable numerically. Results of test cases are given to study the effect of adding noise and of different distributions for the pulse amplitudes. A preliminary analysis of the light curve of the quasar 3C 272 is considered as an example.

Scargle, J. D.↗
NASA NTRS1984

Predictive thermal-vacuum modeling by variable coefficients

In computer-assisted thermal-vacuum testing, an accurate mathematical thermal model is necessary. This model can indicate what effect a given change in control temperature will produce on the test article. A three-node thermal model is used to illustrate some numerical methods for increasing the accuracy by generating more accurate coefficients for the model.

Chuvala, J. T.↗
NASA NTRS1984

Direct numerical simulation of curved turbulent channel flow

Low Reynolds number, mildly curved, turbulent channel flow has been simulated numerically without subgrid scale models. A new spectral numerical method developed for this problem was used, and the computations were performed with 2 million degrees of freedom. A variety of statistical and structural information has been extracted from the computed flow fields. These include mean velocity, turbulence stresses, velocity skewness, and flatness factors, space time correlations and spectra, all the terms in the Reynolds stress balance equations, and contour and vector plots of instantaneous velocity fields. The effects of curvature on this flow were determined by comparing the concave and convex sides of the channel. The observed effects are consistent with experimental observations for mild curvature. The most significant difference in the turbulence statistics between the concave and convex sides was in the Reynolds shear stress. This was accompanied by significant differences in the terms of the Reynolds shear stress balance equations. In addition, it was found that stationary Taylor-Gortler vortices were present and that they had a significant effect on the flow by contributing to the mean Reynolds shear stress, and by affecting the underlying turbulence.

Moser, R. D.↗
NASA NTRS1985

Particle behavior in the magnetosphere

The Rice Convection Model deals with large-scale processes in the earth's inner and middle magnetosphere, including coupling to the ionosphere. Starting from appropriate initial and boundary conditions, the model computes the following physical parameters: ionospheric electric fields and currents; magnetospheric particle distributions, electric fields, and electric currents; and magnetic-field-aligned (Birkeland) currents connecting the two regions. This paper evaluates work on the model, with emphasis on the assumptions made, the basic equations, and the numerical methods. The theoretical basis of the model is compared and contrasted with standard magnetohydrodynamics. The limitations imposed by the major assumptions are discussed. Model inputs and boundary conditions are listed, and the methods of specifying them discussed. Some physical conclusions and insights that have been gained from the model are listed and described very briefly. References are given to published discussions of the major points of physics.

Wolf, R. A.↗