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 109 records · Page 6

Effects of forward velocity on turbulent jet mixing noise

Flight simulation experiments were conducted in an anechoic free jet facility over a broad range of model and free jet velocities. The resulting scaling laws were in close agreement with scaling laws derived from theoretical and semiempirical considerations. Additionally, measurements of the flow structure of jets were made in a wind tunnel by using a laser velocimeter. These tests were conducted to describe the effects of velocity ratio and jet exit Mach number on the development of a jet in a coflowing stream. These turbulence measurements and a simplified Lighthill radiation model were used in predicting the variation in radiated noise at 90 deg to the jet axis with velocity ratio. Finally, the influence of forward motion on flow-acoustic interactions was examined through a reinterpretation of the 'static' numerical solutions to the Lilley equation.

Plumblee, H. E., Jr.↗

Application of turbulence modeling to predict surface heat transfer in stagnation flow region of circular cylinder

A theoretical analysis and numerical calculations for the turbulent flow field and for the effect of free-stream turbulence on the surface heat transfer rate of a stagnation flow are presented. The emphasis is on the modeling of turbulence and its augmentation of surface heat transfer rate. The flow field considered is the region near the forward stagnation point of a circular cylinder in a uniform turbulent mean flow. The free stream is steady and incompressible with a Reynolds number of the order of 10 to the 5th power and turbulence intensity of less than 5 percent. For this analysis, the flow field is divided into three regions: (1) a uniform free-stream region where the turbulence is homogeneous and isotropic; (2) an external viscid flow region where the turbulence is distorted by the variation of the mean flow velocity; and, (3) an anisotropic turbulent boundary layer region over the cylinder surface. The turbulence modeling techniques used are the kappa-epsilon two-equation model in the external flow region and the time-averaged turbulence transport equation in the boundary layer region. The turbulence double correlations, the mean velocity, and the mean temperature within the boundary layer are solved numerically from the transport equations. The surface heat transfer rate is calculated as functions of the free-stream turbulence longitudinal microlength scale, the turbulence intensity, and the Reynolds number.

Wang, Chi R.↗

Ozone vertical distribution retrieval from ground-based high resolution infrared solar spectra

A practical procedure for the retrieval of ozone vertical profiles from ground-based high resolution Fourier transform infrared solar spectra has been developed. The analysis is based on a multilayer line-by-line forward model and a semi-empirical version of the optimal estimation inversion method of Rodgers. The 1002.6-1003.2 cm(exp -1) spectral interval has been selected for the analysis on the basis of synthetic spectrum calculations. This interval contains numerous ozone lines covering a range of intensities and providing retrieval sensitivity from ground level to about 35 km. Characterization of the method and an error analysis have been performed. For a spectral resolution of 0.05-0.01 cm(exp -1) and a signal-to-noise ratio greater than or equal to 100 the retrieval is stable with a vertical resolution of approximately 5 km attainable near the surface degrading to approximately 10 km in the stratosphere. Synthetic spectra studies show that the a priori profile and weak constraints selected for the retrievals introduce no significant biases for a wide range of ozone profiles.

Pougatchev, N. S.↗

Interrogating Solar Flare Loop Models with IRIS Observations 2: Plasma Properties, Energy Transport, and Future Directions

During solar flares a tremendous amount of magnetic energy is released and transported through the Sun’s atmosphere and out into the heliosphere. Despite over a century of study, many unresolved questions surrounding solar flares are still present. Among those are how does the solar plasma respond to flare energy deposition, and what are the important physical processes that transport that energy from the release site in the corona through the transition region and chromosphere? Attacking these questions requires the concert of advanced numerical simulations and high spatial-, temporal-, and spectral-resolution observations. While flares are 3D phenomenon, simulating the NLTE flaring chromosphere in 3D and performing parameter studies of 3D models is largely outwith our current computational capabilities. We instead rely on state-of-the-art 1D field-aligned simulations to study the physical processes that govern flares. Over the last decade, data from the Interface Region Imaging Spectrograph (IRIS) have provided the crucial observations with which we can critically interrogate the predictions of those flare loop models. Here in Paper 2 of a two-part review of IRIS and flare loop models, I discuss how forward modelling flares can help us understand the observations from IRIS, and how IRIS can reveal where our models do well and where we are likely missing important processes, focussing in particular on the plasma properties, energy transport mechanisms, and future directions of flare modelling.

solar flares↗

Leaf bidirectional reflectance and transmittance in corn and soybean

