Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Direct numerical simulation”

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 1,081 records · Page 60

Identifying Bearing Rotordynamic Coefficients using an Extended Kalman Filter

An Extended Kalman Filter is developed to estimate the linearized direct and indirect stiffness and damping force coefficients for bearings in rotor-dynamic applications from noisy measurements of the shaft displacement in response to imbalance and impact excitation. The bearing properties are modeled as stochastic random variables using a Gauss-Markov model. Noise terms are introduced into the system model to account for all of the estimation error, including modeling errors and uncertainties and the propagation of measurement errors into the parameter estimates. The system model contains two user-defined parameters that can be tuned to improve the filter s performance; these parameters correspond to the covariance of the system and measurement noise variables. The filter is also strongly influenced by the initial values of the states and the error covariance matrix. The filter is demonstrated using numerically simulated data for a rotor-bearing system with two identical bearings, which reduces the number of unknown linear dynamic coefficients to eight. The filter estimates for the direct damping coefficients and all four stiffness coefficients correlated well with actual values, whereas the estimates for the cross-coupled damping coefficients were the least accurate.

Miller, Bard A.↗

Identifying Bearing Rotordynamic Coefficients using an Extended Kalman Filter

An Extended Kalman Filter is developed to estimate the linearized direct and indirect stiffness and damping force coefficients for bearings in rotor-dynamic applications from noisy measurements of the shaft displacement in response to imbalance and impact excitation. The bearing properties are modeled as stochastic random variables using a Gauss-Markov model. Noise terms are introduced into the system model to account for all of the estimation error, including modeling errors and uncertainties and the propagation of measurement errors into the parameter estimates. The system model contains two user-defined parameters that can be tuned to improve the filter s performance; these parameters correspond to the covariance of the system and measurement noise variables. The filter is also strongly influenced by the initial values of the states and the error covariance matrix. The filter is demonstrated using numerically simulated data for a rotor-bearing system with two identical bearings, which reduces the number of unknown linear dynamic coefficients to eight. The filter estimates for the direct damping coefficients and all four stiffness coefficients correlated well with actual values, whereas the estimates for the cross-coupled damping coefficients were the least accurate.

Miller, Brad A.↗

Hypersonic Magneto-Fluid-Dynamic Compression in Cylindrical Inlet

Hypersonic magneto-fluid-dynamic interaction has been successfully performed as a virtual leading-edge strake and a virtual cowl of a cylindrical inlet. In a side-by-side experimental and computational study, the magnitude of the induced compression was found to be depended on configuration and electrode placement. To better understand the interacting phenomenon the present investigation is focused on a direct current discharge at the leading edge of a cylindrical inlet for which validating experimental data is available. The present computational result is obtained by solving the magneto-fluid-dynamics equations at the low magnetic Reynolds number limit and using a nonequilibrium weakly ionized gas model based on the drift-diffusion theory. The numerical simulation provides a detailed description of the intriguing physics. After validation with experimental measurements, the computed results further quantify the effectiveness of a magnet-fluid-dynamic compression for a hypersonic cylindrical inlet. At a minuscule power input to a direct current surface discharge of 8.14 watts per square centimeter of electrode area produces an additional compression of 6.7 percent for a constant cross-section cylindrical inlet.

Shang, Joseph S.↗

Evaluation of Advanced Stirling Convertor Net Heat Input Correlation Methods Using a Thermal Standard

The U.S. Department of Energy (DOE) and Lockheed Martin Space Systems Company (LMSSC) have been developing the Advanced Stirling Radioisotope Generator (ASRG) for use as a power system for space science missions. This generator would use two high-efficiency Advanced Stirling Convertors (ASCs), developed by Sunpower Inc. and NASA Glenn Research Center (GRC). The ASCs convert thermal energy from a radioisotope heat source into electricity. As part of ground testing of these ASCs, different operating conditions are used to simulate expected mission conditions. These conditions require achieving a particular operating frequency, hot end and cold end temperatures, and specified electrical power output for a given net heat input. In an effort to improve net heat input predictions, numerous tasks have been performed which provided a more accurate value for net heat input into the ASCs, including testing validation hardware, known as the Thermal Standard, to provide a direct comparison to numerical and empirical models used to predict convertor net heat input. This validation hardware provided a comparison for scrutinizing and improving empirical correlations and numerical models of ASC-E2 net heat input. This hardware simulated the characteristics of an ASC-E2 convertor in both an operating and non-operating mode. This paper describes the Thermal Standard testing and the conclusions of the validation effort applied to the empirical correlation methods used by the Radioisotope Power System (RPS) team at NASA Glenn.

