Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “time integrator”

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 955 records · Page 53
All sourcesNASA NTRSDOE OSTI
NASA NTRS1995

Robust integration schemes for generalized viscoplasticity with internal-state variables. Part 1: Theoretical developments and applications

This two-part report is concerned with the development of a general framework for the implicit time-stepping integrators for the flow and evolution equations in generalized viscoplastic models. The primary goal is to present a complete theoretical formulation, and to address in detail the algorithmic and numerical analysis aspects involved in its finite element implementation, as well as to critically assess the numerical performance of the developed schemes in a comprehensive set of test cases. On the theoretical side, the general framework is developed on the basis of the unconditionally-stable, backward-Euler difference scheme as a starting point. Its mathematical structure is of sufficient generality to allow a unified treatment of different classes of viscoplastic models with internal variables. In particular, two specific models of this type, which are representative of the present start-of-art in metal viscoplasticity, are considered in applications reported here; i.e., fully associative (GVIPS) and non-associative (NAV) models. The matrix forms developed for both these models are directly applicable for both initially isotropic and anisotropic materials, in general (three-dimensional) situations as well as subspace applications (i.e., plane stress/strain, axisymmetric, generalized plane stress in shells). On the computational side, issues related to efficiency and robustness are emphasized in developing the (local) interative algorithm. In particular, closed-form expressions for residual vectors and (consistent) material tangent stiffness arrays are given explicitly for both GVIPS and NAV models, with their maximum sizes 'optimized' to depend only on the number of independent stress components (but independent of the number of viscoplastic internal state parameters). Significant robustness of the local iterative solution is provided by complementing the basic Newton-Raphson scheme with a line-search strategy for convergence. In the present first part of the report, we focus on the theoretical developments, and discussions of the results of numerical-performance studies using the integration schemes for GVIPS and NAV models.

Saleeb, Atef F.↗
NASA NTRS1995

Robust integration schemes for generalized viscoplasticity with internal-state variables. Part 2: Algorithmic developments and implementation

This two-part report is concerned with the development of a general framework for the implicit time-stepping integrators for the flow and evolution equations in generalized viscoplastic models. The primary goal is to present a complete theoretical formulation, and to address in detail the algorithmic and numerical analysis aspects involved in its finite element implementation, as well as to critically assess the numerical performance of the developed schemes in a comprehensive set of test cases. On the theoretical side, the general framework is developed on the basis of the unconditionally-stable, backward-Euler difference scheme as a starting point. Its mathematical structure is of sufficient generality to allow a unified treatment of different classes of viscoplastic models with internal variables. In particular, two specific models of this type, which are representative of the present start-of-art in metal viscoplasticity, are considered in applications reported here; i.e., fully associative (GVIPS) and non-associative (NAV) models. The matrix forms developed for both these models are directly applicable for both initially isotropic and anisotropic materials, in general (three-dimensional) situations as well as subspace applications (i.e., plane stress/strain, axisymmetric, generalized plane stress in shells). On the computational side, issues related to efficiency and robustness are emphasized in developing the (local) interative algorithm. In particular, closed-form expressions for residual vectors and (consistent) material tangent stiffness arrays are given explicitly for both GVIPS and NAV models, with their maximum sizes 'optimized' to depend only on the number of independent stress components (but independent of the number of viscoplastic internal state parameters). Significant robustness of the local iterative solution is provided by complementing the basic Newton-Raphson scheme with a line-search strategy for convergence. In the present second part of the report, we focus on the specific details of the numerical schemes, and associated computer algorithms, for the finite-element implementation of GVIPS and NAV models.

Li, Wei↗
NASA NTRS2005

Rotordynamic Analysis and Feasibility Study of a Disk Spin Test Facility for Rotor Health Monitoring

