Search NASA⌕ Search

SEARCH · Search NASA

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 523 records · Page 29

Non-standard analysis and embedded software

One model for computing in the future is ubiquitous, embedded computational devices analogous to embedded electrical motors. Many of these computers will control physical objects and processes. Such hidden computerized environments introduce new safety and correctness concerns whose treatment go beyond present Formal Methods. In particular, one has to begin to speak about Real Space software in analogy with Real Time software. By this we mean, computerized systems which have to meet requirements expressed in the real geometry of space. How to translate such requirements into ordinary software specifications and how to carry out proofs is a major challenge. In this talk we propose a research program based on the use of no-standard analysis. Much detail remains to be carried out. The purpose of the talk is to inform the Formal Methods community that Non-Standard Analysis provides a possible avenue to attack which we believe will be fruitful.

Platek, Richard↗

Visualizing turbulent mixing of gases and particles

A physical model and interactive computer graphics techniques have been developed for the visualization of the basic physical process of stochastic dispersion and mixing from steady-state CFD calculations. The mixing of massless particles and inertial particles is visualized by transforming the vector field from a traditionally Eulerian reference frame into a Lagrangian reference frame. Groups of particles are traced through the vector field for the mean path as well as their statistical dispersion about the mean position by using added scalar information about the root mean square value of the vector field and its Lagrangian time scale. In this way, clouds of particles in a turbulent environment are traced, not just mean paths. In combustion simulations of many industrial processes, good mixing is required to achieve a sufficient degree of combustion efficiency. The ability to visualize this multiphase mixing can not only help identify poor mixing but also explain the mechanism for poor mixing. The information gained from the visualization can be used to improve the overall combustion efficiency in utility boilers or propulsion devices. We have used this technique to visualize steady-state simulations of the combustion performance in several furnace designs.

Ma, Kwan-Liu↗

The NASA Lewis Research Center: An Economic Impact Study

The NASA Lewis Research Center (LeRC), established in 1941, is one of ten NASA research centers in the country. It is situated on 350 acres of land in Cuyahoga County and occupies more than 140 buildings and over 500 specialized research and test facilities. Most of LeRC's facilities are located in the City of Cleveland; some are located within the boundaries of the cities of Fairview Park and Brookpark. LeRC is a lead center for NASA's research, technology, and development in the areas of aeropropulsion and selected space applications. It is a center of excellence for turbomachinery, microgravity fluid and combustion research, and commercial communication. The base research and technology disciplines which serve both aeronautics and space areas include materials and structures, instrumentation and controls, fluid physics, electronics, and computational fluid dynamics. This study investigates LeRC's economic impact on Northeast Ohio's economy. It was conducted by The Urban Center's Economic Development Program in Cleveland State University's Levin College of Urban Affairs. The study measures LeRC's direct impact on the local economy in terms of jobs, output, payroll, and taxes, as well as the indirect impact of these economic activities when they 'ripple' throughout the economy. To fully explain LeRC's overall impact on the region, its contributions in the areas of technology transfer and education are also examined. The study uses a highly credible and widely accepted research methodology. First, regional economic multipliers based on input-output models were used to estimate the effect of LERC spending on the Northeast Ohio economy. Second, the economic models were complemented by interviews with industrial, civic, and university leaders to qualitatively assess LeRC's impact in the areas of technology transfer and education.

Austrian, Ziona↗

Computational Modeling Method for Superalloys

Computer modeling based on theoretical quantum techniques has been largely inefficient due to limitations on the methods or the computer needs associated with such calculations, thus perpetuating the notion that little help can be expected from computer simulations for the atomistic design of new materials. In a major effort to overcome these limitations and to provide a tool for efficiently assisting in the development of new alloys, we developed the BFS method for alloys, which together with the experimental results from previous and current research that validate its use for large-scale simulations, provide the ideal grounds for developing a computationally economical and physically sound procedure for supplementing the experimental work at great cost and time savings.

Bozzolo, Guillermo↗

