SEARCH · Search NASA

Results for “Finite expression method”

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

All sources NASA NTRS DOE OSTI

NASA NTRS1982

A theoretical study of a two-wavelength lidar technique for the measurement of atmospheric temperature profiles

The theory of differential absorption lidar measurements for lines with a Voigt profile is given and applied to a two-wavelength technique for measuring the atmospheric temperature profile using a high J line in the oxygen A band. Explicit expressions for the temperature and pressure dependence of the absorption coefficient are developed for lines with a Voigt profile. An iteration procedure for calculating the temperature for narrow laser bandwidths is described which has an accuracy better than 0.2 K for bandwidths less than 0.01/cm. To reduce the errors in lidar measurements due to uncertainties in pressure, a method for estimating the pressure from the temperature profile is described. A procedure for extending the differential absorption technique to the case of finite laser bandwidth with good accuracy is also described. Simulation results show that a knowledge of the laser frequency is needed to 0.005/cm for accurate temperature measurements. Evaluation of the sensitivity for both ground- and Shuttle-based measurements shows accuracies generally better than 1 K. This technique allows up to an order of magnitude improvement in sensitivity compared to other differential absorption lidar techniques.

Korb, C. L.↗

NASA NTRS1987

Exact image theory for the problem of dielectric/magnetic slab

Exact image method, recently introduced for the exact solution of electromagnetic field problems involving homogeneous half spaces and microstrip-like geometries, is developed for the problem of homogeneous slab of dielectric and/or magnetic material in free space. Expressions for image sources, creating the exact reflected and transmitted fields, are given and their numerical evaluation is demonstrated. Nonradiating modes, guided by the slab and responsible for the loss of convergence of the image functions, are considered and extracted. The theory allows, for example, an analysis of finite ground planes in microstrip antenna structures.

Lindell, I. V.↗

NASA NTRS1987

Effect of current on spectrum of breaking waves in water of finite depth

This paper presents an approximate method to compute the mean value, the mean square value and the spectrum of waves in water of finite depth taking into account the effect of wave breaking with or without the presence of current. It is assumed that there exists a linear and Gaussian ideal wave train whose spectrum is first obtained using the wave energy flux balance equation without considering wave breaking. The Miche wave breaking criterion for waves in finite water depth is used to limit the wave elevation and establish an expression for the breaking wave elevation in terms of the elevation and its second time derivative of the ideal waves. Simple expressions for the mean value, the mean square value and the spectrum are obtained. These results are applied to the case in which a deep water unidirectional wave train, propagating normally towards a straight shoreline over gently varying sea bottom of parallel and straight contours, encounters an adverse steady current whose velocity is assumed to be uniformly distributed with depth. Numerical results are obtained and presented in graphical form.

Tung, C. C.↗

NASA NTRS2010

Repeatability and Reproducibility of Compression Strength Measurements Conducted According to ASTM E9

Ten commercial laboratories participated in an interlaboratory study to establish the repeatability and reproducibility of compression strength tests conducted according to ASTM International Standard Test Method E9. The test employed a cylindrical aluminum AA2024-T351 test specimen. Participants measured elastic modulus and 0.2 % offset yield strength, YS(0.2 % offset), using an extensometer attached to the specimen. The repeatability and reproducibility of the yield strength measurement, expressed as coefficient of variations were cv(sub r)= 0.011 and cv(sub R)= 0.020 The reproducibility of the test across the laboratories was among the best that has been reported for uniaxial tests. The reported data indicated that using diametrically opposed extensometers, instead of a single extensometer doubled the precision of the test method. Laboratories that did not lubricate the ends of the specimen measured yield stresses and elastic moduli that were smaller than those measured in laboratories that lubricated the specimen ends. A finite element analysis of the test specimen deformation for frictionless and perfect friction could not explain the discrepancy, however. The modulus measured from stress-strain data were reanalyzed using a technique that finds the optimal fit range, and applies several quality checks to the data. The error in modulus measurements from stress-strain curves generally increased as the fit range decreased to less than 40 % of the stress range.

Luecke, William E.↗

NASA NTRS1982

Nonlinear Green's function method for unsteady transonic flows

