Search NASA⌕ Search

SEARCH · Search NASA

Results for “numerical forward model”

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 37 records · Page 2

Tip Vortex and Wake Characteristics of a Counterrotating Open Rotor

One of the primary noise sources for Open Rotor systems is the interaction of the forward rotor tip vortex and blade wake with the aft rotor. NASA has collaborated with General Electric on the testing of a new generation of low noise, counterrotating Open Rotor systems. Three-dimensional particle image velocimetry measurements were acquired in the intra-rotor gap of the Historical Baseline blade set. The velocity measurements are of sufficient resolution to characterize the tip vortex size and trajectory as well as the rotor wake decay and turbulence character. The tip clearance vortex trajectory is compared to results from previously developed models. Forward rotor wake velocity profiles are shown. Results are presented in a form as to assist numerical modeling of Open Rotor system aerodynamics and acoustics.

VanZante, Dale E.↗

Tip Vortex and Wake Characteristics of a Counterrotating Open Rotor

One of the primary noise sources for Open Rotor systems is the interaction of the forward rotor tip vortex and blade wake with the aft rotor. NASA has collaborated with General Electric on the testing of a new generation of low noise, counterrotating Open Rotor systems. Three-dimensional particle image velocimetry measurements were acquired in the intra-rotor gap of the Historical Baseline blade set. The velocity measurements are of sufficient resolution to characterize the tip vortex size and trajectory as well as the rotor wake decay and turbulence character. The tip clearance vortex trajectory is compared to results from previously developed models. Forward rotor wake velocity profiles are shown. Results are presented in a form as to assist numerical modeling of Open Rotor system aerodynamics and acoustics.

Van Zante, Dale E.↗

(GO)2-SIM: a GCM-Oriented Ground-Observation Forward-Simulator Framework for Objective Evaluation of Cloud and Precipitation Phase

General circulation model (GCM) evaluation using ground-based observations is complicated by inconsistencies in hydrometeor and phase definitions. Here we describe (GO)2-SIM, a forward simulator designed for objective hydrometeor-phase evaluation, and assess its performance over the North Slope of Alaska using a 1-year GCM simulation. For uncertainty assessment, 18 empirical relationships are used to convert model grid-average hydrometeor (liquid and ice, cloud, and precipitation) water contents to zenith polarimetric micropulse lidar and Ka-band Doppler radar measurements, producing an ensemble of 576 forward-simulation realizations. Sensor limitations are represented in forward space to objectively remove from consideration model grid cells with undetectable hydrometeor mixing ratios, some of which may correspond to numerical noise.Phase classification in forward space is complicated by the inability of sensors to measure ice and liquid signals distinctly. However, signatures exist in lidar–radar space such that thresholds on observables can be objectively estimated and related to hydrometeor phase. The proposed phase-classification technique leads to misclassification in fewer than 8% of hydrometeor-containing grid cells. Such misclassifications arise because, while the radar is capable of detecting mixed-phase conditions, it can mistake water- for ice-dominated layers. However, applying the same classification algorithm to forward-simulated and observed fields should generate hydrometeor-phase statistics with similar uncertainty. Alternatively, choosing to disregard how sensors define hydrometeor phase leads to frequency of occurrence discrepancies of up to 40%. So, while hydrometeor-phase maps determined in forward space are very different from model "reality" they capture the information sensors can provide and thereby enable objective model evaluation.

Lamer, K.↗

Helicopter rotor wake geometry and its influence in forward flight. Volume 1: Generalized wake geometry and wake effect on rotor airloads and performance

