Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “numerical continuation”

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

The NASA Ames Research Center one- and two-dimensional stratospheric models. Part 1: The one-dimensional model

A one-dimensional model of stratospheric trace constituents is described in detail. Specifically, the numerical solution of the species continuity equations, including a technique for treating the stiff differential equations representing the chemical kinetic terms, and an appropriate method for simulating the diurnal variations of the species concentrations, are discussed. A specialized treatment of atmospheric photodissociation rates is outlined in the text. The choice of a vertical eddy diffusivity profile and its success in predicting the vertical tracer distributions (carbon 14, methane, and nitrous oxide) are also discussed.

Turco, R. P.↗

The NASA Ames Research Center one- and two-dimensional stratospheric models. Part 2: The two-dimensional model

The two-dimensional model of stratospheric constituents is presented in detail. The derivation of pertinent transport parameters and the numerical solution of the species continuity equations, including a technique for treating the stiff differential equations that represent the chemical kinetic terms, and appropriate methods for simulating the diurnal variations of the solar zenith angle and species concentrations are discussed. Predicted distributions of tracer constituents (ozone, carbon 14, nitric acid) are compared with observed distributions.

Whitten, R. C.↗

An analytic formulation for heating source memory in the thermospheric composition

Numerical calculations from a spectral circulation model are utilized to construct an analytic Green's function formulation describing the meridional, time-dependent thermospheric composition and temperature response during magnetic storms. The purpose is to develop a formulation that embodies source memory while being sufficiently simple to serve as a heuristic guide for empirical modeling. By passing from the discrete Fourier series representation, utilized for the numerical circulation model, to a continuous Fourier integral representation, explicit waves are obtained for the thermospheric response times. The response times are altitude and species dependent and can exceed two days below 200 km. Thus, for certain storm scenarios, pronounced source memory signatures for the composition and temperature are predicted. Response times obtained from the formulation are shown to give a response consistent with previously published neutral composition data from AE-C for the February 1974 storm when an ap dependent heat source is employed.

Porter, H. S.↗

A method for extracting meridonal winds from ionosonde measurements by using ionospheric models

There has been great progress in modelling and measuring the dynamics of the neutral upper atmosphere in recent years. However, future progress will depend on the availability of global measurements of neutral winds. Attention is drawn to a relatively cheap means of supplementing the data base of neutral winds provided by radar and optical measurements with data obtained by ionosondes. Rishbeth in his review of F-region dynamics, derived the relationship between the height of the F2 layer and the component of the neutral wind parallel to the magnetic field of the Earth. The sensitivity is examined of the height and density of the F2 layer over Boulder, Colorado on 30 July 1982 to changes in meridional wind speed using a comprehensive interhemispheric numerical model that solves the continuity and momentum equations for H+ and O+, the energy equations for Te and Ti, and the 2-stream photoelectron equations to obtain electron heating rates. For the neutral atmosphere temperature and densities, the mass spectrometer incoherent scatter radar model of Hedin was used.

Richards, P. G.↗

Aspects of the determination of winds by means of scatterometry and of the utilization of vector wind data for meteorological forecasts

The present paper provides a description of four aspects of scatterometer winds and their uses. The theory of wave generation by the wind is considered along with an analysis of the properties of superobservations, and studies of intermittent versus continuous data assimilation methods for numerical weather predictions which use remotely sensed data. A comparison of the sum of squares versus the maximum likelihood method for recovering the vector winds is also conducted. Questions regarding wind speed, friction velocity, or normal stress are discussed and synoptic scale fields from Seasat-SASS data are examined.

Pierson, W. J., Jr.↗

A mathematical model of the middle and high latitude ionosphere

