Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Algebraic models”

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 829 records · Page 46

Reduced basis technique for calculating sensitivity coefficients of nonlinear structural response

An efficient reduced basis technique is presented for calculating the sensitivity of nonlinear structural response to variations in the design variables. The structure is discretized by using two-field mixed finite element models. The vector of structural response and its sensitivity coefficients (derivatives with respect to design variables) are each expressed as a linear combination of a small number of basis (or global approximation) vectors. The Bubnov-Galerkin technique is then used to approximate each of the finite element equations governing the response and the sensitivity coefficients by a small number of algebraic equations in the amplitudes of these vectors. The path derivatives (derivatives of the response vector with respect to path parameters, e.g., load parameters) are used as basis vectors for approximating the response. A combination of the path derivatives and their derivatives with respect to the design variables is used for approximating the sensitivity coefficients. The potential of the proposed technique is discussed and its effectiveness is demonstrated by means of numerical examples of laminated composite plates subjected to mechanical and thermal loads.

Noor, Ahmed K.↗

Linear decentralized systems with special structure

Certain fundamental structures associated with linear systems having internal symmetries are outlined. It is shown that the theory of finite-dimensional algebras and their representations are closely related to such systems. It is also demonstrated that certain problems in the decentralized control of symmetric systems are equivalent to long-standing problems of linear systems theory. Even though the structure imposed arose in considering the problems of twin-lift helicopters, any large system composed of several identical intercoupled control systems can be modeled by a linear system that satisfies the constraints imposed. Internal symmetry can be exploited to yield new system-theoretic invariants and a better understanding of the way in which the underlying structure affects overall system performance.

Martin, C. F.↗

Coupled Riccati equations for complex plane constraint

A new Linear Quadratic Gaussian design method is presented which provides prescribed imaginary axis pole placement for optimal control and estimation systems. This procedure contributes another degree of design freedom to flexible spacecraft control. Current design methods which interject modal damping into the system tend to have little affect on modal frequencies, i.e., they predictably shift open plant poles horizontally in the complex plane to form the closed loop controller or estimator pole constellation, but make little provision for vertical (imaginary axis) pole shifts. Imaginary axis shifts which reduce the closed loop model frequencies (the bandwidths) are desirable since they reduce the sensitivity of the system to noise disturbances. The new method drives the closed loop modal frequencies to predictable (specified) levels, frequencies as low as zero rad/sec (real axis pole placement) can be achieved. The design procedure works through rotational and translational destabilizations of the plant, and a coupling of two independently solved algebraic Riccati equations through a structured state weighting matrix. Two new concepts, gain transference and Q equivalency, are introduced and their use shown.

Strong, Kristin M.↗

Modeling Sound Propagation Through Non-Axisymmetric Jets

A method for computing the far-field adjoint Green's function of the generalized acoustic analogy equations under a locally parallel mean flow approximation is presented. The method is based on expanding the mean-flow-dependent coefficients in the governing equation and the scalar Green's function in truncated Fourier series in the azimuthal direction and a finite difference approximation in the radial direction in circular cylindrical coordinates. The combined spectral/finite difference method yields a highly banded system of algebraic equations that can be efficiently solved using a standard sparse system solver. The method is applied to test cases, with mean flow specified by analytical functions, corresponding to two noise reduction concepts of current interest: the offset jet and the fluid shield. Sample results for the Green's function are given for these two test cases and recommendations made as to the use of the method as part of a RANS-based jet noise prediction code.

Leib, Stewart J.↗

Inverse methods for assessing ship-of-opportunity networks and estimating circulation and winds from tropical expendable bathythermograph data

Inverse methods for estimating the surface ciculation of the equatorial Pacific by combining a linear reduced-gravity shallow-water model with the Tropical Ocean-Global Atmosphere ship-of-opportunity expendable bathythermograph (TOGA SOP XBT) observing program are examined. It is demonstrated that a simple linear model of the upper circulation of the equatorial Pacific can be successfully used as a weak constraint when smoothing the TOGA SOP XBT data. A circulation is sought as the weighted least squares fit to the dynamics and the data. The solution method is an expansion in representer functions, and the generalized inverse problem is thereby reduced from a functional problem to an algebraic problem for the coefficients of the representer. A specific inverse calculation using synthetic forcing and data is presented.