Briggs, Maxwell H.↗

Simulation of the Impact of New Aircraft- and Satellite-Based Ocean Surface Wind Measurements on H*Wind Analyses and Numerical Forecasts

The Hurricane Imaging Radiometer (HIRAD) is a new airborne microwave remote sensor for hurricane observations that is currently under development by NASA Marshall Space Flight Center, NOAA Hurricane Research Division, the University of Central Florida and the University of Michigan. HIRAD is being designed to enhance the realtime airborne ocean surface winds observation capabilities of NOAA and USAF Weather Squadron hurricane hunter aircraft using the operational airborne Stepped Frequency Microwave Radiometer (SFMR). Unlike SFMR, which measures wind speed and rain rate along the ground track directly beneath the aircraft, HIRAD will provide images of the surface wind and rain field over a wide swath (~ 3 x the aircraft altitude). The present paper describes a set of Observing System Simulation Experiments (OSSEs) in which measurements from the new instrument as well as those from existing instruments (air, surface, and space-based) are simulated from the output of a detailed numerical model, and those results are used to construct H*Wind analyses. The H*Wind analysis, a product of the Hurricane Research Division of NOAA s Atlantic Oceanographic and Meteorological Laboratory, brings together wind measurements from a variety of observation platforms into an objective analysis of the distribution of wind speeds in a tropical cyclone. This product is designed to improve understanding of the extent and strength of the wind field, and to improve the assessment of hurricane intensity. See http://www.aoml.noaa.gov/hrd/data_sub/wind.html. Evaluations will be presented on the impact of the HIRAD instrument on H*Wind analyses, both in terms of adding it to the full suite of current measurements, as well as using it to replace instrument(s) that may not be functioning at the future time the HIRAD instrument is deployed. Plans to demonstrate the potential for HIRAD to improve numerical weather prediction of hurricanes will also be presented.

Miller, Timothy↗

Simulation of the Impact of New Aircraft- and Satellite-Based Ocean Surface Wind Measurements on H*Wind Analyses and Numerical Forecasts

The Hurricane Imaging Radiometer (HIRAD) is a new airborne microwave remote sensor for hurricane observations that is currently under development by NASA Marshall Space Flight Center, NOAA Hurricane Research Division, the University of Central Florida and the University of Michigan. HIRAD is being designed to enhance the realtime airborne ocean surface winds observation capabilities of NOAA and USAF Weather Squadron hurricane hunter aircraft using the operational airborne Stepped Frequency Microwave Radiometer (SFMR). Unlike SFMR, which measures wind speed and rain rate along the ground track directly beneath the aircraft, HIRAD will provide images of the surface wind and rain field over a wide swath (~ 3 x the aircraft altitude). The present paper describes a set of Observing System Simulation Experiments (OSSEs) in which measurements from the new instrument as well as those from existing instruments (air, surface, and space-based) are simulated from the output of a detailed numerical model, and those results are used to construct H*Wind analyses. The H*Wind analysis, a product of the Hurricane Research Division of NOAA s Atlantic Oceanographic and Meteorological Laboratory, brings together wind measurements from a variety of observation platforms into an objective analysis of the distribution of wind speeds in a tropical cyclone. This product is designed to improve understanding of the extent and strength of the wind field, and to improve the assessment of hurricane intensity. See http://www.aoml.noaa.gov/hrd/data_sub/wind.html. Evaluations will be presented on the impact of the HIRAD instrument on H*Wind analyses, both in terms of adding it to the full suite of current measurements, as well as using it to replace instrument(s) that may not be functioning at the future time the HIRAD instrument is implemented. Also shown will be preliminary results of numerical weather prediction OSSEs in which the impact of the addition of HIRAD observations to the initial state on numerical forecasts of the hurricane intensity and structure is assessed.

Miller, Timothy↗

Simulation of the Impact of New Aircraft-and Satellite-based Ocean Surface Wind Measurements on Wind Analyses and Numerical Forecasts

