Computation of normal modes from identified complex modes
(Previously announced in STAR as N82-75431)
SEARCH · Search NASA
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.
(Previously announced in STAR as N82-75431)
Postprocessing algorithms devised to facilitate vibrational-mode analyses of dynamics of complicated structures. Yields inertia matrices and elastic/rigid-coupling matrices. Such analyses important in simulation and control in active suppression of vibrations in large building or in precise aiming of large antenna.
An existing computer program, used for predicting the natural frequencies and mode shapes of helicopter rotor blades, was refined to improve program accuracy and versatility. The program is based on the Holzer-Myklestad approach adapted for rotating beams. Coupled vertical (out-of-plane), horizontal (in-plane), and torsional mode characteristics were determined for a variety of hub and blade configurations. The resulting program is documented by presenting the recursion equations and techniques for determining natural frequencies and mode shapes, input data requirements, and descriptions of various program outputs. The accuracy of the program is demonstrated by comparing computed results with exact solutions to classical problems and experimental data.
A simulation of time-optimal intercept algorithms for on-board computation of control commands is described. The effects of three different display modes and two different computation modes on the pilots' ability to intercept a moving target in minimum time were tested. Both computation modes employed singular perturbation theory to help simplify the two-point boundary value problem associated with trajectory optimization. Target intercept time was affected by both the display and computation modes chosen, but the display mode chosen was the only significant influence on the miss distance.
This User's Manual was prepared to provide the engineer with the information required to run the coupled mode version of the Normal Modes Rotor Aeroelastic Analysis Computer Program. The manual provides a full set of instructions for running the program, including calculation of blade modes, calculations of variable induced velocity distribution and the calculation of the time history of the response for either a single blade or a complete rotor with an airframe (the latter with constant inflow).
High-frequency pressure-fluctuation measurements were made in AEDC Tunnel 9 at Mach 10 and the NASA Langley 15-Inch Mach 6 and 31-Inch Mach 10 tunnels. Measurements were made on a 7deg-half-angle cone model. Pitot measurements of freestream pressure fluctuations were also made in Tunnel 9 and the Langley Mach-6 tunnel. For the first time, second-mode waves were measured in all of these tunnels, using 1-MHz-response pressure sensors. In Tunnel 9, second-mode waves could be seen in power spectra computed from records as short as 80 micro-s. The second-mode wave amplitudes were observed to saturate and then begin to decrease in the Langley tunnels, indicating wave breakdown. Breakdown was estimated to occur near N approx. equals 5 in the Langley Mach-10 tunnel. The unit-Reynolds-number variations in the data from Tunnel 9 were too large to see the same processes. In Tunnel 9, the measured transition locations were found to be at N = 4.5 using thermocouples, and N = 5.3 using 50-kHz-response pressure sensors. What appears to be a very long transitional region was observed at a unit Reynolds number of 13.5 million per meter in Tunnel 9. These results were consistent with the high-frequency pressure fluctuation measurements. High-frequency pressure fluctuation measurements indicated that transition did occur in the Langley Mach-6 tunnel, but the location of transition was not precisely determined. Unit Reynolds numbers in the Langley Mach-10 tunnel were too low to observe transition. More analysis of this data set is expected in the future.
Using the normal modes of the GLAS 4th Order Model, a Machenhauer nonlinear normal mode initialization (NLNMI) was carried out for the external vertical mode using the GLAS 4th Order shallow water equations model for an equivalent depth corresponding to that associated with the external vertical mode. A simple procedure was devised which was directed at identifying computational modes by following the rate of increase of BAL sub M, the partial (with respect to the zonal wavenumber m) sum of squares of the time change of the normal mode coefficients (for fixed vertical mode index) varying over the latitude index L of symmetric or antisymmetric gravity waves. A working algorithm is presented which speeds up the convergence of the iterative Machenhauer NLNMI. A 24 h integration using the NLNMI state was carried out using both Matsuno and leap-frog time-integration schemes; these runs were then compared to a 24 h integration starting from a non-initialized state. The maximal impact of the nonlinear normal mode initialization was found to occur 6-10 hours after the initial time.
Progress is reported in obtaining valuable results needed for the full interpretation of the spectral reflectance curves of solar system objects. The degree to which photometric geometry affects spectral reflectance curves was demonstrated. Various forms of photometric functions were compared and a function adequate for describing the scattering properties of low and moderately reflecting materials was developed and applied in a study of the phase coefficients of various materials, as well as in a study of how the shape of a body affects the spectral reflectance properties. The adequacy of the photometric function for Mars-like analogs was studied. The goniometer system is being converted to a computer driven mode. As soon as computer controls are integrated in the goniometer, the phase dependence 0.95 micron feature in meteorite spectra is scheduled to begin.
Two computational procedures for analyzing complex structural systems for their natural modes and frequencies of vibration are presented. Both procedures are based on a substructures methodology and both employ the finite-element stiffness method to model the constituent substructures. The first procedure is a direct method based on solving the eigenvalue problem associated with a finite-element representation of the complete structure. The second procedure is a component-mode synthesis scheme in which the vibration modes of the complete structure are synthesized from modes of substructures into which the structure is divided. The analytical basis of the methods contains a combination of features which enhance the generality of the procedures. The computational procedures exhibit a unique utilitarian character with respect to the versatility, computational convenience, and ease of computer implementation. The computational procedures were implemented in two special-purpose computer programs. The results of the application of these programs to several structural configurations are shown and comparisons are made with experiment.
Computational Fluid Dynamics (CFD) analysis results are compared with benchmark quality test data from the Propulsion Engineering Research Center's (PERC) Rocket Based Combined Cycle (RBCC) experiments to verify fluid dynamic code and application procedures. RBCC engine flowpath development will rely on CFD applications to capture the multi-dimensional fluid dynamic interactions and to quantify their effect on the RBCC system performance. Therefore, the accuracy of these CFD codes must be determined through detailed comparisons with test data. The PERC experiments build upon the well-known 1968 rocket-ejector experiments of Odegaard and Stroup by employing advanced optical and laser based diagnostics to evaluate mixing and secondary combustion. The Finite Difference Navier Stokes (FDNS) code was used to model the fluid dynamics of the PERC RBCC ejector mode configuration. Analyses were performed for both Diffusion and Afterburning (DAB) and Simultaneous Mixing and Combustion (SMC) test conditions. Results from both the 2D and the 3D models are presented.
Computer program for determining transverse natural frequencies and mode shapes of uniform beams
NASA Glenn Research Center's Computational Materials Group has integrated a graphical user interface with in-house-developed surface modeling capabilities, with the goal of using computationally efficient atomistic simulations to aid the development of advanced aerospace materials, through the modeling of alloy surfaces, surface alloys, and segregation. The software is also ideal for modeling nanomaterials, since surface and interfacial effects can dominate material behavior and properties at this level. Through the combination of an accurate atomistic surface modeling methodology and an efficient computational engine, it is now possible to directly model these types of surface phenomenon and metallic nanostructures without a supercomputer. Fulfilling a High Operating Temperature Propulsion Components (HOTPC) project level-I milestone, a graphical user interface was created for a suite of quantum approximate atomistic materials modeling Fortran programs developed at Glenn. The resulting "Alloy Design Workbench-Surface Modeling Package" (ADW-SMP) is the combination of proven quantum approximate Bozzolo-Ferrante-Smith (BFS) algorithms (refs. 1 and 2) with a productivity-enhancing graphical front end. Written in the portable, platform independent Java programming language, the graphical user interface calls on extensively tested Fortran programs running in the background for the detailed computational tasks. Designed to run on desktop computers, the package has been deployed on PC, Mac, and SGI computer systems. The graphical user interface integrates two modes of computational materials exploration. One mode uses Monte Carlo simulations to determine lowest energy equilibrium configurations. The second approach is an interactive "what if" comparison of atomic configuration energies, designed to provide real-time insight into the underlying drivers of alloying processes.
Design and performance details of the advanced information processing system (AIPS) for fault and damage tolerant data processing on aircraft and spacecraft are presented. AIPS comprises several computers distributed throughout the vehicle and linked by a damage tolerant data bus. Most I/O functions are available to all the computers, which run in a TDMA mode. Each computer performs separate specific tasks in normal operation and assumes other tasks in degraded modes. Redundant software assures that all fault monitoring, logging and reporting are automated, together with control functions. Redundant duplex links and damage-spread limitation provide the fault tolerance. Details of an advanced design of a laboratory-scale proof-of-concept system are described, including functional operations.
This thesis describes the investigation of superconducting multi-mode cavity architecture for superconducting transmon-based quantum computing. The dissertation highlights useful features of radio-frequency cavities used for quantum computing, with a brief discussion on the advantages of superconducting cavities. In the subsequent section, the concept of transmon is introduced and the mechanism of coupling with a cavity is analyzed. Once the foundations of the topic are laid down, the thesis focuses on optimizing a multi-mode superconducting rf cavity design originally developed for high-energy physics applications. Such section articulates in two main parts, the first part concerning the "bare" cavity remake via finite-elements eigenmode simulations using a computer-aided design software called CST Studio Suite®. In the second part, a transmon qubit is physically inserted into the modified cavity to assess the qubit-cavity coupling of the new design. In evaluating said coupling, two distinct analysis methods are used, namely the black-box quantization method and the energy participation ratio method, both implemented using Ansys® High-Frequency Electromagnetic-Field Simulator, or HFSS™. Results from the two evaluations, compared together, show that the optimized design meets the requisites to be used for quantum computing purposes.
The numerical modeling of the normal modes of the global oceans is addressed. The results of such modeling could be expected to serve as a guide in the analysis of observations and measurements intended to detect these modes. The numerical computation of normal modes of the global oceans is a field in which several investigations have obtained results during the past 15 years. The results seem to be model-dependent to an unsatisfactory extent. Some modeling areas, such as higher resolution of the bathymetry, inclusion of self-attraction and loading, the role of the Arctic Ocean, and systematic testing by means of diagnostic models are addressed. The results show that the present state of the art is such that a final solution to the normal mode problem still lies in the future. The numerical experiments show where some of the difficulties are and give some insight as to how to proceed in the future.
Truncation errors in natural frequencies as computed by method of component mode synthesis
GPU Monte Carlo Simulation Code with a taste of RASPA We present enhancements in Monte Carlo simulation speed and functionality within an open-source code, gRASPA, which uses graphical processing units (GPUs) to achieve significant performance improvements compared to serial, CPU implementations of Monte Carlo. The code supports a wide range of Monte Carlo simulations, including canonical ensemble (NVT), grand canonical, NVT Gibbs, Widom test particle insertions, and continuous-fractional component Monte Carlo. Implementation of grand canonical transition matrix Monte Carlo (GC-TMMC) and a novel feature to allow different moves for the different components of metal-organic framework (MOF) structures exemplify the capabilities of gRASPA for precise free energy calculations and enhanced adsorption studies, respectively. The introduction of a High-Throughput Computing (HTC) mode permits many Monte Carlo simulations on a single GPU device for accelerated materials discovery. The code can incorporate machine learning (ML) potentials. The open-source nature of gRASPA promotes reproducibility and openness in science, and users may add features to the code and optimize it for their own purposes. The code is written in CUDA/C++ and SYCL/C++ to support different GPU vendors. The gRASPA code is publicly available at https://github.com/snurr-group/gRASPA.
A perturbation technique is used to compute the rotational corrections to the nonradial oscillation spectrum of a realistic neutron star model. The l = l(0) oscillations are coupled to l = l(0) + or - 1 oscillations by the Coriolis force. To calculate the rotational corrections to the quadrupole, l(0) = 2, oscillation modes, previous calculations of the nonradial oscillation spectrum of this neutron star model are extended to include l(0) = 1 and l(0) = 3 modes, as well as many higher overtones at each value of l(o). As an example of this technique, the rotational corrections to two quadrupole toroidal modes are computed. For these modes the coupling to the spheroidal class of oscillations produced by the Coriolis force introduces a nonzero radial component to the velocity field, as well as a nonvanishing Lagrangian density perturbation. This result is used to compute the neutrino damping rates for these toroidal modes of a rotating neutron star. The neutrino damping time can approach the gravitational radiation damping time in rotating neutron stars if the central temperature is high enough. The rotationally induced coupling of spheroidal oscillations to toroidal modes can also produce significant displacements at the stellar surface in some of the toroidal modes. The damping rates due to this process for the rotationally corrected toroidal modes are estimated. The coupling of surface motions to the toroidal modes also have interesting implications for channeling energy, e.g., that associated with a glitch in the crust, to the surface of the star. Perhaps this might produce observable effects in the pulsar emission process or a gamma-ray burst event.