SEARCH · Engineering Papers

Results for “Numerical approximations”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,261 records · Page 70

All sources NASA NTRS DOE OSTI

NASA NTRS1990

Elimination of Gibbs' phenomena from error analysis of finite element results

This paper is one of a series on error analysis and correction of finite element solutions for plates and shells. The error analysis in the earlier papers used half-range double Fourier sine series for numerical harmonic analysis. The half-range formulas are simple to apply, but they can be inaccurate near the ends of the ranges of the independent variables. The Gibbs' phenomenon exhibited by half-range sine series in one independent variable has a two-dimensional analog; a classic example is the Navier solution in a double half-range sine series for the simply supported plate under a uniform load. A simple change of variables is introduced in the paper to improve the accuracy of the double Fourier sine series without adding complexity to the numerical analysis. The change of variables is applied to the problem of approximating a transverse load that is tabulated on a rectangular grid. A solution based on the change of variables is compared with results from the Navier solution for the simply supported plate problem and finite element results for the same problem.

Thurston, Gaylen A.↗

NASA NTRS1990

Resonance absorption of propagating fast waves in a cold plasma

Absorption of propagating waves impinging on a surface in which the plasma and magnetic field may change is investigated by examining in depth the problem of a combination of cold plasma, uniform magnetic field and a surface density which varies linearly from zero at the left end to some finite value at the right end, beyond which the density is constant. Two cases are considered: one in which the plasma is a vacuum everywhere to the left of the surface (which may correspond to coronal conditions) and one in which the plasma density jumps to a very large value to the left of the surface (which may mimic the magnetosphere with the dense region at the left corresponding to the plasmasphere). A complete discussion of the resonance absorption of propagating fast waves for the case considered by Kiveloson and Southwood (1986) is presented, emphasizing approximate analytical results whenever possible; these results are then compared with exact numerical solutions.

Hollweg, Joseph V.↗

NASA NTRS1991

Estimation of optical flow in airborne electro-optical sensors by stochastic approximation

The essence of motion or range estimation by passive electrooptical means is the ability to determine the correspondence of picture elements in pairs of image frames and to estimate their coordinates and their disparity (relative shifts) in the image plane of an electrooptical imaging sensor. The disparity can be in successive frames due to self-motion or in simultaneous frames of a stereo pair. A key issue is to provide these estimates on-line. This paper describes the theoretical background of such an interframe shift estimator. It is based on a stochastic gradient algorithm, specifically implementing a form of stochastic approximation, which can achieve rapid convergence of the shift estimate. Analytical and numerical simulation examples for random texture and isolated features validate the feasibility and the effectiveness of the estimator.

Merhav, S. J.↗

NASA NTRS1993

Recent developments in equivalent plate modeling for wing shape optimization

A new technique for structural modeling of airplane wings is presented taking transverse shear effects into account. The kinematic assumptions of first order shear deformation plate theory in combination with numerical analysis based on simple polynomials which define geometry, construction and displacement approximations lead to analytical expressions for elements of the stiffness and mass matrices and load vector. Contributions from the cover skins, spar and rib caps and spar and rib webs are included as well as concentrated springs and concentrated masses. Limitations of current equivalent plate wing modeling techniques based on classical plate theory are discussed, and the improved accuracy of the new equivalent plate technique is demonstrated through comparison to finite element analysis and test results. Analytical derivatives of stiffness, mass and load terms with respect to wing shape lead to analytic sensitivities of displacements, stresses and natural modes with respect to planform shape and depth distribution. This makes the new capability an effective structural tool for wing shape optimization.

Livne, Eli↗

NASA NTRS1992

Emissions of sulfur gases from marine and freshwater wetlands of the Florida Everglades: Rates and extrapolation using remote sensing