The Hurricane Imaging Radiometer (HIRAD) is a new airborne microwave remote sensor for hurricane observations that is currently under development by NASA Marshall Space Flight Center, NOAA Hurricane Research Division, the University of Central Florida and the University of Michigan. HIRAD is being designed to enhance the realtime airborne ocean surface winds observation capabilities of NOAA and USAF Weather Squadron hurricane hunter aircraft currently using the operational airborne Stepped Frequency Microwave Radiometer (SFMR). Unlike SFMR, which measures wind speed and rain rate along the ground track directly beneath the aircraft, HIRAD will provide images of the surface wind and rain field over a wide swath (approx. 3 x the aircraft altitude). The present paper describes a set of Observing System Simulation Experiments (OSSEs) in which measurements from the new instrument as well as those from existing instruments (air, surface, and space-based) are simulated from the output of a detailed numerical model, and those results are used to construct H*Wind analyses, a product of the Hurricane Research Division of NOAA s Atlantic Oceanographic and Meteorological Laboratory. Evaluations will be presented on the impact of the HIRAD instrument on H*Wind analyses, both in terms of adding it to the full suite of current measurements, as well as using it to replace instrument(s) that may not be functioning at the future time the HIRAD instrument is implemented. Also shown will be preliminary results of numerical weather prediction OSSEs in which the impact of the addition of HIRAD observations to the initial state on numerical forecasts of the hurricane intensity and structure is assessed.

Miller, TImothy↗

Topological structure of the magnetotail as a function of interplanetary magnetic field direction

Magnetic reconnection between the interplanetary magnetic field (IMF) and the geomagnetic field is thought to play a major role in the transfer of solar wind momentum and energy to the magnetosphere. As the angle between the IMF and the geomagnetic field is changed at the bow of the magnetosphere, the topological record of the location of the reconnection region should be recorded in the magnetosheath and on the magnetopause along the flanks of the tail, because the super fast flow freezes strong magnetic gradients formed in the bow reconnection regions into the plasma downstream. In this report, we present results from a three-dimensional, magnetohydrodynamic (MHD), global numerical simulation code for the location of the separatrix between unconnected IMF magnetosheath field lines and reconnected field lines which penetrate the magnetopause and connect to the polar ionosphere. The angle between the IMF direction and the line where the separatrix crosses the magnetopause is shown to be a sensitive function of the IMF clock angle. We also explain how this behavior can be used to derive an approximate relation for the dependence of the cross-polar voltage on the IMF clock angle. We conclude with a note of caution concerning the importance of physical boundary conditions in magnetoplasma simulations.

Fedder, J. A.↗

An analytical transformation for remote sensing of clear-column atmospheric temperature profiles

An exact analytical transformation is presented for the remote sensing of atmospheric temperature profiles in the presence of clouds. The transformation permits direct retrieval of clear-column temperature profiles without the need to predetermine the corresponding clear-column radiances. A numerical illustration is given for simulated observations in the 15 micron CO2 band in the terrestrial atmosphere. The resulting method is numerically very fast and is especially suitable for handling the massive amount of data needed for numerical weather prediction.

Chahine, M. T.↗

Simulated traveling interplanetary disturbances initiated by various solar phenomena

Earlier studies have shown that interplanetary disturbances observed by spacecraft at large heliocentric distances from the sun can be directly associated with various solar activities such as surges, sprays, eruptive prominences, and flares. A digital computer code to simulate these phenomena is developed. This numerical code solves a set of conservation equations for an adiabatic time-dependent spherically-symmetric fluid flow. The solar events are simulated by perturbations in the dependent variables at the lower boundary of the quiet sun corona, and the time-dependent numerical solution simulates the coronal response. Three physical cases are simulated by using the proposed numerical code: (1) part of a surge, spray and prominence; (2) a single flare; and (3) two flares whose interplanetary consequences interact in space. Results for the first two cases are compared with both ground-based and space probe observations. Fairly good agreement with observations is found. A comparison of the third case with observations is presently under way.

Dryer, M.↗

A homotopy continuation method for general preliminary orbit determination and special application to the Tracking and Data Relay Satellite System