Coseismic Excitation of the Earth's Polar Motion

Apart from the "shaking" near the epicenter that is the earthquake, a seismic event creates a permanent field of dislocation in the entire Earth. This redistribution of mass changes (slightly) the Earth's inertia tensor; and the Earth's rotation will change in accordance with the conservation of angular momentum. Similar to this seismic excitation of Earth rotation variations, the same mass redistribution causes (slight) changes in the Earth's gravitational field expressible in terms of changes in the Stokes coefficients of its harmonic expansion. In this paper, we give a historical background of the subject and discuss the related physics; we then compute the geodynamic effects caused by earthquakes based on normal-mode summation scheme. The effects are computed using the centroid moment tensor (CMT) solutions for 15,814 major earthquakes from Jan., 1977, through Feb., 1999, as provided in the Harvard CMT catalog. The computational results further strengthens these findings and conclusions: (i) the strong tendency for earthquakes to make the Earth rounder and more compact (however slightly) continues; (ii) so does the trend in the seismic "nudging" of the rotation pole toward the general direction of approx. 140 E, roughly opposite to that of the observed polar drift, but two orders of magnitude smaller in drift speed.

Chao, B. F.↗

2D Quantum Transport Modeling in Nanoscale MOSFETs

We have developed physical approximations and computer code capable of realistically simulating 2-D nanoscale transistors, using the non-equilibrium Green's function (NEGF) method. This is the most accurate full quantum model yet applied to 2-D device simulation. Open boundary conditions, oxide tunneling and phase-breaking scattering are treated on an equal footing. Electron bandstructure is treated within the anisotropic effective mass approximation. We present the results of our simulations of MIT 25 and 90 nm "well-tempered" MOSFETs and compare them to those of classical and quantum corrected models. The important feature of quantum model is smaller slope of Id-Vg curve and consequently higher threshold voltage. These results are consistent with 1D Schroedinger-Poisson calculations. The effect of gate length on gate-oxide leakage and subthreshold current has been studied. The shorter gate length device has an order of magnitude smaller leakage current than the longer gate length device without a significant trade-off in on-current.

Svizhenko, Alexei↗

Supporting Scientific Analysis within Collaborative Problem Solving Environments

Collaborative problem solving environments for scientists should contain the analysis tools the scientists require in addition to the remote collaboration tools used for general communication. Unfortunately, most scientific analysis tools have been designed for a "stand-alone mode" and cannot be easily modified to work well in a collaborative environment. This paper addresses the questions, "What features are desired in a scientific analysis tool contained within a collaborative environment?", "What are the tool design criteria needed to provide these features?", and "What support is required from the architecture to support these design criteria?." First, the features of scientific analysis tools that are important for effective analysis in collaborative environments are listed. Next, several design criteria for developing analysis tools that will provide these features are presented. Then requirements for the architecture to support these design criteria are listed. Sonic proposed architectures for collaborative problem solving environments are reviewed and their capabilities to support the specified design criteria are discussed. A deficiency in the most popular architecture for remote application sharing, the ITU T. 120 architecture, prevents it from supporting highly interactive, dynamic, high resolution graphics. To illustrate that the specified design criteria can provide a highly effective analysis tool within a collaborative problem solving environment, a scientific analysis tool that contains the specified design criteria has been integrated into a collaborative environment and tested for effectiveness. The tests were conducted in collaborations between remote sites in the US and between remote sites on different continents. The tests showed that the tool (a tool for the visual analysis of computer simulations of physics) was highly effective for both synchronous and asynchronous collaborative analyses. The important features provided by the tool (and made possible by the specified design criteria) are: 1. The tool provides highly interactive, dynamic, high resolution, 3D graphics. 2. All remote scientists can view the same dynamic, high resolution, 3D scenes of the analysis as the analysis is being conducted. 3. The responsiveness of the tool is nearly identical to the responsiveness of the tool in a stand-alone mode. 4. The scientists can transfer control of the analysis between themselves. 5. Any analysis session or segment of an analysis session, whether done individually or collaboratively, can be recorded and posted on the Web for other scientists or students to download and play in either a collaborative or individual mode. 6. The scientist or student who downloaded the session can, individually or collaboratively, modify or extend the session with his/her own "what if" analysis of the data and post his/her version of the analysis back onto the Web. 7. The peak network bandwidth used in the collaborative sessions is only 1K bit/second even though the scientists at all sites are viewing high resolution (1280 x 1024 pixels), dynamic, 3D scenes of the analysis. The links between the specified design criteria and these performance features are presented.

