SEARCH · Engineering Papers
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.
A computer program for redesigning structural modes to reduce response has been initiated. The linear regulator approach in modal coordinates has been implemented. It is noted that the transformation of solution to physical structure is a major problem. It is concluded that the solution of stiffness equations and damping equations can be done separately as NXN set of (matrix Riccati) equations.
An 'electrograph' instrument designed for measurement of macroscopic electric fields in solar plasmas is described, using the polarization dependence of line width in Stark-broadened hydrogen Paschen emission lines. Observations of quiescent prominences and limb chromosphere with the electrograph at the NSO/Sac Peak Evans Coronal Facility provide upper limits of 5-10 V/cm for transverse macroscopic electric fields in these structures, averaged over an area of about 5 x 7 arcsec. Random thermal motions of hydrogen ions across magnetic field lines generate a quasi-static electric field, which should be distinguishable from pressure broadening in the intensely magnetized chromosphere over a sunspot, given an electrograph sensitivity a factor 2-3 better than that achieved here. Future electrograph measurements of limb flares, post-flare loops, and eruptive prominences, even at 5-V/cm sensitivity, could provide a useful new test of reconnection and discharge effects in such dynamic structures.
In this paper we investigate the thermal and dynamic properties of dynamic structures in and around a prominence channel observed on the limb on 17 April 2003. Observations were taken with the Solar and Heliospheric Observatory's Solar Ultraviolet Measurements of Emitted Radiation (SOHO/SUMER) in lines formed at temperatures from 80,000 to 1.6 MK. The instrument was pointed to a single location and took a series of 90 s exposures. Two-dimensional context was provided by the Transition Region and Coronal Explorer (TRACE) in the UV and EUV and the Kanzelhohe Solar Observatory in H-alpha. Two dynamic features were studied in depth: an activated prominence and repeated motions in a loop near the prominence. We calculated three-dimensional geometries and trajectories, differential emission measure, and limits on the mass, pressure, average density, and kinetic and thermal energies. These observations provide important tests for models of dynamics in prominences and cool (approx. 10(exp 5) K)loops, which will ultimately lead to a better understanding the mechanism(s) leading to energy and mass flow in these solar features.
This paper reports on an initial assessment of using a Field-Programmable Gate Array (FPGA) computational device as a new tool for solving structural mechanics problems. A FPGA is an assemblage of binary gates arranged in logical blocks that are interconnected via software in a manner dependent on the algorithm being implemented and can be reprogrammed thousands of times per second. In effect, this creates a computer specialized for the problem that automatically exploits all the potential for parallel computing intrinsic in an algorithm. This inherent parallelism is the most important feature of the FPGA computational environment. It is therefore important that if a problem offers a choice of different solution algorithms, an algorithm of a higher degree of inherent parallelism should be selected. It is found that in structural analysis, an 'analog computer' style of programming, which solves problems by direct simulation of the terms in the governing differential equations, yields a more favorable solution algorithm than current solution methods. This style of programming is facilitated by a 'drag-and-drop' graphic programming language that is supplied with the particular type of FPGA computer reported in this paper. Simple examples in structural dynamics and statics illustrate the solution approach used. The FPGA system also allows linear scalability in computing capability. As the problem grows, the number of FPGA chips can be increased with no loss of computing efficiency due to data flow or algorithmic latency that occurs when a single problem is distributed among many conventional processors that operate in parallel. This initial assessment finds the FPGA hardware and software to be in their infancy in regard to the user conveniences; however, they have enormous potential for shrinking the elapsed time of structural analysis solutions if programmed with algorithms that exhibit inherent parallelism and linear scalability. This potential warrants further development of FPGA-tailored algorithms for structural analysis.
We will discuss the thermal and dynamic properties of dynamic structures in and around a prominence channel observed on the limb on 17 April 2003. Observations were taken with the Solar and Heliospheric Observatory's Solar Ultraviolet Measurements of Emitted Radiation (SOHO/SUMER) in lines formed at temperatures from 80,000 to 1.6 million K. The instrument was pointed to a single location and took a series of 90 s exposures. Two-dimensional context was provided by the Transition Region and Coronal Explorer (TRACE) in the UV and EUV and by the Kanzelhohe Solar Observatory in H-alpha. Two dynamic features were studied in depth: a prominence activation and repeated motions in a loop nearby the prominence. We were able to calculate three-dimensional geometries and trajectories, differential emission measure, and limits on the average density, kinetic and thermal energies. These observations provide important tests for models of dynamics in prominences and cool (approx.10(exp 5) K) loops.
The turbine under study incorporates the combination of two features: the application of blade cyclic pitch variation; and the use of yaw angle control for rotor speed and torque regulation. Due to its emasculation by passive cyclic pitch variation the rotor can be rapidly yawed without encountering gyroscopic and aerodynamic hub moments and without noticeable out of plane excursions. The two bladed upwind rotor is vane stabilized and of very simple and rugged design. The principle was first checked out with a small scale wind tunnel model and then tested in the atmosphere with a 7.6 meter diameter experimental fully instrumented wind turbine driving a 3 phase alternator. The test results are summarized with respect to structural dynamics and yaw dynamics.
The addition of a dynamic loads computation capability to the Unified Flexible Spacecraft Simulation (UFSS) program is discussed. The added capability provides a means for determining the internal member loads due to the time-varying external loading conditions experienced by an orbiting spacecraft/cluster. The flexible bodies are modeled as a system of joints or nodes which are interconnected by weightless finite element members. All masses are lumped at the joints. The orthogonal functions used to describe the spatial deformation of the bodies are normally taken to be the orthonormal cantilever modes produced by a standard structural dynamics program.
Condensate storage tanks (CSTs) in nuclear power plants (NPPs) are classified as critical equipment capable of surviving strong shaking in a design basis earthquake to assure the ability to subsequently provide cooling water. Structural modeling and dynamic analysis of CSTs are complicated due to fluid-structure interaction (FSI) and coupling issues, which could have significant effects on seismic response. The capability of a number of approximate models to capture seismic response is investigated using dynamic analysis results from detailed finite element models. Modal analyses and time history analyses are carried out using both 2D and 3D models to predict the dynamic behavior of CSTs as a function of ground motion intensity. Although 2D simplified models can be used to quickly evaluate the dynamic response of CSTs when there is substantial margin to failure, a detailed 3D model is required when it is necessary to examine a limit state associated with the failure mode of the tank. The findings show good agreement between the natural frequencies of convective modes determined from the 2D and 3D CST models. The difference between the frequencies of impulsive modes were found to be equal to or less than 25%.
Carbon dioxide electroreduction facilitates the sustainable synthesis of fuels and chemicals. Although Cu enables CO 2 -to-multicarbon product (C 2+ ) conversion, the nature of the active sites under operating conditions remains elusive. Importantly, identifying active sites of high-performance Cu nanocatalysts necessitates nanoscale, time-resolved operando techniques. Here, we present a comprehensive investigation of the structural dynamics during the life cycle of Cu nanocatalysts. A 7 nm Cu nanoparticle ensemble evolves into metallic Cu nanograins during electrolysis before complete oxidation to single-crystal Cu 2 O nanocubes following post-electrolysis air exposure. Operando analytical and four-dimensional electrochemical liquid-cell scanning transmission electron microscopy shows the presence of metallic Cu nanograins under CO 2 reduction conditions. Correlated high-energy-resolution time-resolved X-ray spectroscopy suggests that metallic Cu, rich in nanograin boundaries, supports undercoordinated active sites for C–C coupling. Quantitative structure–activity correlation shows that a higher fraction of metallic Cu nanograins leads to higher C 2+ selectivity. A 7 nm Cu nanoparticle ensemble, with a unity fraction of active Cu nanograins, exhibits sixfold higher C 2+ selectivity than the 18 nm counterpart with one-third of active Cu nanograins. Importantly, the correlation of multimodal operando techniques serves as a powerful platform to advance our fundamental understanding of the complex structural evolution of nanocatalysts under electrochemical conditions.
The basic threat to the more widespread use of composites in aircraft primary structure is the problem of residual strength in the presence of damage. Specific problem areas that need to be addressed include the following: (1) modeling/understanding composite failure (impact: relate strain energy release rate to strain-to-failure, open hole compression: relate shear crippling to strain-to-failure, and account for variables such as thickness and stacking sequence); (2) micromechanics models; (3) consider the use of hybrid combinations, both interply and intraply, using graphite, Kevlar, and glass; and (4) modeling/understanding the role of the interface. Solving these problems involves a multidisciplinary approach including dynamics, structural stability, composite mechanics, and fracture mechanics.
Incompressible magnetohydrodynamic (MHD) turbulence and magnetic dynamos, which occur in magnetofluids with large fluid and magnetic Reynolds numbers, will be discussed. When Reynolds numbers are large and energy decays slowly, the distribution of energy with respect to length scale becomes quasi-stationary and MHD turbulence can be described statistically. In the limit of infinite Reynolds numbers, viscosity and resistivity become zero and if these values are used in the MHD equations ab initio, a model system called ideal MHD turbulence results. This model system is typically confined in simple geometries with some form of homogeneous boundary conditions, allowing for velocity and magnetic field to be represented by orthogonal function expansions. One advantage to this is that the coefficients of the expansions form a set of nonlinearly interacting variables whose behavior can be described by equilibrium statistical mechanics, i.e., by a canonical ensemble theory based on the global invariants (energy, cross helicity and magnetic helicity) of ideal MHD turbulence. Another advantage is that truncated expansions provide a finite dynamical system whose time evolution can be numerically simulated to test the predictions of the associated statistical mechanics. If ensemble predictions are the same as time averages, then the system is said to be ergodic; if not, the system is nonergodic. Although it had been implicitly assumed in the early days of ideal MHD statistical theory development that these finite dynamical systems were ergodic, numerical simulations provided sufficient evidence that they were, in fact, nonergodic. Specifically, while canonical ensemble theory predicted that expansion coefficients would be (i) zero-mean random variables with (ii) energy that decreased with length scale, it was found that although (ii) was correct, (i) was not and the expected ergodicity was broken. The exact cause of this broken ergodicity was explained, after much investigation, by greatly extending the statistical theory of ideal MHD turbulence. The mathematical details of broken ergodicity, in fact, give a quantitative explanation of how coherent structure, dynamic alignment and force-free states appear in turbulent magnetofluids. The relevance of these ideal results to real MHD turbulence occurs because broken ergodicity is most manifest in the ideal case at the largest length scales and it is in these largest scales that a real magnetofluid has the least dissipation, i.e., most closely approaches the behavior of an ideal magnetofluid. Furthermore, the effects grow stronger when cross and magnetic helicities grow large with respect to energy, and this is exactly what occurs with time in a real magnetofluid, where it is called selective decay. The relevance of these results found in ideal MHD turbulence theory to the real world is that they provide at least a qualitative explanation of why confined turbulent magnetofluids, such as the liquid iron that fills the Earth's outer core, produce stationary, large-scale magnetic fields, i.e., the geomagnetic field. These results should also apply to other planets as well as to plasma confinement devices on Earth and in space, and the effects should be manifest if Reynolds numbers are high enough and there is enough time for stationarity to occur, at least approximately. In the presentation, details will be given for both theoretical and numerical results, and references will be provided.
Incompressible magnetohydrodynamic (MHD) turbulence and magnetic dynamos, which occur in magnetofluids with large fluid and magnetic Reynolds numbers, will be discussed. When Reynolds numbers are large and energy decays slowly, the distribution of energy with respect to length scale becomes quasi-stationary and MHD turbulence can be described statistically. In the limit of infinite Reynolds numbers, viscosity and resistivity become zero and if these values are used in the MHD equations ab initio, a model system called ideal MHD turbulence results. This model system is typically confined in simple geometries with some form of homogeneous boundary conditions, allowing for velocity and magnetic field to be represented by orthogonal function expansions. One advantage to this is that the coefficients of the expansions form a set of nonlinearly interacting variables whose behavior can be described by equilibrium statistical mechanics, i.e., by a canonical ensemble theory based on the global invariants (energy, cross helicity and magnetic helicity) of ideal MHD turbulence. Another advantage is that truncated expansions provide a finite dynamical system whose time evolution can be numerically simulated to test the predictions of the associated statistical mechanics. If ensemble predictions are the same as time averages, then the system is said to be ergodic; if not, the system is nonergodic. Although it had been implicitly assumed in the early days of ideal MHD statistical theory development that these finite dynamical systems were ergodic, numerical simulations provided sufficient evidence that they were, in fact, nonergodic. Specifically, while canonical ensemble theory predicted that expansion coefficients would be (i) zero-mean random variables with (ii) energy that decreased with length scale, it was found that although (ii) was correct, (i) was not and the expected ergodicity was broken. The exact cause of this broken ergodicity was explained, after much investigation, by greatly extending the statistical theory of ideal MHD turbulence. The mathematical details of broken ergodicity, in fact, give a quantitative explanation of how coherent structure, dynamic alignment and force-free states appear in turbulent magnetofluids. The relevance of these ideal results to real MHD turbulence occurs because broken ergodicity is most manifest in the ideal case at the largest length scales and it is in these largest scales that a real magnetofluid has the least dissipation, i.e., most closely approaches the behavior of an ideal magnetofluid. Furthermore, the effects grow stronger when cross and magnetic helicities grow large with respect to energy, and this is exactly what occurs with time in a real magnetofluid, where it is called selective decay. The relevance of these results found in ideal MHD turbulence theory to the real world is that they provide at least a qualitative explanation of why confined turbulent magnetofluids, such as the liquid iron that fills the Earth's outer core, produce stationary, large-scale magnetic fields, i.e., the geomagnetic field. These results should also apply to other planets as well as to plasma confinement devices on Earth and in space, and the effects should be manifest if Reynolds numbers are high enough and there is enough time for stationarity to occur, at least approximately. In the presentation, details will be given for both theoretical and numerical results, and references will be provided.
Abstract This paper presents the development of a multi‐fidelity digital twin structural model (virtual model) of an as‐built wind turbine blade. The goal is to develop and demonstrate an approach to produce an accurate and detailed model of the as‐built blade for use in verifying the performance of the operating two‐bladed, downwind rotor. The digital twin model development methodology, presented herein, involves a novel calibration process to integrate a wide range of information including design specifications, manufacturing information, and structural testing data (modal and static) to produce a multi‐fidelity digital twin structural model: a detailed high‐fidelity model (i.e., 3D finite element analysis [FEA]) and consistent beam‐type models for aeroelastic simulation. A key element is that the multi‐fidelity structural digital twin method follows the rotor from the stages of design, to manufacturing, then to the ground testing and field operation. The result of this comprehensive approach is an accurate multi‐fidelity digital twin structural model for the geometric, structural, and structural dynamic properties of the as‐built blade within a 1% match in mass properties, 3.2% in blade frequencies, and 6% in deflection. The different stages of processing this information within the methodology are discussed. The rotor examined is the SUMR‐Demonstrator (SUMR‐D), which was installed on the Controls Advanced Research Testbed (CART‐2) wind turbine at the National Wind Technology Center. The digital twin model developed here was utilized to design controllers to safely operate SUMR‐D in field tests, which are providing additional data for further evaluation and development of the multi‐fidelity digital twin structural model.
The goal of flight flutter testing is to detect possibly destructive modes of aircraft vibration which may arise during flight from interaction of aerodynamic forces with structural dynamic properties of the airframe. This is typically accomplished by exciting the airframe with a time varying force and monitoring the response of the aircraft throughout its flight envelope. The data generated must be analyzed and presented so that the frequency and time of occurrence of excited modes are clearly and unambiguously displayed. Processing and display in near real time is also desirable. Display of data in the time-frequency plane is a natural choice because it is a familiar and intuitive framework. The Matching Pursuit algorithm provides a time-frequency analysis with good adaptability to signal structure and good signal representation in the time-frequency plane. Improvements in efficiency are needed before the algorithm can be used in real time, however.
Under the direction of the NASA In-Space Propulsion Technology Office, the team of L Garde, NASA Jet Propulsion Laboratory, Ball Aerospace, and NASA Langley Research Center has been developing a scalable solar sail configuration to address NASA s future space propulsion needs. Prior to a flight experiment of a full-scale solar sail, a comprehensive test program was implemented to advance the technology readiness level of the solar sail design. These tests consisted of solar sail component, subsystem, and sub-scale system ground tests that simulated the aspects of the space environment such as vacuum and thermal conditions. In July 2005, a 20-m four-quadrant solar sail system test article was tested in the NASA Glenn Research Center s Space Power Facility to measure its static and dynamic structural responses. Key to the maturation of solar sail technology is the development of validated finite element analysis (FEA) models that can be used for design and analysis of solar sails. A major objective of the program was to utilize the test data to validate the FEA models simulating the solar sail ground tests. The FEA software, ABAQUS, was used to perform the structural analyses to simulate the ground tests performed on the 20-m solar sail test article. This paper presents the details of the FEA modeling, the structural analyses simulating the ground tests, and a comparison of the pretest and post-test analysis predictions with the ground test results for the 20-m solar sail system test article. The structural responses that are compared in the paper include load-deflection curves and natural frequencies for the beam structural assembly and static shape, natural frequencies, and mode shapes for the solar sail membrane. The analysis predictions were in reasonable agreement with the test data. Factors that precluded better correlation of the analyses and the tests were unmeasured initial conditions in the test set-up.
A spline-based method for approximating thin shell dynamics is presented here. While the method is developed in the context of the Donnell-Mushtari thin shell equations, it can be easily extended to the Byrne-Flugge-Lur'ye equations or other models for shells of revolution as warranted by applications. The primary requirements for the method include accuracy, flexibility and efficiency in smart material applications. To accomplish this, the method was designed to be flexible with regard to boundary conditions, material nonhomogeneities due to sensors and actuators, and inputs from smart material actuators such as piezoceramic patches. The accuracy of the method was also of primary concern, both to guarantee full resolution of structural dynamics and to facilitate the development of PDE-based controllers which ultimately require real-time implementation. Several numerical examples provide initial evidence demonstrating the efficacy of the method.
The interactions of control system and distributed flexible structural dynamics is explored for mechanical arms. A modeling process using 4 x 4 transfer matrices is described which permits the closed loop response of many current arm configurations to be evaluated. Root locus, frequency response, modal shapes, and time impulse response have all been obtained from the digital computer implementation of this model, which is oriented to arm design and allows for easy variation of the arm configuration through data cards. The model corresponds with experimentally observed natural frequencies with an average error of less than 5% in the first three flexible modes in the seven cases considered. The model was used to explore the limits imposed by structural flexibility on a nondimensionalized two link arm with one and two joints for planar motion.
The paper is concerned with the task of estimating modal parameters from system response measurement in aircraft flutter testing. A frequency-domain derivation of an instrumental-variables algorithm is presented for a linear time-invariant dynamic system of order n. Basically, this algorithm fits a set of poles and zeros to the measured transfer function. An illustrative example is provided regarding the application of the algorithm to aeroelasticity testing. It is shown that the algorithm can be implemented for on-line data reduction with a microcomputer-based analysis system. By using instrumental variables the sensitivity of the modal parameter estimates to noise in the system-response measurements is reduced greatly. The algorithm is expected to be a powerful and valuable tool for on-line estimation of modal parameters in flutter testing and should be useful in control system and structural dynamics tests.