Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “linear equations”

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 145 records · Page 8
All sourcesNASA NTRSDOE OSTI
NASA NTRS1998

WARP3D-Release 10.8: Dynamic Nonlinear Analysis of Solids using a Preconditioned Conjugate Gradient Software Architecture

This report describes theoretical background material and commands necessary to use the WARP3D finite element code. WARP3D is under continuing development as a research code for the solution of very large-scale, 3-D solid models subjected to static and dynamic loads. Specific features in the code oriented toward the investigation of ductile fracture in metals include a robust finite strain formulation, a general J-integral computation facility (with inertia, face loading), an element extinction facility to model crack growth, nonlinear material models including viscoplastic effects, and the Gurson-Tver-gaard dilatant plasticity model for void growth. The nonlinear, dynamic equilibrium equations are solved using an incremental-iterative, implicit formulation with full Newton iterations to eliminate residual nodal forces. The history integration of the nonlinear equations of motion is accomplished with Newmarks Beta method. A central feature of WARP3D involves the use of a linear-preconditioned conjugate gradient (LPCG) solver implemented in an element-by-element format to replace a conventional direct linear equation solver. This software architecture dramatically reduces both the memory requirements and CPU time for very large, nonlinear solid models since formation of the assembled (dynamic) stiffness matrix is avoided. Analyses thus exhibit the numerical stability for large time (load) steps provided by the implicit formulation coupled with the low memory requirements characteristic of an explicit code. In addition to the much lower memory requirements of the LPCG solver, the CPU time required for solution of the linear equations during each Newton iteration is generally one-half or less of the CPU time required for a traditional direct solver. All other computational aspects of the code (element stiffnesses, element strains, stress updating, element internal forces) are implemented in the element-by- element, blocked architecture. This greatly improves vectorization of the code on uni-processor hardware and enables straightforward parallel-vector processing of element blocks on multi-processor hardware.

Koppenhoefer, Kyle C.↗
NASA NTRS2001

The Detection of Radiated Modes from Ducted Fan Engines

The bypass duct of an aircraft engine is a low-pass filter allowing some spinning modes to radiate outside the duct. The knowledge of the radiated modes can help in noise reduction, as well as the diagnosis of noise generation mechanisms inside the duct. We propose a nonintrusive technique using a circular microphone array outside the engine measuring the complex noise spectrum on an arc of a circle. The array is placed at various axial distances from the inlet or the exhaust of the engine. Using a model of noise radiation from the duct, an overdetermined system of linear equations is constructed for the complex amplitudes of the radial modes for a fixed circumferential mode. This system of linear equations is generally singular, indicating that the problem is illposed. Tikhonov regularization is employed to solve this system of equations for the unknown amplitudes of the radiated modes. An application of our mode detection technique using measured acoustic data from a circular microphone array is presented. We show that this technique can reliably detect radiated modes with the possible exception of modes very close to cut-off.

Farassat, F.↗
NASA NTRS1994

On Dirac equations for linear magnetoacoustic waves propagating in an isothermal atmosphere

A new analytical approach to study linear magnetoacoustic waves propagating in an isothermal, stratified, and uniformly magnetized atmosphere is presented. The approach is based on Dirac equations, and the theory of Sturm-Liouville operators is used to investigate spectral properties of the obtained Dirac Hamiltonians. Two cases are considered: (1) the background magnetic field is vertical, and the waves are separated into purely magnetic (transverse) and purely acoustic (longitudinal) modes; and (2) the field is tilted with respect to the vertical direction and the magnetic and acoustic modes become coupled giving magnetoacoustic waves. For the first case, the Dirac Hamiltonian possesses either a discrete spectrum, which corresponds to standing magnetic waves, or a continuous spectrum, which can be clearly identified with freely propagating acoustic waves. For the second case, the quantum mechanical perturbation calculus is used to study coupling and energy exchange between the magnetic and acoustic components of magnetoacoustic waves. It is shown that this coupling may efficiently prevent trapping of magnetoacoustic waves instellar atmospheres.

Alicki, R.↗
NASA NTRS2019

The Impact of Neglecting Ice Phase on Cloud Optical Depth Retrievals from AERONET Cloud Mode Observations