Watson, Velvin R.↗

The Turbulent/Non-Turbulent Interface Bounding a Far-Wake

The velocity fields of a turbulent wake behind a flat plate obtained from the direct numerical simulations of Moser et al. are used to study the structure of the flow in the intermittent zone where there are, alternately, regions of fully turbulent flow and non-turbulent velocity fluctuations either side of a thin randomly moving interface. Comparisons are made with a wake that is 'forced' by amplifying initial velocity fluctuations. There is also a random temperature field T in the flow; T varies between constant values of 0.0 and 1.0 on the sides of the wake. The value of the Reynolds number based on the centreplane mean velocity defect and halfwidth b of the wake is Re approx. = 2000. It is found that the thickness of the continuous interface is about equal to 0.07b, whereas the amplitude of fluctuations of the instantaneous interface displacement y(sub I)(t) is an order of magnitude larger, being about 0.5b. This explains why the mean statistics of vorticity in the intermittent zone can be calculated in terms of the probability distribution of y(sub I) and the instantaneous discontinuity in vorticity across the interface. When plotted as functions of y - y(sub I), the conditional mean velocity (U) and temperature (T) profiles show sharp jumps Delta(U) and Delta(T) at the interface adjacent to a thick zone where (U) and (T) vary much more slowly. Statistics for the vorticity and velocity variances, available in such detail only from DNS data, show how streamwise and spanwise components of vorticity are generated by vortex stretching in the bulges of the interface. Flow fields around the interface, analyzed in terms of the local streamline pattern, confirm previous results that the advancement of the vortical interface into the irrotational flow is driven by large-scale eddy motion. It is argued that because this is an inviscid mechanism the entrainment process is not sensitive to the value of Re, and that small-scale nibbling only plays a subsidiary role. While mean Reynolds stresses decrease gradually in the intermittent zone, conditional stresses are found to decrease sharply towards zero at the interface. Using one-point turbulence models applied to either unconditional or conditional statistics for the turbulent region and then averaged, the entrainment rate E(sub b) would, if calculated exactly, be zero. But if computed with standard computational methods, E(sub b) would be non-zero because of numerical diffusion. It is concluded that the current practice in statistical models of approximating entrainment by a diffusion process is computationally arbitrary and physically incorrect. An analysis shows how E(sub b) is related to Delta(U) and the jump in shear stress at the interface, and correspondingly to Delta(T) and the heat flux.

Bisset, David K.↗

2D Quantum Mechanical Study of Nanoscale MOSFETs

With the onset of quantum confinement in the inversion layer in nanoscale MOSFETs, behavior of the resonant level inevitably determines all device characteristics. While most classical device simulators take quantization into account in some simplified manner, the important details of electrostatics are missing. Our work addresses this shortcoming and provides: (a) a framework to quantitatively explore device physics issues such as the source-drain and gate leakage currents, DIBL, and threshold voltage shift due to quantization, and b) a means of benchmarking quantum corrections to semiclassical models (such as density-gradient and quantum-corrected MEDICI). We have developed physical approximations and computer code capable of realistically simulating 2-D nanoscale transistors, using the non-equilibrium Green's function (NEGF) method. This is the most accurate full quantum model yet applied to 2-D device simulation. Open boundary conditions and oxide tunneling are treated on an equal footing. Electrons in the ellipsoids of the conduction band are treated within the anisotropic effective mass approximation. We present the results of our simulations of MIT 25, 50 and 90 nm "well-tempered" MOSFETs and compare them to those of classical and quantum corrected models. The important feature of quantum model is smaller slope of Id-Vg curve and consequently higher threshold voltage. Surprisingly, the self-consistent potential profile shows lower injection barrier in the channel in quantum case. These results are qualitatively consistent with ID Schroedinger-Poisson calculations. The effect of gate length on gate-oxide leakage and subthreshold current has been studied. The shorter gate length device has an order of magnitude smaller current at zero gate bias than the longer gate length device without a significant trade-off in on-current. This should be a device design consideration.

