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 289 records · Page 16

Theory of Electrostatic Fields in the Ionosphere at Equatorial Latitudes

The properties of the elongated electrostatic fields that are required to provide the coupling mechanism in the dynamo-motor concept of the E and F regions of the ionosphere are examined theoretically for the conditions that prevail in equatorial latitudes. The analysis is developed for an electrostatic field of arbitrary horizontal scale in a horizontally stratified partly ionized gas subject to an imposed magnetic field having the form of a parabolic arch over the equator. The anisotropic character and continuous variation with height of the conductivity are retained throughout , and numerical solutions are determined for the attenuation of the electric field with distance along the field line. The results are similar qualitatively to those found previously upon analysis of the corresponding problem for middle latitudes, but the attenuation of the electrostatic field with height is considerably greater. It develops, in particular, that the coupling between E and F regions is very small for fields having horizontal wavelengths of a few kilometers. It follows that the dynamo-motor concept could not be used to account for the presence of irregularities of this scale. On the other hand, it is found that almost all the attenuation occurs at heights near that of the E region. If an electrostatic field having a horizontal wavelength of a few kilometers could be produced at heights of, say 200 km or greater, it follows that it would be very elongated and extend with little change in amplitude from hemisphere to hemisphere. The results also indicate that effective coupling could be achieved for fields having horizontal wavelengths of the order of tens of kilometers, but the assumptions introduced to simplify the analysis may impair the quantitative reliability of the results for fields of this scale.

Briggs, Benjamin R.↗

Computations and turbulent flow modeling in support of helicopter rotor technology

Efforts continue to investigate the applicability of the Deiwert's time dependent numerical airfoil calculation code to the simulation of two dimensional airfoil flows with large amounts of separation present. A modified logarithmic region mixing length term produced substantial differences in the calculated flow fields for the large separation of the M=0.4, high angle of attack cases. Reasons are given for interest in applying that model to previously investigated cases of the 64A010 airfoil section at Mach number 0.8 for angles of attack of 4 deg and 6.2 deg. The effects of these turbulence models at alterations and the time dependent behavior from impulsively started conditions are discussed.

Rose, W. C.↗

Aeroassist flight experiment aerodynamics and aerothermodynamics

The problem is to determine the transitional flow aerodynamics and aerothermodynamics, including the base flow characteristics, of the Aeroassist Flight Experiment (AFE). The justification for the computational fluid dynamic (CFD) Application stems from MSFC's system integration responsibility for the AFE. To insure that the AFE objectives are met, MSFC must understand the limitations and uncertainties of the design data. Perhaps the only method capable of handling the complex physics of the rarefied high energy AFE trajectory is Bird's Direct Simulation Monte Carlo (DSMC) technique. The 3-D code used in this analysis is applicable only to the AFE geometry. It uses the Variable Hard Sphere (VHS) collision model and five specie chemistry model available from Langley Research Center. The code is benchmarked against the AFE flight data and used as an Aeroassisted Space Transfer Vehicle (ASTV) design tool. The code is being used to understand the AFE flow field and verify or modify existing design data. Continued application to lower altitudes is testing the capability of the Numerical Aerodynamic Simulation Facility (NASF) to handle 3-D DSMC and its practicality as an ASTV/AFE design tool.

Brewer, Edwin B.↗

Low frequency electric and magnetic fields

Following preliminary investigations of the low frequency electric and magnetic fields that may exists in the Earth-ionospheric cavity, measurements were taken with state-of-the art spectrum analyzers. As a follow up to this activity, an investigation was initiated to determine sources and values for possible low frequency signal that would appear in the cavity. The lowest cavity resonance is estimated at about 8 Hz, but lower frequencies may be an important component of our electromagnetic environment. The potential field frequencies produced by the electron were investigated by a classical model that included possible cross coupling of the electric and gravitation fields. During this work, an interesting relationship was found that related the high frequency charge field with the extremely low frequency of the gravitation field. The results of numerical calculations were surprisingly accurate and this area of investigation is continuing. The work toward continued development of a standardized monitoring facility is continuing with the potential of installing the prototype at West Virginia State College early in 1990. This installation would be capable of real time monitoring of ELF signals in the Earth-ionoshpere cavity and would provide some directional information. A high gain, low noise, 1/f frequency corrected preamplifier was designed and tested for the ferrite core magnetic sensor. The potential application of a super conducting sensor for the ELF magnetic field detection is under investigation. It is hoped that a fully operational monitoring network could pinpoint the location of ELF signal sources and provide new information on where these signals originate and what causes them, assuming that they are natural in origin.

Spaniol, Craig↗

Numerical methods for the analysis of sampled-data systems and for the computation of system zeros

MARSYAS is a computer aided control system design package for the simulation and analysis of open loop and closed loop dynamic systems. Outlined here is the numerical and theoretical basis of the MARSYAS functions developed during the summer of 1991. In particular, the numerical computation of the matrix exponential e(exp A) = I + A + A(exp 2)/2! + A(exp 3)/3! + ... and the numerical computation of the finite system zeros of a dynamic system (continuous time or discrete time) are presented.

