Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Matrix elements”

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 1,585 records · Page 88

Advanced composites for large Navy spacecraft

An overview is given of work conducted on contract for the Naval Sea Systems Command. The objective of this contract was to provide direction for the development of high modulus graphite reinforced metal matrix composites. These advanced materials can have a significant effect on the performance of a spacecraft before, during and after an evasive maneuver. The work conducted on this program was organized into seven technical tasks. Task 1 was development of a generic Navy spacecraft model. Finite element models of candidate structural designs were developed. In Task 2, the finite-element model(s) of the structure were used to conduct analytical assessments involving conventional materials, resin matrix composites and metal matrix composites (MMC). In Task 3 and 4, MMC material design, fabrication and evaluation was conducted. This consisted of generating material designs and developing a data base for a broad range of graphite reinforced MMC materials. All material was procured according to specifications which set material quality and material property standards. In Task 5, a set of evasive maneuvering requirements were derived and used in Task 6 to conduct analytical simulations. These analytical simulations used current SOA material properties and projected material properties to provide an indication of key payoffs for material development. In Task 7, a set of material development recommendations was generated.

Davis, William E.↗

General purpose program to generate compatibility matrix for the integrated force method

An efficient procedure for obtaining the compatibility conditions of finite-element models involves the generation of both field and compatibility conditions from deformation-displacement relations, using (1) the compatibility bandwidth, and (2) the node-determinacy concept. A computer program thus structured will generate sparse and banded compatibility conditions for a structure that is idealized by the finite elements.

Nagabhusanam, J.↗

Elastic plate spallation

The dynamic finite element method was used to investigate elastic stress waves in a plate. Strain displacement and stress strain relations are discussed along with the stiffness and mass matrix. The results of studying point load, and distributed load over small, intermediate, and large radii are reported. The derivation of finite element matrices, and the derivation of lumped and consistent matrices for one dimensional problems with Laplace transfer solutions are included. The computer program JMMSPALL is also included.

Oline, L.↗

Ocean Energy Systems Wave Energy Modeling Task 10.4: Numerical Modeling of a Fixed Oscillating Water Column

This paper reports on an ongoing international effort to establish guidelines for numerical modeling of wave energy converters, initiated by the International Energy Agency Technology Collaboration Program for Ocean Energy Systems. Initial results for point absorbers were presented in previous work, and here we present results for a breakwater-mounted Oscillating Water Column (OWC) device. The experimental model is at scale 1:4 relative to a full-scale installation in a water depth of 12.8 m. The power-extracting air turbine is modeled by an orifice plate of 1–2% of the internal chamber surface area. Measurements of chamber surface elevation, air flow through the orifice, and pressure difference across the orifice are compared with numerical calculations using both weakly-nonlinear potential flow theory and computational fluid dynamics. Both compressible- and incompressible-flow models are considered, and the effects of air compressibility are found to have a significant influence on the motion of the internal chamber surface. Recommendations are made for reducing uncertainties in future experimental campaigns, which are critical to enable firm conclusions to be drawn about the relative accuracy of the numerical models. It is well-known that boundary element method solutions of the linear potential flow problem (e.g., WAMIT) are singular at infinite frequency when panels are placed directly on the free surface. This is problematic for time-domain solutions where the value of the added mass matrix at infinite frequency is critical, especially for OWC chambers, which are modeled by zero-mass elements on the free surface. A straightforward rational procedure is described to replace ad-hoc solutions to this problem that have been proposed in the literature.

16 TIDAL AND WAVE POWER↗

A Framework for Integrated Component and System Analyses of Instabilities

Instabilities associated with fluid handling and operation in liquid rocket propulsion systems and test facilities usually manifest themselves as structural vibrations or some form of structural damage. While the source of the instability is directly related to the performance of a component such as a turbopump, valve or a flow control element, the associated pressure fluctuations as they propagate through the system have the potential to amplify and resonate with natural modes of the structural elements and components of the system. In this paper, the authors have developed an innovative multi-level approach that involves analysis at the component and systems level. The primary source of the unsteadiness is modeled with a high-fidelity hybrid RANS/LES based CFD methodology that has been previously used to study instabilities in feed systems. This high fidelity approach is used to quantify the instability and understand the physics associated with the instability. System response to the driving instability is determined through a transfer matrix approach wherein the incoming and outgoing pressure and velocity fluctuations are related through a transfer (or transmission) matrix. The coefficients of the transfer matrix for each component (i.e. valve, pipe, orifice etc.) are individually derived from the flow physics associated with the component. A demonstration case representing a test loop/test facility comprised of a network of elements is constructed with the transfer matrix approach and the amplification of modes analyzed as the instability propagates through the test loop.

Ahuja, Vineet↗

hr-Adaptivity for nonconforming high-order meshes with the target matrix optimization paradigm

