Search NASA⌕ Search

SEARCH · Search NASA

Results for “Physics 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 451 records · Page 25

Plant canopy transpiration in bioregenerative life support systems - The link between mechanistic and empirical models

A model of water transpiration in a plant canopy that combines two approaches is presented. The first approach is to account for underlying physical processes, while the second is to empirically incorporate transpiration data now being generated at the Johnson Center Variable Pressure Growth Chamber. The two approaches, physical modeling and data analysis, make it possible to produce a model that is more robust than either the standard first-principles model or a straightforward empirical model. It is shown that the present transpiration model is able to efficiently capture the dynamic behavior of the plant canopy over the entire range of environmental parameters now envisioned to be important in an operating controlled ecological life support system (CELSS). Examples of the use of this model in assessing plant canopy dynamics and CELSS design options are also presented.

Sirko, Robert J.↗

Code-to-Code Comparison, and Material Response Modeling of Stardust and MSL using PATO and FIAT

This report provides a code-to-code comparison between PATO, a recently developed high fidelity material response code, and FIAT, NASA's legacy code for ablation response modeling. The goal is to demonstrates that FIAT and PATO generate the same results when using the same models. Test cases of increasing complexity are used, from both arc-jet testing and flight experiment. When using the exact same physical models, material properties and boundary conditions, the two codes give results that are within 2% of errors. The minor discrepancy is attributed to the inclusion of the gas phase heat capacity (cp) in the energy equation in PATO, and not in FIAT.

ablative materials↗

Acoustic space occupancy: Combining ecoacoustics and lidar to model biodiversity variation and detection bias across heterogeneous landscapes

There is global interest in quantifying changing biodiversity in human-modified landscapes. Ecoacoustics may offer a promising pathway for supporting multi-taxa monitoring, but its scalability has been hampered by the sonic complexity of biodiverse ecosystems and the imperfect detectability of animal-generated sounds. The acoustic signature of a habitat, or soundscape, contains information about multiple taxa and may circumvent species identification, but robust statistical technology for characterizing community-level attributes is lacking. Here, we present the Acoustic Space Occupancy Model, a flexible hierarchical framework designed to account for detection artifacts from acoustic surveys in order to model biologically relevant variation in acoustic space use among community assemblages. We illustrate its utility in a biologically and structurally diverse Amazon frontier forest landscape, a valuable test case for modeling biodiversity variation and acoustic attenuation from vegetation density. We use complementary airborne lidar data to capture aspects of 3D forest structure hypothesized to influence community composition and acoustic signal detection. Our novel analytic framework permitted us to model both the assembly and detectability of soundscapes using lidar-derived estimates of forest structure. Our empirical predictions were consistent with physical models of frequency-dependent attenuation, and we estimated that the probability of observing animal activity in the frequency channel most vulnerable to acoustic attenuation varied by over 60%, depending on vegetation density. There were also large differences in the biotic use of acoustic space predicted for intact and degraded forest habitats, with notable differences in the soundscape channels predominantly occupied by insects. This study advances the utility of ecoacoustics by providing a robust modeling framework for addressing detection bias from remote audio surveys while preserving the rich dimensionality of soundscape data, which may be critical for inferring biological patterns pertinent to multiple taxonomic groups in the tropics. Our methodology paves the way for greater integration of remotely sensed observations with high-throughput biodiversity data to help bring routine, multi-taxa monitoring to scale in dynamic and diverse landscapes.

Airborne lidar↗

Transfer function models for using empirical and physics-based simulation signal response data

In many situations, real or induced flaws such as tight cracks with known morphology cannot be manufactured in part geometry specimens or in real parts. Typically, surface fatigue cracks are manufactured in simple geometry specimens such as flat plates, dog-bone shaped flat or cylindrical specimens. If a nondestructive evaluation (NDE)technique is required to provide a reliably detectable flaw size, denoted asa90/95, for induced flaws in a part, then a direct method for qualifying the NDE procedure is to use appropriate induced flaw specimens and perform NDE procedure demonstration on the specimens. Probability of detection (POD)analysis of the empirical data may provide estimation of a90/95. This approach is described as direct POD demonstration testing, which may follow guidelines of MIL-HDBK-1823. This paper considers a case, where embedded tight cracklike induced flaws are to be detected reliably using a signal response based NDE procedure. Here, it is assumed that it is not practical to make surface or embedded induced flaw specimens in part geometry or configuration. Therefore, a direct POD demonstration testing cannot be undertaken. It is also assumed that simulation of signal response is possible for both surface and embedded induced flaws in part geometry specimens using a physics-based model. The proposed approach for NDE procedure qualification uses artificial flaws in simple geometry and part geometry specimens, and induced flaws in the same type of simple geometry specimens. Signal response data is taken on all sets of artificial and induced flaws in simple geometry and part geometry specimens. Moreover, simulated signal response data is generated for surface and embedded flaws. Thus, a case of five signal response versus flaw size datasets is considered. Three of the datasets are empirical and two datasets are physics model-based simulation datasets. A method of devising and using transfer function calculation dataset blocks to estimate the either the reliably detectable flaw size or the demonstration flaw size is provided.