A time-dependent three-dimensional model of the middle and high latitude ionosphere is described. The density distributions of six ion species NO(+), O(2+), N(2+), O(+), N(+), He(+), and the electron and ion temperatures are obtained from a numerical solution of the appropriate continuity, momentum, and energy equations. The equations are solved as a function of height for an inclined magnetic field at E and F region altitudes. The three-dimensional nature of the model is obtained by following flux tubes of plasma as they convect or corotate through a moving neutral atmosphere. The model takes account of field-aligned diffusion, cross-field electrodynamic drifts, thermospheric winds, polar wind escape, energy-dependent chemical reactions, neutral composition changes, ion production due to solar EUV radiation and auroral precipitation, thermal conduction, diffusion-thermal heat flow and local heating and cooling processes. The model also takes account of the offset between the geomagnetic and geographic poles.

Schunk, R. W.↗

A new class of solutions for interstellar magnetohydrodynamic shock waves

An analysis is presented of the equations of motion for steady MHD shock waves proopagating in interstellar clouds, for boundary conditions that preclude C shocks. In addition to J shocks, in which the neutral fluid component becomes subsonic at an adiabatic jump front, the equations admit a new class of solutions, called C-asterisk shocks, in which the transition to subsonic flow occurs continuously at a sonic point. Numerical methods are developed for computing the structure of J and C-asterisk shocks propagating in diffuse interstellar clouds. The effects of chemical, ionization, and recombination processes are included in this treatment. An alternative numerical method, which uses artificial viscosity to facilitate integration through sonic points, is analyzed and shown to be invalid. A set of exemplary solutions, computed for realistic shock parameters, shows that C-asterisk shocks occur for a broad range of conditions relevant to diffuse interstellar clouds.

Roberge, W. G.↗

Penetration experiments in aluminum and Teflon targets of widely variable thickness

The morphologies and detailed dimensions of hypervelocity craters and penetration holes on space-exposed surfaces faithfully reflect the initial impact conditions. However, current understanding of this postmortem evidence and its relation to such first-order parameters as impact velocity or projectile size and mass is incomplete. While considerable progress is being made in the numerical simulation of impact events, continued impact simulations in the laboratory are needed to obtain empirical constraints and insights. This contribution summarizes such experiments with Al and Teflon targets that were carried out in order to provide a better understanding of the crater and penetration holes reported from the Solar Maximum Mission (SMM) and the Long Duration Exposure Facility (LDEF) satellites. A 5-mm light gas gun was used to fire spherical soda-lime glass projectiles from 50 to 3175 microns in diameter (D(sub P)), at a nominal 6 km/s, into Al (1100 series; annealed) and Teflon (Teflon(sup TFE)) targets. Targets ranged in thickness (T) from infinite halfspace targets (T approx. equals cm) to ultrathin foils (T approx. equals micron), yielding up to 3 degrees of magnitude variation in absolute and relative (D(sub P)/T) target thickness. This experimental matrix simulates the wide range in D(sub P)/T experienced by a space-exposed membrane of constant T that is being impacted by projectiles of widely varying sizes.

Hoerz, F.↗

The Axisymmetric Pulsar Magnetosphere

We present the structure of the axisymmetric force-free magnetosphere of an aligned rotating magnetic dipole, in the case in which there exists a sufficiently large charge density (whose origin we do not question) to satisfy the ideal MHD condition, E (raised dot) B = 0, everywhere. The unique distribution of electric current along the open magnetic field lines which is required for the solution to be continuous and smooth is obtained numerically. We expect that our solution will be useful as the starting point for detailed studies of pulsar magnetospheres under more general conditions, namely when either the force-free and/or the ideal MHD condition E (raised dot) B = 0 are not valid in the entire magnetosphere. Based on our solution, we consider that the most likely positions of such an occurrence are the polar cap, the crossings of the zero space charge surface by open field lines, and the return current boundary, but not the light cylinder.

Contopoulos, J.↗

Investigating the Effects of Temperature on the Signatures of Shocks Propagated Through Impacts into Minerals Found in Comets and Asteroids