Clouds present many challenges to climate modelling. To develop and verify the parameterisations needed to allow climate models to represent cloud structure and processes, there is a need for high-quality observations of cloud optical depth from locations around the world. Retrievals of cloud optical depth are obtainable from radiances measured by Aerosol Robotic Network (AERONET) radiometers in “cloud mode” using a two-wavelength retrieval method. However, the method is unable to detect cloud phase, and hence assumes that all of the cloud in a profile is liquid. This assumption has the potential to introduce errors into long term statistics of retrieved optical depth for clouds that also contain ice. Using a set of idealized cloud profiles we find that, for optical depths above 20, the fractional error in retrieved optical depth is a linear function of the fraction of the optical depth that is due to the presence of ice cloud (“ice fraction”). Clouds that are entirely ice have positive errors with magnitudes of the order of 55% to 70 %. We derive a simple linear equation that can be used as a correction at AERONET sites where ice fraction can be independently estimated. Using this linear equation, we estimate the magnitude of the error for a set of cloud profiles from five sites of the Atmospheric Radiation Measurement programme. The dataset contains separate retrievals of ice and liquid retrievals; hence ice fraction can be estimated. The magnitude of the error at each location was related to the relative frequencies of occurrence in thick frontal cloud at the mid-latitude sites and of deep convection at the tropical sites – that is, of deep cloud containing both ice and liquid particles. The long-term mean optical depth error at the five locations spans the range 2– 4, which we show to be small enough to allow calculation of top-of-atmosphere flux to within 10% and surface flux to about 15 %.

Shonk, Jonathan K. P.↗
NASA NTRS1993

Methodology for sensitivity analysis, approximate analysis, and design optimization in CFD for multidisciplinary applications

In this study involving advanced fluid flow codes, an incremental iterative formulation (also known as the delta or correction form) together with the well-known spatially-split approximate factorization algorithm, is presented for solving the very large sparse systems of linear equations which are associated with aerodynamic sensitivity analysis. For smaller 2D problems, a direct method can be applied to solve these linear equations in either the standard or the incremental form, in which case the two are equivalent. Iterative methods are needed for larger 2D and future 3D applications, however, because direct methods require much more computer memory than is currently available. Iterative methods for solving these equations in the standard form are generally unsatisfactory due to an ill-conditioning of the coefficient matrix; this problem can be overcome when these equations are cast in the incremental form. These and other benefits are discussed. The methodology is successfully implemented and tested in 2D using an upwind, cell-centered, finite volume formulation applied to the thin-layer Navier-Stokes equations. Results are presented for two sample airfoil problems: (1) subsonic low Reynolds number laminar flow; and (2) transonic high Reynolds number turbulent flow.

Taylor, Arthur C., III↗
NASA NTRS1978

Non-linear effects on sound propagation in a variable area duct carrying high subsonic Mach number flows

The general one-dimensional non-linear equation for the acoustic velocity potential in a variable area duct carrying high subsonic Mach number flows is presented and solved numerically using an implicit finite difference scheme. For the linearized equation with no flow, the present scheme is compared with the exact solution and the fourth-order Runge-Kutta method with excellent agreement for dimensionless time periods. Non-linear solutions are more sensitive to Mach numbers and exciting amplitudes and less sensitive to exciting frequencies. In general, non-linear effects can be safely neglected for low Mach number flows but must be accounted for when high Mach number flows are encountered.

Lumsdaine, E.↗
NASA NTRS1991

Direct time integration of Maxwell's equations in linear dispersive media with absorption for scattering and propagation of femtosecond electromagnetic pulses

The initial results for femtosecond pulse propagation and scattering interactions for a Lorentz medium obtained by a direct time integration of Maxwell's equations are reported. The computational approach provides reflection coefficients accurate to better than 6 parts in 10,000 over the frequency range of dc to 3 x 10 to the 16th Hz for a single 0.2-fs Gaussian pulse incident upon a Lorentz-medium half-space. New results for Sommerfeld and Brillouin precursors are shown and compared with previous analyses. The present approach is robust and permits 2D and 3D electromagnetic pulse propagation directly from the full-vector Maxwell's equations.

Joseph, Rose M.↗
NASA NTRS1959

Investigation of the Stability of Very Flat Spins and Analysis of Effects of Applying various Moments Utilizing the Three Moment Equations of Motion