An analytic investigation to generalize wake geometry of a helicopter rotor in steady level forward flight and to demonstrate the influence of wake deformation in the prediction of rotor airloads and performance is described. Volume 1 presents a first level generalized wake model based on theoretically predicted tip vortex geometries for a selected representative blade design. The tip vortex distortions are generalized in equation form as displacements from the classical undistorted tip vortex geometry in terms of vortex age, blade azimuth, rotor advance ratio, thrust coefficient, and number of blades. These equations were programmed to provide distorted wake coordinates at very low cost for use in rotor airflow and airloads prediction analyses. The sensitivity of predicted rotor airloads, performance, and blade bending moments to the modeling of the tip vortex distortion are demonstrated for low to moderately high advance ratios for a representative rotor and the H-34 rotor. Comparisons with H-34 rotor test data demonstrate the effects of the classical, predicted distorted, and the newly developed generalized wake models on airloads and blade bending moments. Use of distorted wake models results in the occurrence of numerous blade-vortex interactions on the forward and lateral sides of the rotor disk. The significance of these interactions is related to the number and degree of proximity to the blades of the tip vortices. The correlation obtained with the distorted wake models (generalized and predicted) is encouraging.

Egolf, T. A.↗

Numerical simulation and comparison of symmetrical/supercritical airfoils for the near tip region of a helicopter in forward flight

Aerodynamic loads on a multi-bladed helicopter rotor in forward flight at transonic tip conditions are calculated. The unsteady, three-dimensional, time-accurate compressible Reynolds-averaged thin layer Navier-Stokes equations are solved in a rotating coordinate system on a body-conformed, curvilinear grid of C-H topology. Detailed boundary layer and global numerical comparisons of NACA-0012 symmetrical and CAST7-158 supercritical airfoils are made under identical forward flight conditions. The rotor wake effects are modeled by applying a correction to the geometric angle of attack of the blade. This correction is obtained by computing the local induced downwash velocity with a free wake analysis program. The calculations are performed on the Numerical Aerodynamic Simulation Cray 2 and the VPS32 (a derivative of a Cyber 205 at the Langley Research Center) for a model helicopter rotor in forward flight.

Badavi, F. F.↗

Program Computes Flows Of Fluids And Heat

SINDA'85/FLUINT incorporates lumped-parameter-network and one-dimensional-flow mathematical models. System enables analysis of mutual influences of thermal and flow phenomena. Offers two finite-difference numerical solution techniques: forward-difference explicit approximation and Crank-Nicholson approximation. Enables simulation of nonuniform heating and facilitates mathematical modeling of thin-walled heat exchangers. Ability to model nonequilibrium behavior within two-phase volumes included. Recent changes in program improve modeling of real evaporator pumps and other capillary-assist evaporators. Written in FORTRAN 77.

Cullimore, Brent↗

Solar Cycle Predictions Near Sunspot Minimum

Observations of solar magnetic activity and the dynamics of the solar convection zone have produced severe constraints on models of the Sun's magnetic dynamo. These constraints are so severe that, at present, we do not have numerical models that can accept the current conditions and then march forward in time to predict future activity. Given this state of solar dynamo theory we are forced to examine previous behavior to discover patterns and trends that afford us some measure of predictability. Here we examine the behavior of several indicators of solar activity near solar minimum that are well correlated with the amplitude of the following solar maximum to predict the level of solar activity over cycle 23. Sunspot numbers, areas, and positions are useful for characterizing solar cycle behavior due to the extent of the data (12 cycles or more). These data exhibit several patterns that relate future activity to past behavior. With the Odd-Even effect the odd numbered cycles have been larger than their even numbered predecessors for each of the last six cycle pairs. With the Amplitude-Period effect short period cycles have been followed by large amplitude cycles and long period cycles have been followed by small amplitude cycles for 10 of the last 13 cycles. With the Maximum-Minimum effect the sunspot number at minimum is directly correlated with the sunspot number at maximum for a given cycle. The geomagnetic indices aa and Ap are also related to solar activity by the connections between disturbances in the solar wind and variations in the Earth's magnetic field. Like the Maximum-Minimum effect for sunspots, the size of the aa and Ap indices at minimum are directly related to the amplitude of the following maximum. The number of geomagnetically disturbed days (days with Ap >= 25) over the course of a cycle is another indicator for the size of the next cycle. The aa and Ap indices can each be separated into a component in phase with the current sunspot cycle and an "extended cycle" component that is directly related to the following cycle. Individually each of these indicators provides some measure of predictive ability for the amplitudes of previous cycles. Taken in combination we find that the predictions are much more accurate. Our multi-variate analysis indicates that a maximum sunspot number of 160 +/- 20 should occur near the end of 1999 for cycle 23.