Nondestructive evaluation↗

How Consistent are Recent Variations in the Tropical Energy and Water Cycle Resolved by Satellite Measurements?

One notable aspect of Earth's climate is that although the planet appears to be very close to radiative balance at top-of-atmosphere (TOA), the atmosphere itself and underlying surface are not. Profound exchanges of energy between the atmosphere and oceans, land and cryosphere occur over a range of time scales. Recent evidence from broadband satellite measurements suggests that even these TOA fluxes contain some detectable variations. Our ability to measure and reconstruct radiative fluxes at the surface and at the top of atmosphere is improving rapidly. One question is 'How consistent, physically, are these diverse remotely-sensed data sets'? The answer is of crucial importance to understanding climate processes, improving physical models, and improving remote sensing algorithms. In this work we will evaluate two recently released estimates of radiative fluxes, focusing primarily on surface estimates. The International Satellite Cloud Climatology Project 'FD' radiative flux profiles are available from mid-1983 to near present and have been constructed by driving the radiative transfer physics from the Goddard Institute for Space Studies (GISS) global model with ISCCP clouds and TOVS (TIROS Operational Vertical Sounder)thermodynamic profiles. Full and clear sky SW and LW fluxes are produced. A similar product from the NASA/GEWEX Surface Radiation Budget Project using different radiative flux codes and thermodynamics from the NASA/Goddard Earth Observing System (GEOS-1) assimilation model makes a similar calculation of surface fluxes. However this data set currently extends only through 1995. We also employ precipitation measurements from the Global Precipitation Climatology Project (GPCP) and the Tropical Rainfall Measuring Mission (TRMM). Finally, ocean evaporation estimates from the Special Sensor Microwave Imager (SSM/I) are considered as well as derived evaporation from the NCAR/NCEP Reanalysis. Additional information is included in the original extended abstract.

Robertson, F. R.↗

U.S. Eastern Continental Shelf Carbon Cycling (USECoS): Modeling, Data Assimilation, and Analysis

Although the oceans play a major role in the uptake of fossil fuel CO2 from the atmosphere, there is much debate about the contribution from continental shelves, since many key shelf fluxes are not yet well quantified: the exchange of carbon across the land-ocean and shelf-slope interfaces, air-sea exchange of CO2, burial, and biological processes including productivity. Our goal is to quantify these carbon fluxes along the eastern U.S. coast using models quantitatively verified by comparison to observations, and to establish a framework for predicting how these fluxes may be modified as a result of climate and land use change. Our research questions build on those addressed with previous NASA funding for the USECoS (U.S. Eastern Continental Shelf Carbon Cycling) project. We have developed a coupled biogeochemical ocean circulation model configured for this study region and have extensively evaluated this model with both in situ and remotely-sensed data. Results indicate that to further reduce uncertainties in the shelf component of the global carbon cycle, future efforts must be directed towards 1) increasing the resolution of the physical model via nesting and 2) making refinements to the biogeochemical model and quantitatively evaluating these via the assimilation of biogeochemical data (in situ and remotely-sensed). These model improvements are essential for better understanding and reducing estimates of uncertainties in current and future carbon transformations and cycling in continental shelf systems. Our approach and science questions are particularly germane to the carbon cycle science goals of the NASA Earth Science Research Program as well as the U.S. Climate Change Research Program and the North American Carbon Program. Our interdisciplinary research team consists of scientists who have expertise in the physics and biogeochemistry of the U.S. eastern continental shelf, remote-sensing data analysis and data assimilative numerical models.

Mannino, Antonio↗

Development of a three-dimensional time-dependent flow field model

