Search NASA⌕ Search

SEARCH · Search NASA

Results for “Higher order calculations”

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 181 records · Page 10
All sourcesNASA NTRSDOE OSTI
NASA NTRS1991

The propagation characteristics of the plate modes of acoustic emission waves in thin aluminum plates and thin graphite/epoxy composite plates and tubes

Acoustic emission was interpreted as modes of vibration in plates. Classical plate theory was used to predict dispersion curves for the two fundamental modes and to calculate the shapes of flexural waveforms produced by vertical step function loading. There was good agreement between theoretical and experimental results for aluminum. Composite materials required the use of a higher order plate theory (Reissner-Mindlin) to get good agreement with the measured velocities. Four composite plates with different laminate stacking sequences were studied. The dispersion curves were determined from phase spectra of the time dependent waveforms. Plate modes were shown to be useful for determining the direction of source motion. Aluminum plates were loaded by breaking a pencil lead against their surface. By machining slots at angles to the plane of a plate, the direction in which the force acted was varied. Changing the source motion direction produced regular variations in the waveforms. To demonstrate applicability beyond simple plates, waveforms produced by lead breaks on a thin walled composite tube were also shown to be interpretable as plate modes. The tube design was based on the type of struts proposed for Space Station Freedom's trussed structures.

Prosser, William H.↗
NASA NTRS2009

Measuring Massive Black Hole Binaries with LISA

The coalescence of two massive black holes produces gravitational waves (GWs) which can be detected by the space-based detector LISA. By measuring these waves, LISA can determine the various parameters which characterize the source. Measurements of the black hole masses and spins will provide information about the growth of black holes and their host galaxies over time. Measurements of a source's sky position and distance may help astronomers identify an electromagnetic counterpart to the GW event. The counterpart's redshift, combined with the GW-measured luminosity distance, can then be used to measure the Hubble constant and the dark energy parameter $w$. Because the potential science output is so high, it is useful to know in advance how well LISA can measure source parameters for a wide range of binaries. We calculate expected parameter estimation errors using the well-known Fisher matrix method. Our waveform model includes the physics of spin precession, as well as subleading harmonics. When these higher-order effects are not included, strong degeneracies between some parameters cause them to be poorly determined by a GW measurement. When precession and subleading harmonics are properly included, the degeneracies are broken, reducing parameter errors by one to several orders of magnitude.

Lang, Ryan N.↗
NASA NTRS1980

Higher order farfield drag minimization for a subcritical wing design code

A higher order Trefftz plane model of the undistorted interacting wing wakes of two symmetric subsonic planforms has been developed and integrated with an existing vortex lattice wing design code. The Trefftz plane calculation provides continuous, piecewise quadratically varying bound circulations, which are then interpolated in the nearfield vortex lattice representation of the wing. Integration of the surface slopes for a specified chord loading function provides the design wing camber surfaces. This paper outlines the theoretical development and presents comparisons between results obtained using the present theory and (1) previous exact solutions for minimum induced drag and (2) previous wing designs obtained using a discrete vortex wake moded.

Kuhlman, J. M.↗
NASA NTRS1981

The dynamics and control of large flexible space structures-IV

The effects of solar radiation pressure as the main environmental disturbance torque were incorporated into the model of the rigid orbiting shallow shell and computer simulation results indicate that within the linear range the rigid modal amplitudes are excited in proportion to the area to mass ratio. The effect of higher order terms in the gravity-gradient torque expressions previously neglected was evaluated and found to be negligible for the size structures under consideration. A graph theory approach was employed for calculating the eigenvalues of a large flexible system by reducing the system (stiffness) matrix to lower ordered submatrices. The related reachability matrix and term rank concepts are used to verify controllability and can be more effective than the alternate numerical rank tests. Control laws were developed for the shape and orientation control of the orbiting flexible shallow shell and numerical results presented.

Bainum, P. M.↗
NASA NTRS1984

The application of a second generation low-order panel method - program 'Vsaero' - to powerplant installation studies

Results from the application of a second generation low-order panel method to a number of powerplant installations are presented. The method, which retains the advantages of ease of use and low operating cost of earlier panel methods while retaining the aerodynamic modelling rigor of higher-order methods was used to predict the flow in and around typical nacelles and aircraft for a wide range of conditions. The ability of the second generation low-order analysis to calculate internal flows without the leakage problems associated with earlier programs is demonstrated and correlation between calculated and measured surface pressures and flow behavior, including predicting streamlines and locating regions of separated flow, is presented.

Clark, D. R.↗
NASA NTRS1989