A recent mathematical technique for solving systems of equations is applied in a very general way to the orbit determination problem. The use of this technique, the homotopy continuation method, was motivated by the possible need to perform preliminary orbit determination in the Tracking and Data Relay Satellite System (TDRSS), using range and Doppler tracking alone. In the formulation studied here, a set of six tracking observations is continuously transformed from a set with known solutions to the given set with unknown solutions, and the corresponding orbit state vector is followed from the a priori orbit state to the solution orbit states. A modular numerical algorithm is developed to follow the state vector. The technique is formulated and then tested, using both real and simulated tracking, with positive results. Direct comparisons show that the basic method has a much larger radius of convergence than standard or simply modified Newton-Raphson methods. An extension of the basic method is formulated that is expected to produce a global method.

Smith, R. L.↗

Compositional redistribution during casting of Hg sub 0.8 Cd sub 0.2 Te alloys

A series of Hg(0.8)Cd(0.2)Te ingots was cast both vertically and horizontally under well-defined thermal conditions by using a two-zone furnace with isothermal heat-pipe liners. The main objective of the experiments was to establish correlations between casting parameters and compositional redistribution and to develop ground-based data for a proposed flight experiment of casting of Hg(1-x)Cd(x)Te alloys under reduced gravity conditions. The compositional variations along the axial and radial directions were determined by precision density measurements, infrared transmission spectra, and X-ray energy dispersion spectrometry. Comparison between the experimental results and a numerical simulation of the solidification process of Hg(0.8)Cd(0.2)Te is described.

Su, Ching-Hua↗

Compositional redistribution during casting of Hg(0.8)Cd(0.2)Te alloys

A series of Hg(0.8)Cd(0.2)Te ingots was cast both vertically and horizontally under well-defined thermal conditions by using a two-zone furnace with isothermal heat-pipe liners. The main objective of the experiments was to establish correlations between casting parameters and compositional redistribution and to develop ground-based data for a proposed flight experiment of casting of Hg(1-x)Cd(x)Te alloys under reduced gravity conditions. The compositional variations along the axial and radial directions were determined by precision density measurements, infrared transmission spectra, and X-ray energy dispersion spectrometry. Comparison between the experimental results and a numerical simulation of the solidification process of Hg(0.8)Cd(0.2)Te is described.

Su, Ching-Hua↗

Stability and transition in supersonic boundary layers

The three-dimensional time-dependent, compressible Navier-Stokes equations are numerically solved by a Fourier-Chebyshev collocation algorithm to study the stability of a Mach 4.5 flow over a flat plate. Several nonlinear direct simulations suggest the existence of a secondary instability which might provide a possible route to transition. Pertinent differences in the energy content of the various Fourier modes between this instability and the more common incompressible K-type instabilities are pointed out.

Erlebacher, Gordon↗

Banded surface flow maintained by convection in a model of the rapidly rotating giant planets

In three-dimensional numerical simulations of a rapidly rotating Boussinesq fluid shell, thermally driven convection in the form of columns parallel to the rotation axis generates an alternately directed mean zonal flow with a cylindrical structure. The mean structure at the outer spherical surface consists of a broad eastward flow at the equator and alternating bands of westward and eastward flows at higher latitudes in both hemispheres. The banded structure persists even though the underlying convective motions are time-dependent. These results, although still far from the actual motions seen on Jupiter and Saturn, provide support for theoretical suggestions that thermal convection can account for the remarkable banded flow structures on these planets.

Sun, Zi-Ping↗

Computational Investigation of a Boundary-Layer Ingesting Propulsion System for the Common Research Model

The present paper examines potential propulsive and aerodynamic benefits of integrating a Boundary-Layer Ingestion (BLI) propulsion system into a typical commercial aircraft using the Common Research Model (CRM) geometry and the NASA Tetrahedral Unstructured Software System (TetrUSS). The Numerical Propulsion System Simulation (NPSS) environment is used to generate engine conditions for CFD analysis. Improvements to the BLI geometry are made using the Constrained Direct Iterative Surface Curvature (CDISC) design method. Previous studies have shown reductions of up to 25% in terms of propulsive power required for cruise for other axisymmetric geometries using the BLI concept. An analysis of engine power requirements, drag, and lift coefficients using the baseline and BLI geometries coupled with the NPSS model are shown. Potential benefits of the BLI system relating to cruise propulsive power are quantified using a power balance method, and a comparison to the baseline case is made. Iterations of the BLI geometric design are shown and any improvements between subsequent BLI designs presented. Simulations are conducted for a cruise flight condition of Mach 0.85 at an altitude of 38,500 feet and an angle of attack of 2 deg for all geometries. A comparison between available wind tunnel data, previous computational results, and the original CRM model is presented for model verification purposes along with full results for BLI power savings. Results indicate a 14.4% reduction in engine power requirements at cruise for the BLI configuration over the baseline geometry. Minor shaping of the aft portion of the fuselage using CDISC has been shown to increase the benefit from Boundary-Layer Ingestion further, resulting in a 15.6% reduction in power requirements for cruise as well as a drag reduction of eighteen counts over the baseline geometry.