Svizhenko, Alexei↗

2D Quantum Transport Modeling in Nanoscale MOSFETs

With the onset of quantum confinement in the inversion layer in nanoscale MOSFETs, behavior of the resonant level inevitably determines all device characteristics. While most classical device simulators take quantization into account in some simplified manner, the important details of electrostatics are missing. Our work addresses this shortcoming and provides: (a) a framework to quantitatively explore device physics issues such as the source-drain and gate leakage currents, DIBL, and threshold voltage shift due to quantization, and b) a means of benchmarking quantum corrections to semiclassical models (such as density- gradient and quantum-corrected MEDICI). We have developed physical approximations and computer code capable of realistically simulating 2-D nanoscale transistors, using the non-equilibrium Green's function (NEGF) method. This is the most accurate full quantum model yet applied to 2-D device simulation. Open boundary conditions, oxide tunneling and phase-breaking scattering are treated on equal footing. Electrons in the ellipsoids of the conduction band are treated within the anisotropic effective mass approximation. Quantum simulations are focused on MIT 25, 50 and 90 nm "well- tempered" MOSFETs and compared to classical and quantum corrected models. The important feature of quantum model is smaller slope of Id-Vg curve and consequently higher threshold voltage. These results are quantitatively consistent with I D Schroedinger-Poisson calculations. The effect of gate length on gate-oxide leakage and sub-threshold current has been studied. The shorter gate length device has an order of magnitude smaller current at zero gate bias than the longer gate length device without a significant trade-off in on-current. This should be a device design consideration.

Svizhenko, Alexei↗

2-D Modeling of Nanoscale MOSFETs: Non-Equilibrium Green's Function Approach

We have developed physical approximations and computer code capable of realistically simulating 2-D nanoscale transistors, using the non-equilibrium Green's function (NEGF) method. This is the most accurate full quantum model yet applied to 2-D device simulation. Open boundary conditions and oxide tunneling are treated on an equal footing. Electrons in the ellipsoids of the conduction band are treated within the anisotropic effective mass approximation. Electron-electron interaction is treated within Hartree approximation by solving NEGF and Poisson equations self-consistently. For the calculations presented here, parallelization is performed by distributing the solution of NEGF equations to various processors, energy wise. We present simulation of the "benchmark" MIT 25nm and 90nm MOSFETs and compare our results to those from the drift-diffusion simulator and the quantum-corrected results available. In the 25nm MOSFET, the channel length is less than ten times the electron wavelength, and the electron scattering time is comparable to its transit time. Our main results are: (1) Simulated drain subthreshold current characteristics are shown, where the potential profiles are calculated self-consistently by the corresponding simulation methods. The current predicted by our quantum simulation has smaller subthreshold slope of the Vg dependence which results in higher threshold voltage. (2) When gate oxide thickness is less than 2 nm, gate oxide leakage is a primary factor which determines off-current of a MOSFET (3) Using our 2-D NEGF simulator, we found several ways to drastically decrease oxide leakage current without compromising drive current. (4) Quantum mechanically calculated electron density is much smaller than the background doping density in the poly silicon gate region near oxide interface. This creates an additional effective gate voltage. Different ways to. include this effect approximately will be discussed.

Svizhenko, Alexei↗

Fusion for Space Propulsion