Recently, National Aeronautics and Space Administration (NASA) initiated a program to achieve the significant improvement in aviation safety. One of the technical challenges is the design and development of accelerated experiments that mimic critical damage cases encountered in engine components. The Nondestructive Evaluation (NDE) Group at the NASA Glenn Research Center (GRC) is currently addressing the goal concerning propulsion health management and the development of propulsion system specific technologies intended to detect potential failures prior to catastrophe. For this goal the unique disk spin simulation system was assembled at NASA GRC, which allows testing of rotors with the spinning speeds up to 10K RPM, and at the elevated temperature environment reaching 540 C (1000 F). It is anticipated that the facility can be employed for detection of Low Cycle Fatigue disk cracking and further High Cycle Fatigue blade vibration. The controlled crack growth studies at room and elevated temperatures can be conducted on the turbine wheels, and various NDE techniques can be integrated and assessed as in-situ damage monitoring tools. Critical rotating parts in advanced gas turbine engines such as turbine disks frequently operate at high temperature and stress for long periods of time. The integrity of these parts must be proven by non-destructive evaluation (NDE) during various machining steps ranging from forging blank to finished shape, and also during the systematic overhaul inspections. Conventional NDE methods, however, have unacceptable limits. Some of these techniques are time-consuming and inconvenient for service aircraft testing. Almost all of these techniques require that the vicinity of the damage is known in advance. These experimental techniques can provide only local information and no indication of the structural strength at a component and/or system level. The shortcomings of currently available NDE methods lead to the requirement of new damage detection techniques that can provide global information on the rotating components/system, and, in addition, they do not require direct human access to the operating system. During this period of research considerable effort was directed towards the further development of experimental facility and development of the vibration-based crack detection methodology for rotating disks and shafts. A collection of papers and reports were written to describe the results of this work. The attached captures that effort and represents the research output during the grant period.

Sawicki, Jerzy T.↗
NASA NTRS2005

Afterbody Heating Predictions for a Mars Science Laboratory Entry Vehicle

The Mars Science Laboratory mission intends to deliver a large rover to the Martian surface within 10 km of its target site. One candidate entry vehicle aeroshell consists of a 3.75-m diameter, 70-deg sphere-cone forebody and a biconic afterbody similar to that of Viking. This paper presents computational fluid dynamics predictions of laminar afterbody heating rates for this configuration and a 2010 arrival at Mars. Computational solutions at flight conditions used an 8-species Mars gas model in chemical and thermal non-equilibrium. A grid resolution study examined the effects of mesh spacing on afterbody heating rates and resulted in grids used for heating predictions on a reference entry trajectory. Afterbody heating rate reaches its maximum value near 0.6 W/sq cm on the first windward afterbody cone at the time of peak freestream dynamic pressure. Predicted afterbody heating rates generally are below 3% of the forebody laminar nose cap heating rate throughout the design trajectory. The heating rates integrated over time provide total heat load during entry, which drives thermal protection material thickness.

Edquist, Karl T.↗
NASA NTRS2005

Flight Test Comparison Between Enhanced Vision (FLIR) and Synthetic Vision Systems

Limited visibility and reduced situational awareness have been cited as predominant causal factors for both Controlled Flight Into Terrain (CFIT) and runway incursion accidents. NASA s Synthetic Vision Systems (SVS) project is developing practical application technologies with the goal of eliminating low visibility conditions as a causal factor to civil aircraft accidents while replicating the operational benefits of clear day flight operations, regardless of the actual outside visibility condition. A major thrust of the SVS project involves the development/demonstration of affordable, certifiable display configurations that provide intuitive out-the-window terrain and obstacle information with advanced pathway guidance. A flight test evaluation was conducted in the summer of 2004 by NASA Langley Research Center under NASA s Aviation Safety and Security, Synthetic Vision System - Commercial and Business program. A Gulfstream G-V aircraft, modified and operated under NASA contract by the Gulfstream Aerospace Corporation, was flown over a 3-week period at the Reno/Tahoe International Airport and an additional 3-week period at the NASA Wallops Flight Facility to evaluate integrated Synthetic Vision System concepts. Flight testing was conducted to evaluate the performance, usability, and acceptance of an integrated synthetic vision concept which included advanced Synthetic Vision display concepts for a transport aircraft flight deck, a Runway Incursion Prevention System, an Enhanced Vision Systems (EVS), and real-time Database Integrity Monitoring Equipment. This paper focuses on comparing qualitative and subjective results between EVS and SVS display concepts.

Arthur, Jarvis J., III↗
NASA NTRS2015

SOAREX-8 Suborbital Experiments 2015 - A New Paradigm for Small Spacecraft Communication

