Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “computational structural dynamics”

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 991 records · Page 55

Development of Unsteady Aerodynamic and Aeroelastic Reduced-Order Models Using the FUN3D Code

Recent significant improvements to the development of CFD-based unsteady aerodynamic reduced-order models (ROMs) are implemented into the FUN3D unstructured flow solver. These improvements include the simultaneous excitation of the structural modes of the CFD-based unsteady aerodynamic system via a single CFD solution, minimization of the error between the full CFD and the ROM unsteady aero- dynamic solution, and computation of a root locus plot of the aeroelastic ROM. Results are presented for a viscous version of the two-dimensional Benchmark Active Controls Technology (BACT) model and an inviscid version of the AGARD 445.6 aeroelastic wing using the FUN3D code.

Silva, Walter A.↗

Investigation of Aerodynamic Matching Sensitivities for Conceptual Design

A methodology for aerodynamic matching was implemented in the Higher fidelity Conceptual Design and structural optimization (HCDstruct) tool to support transonic flutter simulations. The Enhanced Factor Correction Technique was implemented in HCDstruct to correct the Nastran doublet lattice method using data from corresponding computational fluid dynamics simulations. In this paper, sensitivities associated with this aerodynamic matching technique were investigated using HCDstruct and a Boeing 737-200-like aircraft configuration. Specifically, sensitivities of the aerodynamic matching method were investigated with respect to changes in flight conditions, such as changes in Mach number and angle of attack, and also in key design variables for the main wing, such as changes in aspect ratio and sweep angle. The variations in flight conditions and design parameters were mapped to estimated prediction error in sectional lift for the main wing, which helps to demonstrate the extent to which a user may expect a given aerodynamic correction matrix to provide suitable aerodynamic predictions for a realistic aircraft wing.

Jesse R. Quinlan↗

Investigation of Aerodynamic Matching Sensitivities for Conceptual Design

A methodology for aerodynamic matching was implemented in the Higher fidelity Conceptual Design and structural optimization (HCDstruct) tool to support transonic flutter simulations. The Enhanced Factor Correction Technique was implemented in HCDstruct to correct the Nastran doublet lattice method using data from corresponding computational fluid dynamics simulations. In this paper, sensitivities associated with this aerodynamic matching technique were investigated using HCDstruct and a Boeing 737-200-like aircraft configuration. Specifically, sensitivities of the aerodynamic matching method were investigated with respect to changes in flight conditions, such as changes in Mach number and angle of attack, and also in key design variables for the main wing, such as changes in aspect ratio and sweep angle. The variations in flight conditions and design parameters were mapped to estimated prediction error in sectional lift for the main wing, which helps to demonstrate the extent to which a user may expect a given aerodynamic correction matrix to provide suitable aerodynamic predictions for a realistic aircraft wing.

Jesse Ray Quinlan↗

Computational methods for the identification of spatially varying stiffness and damping in beams

A numerical approximation scheme for the estimation of functional parameters in Euler-Bernoulli models for the transverse vibration of flexible beams with tip bodies is developed. The method permits the identification of spatially varying flexural stiffness and Voigt-Kelvin viscoelastic damping coefficients which appear in the hybrid system of ordinary and partial differential equations and boundary conditions describing the dynamics of such structures. An inverse problem is formulated as a least squares fit to data subject to constraints in the form of a vector system of abstract first order evolution equations. Spline-based finite element approximations are used to finite dimensionalize the problem. Theoretical convergence results are given and numerical studies carried out on both conventional (serial) and vector computers are discussed.

Banks, H. T.↗

Simultaneous expansion and orthogonalization of measured modes for structure identification

Tests of large structures on-orbit will be performed with measurements at a relatively few structure points. Values for the unmeasured degrees of freedom (dofs) can be estimated based on measured dofs and analytical model dynamic information. These 'expanded' mode shapes are useful for optimal-update identification and damage location as well as test/analysis correlation. A new method of expansion for test mode shape vectors is developed from the orthogonal Procrustes problem from computational linear algebra. A subspace defined by the set of measured dofs is compared to a subspace defined by mode shapes from an analytical model of the structure. The method simultaneously expands and orthogonalizes the mode shape vectors. Two demonstration problems are used to compare the new method to current expansion techniques. One demonstration uses test data from a laboratory scale-model truss structure. Performance of the new method is comparable or superior to that of the previous expansion methods which require separate orthogonalization.

Smith, Suzanne Weaver↗

Numerical simulation of compressor endwall and casing treatment flow phenomena

