Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Computational Physics”

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 955 records · Page 53

Turbulent excitation of spontaneous reconnection

The long-term nonlinear evolution of a tearing-mode-unstable sheared-field plasma in a turbulent environment is explored. Two different physical configurations are modeled, and a different computational system is used for each. Results of both sets of calculations show that magnetic tearing arises spontaneously provided that the initial turbulence energy level is below the natural saturation level of the tearing instability. The relationship between these results and those of previous calculations are briefly discussed, concluding that there are no significant unexplainable disagreements.

Deeds, D.↗

Fault diagnosis

The objective of the research in this area of fault management is to develop and implement a decision aiding concept for diagnosing faults, especially faults which are difficult for pilots to identify, and to develop methods for presenting the diagnosis information to the flight crew in a timely and comprehensible manner. The requirements for the diagnosis concept were identified by interviewing pilots, analyzing actual incident and accident cases, and examining psychology literature on how humans perform diagnosis. The diagnosis decision aiding concept developed based on those requirements takes abnormal sensor readings as input, as identified by a fault monitor. Based on these abnormal sensor readings, the diagnosis concept identifies the cause or source of the fault and all components affected by the fault. This concept was implemented for diagnosis of aircraft propulsion and hydraulic subsystems in a computer program called Draphys (Diagnostic Reasoning About Physical Systems). Draphys is unique in two important ways. First, it uses models of both functional and physical relationships in the subsystems. Using both models enables the diagnostic reasoning to identify the fault propagation as the faulted system continues to operate, and to diagnose physical damage. Draphys also reasons about behavior of the faulted system over time, to eliminate possibilities as more information becomes available, and to update the system status as more components are affected by the fault. The crew interface research is examining display issues associated with presenting diagnosis information to the flight crew. One study examined issues for presenting system status information. One lesson learned from that study was that pilots found fault situations to be more complex if they involved multiple subsystems. Another was pilots could identify the faulted systems more quickly if the system status was presented in pictorial or text format. Another study is currently under way to examine pilot mental models of the aircraft subsystems and their use in diagnosis tasks. Future research plans include piloted simulation evaluation of the diagnosis decision aiding concepts and crew interface issues. Information is given in viewgraph form.

Abbott, Kathy↗

Deriving Earth science products from SSM/I

There are 3 components to this project. The first is the production of a quality controlled Level-1 data product for the Special Sensor Microwave Imager (SSM/I). Second is the generation of research quality ocean products. The 3rd component is studying the feasibility of obtaining both wind speed and wind direction from satellite microwave radiometers. The status of the three studies are discussed briefly. Time series of the SSM/I physical temperatures, receiver gains, and noise figures were studied in order to access the stability and integrity of the sensor data. A physically based retrieval algorithm was used to compute ocean products from the SSM/I Level-1 data. It simultaneously finds the near surface wind speed, the columnar water vapor, and the columnar liquid water. A wind direction signal was found in the SSM/I brightness temperature.

Wentz, Frank J.↗

Vortex ring bubbles

Toroidal bubbles with circulation are studied numerically and by means of a physically motivated model equation. Two series of computations are performed by a boundary-integral method. One set shows the starting motion of an initially spherical bubble as a gravitationally driven liquid jet penetrates through the bubble from below causing a toroidal geometry to develop. The jet becomes broader as surface tension increases and fails to penetrate if surface tension is too large. The dimensionless circulation that develops is not very dependent on the surface tension. The second series of computations starts from a toroidal geometry, with circulation determined from the earlier series, and follows the motion of the rising and spreading vortex ring. Some modifications to the boundary-integral formulation were devised to handle the multiply connected geometry. The computations uncovered some unexpected rapid oscillations of the ring radius. These oscillations and the spreading of the ring are explained by the model equation which provides a more general description of vortex ring bubbles than previously available.

Lundgren, T. S.↗

Mass and energy flow in prominences

Mass and energy flow in quiescent prominences is considered based on the hypothesis that active region prominences have a different structure and thus different mass and energy flow characteristics. Several important physical parameters have been plotted using the computational model, representing the evolutionary process after the prominence formation. The temperature, velocity, conductive flux, and enthalpy flux are plotted against distance from the highest point in the loop to the coolest part of the prominence. It is shown that the maximum velocity is only about 5 km/s. The model calculations indicate that the transition region of prominences is dominated by complex processes. It is necessary to take into account mass flow at temperatures below 200,000 K, and both mass flow and optical depth effects in hydrogen at temperatures below 30,000 K. Both of these effects lead to a less steep temperature gradient through the prominence corona interface than can be obtained from the conduction alone.