In 2015 NASA plans to launch a payload to 280 Km altitude on a sounding rocket from the Wallops Flight Facility. This payload will contain several novel technologies that work together to demonstrate methodologies for space sample return missions and for nanosatellite communications in general. The payload will deploy and test an Exo-Brake, which slows the payload aerodynamically, providing eventual de-orbit and recovery of future ISS samples through a Small Payload Quick Return project. In addition, this flight addresses future Mars mission entry technology, space-to-space communications using the Iridium Short Messaging Service (SMS), GPS tracking, and wireless sensors using the ZigBee protocol. SOAREX-8 is being assembled and tested at Ames Research Center (ARC) and the NASA Engineering and Safety Center (NESC) is funding sensor and communications work. Open source Arduino technology and software are used for system control. The ZigBee modules used are XBee units that connect analog sensors for temperature, air pressure and acceleration measurement wirelessly to the payload telemetry system. Our team is developing methods for power distribution and module mounting, along with software for sensor integration, data assembly and downlink. We have demonstrated relaying telemetry to the ground using the Iridium satellite constellation on a previous flight, but the upcoming flight will be the first time we integrate useful flight test data from a ZigBee wireless sensor network. Wireless sensor data will measure the aerodynamic efficacy of the Exo-Brake permitting further on orbit flight tests of improved designs. The Exo-Brake is 5 sq m in area and will be stored in a container and deployed during ascent once the payload is jettisoned from the launch vehicle. We intend to further refine the hardware and continue testing on balloon launches, future sounding rocket flights and on nanosatellite missions. The use of standards-based and open source hardware/software has allowed for this project to be completed with a very modest budget and a challenging schedule. There is a wealth of hardware and software available for both the Arduino platform and the XBee, all low-cost or open-source. Along with the Exo-Brake hardware and deployment discussion, this paper will describe in detail the system architecture emphasizing the successful use of open source hardware and software to minimize effort and cost. Testing procedures, radio frequency interference (RFI) mitigation, success criteria and expected results will also be discussed. The use of Iridium short messaging capability for space-to-space links, standards-based wireless sensor networks, and other innovative communications technology are also presented.

aerobrake↗
NASA NTRS2018

Practical Aspects of Real-Time Modeling for the Learn-To-Fly Concept

Practical aspects of autonomous onboard real-time modeling for the NASA Learn-to-Fly concept are examined using flight data from two subscale test aircraft. A practical real-time global nonlinear aerodynamic modeling method is developed and explained, along with the multiple-input excitation needed for effective real-time modeling. Critical issues for integrating real-time global nonlinear aerodynamic modeling and local linear modeling with the real-time learning adaptive control and guidance elements of the Learn-to-Fly algorithm are identified and discussed. Real-time modeling results from NASA Learn-to-Fly flight demonstrations are presented and evaluated using model fit metrics and prediction tests.

Morelli, Eugene A.↗
NASA NTRS1979

Evaluation of MOSTAS computer code for predicting dynamic loads in two bladed wind turbines

Calculated dynamic blade loads were compared with measured loads over a range of yaw stiffnesses of the DOE/NASA Mod-O wind turbine to evaluate the performance of two versions of the MOSTAS computer code. The first version uses a time-averaged coefficient approximation in conjunction with a multi-blade coordinate transformation for two bladed rotors to solve the equations of motion by standard eigenanalysis. The second version accounts for periodic coefficients while solving the equations by a time history integration. A hypothetical three-degree of freedom dynamic model was investigated. The exact equations of motion of this model were solved using the Floquet-Lipunov method. The equations with time-averaged coefficients were solved by standard eigenanalysis.

Kaza, K. R. V.↗
NASA NTRS1995

Coronal mass ejections and the evolution of the large-scale corona

Coronal mass ejections (CMEs) are dynamic, large-scale events in the solar corona that expel plasma and magnetic fields into the solar wind. The structure of the large-scale corona influences CME occurrence; for example, CMEs are most frequently observed to be disruptions of coronal or helmet streamer configurations on the Sun. We have investigated the evolution of the large-scale corona in the presence of differential rotation, using time-dependent magnetohydrodynamic (MHD) simulations in two and three dimensions. An equilibrium coronal configuration is first computed by superimposing a Parker solar wind flow on an initial potential field, the equations are integrated in time until a steady-state is reached. When differential rotation is introduced, we find that the configuration does not reach a steady state; the closed field regions (helmet streamers) disrupt recurrently, ejecting plasmoids into the solar wind. Our results suggest that differential rotation may be one of the mechanisms by which mass ejections are initiated. We will describe the evolution of the polarization brightness (pB) for the resulting configurations, as well as the propagation of the ejected plasmoids in the solar wind, and we will compare the effects of differential rotation with other possible initiation mechanisms, such as the emergence of new magnetic flux.

Linker, J. A.↗
NASA NTRS2000

Supercomputing Aspects for Simulating Incompressible Flow