A three-dimensional, time-dependent mathematical model to represent Mobile Bay was developed. Computer programs were developed which numerically solve the appropriate conservation equations for predicting bay and estuary flow fields. The model is useful for analyzing the dispersion of sea water into fresh water and the transport of sediment, and for relating field and physical model data.

Farmer, R. C.↗

Calculations of combustion response profiles and oscillations

The theory and procedures for determining the characteristics of pressure oscillations in rocket engines with prescribed burning rate oscillations are presented. Pressure and velocity oscillations calculated using this procedure are presented for the Space Shuttle Main Engine (SSME) to show the influence of baffles and absorbers on the burning rate oscillations required to achieve neutral stability. Results of calculations to determine local combustion responses using detailed physical models for injection, atomization, and vaporization with gas phase oscillations in baffled and unbaffled SSME combustors are presented. The contributions of the various physical phenomena occurring in a combustor to oscillations in combustion response were determined.

Priem, Richard J.↗

Computational Modeling of the Ames 11-Ft Transonic Wind Tunnel in Conjunction with IofNEWT

Technical advances in Computational Fluid Dynamics have now made it possible to simulate complex three-dimensional internal flows about models of various size placed in a Transonic Wind Tunnel. TWT wall interference effects have been a source of error in predicting flight data from actual wind tunnel measured data. An advantage of such internal CFD calculations is to directly compare numerical results with the actual tunnel data for code assessment and tunnel flow analysis. A CFD capability has recently been devised for flow analysis of the NASA/Ames 11-Ft TWT facility. The primary objectives of this work are to provide a CFD tool to study the NASA/Ames 11-Ft TWT flow characteristics, to understand the slotted wall interference effects, and to validate CFD codes. A secondary objective is to integrate the internal flowfield calculations with the Pressure Sensitive Paint data, a surface pressure distribution capability in Ames' production wind tunnels. The effort has been part of the Ames IofNEWT, Integration of Numerical and Experimental Wind Tunnels project, which is aimed at providing further analytical tools for industrial application. We used the NASA/Ames OVERFLOW code to solve the thin-layer Navier-Stokes equations. Viscosity effects near the model are captured by Baldwin-Lomax or Baldwin-Barth turbulence models. The solver was modified to model the flow behavior in the vicinity of the tunnel longitudinal slotted walls. A suitable porous type wall boundary condition was coded to account for the cross-flow through the test section. Viscous flow equations were solved in generalized coordinates with a three-factor implicit central difference scheme in conjunction with the Chimera grid procedure. The internal flow field about the model and the tunnel walls were descretized by the Chimera overset grid system. This approach allows the application of efficient grid generation codes about individual components of the configuration; separate minor grids were developed to resolve the model and overset onto a main grid which discretizes the interior of the tunnel test section. Individual grid components axe not required to have mesh boundaries joined in any special way to each other or to the main tunnel grid. Programs have been developed to rotate the model about the tunnel pivot point and rotation axis, similar to that of the tunnel turntable mechanism for adjusting the pitch of the physical model in the test section.

Djomehri, M. Jahed↗

Earth Satellite Population Instability: Underscoring the Need for Debris Mitigation