Bennett, Andrew F.↗

Mixed isoparametric finite element models of laminated composite shells

Mixed shear-flexible isoparametric elements are presented for the stress and free vibration analysis of laminated composite shallow shells. Both triangular and quadrilateral elements are considered. The 'generalized' element stiffness, consistent mass, and consistent load coefficients are obtained by using a modified form of the Hellinger-Reissner mixed variational principle. Group-theoretic techniques are used in conjunction with computerized symbolic integration to obtain analytic expressions for the stiffness, mass and load coefficients. A procedure is outlined for efficiently handling the resulting system of algebraic equations. The accuracy of the mixed isoparametric elements developed is demonstrated by means of numerical examples, and their advantages over commonly used displacement elements are discussed.

Noor, A. K.↗

Nonlinear evolution of magnetopause tearing modes

Since the magnetosheath plasma is highly turbulent, reconnection at the dayside magnetopause is likely to be temporally unsteady. The tearing mode can be viewed as a model for the unsteady development of a reconnecting magnetic topology. Magnetopause tearing occurs in the guide-field limit and has a wave packet spatial structure in the east-west direction. This paper solves for the nonlinear evolution of a single wavelength guide-field tearing mode including the effects of finite transit time on the Landau resonant electrons. Short wavelength modes evolve algebraically in time with perturbation amplitudes proportional to t-squared. Long wavelength modes are fully nonlinear, and the amplitude grows linearly in time.

Coroniti, F. V.↗

Calculation of plasma properties in ion sources

Simple algebraic equations are developed that can be used to calculate the average primary electron density, the average primary to Maxwellian density ratio, and the average Maxwellian electron temperature in the discharge chamber of a cusped magnetic field ion source. These properties may be calculated as functions of the propellant gas flow rate and the propellant utilization efficiency. Each of the plasma properties calculated using the model agrees well with the corresponding experimental data for both argon and xenon gases over a wide range of propellant flow rates and utilization efficiencies.

Brophy, J. R.↗

Time-optical spinup maneuvers of flexible spacecraft

Attitude controllers for spacecraft have been based on the assumption that the bodies being controlled are rigid. Future spacecraft, however, may be quite flexible. Many applications require spinning up/down these vehicles. In this work the minimum time control of these maneuvers is considered. The time-optimal control is shown to possess an important symmetry property. Taking advantage of this property, the necessary and sufficient conditions for optimality are transformed into a system of nonlinear algebraic equations in the control switching times during one half of the maneuver, the maneuver time, and the costates at the mid-maneuver time. These equations can be solved using a homotopy approach. Control spillover measures are introduced and upper bounds on these measures are obtained. For a special case these upper bounds can be expressed in closed form for an infinite dimensional evaluation model. Rotational stiffening effects are ignored in the optimal control analysis. Based on a heuristic argument a simple condition is given which justifies the omission of these nonlinear effects. This condition is validated by numerical simulation.

Singh, G.↗

A generalization of algebraic surface drawing

An implicit surface mathematical description of three-dimensional space is defined in terms of all points which satisfy some equation F(x, y, z) equals 0. This form is ideal for space-shaded picture drawing, where the coordinates are substituted for x and y and the equation is solved for z. A new algorithm is presented which is applicable to functional forms other than those of first- and second-order polynomial functions, such as the summation of several Gaussian density distributions. The algorithm was created in order to model electron density maps of molecular structures, but is shown to be capable of generating shapes of esthetic interest.

Blinn, J. F.↗

Parametric study of stress-intensity factors in bonded composite stringer panels