Hathaway, David H.↗

Response of multi-panel assembly to noise from a jet in forward motion

A model of the interaction of the noise from a spreading subsonic jet with a 4 panel assembly is studied numerically in two dimensions. The effect of forward motion of the jet is accounted for by considering a uniform flow field superimposed on a mean jet exit profile. The jet is initially excited by a pulse-like source inserted into the flow field. The pulse triggers instabilities associated with the inviscid instability of the jet shear layer. These instabilities generate sound which in turn serves to excite the panels. We compare the sound from the jet, the responses of the panels and the resulting acoustic radiation for the static jet and the jet in forward motion. The far field acoustic radiation, the panel response and sound radiated from the panels are all computed and compared to computations of a static jet. The results demonstrate that for a jet in forward motion there is a reduction in sound in downstream directions and an increase in sound in upstream directions in agreement with experiments. Furthermore, the panel response and radiation for a jet in forward motion exhibits a downstream attenuation as compared with the static case.

Bayliss, A.↗

Forward Monte Carlo Computations of Polarized Microwave Radiation

Microwave radiative transfer computations continue to acquire greater importance as the emphasis in remote sensing shifts towards the understanding of microphysical properties of clouds and with these to better understand the non linear relation between rainfall rates and satellite-observed radiance. A first step toward realistic radiative simulations has been the introduction of techniques capable of treating 3-dimensional geometry being generated by ever more sophisticated cloud resolving models. To date, a series of numerical codes have been developed to treat spherical and randomly oriented axisymmetric particles. Backward and backward-forward Monte Carlo methods are, indeed, efficient in this field. These methods, however, cannot deal properly with oriented particles, which seem to play an important role in polarization signatures over stratiform precipitation. Moreover, beyond the polarization channel, the next generation of fully polarimetric radiometers challenges us to better understand the behavior of the last two Stokes parameters as well. In order to solve the vector radiative transfer equation, one-dimensional numerical models have been developed, These codes, unfortunately, consider the atmosphere as horizontally homogeneous with horizontally infinite plane parallel layers. The next development step for microwave radiative transfer codes must be fully polarized 3-D methods. Recently a 3-D polarized radiative transfer model based on the discrete ordinate method was presented. A forward MC code was developed that treats oriented nonspherical hydrometeors, but only for plane-parallel situations.

Battaglia, A.↗

Models for Water Isotopes Constrained with Data from Crystal Face

During the year covered by this proposal we conducted work on several different topics, as reflected by our publications listed below. One major activity was to work with a group of about 10 scientists from around the country to prepare a science-planning document (Tropical Composition, Cloud and Climate Coupling Experiment (TC4)) that outlined the rationale, locations, strategy to accomplish the goals, and possible payloads for a set of three tropical missions. We also prepared background materials for various NRAs being prepared at NASA Headquarters for missions in Costa Rica, Darwin and Guam. Unfortunately budgetary constraints prevented these missions from moving forward. In conjunction with the group NASA Ames we built a new numerical model for deep convection and have applied that model to simulate the CRYSTAL isotope data. Our goal in particular has been to better understanding how convection distributes water vapor isotopes. CRYSTAL observations of water isotopes are very different from those suggested by previous workers who assumed the isotopes would obey Rayleigh fractionation. The water isotope study has several implications. First it is a check on the realism of the deep convection model. Second, the isotopes are a measure of the precipitation removal in the atmosphere. Hence they provide a constraint on a parameter that is difficult to otherwise measure. Finally it has been suggested that isotopes may be the key to unraveling the water transport into the stratosphere and upper troposphere. Such transport is critical both for the radiation balance and for stratospheric chemistry. Ours is the first model that is able to treat this transport. Our initial results are now in press in Geophys. Res. Lett. Essentially we are able to explain the vertical profiles of isotopes in the tropical tropopause transition layer. We are also able to account for stratospheric humidity ana isotope abundances with this model. The data suggest that isotopes do not provide a clear constraint on the mechanism by which water enters the stratosphere- whether by convection, or by slow ascent. Our work is relevant for the water isotope comparison experiments recently done by NASA. We are conducting numerical experiments related to this project to help understand that data.