A recent study by NASA indicates that the implementation of international orbital debris mitigation measures alone will not prevent a significant increase in the artificial Earth satellite population, beginning in the second half of this century. Whereas the focus of the aerospace community for the past 25 years has been on the curtailment of the generation of long-lived orbital debris, active remediation of the current orbital debris population should now be reconsidered to help preserve near-Earth space for future generations. In particular, we show in this paper that even if launch operations were to cease today, the population of space debris would continue to grow. Further, proposed remediation techniques do not appear to offer a viable solution. We therefore recommend that, while the aerospace community maintains the current debris-limiting mission regulations and postmission disposal procedures, future emphasis should be placed on finding new remediation technologies for solving this growing problem. Since the launch of Sputnik 1, space activities have created an orbital debris environment that poses increasing impact risks to existing space systems, including human space flight and robotic missions (1, 2). Currently, more than 9,000 Earth orbiting man-made objects (including many breakup fragments), with a combined mass exceeding 5 million kilograms, are tracked by the US Space Surveillance Network and maintained in the US satellite catalog (3-5). Three accidental collisions between cataloged satellites during the period from late 1991 to early 2005 have already been documented (6), although fortunately none resulted in the creation of large, trackable debris clouds. Several studies conducted during 1991-2001 demonstrated, with assumed future launch rates, the unintended growth potential of the Earth satellite population, resulting from random, accidental collisions among resident space objects (7-13). In some low Earth orbit (LEO) altitude regimes where the number density of satellites is above a critical spatial density, the production rate of new satellites (i.e., debris) due to collisions exceeds the loss of objects due to orbital decay. NASA s evolutionary satellite population model LEGEND (LEO-to-GEO Environment Debris model), developed by the Orbital Debris Program Office at the NASA Lyndon B. Johnson Space Center, is a high fidelity three-dimensional physical model that is capable of simulating the historical satellite environment, as well as the evolution of future debris populations (14, 15). The subject study assumed no rocket bodies and spacecraft were launched after December 2004, and no future disposal maneuvers were allowed for existing spacecraft, few of which currently have such a capability. The rate of satellite explosions would naturally decrease to zero within a few decades as the current satellite population ages. The LEGEND future projection adopts a Monte Carlo approach to simulate future on-orbit explosions and collisions. Within a given projection time step, once the explosion probability is estimated for an intact object, a random number is drawn and compared with the probability to determine if an explosion would occur. A similar procedure is applied to collisions for each pair of target and projectile involved within the same time step. Due to the nature of the Monte Carlo process, multiple projection runs must be performed and analyzed before one can draw reliable and meaningful conclusions from the outcome. A total of fifty, 200-year future projection Monte Carlo simulations were executed and evaluated (16).

Liou, Jer-chyi↗

Timing of the Accreting Millisecond Pulsar IGR J17591-2342: Evidence of Spin-down during Accretion

We report on the phase-coherent timing analysis of the accreting millisecond X-ray pulsar IGR J17591–2342, using Neutron Star Interior Composition Explorer (NICER) data taken during the outburst of the source between 2018 August 15 and 2018 October 17. We obtain an updated orbital solution of the binary system. We investigate the evolution of the neutron star spin frequency during the outburst, reporting a refined estimate of the spin frequency and the first estimate of the spin frequency derivative ( ˙ν ∼ −7 × 10−14 Hz s−1), confirmed independently from the modelling of the fundamental frequency and its first harmonic. We further investigate the evolution of the X-ray pulse phases adopting a physical model that accounts for the accretion material torque as well as the magnetic threading of the accretion disc in regions where the Keplerian velocity is slower than the magnetosphere velocity. From this analysis we estimate the neutron star magnetic field Beq = 2.8(3) × 108 G. Finally, we investigate the pulse profile dependence on energy finding that the observed behaviour of the pulse fractional amplitude and lags as a function of energy is compatible with the down-scattering of hard X-ray photons in the disc or the neutron star surface.

A Sanna↗

NASA GPM GV Science Implementation

Pre-launch algorithm development & post-launch product evaluation: The GPM GV paradigm moves beyond traditional direct validation/comparison activities by incorporating improved algorithm physics & model applications (end-to-end validation) in the validation process. Three approaches: 1) National Network (surface): Operational networks to identify and resolve first order discrepancies (e.g., bias) between satellite and ground-based precipitation estimates. 2) Physical Process (vertical column): Cloud system and microphysical studies geared toward testing and refinement of physically-based retrieval algorithms. 3) Integrated (4-dimensional): Integration of satellite precipitation products into coupled prediction models to evaluate strengths/limitations of satellite precipitation producers.

Petersen, W. A.↗

A model of the origin of the Jovian ring

Assuming that the micron-sized particles making up the bright Jovian ring are fragments of erosive collisions between micrometeoroid projectiles and large parent bodies, a physical model of the ring is calculated. This leads to a well-defined size distribution for the ejecta, whose optical properties can be compared with observation. The (most likely silicate) ejecta material maximum diameter is estimated to be 0.1 micron, and most likely the result of Io volcanic activity. The impact model's determination of ejecta size distribution in turn determines the structure of the ring, with the largest ejecta forming the bright ring, medium-sized ejecta forming a disk that extends to the Jovian atmosphere, and small ejecta forming a faint halo whose structure is dominated by electromagnetic forces.

Gruen, E.↗

The X-ray surface brightness distribution and spectral properties of six early-type galaxies