Nonlinear interactions in mixing layers and compressible heated round jets

The nonlinear interactions between a fundamental instability mode and both its harmonics and the changing mean flow are studied using the weakly nonlinear stability theory of Stuart and Watson, and numerical solutions of coupled nonlinear partial differential equations. The first part focuses on incompressible cold (or isothermal; constant temperature throughout) mixing layers, and for these, the first and second Landau constants are calculated as functions of wavenumber and Reynolds number. It is found that the dominant contribution to the Landau constants arises from the mean flow changes and not from the higher harmonics. In order to establish the range of validity of the weakly nonlinear theory, the weakly nonlinear and numerical solutions are compared and the limitation of each is discussed. At small amplitudes and at low-to-moderate Reynolds numbers, the two results compare well in describing the saturation of the fundamental, the distortion of the mean flow, and the initial stages of vorticity roll-up. At larger amplitudes, the interaction between the fundamental, second harmonic, and the mean flow is strongly nonlinear and the numerical solution predicts flow oscillations, whereas the weakly nonlinear theory yields saturation. In the second part, the weakly nonlinear theory is extended to heated (or nonisothermal; mean temperature distribution) subsonic round jets where quadratic and cubic nonlinear interactions are present, and the Landau constants also depend on jet temperature ratio, Mach number and azimuthal mode number. Under exponential growth and nonlinear saturation, it is found that heating and compressibility suppress the growth of instability waves, that the first azimuthal mode is the dominant instability mode, and that the weakly nonlinear solution describes the early stages of the roll-up of an axisymmetric shear layer. The receptivity of a typical jet flow to pulse type input disturbance is also studied by solving the initial value problem and then examining the behavior of the long-time solution.

Jarrah, Yousef Mohd↗
NASA NTRS1992

A high-frequency ray analysis of the electromagnetic backscattering by a two-dimensional ogive

The electromagnetic backscattering from a two-dimensional ogive is examined, including a mechanism-by-mechanism account for the scattered field. The uniform theory of diffraction is used in conjunction with a simple creeping wave representation to yield a solution which provides an accurate and continuous backscattering result for a wide range of ogive geometry parameters. Although not a rigorous derivation of higher-order diffraction ceofficients, this solution provides insight into the relative contributions of the various creeping waves which are generated and propagate around the two-dimensional ogive. Comparisons with moment method calculations are used to demonstrate the range of applicability of this high-frequency asymptotic solution.

Shamansky, H.↗
NASA NTRS1979

Interaction between a normal shock wave and a turbulent boundary layer at high transonic speeds. Part 1: Pressure distribution

Analytical solutions are derived which incorporate additional physical effects as higher order terms for the case when the sonic line is very close to the wall. The functional form used for the undisturbed velocity profile is described to indicate how various parameters will be calculated for later comparison with experiment. The basic solutions for the pressure distribution are derived. Corrections are added for flow along a wall having longitudinal curvature and for flow in a circular pipe, and comparisons with available experimental data are shown.

Messiter, A. F.↗
NASA NTRS2021

Temperature Fluctuations and Boundary Layer Turbulence as seen by Mars Exploration Rovers Miniature Thermal Emission Spectrometer

The Mini-TES spectrometer carried by Spirit and Opportunity is a unique instrument that enabled vertical temperature profiles to be retrieved at high frequency (every 2 s, or “ICK”). Observations obtained from this instrument provide information on boundary layer activity. As an extension of previous work, we have retrieved ICK-by-ICK vertical temperature profiles to characterize turbulent behavior at each rover site for up to two Mars Years. We have also retrieved ICK-by-ICK near-surface (1.1 m) atmospheric temperatures and surface temperatures from dedicated atmospheric observations as well as geologic observations that extend the coverage of near-surface observations. The extension and increased frequency of retrievals as well as the additional observations provide a more complete view of the vertical and temporal variation of temperature fluctuations within the PBL. Temperature fluctuations are assessed using detrended temperature time-series data. The results show fluctuations can be as large as 10 K at the surface and 8 K aloft with timescales on the order of minutes. We have calculated the peak-to-peak amplitudes of temperature fluctuations and the corresponding timescales over which they occur. Spirit and Opportunity temperature fluctuation amplitudes decrease in magnitude vertically, where higher atmospheric levels show smaller fluctuations in amplitudes than at the surface. The magnitude of temperature fluctuations changes with the difference in surface and air temperature. This difference varies little as a function of season at both rover sites. Fluctuations are observed to quickly decrease in the late afternoon reaching a minimum near 17:00 LTST followed by moderately increased activity near sunset caused by forced convection. There is no significant change in the observed convective activity as a function of dust optical depth during regional dust storms. The few large temperature fluctuations that are concurrent with dust events are likely the result of large-scale systems.

