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 163 records · Page 9

Technical Challenges Associated with In-Air Wingtip Docking of Aircraft in Forward Flight

Autonomous in-air wingtip docking of aircraft offers significant opportunity for system level performance gains for numerous aircraft applications. Several of the technical challenges facing wingtip docking of fixed-wing aircraft are addressed in this paper, including: close proximity aerodynamic coupling; mechanisms and operations for robust docking; and relative state estimation methods. A simulation framework considering the aerodynamics, rigid-body dynamics, and vehicle controls is developed and used to perform docking sensitivity studies for a system of two 5.5% scale NASA Generic Transport Model aircraft. Additionally, proof of- concept testing of a candidate docking mechanism designed to move the primary wingtip vortex inboard suggests the viability of such an approach for achieving robust docking.

Quinlan, Jesse R.↗

Bayesian Approach to the Joint Inversion of Gravity and Magnetic Data, with Application to the Ismenius Area of Mars

This viewgraph presentation reviews a Bayesian approach to the inversion of gravity and magnetic data with specific application to the Ismenius Area of Mars. Many inverse problems encountered in geophysics and planetary science are well known to be non-unique (i.e. inversion of gravity the density structure of a body). In hopes of reducing the non-uniqueness of solutions, there has been interest in the joint analysis of data. An example is the joint inversion of gravity and magnetic data, with the assumption that the same physical anomalies generate both the observed magnetic and gravitational anomalies. In this talk, we formulate the joint analysis of different types of data in a Bayesian framework and apply the formalism to the inference of the density and remanent magnetization structure for a local region in the Ismenius area of Mars. The Bayesian approach allows prior information or constraints in the solutions to be incorporated in the inversion, with the "best" solutions those whose forward predictions most closely match the data while remaining consistent with assumed constraints. The application of this framework to the inversion of gravity and magnetic data on Mars reveals two typical challenges - the forward predictions of the data have a linear dependence on some of the quantities of interest, and non-linear dependence on others (termed the "linear" and "non-linear" variables, respectively). For observations with Gaussian noise, a Bayesian approach to inversion for "linear" variables reduces to a linear filtering problem, with an explicitly computable "error" matrix. However, for models whose forward predictions have non-linear dependencies, inference is no longer given by such a simple linear problem, and moreover, the uncertainty in the solution is no longer completely specified by a computable "error matrix". It is therefore important to develop methods for sampling from the full Bayesian posterior to provide a complete and statistically consistent picture of model uncertainty, and what has been learned from observations. We will discuss advanced numerical techniques, including Monte Carlo Markov

data analysis↗

Electromagnetic Scattering Analysis of Large Size Asteroids/Comets for Reflection/Transmission Tomography (RTT)

A precise knowledge of the interior structure of asteroids, comets, and Near Earth Objects (NEO) is important to assess the consequences of their impacts with the Earth and develop efficient mitigation strategies. Knowledge of their interior structure also provides opportunities for extraction of raw materials for future space activities. Low frequency radio sounding is often proposed for investigating interior structures of asteroids and NEOs. For designing and optimizing radio sounding instrument it is advantageous to have an accurate and efficient numerical simulation model of radio reflection and transmission through large size bodies of asteroid shapes. In this presentation we will present electromagnetic (EM) scattering analysis of electrically large size asteroids using (1) a weak form formulation and (2) also a more accurate hybrid finite element method/method of moments (FEM/MOM) to help estimate their internal structures. Assuming the internal structure with known electrical properties of a sample asteroid, we first develop its forward EM scattering model. From the knowledge of EM scattering as a function of frequency and look angle we will then present the inverse scattering procedure to extract its interior structure image. Validity of the inverse scattering procedure will be presented through few simulation examples.

Deshpande, Manohar↗

A numerical study of the unsteady leading edge separation bubble on an oscillating airfoil