Advantages to employing Green's function in describing unsteady three-dimensional transonic flows are explored. The development of the function for application to linear subsonic and supersonic unsteady aerodynamics is reviewed. It is shown that unique solutions are possible for external flows, with all functional expressions being defined in Prandtl-Glauert space. The development of methods of using the Green's function for transonic flows is traced, noting the necessity of including the effects of significant nonlinear terms. The steady-state problem is considered to demonstrate the shock-capturing ability of the method and the usefulness of the function in the incompressible, subsonic, transonic, and supersonic areas of potential unsteady three-dimensional flows around complex configurations. Computational time is asserted to be an order of magnitude less than with finite difference methods.

Tseng, K.↗

NASA NTRS1984

On the Kelvin-Helmholtz instabilities of supersonic shear layers

Analytical consideration is given to a Kelvin-Helmholtz instability in a finite width supersonic shear layer with linear and sinusoidal velocity profiles in the transition region. An analogy is explored between the pressure perturbation and the bound and scattering wave states of the one-dimensional Schroedinger equation. A numerical model is obtained for generating contour maps of the real and imaginary parts of the perturbation frequency in the (A,B)-plane. A comparison is made between the numerical results and analytical expressions derived for very small Mach numbers. Two unstable regimes are characterized in the (A,B)-plane. The method is concluded useful for both laboratory and astrophysics problems.

Choudhury, S. R.↗

NASA NTRS1994

Spatial-Operator Algebra For Flexible-Link Manipulators

Method of computing dynamics of multiple-flexible-link robotic manipulators based on spatial-operator algebra, which originally applied to rigid-link manipulators. Aspects of spatial-operator-algebra approach described in several previous articles in NASA Tech Briefs-most recently "Robot Control Based on Spatial-Operator Algebra" (NPO-17918). In extension of spatial-operator algebra to manipulators with flexible links, each link represented by finite-element model: mass of flexible link apportioned among smaller, lumped-mass rigid bodies, coupling of motions expressed in terms of vibrational modes. This leads to operator expression for modal-mass matrix of link.

Jain, Abhinandan↗

NASA NTRS1977

A first-principles derivation of Doppler noise expected from solar wind density fluctuations

The level of Doppler noise (DN) expected from solar wind (SW) density fluctuations (DF) is derived beginning with the expression for refractive index variations. The calculation takes account of up and downlink paths and of the method actually used to produce the DN values. The usual assumptions that the DF are frozen in, that the large-scale radial variation can be separated from the DF, that the DF power spectrum is a power law with outer scale k sub 0 and that the DF are homogeneous on scales less than 2c delta t, delta t = sample time, are made. The result agrees quite well with previous observations of DN. Corrections for the finite number of points used in the actual algorithm are discussed.

Callahan, P. S.↗

NASA NTRS1994

Higher harmonic control analysis for vibration reduction of helicopter rotor systems

An advanced higher harmonic control (HHC) analysis has been developed and applied to investigate its effect on vibration reduction levels, blade and control system fatigue loads, rotor performance, and power requirements of servo-actuators. The analysis is based on a finite element method in space and time. A nonlinear time domain unsteady aerodynamic model, based on the indicial response formulation, is used to calculate the airloads. The rotor induced inflow is computed using a free wake model. The vehicle trim controls and blade steady responses are solved as one coupled solution using a modified Newton method. A linear frequency-domain quasi-steady transfer matrix is used to relate the harmonics of the vibratory hub loads to the harmonics of the HHC inputs. Optimal HHC is calculated from the minimization of the vibratory hub loads expressed in term of a quadratic performance index. Predicted vibratory hub shears are correlated with wind tunnel data. The fixed-gain HHC controller suppresses completely the vibratory hub shears for most of steady or quasi-steady flight conditions. HHC actuator amplitudes and power increase significantly at high forward speeds (above 100 knots). Due to the applied HHC, the blade torsional stresses and control loads are increased substantially. For flight conditions where the blades are stalled considerably, the HHC input-output model is quite nonlinear. For such cases, the adaptive-gain controller is effective in suppressing vibratory hub loads, even though HHC may actually increase stall areas on the rotor disk. The fixed-gain controller performs poorly for such flight conditions. Comparison study of different rotor systems indicates that a soft-inplane hingeless rotor requires less actuator power at high speeds (above 130 knots) than an articulated rotor, and a stiff-inplane hingeless rotor generally requires more actuator power than an articulated or a soft-inplane hingeless rotor. Parametric studies for a hingeless rotor operating in a transition flight regime and for an articulated rotor operating at the level-flight boundary (high speed and high thrust conditions) indicate that blade parameters including flap, lag, torsion stiffness distributions, linear pretwist, chordwise offset of center-of-mass from elastic axis and chordwise offset of elastic axis from aerodynamic center can be selected to minimize the actuator power requirements for HHC.