Poland, Arthur I.↗

Tomography, the geoid and plate motions

The dynamics of the earth's mantle and its relationship to mantle structure as revealed by seismic tomography, the geoid, and plate motions are discussed in a critical review of U.S. research from the period 1987-1990. Sections are devoted to plates and mantle convection; seismic anomalies, flow, and the geoid; subducted slabs and the 670-km discontinuity; the physical properties of the mantle; plumes; and computational and experimental modeling efforts. Also included is a comprehensive bibliography for the period.

Gable, Carl W.↗

Extension of convection modeling into the high-latitude ionosphere - Some theoretical difficulties

The Rice Convection Model (RCM) is extended and merged with empirical models so as to cover the entire high-latitude ionosphere with the aim of providing precipitation and electric field inputs for ionosphere and thermosphere modelers and producing a model in which the boundaries of the precipitation and electric field patterns maintain physically consistent relationships to each other. The computed auroral electron energy flux, plotted as a function of latitude, exhibited an exaggerated two-peak structure. When no floor was placed under the precipitation rate, the minimum between the two peaks was much too deep to be consistent with typical observations. The regions of excessively weak precipitation map to equatorial distances of 15-35 RE and thus to the regions of the plasma sheet that were not included in previous self-consistent convection calculations.

Wolf, R. A.↗

The OH distribution in cometary atmospheres - A collisional Monte Carlo model for heavy species

The study presents an extension of the cometary atmosphere Monte Carlo particle trajectory model formalism which makes it both physically correct for heavy species and yet computationally reasonable. The derivation accounts for the collision path and scattering redirection of a heavy radical traveling through a fluid coma with a given radial distribution in outflow speed and temperature. The revised model verifies that the earlier fast-H atom approximations used in earlier work are valid, and it is applied to a case where the heavy radical formalism is necessary: the OH distribution. It is found that a steeper variation of water production rate with heliocentric distance is required for a water coma which is consistent with the velocity-resolved observations of Comet P/Halley.

Combi, Michael R.↗

FORTH direct execution processors in the Hopkins ultraviolet telescope

Engineers at the Johns Hopkins University Applied Physics Laboratory have designed and built two computers which will control and monitor the Hopkins Ultraviolet Telescope, a 1986 Space Shuttle experiment. These computers implement a microcoded FORTH nucleus in a word-addressed AMD 2900 series bit-slice architecture. All programs for them are written in FORTH, which takes the place of assembly language in this architecture.

Ballard, Ben↗

Development of deterministic transport methods for low energy neutrons for shielding in space

Transport of low energy neutrons associated with the galactic cosmic ray cascade is analyzed in this dissertation. A benchmark quality analytical algorithm is demonstrated for use with BRYNTRN, a computer program written by the High Energy Physics Division of NASA Langley Research Center, which is used to design and analyze shielding against the radiation created by the cascade. BRYNTRN uses numerical methods to solve the integral transport equations for baryons with the straight-ahead approximation, and numerical and empirical methods to generate the interaction probabilities. The straight-ahead approximation is adequate for charged particles, but not for neutrons. As NASA Langley improves BRYNTRN to include low energy neutrons, a benchmark quality solution is needed for comparison. The neutron transport algorithm demonstrated in this dissertation uses the closed-form Green's function solution to the galactic cosmic ray cascade transport equations to generate a source of neutrons. A basis function expansion for finite heterogeneous and semi-infinite homogeneous slabs with multiple energy groups and isotropic scattering is used to generate neutron fluxes resulting from the cascade. This method, called the FN method, is used to solve the neutral particle linear Boltzmann transport equation. As a demonstration of the algorithm coded in the programs MGSLAB and MGSEMI, neutron and ion fluxes are shown for a beam of fluorine ions at 1000 MeV per nucleon incident on semi-infinite and finite aluminum slabs. Also, to demonstrate that the shielding effectiveness against the radiation from the galactic cosmic ray cascade is not directly proportional to shield thickness, a graph of transmitted total neutron scalar flux versus slab thickness is shown. A simple model based on the nuclear liquid drop assumption is used to generate cross sections for the galactic cosmic ray cascade. The ENDF/B V database is used to generate the total and scattering cross sections for neutrons in aluminum. As an external verification, the results from MGSLAB and MGSEMI were compared to ANISN/PC, a routinely used neutron transport code, showing excellent agreement. In an application to an aluminum shield, the FN method seems to generate reasonable results.