A numerical study is presented of the flow in the endwall region of a compressor blade row, in conditions of operation with both smooth and grooved endwalls. The computations are first compared to velocity field measurements in a cantilevered stator/rotating hub configuration to confirm that the salient features are captured. Computations are then interrogated to examine the tip leakage flow structure since this is a dominant feature of the endwall region. In particular, the high blockage that can exist near the endwalls at the rear of a compressor blade passage appears to be directly linked to low total pressure fluid associated with the leakage flow. The fluid dynamic action of the grooved endwall, representative of the casing treatments that have been most successful in suppressing stall, is then simulated computationally and two principal effects are identified. One is suction of the low total pressure, high blockage fluid at the rear of the passage. The second is energizing of the tip leakage flow, most notably in the core of the leakage vortex, thereby suppressing the blockage at its source.

Crook, A. J.↗

Reduced-Order Aerodynamic Modeling Based on CFD Frequency Responses from Multisine Inputs

A system identification analysis was performed to determine reduced-order models of a computational fluid dynamics (CFD) solver for linear aeroelastic analysis and control design. The application was to the FUN3D code and the flexible half-span wind tunnel test article, in transonic flow conditions, used in the NASA-Boeing collaboration called the Integrated Adaptive Wing Technology Maturation (IAWTM) project. Multiple inputs (structural mode displacements and control surface deflections) were simultaneously excited with orthogonal phase-optimized multisines and multiple outputs (generalized aerodynamic forces) were recorded, from which the matrix of frequency responses were computed using a single CFD run. A state-space model was then fit to the frequency response data using a maximum-likelihood estimator.

System identification↗

Structure and Dynamics of Bleed-Controlled Impinging Shock/Turbulent-Boundary-Layer Interactions

Porous boundary-layer bleed is often employed to alleviate the detrimental effects of shock/turbulent-boundary-layer interactions (STBLI) in high-speed vehicles. However, outstanding questions remain on the underlying mechanisms through which bleed influences STBLI dynamics. This work addresses this knowledge gap by performing a scale-resolved Large-Eddy Simulation of an impinging STBLI controlled by discrete bleed holes. The simulation has a freestream Mach number of 2.5, an oblique shock generator angle of 8 degrees, and a Reynolds number based on the incoming momentum thickness of 2800. Bleed substantially reduces the separation size and upstream influence of the STBLI compared to the canonical, uncontrolled case. Consequently, the characteristic, low-frequency motion of STBLI-induced separation is mitigated, as revealed by wall-pressure spectra. Energy is shifted from the low frequency separation behavior to a higher-frequency signature caused by the bleed shocks and expansions. The scale-resolved simulations also reveal that these small, localized structures severely limit the turbulence production traditionally seen in the uncontrolled interaction.

computational fluid dynamics↗

NERVA dynamic analysis methodology, SPRVIB

The general dynamic computer code called SPRVIB (Spring Vib) developed in support of the NERVA (nuclear engine for rocket vehicle application) program is described. Using normal mode techniques, the program computes kinematical responses of a structure caused by various combinations of harmonic and elliptic forcing functions or base excitations. Provision is made for a graphical type of force or base excitation input to the structure. A description of the required input format and a listing of the program are presented, along with several examples illustrating the use of the program. SPRVIB is written in FORTRAN 4 computer language for use on the CDC 6600 or the IBM 360/75 computers.

Vronay, D. F.↗

Magnetic field and electric current structure in the chromosphere

The three-dimensional vector magnetic field structure in the chromosphere above an active region is deduced by using high-resolution H-alpha filtergrams together with a simultaneous digital magnetogram. An analog model of the field is made with 400 metal wires representing field lines that outline the H-alpha structure. The height extent of the field is determined from vertical field-gradient observations around sunspots, from observed fibril heights, and from an assumption that the sources of the field are largely local. The computed electric currents (typically 10 mA/sq m) are found to flow in patterns not similar to observed features and not parallel to magnetic fields. Force structures correspond to observed solar features; the dynamics to be expected include: downward motion in bipolar areas in the lower chromosphere, an outflow of the outer chromosphere into the corona with radially outward flow above bipolar plage regions, and motion of arch filament systems.

Dravins, D.↗

Dynamic characterization of solid rockets