The primary objective of this research is to support the design of liquid rocket systems for the Advanced Space Transportation System. Since the space launch systems in the near future are likely to rely on liquid rocket engines, increasing the efficiency and reliability of the engine components is an important task. One of the major problems in the liquid rocket engine is to understand fluid dynamics of fuel and oxidizer flows from the fuel tank to plume. Understanding the flow through the entire turbo-pump geometry through numerical simulation will be of significant value toward design. One of the milestones of this effort is to develop, apply and demonstrate the capability and accuracy of 3D CFD methods as efficient design analysis tools on high performance computer platforms. The development of the Message Passage Interface (MPI) and Multi Level Parallel (MLP) versions of the INS3D code is currently underway. The serial version of INS3D code is a multidimensional incompressible Navier-Stokes solver based on overset grid technology, INS3D-MPI is based on the explicit massage-passing interface across processors and is primarily suited for distributed memory systems. INS3D-MLP is based on multi-level parallel method and is suitable for distributed-shared memory systems. For the entire turbo-pump simulations, moving boundary capability and efficient time-accurate integration methods are built in the flow solver, To handle the geometric complexity and moving boundary problems, an overset grid scheme is incorporated with the solver so that new connectivity data will be obtained at each time step. The Chimera overlapped grid scheme allows subdomains move relative to each other, and provides a great flexibility when the boundary movement creates large displacements. Two numerical procedures, one based on artificial compressibility method and the other pressure projection method, are outlined for obtaining time-accurate solutions of the incompressible Navier-Stokes equations. The performance of the two methods is compared by obtaining unsteady solutions for the evolution of twin vortices behind a flat plate. Calculated results are compared with experimental and other numerical results. For an unsteady flow, which requires small physical time step, the pressure projection method was found to be computationally efficient since it does not require any subiteration procedure. It was observed that the artificial compressibility method requires a fast convergence scheme at each physical time step in order to satisfy the incompressibility condition. This was obtained by using a GMRES-ILU(0) solver in present computations. When a line-relaxation scheme was used, the time accuracy was degraded and time-accurate computations became very expensive.

Kwak, Dochan↗
NASA NTRS2020

System Design for Near-Global Imaging of Triton

We propose a near-term low-cost mission to perform imaging of over 80% of the surface of Neptune’s moon Triton during a single flyby. To accomplish this we take advantage of Triton’s 141 hr orbit period and design narrow-angle and wide-angle cameras to image over the wide range of distances and lighting con-ditions. Of particular interest in this paper is the process for high-resolution im-aging as the range to the surface comes toward closest approach. To perform this, we traded a time-delay integration camera design with a slew-and-settle full-frame imager, and decided on the latter implementation. This required a high-performance attitude control system to cover the full disk in the allotted time, so a number of mosaic designs were iterated with simulated ADCS capa-bilities until a system design was established that showed appropriate margins. This paper describes the system design, including mission design summary, camera design summary, mosaicking pattern options, timelines, and detailed ADCS performance simulations showing slew and settle and pointing stability.

Frazier, William↗
NASA NTRS1985

Optimal generalized multistep integration formulae for real-time digital simulation

The problem of discretizing a dynamical system for real-time digital simulation is considered. Treating the system and its simulation as stochastic processes leads to a statistical characterization of simulator fidelity. A plant discretization procedure based on an efficient matrix generalization of explicit linear multistep discrete integration formulae is introduced, which minimizes a weighted sum of the mean squared steady-state and transient error between the system and simulator outputs.

Moerder, D. D.↗
NASA NTRS1988

A numerical method for computing unsteady 2-D boundary layer flows

A numerical method for computing unsteady two-dimensional boundary layers in incompressible laminar and turbulent flows is described and applied to a single airfoil changing its incidence angle in time. The solution procedure adopts a first order panel method with a simple wake model to solve for the inviscid part of the flow, and an implicit finite difference method for the viscous part of the flow. Both procedures integrate in time in a step-by-step fashion, in the course of which each step involves the solution of the elliptic Laplace equation and the solution of the parabolic boundary layer equations. The Reynolds shear stress term of the boundary layer equations is modeled by an algebraic eddy viscosity closure. The location of transition is predicted by an empirical data correlation originating from Michel. Since transition and turbulence modeling are key factors in the prediction of viscous flows, their accuracy will be of dominant influence to the overall results.

Krainer, Andreas↗
NASA NTRS1984

Addendum: Development of a preprototype times wastewater recovery subsystem

The results of the second generation operational improvements and the TIMES (Thermoelectric Integrated Membrane Evaporation Subsystem) 2 study are covered. Areas covered in the second generation operational improvements are improved temperature control, water quality improvements, subsytem operational improvements, solid handling improvements, wastewater pretreatment optimization, and membrane rejuvenation concepts. The task for the TIMES 2 study are thermoelectric regenerator improvement, recycle loop pH operational criteria, recycle loop component optimization, and hollow fiber membrane evaporator improvement. Results are presented and conclusions are drawn from both studies.

Dehner, G. F.↗
NASA NTRS1986

Long-Term Electronic Timer