Rates of emissions of the biogenic sulfur (S) gases carbonyl sulfide (COS), methyl mercaptan (MSH), dimethyl sulfide (DMS), and carbon disulfide (CS2) were measured in a variety of marine and freshwater wetland habitats in the Florida Everglades during a short duration period in October using dynamic chambers, cryotrapping techniques, and gas chromatography. The most rapid emissions of greater than 500 nmol/m(sup -2)h(sup -1) occurred in red mangrove-dominated sites that were adjacent to open seawater and contained numerous crab burrows. Poorly drained red mangrove sites exhibited lower fluxes of approximately 60 nmol/m(sup -2)h(sup -1) which were similar to fluxes from the black mangrove areas which dominated the marine-influenced wetland sites in the Everglades. DMS was the dominant organo-S gas emitted especially in the freshwater areas. Spectral data from a scene from the Landsat thematic mapper were used to map habitats in the Everglades. Six vegetation categories were delineated using geographical information system software and S gas emission were extrapolated for the entire Everglades National Park. The black mangrove-dominated areas accounted for the largest portion of S gas emissions to the area. The large area extent of the saw grass communities (42 percent) accounted for approximately 24 percent of the total S emissions.

Hines, Mark E.↗

NASA NTRS1981

Free-Vibration Analysis of Rotating Beams by a Variable-Order Finite-Element Method

The free vibration of rotating beams is analyzed by means of a finite-element method of variable order. This method entails displacement functions that are a complete power series of a variable number of terms. The terms are arranged so that the generalized coordinates are composed of displacements and slopes at the element extremities and, additionally, displacements at certain points within the element. The displacement is assumed to be analytic within an element and thus can be approximated to any degree of accuracy desired by a complete power series. Numerical results are presented for uniform beams with zero and nonzero hub radii, tapered beams, and a nonuniform beam with discontinuities. Since the present method reduces to a conventional beam finite-element method for a cubic displacement function, the results are compared and found to be superior to the conventional results in terms of accuracy for a given number of degrees of freedom. Indeed, essentially exact eigenvalues and eigenvectors are obtained with this technique, which is far more rapidly convergent than other approaches in the literature.

Hodges, Dewey H.↗

NASA NTRS1998

Preliminary Work for Modeling the Propellers of an Aircraft as a Noise Source in an Acoustic Boundary Element Analysis

An algorithm for generating appropriate velocity boundary conditions for an acoustic boundary element analysis from the kinematics of an operating propeller is presented. It constitutes the initial phase of Integrating sophisticated rotorcraft models into a conventional boundary element analysis. Currently, the pressure field is computed by a linear approximation. An initial validation of the developed process was performed by comparing numerical results to test data for the external acoustic pressure on the surface of a tilt-rotor aircraft for one flight condition.

Vlahopoulos, Nickolas↗

NASA NTRS1999

Finite Differences and Collocation Methods for the Solution of the Two Dimensional Heat Equation

In this paper we combine finite difference approximations (for spatial derivatives) and collocation techniques (for the time component) to numerically solve the two dimensional heat equation. We employ respectively a second-order and a fourth-order schemes for the spatial derivatives and the discretization method gives rise to a linear system of equations. We show that the matrix of the system is non-singular. Numerical experiments carried out on serial computers, show the unconditional stability of the proposed method and the high accuracy achieved by the fourth-order scheme.

Kouatchou, Jules↗

NASA NTRS1995

An Investigation of Multipath Effects on the GPS System During Auto-Rendezvous and Capture

The proposed use of a Cargo Transport Vehicle (CTV) to carry hardware to the Space Station Freedom (SSF) during the construction phase of the SSF project requires remote maneuvering of the CTV. The CTV is not a manned vehicle. Obtaining the relative positions of the CTV and SSF for remote auto-rendezvous and capture (AR&C) scenarios will rely heavily on the Global Positioning System (GPS). The GPS system is expected to guide the CTV up to a distance of 100 to 300 meters from the SSF. At some point within this range, an optical docking system will take over the remote guidance for capture. During any remote guidance by GPS it is possible that significant multipath signals may be caused by large objects in the vicinity of the module being remotely guided. This could alter the position obtained by the GPS system from the actual position. Due to the nature of the GPS signals, it has been estimated that if the difference in distance between the Line of Sight (LOS) path and the multipath is greater than 300 meters, the GPS system is capable of discriminating between the direct signal and the reflected (or multipath) signal. However, if the path difference is less than 300 meters, one must be concerned. This report details the work accomplished by the Electromagnetic Simulations Laboratory at Marquette University over the period December 1993 to May 1995. This work is an investigation of the strength and phase of a multipath signal arriving at the CTV relative to the direct or line of sight (LOS) signal. The signal originates at a GPS satellite in half geo-stationary orbit and takes two paths to the CTV: (1) the direct or LOS path from the GPS satellite to the CTV; and (2) a scattered path from the GPS satellite to the SSF module and then to the CTV. The scattering from a cylinder has been computed using the physical optics approximation for the current. No other approximations or assumptions have been made including no assumptions regarding the far field or Fresnel field approximations. The integrations required to obtain the scattered field have been computed numerically using an N dimensional Romberg integration. The total scattered electric field is then projected onto the RCP component in the direction of propagation only. The direct or line of sight signal is then used to compute the relative strength and phase of the scattered field. The trajectory of the CTV has been parameterized into 4,214 points that are calculated for each of the geometries investigated. The motion of the CTV between points is small enough for the magnitude data (dB down from direct signal) to appear very smooth; however, because of the distances and wavelengths involved, the phase of the scattered field relative to the direct signal varies very rapidly.