The need for fusion propulsion for interplanetary flights is discussed. For a propulsion system, there are three important system attributes: (1) The absolute amount of energy available, (2) the propellant exhaust velocity, and (3) the jet power per unit mass of the propulsion system (specific power). For efficient and affordable human exploration of the solar system, propellant exhaust velocity in excess of 100 km/s and specific power in excess of 10 kW/kg are required. Chemical combustion obviously cannot meet the requirement in propellant exhaust velocity. Nuclear fission processes typically result in producing energy in the form of heat that needs to be manipulated at temperatures limited by materials to about 2,800 K. Using the fission energy to heat a low atomic weight propellant produces propellant velocity of the order of 10 kinds. Alternatively the fission energy can be converted into electricity that is used to accelerate particles to high exhaust velocity. However, the necessary power conversion and conditioning equipment greatly increases the mass of the propulsion system. Fundamental considerations in waste heat rejection and power conditioning in a fission electric propulsion system place a limit on its jet specific power to the order of about 0.2 kW/kg. If fusion can be developed for propulsion, it appears to have the best of all worlds - it can provide the largest absolute amount of energy, the propellant exhaust velocity (> 100 km/s), and the high specific jet power (> 10 kW/kg). An intermediate step towards fusion propulsion might be a bimodal system in which a fission reactor is used to provide some of the energy to drive a fusion propulsion unit. There are similarities as well as differences between applying fusion to propulsion and to terrestrial electrical power generation. The similarities are the underlying plasma and fusion physics, the enabling component technologies, the computational and the diagnostics capabilities. These physics and engineering capabilities have been demonstrated for a fusion reactor gain (Q) of the order of unity (TFTR: 0.25, JET: 0.65, JT-60: Q(sub eq) approx. 1.25). These technological advances made it compelling for considering fusion for propulsion.

Thio, Y. C. Francis↗

Constitutive Relation in Transitional Granular Flows

To study the constitutive behavior of granular materials, the presence of gravity is detrimental. Although empirical relations have been obtained for engineering designs to control granular flows on Earth, it is not known how well these Earthbound relations can be used in another gravity field. Fundamental understanding must be derived to reliably design for granular flows in space exploration. There are two extremes of granular flows of which significant amount of knowledge is available. One deals with a dense and quasi-static situation where the deformation rate nearly vanishes. The other deals with dilute and rapidly fluctuating grain velocities where particle inertia dominates. This project, funded by the NASA Microgravity Fluid Physics Program, aims to study this transitional regime via physical experiments and computer simulations. A conceptual model has been established as described below. There are two natural time scales in a granular flow. One is the travel time between two consecutive collisions and the other is the duration of a collision contact. At a very low shear-rate, the shear-induced particle velocity is low. Hence the travel time between collisions is longer than the contact time between colliding particles. Binary collisions prevail. As the shear-rate increases, the traveling time between collisions reduces and the probability of multiple collisions goes up. These particle groups disperse shortly after and new groups form. When shear-rate is further increased, clusters grow in size due to an increasing chance for free particles to join before groups have the time to disperse. The maximum cluster size may depend on the global concentration and material properties. As the solid concentration approaches zero, the cluster size goes to one particle diameter. The maximum possible cluster size under any condition is the container size, provided that the shear flow is inside a container. The critical shear-rates that dictate the initiation of the multiple contacts, and the size and lifetime of the collision clusters, are functions of the concentration also. A 'regime' theory has been proposed by Babic et al. This theory suggested that both the solid concentration, C, and the non-dimensional shear-rate, B, are important in determining the regimes of the granular constitutive law. Additional information is included in the original extended abstract.

Shen, Hayley H.↗

Classification of Tropical Oceanic Precipitation Using High Altitude Aircraft Microwave and Electric Field Measurements