Toon, Owen B.↗

Stability of numerical integration techniques for transient rotor dynamics

A finite element model of a rotor bearing system was analyzed to determine the stability limits of the forward, backward, and centered Euler; Runge-Kutta; Milne; and Adams numerical integration techniques. The analysis concludes that the highest frequency mode determines the maximum time step for a stable solution. Thus, the number of mass elements should be minimized. Increasing the damping can sometimes cause numerical instability. For a uniform shaft, with 10 mass elements, operating at approximately the first critical speed, the maximum time step for the Runge-Kutta, Milne, and Adams methods is that which corresponds to approximately 1 degree of shaft movement. This is independent of rotor dimensions.

Kascak, A. F.↗

Inversion of Magnetic Measurements of the CHAMP Satellite Over the Pannonian Basin

The Pannonian Basin is a deep intra-continental basin that formed as part of the Alpine orogeny. In order to study the nature of the crustal basement we used the long-wavelength magnetic anomalies acquired by the CHAMP satellite. The anomalies were distributed in a spherical shell, some 107,927 data recorded between January 1 and December 31 of 2008. They covered the Pannonian Basin and its vicinity. These anomaly data were interpolated into a spherical grid of 0.5 x 0.5, at the elevation of 324 km by the Gaussian weight function. The vertical gradient of these total magnetic anomalies was also computed and mapped to the surface of a sphere at 324 km elevation. The former spherical anomaly data at 425 km altitude were downward continued to 324 km. To interpret these data at the elevation of 324 km we used an inversion method. A polygonal prism forward model was used for the inversion. The minimum problem was solved numerically by the Simplex and Simulated annealing methods; a L2 norm in the case of Gaussian distribution parameters and a L1 norm was used in the case of Laplace distribution parameters. We INTERPRET THAT the magnetic anomaly WAS produced by several sources and the effect of the sable magnetization of the exsolution of hemo-ilmenite minerals in the upper crustal metamorphic rocks.

Kis, K. I.↗

An explicit mixed numerical method for mesoscale model

A mixed numerical method has been developed for mesoscale models. The technique consists of a forward difference scheme for time tendency terms, an upstream scheme for advective terms, and a central scheme for the other terms in a physical system. It is shown that the mixed method is conditionally stable and highly accurate for approximating the system of either shallow-water equations in one dimension or primitive equations in three dimensions. Since the technique is explicit and two time level, it conserves computer and programming resources.

Hsu, H.-M.↗

Low gravity transfer line chilldown

A code has been developed that solves for the transfer line chilldown time and flow and heat transfer characteristics in one-g environment. The code solves the transient, one dimensional, space averaged mass, momentum, and energy conservation equations for liquid-vapor two-phase flow in tubes. The physical configuration solved is that appropriate for bottom coolant injection in a vertically supported heated tube. Four distinct regions are considered consecutively: fully liquid, inverted annular, dispersed, and fully vapor flow. The conservation equations for both the liquid and the vapor are solved in each region separately. Also, in each region the mass and energy transport between each phase as well as the energy and momentum transport between the tube wall and the fluid are accounted for. A finite wall thickness is also considered. The model described above was solved numerically through a mixed finite difference scheme with forward time marching. The inverted annular regime was resolved using a semi-implicit finite differencing while the dispersed regime was solved explicitly. Also, a staggered mesh was used in which the velocity was resolved at mesh boundaries while all other field variables were resolved at the mesh centroids. Different mesh sizes were used depending on the region of solution. A coarse mesh was used in the dispersed flow region while a much finer mesh was used in both the inverted annular flow region and the tube walls.

B. N. Antar↗

Freezable Radiator Model Correlation Improvements and Fluids Study