Hodel, A. Scottedward↗

Lidar Atmospheric Observatory in the Canadian Arctic

During the last decade there have been growing concerns about a broad variety of atmospheric properties. Among these, a depletion of the stratospheric ozone layer has attracted considerable attention from the general public, politicians and scientists due to its vital impact for the entire global biosphere. One of the major warning signs was the discovery of the 'ozone hole' in the Antarctic region where the concentration of the ozone in the stratosphere was significantly reduced. At present the stratospheric ozone layer in this region is being continuously monitored by groups of scientists from around the world and numerous observations of the ozone layer on the global scale have clearly demonstrated the process of ozone depletion. Recent observations by NASA have shown significant depletion in the Arctic region. This paper provides an initial description of two lidars that are planned to be installed in a new observatory for atmospheric studies in the Canadian Arctic. This observatory is being constructed under the supervision of the Atmospheric Environment Services (AES) of Canada as a part of Green Plan - an initiative of the Federal Government of Canada. The station is located at Eureka on Ellesmere Island at a latitude of 80 degrees N and a longitude of 86 degrees W.

Ulitsky, Arkady↗

A boundary integral method for an inverse problem in thermal imaging

An inverse problem in thermal imaging involving the recovery of a void in a material from its surface temperature response to external heating is examined. Uniqueness and continuous dependence results for the inverse problem are demonstrated, and a numerical method for its solution is developed. This method is based on an optimization approach, coupled with a boundary integral equation formulation of the forward heat conduction problem. Some convergence results for the method are proved, and several examples are presented using computationally generated data.

Bryan, Kurt↗

Needs and status of CFD code validation

The two types of Computational Fluid Dynamics code validations, solution-to-solution comparison and solution-to-experiment comparison, are discussed. It is suggested that to develop more detailed experiments the following things are necessary: (1) further development of turbulence models; (2) better methods for numerical validation of CFD codes; (3) evaluation of disagreements; and (4) continued determination of experimental scatter. All data and results are presented in viewgraph form.

Holst, Terry L.↗

Structure and Dynamics of the Solar Corona

Advanced computational techniques were used to study solar coronal heating and coronal mass ejections. A three dimensional, time dependent resistive magnetohydrodynamic code was used to study the dynamic response of a model corona to continuous, slow, random magnetic footpoint displacements in the photosphere. Three dimensional numerical simulations of the response of the corona to simple smooth braiding flows in the photosphere were calculated to illustrate and understand the spontaneous formation of current filaments. Two dimensional steady state helmet streamer configurations were obtained by determining the time asymptotic state of the interaction of an initially one dimensinal transponic solar wind with a spherical potential dipole field. The disruption of the steady state helmet streamer configuration was studied as a response to shearing of the magnetic footpoints of the closed field lines under the helmet.

Schnack, D. D.↗

The birthplace of compact groups of galaxies

We use complete redshift surveys to study the redshift neighborhoods of 38 Hickson compact groups (HCGs). Twenty-nine of these HCGs (76%) are embedded in rich looser systems which we call HCG associations. Analysis of the redshift neighborhood of HCGs outside the CfA survey suggests that most HCGs are embedded in more extended physical systems. Rich loose groups extracted from the CfA survey (Ramella et al. (1994)) have physical properties similar to those of the HCG associations. These rich loose groups often contain compact configurations. N-body experiments (Diaferio (1994)) suggest that compact configurations analogous to HCGs form continually during the collapse of rich loose groups. These observational and numerical results suggest that rich loose groups are the birthplace of HCGs.

Ramella, Massimo↗

Parallel inhomogeneity and the Alfven resonance. 1: Open field lines

In light of a recent demonstration of the general nonexistence of a singularity at the Alfven resonance in cold, ideal, linearized magnetohydrodynamics, we examine the effect of a small density gradient parallel to uniform, open ambient magnetic field lines. To lowest order, energy deposition is quantitatively unaffected but occurs continuously over a thickened layer. This effect is illustrated in a numerical analysis of a plasma sheet boundary layer model with perfectly absorbing boundary conditions. Consequences of the results are discussed, both for the open field line approximation and for the ensuing closed field line analysis.

Hansen, P. J.↗

Regularization for Atmospheric Temperature Retrieval Problems