Based on linearized equations of motion utilizing only the three moment equations and assuming only flat-spin conditions, it appears that contemporary designs (with the moment of inertia about the wing axis I(sub Y) considerably greater than the moment of inertia about the fuselage axis I(sub X) having positive values of C(sub l, sub p) (rolling-moment coefficient due to rolling) or positive values of C(sub l, sub beta) (rolling-moment coefficient due to sideslip) will probably not have a stable spin in the flat-spin region near an angle of attack of 90 deg. If the damping in pitch in flat-spin attitudes is zero, stable flat-spin conditions may not be possible on an airplane having the mass primarily distributed along the wings. The effect of moving ailerons with the spin or the effect of applying a positive pitching moment producing recovery for contemporary fighter designs will be greatest for large negative values of C(sub n, sub beta) (yawing-moment coefficient due to sideslip). In addition, for a certain critical value of positive C(sub n, sub beta), the rolling moment applied by moving ailerons with the spin or the application of a positive pitching moment will have no effect on reducing the spin rate.

Klinar, Walter J.↗
NASA NTRS1982

Discrete Fourier transforms of nonuniformly spaced data

Time series or spatial series of measurements taken with nonuniform spacings have failed to yield fully to analysis using the Discrete Fourier Transform (DFT). This is due to the fact that the formal DFT is the convolution of the transform of the signal with the transform of the nonuniform spacings. Two original methods are presented for deconvolving such transforms for signals containing significant noise. The first method solves a set of linear equations relating the observed data to values defined at uniform grid points, and then obtains the desired transform as the DFT of the uniform interpolates. The second method solves a set of linear equations relating the real and imaginary components of the formal DFT directly to those of the desired transform. The results of numerical experiments with noisy data are presented in order to demonstrate the capabilities and limitations of the methods.

Swan, P. R.↗
NASA NTRS1996

Methodology for Sensitivity Analysis, Approximate Analysis, and Design Optimization in CFD for Multidisciplinary Applications

An incremental iterative formulation together with the well-known spatially split approximate-factorization algorithm, is presented for solving the large, sparse systems of linear equations that are associated with aerodynamic sensitivity analysis. This formulation is also known as the 'delta' or 'correction' form. For the smaller two dimensional problems, a direct method can be applied to solve these linear equations in either the standard or the incremental form, in which case the two are equivalent. However, iterative methods are needed for larger two-dimensional and three dimensional applications because direct methods require more computer memory than is currently available. Iterative methods for solving these equations in the standard form are generally unsatisfactory due to an ill-conditioned coefficient matrix; this problem is overcome when these equations are cast in the incremental form. The methodology is successfully implemented and tested using an upwind cell-centered finite-volume formulation applied in two dimensions to the thin-layer Navier-Stokes equations for external flow over an airfoil. In three dimensions this methodology is demonstrated with a marching-solution algorithm for the Euler equations to calculate supersonic flow over the High-Speed Civil Transport configuration (HSCT 24E). The sensitivity derivatives obtained with the incremental iterative method from a marching Euler code are used in a design-improvement study of the HSCT configuration that involves thickness. camber, and planform design variables.

Taylor, Arthur C., III↗
NASA NTRS2003

Development of Numerical Codes for Modeling Electromagnetic Behavior at High Frequencies Near Large Objects

A study into the problem of determining electromagnetic solutions at high frequencies for problems involving complex geometries, large sizes and multiple sources (e.g. antennas) has been initiated. Typical applications include the behavior of antennas (and radiators) installed on complex conducting structures (e.g. ships, aircrafts, etc..) with strong interactions between antennas, the radiation patterns, and electromagnetic signals is of great interest for electromagnetic compatibility control. This includes the overall performance evaluation and control of all on-board radiating systems, electromagnetic interference, and personnel radiation hazards. Electromagnetic computational capability exists at NASA LaRC, and many of the codes developed are based on the Moment Method (MM). However, the MM is computationally intensive, and this places a limit on the size of objects and structures that can be modeled. Here, two approaches are proposed: (i) a current-based hybrid scheme that combines the MM with Physical optics, and (ii) an Alternating Direction Implicit-Finite Difference Time Domain (ADI-FDTD) method. The essence of a hybrid technique is to split the overall scattering surface(s) into two regions: (a) a MM zone (MMZ) which can be used over any part of the given geometry, but is most essential over irregular and "non-smooth" geometries, and (b) a PO sub-region (POSR). Currents induced on the scattering and reflecting surfaces can then be computed in two ways depending on whether the region belonged to the MMZ or was part of the POSR. For the MMZ, the current calculations proceed in terms of basis functions with undetermined coefficients (as in the usual MM method), and the answer obtained by solving a system of linear equations. Over the POSR, conduction is obtained as a superposition of two contributions: (i) currents due to the incident magnetic field, and (ii) currents produced by the mutual induction from conduction within the MMZ. This effectively leads to a reduction in the size of linear equations from N to N - Npo with N being the total number of segments for the entire surface and Npo the number of segments over the POSR. The scheme would be appropriate for relatively large, flat surfaces, and at high frequencies. The ADI-FDTD scheme provides for both transient and steady state analyses. The restrictive Courant-Friedrich-Levy (CFL) condition on the time-step is removed, and so large time steps can be chosen even though the spatial grids are small. This report includes the problem definition, a detailed discussion of both the numerical techniques, and numerical implementations for simple surface geometries. Numerical solutions have been derived for a few simple situations.

Joshi, R. P.↗
NASA NTRS1971

The responses of rising or falling spherical wind sensors to atmospheric wind perturbations.

Analysis of the responses of rising or falling spherical wind sensors to atmospheric wind perturbations on the wind profile in the vertical, using Fourier transform techniques. The linearized equations of motion of a sensor that is subject to drag and gravitational body forces are developed by perturbing a sensor about an equilibrium uniform motion with wind fluctuations which have vertical variations. The wind environment and sensor velocities are decomposed with stochastic Fourier-Stieltjes integrals, and the linearized equations of motion are used to derive the response functions and phase angles of the sensor motions. The results of the analysis are used to analyze the response properties of the Jimsphere balloon wind sensor. It is shown that, in general, the transfer functions associated with the horizontal sensor motions are smaller than the transfer functions associated with the vertical sensor motions in the omega times T product range from zero to infinity, omega being the wind perturbation frequency and T a time constant of the system. Thus, the sensor is more responsive to vertical than to horizontal air motions.

Fichtl, G. H.↗
NASA NTRS1951

Integrals and integral equations in linearized wing theory

The formulas of subsonic and supersonic wing theory for source, doublet, and vortex distributions are reviewed and a systematic presentation is provided which relates these distributions to the pressure and to the vertical induced velocity in the plane of the wing. It is shown that care must be used in treating the singularities involved in the analysis and that the order of integration is not always reversible. Concepts suggested by the irreversibility of order of integration are shown to be useful in the inversion of singular integral equations when operational techniques are used. A number of examples are given to illustrate the methods presented, attention being directed to supersonic flight speed.

Lomax, Harvard↗
NASA NTRS1976

Sound propagation in choked ducts

The linearized equations describing the propagation of sound in variable area ducts containing flow are shown to be singular when the duct mean flow is sonic. The singularity is removed when previously ignored nonlinear terms are retained. The results of a numerical study, for the case of plane waves propagating in a one-dimensional converging-diverging duct, show that the sound field is adequately described by the linearized equations only when the axial mean flow Mach number at the duct throat M sub th 0.6. For M sub th 0.6, the numerical results showed that acoustic energy flux was not conserved. An attempt was made to extend the study to include the nonlinear behavior of the sound field. Meaningful results were not obtained due, primarily, to numerical difficulties.

Hersh, A. S.↗
NASA NTRS1946

Numerical procedures for the calculation of the stresses in monocoques III : calculation of the bending moments in fuselage frames

This report deals with the calculation of the bending moments in and the distortions of fuselage rings upon which known concentrated and distributed loads are acting. In the procedure suggested, the ring is divided into a number of beams each having a constant radius of curvature. The forces and moments caused in the end sections of the beams by individual unit displacements of the end sections are listed in a table designated as the operations table in conformity with Southwell's nomenclature. The operations table and the external loads are equivalent to a set of linear equations. For their solution the following three procedures are presented: 1) Southwell's method of systematic relaxations. This is a step-by-step approximation procedure guided by the physical interpretation of the changes in the values of the unknown. 2) The growing unit procedure in which the individual beams are combined successively into beams of increasing length until finally the entire ring becomes a single beam. In each step of the procedure a set of not more than three simultaneous linear equations is solved. 3) Solution of the entire set of simultaneous equations by the methods of the matrix calculus. In order to demonstrate the manner in which the calculations may be carried out, the following numerical examples are worked out: 1) Curved beam with both its end sections rigidly fixed. The load is a concentrated force. 2) Egg-shape ring with symmetric concentrated loads. 3) Circular ring with antisymmetric concentrated loads and shear flow (torsion of the fuselage). 4) Same with V-braces incorporated in the ring. 5) Egg-shape ring with antisymmetric concentrated loads and shear flow (torsion of the fuselage). 6) Same with V-braces incorporated in the ring. The results of these calculations are checked, whenever possible, by calculations carried out according to known methods of analysis. The agreement is found to be good. The amount of work necessary for the solution of ring problems by the methods described in the present report is practically independent of the degree of redundancy of the structure. For this reason the methods are recommended for use particularly in problems of rings having one or more internal bracing elements.

LOADS - FUSELAGES↗
NASA NTRS2007

USM3D Unstructured Grid Solutions for CAWAPI at NASA LaRC

In support the Cranked Arrow Wing Aerodynamic Project International (CAWAPI) to improve the Technology Readiness Level of flow solvers by comparing results with measured F-16XL-1 flight data, NASA Langley employed the TetrUSS unstructured grid solver, USM3D, to obtain solutions for all seven flight conditions of interest. A newly available solver version that incorporates a number of turbulence models, including the two-equation linear and non-linear k-epsilon, was used in this study. As a first test, a choice was made to utilize only a single grid resolution with the solver for the simulation of the different flight conditions. Comparisons are presented with three turbulence models in USM3D, flight data for surface pressure, boundary-layer profiles, and skin-friction results, as well as limited predictions from other solvers. A result of these comparisons is that the USM3D solver can be used in an engineering environment to predict flow physics on a complex configuration at flight Reynolds numbers with a two-equation linear k-epsilon turbulence model.

Lamar, John E.↗