In this study, we present an \(hr\) -adaptivity framework for optimization of high-order meshes. This work extends the r -adaptivity method by Dobrev et al. (Comput Fluids, 2020), where we utilized the Target-Matrix Optimization Paradigm (TMOP) to minimize a functional that depends on each element’s current and target geometric parameters: element aspect-ratio , size , skew , and rotation . Since fixed mesh topology limits the ability to achieve the target size and aspect-ratio at each position, in this paper, we augment the r -adaptivity framework with nonconforming adaptive mesh refinement to further reduce the error with respect to the target geometric parameters. The proposed formulation, referred to as \(hr\) -adaptivity, introduces TMOP-based quality estimators to satisfy the aspect-ratio target via anisotropic refinements and size target via isotropic refinements in each element of the mesh. The methodology presented is purely algebraic, extends to both simplices and hexahedra/quadrilaterals of any order, and supports nonconforming isotropic and anisotropic refinements in 2D and 3D. Using a problem with a known exact solution, we demonstrate the effectiveness of \(hr\) -adaptivity over both r - and \(h\) -adaptivity in obtaining similar accuracy in the solution with significantly fewer mesh nodes. We also present several examples that show that \(hr\) -adaptivity can help satisfy geometric targets even when \(r\) -adaptivity fails to do so, due to the topology of the initial mesh.

42 ENGINEERING↗

Degree and plane of polarization of multiple scattered light. 1: Homogeneous cloud layers

The degree of polarization and the direction of the plane of polarization are calculated by a Monte Carlo method for homogeneous layers. Two solar zenith angles and a range of optical thicknesses up to 10 are considered. The results are compared with calculations for single scattered photons. For a given pair of incident and scattered directions, there are only two possible values for the direction of the plane of polarization differing by 90 deg for single scattering from spherical aerosols. The choice between these two values depends only on the sign of the element M(-) in the first row and second column of the scattering matrix in the I, Q, U, V representation. In most cases there is little change in the direction of the plane of polarization when multiple scattering is taken into account, so that this quantity can usually be predicted from a very simple trigonometric relationship to good accuracy. Measurements of the direction of the plane of polarization at appropriately chosen angles provides information about the size distribution of the scattering centers.

Kattawar, G. W.↗

Major and trace element chemistry of Boulder 1 at Station 2, Apollo 17

Twenty-seven samples from Boulder 1 at Station 2 are analyzed for major and trace elements by atomic absorption spectrophotometry and neutron activation analysis. Two types of matrix and several types of clast materials are characterized on the basis of their chemistry. It is shown that one matrix type is a common material at the Apollo 17 site, while the other is probably exotic to that site. The most unusual clast materials found are coarse norite (an old rock no longer found in millimeter fragments at the site) and pigeonite basalt (possibly a highland volcanic rock). It is concluded that the boulder-forming process combined materials from at least two different localities or vertical strata.

Blanchard, D. P.↗

On a numerical solution of the supersonic panel flutter eigenproblem

An automated digital computer procedure is presented in this paper which enables efficient solution of the eigenvalue problem associated with the supersonic panel flutter phenomena. The step-by-step incremental solution procedure is based on an inverse iteration technique which effectively utilizes solution results from the previous step in determining such results during the current solution step. Also, the computations are limited to the determination of a few specific roots only, which are expected to contain the flutter mode, and this is achieved at each step without having to compute any other root. The structural discretization achieved by the finite-element method yields highly banded stiffness, mass, and aerodynamic matrices; the aerodynamic matrix evaluated by the linearized piston theory is real but unsymmetric in nature. The solution algorithm presented in this paper fully exploits the banded form of the associated matrices, and the resulting computer program written in FORTRAN V for the JPL UNIVAC 1108 computer proves to be most efficient and economical when compared to existing procedures of such analysis. Numerical results are presented for a two-dimensional panel flutter problem.

Gupta, K. K.↗

FCAP - A new tool for the performance and structural analysis for complex flexible aircraft with active control

A very general formulation has been developed for the analysis of complex flexible-aircraft configurations in presence of active control. The formulation is incorporated into a computer program, Flight Control Analysis Program (FCAP), which is designed in a modular fashion to incorporate aircraft dynamics and aerodynamics for complex configurations, as well as sensor, actuator and control logic element dynamics. Formulation of the total aircraft dynamic system is accomplished in matrix form by casting equations in state vector format. The system stability and performance are determined in either the frequency or time domain using classical analysis techniques. Fully unsteady aerodynamics is used for the evaluation of the flutter characteristics as well as the gust response of the aircraft.

Morino, L.↗

A time dependent difference theory for sound propagation in ducts with flow

A time dependent numerical solution of the linearized continuity and momentum equation was developed for sound propagation in a two dimensional straight hard or soft wall duct with a sheared mean flow. The time dependent governing acoustic difference equations and boundary conditions were developed along with a numerical determination of the maximum stable time increments. A harmonic noise source radiating into a quiescent duct was analyzed. This explicit iteration method then calculated stepwise in real time to obtain the transient as well as the steady state solution of the acoustic field. Example calculations were presented for sound propagation in hard and soft wall ducts, with no flow and plug flow. Although the problem with sheared flow was formulated and programmed, sample calculations were not examined. The time dependent finite difference analysis was found to be superior to the steady state finite difference and finite element techniques because of shorter solution times and the elimination of large matrix storage requirements.