A physically intuitive and computationally simple precipitation mapping algorithm has been developed for use with the airborne Advanced Microwave Precipitation Radiometer (AMPR). The algorithm is based on microwave emission and scattering properties of precipitation. Specifically, emission by liquid water allows increasing brightness temperatures at low frequencies to be interpreted as increasing rain rates. Scattering by large hydrometeors (particularly graupel and hail) causes relative minima in the brightness temperatures, with progressively larger hydrometeors scattering progressively longer wavelengths. The vigor of convection is therefore ascertained according to which wavelengths are being significantly scattered. The combination of emission and scattering information from four microwave channels is used to assign a precipitation category, which is related to the liquid rain rate, the vertical extent of precipitation, and the vigor of convection. The qualitative precipitation categories output by the passive microwave algorithm have been verified using coincident radar (ER-2 Doppler Radar - EDOP) and electric field measurements (Lightning Instrument Package - LIP). These coincident measurements can subsequently be used to quantify rain rates, hydrometeor contents, and vertical profiles that are typical for each precipitation category. This algorithm has been developed using an airborne platform. Comparisons are being made with other airborne, satellite, and ground-based radar and radiometer data. This technique shows promise both as a research tool and potentially as a real-time analysis tool, which could be applied to either traditional or uninhabited aerial vehicles.

Cecil, Daniel J.↗

Classification of Tropical Oceanic Precipitation Using High Altitude Aircraft Microwave and Electric Field Measurements

A physically intuitive and computationally simple precipitation mapping algorithm has been developed for use with the airborne Advanced Microwave Precipitation Radiometer (AMPR). The algorithm is based on microwave emission and scattering properties of precipitation. Specifically, emission by liquid water allows increasing brightness temperatures at low frequencies to be interpreted as increasing rain rates. Scattering by large hydrometeors (particularly graupel and hail) causes relative minima in the brightness temperatures, with progressively larger hydrometeors scattering progressively longer wavelengths. The vigor of convection is therefore ascertained according to which wavelengths are being significantly scattered. The combination of emission and scattering information from four microwave channels is used to assign a precipitation category, which is related to the liquid rain rate, the vertical extent of precipitation, and the vigor of convection. The qualitative precipitation categories output by the passive microwave algorithm have been verified using coincident radar (ER-2 Doppler Radar - EDOP) and electric field measurements (Lightning Instrument Package - LIP). These coincident measurements can subsequently be used to quantify rain rates, hydrometeor contents, and vertical profiles that are typical for each precipitation category. This algorithm has been developed using an airborne platform. Comparisons are being made with other airborne, satellite, and ground-based radar and radiometer data. This technique shows promise both as a research tool and potentially as a real-time analysis tool, which could be applied to either traditional or uninhabited aerial vehicles.

Cecil, Daniel J.↗

Modeling Code Is Helping Cleveland Develop New Products

Master Builders, Inc., is a 350-person company in Cleveland, Ohio, that develops and markets specialty chemicals for the construction industry. Developing new products involves creating many potential samples and running numerous tests to characterize the samples' performance. Company engineers enlisted NASA's help to replace cumbersome physical testing with computer modeling of the samples' behavior. Since the NASA Lewis Research Center's Structures Division develops mathematical models and associated computation tools to analyze the deformation and failure of composite materials, its researchers began a two-phase effort to modify Lewis' Integrated Composite Analyzer (ICAN) software for Master Builders' use. Phase I has been completed, and Master Builders is pleased with the results. The company is now working to begin implementation of Phase II.

Source record↗

The world problem: on the computability of the topology of 4-manifolds

Topological classification of the 4-manifolds bridges computation theory and physics. A proof of the undecidability of the homeomorphy problem for 4-manifolds is outlined here in a clarifying way. It is shown that an arbitrary Turing machine with an arbitrary input can be encoded into the topology of a 4-manifold, such that the 4-manifold is homeomorphic to a certain other 4-manifold if and only if the corresponding Turing machine halts on the associated input. Physical implications are briefly discussed.

vanMeter, J. R.↗

The Equivalence of Precession Phenomena in Metric Theories of Gravity

The requirement of general covariance imparts to metric theories of gravity, such as general relavity, important structural features. A precise mathematical form results, ensuring that computation of observable physical effects in the theory gives the same answers independently of the chosen system of coordinates.

metric theories gravitation relativity↗