Timing circuit turns on power source on command, then turns it off again after preset interval. In comparison with prior devices, unit consumes little power and is smaller, lighter in weight, and less complicated. Timer includes oscillator and counter in integrated circuit. Timing interval equals oscillator period multiplied by number of cycles to be counted.

Temple, G.↗
NASA NTRS1993

Conical Euler analysis and active roll suppression for unsteady vortical flows about rolling delta wings

A conical Euler code was developed to study unsteady vortex-dominated flows about rolling, highly swept delta wings undergoing either forced motions or free-to-roll motions that include active roll suppression. The flow solver of the code involves a multistage, Runge-Kutta time-stepping scheme that uses a cell-centered, finite-volume, spatial discretization of the Euler equations on an unstructured grid of triangles. The code allows for the additional analysis of the free to-roll case by simultaneously integrating in time the rigid-body equation of motion with the governing flow equations. Results are presented for a delta wing with a 75 deg swept, sharp leading edge at a free-stream Mach number of 1.2 and at 10 deg, 20 deg, and 30 deg angle of attack alpha. At the lower angles of attack (10 and 20 deg), forced-harmonic analyses indicate that the rolling-moment coefficients provide a positive damping, which is verified by free-to-roll calculations. In contrast, at the higher angle of attack (30 deg), a forced-harmonic analysis indicates that the rolling-moment coefficient provides negative damping at the small roll amplitudes. A free-to-roll calculation for this case produces an initially divergent response, but as the amplitude of motion grows with time, the response transitions to a wing-rock type of limit cycle oscillation, which is characteristic of highly swept delta wings. This limit cycle oscillation may be actively suppressed through the use of a rate-feedback control law and antisymmetrically deflected leading-edge flaps. Descriptions of the conical Euler flow solver and the free-to roll analysis are included in this report. Results are presented that demonstrate how the systematic analysis of the forced response of the delta wing can be used to predict the stable, neutrally stable, and unstable free response of the delta wing. These results also give insight into the flow physics associated with unsteady vortical flows about delta wings undergoing forced motions and free-to-roll motions, including the active suppression of the wing-rock type phenomenon. The conical Euler methodology developed is directly extend able to three-dimensional calculations.

Lee-Rausch, Elizabeth M.↗
NASA NTRS1981

Fuel-Conservation Guidance System for Powered-Lift Aircraft

A technique is described for the design of fuel-conservative guidance systems and is applied to a system that was flight tested on board NASA's sugmentor wing jet STOL research aircraft. An important operational feature of the system is its ability to rapidly synthesize fuel-efficient trajectories for a large set of initial aircraft positions, altitudes, and headings. This feature allows the aircraft to be flown efficiently under conditions of changing winds and air traffic control vectors. Rapid synthesis of fuel-efficient trajectories is accomplished in the airborne computer by fast-time trajectory integration using a simplified dynamic performance model of the aircraft. This technique also ensures optimum flap deployment and, for powered-lift STOL aircraft, optimum transition to low-speed flight. Also included in the design is accurate prediction of touchdown time for use in four-dimensional guidance applications. Flight test results have demonstrated that the automatically synthesized trajectories produce significant fuel savings relative to manually flown conventional approaches.

Erzberger, Heinz↗
NASA NTRSDate not supplied

Scale-Resolving Simulations of a Supersonic Retro-Propulsion Concept For Mars Entry, Descent, and Landing

Supersonic retro-propulsion (SRP) is a deceleration technology that could enable largerpayloads and potentially humans to be brought from space to the surface of Mars safely.Accurate and reliable CFD predictions are needed to help design future Mars entry, descent,and landing (EDL) vehicles due to the cost and limitations of wind tunnel tests. TypicalCFD approaches struggle at points along the Mars EDL trajectory with high enough altitude(low enough ambient pressure) that the SRP motor plumes are over-expanded. This workdemonstrates how high-order scale-resolving simulations on Cartesian AMR grids can provideaccurate and reliable time-averaged aerodynamic loads on the vehicle for this challengingflow regime. Greater than second order convergence in the time-averaged integrated axialforce𝐶𝐴is obtained across the three finest mesh resolutions simulated, and the estimateduncertainty is smaller than the unsteady standard deviation – demonstrating the accuracy andreliability of the method. Although costly, this approach can yield higher confidence in𝐶𝐴thanpreviously obtained in a reasonable turnaround time and is within reach for a few points in thetrajectory (i.e. Mach 2.4, 2.52, and 2.78). Further computational performance and numericalmethods improvements are needed to reduce the cost enough to bring this method into the EDLengineering design cycle.

S&M↗