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.
Optimal temperature control design using IBM 7094 digital computer program for thermal environment simulation of stabilized multisurfaced earth- orbiting spacecraft structure
Predictive methods/computer codes for the computational simulation of acoustic fatigue resistance of hot composite structures subjected to acoustic excitation emanating from an adjacent vibrating component are discussed. Select codes developed over the past two decades at the NASA Lewis Research Center are used. The codes include computation of acoustic noise generated from a vibrating component, degradation in material properties of a composite laminate at use temperature, dynamic response of acoustically excited hot multilayered composite structure, degradation in the first ply strength of the excited structure due to acoustic loading, and acoustic fatigue resistance of the excited structure, including the propulsion environment. Effects of the laminate lay-up and environment on the acoustic fatigue life are evaluated. The results show that, by keeping the angled plies on the outer surface of the laminate, a substantial increase in the acoustic fatigue life is obtained. The effect of environment (temperature and moisture) is to relieve the residual stresses leading to an increase in the acoustic fatigue life of the excited panel.
We introduce improved numerical techniques for simulating the structural and compositional evolution of planetary rings due to micrometeoroid bombardment and subsequent ballistic transport of impact ejecta. Our current, robust code is capable of modeling structural changes and pollution transport simultaneously over long times on both local and global scales. In this paper, we describe the methodology based on the original structural code of Durisen et al. (1989, Icarus 80, 136-166) and on the pollution transport code of Cuzzi and Estrada (1998, Icarus 132, 1-35). We provide demonstrative simulations to compare with, and extend upon previous work, as well as examples of how ballistic transport can maintain the observed structure in Saturn's rings using available Cassini occultation optical depth data. In particular, we explicitly verify the claim that the inner B (and presumably A) ring edge can be maintained over long periods of time due to an ejecta distribution that is heavily biased in the prograde direction through a balance between the sharpening effects of ballistic transport and the broadening effects of viscosity. We also see that a "ramp"-like feature forms over time just inside that edge. However, it does not remain linear for the duration of the runs presented here unless a less steep ejecta velocity distribution is adopted. We also model the C ring plateaus and find that their outer edges can be maintained at their observed sharpness for long periods due to ballistic transport. We hypothesize that the addition of a significant component of a retrograde-biased ejecta distribution may help explain the linearity of the ramp and is probably essential for maintaining the sharpness of C ring plateau inner edges. This component would arise for the subset of micrometeoroid impacts which are destructive rather than merely cratering. Such a distribution will be introduced in future work.
Ionic liquids have been proposed as candidate electrolytes for high-energy density, rechargeable batteries. We present an extensive computational analysis supported by experimental comparisons of the bulk and interfacial properties of a representative set of these electrolytes as a function of Li-salt doping. We begin by investigating the bulk electrolyte using quantum chemistry and ab initio molecular dynamics to elucidate the solvation structure of Li(+). MD simulations using the polarizable force field of Borodin and coworkers were then performed, from which we obtain an array of thermodynamic and transport properties. Excellent agreement is found with experiments for diffusion, ionic conductivity, and viscosity. Combining MD simulations with electronic structure computations, we computed the electrochemical window of the electrolytes across a range of Li(+)-doping levels and comment on the role of the liquid environment. Finally, we performed a suite of simulations of these Li-doped electrolytes at ideal electrified interfaces to evaluate the differential capacitance and the equilibrium Li(+) distribution in the double layer. The magnitude of differential capacitance is in good agreement with our experiments and exhibits the characteristic camel-shaped profile. In addition, the simulations reveal Li(+) to be highly localized to the second molecular layer of the double layer, which is supported by additional computations that find this layer to be a free energy minimum with respect to Li(+) translation.
A FORTRAN-operated and an Ada-oriented design for the same system are compared to learn whether an essentially different design was produced using Ada. The designs were produced by an experiment that involves the parallel development of software for a spacecraft dynamics simulator. Design differences are identified in the use of abstractions, system structure, and simulator operations. Although the designs were significantly different, this result may be influenced by some special characteristics discussed.
A FORTRAN oriented and an Ada oriented design for the same system are compared to learn whether an essentially different design was produced using Ada. The designs were produced by an experiment that involves the parallel development of software for a spacecraft dynamics simulator. Design differences are identified in the use of abstractions, system structure, and simulator operations. Although the designs were vastly different, this result may be influenced by some special characteristics discussed.
VULCAN-CFD offers a comprehensive set of capabilities to enable the simulation of continuum flowfields from subsonic to hypersonic conditions. The governing equations that are employed include allowances for both chemical and thermal nonequilibrium processes, coupled with a wide variety of turbulence models for both Reynolds-averaged and large eddy simulations. The software package can simulate two-dimensional, axisymmetric, or three-dimensional problems on structured multiblock meshes or fully unstructured meshes. A parabolic (i.e., space-marching) treatment can also be used for any subset of a structured mesh that can accommodate this solution strategy. The flow solver provides a significant level of geometric flexibility for structured grid simulations by allowing for arbitrary face-to-face C(0) continuous and non-C(0) continuous block interface connectivities. The unstructured grid paradigm allows for mixed element unstructured meshes that contain any combination of tetrahedral, prismatic, pyramidal, and hexahedral cell elements. The flow solver is also fully parallelized using MPI (Message Passing Interface) libraries in a data-parallel fashion, allowing for efficient simulations on modern High Performance Computing (HPC) systems. This document provides information related to the installation and execution of the VULCAN-CFD software package. A detailed description of the physical and numerical models available in the software are provided in the VULCAN-CFD Theory Manual.
A computational simulation procedure is described for multidisciplinary analysis and tailoring of layered multi-material hot composite engine structural components subjected to simultaneous multiple discipline-specific thermal, structural, vibration, and acoustic loads. The effect of aggressive environments is also simulated. The simulation is based on a three-dimensional finite element analysis technique in conjunction with structural mechanics codes, thermal/acoustic analysis methods, and tailoring procedures. The integrated multidisciplinary simulation procedure is general-purpose including the coupled effects of nonlinearities in structure geometry, material, loading, and environmental complexities. The composite material behavior is assessed at all composite scales, i.e., laminate/ply/constituents (fiber/matrix), via a nonlinear material characterization hygro-thermo-mechanical model. Sample tailoring cases exhibiting nonlinear material/loading/environmental behavior of aircraft engine fan blades, are presented. The various multidisciplinary loads lead to different tailored designs, even those competing with each other, as in the case of minimum material cost versus minimum structure weight and in the case of minimum vibration frequency versus minimum acoustic noise.
Results from numerical simulations of auroral current sheet equilibria and associated V-shaped potential structures are presented. It is shown that with allowance for both hot magnetospheric ion and cold ionospheric ion populations, the perpendicular potential drop, Phi(m), associated with a non-neutral auroral current sheet is critically controlled by the temperature of the 'heated' ionospheric ions. The heating is caused by the wave turbulence excited by the auroral current sheet. In the presence of heated ionospheric ions, a relatively large variation in the temperature of the hot magnetospheric ion population causes a very small variation in the potential drop Phi(m). The perpendicular potential drop acts to produce a V-shaped double layer with multiple potential steps parallel to the magnetic field when a zero potential boundary condition is imposed at the ionospheric boundary. Outside the V-shaped potential structure, ionospheric return currents develop self-consistently.
A relatively simple approach to the analysis of linear spacecraft control/structure interaction problems is presented. The approach uses a commercially available structural system dynamic analysis package for both controller and plant dynamics, thus obviating the need to transfer data between separate programs. The unilateral coupling between components in the control system block diagram is simulated using sparse matrix stiffness and damping elements available in the structural dynamic code. The approach is illustrated with a series of simple tutorial examples of a rigid spacecraft core with flexible appendages.
We have simulated plasma transport processes in the presence of a quasi-two-dimensional current filament, that generated kV potential structure in the auroral region. The simulation consists of a set of one-dimensional flux tube simulations with different imposed time-dependent, field-aligned currents. The model uses the 16 moment system of equations and simultaneously solves coupled continuity and momentum equations and equations describing the transport along the magnetic field lines of parallel and perpendicular thermal energy and heat flows for each species. The lower end of the simulation is at an altitude of 800 km, in the collisional topside ionosphere, while the upper end is at 10 R(sub E) in the magnetosphere. The plasma consists of hot electrons and protons of magnetospheric origin and low-energy electrons, protons, and oxygen ions of ionospheric origin. The dynamical interaction of the individual current filaments with ionospheric and magnetospheric plasma generates a potential structure in the horizontal direction and kilovolt field-aligned potential drops along the field lines. The side-by-side display exhibits the evolution of the implied potential structure in the horizontial direction. In the presence of this potential structure and parallel electric field ionospheric plasma density is depleted and velocity is reduced, while density enhancement and increased velocity is observed in magnetospheric plasma. The ionospheric and magnetospheric electron temperatures increase below 2 R(sub E) due to magnetic mirror force on converging geomagnetic field lines. The primary cross-field motion produced by the horizontal E field (E x B drift) is perpendicular to both of the significant spatial directions and is thus ignorable in this geometry. The effects of other cross-field drift processes are discussed. The simulation thus provides insight into the dynamical evolution of two-dimensional potential structures driven by an imposed finite width, field-aligned current profile.
The flow past the Boeing Transonic Truss-Braced Wing (TTBW) aircraft is simulated with a Hybrid Reynolds-Averaged Navier Stokes Large-Eddy Simulations (HRLES) turbulence modeling approach to study transonic buffet onset and evaluate the predictive capability of the numerical approach. All simulations used structured overset curvilinear grids for the free-air configuration, with all simulations performed using the Launch, Ascent, and Vehicle Aerodynamics (LAVA) computational framework. Aerodynamic loads and surface pressure obtained from the HRLES solutions are compared to results from unsteady RANS (URANS) and experiments conducted at NASA Ames 11-by 11-Foot Transonic Wind Tunnel. Overall good agreement is obtained in the predicted loads and surface pressure. Unsteady pressure data from HRLES and URANS simulations is utilized to compute power spectral density (PSD) to predict transonic buffet onset and compared to the Kulite data obtained from the experiment. Comparisons of the PSD spectrum also show reasonable agreement with the wind tunnel data, with HRLES and URANS predicting buffet onset at a slightly earlier angle-of-attack than the experiment.
Software aspects of control-structure interaction (CSI) analysis are discussed. The following subject areas are covered: (1) implementation of a partitioned algorithm for simulation of large CSI problems; (2) second-order discrete Kalman filtering equations for CSI simulations; and (3) parallel computations and control of adaptive structures.
Draper Station Analysis Tool (DSAT) is a computer program, built on commercially available software, for simulating and analyzing complex dynamic systems. Heretofore used in designing and verifying guidance, navigation, and control systems of the International Space Station, DSAT has a modular architecture that lends itself to modification for application to spacecraft or terrestrial systems. DSAT consists of user-interface, data-structures, simulation-generation, analysis, plotting, documentation, and help components. DSAT automates the construction of simulations and the process of analysis. DSAT provides a graphical user interface (GUI), plus a Web-enabled interface, similar to the GUI, that enables a remotely located user to gain access to the full capabilities of DSAT via the Internet and Webbrowser software. Data structures are used to define the GUI, the Web-enabled interface, simulations, and analyses. Three data structures define the type of analysis to be performed: closed-loop simulation, frequency response, and/or stability margins. DSAT can be executed on almost any workstation, desktop, or laptop computer. DSAT provides better than an order of magnitude improvement in cost, schedule, and risk assessment for simulation based design and verification of complex dynamic systems.
Analysis of interaction occurring between space debris and orbiting structures is of great interest to the planning and survivability of space assets. Computer simulation of the impact events using hydrodynamic codes can provide some understanding of the processes but the problems involved with this fundamental approach are formidable. First, any realistic simulation is necessarily three-dimensional, e.g., the impact and breakup of a satellite. Second, the thickness of important components such as satellite skins or bumper shields are small with respect to the dimension of the structure as a whole, presenting severe zoning problems for codes. Thirdly, the debris cloud produced by the primary impact will yield many secondary impacts which will contribute to the damage and possible breakup of the structure. The problem was approached by choosing a relatively new computational technique that has virtues peculiar to space impacts. The method is called Smoothed Particle Hydrodynamics.
Airframe noise is an important issue during the approach and landing phases of flight because engines are throttled back, and airframe sources become an equal or greater contributor to environmental noise. The leading-edge slat component of the high-lift system for transport aircraft has been identified as a prominent source of airframe noise. Noise production mechanisms and potential treatments to mitigate airframe noise have been areas of active research for nearly three decades. Two promising technologies for slat noise reduction are the slat-gap filler (SGF) and slat-cove filler (SCF). Previous work was done using a 2D section of a transport class aircraft wing to demonstrate their feasibility via benchtop and computational models. A wind tunnel test campaign was recently undertaken to test the SGF and SCF technologies in a flight-like environment by installing dynamically scaled treatments on the 10%-scale High-Lift Common Research Model (CRM-HL), which is representative of a commercial twin-aisle transport aircraft. Advancement of the SGF and SCF technologies to overcome challenges of integration with 3D airframe structures, computational simulations of the CRM-HL predicting the aerodynamic and aeroacoustic performance and structural loads, and design of wind tunnel model hardware with features accommodating noise treatments were conducted in parallel. Details of innovations in the SGF and SCF technologies for 3D airframe feasibility will be described. The loads the treatments must sustain in flight and the influence of those loads on the design is also discussed. Finally, the design of wing leading-edge assemblies to treat the 10% CRM-HL with dynamically scaled SGF and SCF components is presented. The treatments were successfully tested in 2020/2021 providing a demonstration of leading-edge noise reduction using components that would allow slat articulation.
Graphite/epoxy composite thin shell structures were simulated to investigate damage and fracture progression due to internal pressure and axial loading. Defective and defect-free structures (thin cylinders) were examined. The three different laminates examined had fiber orientations of (90/0/+/-0)(sub s), where 0 is 45, 60, and 75 deg. CODSTRAN, an integrated computer code that scales up constituent level properties to the structural level and accounts for all possible failure modes, was used to simulate composite degradation under loading. Damage initiation, growth, accumulation, and propagation to fracture were included in the simulation. Burst pressures for defective and defect-free shells were compared to evaluate damage tolerance. The results showed that damage initiation began with matrix failure whereas damage and/or fracture progression occurred as a result of additional matrix failure and fiber fracture. In both thin cylinder cases examined (defective and defect-free), the optimum layup configuration was (90/0/+/-60)(sub s) because it had the best damage tolerance with respect to the burst pressure.
The NESSUS probabilistic structural analysis computer program combines state-of-the-art probabilistic algorithms with general purpose structural analysis methods to compute the probabilistic response and the reliability of engineering structures. Uncertainty in loading, material properties, geometry, boundary conditions and initial conditions can be simulated. The structural analysis methods include nonlinear finite element and boundary element methods. Several probabilistic algorithms are available such as the advanced mean value method and the adaptive importance sampling method. The scope of the code has recently been expanded to include probabilistic life and fatigue prediction of structures in terms of component and system reliability and risk analysis of structures considering cost of failure. The code is currently being extended to structural reliability considering progressive crack propagation. Several examples are presented to demonstrate the new capabilities.