Ganapol, Barry↗

Advanced interactive display formats for terminal area traffic control

This report describes the basic design considerations for perspective air traffic control displays. A software framework has been developed for manual viewing parameter setting (MVPS) in preparation for continued, ongoing developments on automated viewing parameter setting (AVPS) schemes. Two distinct modes of MVPS operations are considered, both of which utilize manipulation pointers imbedded in the three-dimensional scene: (1) direct manipulation of the viewing parameters -- in this mode the manipulation pointers act like the control-input device, through which the viewing parameter changes are made. Part of the parameters are rate controlled, and part of them position controlled. This mode is intended for making fast, iterative small changes in the parameters. (2) Indirect manipulation of the viewing parameters -- this mode is intended primarily for introducing large, predetermined changes in the parameters. Requests for changes in viewing parameter setting are entered manually by the operator by moving viewing parameter manipulation pointers on the screen. The motion of these pointers, which are an integral part of the 3-D scene, is limited to the boundaries of the screen. This arrangement has been chosen in order to preserve the correspondence between the spatial lay-outs of the new and the old viewing parameter setting, a feature which contributes to preventing spatial disorientation of the operator. For all viewing operations, e.g. rotation, translation and ranging, the actual change is executed automatically by the system, through gradual transitions with an exponentially damped, sinusoidal velocity profile, in this work referred to as 'slewing' motions. The slewing functions, which eliminate discontinuities in the viewing parameter changes, are designed primarily for enhancing the operator's impression that he, or she, is dealing with an actually existing physical system, rather than an abstract computer-generated scene. The proposed, continued research efforts will deal with the development of automated viewing parameter setting schemes. These schemes employ an optimization strategy, aimed at identifying the best possible vantage point, from which the air traffic control scene can be viewed for a given traffic situation. They determine whether a change in viewing parameter setting is required and determine the dynamic path along which the change to the new viewing parameter setting should take place.

Grunwald, Arthur J.↗

Simulations of vortex generators

We are interested in the study, via direct numerical simulations, of active vortex generators. Vortex generators may be used to modify the inner part of the boundary layer or to control separation thus enhancing the performance and maneuverability of aerodynamic configurations. We consider generators that consist of a surface cavity elongated in the stream direction and partially covered with a moving lid that at rest lies flush with the boundary. Streamwise vorticity is generated and ejected due to the oscillatory motion of the lid. The present simulations complement relevant experimental investigations of active vortex generators at NASA Ames and Stanford University (Saddoughi, 1994, and Jacobson and Reynolds, 1993). Jacobson and Reynolds (1993) used a piezoelectric device in water, allowing for small amplitude high frequency oscillations. They placed the lid asymmetrically on the cavity and observed a strong outward velocity at the small gap of the cavity. Saddoughi used a larger mechanically driven device in air to investigate this flow and he observed a jet emerging from the wide gap of the configuration, contrary to the findings of Jacobson and Reynolds. Our task is to simulate the flows generated by these devices and to conduct a parametric study that would help us elucidate the physical mechanisms present in the flow. Conventional computational schemes encounter difficulties when simulating flows around complex configurations undergoing arbitrary motions. Here we present a formulation that achieves this task on a purely Lagrangian frame by extending the formulation presented by Koumoutsakos, Leonard and Pepin (1994). The viscous effects are taken into account by modifying the strength of the particles, whereas fast multipole schemes employing hundreds of thousands of particles allow for high resolution simulations. The results of the present simulations would help us assess some of the effects of three-dimensionality in experiments and investigate the role of two-dimensional vortex generation due to an oscillating lid.

Koumoutsakos, P.↗

Multiwavelength Modeling of Nova Atmospheres