Richie, James E.↗

NASA NTRS2000

Exchange Forces in Dispersion Relations Investigated Using Circuit Relations

We propose a novel method to compute in an exact manner the left-hand cut discontinuity of the electron-Atom partial wave scattering amplitude in the complex energy plane within the static stage approximation. Zero energy dispersion relations for electron-Hydrogen scattering are computed numerically for illustration.

Vrinceanu, D.↗

NASA NTRS2000

Accuracy Improvement in Magnetic Field Modeling for an Axisymmetric Electromagnet

This paper examines the accuracy and calculation speed for the magnetic field computation in an axisymmetric electromagnet. Different numerical techniques, based on an adaptive nonuniform grid, high order finite difference approximations, and semi-analitical calculation of boundary conditions are considered. These techniques are being applied to the modeling of the Variable Specific Impulse Magnetoplasma Rocket. For high-accuracy calculations, a fourth-order scheme offers dramatic advantages over a second order scheme. For complex physical configurations of interest in plasma propulsion, a second-order scheme with nonuniform mesh gives the best results. Also, the relative advantages of various methods are described when the speed of computation is an important consideration.

Ilin, Andrew V.↗

NASA NTRS2000

Global MHD Simulations of Space Plasma Environments: Heliosphere, Comets, Magnetospheres of Plants and Satellites

Magnetohydrodynamics (MHD) provides an approximate description of a great variety of processes in space physics. Accurate numerical solutions of the MHD equations are still a challenge, but in the past decade a number of robust methods have appeared. Once these techniques made the direct solution of MHD equations feasible, a number of global three-dimensional models were designed and applied to many space physics objects. The range of these objects is truly astonishing, including active galactic nuclei, the heliosphere, the solar corona, and the solar wind interaction with planets, satellites, and comets. Outside the realm of space physics, MHD theory has been applied to such diverse problems as laboratory plasmas and electromagnetic casting of liquid metals. In this paper we present a broad spectrum of models of different phenomena in space science developed in the recent years at the University of Michigan. Although the physical systems addressed by these models are different, they all use the MHD equations as a unifying basis.

Kabin, K.↗

NASA NTRS2005

Aspects of Numerical Simulation of Circulation Control Airfoils

The mass-averaged compressible Navier-Stokes equations are solved for circulation control airfoils. Numerical solutions are computed with a multigrid method that uses an implicit approximate factorization smoother. The effects of flow conditions (e.g., free-stream Mach number, angle of attack, momentum coefficient) and mesh on the prediction of circulation control airfoil flows are considered. In addition, the impact of turbulence modeling, including curvature effects and modifications to reduce eddy viscosity levels in the wall jet (i.e., Coanda flow), is discussed. Computed pressure distributions are compared with available experimental data.

Swanson, R. C.↗

NASA NTRS2006

A Historical Search for the Occurrence of Habitable Ground Ice at the Phoenix Landing Site

A numerical model of the thermal history of Martian ground ice at the approximate location of the planned Phoenix landing site has been developed and used to identify instances of relatively warm ground ice over the last 10 Ma. Many terrestrial organisms are adapted to life at or below the freezing temperature of water, and we will use the approximate doubling time of terrestrial microbial populations as a function of temperature, is used as a metric against which to assess the "habitability" of Martian ground ice.