A preliminary study of the unsteady viscous flow in the region of an airfoil leading edge was performed, in which the interaction between the viscous and inviscid flow fields is neglected. The solution method uses the finite difference form of the governing equations throughout the separated flow field and incorporates a transition model based on the integral turbulence kinetic energy equation. The validity of the numerical procedure is verified by making comparisons with analytical solutions to several test problems, including unsteady flow over a plate oscillating in its own plane. The method was then applied to the problem of unsteady viscous flow over a NACA 0012 airfoil oscillating sinusoidally in pitch. The flow field characteristics were in qualitative agreement with experimental results. The bubble moved forward on the airfoil and decreased in size as incidence was increased. Viscous flow in the leading edge region was found to be quasi-steady, while bubble height varied inversely with Reynolds number.

Kreskovsky, J. P.↗

Spatially random models, estimation theory, and robot arm dynamics

Spatially random models provide an alternative to the more traditional deterministic models used to describe robot arm dynamics. These alternative models can be used to establish a relationship between the methodologies of estimation theory and robot dynamics. A new class of algorithms for many of the fundamental robotics problems of inverse and forward dynamics, inverse kinematics, etc. can be developed that use computations typical in estimation theory. The algorithms make extensive use of the difference equations of Kalman filtering and Bryson-Frazier smoothing to conduct spatial recursions. The spatially random models are very easy to describe and are based on the assumption that all of the inertial (D'Alembert) forces in the system are represented by a spatially distributed white-noise model. The models can also be used to generate numerically the composite multibody system inertia matrix. This is done without resorting to the more common methods of deterministic modeling involving Lagrangian dynamics, Newton-Euler equations, etc. These methods make substantial use of human knowledge in derivation and minipulation of equations of motion for complex mechanical systems.

Rodriguez, G.↗

Numerical simulation of corotating interaction regions

We report on a numerical simulation of Corotating Interaction Regions (CIRs) in the solar wind. We employ an MHD calculation with a tilted dipole model: slow wind is assumed around the tilted current sheet. As a simplifying assumption the solar wind is kept radial, meridional and azimuthal deflections are neglected. Despite of this simplification the model reproduces several observed features of CIR-s. CIR-s turn out most prominent at mid-latitudes. Toward increasing heliographic latitudes, reverse shocks become dominant, forward shocks tend to become weaker then disappear. The frozen in magnetic field is calculated and applied for the simulation of the transport of cosmic rays in the heliosphere. Some consequences on recurrent cosmic ray phenomena and on the acceleration of energetic particles will be discussed.

Kota, J.↗

Investigation B: Laminar confined coaxial entrance flow with heat generation

The results of a parametric study on the entrance flow region in a gas core nuclear reactor are presented. The physical system is modeled as laminar confined, coaxial flow with heat generation in the inner fluid. The governing equations include the boundary layer approximations and the assumptions of only radial radiative transport of energy represented as an energy diffusion term. The Von Mises transformation and a zeta transformation are used to transform the equations into nonlinear nonhomogeneous convective-diffusion equations. A unique combination of forward and backward difference equations which yields accurate results at moderate computational times, is used in the numerical method. Results show that the rapidly accelerating, heat generating inner stream actually shrinks in radius as it expands axially.

Bobba, G. K. M.↗

A theoretical analysis of the extinction limits of a methane-air opposed-jet diffusion flame

A theoretical analysis is described for a methane-air diffusion flame stabilized in the forward stagnation region of a porous metal cylinder in a forced convective flow. The analysis includes effects of radiative heat loss from the porous metal surface and finite rate kinetics but neglects the effects of gravity. The theoretically predicted extinction limits compare well with experimentally observed extinction limits from the literature. After the predicted limits compared well with the experimental limits, a parametric study of the effect of fuel surface emissivity and Lewis number was conducted with the numerical model. It was found that the computed blowoff limit is independent of radiative heat loss for high fuel blowing velocities but is a strong function of Lewis number. At low fuel blowing velocities, the extinction limit varies with both radiative heat loss and Lewis number. It is discovered, however, that even if thermal losses from the fuel surface are absent, the flame can extinguish at the fuel surface independently of Lewis number due to excessive reaction zone thinning.

Olson, S. L.↗

Formulation and solution of rotary-wing aeroelastic stability and response problems

The state of the art in the formulation and solution of rotary-wing aeroelastic stability and response problems is reviewed in detail. The approximations used in the structural, inertia and aerodynamic operators are discussed. The important role of geometric nonlinearities, due to moderate deflections, and aerodynamic stall in the aeroelastic stability and response problem are identified. It is also shown that geometric nonlinearities are of primary importance in aeroelastic stability calculations, and have a more limited, though important, role in response calculations. Next, formulation of coupled rotor/fuselage problems is described, for both air and ground resonance type problems. Both topics, the isolated blade problem and the coupled rotor/fuselage problem, are treated for both hover and forward flight. Solution of aeroelastic stability and response problems proceeds in two stages. First, the spatial dependence is eliminated by using Galerkin's method, or by using the finite element method. Next the nonlinear, or linear, ordinary differential equation with periodic coefficients have to be solved for stability or response. Efficient numerical methods for accomplishing these objectives are presented in a comprehensive manner. The paper contains a number of illustrative numerical results which are intended to clarify various aspects of the modeling process and serve as representative results for both aeroelastic stability and response calculations for a variety of blade and rotor configurations.

Friedmann, P. P.↗

Numerical simulation of helicopter multi-bladed rotor flow

A numerical method for solving the Euler equations for multiblade rotors has been developed. The computed results, including one nonlifting forward-flight and several lifting hovering cases, are compared with experimental data. For forward-flight case, the unsteady growth and decay of the shock agree well with experimental results. For hovering case, at low and moderate transonic tip speed, the results show good agreement with experimental data in the tip region. The near-wake effects can be captured without wake modeling.

Chen, C. L.↗

Highly Variable Cycle Nozzle Concept: Validation of Flow and Noise Predictions

Results from experimental and numerical studies of highly Variable Cycle (HVC) exhaust model were presented. The model was designed and fabricated under a Supersonics NRA awarded to Rolls-Royce. The model had a lobed mixer for the core stream nozzle, and elliptic fan stream nozzle, and an ejector. Experiments included far-field acoustic array, phased array, and Particle Image Velocimetry (PIV) measurements. Numerical studies included flow simulations using the WIND-US code and far-field acoustic solutions using an acoustic analogy developed by Goldstein (2003) and Leib and Goldstein (2011). Far-field acoustic measurements showed increased noise levels over the round baseline nozzle when using non-static forward flight conditions. Phased array measurements showed noise sources near the ejector doors when tones were produced for small ejector door positions. Ejector door separation identified in the experiments was reproduced in the numerical flow simulations. Acoustic solutions were unable to match levels measured in the peak jet noise direction indicating additional development work is needed to predict noise from highly three-dimensional flows.

Halbig, Michael C.↗

The extinction properties of forest components

The effect of each forest component on the extinction of electromagnetic waves is investigated by modeling the branches with finite cylinders, deciduous leaves with elliptic disks, and coniferous leaves with needles. The inner field is estimated by the field inside an infinitely long cylinder of similar properties for the branches, and by the Shifrin approximation for the leaves. For each forest component analytic expressions were derived for the extinction cross section via the forward scattering theorem and for ohmic and scattered losses. For branches, the variation of the extinction cross section obtained via the forward scattering theorem is illustrated numerically as a function of the branch radius and the imaginery part of its dielectric constant. It is compared with the measurements from a single branch. For the leaves, the forward scattering theorem gives value for the extinction cross section equal to the ohmic cross section.

Karam, M. A.↗

A fully Galerkin method for the recovery of stiffness and damping parameters in Euler-Bernoulli beam models

A fully Sinc-Galerkin method for recovering the spatially varying stiffness and damping parameters in Euler-Bernoulli beam models is presented. The forward problems are discretized with a sinc basis in both the spatial and temporal domains thus yielding an approximate solution which converges exponentially and is valid on the infinite time interval. Hence the method avoids the time-stepping which is characteristic of many of the forward schemes which are used in parameter recovery algorithms. Tikhonov regularization is used to stabilize the resulting inverse problem, and the L-curve method for determining an appropriate value of the regularization parameter is briefly discussed. Numerical examples are given which demonstrate the applicability of the method for both individual and simultaneous recovery of the material parameters.

Smith, R. C.↗

The Impact Trajectory of Asteroid 2008 TC3

Asteroid 2008 TC3 was the rst asteroid ever discovered before reaching Earth. By using the almost 900 astrometric observations acquired prior to impact we estimate the trajectory of 2008 TC3 and the ground-track of the impact location as a function of altitude. For a reference altitude of 100 km the impact location 3- formal uncertainty is a 1.4 km 0.15 km ellipse with a semimajor axis azimuth of 105. We analyze the contribution of modeling errors and nd that the second-order zonal harmonics of the Earth gravity eld moves the ground-track by more than 1 km and the location along the ground-track by more than 2 km. Non-zonal and higher order harmonics only change the impact prediction by less than 20 m. The contribution of the atmospheric drag to the trajectory of 2008 TC3 is at the numerical integration error level, a few meters, down to an altitude of 50 km. Integrating forward to lower altitudes and ignoring the break-up of 2008 TC3, the atmospheric drag causes an along-track error that can be as large as a few kilometers at sea level. The locations of the recovered meteorites is consistent with the computed ground-track.

Trajectory↗

Multiple scattering of light in a spherical cometary atmosphere with an axisymmetric dust jet. II - Image simulation

A numerical solution for the multiple light scattering in spherical axisymmetric geometry is applied to the simulation of images of a coma as it would appear to a near-flying satellite such as Giotto. The appearance of symmetric comas and dust jets is examined in detail; the nucleus visibility is studied; the effect of forward scattering is considered; and single and multiple scattering effects are quantified. Attention is given to simulated images of a coma with a hollow cone of dust, as predicted by dust-gas hydrodynamic modeling. The cone's appearance is very similar to the northern area of activity on Comet Halley, observed by the Giotto HMC.

Chick, Kenneth M.↗

A fully Sinc-Galerkin method for Euler-Bernoulli beam models

A fully Sinc-Galerkin method in both space and time is presented for fourth-order time-dependent partial differential equations with fixed and cantilever boundary conditions. The Sinc discretizations for the second-order temporal problem and the fourth-order spatial problems are presented. Alternate formulations for variable parameter fourth-order problems are given which prove to be especially useful when applying the forward techniques to parameter recovery problems. The discrete system which corresponds to the time-dependent partial differential equations of interest are then formulated. Computational issues are discussed and a robust and efficient algorithm for solving the resulting matrix system is outlined. Numerical results which highlight the method are given for problems with both analytic and singular solutions as well as fixed and cantilever boundary conditions.

Smith, R. C.↗

Recent Developments on Airborne Forward Looking Interferometer for the Detection of Wake Vortices

A goal of these studies was development of the measurement methods and algorithms necessary to detect wake vortex hazards in real time from either an aircraft or ground-based hyperspectral Fourier Transform Spectrometer (FTS). This paper provides an update on research to model FTS detection of wake vortices. The Terminal Area Simulation System (TASS) was used to generate wake vortex fields of 3-D winds, temperature, and absolute humidity. These fields were input to the Line by Line Radiative Transfer Model (LBLRTM), a hyperspectral radiance model in the infrared, employed for the FTS numerical modeling. An initial set of cases has been analyzed to identify a wake vortex IR signature and signature sensitivities to various state variables. Results from the numerical modeling case studies will be presented. Preliminary results indicated that an imaging IR instrument sensitive to six narrow bands within the 670 to 3150 per centimeter spectral region would be sufficient for wake vortex detection. Noise floor estimates for a recommended instrument are a current research topic.

Daniels, Taumi S.↗

Numerical Simulation of Dual-Mode Scramjet Combustors

Results of a numerical investigation of a three-dimensional dual-mode scramjet isolator-combustor flow-field are presented. Specifically, the effect of wall cooling on upstream interaction and flow-structure is examined for a case assuming jet-to-jet symmetry within the combustor. Comparisons are made with available experimental wall pressures. The full half-duct for the isolator-combustor is then modeled in order to study the influence of side-walls. Large scale three-dimensionality is observed in the flow with massive separation forward on the side-walls of the duct. A brief review of convergence-acceleration techniques useful in dual-mode simulations is presented, followed by recommendations regarding the development of a reliable and unambiguous experimental data base for guiding CFD code assessments in this area.

Rodriguez, C. G.↗