Baumeister, K. J.↗

A time dependent difference theory for sound propagation in ducts with flow

A time dependent numerical solution of the linearized continuity and momentum equation is developed for sound propagation in a two-dimensional straight hard or soft wall duct with a sheared mean flow. The time dependent governing acoustic-difference equations and boundary conditions are developed along with a numerical determination of the maximum stable time increments. The analysis begins with a harmonic noise source radiating into a quiescent duct. This explicit iteration method then calculates stepwise in real time to obtain the transient as well as the 'steady' state solution of the acoustic field. Example calculations are presented for sound propagation in hard and soft wall ducts, with no flow and with plug flow. Although the problem with sheared flow has been formulated and programmed, sample calculations have not yet been examined. So far, the time dependent finite difference analysis has been found to be superior to the steady state finite difference and finite element techniques because of shorter solution times and the elimination of large matrix storage requirements.

Baumeister, K. J.↗

Biaxial tests of flat graphite/epoxy laminates

The influence of biaxially applied loads on the strength of composite materials containing holes was analyzed. The analysis was performed through the development of a three dimensional, finite element computer program that is capable of evaluating fiber breakage, delamination, and matrix failure. Realistic failure criteria were established for each of the failure modes, and the influence of biaxial loading on damage accumulation under monotonically increasing loading was examined in detail. Both static and fatigue testing of specially designed biaxial specimens containing central holes were performed. Static tests were performed to obtain an understanding of the influence of biaxial loads on the fracture strength of composite materials and to provide correlation with the analytical predictions. The predicted distributions and types of damage are in reasonable agreement with the experimental results. A number of fatigue tests were performed to determine the influence of cyclic biaxial loads on the fatigue life and residual strength of several composite laminates.

Liebowitz, H.↗

Space applications of Automation, Robotics and Machine Intelligence Systems (ARAMIS). Volume 2: Space projects overview

Applications of automation, robotics, and machine intelligence systems (ARAMIS) to space activities, and their related ground support functions are studied so that informed decisions can be made on which aspects of ARAMIS to develop. The space project breakdowns, which are used to identify tasks ('functional elements'), are described. The study method concentrates on the production of a matrix relating space project tasks to pieces of ARAMIS.

Miller, R. H.↗

Triangular system solutions on an optical systolic processor

It is noted that the simplified system of equations presented by Casasent and Ghosh (1983) can also be solved optically. One realization of this is described, with attention given to data flow and speed for the case of a lower triangular system of equations. The solution of an upper triangular system of equations can be handled in analogous fashion. The problem requires a matrix-vector processor. Only one carrier frequency input to the acoustooptic cell and only one output detector are used. To multiply a matrix by a vector in this system, the time-multiplexed vector elements are fed into the acoustooptic cell.

Ghosh, A.↗

Extended aeroelastic analysis for helicopter rotors with prescribed hub motion and blade appended penduluum vibration absorbers

The mathematical development for the expanded capabilities of the G400 rotor aeroelastic analysis was examined. The G400PA expanded analysis simulates the dynamics of all conventional rotors, blade pendulum vibration absorbers, and the higher harmonic excitations resulting from prescribed vibratory hub motions and higher harmonic blade pitch control. The methodology for modeling the unsteady stalled airloads of two dimensional airfoils is discussed. Formulations for calculating the rotor impedance matrix appropriate to the higher harmonic blade excitations are outlined. This impedance matrix, and the associated vibratory hub loads, are the rotor dynamic characteristic elements for use in the simplified coupled rotor/fuselage vibration analysis (SIMVIB). Updates to the development of the original G400 theory, program documentation, user instructions and information are presented.

Bielawa, R. L.↗

Computational technology for flight vehicles; Proceedings of the Symposium, Washington, DC, Nov. 5-7, 1990

Recent advances in computational fluid mechanics are discussed in reviews and reports. Sections are devoted to (1) the modeling of local phenomena and edge effects in solids, (2) stochastic modeling and simulation of fracture toughness, and (3) partitioning strategy and new finite elements. Particular attention is given to global and local finite-element/spectral-boundary-element techniques for failure analysis, simulations of microfracture in metal-matrix composites, fatigue analysis of cracked anisotropic plates under stochastic loading, mathematical modeling for the analysis of nonlinear aircraft dynamics, physical and mathematical modeling of wave propagation in the Ariane 5 VEB structure, partitioning based on symmetry transformations, an FEM approach to adaptive reliability assurance, and time-domain FEMs for the large rotational dynamics of multibody systems.

Noor, Ahmed K.↗