We have analyzed the early optically thick ultraviolet spectra of Nova OS And 1986 using a grid of spherically symmetric, non-LTE, line-blanketed, expanding model atmospheres and synthetic spectra with the following set of parameters: 5,000 less than or equal to T(sub model) less than or equal to 60,000K, solar abundances, (rho)(alpha) r(sup -3), v(sub max) = 2000 km/s, L = 6 x 10(exp 4) solar luminosity, and a statistical or microturbulent velocity of 50 km/s. We used the synthetic spectra to estimate the model parameters corresponding to the observed IUE spectra. The fits to the observations were then iteratively improved by changing the parameters of the model atmospheres, in particular T(sub model) and the abundances, to arrive at the best fits to the optically thick pseudo-continuum and the features found in the IUE spectra. The IUE spectra show two different optically thick subphases. The earliest spectra, taken a few days after maximum optical light, show a pseudo-continuum created by overlapping absorption lines. The later observations, taken approximately 3 weeks after maximum light, show the simultaneous presence of allowed, semi-forbidden, and forbidden lines in the observed spectra. Analysis of these phases indicate that OS And 86 had solar metallicities except for Mg which showed evidence of being underabundant by as much as a factor of 10. We determine a distance of 5.1 kpc to OS And 86 and derive a peak bolometric luminosity of approximately 5 x 10(exp 4) solar luminosity. The computed nova parameters provide insights into the physics of the early outburst and explain the spectra seen by IUE. Lastly, we find evidence in the later observations for large non-LTE effects of Fe II which, when included, lead to much better agreement with the observations.

Hauschildt, P. H.↗

Advanced Interactive Display formats for Terminal Area Traffic Control

This report describes the basic design considerations for perspective Air Traffic Control displays. A software framework has been developed for manual viewing parameter setting (MVPS) in preparation for continued, ongoing developments on automated viewing parameter setting (AVPS) schemes. Two distinct modes of MVPS operations are considered, both of which utilize manipulation pointers imbedded in the three-dimensional scene: (1) direct manipulation of the viewing parameters; in this mode the manipulation pointers act like the control-input device, through which the viewing parameter changes are made. Part of the parameters are rate controlled, and part of them position controlled. This mode is intended for making fast, iterative small changes in the parameters. (2) indirect manipulation of the viewing parameters. This mode is intended primarily for introducing large, predetermined changes in the parameters. Requests for changes in viewing parameter setting are entered manually by the operator by moving viewing parameter manipulation pointers on the screen. The motion of these pointers, which are an integral part of the 3-D scene, is limited to the boundaries of screen. This arrangement has been chosen, in order to preserve the correspondence between the spatial lay-outs of the new and the old viewing parameter setting, a feature which contributes to preventing spatial disorientation of the operator. For all viewing operations, e.g. rotation, translation and ranging, the actual change is executed automatically by the system, through gradual transitions with an exponentially damped, sinusoidal velocity profile, in this work referred to as 'slewing' motions. The slewing functions, which eliminate discontinuities in the viewing parameter changes, are designed primarily for enhancing the operator's impression that he, or she, is dealing with an actually existing physical system, rather than an abstract computer-generated scene, The proposed, continued research efforts will deal with the development of automated viewing parameter setting schemes. These schemes employ an optimization strategy, aimed at identifying the best possible vantage point, from which the Air Traffic Control scene can be viewed for a given traffic situation. They determine whether a change in viewing parameter setting is required and determine the dynamic path along which the change to the new viewing parameter setting should take place.

Grunwald, Arthur J.↗

Cranked Arrow Wing (F-16XL-1) Flight Flow Physics with CFD Predictions at Subsonic and Transonic Speeds

The computational fluid dynamics (CFD) modeling used has produced reasonably good global upper-surface pressure coefficient comparisons with measured flight data at both transonic and subsonic speeds at the angles of attack presented. Boundary layer comparisons showed the profiles to be reasonably well predicted inboard and under the primary vortex system. However, the secondary vortex profile was not well predicted either at the anticipated separation point or under the secondary vortex. Moreover, the flight data showed there to be a vortex/boundary-layer interaction that occurred in the vicinity of the secondary vortex. The spanwise distribution of local skin friction measured data was reasonably well predicted, especially away from the wing leading-edge. Lastly, predicted and measured flight-pressures, as well as flight-image data, for the F-16XL-1 airplane are now available via the World Wide Web.

Lamar, John E.↗

Simulations of Active Vortex Generators