Comets and asteroids are subjected to extremely cold conditions throughout their lifetimes. During their sojourns in the solar system, they are subjected to collisions at speeds that are easily capable of generating shock waves in their constituent materials. In addition to ices, more common silicate minerals such as olivines and pyroxenes are important components of these objects. The collision-induced shocks could affect the spectral signatures of those mineral components, which could in turn be detected telescopically. We have embarked on a project to determine how impact-generated shock might affect the reflectance spectra and structures of select silicates as both impact speed and target temperature are varied systematically. While the effects of impact speed (in the form of shock stress) on numerous materials have been and continue to be studied, the role of target temperature has received comparatively little attention, presumably because of the operational difficulties it can introduce to experimentation. Our experiments were performed with the vertical gun in the Experimental Impact Laboratory of the Johnson Space Center. A liquid-nitrogen system was plumbed to permit cooling of the target container and its contents under vacuum to temperatures as low as -100 C (173 K). Temperatures were monitored by thermocouples mounted on the outside of the target container. Because those sensors were not in contact with the target material at impact, the measured temperatures are treated as lower limits for the actual values. Peridot (Mg-rich olivine) and enstatite (Mg-rich orthopyroxene) were used as targets, which involved the impact of alumina (Al2O3) spheres at speeds of 2.0 - 2.7 km s(exp -1) and temperatures covering 25 C to -100 C (298 K to 173 K). We have begun collecting and analyzing data in the near to mid-IR with a Fourier-transform infrared spectrometer, and preliminary analyses show that notable differences in absorption-band strength and position occur as functions of both impact speed (peak shock stress) and initial temperature.

Lederer, Susan M.↗

Multidisciplinary Optimization of a Transport Aircraft Wing using Particle Swarm Optimization

The purpose of this paper is to demonstrate the application of particle swarm optimization to a realistic multidisciplinary optimization test problem. The paper's new contributions to multidisciplinary optimization is the application of a new algorithm for dealing with the unique challenges associated with multidisciplinary optimization problems, and recommendations as to the utility of the algorithm in future multidisciplinary optimization applications. The selected example is a bi-level optimization problem that demonstrates severe numerical noise and has a combination of continuous and truly discrete design variables. The use of traditional gradient-based optimization algorithms is thus not practical. The numerical results presented indicate that the particle swarm optimization algorithm is able to reliably find the optimum design for the problem presented here. The algorithm is capable of dealing with the unique challenges posed by multidisciplinary optimization as well as the numerical noise and truly discrete variables present in the current example problem.

Sobieszczanski-Sobieski, Jaroslaw↗

Solar Thermal Concept Evaluation

Concentrated solar thermal energy can be utilized in a variety of high temperature applications for both terrestrial and space environments. In each application, knowledge of the collector and absorber's heat exchange interaction is required. To understand this coupled mechanism, various concentrator types and geometries, as well as, their relationship to the physical absorber mechanics were investigated. To conduct experimental tests various parts of a 5,000 watt, thermal concentrator, facility were made and evaluated. This was in anticipation at a larger NASA facility proposed for construction. Although much of the work centered on solar thermal propulsion for an upper stage (less than one pound thrust range), the information generated and the facility's capabilities are applicable to material processing, power generation and similar uses. The numerical calculations used to design the laboratory mirror and the procedure for evaluating other solar collectors are presented here. The mirror design is based on a hexagonal faceted system, which uses a spherical approximation to the parabolic surface. The work began with a few two dimensional estimates and continued with a full, three dimensional, numerical algorithm written in FORTRAN code. This was compared to a full geometry, ray trace program, BEAM 4, which optimizes the curvatures, based on purely optical considerations. Founded on numerical results, the characteristics of a faceted concentrator were construed. The numerical methodologies themselves were evaluated and categorized. As a result, the three-dimensional FORTRAN code was the method chosen to construct the mirrors, due to its overall accuracy and superior results to the ray trace program. This information is being used to fabricate and subsequently, laser map the actual mirror surfaces. Evaluation of concentrator mirrors, thermal applications and scaling the results of the 10 foot diameter mirror to a much larger concentrator, were studied. Evaluations, recommendations and pit falls regarding the structure, materials and facility design are presented.