Blumenthal, Brennan T.↗

Investigation of the Polytropic Relationship Between Density and Temperature Within Interplanetary Coronal Mass Ejections Using Numerical Simulations

Single-point spacecraft measurements within coronal mass ejections (CMEs) often exhibit a negative correlation between electron density and temperature. At least two opposing interpretations have been suggested for this relationship. If, on one hand, these single spacecraft observations provide direct measures of the polytropic properties of the plasma, then they imply that the polytropic index for the electrons gamma(sub e) is often < 1. Moreover, since the electrons carry the bulk of the pressure (via their significantly higher temperature), this further implies that the dynamics of CME evolution are dominated by an effective polytropic index gamma(sub e)ff < 1. On the other hand, gamma < 1 implies that as the ejecta propagate away from the Sun and expand, they also heat up; a result clearly at odds with in situ observations. In contrast to these CME intervals, many studies have shown that the quiescent solar wind exhibits a positive correlation between electron density and temperature, suggesting that gamma(sub e) > 1. In this study we simulate the evolution of a variety of CME-like disturbances in the solar wind using a one-dimensional, single-fluid model, to address the interpretation of the relationship between electron density and temperature within CMEs at fixed locations in space. Although we strictly impose a polytropic relationship (with gamma = constant) throughout our simulations, we demonstrate that a variety of correlations can exist between density and temperature at fixed points. Furthermore, we demonstrate that the presence of only local uncorrelated random fluctuations in density and temperature can produce a negative correlation. Consequently, we conclude that these single-point observations of negative correlations between electron density and temperature cannot be used to infer the value of gamma(sub e). Instead, we suggest that entropy variations, together with the plasma's tendency to achieve pressure balance with its surroundings, are responsible for the observed profiles.

Riley, Pete↗

Comparison of remote sensing measurements with a two-scale polarimetric emission and scattering model for sea surfaces

Recently, it has been observed that the brightness temperatures of sea surfaces correlate with the azimuth angle of the ocean wind vector, indicating that ocean wind direction can not only be retrieved from the microwave backscatter, but also from the brightness temperature measurements. In this paper, comparison of the theoretical result calculated from a two-scale emission and scattering model with the Seasat and SSM/I model functions and aircraft measurements is presented and potential applications of microwave polarimetry are discussed. In our two-scale model, the modified reflectivities of large scale surfaces are calculated by extending the small perturbation method to the second order for small scale perturbation with anisotropic directional spectrum. It was found that the modified reflectivities derived from the second-order scattered field agree excellently well with the results obtained from a Monte Carlo simulation technique which numerically calculates the polarimetric reflectivities of one-dimensional random rough surfaces with a power-law spectrum. Without the second-order correction, the modified reflectivities of the rough surfaces are significantly over-estimated and sign errors are observed in the third and fourth Stokes parameters for thermal emissions. The surface spectrum parameters and two-scale cutoff are selected so that the calculated scattering and emission signatures agree with the reported model functions for Seasat and SSM/I. Subsequently, the polarimetric signatures of sea surfaces are illustrated to indicate the possibility of reducing the number of azimuthal looks required for spaceborne sensors in the remote sensing of ocean wind by using the polarimetric information. Furthermore, it is found that contrary to the dependence of backscattering coefficients on incidence angles, polarimetric brightness temperatures display a stronger wind direction dependence in the near nadir-looking direction than away from nadir. Finally, we discuss the potential of a single-look polarimetric radiometer for wind retrieval.

Yueh, S. H.↗