Freezable radiators offer an attractive solution to the issue of thermal control system scalability. As thermal environments change, a freezable radiator will effectively scale the total heat rejection it is capable of as a function of the thermal environment and flow rate through the radiator. Scalable thermal control systems are a critical technology for spacecraft that will endure missions with widely varying thermal requirements. These changing requirements are a result of the space craft s surroundings and because of different thermal rejection requirements during different mission phases. However, freezing and thawing (recovering) a radiator is a process that has historically proven very difficult to predict through modeling, resulting in highly inaccurate predictions of recovery time. To attempt to improve this, tests were conducted in 2009 to determine whether the behavior of a simple stagnating radiator could be predicted or emulated in a Thermal Desktop(trademark) numerical model. A 50-50 mixture of DowFrost HD and water was used as the working fluid. Efforts to scale this model to a full scale design, as well as efforts to characterize various thermal control fluids at low temperatures are also discussed. Previous testing and modeling efforts showed that freezable radiators could be operated as intended, and be fairly, if not perfectly predicted by numerical models. This paper documents the improvements made to the numerical model, and outcomes of fluid studies that were determined necessary to go forward with further radiator testing.

Lillibridge, Sean↗

Numerical evaluation of the jet noise source distribution from far-field cross correlations

This paper contains the development of techniques to determine the relationship between the unknown source correlation function to the correlation of scattered amplitudes in a jet. This study has application to the determination of forward motion effects. The technique has been developed and tested on a model jet of high subsonic flow. Numerical solution was obtained by solving the Fredholm integral equation of the first kind. Interpretation of the apparent source distribution and its application to flight testing are provided.

Maestrello, L.↗

Can Aerosol Direct Radiative Effects Account for Analysis Increments of Temperature in the Tropical Atlantic?

In the late 1990's, prior to the launch of the Terra satellite, atmospheric general circulation models (GCMs) did not include aerosol processes because aerosols were not properly monitored on a global scale and their spatial distributions were not known well enough for their incorporation in operational GCMs. At the time of the first GEOS Reanalysis (Schubert et al. 1993), long time series of analysis increments (the corrections to the atmospheric state by all available meteorological observations) became readily available, enabling detailed analysis of the GEOS-1 errors on a global scale. Such analysis revealed that temperature biases were particularly pronounced in the Tropical Atlantic region, with patterns depicting a remarkable similarity to dust plumes emanating from the African continent as evidenced by TOMS aerosol index maps. Yoram Kaufman was instrumental encouraging us to pursue this issue further, resulting in the study reported in Alpert et al. (1998) where we attempted to assess aerosol forcing by studying the errors of a the GEOS-1 GCM without aerosol physics within a data assimilation system. Based on this analysis, Alpert et al. (1998) put forward that dust aerosols are an important source of inaccuracies in numerical weather-prediction models in the Tropical Atlantic region, although a direct verification of this hypothesis was not possible back then. Nearly 20 years later, numerical prediction models have increased in resolution and complexity of physical parameterizations, including the representation of aerosols and their interactions with the circulation. Moreover, with the advent of NASA's EOS program and subsequent satellites, atmospheric aerosols are now monitored globally on a routine basis, and their assimilation in global models are becoming well established. In this talk we will reexamine the Alpert et al. (1998) hypothesis using the most recent version of the GEOS-5 Data Assimilation System with assimilation of aerosols. We will explicitly calculate the impact of aerosols on the temperature analysis increments in the tropical Atlantic and assess the extent to which inclusion of atmospheric aerosols have reduced these increments.

General Circulation Models↗

Numerical study of three-dimensional turbulent flow interactions between blockage models and wind tunnels including longitudinally slotted test sections

A spatial forward-marching approach is applied to compute three-dimensional turbulent flows for several blockage models in free flight, in a solid wall wind tunnel, and in a wind tunnel with longitudinal slots in the test section. The effects of area blockage in the tunnel, model growth and tunnel wall boundary layers, and of the slots are included. The large blockage models are found to have significant wall interference effects which can be reduced by the slots. The effects of the latter are confined to the region near the tunnel wall. Model/wall interference effects are not limited to the effects of area blockage; in particular, boundary layer profile shapes for a wind tunnel model in a tunnel are different from shapes for a model in free flight even when slots are used. This indicates that the flow responds differently in these two cases with the same pressure gradient.

Lin, S.-J.↗