Mars, atmosphere↗
NASA NTRS1957

Comparison of several methods for obtaining the time response of linear systems to either a unit impulse or arbitrary input from frequency-response data

Several methods of obtaining the time response of Linear systems to either a unit impulse or an arbitrary input from frequency-response data are described and compared. Comparisons indicate that all the methods give good accuracy when applied to a second-order system; the main difference is the required computing time. The methods generally classified as inverse Laplace transform methods were found to be most effective in determining the response to a unit impulse from frequency-response data of higher order systems. Some discussion and examples are given of the use of such methods as flight-data-analysis techniques in predicting loads and motions of a flexible aircraft on the basis of simple calculations when the aircraft frequency response is known.

Donegan, James J↗
NASA NTRS2017

Metallic Junction Thermoelectric Device Simulations

Thermoelectric junctions made of semiconductors have existed in radioisotope thermoelectric generators (RTG) for deep space missions, but are currently being adapted for terrestrial energy harvesting. Unfortunately, these devices are inefficient, operating at only 7% efficiency. This low efficiency has driven efforts to make high-figure-of-merit thermoelectric devices, which require a high electrical conductivity but a low thermal conductivity, a combination that is difficult to achieve. Lowered thermal conductivity has increased efficiency, but at the cost of power output. An alternative setup is to use metallic junctions rather than semiconductors as thermoelectric devices. Metals have orders of magnitude more electrons and electronic conductivities higher than semiconductors, but thermal conductivity is higher as well. To evaluate the viability of metallic junction thermoelectrics, a two dimensional heat transfer MATLAB simulation was constructed to calculate efficiency and power output. High Seebeck coefficient alloys, Chromel (90%Ni-10%Cr) and Constantan (55%Cu-45%Ni), produced efficiencies of around 20-30%. Parameters such as the number of layers of junctions, lateral junction density, and junction sizes for both series- and parallel-connected junctions were explored.

Duzik, Adam J.↗
NASA NTRS1989

Application of direct solvers to unstructured meshes for the Euler and Navier-Stokes equations using upwind schemes

The application of Newton iteration to inviscid and viscous airfoil calculations on unstructured meshes is examined. A cell-centered finite volume scheme is employed on an unstructured mesh consisting of triangles. Roe's flux difference splitting scheme is used to compute the inviscid fluxes. Higher order accuracy is achieved by an interpolation procedure that makes use of auxiliary gradients. The efficient solution of the sparse linear system of equations which arises upon linearization in time is addressed. Results are presented for inviscid and viscous test cases. The complications which arise due to the introduction of nonlinear limiters are addressed.

Venkatakrishnan, V.↗
NASA NTRS1989

Advanced development of the boundary element method for steady-state heat conduction

Considerable progress has been made in recent years toward advancing the state-of-the-art in solid mechanics boundary element technology. In the present work, much of this new technology is applied in the development of a general-purpose boundary element method (BEM) for steady-state heat conduction. In particular, the BEM implementation involves the use of higher-order conforming elements, self-adaptive integration and multi-region capability. Two- and three-dimensional, as well as axisymmetric analysis, are incorporated within a unified framework. In addition, techniques are introduced for the calculation of boundary flux, and for the inclusion of thermal resistance across interfaces. As a final extension, an efficient formulation is developed for the analysis of solid three-dimensional bodies with embedded holes. For this last class of problems, the new BEM formulation is particularly attractive, since use of the alternatives (i.e. finite element or finite difference methods) is not practical. A number of detailed examples illustrate the suitability and robustness of the present approach for steady-state heat conduction.

Dargush, G. F.↗
NASA NTRS2001

2D Quantum Transport Modeling in Nanoscale MOSFETs

We have developed physical approximations and computer code capable of realistically simulating 2-D nanoscale transistors, using the non-equilibrium Green's function (NEGF) method. This is the most accurate full quantum model yet applied to 2-D device simulation. Open boundary conditions, oxide tunneling and phase-breaking scattering are treated on an equal footing. Electron bandstructure is treated within the anisotropic effective mass approximation. We present the results of our simulations of MIT 25 and 90 nm "well-tempered" MOSFETs and compare them to those of classical and quantum corrected models. The important feature of quantum model is smaller slope of Id-Vg curve and consequently higher threshold voltage. These results are consistent with 1D Schroedinger-Poisson calculations. The effect of gate length on gate-oxide leakage and subthreshold current has been studied. The shorter gate length device has an order of magnitude smaller leakage current than the longer gate length device without a significant trade-off in on-current.