Zent, Aaron P.↗

NASA NTRS2008

Unsteady Analysis of Blade and Tip Heat Transfer as Influenced by the Upstream Momentum and Thermal Wakes

The effect of the upstream wake on the blade heat transfer has been numerically examined. The geometry and the flow conditions of the first stage turbine blade of GE s E3 engine with a tip clearance equal to 2 percent of the span was utilized. Based on numerical calculations of the vane, a set of wake boundary conditions were approximated, which were subsequently imposed upon the downstream blade. This set consisted of the momentum and thermal wakes as well as the variation in modeled turbulence quantities of turbulence intensity and the length scale. Using a one-blade periodic domain, the distributions of unsteady heat transfer rate on the turbine blade and its tip, as affected by the wake, were determined. Such heat transfer coefficient distribution was computed using the wall heat flux and the adiabatic wall temperature to desensitize the heat transfer coefficient to the wall temperature. For the determination of the wall heat flux and the adiabatic wall temperatures, two sets of computations were required. The results were used in a phase-locked manner to compute the unsteady or steady heat transfer coefficients. It has been found that the unsteady wake has some effect on the distribution of the time averaged heat transfer coefficient on the blade and that this distribution is different from the distribution that is obtainable from a steady computation. This difference was found to be as large as 20 percent of the average heat transfer on the blade surface. On the tip surface, this difference is comparatively smaller and can be as large as four percent of the average.

Ameri, Ali A.↗

NASA NTRS2008

Direct Solve of Electrically Large Integral Equations for Problem Sizes to 1M Unknowns

Matrix methods for solving integral equations via direct solve LU factorization are presently limited to weeks to months of very expensive supercomputer time for problems sizes of several hundred thousand unknowns. This report presents matrix LU factor solutions for electromagnetic scattering problems for problem sizes to one million unknowns with thousands of right hand sides that run in mere days on PC level hardware. This EM solution is accomplished by utilizing the numerical low rank nature of spatially blocked unknowns using the Adaptive Cross Approximation for compressing the rank deficient blocks of the system Z matrix, the L and U factors, the right hand side forcing function and the final current solution. This compressed matrix solution is applied to a frequency domain EM solution of Maxwell's equations using standard Method of Moments approach. Compressed matrix storage and operations count leads to orders of magnitude reduction in memory and run time.

Shaeffer, John↗

NASA NTRS2011

AGN Obscuration Through Dusty Infrared Dominated Flows: Radiation-Hydrodynamics Solution for the Wind - 1

We construct a radiation-hydrodynamics model for the obscuring toroidal structure in active galactic nuclei. In this model the obscuration is produced at parsec scale by a dense, dusty wind which is supported by infrared radiation pressure on dust grains. To find the distribution of radiation pressure, we numerically solve the 2D radiation transfer problem in a flux limited diffusion approximation. We iteratively couple the solution with calculations of stationary 1D models for the wind, and obtain the z-component of the velocity. Our results demonstrate that for AGN luminosities greater than 0.1 L(sub edd) external illumination can support a geometrically thick obscuration via outflows driven by infrared radiation pressure. The terminal velocity of marginally Compton-thin models (0.2 < tau(sub T) < 0.6), is comparable to or greater than the escape velocity. In Compton thick models the maximum value of the vertical component of the velocity is lower than the escape velocity, suggesting that a significant part of our torus is in the form of failed wind. The results demonstrate that obscuration via normal or failed infrared-driven winds is a viable option for the AGN torus problem and AGN unification models. Such winds can also provide an important channel for AGN feedback.

Dorodnitsyn, A.↗

NASA NTRS2012

Infrared Dielectric Properties of Low-stress Silicon Nitride

Silicon nitride thin films play an important role in the realization of sensors, filters, and high-performance circuits. Estimates of the dielectric function in the far- and mid-IR regime are derived from the observed transmittance spectra for a commonly employed low-stress silicon nitride formulation. The experimental, modeling, and numerical methods used to extract the dielectric parameters with an accuracy of approximately 4% are presented.

dielectric parameters↗