Hawk, Clark W.↗

Recent Advances In Science Support For Isolated Droplet Combustion Experiments

In a joint program involving Prof. F.A. Williams of the University of California, San Diego and Dr. V. Nayagam of the National Center for Microgravity Research, the combustion characteristics of isolated liquid fuel droplets of n-heptane, n-decane, methanol, methanol-water, ethanol and ethanol-water having initial diameters between about 1 mm and 6 mm continues to be investigated. The objectives of the work are to improve fundamental knowledge of droplet combustion dynamics for pure fuels and fuel-water mixtures through microgravity experiments and theoretical analyses. The Princeton contributions support the engineering design, data analysis, and data interpretation requirements for the study of initially single component, spherically symmetric, isolated droplet combustion studies through experiments and numerical modeling. UCSD contributions are described in a companion communication in this conference. The Princeton effort also addresses the analyses of Fiber Supported Droplet Combustion (FSDC) experiments conducted with the above fuels and collaborative work with others who are investigating droplet combustion in the presence of steady convection. A thorough interpretation of droplet burning behavior for n-heptane and n-decane over a relatively wide range of conditions also involves the influences of sooting on the combustion behavior, and this particular aspect on isolated burning of droplets is under consideration in a collaborative program underway with Drexel University. This collaboration is addressed in another communication at this conference. The one-dimensional, time-dependent, numerical modeling approach that we have continued to evolve for analyzing isolated, quiescent droplet combustion data has been further applied to investigate several facets of isolated droplet burning of simple alcohols, n-heptane, and n-decane. Some of the new results are described below.

Dryer, F. L.↗

Effects of sounding temperature assimilation on weather forecasting - Model dependence studies

In comparing various methods for the assimilation of remote sounding information into numerical weather prediction (NWP) models, the problem of model dependence for the different results obtained becomes important. The paper investigates two aspects of the model dependence question: (1) the effect of increasing horizontal resolution within a given model on the assimilation of sounding data, and (2) the effect of using two entirely different models with the same assimilation method and sounding data. Tentative conclusions reached are: first, that model improvement as exemplified by increased resolution, can act in the same direction as judicious 4-D assimilation of remote sounding information, to improve 2-3 day numerical weather forecasts. Second, that the time continuous 4-D methods developed at GLAS have similar beneficial effects when used in the assimilation of remote sounding information into NWP models with very different numerical and physical characteristics.

Ghil, M.↗

A theoretical framework for convergence and continuous dependence of estimates in inverse problems for distributed parameter systems

Numerical techniques for parameter identification in distributed-parameter systems are developed analytically. A general convergence and stability framework (for continuous dependence on observations) is derived for first-order systems on the basis of (1) a weak formulation in terms of sesquilinear forms and (2) the resolvent convergence form of the Trotter-Kato approximation. The extension of this framework to second-order systems is considered.

Banks, H. T.↗

Dual-Mode Combustion of Hydrogen in a Mach 5, Continuous-Flow Facility

Results of an experimental and numerical study of a dual-mode scramjet combustor are reported. The experiment consisted of a direct-connect test of a Mach 2 hydrogen-air combustor with a single unswept-ramp fuel injector. The flow stagnation enthalpy simulated a flight Mach number of 5. Measurements were obtained using conventional wall instrumentation and a particle-imaging laser diagnostic technique. The particle imaging was enabled through the development of a new apparatus for seeding fine silicon dioxide particles into the combustor fuel stream. Numerical simulations of the combustor were performed using the GASP code. The modeling, and much of the experimental work, focused on the supersonic combustion mode. Reasonable agreement was observed between experimental and numerical wall pressure distributions. However, the numerical model was unable to predict accurately the effects of combustion on the fuel plume size, penetration, shape, and axial growth.

Goyne, C. P.↗