We are interested in the study, via numerical simulations, of active vortex generators. Vortex generators may be used to modify the inner part of the boundary layer or to control separation thus enhancing the performance and maneuverability of aerodynamic configurations. we consider generators that consist of a surface cavity elongated in the streamwise direction and partially covered with a moving lid that at rest lies flush with the boundary. Streamwise vorticity is generated and ejected due to the oscillatory motion of the lid. The present simulations complement relevant experimental investigations of active vortex generators that have been conducted at NASA Ames Research Center and Stanford University used a piezoelectric device in water, allowing for small amplitude high frequency oscillations. They placed the lid asymmetrically on the cavity and observed a strong outward velocity at the small gap of the cavity. Saddoughi used a larger mechanically driven device in air to investigate this flow and observed a jet emerging from the wide gap of the configuration, contrary to the findings of Jacobson and Reynolds We are simulating the flows generated by these devices and we are conducting a parametric study that would help us elucidate the physical mechanisms present in the flow. Conventional computational schemes encounter difficulties when simulating flows around complex configurations undergoing arbitrary motions. Here we present a formulation that achieves this task on a purely Lagrangian frame by extending the formulation presented by Koumoutsakos, Leonard and Pepin (1994). The viscous effects are taken into account by modifying the strength Of the particles, whereas fast multipole schemes employing hundreds of thousands of particles allow for high resolution simulations. We shall present simulation results of an oscillating plate at various Reynolds numbers and Strouhal frequencies.

Mansour, N. N.↗

Simulations of Active Vortex Generators

We are interested in the study, via numerical simulations, of active vortex generators. Vortex generators may be used to modify the inner part of the boundary layer or to control separation thus enhancing the performance and maneuverability of aerodynamic configurations. We consider generators that consist of a surface cavity elongated in the streamwise direction and partially covered with a moving lid that at rest lies flush with the boundary. Streamwise voracity is generated and ejected due to the oscillatory motion of the lid. The present simulations c Implement relevant experimental investigations of active vortex generators that have been conducted at NASA Ames Research Center and Stanford University. Jacobson and Reynolds used a piezoelectric device in water, allowing for small amplitude high frequency oscillations. They placed the lid asymmetrically on the cavity and observed a strong outward velocity at the small gap of the cavity. Saddoughi used a larger mechanically driven device in air to investigate this flow and observed a jet emerging from the wide gap of the configuration, contrary to the findings of Jacobson and Reynolds. More recently, Lachowiez and Wlezien are investigating the flow generated by an electro-mechanically driven lid to be used for assertion control in aerodynamic applications. We are simulating the flows generated by these devices and we are conducting a parametric study that would help us elucidate the physical mechanisms present in the flow. Conventional computational schemes encounter difficulties when simulating flows around complex configurations undergoing arbitrary motions. Here we present a formulation that achieves this task on a purely Lagrangian frame by extending the formulation presented by Koumoutsakos, Leonard and Pepin. The viscous effects are taken into account by modifying the strength of the particles, whereas fast multipole schemes employing hundreds of thousands ol'particle's allow for high resolution simulations. We shall present simulation results of an oscillating plate at various Reynolds numbers and Strouhal frequencies. Estimates of the forces needed to drive the devices will also be presented.

Mansour, N. N.↗

A Response Surface Methodology for Bi-Level Integrated System Synthesis (BLISS)

The report describes a new method for optimization of engineering systems such as aerospace vehicles whose design must harmonize a number of subsystems and various physical phenomena, each represented by a separate computer code, e.g., aerodynamics, structures, propulsion, performance, etc. To represent the system internal couplings, the codes receive output from other codes as part of their inputs. The system analysis and optimization task is decomposed into subtasks that can be executed concurrently, each subtask conducted using local state and design variables and holding constant a set of the system-level design variables. The subtasks results are stored in form of the Response Surfaces (RS) fitted in the space of the system-level variables to be used as the subtask surrogates in a system-level optimization whose purpose is to optimize the system objective(s) and to reconcile the system internal couplings. By virtue of decomposition and execution concurrency, the method enables a broad workfront in organization of an engineering project involving a number of specialty groups that might be geographically dispersed, and it exploits the contemporary computing technology of massively concurrent and distributed processing. The report includes a demonstration test case of supersonic business jet design.

Altus, Troy David↗