Bidirectional optical properties of leaves must be adequately characterized to develop comprehensive and reliably predictive canopy radiative-transfer models. Directional reflectance and transmittance factors of individual corn and soybean leaves were measured at source incidence angles (SIAs) 20, 45, and 70 deg and numerous view angles in the visible and NIR. Bidirectional reflectance distributions changed with increasing SIA, with forward scattering most pronounced at 70 deg. Directional-hemispherical reflectance generally increased and transmittance decreased with increased SIA. Directional-hemispherical reflectance factors were higher and transmittances were lower than the nadir-viewed reflectance component.

Walter-Shea, E. A.↗

Thermal and Nonthermal X-ray Emission from the Forward Shock in Tycho's Supernova Remnant

We present Chandra CCD images of Tycho's supernova remnant that delineate its outer shock, seen as a thin, smooth rim along the straight northeastern edge and most of the circular western half. The images also show that the Si and S ejecta are highly clumpy, and have reached the forward shock at numerous locations. Most of the X-ray spectra that we examine along the rim show line emission from Si and S, which in some cases must come from ejecta; the continuum is well represented by either thermal or nonthermal models. In the case that the continuum is assumed to be thermal, the temperatures at the rim are all similar at about 2 keV, and the ionization ages are very low because of the overall weakness of the line emission. Assuming shock velocities inferred from radio and X-ray expansion measurements, these temperatures are substantially below those expected for equilibration of the electron and ion temperatures; electron to mean temperature ratios of approximately less than 0.1 - 0.2 indicate at most modest collisionless heating of the electrons at the shock. The nonthermal contribution to these spectra may be important, however, and may account for as many as half of the counts in the 4-6 keV energy range, based on an extrapolation of the hard X-ray spectrum above 10 keV.

Hwang, Una↗

Towards developing robust algorithms for solving partial differential equations on MIMD machines

Methods for efficient computation of numerical algorithms on a wide variety of MIMD machines are proposed. These techniques reorganize the data dependency patterns to improve the processor utilization. The model problem finds the time-accurate solution to a parabolic partial differential equation discretized in space and implicitly marched forward in time. The algorithms are extensions of Jacobi and SOR. The extensions consist of iterating over a window of several timesteps, allowing efficient overlap of computation with communication. The methods increase the degree to which work can be performed while data are communicated between processors. The effect of the window size and of domain partitioning on the system performance is examined both by implementing the algorithm on a simulated multiprocessor system.

Saltz, J. H.↗

Towards developing robust algorithms for solving partial differential equations on MIMD machines

Methods for efficient computation of numerical algorithms on a wide variety of MIMD machines are proposed. These techniques reorganize the data dependency patterns to improve the processor utilization. The model problem finds the time-accurate solution to a parabolic partial differential equation discretized in space and implicitly marched forward in time. The algorithms are extensions of Jacobi and SOR. The extensions consist of iterating over a window of several timesteps, allowing efficient overlap of computation with communication. The methods increase the degree to which work can be performed while data are communicated between processors. The effect of the window size and of domain partitioning on the system performance is examined both by implementing the algorithm on a simulated multiprocessor system.

Saltz, Joel H.↗

A new methodology for free wake analysis using curved vortex elements

A method using curved vortex elements was developed for helicopter rotor free wake calculations. The Basic Curve Vortex Element (BCVE) is derived from the approximate Biot-Savart integration for a parabolic arc filament. When used in conjunction with a scheme to fit the elements along a vortex filament contour, this method has a significant advantage in overall accuracy and efficiency when compared to the traditional straight-line element approach. A theoretical and numerical analysis shows that free wake flows involving close interactions between filaments should utilize curved vortex elements in order to guarantee a consistent level of accuracy. The curved element method was implemented into a forward flight free wake analysis, featuring an adaptive far wake model that utilizes free wake information to extend the vortex filaments beyond the free wake regions. The curved vortex element free wake, coupled with this far wake model, exhibited rapid convergence, even in regions where the free wake and far wake turns are interlaced. Sample calculations are presented for tip vortex motion at various advance ratios for single and multiple blade rotors. Cross-flow plots reveal that the overall downstream wake flow resembles a trailing vortex pair. A preliminary assessment shows that the rotor downwash field is insensitive to element size, even for relatively large curved elements.

Bliss, Donald B.↗

An analysis of the flow field in the region of the ASRM field joints

