An absorption integral for planetary atmospheres.
Absorption integral for isothermal planetary atmospheres with exponential pressure distribution
SEARCH · Search NASA
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.
Absorption integral for isothermal planetary atmospheres with exponential pressure distribution
Problems in the design of feedback controllers for chaotic dynamical systems are considered theoretically, focusing on two cases where chaos arises only when a nonchaotic continuous-time system is discretized into a simpler discrete-time systems (exponential discretization and pseudo-Euler integration applied to Lotka-Volterra competition and prey-predator systems). Numerical simulation results are presented in extensive graphs and discussed in detail. It is concluded that care must be taken in applying standard dynamical-systems methods to control systems that may be discontinuous or nondifferentiable.
When incorporating formulations of instantaneous solar heating or photolytic rates as functions of altitude and sun angle into long range forecasting models, it may be desirable to replace the time integrals by daily average rates that are simple functions of latitude and season. This replacement is accomplished by approximating the integral over the solar day by a pure exponential. This gives a daily average rate as a multiplication factor times the instantaneous rate evaluated at an appropriate sun angle. The accuracy of the exponential approximation is investigated by a sample calculation using an instantaneous ozone heating formulation available in the literature.
Adopting the formalism of Parsons and Daly (1983), analytical integral equations (Green's function integrals) are derived which relate gravity anomalies and dynamic boundary topography with temperature as a function of wavenumber for a fluid layer whose viscosity varies exponentially with depth. In the earth, such a viscosity profile may be found in the asthenosphere, where the large thermal gradient leads to exponential decrease of viscosity with depth, the effects of a pressure increase being small in comparison. It is shown that, when viscosity varies rapidly, topography kernels for both the surface and bottom boundaries (and hence the gravity kernel) are strongly affected at all wavelengths.
A new time integrator for large, stiff systems of linear and nonlinear coupled differential equations is described. For linear systems, the method consists of forming a small (5-15-term) Krylov space using the Jacobian of the system and carrying out exact exponential propagation within this space. Nonlinear corrections are incorporated via a convolution integral formalism; the integral is evaluated via approximate Krylov methods as well. Gains in efficiency ranging from factors of 2 to 30 are demonstrated for several test problems as compared to a forward Euler scheme and to the integration package LSODE.
A computer-oriented technique is presented for performing a nonlinear exponential regression analysis on decay-type experimental data. The technique involves the least squares procedure wherein the nonlinear problem is linearized by expansion in a Taylor series. A linear curve fitting procedure for determining the initial nominal estimates for the unknown exponential model parameters is included as an integral part of the technique. A correction matrix was derived and then applied to the nominal estimate to produce an improved set of model parameters. The solution cycle is repeated until some predetermined criterion is satisfied.
Axial and off-axis radiative cooling of cylindrical shock-heated hydrogen-helium plasmas is investigated theoretically and experimentally. The coupled fluid dynamic-radiative transfer equations are solved by a combination of approximation techniques aimed at simplifying the computation of the flux divergence term, namely, the quasi-isothermal approximation and the exponential approximation developed for the solid angle integration. The accuracy of the approximation schemes has been assessed and found acceptable for applying the methods to the rapid computation of the radiatively coupled flow problem. Radiative cooling experiments were conducted in a 6-inch annular arc accelerator shock tube (ANAA) for an initial pressure of 1 torr and shock speeds from 35 to 45 Km/sec. The results indicate that the lateral cooling is small compared with the axial cooling, and that better agreement is achieved between the data and the theoretical results by inclusion of the lateral temperature gradient.
Transonic aeroelastic solutions based upon the transonic small perturbation potential equation were studied. Time-marching transient solutions of plunging and pitching airfoils were analyzed using a complex exponential modal identification technique, and seven alternative integration techniques for the structural equations were evaluated. The HYTRAN2 code was used to determine transonic flutter boundaries versus Mach number and angle-of-attack for NACA 64A010 and MBB A-3 airfoils. In the code, a monotone differencing method, which eliminates leading edge expansion shocks, is used to solve the potential equation. When the effect of static pitching moment upon the angle-of-attack is included, the MBB A-3 airfoil can have multiple flutter speeds at a given Mach number.
Transonic aeroelastic solutions based upon the transonic small perturbation potential equation are studied. Time-marching transient solutions of plunging and pitching airfoils are analyzed using a complex exponential modal identification technique, and seven alternative integration techniques for the structural equations are evaluated. The HYTRAN2 code is used to determine transonic flutter boundaries versus Mach number and angle-of-attack for NACA 64A010 and MBB A-3 airfoils. In that code, a monotone differencing method, which eliminates leading edge expansion shocks, is used to solve the potential equation. When the effect of static pitching moment upon the angle-of-attack is included, the MBB A-3 airfoil can have multiple flutter speeds at a given Mach number.
Verification and Validation (V&V) at JPL is traditionally performed on flight or flight-like hardware running flight software. For some time, the complexity of avionics has increased exponentially while the time allocated for system integration and associated V&V testing has remained fixed. There is an increasing need to perform comprehensive system level V&V using modeling and simulation, and to use scarce hardware testing time to validate models; the norm for thermal and structural V&V for some time. Our approach extends model-based V&V to electronics and software through functional and structural models implemented in SysML. We develop component models of electronics and software that are validated by comparison with test results from actual equipment. The models are then simulated enabling a more complete set of test cases than possible on flight hardware. SysML simulations provide access and control of internal nodes that may not be available in physical systems. This is particularly helpful in testing fault protection behaviors when injecting faults is either not possible or potentially damaging to the hardware. We can also model both hardware and software behaviors in SysML, which allows us to simulate hardware and software interactions. With an integrated model and simulation capability we can evaluate the hardware and software interactions and identify problems sooner. The primary missing piece is validating SysML model correctness against hardware; this experiment demonstrated such an approach is possible.
Hard X-ray measurements of six gamma-ray bursts observed by detectors on the OGO-5 and OSO-6 satellites during the period Oct. 1969 to Apr. 1971 are presented. Spectra for five of the six bursts were determined using measurements from both satellites in order to reduce ambiguities due to uncertain source locations. A significant fraction, 20 to 60%, of the energy of the bursts fall in the hard X-ray range (20 to 130 keV). The time-integrated spectra were fitted by power-law, exponential, and thermal bremsstrahlung functions. They are consistent with power-laws which steepen at energies approx. above 150 keV, as reported earlier for two other bursts. Evidence for spectral variability from event to event in the hard X-ray region is presented. The hard X-ray spectra of the gamma-ray bursts differ from those of the recently discovered 1 to 15 keV bursts.
The transport properties of the oxygen system, O2, O, O(+), and the electron, e, are calculated and their applications are discussed. The possible two-body interactions in the oxygen system are examined and the exponential repulsive potential is used to determine collision integrals for these interactions. Values are presented for the viscosity of the pure gases, the self-diffusion coefficient of the pure gases at 1 atm pressure, the translational thermal conductivity of the pure gases, and the internal thermal conductivity of O2.
Numerical cloud models have been developed and applied extensively to study cloud-scale and mesoscale processes during the past four decades. The distinctive aspect of these cloud models is their ability to treat explicitly (or resolve) cloud-scale dynamics. This requires the cloud models to be formulated from the non-hydrostatic equations of motion that explicitly include the vertical acceleration terms since the vertical and horizontal scales of convection are similar. Such models are also necessary in order to allow gravity waves, such as those triggered by clouds, to be resolved explicitly. In contrast, the hydrostatic approximation, usually applied in global or regional models, does allow the presence of gravity waves. In addition, the availability of exponentially increasing computer capabilities has resulted in time integrations increasing from hours to days, domain grids boxes (points) increasing from less than 2000 to more than 2,500,000 grid points with 500 to 1000 m resolution, and 3-D models becoming increasingly prevalent. The cloud resolving model is now at a stage where it can provide reasonably accurate statistical information of the sub-grid, cloud-resolving processes poorly parameterized in climate models and numerical prediction models.
In a self-gravitating system of point particles such as a spherical star cluster, small disturbances to an orbit grow exponentially on a time-scale comparable with the crossing time. The results of N-body integrations are therefore extremely sensitive to numerical errors: in practice it is almost impossible to follow orbits of individual particles accurately for more than a few crossing times. We demonstrate that numerical orbits in the gravitational N-body problem are often shadowed by true orbits for many crossing times. This result enhances our confidence in the use of N-body integrations to study the evolution of stellar systems.
It is known that the exact analytic solutions of wave scattering by a circular cylinder, when they exist, are not in a closed form but in infinite series which converges slowly for high frequency waves. In this paper, we present a fast number solution for the scattering problem in which the boundary integral equations, reformulated from the Helmholtz equation, are solved using a Fourier spectral method. It is shown that the special geometry considered here allows the implementation of the spectral method to be simple and very efficient. The present method differs from previous approaches in that the singularities of the integral kernels are removed and dealt with accurately. The proposed method preserves the spectral accuracy and is shown to have an exponential rate of convergence. Aspects of efficient implementation using FFT are discussed. Moreover, the boundary integral equations of combined single and double-layer representation are used in the present paper. This ensures the uniqueness of the numerical solution for the scattering problem at all frequencies. Although a strongly singular kernel is encountered for the Neumann boundary conditions, we show that the hypersingularity can be handled easily in the spectral method. Numerical examples that demonstrate the validity of the method are also presented.
Hard X-ray measurements of six gamma-ray bursts observed during the period from October 1969 to April 1971 are presented. The measurements were made with detectors on the OGO-5 and OSO-6 satellites. Spectra for five of the six bursts have been determined using measurements from both satellites in order to reduce ambiguities due to uncertain source locations. A significant fraction, about 20-60%, of the energy of the bursts falls in the hard X-ray range (20-130 keV). The time-integrated spectra have been fitted by power-law, exponential, and thermal-bremsstrahlung functions. They are consistent with power laws which steepen at energies of at least 150 keV, as reported earlier for two other bursts. Evidence for spectral variability from event to event in the hard X-ray region is presented. For a power-law representation, the power-law index has values ranging from approximately unity to 2.5. The hard X-ray spectra of the gamma-ray bursts differ significantly from those of the recently discovered 1-15-keV X-ray bursts.
In recent years, exponentially increasing computer power extended Cloud Resolving Model (CRM) integrations from hours to months, the number of computational grid points from less than a thousand to close to ten million. Three-dimensional models are now more prevalent. Much attention is devoted to precipitating cloud systems where the crucial 1-km scales are resolved in horizontal domains as large as 10,000 km in two-dimensions, and 1,000 x 1,000 sq km in three-dimensions. Cloud resolving models now provide statistical information useful for developing more realistic physically based parameterizations for climate models and numerical weather prediction models. It is also expected that NWP and mesoscale models can be run in grid size similar to cloud resolving models through nesting technique. Recently, a multi-scale modeling system with unified physics was developed at NASA Goddard. It consists of (1) a cloud-resolving model (Goddard Cumulus Ensemble model, GCE model). (2) a regional scale model (a NASA unified weather research and forecast, W8F). (3) a coupled CRM and global model (Goddard Multi-scale Modeling Framework, MMF), and (4) a land modeling system. The same microphysical processes, long and short wave radiative transfer and land processes and the explicit cloud-radiation and cloud-land surface interactive processes are applied in this multi-scale modeling system. This modeling system has been coupled with a multi-satellite simulator to use NASA high-resolution satellite data to identify the strengths and weaknesses of cloud and precipitation processes simulated by the model. In this talk, a review of developments and applications of the multi-scale modeling system will be presented. In particular, the results from using multi-scale modeling systems to study the interactions between clouds, precipitation, and aerosols will be presented. Also how to use the multi-satellite simulator to improve precipitation processes will be discussed.
An integral equation is developed for application to time-dependent laboratory experiments in which partial redistribution effects are important. The equation of transport with the Heasly-Kneer (1976) emission coefficient and the equation of statistical equilibrium lead to a time-dependent redistribution function containing an absorption - reemission term which decays exponentially in time and a scattering term which is instantaneous. This integral equation does not agree with an equation written by Payne et al. (1974) that has been used to compare theory with experiments. The difference between the Payne equation and the equation developed here needs to be examined in detail, since it might under some circumstances be on the same order as the difference between partial and complete redistribution.