Stress-intensity factors are determined for an infinite cracked orthotropic sheet adhesively bonded to an orthotropic stringer. Since the stringer is modeled as a semi-infinite sheet, the solution is most appropriate for a crack tip located near a stringer edge. Both adherends are treated as homogeneous, orthotropic media which are representative of many fiber-reinforced composite materials. The complex variable theory of elasticity was used to obtain a set of integral equations describing the problem. The integral equations are replaced by an equivalent sheet of algebraic equations, which are solved to obtain the shear stress distribution in the adhesive layer. From these adhesive stresses, the stress-intensity factors are found. A parametric study is conducted to determine the sensitivity of the system to material properties and specimen configuration. Unless the crack tip is very close to or under the stringer, the stress-intensity factor is approximately that of the unstiffened sheet. However, as the crack propagates beneath the stringer, the stress-intensity factor decreases significantly. Increasing the stringer stiffness or the adhesive stiffness also decreases the stress-intensity factor.

Bigelow, C. A.↗

Approach for Propagating Radiometric Data Uncertainties Through NASA Ocean Color Algorithms

Spectroradiometric satellite observations of the ocean are commonly referred to as “ocean color” remote sensing. NASA has continuously collected, processed, and distributed ocean color datasets since the launch of the Sea-viewing Wide-field-of-view Sensor (SeaWiFS) in 1997. While numerous ocean color algorithms have been developed in the past two decades that derive geophysical data products from sensor-observed radiometry, few papers have clearly demonstrated how to estimate measurement uncertainty in derived data products. As the uptake of ocean color data products continues to grow with the launch of new and advanced sensors, it is critical that pixel-by-pixel data product uncertainties are estimated during routine data processing. Knowledge of uncertainties can be used when studying long-term climate records, or to assist in the development and performance appraisal of bio-optical algorithms. In this method paper we provide a comprehensive overview of how to formulate first-order first-moment (FOFM) calculus for propagating radiometric uncertainties through a selection of bio-optical models. We demonstrate FOFM uncertainty formulations for the following NASA ocean color data products: chlorophyll-a pigment concentration (Chl), the diffuse attenuation coefficient at 490 nm (K(sub d,490)), particulate organic carbon (POC), normalized fluorescent line height (nflh), and inherent optical properties (IOPs). Using a quality-controlled in situ hyperspectral remote sensing reflectance (R(sub rs,i)) dataset, we show how computationally inexpensive, yet algebraically complex, FOFM calculations may be evaluated for correctness using the more computationally expensive Monte Carlo approach. We compare bio-optical product uncertainties derived using our test R(sub rs) dataset assuming spectrally-flat, uncorrelated relative uncertainties of 1, 5, and 10%. We also consider spectrally dependent, uncorrelated relative uncertainties in R(sub rs). The importance of considering spectral covariances in R(sub rs), where practicable, in the FOFM methodology is highlighted with an example SeaWiFS image. We also present a brief case study of two POC algorithms to illustrate how FOFM formulations may be used to construct measurement uncertainty budgets for ecologically-relevant data products. Such knowledge, even if rudimentary, may provide useful information to end-users when selecting data products or when developing their own algorithms.

Bio-optics↗

Maximum likelihood identification of aircraft parameters with unsteady aerodynamic modelling

A simplified aerodynamic force model based on the physical principle of Prandtl's lifting line theory and trailing vortex concept has been developed to account for unsteady aerodynamic effects in aircraft dynamics. Longitudinal equations of motion have been modified to include these effects. The presence of convolution integrals in the modified equations of motion led to a frequency domain analysis utilizing Fourier transforms. This reduces the integro-differential equations to relatively simple algebraic equations, thereby reducing computation time significantly. A parameter extraction program based on the maximum likelihood estimation technique is developed in the frequency domain. The extraction algorithm contains a new scheme for obtaining sensitivity functions by using numerical differentiation. The paper concludes with examples using computer generated and real flight data

Keskar, D. A.↗

An Eigenvalue Analysis of finite-difference approximations for hyperbolic IBVPs