The flow field in the region of a solid rocket motor field joint is very important since fluid dynamic and mechanical propellant stresses can couple to cause a motor failure at a joint. Presented here is an examination of the flow field in the region of the Advanced Solid Rocket Motor (ASRM) field joints. The analyses were performed as a first step in assessing the design of the ASRM forward and aft field joints in order to assure the proper operation of the motor prior to further development of test firing. The analyses presented here were performed by employing a two-dimensional axisymmetric assumption. Fluent/BFC, a three dimensional full Navier-Stokes flow field code, was used to make the numerical calculations. This code utilizes a staggered grid formulation along with the SIMPLER numerical algorithm. Wall functions are used to determine the character of the laminar sublayer, and a standard kappa-epsilon turbulence model is used to close the fluid dynamic equations. The analyses performed to this date verify that the ASRM field joint design operates properly. The fluid dynamic stresses at the field joints are small due to the inherent design of the field joints. A problem observed in some other solid rocket motors is that large fluid dynamic stresses are generated at the motor joint on the downstream propellant grain due to forward facing step geometries. The design of the ASRM field joints are such that this is not a problem as shown by the analyses. Also, the analyses of the inhibitor stub left protruding into the port flow from normal propellant burn back show that more information is necessary to complete these analyses. These analyses were performed as parametric analyses in relation to the height of the inhibitor stub left protruding into the motor port. A better estimate of the amount of the inhibitor stub remaining at later burn times must be determined since the height which the inhibitor stub protrudes into the port flow drastically affects the fluid dynamics induced stresses on the propellant grain at the field joints.

Dill, Richard A.↗

Experimental Measurement of RCS Jet Interaction Effects on a Capsule Entry Vehicle

An investigation was made in NASA Langley Research Center s 31-Inch Mach 10 Tunnel to determine the effects of reaction-control system (RCS) jet interactions on the aft-body of a capsule entry vehicle. The test focused on demonstrating and improving advanced measurement techniques that would aid in the rapid measurement and visualization of jet interaction effects for the Orion Crew Exploration Vehicle while providing data useful for developing engineering models or validation of computational tools used to assess actual flight environments. Measurements included global surface imaging with pressure and temperature sensitive paints and three-dimensional flow visualization with a scanning planar laser induced fluorescence technique. The wind tunnel model was fabricated with interchangeable parts for two different aft-body configurations. The first, an Apollo-like configuration, was used to focus primarily on the forward facing roll and yaw jet interactions which are known to have significant aft-body heating augmentation. The second, an early Orion Crew Module configuration (4-cluster jets), was tested blowing only out of the most windward yaw jet, which was expected to have the maximum heating augmentation for that configuration. Jet chamber pressures and tunnel flow conditions were chosen to approximate early Apollo wind tunnel test conditions. Maximum heating augmentation values measured for the Apollo-like configuration (>10 for forward facing roll jet and 4 for yaw jet) using temperature sensitive paint were shown to be similar to earlier experimental results (Jones and Hunt, 1965) using a phase change paint technique, but were acquired with much higher surface resolution. Heating results for the windward yaw jet on the Orion configuration had similar augmentation levels, but affected much less surface area. Numerical modeling for the Apollo-like yaw jet configuration with laminar flow and uniform jet outflow conditions showed similar heating patterns, qualitatively, but also showed significant variation with jet exit divergence angle, with as much as 25 percent variation in heat flux intensity for a 10 degree divergence angle versus parallel outflow. These results along with the fabrication methods and advanced measurement techniques developed will be used in the next phase of testing and evaluation for the updated Orion RCS configuration.

Buck, Gregory M.↗

The Role of Model and Initial Condition Error in Numerical Weather Forecasting Investigated with an Observing System Simulation Experiment

A series of experiments that explore the roles of model and initial condition error in numerical weather prediction are performed using an observing system simulation experiment (OSSE) framework developed at the National Aeronautics and Space Administration Global Modeling and Assimilation Office (NASA/GMAO). The use of an OSSE allows the analysis and forecast errors to be explicitly calculated, and different hypothetical observing networks can be tested with ease. In these experiments, both a full global OSSE framework and an 'identical twin' OSSE setup are utilized to compare the behavior of the data assimilation system and evolution of forecast skill with and without model error. The initial condition error is manipulated by varying the distribution and quality of the observing network and the magnitude of observation errors. The results show that model error has a strong impact on both the quality of the analysis field and the evolution of forecast skill, including both systematic and unsystematic model error components. With a realistic observing network, the analysis state retains a significant quantity of error due to systematic model error. If errors of the analysis state are minimized, model error acts to rapidly degrade forecast skill during the first 24-48 hours of forward integration. In the presence of model error, the impact of observation errors on forecast skill is small, but in the absence of model error, observation errors cause a substantial degradation of the skill of medium range forecasts.

Model Error↗

The influence of dynamic inflow and torsional flexibility on rotor damping in forward flight from symbolically generated equations

