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 1,531 records · Page 85

Atomic-scale mechanism of carbon nucleation from a deep crustal fluid by replica exchange reactive molecular dynamics simulation

Here we present a mechanistic model of carbon nucleation and growth from a fluid at elevated temperature (T) and pressure conditions, typical of those found in the shallow Earth’s lithosphere. Our model uses a replica exchange reactive molecular dynamics framework in which molecular configurations are swapped between adjacent T replica at regular intervals according to underlying statistical mechanics. This framework allows predicting complex molecular structures and thermodynamics while remaining computationally efficient. Here we simulate the reactivity of an unstable mixture of CO 2 and CH 4 at 1000 K and 1 GPa. We find that the path to thermodynamic equilibrium is initially entropy-driven, producing a diversity of short-lived species, including various alcohols with intermediate carbon oxidation states. Cyclic and polycyclic radicals that are sometimes resonance-stabilized form next and set the stage for carbon nucleation. The carbon exsolution process releases abundant water, is exothermic and starts with the nucleation of a large aggregate of hydrogenated graphene flakes from covalently bonded polycyclic units. The carbon backbone of this nucleus subsequently grows into a hydrogen-depleted fullerene-like structure, before evolving toward a partially bilayered graphene layer. Overall, our results show that the mechanism of graphitic C formation is certainly not bimolecular, and that it may involve a combination of key condensation and radical chain reactions. This will help understand the isotopic, and reactive characteristics of carbon-bearing fluids during their upward transit through the Earth’s mantle and crust. Moreover, the mechanistic insights outlined here present intriguing similarities with the process of soot and interstellar dust formation, which suggests that the widespread distribution of abiotic polyaromatic and graphitic material on Earth and beyond may reflect the prevalence of a fundamental chemical pathway.

58 GEOSCIENCES↗

Gas Phase Synthesis of Benzene via the Propargyl Radical Self-Reaction

Polycyclic Aromatic Hydrocarbons (PAHs) have been invoked in fundamental molecular mass growth processes in our galaxy. We provide compelling evidence of the formation of the very first ringed aromatic and building block of PAHs - benzene - via the self-recombination of two resonantly stabilized propargyl (C 3 H 3 ) radicals in dilute environments using isomer selective synchrotron-based mass spectrometry coupled to theoretical calculations. Along with benzene, three other structural isomers (1,5-hexadiyne, fulvene, 2-ethynyl-1,3-butadiene) and o-benzyne are detected and their branching ratios quantified experimentally and verified with the aid of computational fluid dynamics and kinetic simulations. Furthermore, these results uncover molecular growth pathways not only in interstellar, circumstellar, and solar systems environments, but also in combustion systems gaining us a better understanding of the hydrocarbon chemistry of our universe.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