The eigenvalue spectrum associated with a linear finite-difference approximation plays a crucial role in the stability analysis and in the actual computational performance of the discrete approximation. The eigenvalue spectrum associated with the Lax-Wendroff scheme applied to a model hyperbolic equation was investigated. For an initial-boundary-value problem (IBVP) on a finite domain, the eigenvalue or normal mode analysis is analytically intractable. A study of auxiliary problems (Dirichlet and quarter-plane) leads to asymptotic estimates of the eigenvalue spectrum and to an identification of individual modes as either benign or unstable. The asymptotic analysis establishes an intuitive as well as quantitative connection between the algebraic tests in the theory of Gustafsson, Kreiss, and Sundstrom and Lax-Richtmyer L(sub 2) stability on a finite domain.

Warming, Robert F.↗

An eigenvalue analysis of finite-difference approximations for hyperbolic IBVPs

The eigenvalue spectrum associated with a linear finite-difference approximation plays a crucial role in the stability analysis and in the actual computational performance of the discrete approximation. The eigenvalue spectrum associated with the Lax-Wendroff scheme applied to a model hyperbolic equation was investigated. For an initial-boundary-value problem (IBVP) on a finite domain, the eigenvalue or normal mode analysis is analytically intractable. A study of auxiliary problems (Dirichlet and quarter-plane) leads to asymptotic estimates of the eigenvalue spectrum and to an identification of individual modes as either benign or unstable. The asymptotic analysis establishes an intuitive as well as quantitative connection between the algebraic tests in the theory of Gustafsson, Kreiss, and Sundstrom and Lax-Richtmyer L (sub 2) stability on a finite domain.

Warming, Robert F.↗

A tire contact solution technique

An efficient method for calculating the contact boundary and interfacial pressure distribution was developed. This solution technique utilizes the discrete Fourier transform to establish an influence coefficient matrix for the portion of the pressurized tire surface that may be in the contact region. This matrix is used in a linear algebra algorithm to determine the contact boundary and the array of forces within the boundary that are necessary to hold the tire in equilibrium against a specified contact surface. The algorithm also determines the normal and tangential displacements of those points on the tire surface that are included in the influence coefficient matrix. Displacements within and outside the contact region are calculated. The solution technique is implemented with a finite-element tire model that is based on orthotropic, nonlinear shell of revolution elements which can respond to nonaxisymmetric loads. A sample contact solution is presented.

Tielking, J. T.↗

A cracked orthotropic sheet stiffened by a semi-infinite orthotropic sheet

The stress intensity factor is determined for a cracked orthotropic sheet adhesively bonded to an orthotropic stringer. Since the stringer is modeled as a semiinfinite sheet, the solution is most appropriate for a crack tip located near a stringer edge. Both adherends are treated as homogeneous, orthotropic media. It is assumed they are in plane stress and the adhesive is in pure shear. From Green's functions and the complex variable theory of orthotropic elasticity developed by Lekhnitskiy, a set of integral equations is obtained. The integral equations are replaced by an equivalent set of algebraic equations which is solved to obtain the shear stress distribution in the adhesive layer. With these stresses, equations for the stress intensity factors at both crack tips are found. A parametric study is conducted to determine the sensitivity of the system to material properties and specimen configuration. Unless the crack tip is very close to or under the stringer the stress intensity factor is approximately that of an unstiffened sheet. However, as the crack propagates beneath the stringer, the stress intensity factor decreases significantly. Increasing the stiffness of the stringer or the adhesive also results in a decrease in the stress intensity factor.

Bigelow, C. A.↗

Flow solution on a dual-block grid around an airplane

The compressible flow around a complex fighter-aircraft configuration (fuselage, cranked delta wing, canard, and inlet) is simulated numerically using a novel grid scheme and a finite-volume Euler solver. The patched dual-block grid is generated by an algebraic procedure based on transfinite interpolation, and the explicit Runge-Kutta time-stepping Euler solver is implemented with a high degree of vectorization on a Cyber 205 processor. Results are presented in extensive graphs and diagrams and characterized in detail. The concentration of grid points near the wing apex in the present scheme is shown to facilitate capture of the vortex generated by the leading edge at high angles of attack and modeling of its interaction with the canard wake.

Eriksson, Lars-Erik↗