Nguyen, Khanh Q.↗

NASA NTRS2010

Methods for In-Flight Wing Shape Predictions of Highly Flexible Unmanned Aerial Vehicles: Formulation of Ko Displacement Theory

The Ko displacement theory is formulated for a cantilever tubular wing spar under bending, torsion, and combined bending and torsion loading. The Ko displacement equations are expressed in terms of strains measured at multiple sensing stations equally spaced on the surface of the wing spar. The bending and distortion strain data can then be input to the displacement equations to calculate slopes, deflections, and cross-sectional twist angles of the wing spar at the strain-sensing stations for generating the deformed shapes of flexible aircraft wing spars. The displacement equations have been successfully validated for accuracy by finite-element analysis. The Ko displacement theory that has been formulated could also be applied to calculate the deformed shape of simple and tapered beams, plates, and tapered cantilever wing boxes. The Ko displacement theory and associated strain-sensing system (such as fiber optic sensors) form a powerful tool for in-flight deformation monitoring of flexible wings and tails, such as those often employed on unmanned aerial vehicles. Ultimately, the calculated displacement data can be visually displayed in real time to the ground-based pilot for monitoring the deformed shape of unmanned aerial vehicles during flight.

Ko, William L.↗

NASA NTRS2007

Displacement Theories for In-Flight Deformed Shape Predictions of Aerospace Structures

Displacement theories are developed for a variety of structures with the goal of providing real-time shape predictions for aerospace vehicles during flight. These theories are initially developed for a cantilever beam to predict the deformed shapes of the Helios flying wing. The main structural configuration of the Helios wing is a cantilever wing tubular spar subjected to bending, torsion, and combined bending and torsion loading. The displacement equations that are formulated are expressed in terms of strains measured at multiple sensing stations equally spaced on the surface of the wing spar. Displacement theories for other structures, such as tapered cantilever beams, two-point supported beams, wing boxes, and plates also are developed. The accuracy of the displacement theories is successfully validated by finite-element analysis and classical beam theory using input-strains generated by finite-element analysis. The displacement equations and associated strain-sensing system (such as fiber optic sensors) create a powerful means for in-flight deformation monitoring of aerospace structures. This method serves multiple purposes for structural shape sensing, loads monitoring, and structural health monitoring. Ultimately, the calculated displacement data can be visually displayed to the ground-based pilot or used as input to the control system to actively control the shape of structures during flight.

Ko, William L.↗

NASA NTRS2001

Highly parallel structured adaptive mesh refinement using parallel language-based approaches

Adaptive mesh refinement (AMR) calculations carried out on structured meshes play an exceedingly important role in several areas of science and engineering. A strategy for using Fortran 90 in an object-oriented fashion is presented. This permits AMR applications to be expressed in terms of familiar abstractions that are natural to the process of solving AMR hierarchies. The OpenMP features that are useful for parallel processing of AMR hierarchies in a load balanced fashion on multiprocessors is described.

computational↗

NASA NTRS1990

Derivation of generalized transition/boundary conditions for planar multiple-layer structures

Infinite-order generalized impedance boundary conditions (GIBCs) and generalized sheet transition conditions (GSTCs) for planar multilayer configurations are developed via the Taylor series expansion method. The conditions are derived in a matrix product form where each matrix corresponds to a specific layer. An overall composite boundary/transition condition is obtained by making finite-order approximations to the elements of each matrix for the cases of 'low-contrast' and 'high-contrast' material layers. The accuracy of the truncated boundary conditions is examined by comparing their implied reflection and transmission coefficients with the corresponding exact coefficients. Design curves are also given which relate the maximum order of the conditions required to simulate a coating or layer of specific thickness and contrast. Expressions are then derived for the reflection and transmission coefficients of the GIBC/GSTC sheets, and these are compared to exact coefficients to demonstrate the validity of the derived GIBCs/GSTCs.

Ricoy, M. A.↗

NASA NTRS2015

Numerical Computation of a Continuous-thrust State Transition Matrix Incorporating Accurate Hardware and Ephemeris Models