VISCEL: A general-purpose computer program for analysis of linear viscoelastic structures (user's manual), volume 1

This program, an extension of the linear equilibrium problem solver ELAS, is an updated and extended version of its earlier form (written in FORTRAN 2 for the IBM 7094 computer). A synchronized material property concept utilizing incremental time steps and the finite element matrix displacement approach has been adopted for the current analysis. A special option enables employment of constant time steps in the logarithmic scale, thereby reducing computational efforts resulting from accumulative material memory effects. A wide variety of structures with elastic or viscoelastic material properties can be analyzed by VISCEL. The program is written in FORTRAN 5 language for the Univac 1108 computer operating under the EXEC 8 system. Dynamic storage allocation is automatically effected by the program, and the user may request up to 195K core memory in a 260K Univac 1108/EXEC 8 machine. The physical program VISCEL, consisting of about 7200 instructions, has four distinct links (segments), and the compiled program occupies a maximum of about 11700 words decimal of core storage.

Gupta, K. K.↗

CELFE: Coupled Eulerian-Lagrangian Finite Element program for high velocity impact. Part 1: Theory and formulation

A 3-D finite element program capable of simulating the dynamic behavior in the vicinity of the impact point, together with predicting the dynamic response in the remaining part of the structural component subjected to high velocity impact is discussed. The finite algorithm is formulated in a general moving coordinate system. In the vicinity of the impact point contained by a moving failure front, the relative velocity of the coordinate system will approach the material particle velocity. The dynamic behavior inside the region is described by Eulerian formulation based on a hydroelasto-viscoplastic model. The failure front which can be regarded as the boundary of the impact zone is described by a transition layer. The layer changes the representation from the Eulerian mode to the Lagrangian mode outside the failure front by varying the relative velocity of the coordinate system to zero. The dynamic response in the remaining part of the structure described by the Lagrangian formulation is treated using advanced structural analysis. An interfacing algorithm for coupling CELFE with NASTRAN is constructed to provide computational capabilities for large structures.

Lee, C. H.↗

Structural dynamics of shroudless, hollow fan blades with composite in-lays

Structural and dynamic analyses are presented for a shroudless, hollow titanium fan blade proposed for future use in aircraft turbine engines. The blade was modeled and analyzed using the composite blade structural analysis computer program (COBSTRAN); an integrated program consisting of mesh generators, composite mechanics codes, NASTRAN, and pre- and post-processors. Vibration and impact analyses are presented. The vibration analysis was conducted with COBSTRAN. Results show the effect of the centrifugal force field on frequencies, twist, and blade camber. Bird impact analysis was performed with the multi-mode blade impact computer program. This program uses the geometric model and modal analysis from the COBSTRAN vibration analysis to determine the gross impact response of the fan blades to bird strikes. The structural performance of this blade is also compared to a blade of similar design but with composite in-lays on the outer surface. Results show that the composite in-lays can be selected (designed) to substantially modify the mechanical performance of the shroudless, hollow fan blade.

Aiello, R. A.↗

Conceptual design and analysis of a large antenna utilizing electrostatic membrane management

Conceptual designs and associated technologies for deployment 100 m class radiometer antennas were developed. An electrostatically suspended and controlled membrane mirror and the supporting structure are discussed. The integrated spacecraft including STS cargo bay stowage and development were analyzed. An antenna performance evaluation was performed as a measure of the quality of the membrane/spacecraft when used as a radiometer in the 1 GHz to 5 GHz region. Several related LSS structural dynamic models differing by their stiffness property (and therefore, lowest modal frequencies) are reported. Control system whose complexity varies inversely with increasing modal frequency regimes are also reported. Interactive computer-aided-design software is discussed.

Brooks, A. L.↗

Computing the transmission zeros of large space structures

The transmission zeros of a large space structure can be computed by the general-purpose algorithm of A. Emami-Naeini and P. Van Dooren (1982). However, careful use of the special form of the equations of motion of structural dynamics leads to a new method that is about twice as fast as theirs when applied to a damped structure, and at least 60 times as fast for an undamped one.

Williams, Trevor↗

Beam rider for an Articulated Robot Manipulator (ARM) accurate positioning of long flexible manipulators

Laser beam positioning and beam rider modules were incorporated into the long hollow flexible segment of an articulated robot manipulator (ARM). Using a single laser beam, the system determined the position of the distal ARM endtip, with millimetric precision, in six degrees of freedom, at distances of up to 10 meters. Preliminary designs, using space rated technology for the critical systems, of a two segmented physical ARM, with a single and a dual degree of freedom articulation, were developed, prototyped, and tested. To control the positioning of the physical ARM, an indirect adaptive controller, which used the mismatch between the position of the laser beam under static and dynamic conditions, was devised. To predict the behavior of the system and test the concept, a computer simulation model was constructed. A hierarchical artificially intelligent real time ADA operating system program structure was created. The software was designed for implementation on a dedicated VME bus based Intel 80386 administered parallel processing multi-tasking computer system.

Malachowski, M. J.↗

Freezing in confined geometries

Results of detailed structural studies, using elastic neutron scattering, of the freezing of liquid O2 and D2 in porous vycor glass, are presented. The experimental studies have been complemented by computer simulations of the dynamics of freezing of a Lennard-Jones liquid in narrow channels bounded by molecular walls. Results point to a new simple physical interpretation of freezing in confined geometries.

Sokol, P. E.↗

Direct adaptive impedance control of manipulators

An adaptive scheme for controlling the end-effector impedance of robot manipulators is presented. The proposed control system consists of three subsystems: a simple filter which characterizes the desired dynamic relationship between the end-effector position error and the end-effector/environment contact force, an adaptive controller which produces the Cartesian-space control input required to provide this desired dynamic relationship, and an algorithm for mapping the Cartesian-space control input to a physically realizable joint-space control torque. The controller does not require knowledge of either the structure or the parameter values of the robot dynamics, and it is implemented without calculation of the robot inverse kinematic transformation. As a result, the scheme represents a very general and computationally efficient approach to controlling the impedance of both nonredundant and redundant manipulators. Furthermore, the method can be applied directly to trajectory tracking in free-space motion by removing the impedance filter.

Colbaugh, R.↗

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↗

Tutorial for Non-Orthogonal Configuration Interaction Fragment Studies Using GronOR

A brief tutorial is presented for the GronOR non-orthogonal configuration interaction application illustrating the download and build process, setting up calculations with OpenMolcas and GronOR, and a description of the file structure and input options. An overview of the computational workflow includes the generation of fragment wave functions and integrals with an adapted version of OpenMolcas, running GronOR calculations, and accounting for dynamic correlation corrections. Three end-to-end use cases are provided, representing a molecular dimer, an ion pair, and overlapping molecular fragments.

Straatsma, T.P.↗

Discovering causal structure with reproducing-kernel Hilbert space ε -machines

We merge computational mechanics’ definition of causal states (predictively equivalent histories) with reproducing-kernel Hilbert space (RKHS) representation inference. The result is a widely applicable method that infers causal structure directly from observations of a system’s behaviors whether they are over discrete or continuous events or time. A structural representation—a finite- or infinite-state kernel ϵ-machine—is extracted by a reduced-dimension transform that gives an efficient representation of causal states and their topology. In this way, the system dynamics are represented by a stochastic (ordinary or partial) differential equation that acts on causal states. We introduce an algorithm to estimate the associated evolution operator. Paralleling the Fokker–Planck equation, it efficiently evolves causal-state distributions and makes predictions in the original data space via an RKHS functional mapping. We demonstrate these techniques, together with their predictive abilities, on discrete-time, discrete-value infinite Markov-order processes generated by finite-state hidden Markov models with (i) finite or (ii) uncountably infinite causal states and (iii) continuous-time, continuous-value processes generated by thermally driven chaotic flows. The method robustly estimates causal structure in the presence of varying external and measurement noise levels and for very high-dimensional data.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

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.↗

Constant-potential molecular dynamics simulations of molten salt double layers for FLiBe and FLiNaK

We report the results of constant-potential molecular dynamics simulations of the double layer interface between molten 2LiF–BeF 2 (FLiBe) and 23LiF–6NaF–21KF (FLiNaK) fluoride mixtures and idealized solid electrodes. Employing methods similar to those used in studies of chloride double layers, we compute the structure and differential capacitance of molten fluoride electric double layers as a function of applied voltage. The role of molten salt structure is probed through comparisons between FLiBe and FLiNaK, which serve as models for strong and weak associate-forming salts, respectively. In FLiBe, screening involves changes in Be–F–Be angles and alignment of the oligomers parallel to the electrode, while in FLiNaK, the electric field is screened mainly by rearrangement of individual ions, predominantly the polarizable potassium cation.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