Detailed analysis is presented of the Einstein X-ray observations of six early-type galaxies. The results show that effective cooling is probably present in these systems, at least in the innermost regions. Interaction with the surrounding medium has a major effect on the X-ray surface brightness distribution at large radii, at least for galaxies in clusters. The data do not warrant the general assumptions of isothermality and gravitational hydrostatic equilibrium at large radii. Comparison of the X-ray surface brightness profiles with model predictions indicate that 1/r-squared halos with masses of the order of 10 times the stellar masses are required to match the data. The physical model of White and Chevalier (1984) for steady cooling flows in a King law potential with no heavy halo gives a surface brightness distribution that resembles the data if supernovae heating is present.

Trinchieri, G.↗

Aerothermal Analysis of the Rocket Lab Venus Probe Heatshield

This document provides an overview of the aerothermodynamic analyses performed by the Aerothermodynamics Branch at NASA Langley Research Center for the Rocket Lab Venus Probe (RLVP). In addition to defining the baseline heating environment to the heatshield, this document pursues the experimental validation of key physical models at RLVP-relevant conditions. This experimental validation analysis, which captures the model form uncertainty, is used as one of two primary components of the margin assessment, where the other component is the parametric uncertainty. These model form (experimental) and parametric uncertainty components are used to construct a spatial and time varying margin for the heating to the RLVP heatshield. The margin is evaluated as the sum of the parametric and model form uncertainty components. The model form uncertainty is defined as the difference between the RLVP-relevant measurements and their simulations, using the upper limit uncertainty bounds for both the measurements and simulations in the comparisons. The differences in the dominant physics in the stagnation region and flank lead to the separate RLVP-relevant measurements for assessing the model form uncertainty in these two regions. These regions are addressed as follows: Stagnation Region Heating Environment: For the high-temperature stagnation-region, both the radiative heating and impact of blowing on convective heating are significant, while the impacts of turbulence and roughness are negligible. Coupled radiation and ablation LAURA/HARA solutions with ray-tracing provide the radiative heating over the entire vehicle, including the contributions from the Venus atmosphere and ablation species. Non-ablating LAURA simulations provide the convective heating. During the material-response computation typically used for TPS sizing, this non-ablating convective heating is corrected for the impact of ablation using the blowing correction. Coupled ablation LAURA simulations that capture finite-rate sur-face processes show that this blowing correction may be non-conservative over most of the heatshield. This non-conservatism is due to hydrogen recombination in the finite-rate surface model, which tends to increase the coupled ablation convective heating to near the non-ablating values, therefore making any reduction in the non-ablating value through the blowing correction non-conservative. This non-conservatism due to H catalysis is captured in the parametric component of the margin. The best available ground-test measurements that capture the impact of blowing on stagnation region convective heating, at RLVP-relevant conditions, indicate that the current blowing reduction model is non-conservative by up to 20% at RLVP-relevant blowing rates (the coupled blowing simulations were also non-conservative). Because of the relatively low velocity of the ground tests and the non-Venus atmospheric chemistry, these measurements do not capture the chemistry and therefore do not inform the uncertainty due to H catalysis. However, they do capture the fluid mechanics of blowing. The non-conservatism of the blowing correction implied by these measurements is covered by the model form component of the margin, which leads to total margin values over 50%. For the radiative heating, the shock-tube informed bias approach suggests a model form uncertainty of roughly 20%, while the parametric uncertainty analysis suggests values over 100%. The combined stagnation-point radiation margin of over 100% leads to peak margined radiative heating values of over300 W/cm2, which remains small relative to the peak margined convective heating of nearly 2000 W/cm2. Based on this analysis, at the stagnation point, the peak margined heat rate is 2203 W/cm2 and the margined total heat load is 31.5 kJ/cm2 for the current nominal trajectory. Flank Heating Environment: The forebody flank (and near-shoulder) heating environment is dominated by the impact of turbulence, roughness augmentation, and ablation on the convective heating. An extensive collection of ground test measurements with RLVP-relevant turbulence and roughness is studied to show that the maximum difference between the simulated and measured convective heating is 5%. However, with the exception of the Holden measurements from the 1980s, these measurements do not include roughness elements extending into the supersonic region of the boundary layer, which is likely to occur for RLVP (due to the 45 degree sphere-cone geometry). The interaction between the supersonic flow and roughness could cause convective heating augmentation fundamentally different than for locally subsonic flow. Although these Holden measurements are consistent with the other measurements considered, another path was pursued to assure that the rough-ness height extending into supersonic flow does not fundamentally change the roughness augmentation. This additional path was a computational effort to resolve the roughness elements in the CFD grid, so that the interaction be-tween the roughness elements and locally supersonic flow may be simulated in detail. This roughness-resolved CFD simulation is feasible because of the RLVP forebody TPS’s patterned roughness, which may be approximated analytically, and because of the axisymmetric nominal flow field, which allows a narrow surface region to be simulated and therefore make the computational expense feasible. These grid-resolved roughness simulations, which are performed at actual RLVP flight conditions, result in heating augmentation values that are below the design approach for roughness augmentation. This provides evidence that the design approach for RLVP roughness augmentation is sufficient. Based on this analysis, at this flank or near-shoulder location, the peak margined total heat rate is 2088 W/cm2and the margined total heat load is 26.0 kJ/cm2for the current nominal trajectory. Heat flux, shear, pressure and heat transfer coefficient at the RLVP stagnation point and near shoulder location are evaluated for the entire trajectory, and curved fit to a functional form of F=AρB∞UC∞. These simplified relationships for the nominal and margined aerothermal environments are referred to as aerothermal indicators, and presented at the end of this document.