Svizhenko, Alexei↗
NASA NTRS2021

Evaluation of Higher-order Quadrature Schemes in Improving Computational Efficiency for Orientation-averaged Single-Scattering Properties of Nonspherical Ice Particles

We evaluate several high-order quadrature schemes for accuracy and efficacy in obtaining orientation-averaged single-scattering properties (SSPs). We use the highly efficient MIDAS to perform electromagnetic scattering calculations to evaluate the gain in efficiency from these schemes. MIDAS is shown to be superior to DDSCAT, a popular discrete dipole approximation (DDA) method. This study is motivated by the fact that quality physical precipitation retrievals rely on using accurate orientation-averaged SSPs derived from realistic hydrometeors as input to radiative transfer Models (RTMs). The DDA has been a popular choice for single-scattering calculations, due to its versatility with respect to target geometry. However, being iterative-solver-based (ISB), the most used DDA codes, e.g. DDSCAT and ADDA, must solve the scattering problem for each orientation of the target separately. As the size parameter and geometric anisotropy of the hydrometeor increase, the number of orientations needed to obtain accurate orientation-averages can increase drastically and so does the computation cost incurred by the ISB-DDA methods. MIDAS is a Direct-Solver-Based (DSB) code, its decomposition of the original large matrix with a high rank into multiple more manageable smaller matrices of lower ranks makes it much more computationally efficient while maintaining excellent accuracy. In addition, direct solvers consider all requested orientations at once, giving MIDAS further advantage over popular ISB-DDA methods. MIDAS, when combined with high-order quadrature for orientation averaging, can be greater than three orders of magnitude more efficient in obtaining RTM-ready SSPs of complex-shaped hydrometeors than existing ISB-DDA methods, with the native quadrature schemes they offer.

Ines Fenni↗
NASA NTRS1983

Calculation of potential flow past non-lifting bodies at angle of attack using axial and surface singularity methods

Two different singularity methods have been utilized to calculate the potential flow past a three dimensional non-lifting body. Two separate FORTRAN computer programs have been developed to implement these theoretical models, which will in the future allow inclusion of the fuselage effect in a pair of existing subcritical wing design computer programs. The first method uses higher order axial singularity distributions to model axisymmetric bodies of revolution in an either axial or inclined uniform potential flow. Use of inset of the singularity line away from the body for blunt noses, and cosine-type element distributions have been applied to obtain the optimal results. Excellent agreement to five significant figures with the exact solution pressure coefficient value has been found for a series of ellipsoids at different angles of attack. Solutions obtained for other axisymmetric bodies compare well with available experimental data. The second method utilizes distributions of singularities on the body surface, in the form of a discrete vortex lattice. This program is capable of modeling arbitrary three dimensional non-lifting bodies. Much effort has been devoted to finding the optimal method of calculating the tangential velocity on the body surface, extending techniques previously developed by other workers.

Shu, J. Y.↗
NASA NTRS1971

Test-particle method in kinetic theory of a plasma.

The introduction of a test particle into a system is considered. The system may be described by the Born-Bogoliubov-Green-Kirkwood-Yvon hierarchy. The field particles form a cloud which surrounds the test particle. The cloud is described by a conditional probability function which satisfies a certain equation. A generalization of the superposition principle reported by Rostoker (1964) to higher order correlation functions is discussed. Kinetic equations with the generalized Lenard-Balescu term are obtained, taking into account also diffusion by waves. The characteristics regarding the absorption or emission of waves by particles can be calculated.

Matsuda, K.↗
NASA NTRS1980

Anorthite - Thermal equation of state to high pressures

New shock-wave data are presented for anorthite from which a full high-temperature, high-pressure equation of state is derived. Whereas anorthite has relatively low values of thermal expansion and Grueneisen parameter at zero pressure, it is found that these attain relatively high values in the high density state corresponding to the high-pressure phase Hugoniot but decrease upon compression as expected. It is noted that higher order anharmonic contributions decrease more rapidly with pressure and that the thermal expansion therefore saturates to a high temperature value at pressures above about 100 GPa. Reduction of the Hugoniot data permits shock temperatures to be calculated; it also yields a principal adiabat for the high pressure branch of the Hugoniot. The initial bulk modulus of this adiabat is essentially identical to that of anorthite, whereas the initial density is about 3.40 Mg/cu m.

Jeanloz, R.↗