The structural dynamics of solid rockets in-general was studied. A review is given of the modes of vibration and bending that can exist for a solid propellant rocket, and a NASTRAN computer model is included. Also studied were the dynamic properties of a solid propellant, polybutadiene-acrylic acid-acrylonitrile terpolymer, which may be used in the space shuttle rocket booster. The theory of viscoelastic materials (i.e, Poisson's ratio) was employed in describing the dynamic properties of the propellant. These studies were performed for an eventual booster stage development program for the space shuttle.

Source record↗

Aeronautical propulsion: Present status and future directions

The advancement of aeropropulsion systems continues to provide technology to various portions of the gas turbine field. It is recognized that this area is undergoing considerale change, which will result in substantially improved gas turbine components and systems. These changes are occurring in a number of technical areas including advanced analytical and physical measurement methods, the application of large scientific computers, the dynamic modeling of components and systems, the application of integrated control systems that optimize and improve performance and system condition monitoring, and the development of new and unique materials and structures. As these areas evolve, the ways in which technology will advance, and factors affecting the design and development of new systems, will probably be considerably different than those of today. It is also anticipated that the necessary skilled work force will be different. Certainly there will be changes, but the nature, extent, and rate of those changes can only be surmised at this time.

Hartmann, M. J.↗

Decoupled dynamic analysis of combined systems by iterative determination of interface accelerations

A dynamic analysis technique is presented that can be used to determine the response of a discrete model of a large linear structural system composed of multiple substructures. The technique circumvents the costly computation of the modal characteristics of the combined system. This is accomplished by relying on a predictor-corrector scheme to converge iteratively to the interface accelerations of the combined system, while the equations of motions of the individual structures are integrated separately. In this regard, the temporal slopes of the interface accelerations (jerks) are computed at each time point of integration to predict the interface accelerations at the next time point. The proposed technique is exemplified by conducting a Space Shuttle landing loads analysis; the obtained numerical data demonstrate its reliability and efficiency.

Spanos, P. D.↗

Probabilistic Fatigue Damage Program (FATIG)

FATIG computes fatigue damage/fatigue life using the stress rms (root mean square) value, the total number of cycles, and S-N curve parameters. The damage is computed by the following methods: (a) traditional method using Miner s rule with stress cycles determined from a Rayleigh distribution up to 3*sigma; and (b) classical fatigue damage formula involving the Gamma function, which is derived from the integral version of Miner's rule. The integration is carried out over all stress amplitudes. This software solves the problem of probabilistic fatigue damage using the integral form of the Palmgren-Miner rule. The software computes fatigue life using an approach involving all stress amplitudes, up to N*sigma, as specified by the user. It can be used in the design of structural components subjected to random dynamic loading, or by any stress analyst with minimal training for fatigue life estimates of structural components.

Michalopoulos, Constantine↗

Current Computational Challenges for CMC Processes, Properties, and Structures

In comparison to current state-of-the-art metallic alloys, ceramic matrix composites (CMC) offer a variety of performance advantages, such as higher temperature capability (greater than the approx.2100 F capability for best metallic alloys), lower density (approx.30-50% metal density), and lower thermal expansion. In comparison to other competing high-temperature materials, CMC are also capable of providing significantly better static and dynamic toughness than un-reinforced monolithic ceramics and significantly better environmental resistance than carbon-fiber reinforced composites. Because of these advantages, NASA, the Air Force, and other U.S. government agencies and industries are currently seeking to implement these advanced materials into hot-section components of gas turbine engines for both propulsion and power generation. For applications such as these, CMC are expected to result in many important performance benefits, such as reduced component cooling air requirements, simpler component design, reduced weight, improved fuel efficiency, reduced emissions, higher blade frequencies, reduced blade clearances, and higher thrust. Although much progress has been made recently in the development of CMC constituent materials and fabrication processes, major challenges still remain for implementation of these advanced composite materials into viable engine components. The objective of this presentation is to briefly review some of those challenges that are generally related to the need to develop physics-based computational approaches to allow CMC fabricators and designers to model (1) CMC processes for fiber architecture formation and matrix infiltration, (2) CMC properties of high technical interest such as multidirectional creep, thermal conductivity, matrix cracking stress, damage accumulation, and degradation effects in aggressive environments, and (3) CMC component life times when all of these effects are interacting in a complex stress and service environment. To put these computational issues in perspective, the various modeling needs within these three areas are briefly discussed in terms of their technical importance and their key controlling mechanistic factors as we know them today. Emphasis is placed primarily on the SiC/SiC ceramic composite system because of its higher temperature capability and enhanced development within the CMC industry. A brief summary is then presented concerning on-going property studies aimed at addressing these CMC modeling needs within NASA in terms of their computational approaches and recent important results. Finally an overview perspective is presented on those key areas where further CMC computational studies are needed today to enhance the viability of CMC structural components for high-temperature applications.

DiCarlo, James↗

The LATDYN user's manual

The LATDYN User's Manual presents the capabilities and instructions for the LATDYN (Large Angle Transient DYNamics) computer program. The LATDYN program is a tool for analyzing the controlled or uncontrolled dynamic transient behavior of interconnected deformable multi-body systems which can undergo large angular motions of each body relative other bodies. The program accommodates large structural deformation as well as large rigid body rotations and is applicable, but not limited to, the following areas: (1) development of large flexible space structures; (2) slewing of large space structure components; (3) mechanisms with rigid or elastic components; and (4) robotic manipulations of beam members. Presently the program is limited to two dimensional problems, but in many cases, three dimensional problems can be exactly or approximately reduced to two dimensions. The program uses convected finite elements to affect the large angular motions involved in the analysis. General geometry is permitted. Detailed user input and output specifications are provided and discussed with example runstreams. To date, LATDYN has been configured for CDC/NOS and DEC VAX/VMS machines. All coding is in ANSII-77 FORTRAN. Detailed instructions regarding interfaces with particular computer operating systems and file structures are provided.

Housner, J. M.↗

Multi-Domain Routing in Delay Tolerant Networks

The goal of Delay Tolerant Networking (DTN) is to provide the missing ingredient for the ever-growing collection of communicating nodes in our solar system to become a Solar System Internet (SSI). Great strides have been made in modeling particular types of DTNs, such as schedule- or discovery-based. Now, analogously to the Internet, these smaller DTNs can be considered routing domains which must be stitched together to form the overall SSI. In this paper, we propose a framework for cross-domain routing in DTNs as well as methodologies for detecting these sub-domains. Example time-varying networks are given to demonstrate the techniques proposed. A basic component is the mathematical theory of sheaves, which unifies the underlying model of DTN routing algorithms, by giving rise to routing sheaves – these can be defined for the dynamic and scheduled networks as noted above, and can also be used to define the interfaces between these domains in order to route across them. An immediate application would be routing across discovery-based networks connected by scheduled networks. These DTN subdomains remain elusive, however, and need to become well-defined and properly sized for tractable computability. In particular, a balance must be determined between areas that are too large (i.e. large matrix computations) versus areas that are too small (i.e. “many” single-noded domains). Moreover, the connections between the domains should, at least locally, be chosen to optimize data flow and connectivity: we address this in three ways. First, tools from persistent homology are given to understand underlying structures, reminiscent of hierarchies in the Internet Protocol (IP) addressing. Second, we construct a notion of temporal graph curvature based on network geometry to analyze flows induced by dynamical processes on these networks. Finally, Schrodinger Bridges, a tool arising from statistical physics, are proposed as a method of constructing flows on time-evolving networks with desirable properties such as speed, robustness, and load sensitivity. We construct an approach to temporal hypergraphs to simultaneously model unicast, multicast, and broadcast, using the language of scheme theory, and then consider DTN network coding as a way to achieve network-level computation and organization. The paper concludes with a discussion and ideas for future work.

Alan Hylton↗

Overview of the Predictive Simulation Capability Element of the Plume Surface Interaction Project

As part of the Game Changing Development (GCD) Program, funded by NASA’s Space Technology Mission Directorate (STMD), the development of simulation capability for the prediction of extra-terrestrial Plume Surface Interaction (PSI) environments has been undertaken by the Fluid Dynamics Branch at NASA/MSFC. The GCD PSI Project, planned to be completed over a four year period, contains a Predictive Simulation Capability (PSC) Element focused on creating simulation capability for the reliable and accurate prediction of PSI in Martian (~650 Pa) and Lunar (vacuum) ambient environments. In addition to the PSC Element, the GCD Program also contains a companion Ground Testing Element for development of focused datasets for validation of predictive capability as well as a Flight-focused Instrumentation Element. This paper describes the PSC Element of the PSI Project as well as providing descriptions of recent accomplishments and remaining work. The overall structure of the PSC Element is broken down into four areas of focus. The first area is the Prediction of Plume Flow in low pressure environments. The development approach taken is the augmentation of the existing production-mode computational fluid dynamics (CFD) tool Loci/Chem, with targeted extensions necessary to accurately model rarefied conditions found in both Martian and Lunar applications. Production readiness and validation of predictive capability are the major objectives of this task. The second area is the effect of mixed continuum/rarefied flow on crater development and ejecta sheets. A new CFD application, Loci/GGFS (Gas Granular Flow Solver), is being developed which implements an Eulerian/Eulerian two-phase model of gas- and soil-phases in order to simulate the soil erosion, crater formation, and soil ejecta transport in a fully coupled simulation. This task seeks to verify Loci/GGFS is production-ready as well as perform validation studies to determine the degree of predictive capability achieved by Loci/GGFS. The third area is focused on the details and extension of particle phase modeling of soil. In this task, Discrete Element Modeling (DEM) techniques are used to perform direct simulations of complex soil particles under the action of forcing similar to that to be cause by PSI. The simulation results are then used to construct closures to the Eulerian model of the soil phase used by Loci/GGFS. The fourth area is gas-particle interaction modeling. In this task, experiments are being conducted as well as detailed simulation results are being studied to further understand the complexities of gas-particle interactions in dilute, intermediate, and high soil volume fraction regimes. Improved models of particle drag and the particle turbulent kinetic energy (PTKE) resulting from the interaction of gas flows within particle clouds are the objective of this task.

Jeff West↗