Christopher O Johnston↗

Measurement of material parameters that limit the open-circuit voltage in P-N-junction silicon solar cells

The greatest gains in solar energy conversion efficiency of p-n-junction silicon solar cells come from increasing the open-circuit voltage V sub OC; it is important to understand and characterize the material parameters that limit the V sub OC. Strong experimental evidence exists to support the assertion that either an anomalously large minority carrier charge storage or an anomalously small minority carrier lifetime in the quasi-neutral emitter region limits the open circuit voltage. A method is presented for measuring charge storage and effective lifetime. Static and transient measurements are analyzed using physical models of the solar cell characteristics. This analysis yields the emitter charge storage and life-time, which then can be related to the various physical mechanisms, such as energy band gap shrinkage, that have been proposed earlier as responsible for limiting V sub OC.

Lindholm, F. A.↗

Two-fluid simulations of galaxy formation

We investigate the formation of galaxies and larger structure with a simulation modeling two gravitationally coupled fluids representing dark matter and baryons. The baryon gas dynamics are calculated with a smoothed particle hydrodynamics (SPH) method, and the physics modeled includes thermal pressure, shock heating, and radiative cooling. We simulate a 16 Mpc periodic cube with 64(exp 3) particles in each fluid and 10% baryon mass fraction. We confirm, for the first time experimentally, disk formation as a natural consequence of hierarchical clustering in a large-scale cosmological environment. The majority of isolated galaxies exhibit centrifugally supported disks. A power-law relation between cold baryonic mass and maximum rotation velocity is found, M varies as nu(sub rot)(exp alpha) with alpha = 2.5 after correcting for differential numerical resolution. Both the spatial and velocity distributions of the simulated galaxies are biased with respect to the dark matter. A counts-in-cells analysis indicates that an unphysical degree of merging in the central cluster is likely responsible for the antibias signal in the correlation function. A robust, scale-dependent velocity bias is measured. The ratio of galaxy to dark matter pairwise velocity dispersions on a scale of 1 Mpc is 0.7. The amplitude is only mildly dependent on redshift or mass cutoff and scales with separation as r(exp 0.2). The degree to which these results depend on numerical parameters is discussed. Mass resolution plays a key role in controlling the resulting fraction of cold, dense baryons. The mass fraction associated with galaxies decreases by a factor of approximately greater than 3 when the mass per particle is increased by a factor 8. Photoionization and energy input from supernova will have to be included to determine more carefully the fraction of highly dissipated material and the characteristics of the stellar component of galaxies.

Evrard, August E.↗

Models of the upper atmosphere

Several computer models of the thermosphere are discussed. J70MM, the Jacchia 1970 model with matrix and matrix mean output, and J703X, the Jacchia 1970, 1971, and 1977 models, are of primary interest. The subroutines in these programs were studied in detail, and several mistakes were found and corrected. It is proposed that a weighted average of the three-hour geomagnetic index be used in the models rather than a single index value. The densities (for a given date and time) generated by the J703X program (all three models), using both a single value of the index and weighted averages for three different lag times, are tabulated and discussed. Most of the equations used in the models are the results of empirical curve fits. An attempt was made to generate a theoretical prediction of the thermospheric temperature profile based on a simple physical model of atmospheric heat conduction in a spherically-symmetric shell. The exponential heating theory and the temperature dependence of k are discussed.

Davis, T. M.↗