The combined effects of blade torsion and dynamic inflow on the aeroelastic stability of an elastic rotor blade in forward flight are studied. The governing sets of equations of motion (fully nonlinear, linearized, and multiblade equations) used in this study are derived symbolically using a program written in FORTRAN. Stability results are presented for different structural models with and without dynamic inflow. A combination of symbolic and numerical programs at the proper stage in the derivation process makes the obtainment of final stability results an efficient and straightforward procedure.

Reddy, T. S. R.↗

Navier-Stokes Solutions of Rotor and Rotor-Body Flows

This paper will review the advances made recently in the Navier-Stokes CFD methods to simulate aerodynamics and aeroacoustics of helicopter rotors and rotor-body flows. Although a complete flowfield simulation of full helicopter is currently not feasible with these methods, impressive gains have been made in analyzing individual components of this complex problem in a very detailed manner. The use of the state-of-the-art numerical algorithms in solution methods, in conjunction with powerful supercomputers, like the Cray-2, have enabled noticeable progress to be made in modeling viscous-inviscid interactions, blade-vortex interactions, tip-vortex: simulation and wake effects, as well as high speed impulsive noise in hover and forward flight for isolated rotor blades. This paper will critically evaluate the presently available Euler and Navier-Stokes methods, both finite-difference and finite volume methods using structured and unstructured grids for helicopter applications for accuracy, suitability, and computational efficiency. The review will also include the recent progress made using overset grids to model rotor-body flows. All the material for this review will be drawn from the published material shown below.

Srinivasan, G. R.↗

Turbulence and surface heat transfer near the stagnation point of a circular cylinder in turbulent flow

A turbulent boundary layer flow analysis of the momentum and thermal flow fields near the forward stagnation point due to a circular cylinder in turbulent cross flow is presented. Turbulence modeling length scale, anisotropic turbulence initial profiles and boundary conditions were identified as functions of the cross flow turbulence intensity and the boundary layer flow far field velocity. These parameters were used in a numerical computational procedure to calculate the mean velocity, mean temperature, and turbulence double correlation profiles within the flow field. The effects of the cross flow turbulence on the stagnation region momentum and thermal flow fields were investigated. This analysis predicted the existing measurements of the stagnation region mean velocity and surface heat transfer rate with cross flow Reynolds number and turbulence intensity less than 250,000 and 0.05, respectively.

Wang, C. R.↗

Linearization of manipulator dynamics using spatial operators

Linearized dynamics models for manipulators are useful in robot analysis, motion planning, and control applications. Techniques from the spatial operator algebra are used to obtain closed form operator expressions for two types of linearized dynamics models, the linearized inverse and forward dynamics models. Spatially recursive algorithms of O(n) and O(n-squared) complexity for the computation of the perturbation vector and coefficient matrices for the linearized inverse dynamics model are developed first. Subsequently, operator factorization and inversion identities are used to develop corresponding closed-form expressions for the linearized forward dynamics model (LFDM). Once again, these are used to develop algorithms of O(n) and O(n-squared) complexity for the computation of the perturbation vector and the coefficient matrices. The algorithms for the LFDM do not require the explicit computation of the mass matrix nor its numerical inversion and are also of lower complexity than the conventional O(n-cubed) algorithms.

Jain, A.↗

Electromagnetic wave extinction within a forested canopy

A forested canopy is modeled by a collection of randomly oriented finite-length cylinders shaded by randomly oriented and distributed disk- or needle-shaped leaves. For a plane wave exciting the forested canopy, the extinction coefficient is formulated in terms of the extinction cross sections (ECSs) in the local frame of each forest component and the Eulerian angles of orientation (used to describe the orientation of each component). The ECSs in the local frame for the finite-length cylinders used to model the branches are obtained by using the forward-scattering theorem. ECSs in the local frame for the disk- and needle-shaped leaves are obtained by the summation of the absorption and scattering cross-sections. The behavior of the extinction coefficients with the incidence angle is investigated numerically for both deciduous and coniferous forest. The dependencies of the extinction coefficients on the orientation of the leaves are illustrated numerically.

Karam, M. A.↗

On the Limitations of Variational Bias Correction

Satellite radiances are the largest dataset assimilated into Numerical Weather Prediction (NWP) models, however the data are subject to errors and uncertainties that need to be accounted for before assimilating into the NWP models. Variational bias correction uses the time series of observation minus background to estimate the observations bias. This technique does not distinguish between the background error, forward operator error, and observations error so that all these errors are summed up together and counted as observation error. We identify some sources of observations errors (e.g., antenna emissivity, non-linearity in the calibration, and antenna pattern) and show the limitations of variational bias corrections on estimating these errors.

Moradi, Isaac↗