A significant body of work exists showing that providing a nonlinear programming (NLP) solver with expressions for the problem constraint gradient substantially increases the speed of program execution and can also improve the robustness of convergence, especially for local optimizers. Calculation of these derivatives is often accomplished through the computation of spacecraft's state transition matrix (STM). If the two-body gravitational model is employed as is often done in the context of preliminary design, closed form expressions for these derivatives may be provided. If a high fidelity dynamics model, that might include perturbing forces such as the gravitational effect from multiple third bodies and solar radiation pressure is used then these STM's must be computed numerically. We present a method for the power hardward model and a full ephemeris model. An adaptive-step embedded eight order Dormand-Prince numerical integrator is discussed and a method for the computation of the time of flight derivatives in this framework is presented. The use of these numerically calculated derivatieves offer a substantial improvement over finite differencing in the context of a global optimizer. Specifically the inclusion of these STM's into the low thrust missiondesign tool chain in use at NASA Goddard Spaceflight Center allows for an increased preliminary mission design cadence.

Low Thrust↗

NASA NTRS1978

An efficient, direct finite difference method for computing sound propagation in arbitrarily shaped two-dimensional and axisymmetric ducts without flow

An efficient, direct finite difference method is presented for computing sound propagation in non-stepped two-dimensional and axisymmetric ducts of arbitrarily varying cross section without mean flow. The method is not restricted by axial variation of acoustic impedance of the duct wall linings. The non-uniform two-dimensional or axisymmetric duct is conformally mapped numerically into a rectangular or cylindrical computational domain using a new procedure based on a method of fast direct solution of the Cauchy-Riemann equations. The resulting Helmholtz equation in the computational domain is separable. The solution to the governing equation and boundary conditions is expressed as a linear combination of fundamental solutions. The fundamental solutions are computed only once for each duct shape by means of the fast direct cyclic reduction method for the discrete solution of separable elliptic equations. Numerical results for several examples are presented to show the applicability and efficiency of the method.

Chakravarthy, S.↗

NASA NTRS1991

Eigenvalue sensitivity analysis of planar frames with variable joint and support locations

Two sensitivity equations are derived in this study based upon the continuum approach for eigenvalue sensitivity analysis of planar frame structures with variable joint and support locations. A variational form of an eigenvalue equation is first derived in which all of the quantities are expressed in the local coordinate system attached to each member. Material derivative of this variational equation is then sought to account for changes in member's length and orientation resulting form the perturbation of joint and support locations. Finally, eigenvalue sensitivity equations are formulated in either domain quantities (by the domain method) or boundary quantities (by the boundary method). It is concluded that the sensitivity equation derived by the boundary method is more efficient in computation but less accurate than that of the domain method. Nevertheless, both of them in terms of computational efficiency are superior to the conventional direct differentiation method and the finite difference method.

Chuang, Ching H.↗

NASA NTRS1991

Analytical and experimental study of vibrations in a gear transmission

An analytical simulation of the dynamics of a gear transmission system is presented and compared to experimental results from a gear noise test rig at the NASA Lewis Research Center. The analytical procedure developed couples the dynamic behaviors of the rotor-bearing-gear system with the response of the gearbox structure. The modal synthesis method is used in solving the overall dynamics of the system. Locally each rotor-gear stage is modeled as an individual rotor-bearing system using the matrix transfer technique. The dynamics of each individual rotor are coupled with other rotor stages through the nonlinear gear mesh forces and with the gearbox structure through bearing support systems. The modal characteristics of the gearbox structure are evaluated using the finite element procedure. A variable time steping integration routine is used to calculate the overall time transient behavior of the system in modal coordinates. The global dynamic behavior of the system is expressed in a generalized coordinate system. Transient and steady state vibrations of the gearbox system are presented in the time and frequency domains. The vibration characteristics of a simple single mesh gear noise test rig is modeled. The numerical simulations are compared to experimental data measured under typical operating conditions. The comparison of system natural frequencies, peak vibration amplitudes, and gear mesh frequencies are generally in good agreement.

Choy, F. K.↗

NASA NTRS1979

Numerical results for the thermal scattering functions

A recent formulation in radiative transfer defined the thermal scattering functions that characterize radiative transfer from a general plane-parallel finite medium driven solely by an internal distribution of thermal sources. Exiting diffuse intensities are expressed as space convolutions of the thermal scattering functions with any thermal source distribution. A parametric study is presented to obtain the basic structure of these scattering functions. The independent variables of these azimuthally independent functions are the direction cosine and source location, while the parameters are the single-scattering albedo, total optical depth, and the asymmetry factor in the Henyey-Greenstein phase function. The basic functional trends are discussed by using various parametric plots, and selected results are given to allow numerical checks. The computational method is invariant imbedding.

Cogley, A. C.↗