Passive remote sensing of the atmosphere is used to determine the atmospheric state. A radiometer measures microwave emissions from earth's atmosphere and surface. The radiance measured by the radiometer is proportional to the brightness temperature. This brightness temperature can be used to estimate atmospheric parameters such as temperature and water vapor content. These quantities are of primary importance for different applications in meteorology, oceanography, and geophysical sciences. Depending on the range in the electromagnetic spectrum being measured by the radiometer and the atmospheric quantities to be estimated, the retrieval or inverse problem of determining atmospheric parameters from brightness temperature might be linear or nonlinear. In most applications, the retrieval problem requires the inversion of a Fredholm integral equation of the first kind making this an ill-posed problem. The numerical solution of the retrieval problem requires the transformation of the continuous problem into a discrete problem. The ill-posedness of the continuous problem translates into ill-conditioning or ill-posedness of the discrete problem. Regularization methods are used to convert the ill-posed problem into a well-posed one. In this paper, we present some results of our work in applying different regularization techniques to atmospheric temperature retrievals using brightness temperatures measured with the SSM/T-1 sensor. Simulation results are presented which show the potential of these techniques to improve temperature retrievals. In particular, no statistical assumptions are needed and the algorithms were capable of correctly estimating the temperature profile corner at the tropopause independent of the initial guess.

Velez-Reyes, Miguel↗

Metrics for Linear Kinematic Features in Sea Ice

The treatment of leads as cracks or discontinuities (see Coon et al. presentation) requires some shift in the procedure of evaluation and comparison of lead-resolving models and their validation against observations. Common metrics used to evaluate ice model skills are by and large an adaptation of a least square "metric" adopted from operational numerical weather prediction data assimilation systems and are most appropriate for continuous fields and Eilerian systems where the observations and predictions are commensurate. However, this class of metrics suffers from some flaws in areas of sharp gradients and discontinuities (e.g., leads) and when Lagrangian treatments are more natural. After a brief review of these metrics and their performance in areas of sharp gradients, we present two new metrics specifically designed to measure model accuracy in representing linear features (e.g., leads). The indices developed circumvent the requirement that both the observations and model variables be commensurate (i.e., measured with the same units) by considering the frequencies of the features of interest/importance. We illustrate the metrics by scoring several hypothetical "simulated" discontinuity fields against the lead interpreted from RGPS observations.

Levy, G.↗

Response of the Earth's Magnetosphere to Changes in the Solar Wind

The solar wind couples to the magnetosphere via dynamic pressure and electric field. Pressure establishes the size and shape of the system, while the electric field transfers energy, mass, and momentum to the magnetosphere. When the interplanetary magnetic field (IMF) is antiparallel to the dayside magnetic field, magnetic reconnection connects the IMF to the dipole field. Solar wind transport of the newly opened field lines to the nightside creates an internal convection system. These open field lines must ultimately be closed by reconnection on the nightside. For many decades, it was thought that a magnetospheric substorm was the process for accomplishing this and that all magnetic activity was a consequence of substorms. It is now recognized that there are a variety of modes of response of the magnetosphere to the solar wind. In this paper, we briefly describe these modes and the conditions under which they occur. They include substorms, pseudo-breakups, poleward boundary intensifications (PBI), steady magnetospheric convection (SMC), sawtooth injection events, magnetic storms, high-intensity long-duration continuous AE activities (HILDCAAs), and storm-time activations. There are numerous explanations for these different phenomena, some of which do not involve magnetic reconnection. However, we speculate that it is possible to interpret each mode in terms of differences in the way magnetic reconnection occurs on the nightside.

McPherron, Robert L.↗

Shuttle Propulsion Overview - The Design Challenges

The major elements of the Space Shuttle Main Propulsion System include two reusable solid rocket motors integrated into recoverable solid rocket boosters, an expendable external fuel and oxidizer tank, and three reusable Space Shuttle Main Engines. Both the solid rocket motors and space shuttle main engines ignite prior to liftoff, with the solid rocket boosters separating about two minutes into flight. The external tank separates, about eight and a half minutes into the flight, after main engine shutdown and is safely expended in the ocean. The SSME's, integrated into the Space Shuttle Orbiter aft structure, are reused after post landing inspections. The configuration is called a stage and a half as all the propulsion elements are active during the boost phase, with only the SSME s continuing operation to achieve orbital velocity. Design and performance challenges were numerous, beginning with development work in the 1970's. The solid rocket motors were large, and this technology had never been used for human space flight. The SSME s were both reusable and very high performance staged combustion cycle engines, also unique to the Space Shuttle. The multi body side mount configuration was unique and posed numerous integration and interface challenges across the elements. Operation of the system was complex and time consuming. This paper describes the design challenges and key areas where the design evolved during the program.

Owen, James W.↗

Numerical Simulation of High-Speed Turbulent Reacting Flows

The purpose of this research is to continue our efforts in advancing the state of knowledge in large eddy simulation (LES) methods for computational analysis of high-speed reacting turbulent flows. We have just completed the first year of Phase 3 of this research.

Jaberi, F. A.↗

Turbulence models from the point of view of an industrial user

Since there appears to be little prospect of the development of methods for the solution of the full, time dependent Navier-Stokes equations in the near future, the need continues for turbulence models to approximate the Reynolds stress terms that appear in the time averaged Navier-Stokes equations. Even if methods were available for solving the full equations, this would not necessarily be the optimum choice in all cases. As the cost of numerical computations decreases, the trend toward the use of more complex methods is likely to continue, but there will always be a need for a range of methods, depending on the accuracy and detail required from the calculation.

Birch, S. F.↗