Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “computation time”

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 541 records · Page 30

Implementation of a Hybrid Edge Node-Centroid Node Approach for the Generation of Reduced Thermal Models

Reduced thermal models are often required for delivery to organizations that manage observatory or launch models at the highest levels of assembly. However, the effort to generate reduced models, and verify against their detailed counterparts, is a challenge that has not yet been conclusively solved. Higher level organizations often place a limit on the number of nodes for delivered models with the assumption that smaller models generally result in less computation time. However, the burden of producing and verifying the accuracy of the reduced models is placed primarily on the lower-level organizations, which in turn consumes resources needed to produce these models. Limiting the total number of allowable nodes may also prevent users from taking full advantage of software capabilities that allow for faster generation of models, such as finite elements, which generally require more nodes than older centroid based models. A methodology using Thermal Desktop was described in 2010 which used: (1) finite elements and edge nodes for a conduction matrix, (2) centroid nodes for capacitance and radiative computations, and (3) the super network feature to produce a conduction matrix based only on the centroid nodes. At that time, the methodology was clear, but the implementation would have had to be done manually; however, with the inclusion of the OpenTD API, this methodology can now be implemented programmatically and for the first time, be a viable approach for the generation of reduced models. The approach was implemented and developed at the NASA Goddard Space Flight Center (GSFC) for the Capture, Containment, and Return System (CCRS) payload on the Earth Return Orbiter (ERO) as part of the Mars Sample Return (MSR) mission, resulting in the TCYEE tool. ERO features a spacecraft bus provided by Airbus through the European Space Agency (ESA) with node limitations on the delivered CCRS model provided by GSFC. TCYEE was used to generate the reduced model for this delivery and the predictions compared favorably to the detailed model currently in use for the thermal performance evaluation. Furthermore, TCYEE is being explored for potential use on the Roman Space Telescope (RST) for the generation of reduced models for delivery to the launch provider, which also has node limit requirements on the RST observatory model for use in launch simulation analyses. This paper describes the methodology, its implementation, and compares the performance of reduced models generated by TCYEE to their detailed counterparts.

Hume L. Peabody↗

Implementation of a Hybrid Edge Node-Centroid Node Approach for the Generation of Reduced Thermal Models

Reduced thermal models are often required for delivery to organizations that manage observatory or launch models at the highest levels of assembly. However, the effort to generate reduced models, and verify against their detailed counterparts, is a challenge that has not yet been conclusively solved. Higher level organizations often place a limit on the number of nodes for delivered models with the assumption that smaller models generally result in less computation time. However, the burden of producing and verifying the accuracy of the reduced models is placed primarily on the lower-level organizations, which in turn consumes resources needed to produce these models. Limiting the total number of allowable nodes may also prevent users from taking full advantage of software capabilities that allow for faster generation of models, such as finite elements, which generally require more nodes than older centroid based models. A methodology using Thermal Desktop was described in 2010 which used: (1) finite elements and edge nodes for a conduction matrix, (2) centroid nodes for capacitance and radiative computations, and (3) the super network feature to produce a conduction matrix based only on the centroid nodes. At that time, the methodology was clear, but the implementation would have had to be done manually; however, with the inclusion of the OpenTD API, this methodology can now be implemented programmatically and for the first time, be a viable approach for the generation of reduced models. The approach was implemented and developed at the NASA Goddard Space Flight Center (GSFC) for the Capture, Containment, and Return System (CCRS) payload on the Earth Return Orbiter (ERO) as part of the Mars Sample Return (MSR) mission, resulting in the TCYEE tool. ERO features a spacecraft bus provided by Airbus through the European Space Agency (ESA) with node limitations on the delivered CCRS model provided by GSFC. TCYEE was used to generate the reduced model for this delivery and the predictions compared favorably to the detailed model currently in use for the thermal performance evaluation. Furthermore, TCYEE is being explored for potential use on the Roman Space Telescope (RST) for the generation of reduced models for delivery to the launch provider, which also has node limit requirements on the RST observatory model for use in launch simulation analyses. This paper describes the methodology, its implementation, and compares the performance of reduced models generated by TCYEE to their detailed counterparts.

Thermal desktop↗

Simple numerical method for predicting steady compressible flows

A numerical method for solving the isenthalpic form of the governing equations for compressible viscous and inviscid flows was developed. The method was based on the concept of flux vector splitting in its implicit form. The method was tested on several demanding inviscid and viscous configurations. Two different forms of the implicit operator were investigated. The time marching to steady state was accelerated by the implementation of the multigrid procedure. Its various forms very effectively increased the rate of convergence of the present scheme. High quality steady state results were obtained in most of the test cases; these required only short computational times due to the relative efficiency of the basic method.

Vonlavante, Ernst↗

Upwind algorithm for the parabolized Navier-Stokes equations

A new upwind algorithm based on Roe's scheme has been developed to solve the two-dimensional parabolized Navier-Stokes equations. This method does not require the addition of user-specified smoothing terms for the capture of discontinuities such as shock waves. Thus, the method is easy to use and can be applied without modification to a wide variety of supersonic flowfields. The advantages and disadvantages of this adaptation are discussed in relation to those of the conventional Beam-Warming (1978) scheme in terms of accuracy, stability, computer time and storage requirements, and programming effort. The new algorithm has been validated by applying it to three laminar test cases, including flat-plate boundary-layer flow, hypersonic flow past a 15-deg compression corner, and hypersonic flow into a converging inlet. The computed results compare well with experiment and show a dramatic improvement in the resolution of flowfield details when compared with results obtained using the conventional Beam-Warming algorithm.

Lawrence, Scott L.↗

Variance Decomposition of MEDLI2 Reconstructed Heating Using Neural Networks

The Mars Entry, Descent, and Landing Instrumentation (MEDLI2) sensor suite collected data during entry of the Mars 2020 Perseverance rover into Mars’ atmosphere. An inverse estimation of the backshell and heatshield surface aeroheating was performed, using the data from the MEDLI2 Instrumented Sensor Plugs, a network of thermocouples embedded within the thermal protection system across the aeroshell. Monte Carlo analysis was conducted to assess the sensitivity of the surface heat rate, temperature, and heat load to uncertainties in thermocouple depth and material properties. In this paper, a variance decomposition method using Sobol indices was employed to understand the relative contributions of each uncertainty parameter. Performing this analysis using results from the inverse analysis tool FIAT_Opt was found to require incredibly high computation time, and thus machine learning models were trained and evaluated as a surrogate model for FIAT_Opt. This paper demonstrates that machine learning models can be an efficient, accurate alternative to state-of-the-art inverse analysis tools like FIAT_Opt, especially for computationally-expensive processes. Using these models, the sensitivity analysis showed that uncertainties in heat capacity and thermal conductivity were the main drivers for the overall uncertainty in peak reconstructed heating and heat load.

H S Alpert↗

A Simplified Model of VIPER Thermal Management System. Part I: Loop Heat Pipe

NASA’s Volatiles Investigating Polar Exploration Rover (VIPER) uses a combination of loop heat pipes (LHPs) and heaters as the primary hardware to modulate the temperature of the electronic boxes and payloads instruments. LHPs have been designed and instrumented to be capable of passive and active shutdown. The goal of the passive shutdown, using a thermal control valve (TCV) to limit the flow of the working fluid, is to preserve thermal energy and avoid unnecessary heat leaks during lunar night. The goal of the active shutdown on the other hand, is to forcefully increase the vehicle internal temperature to bakeout the payload instruments. Several numerical modeling tools are available to analyze LHP performance under the complex set of operational requirements and environments levied on the VIPER Thermal Management System (TMS). Yet, while those numerical modeling tools are capable of detailed performance analysis, they tend to take large computational resources and long computational time. In this paper, a model of a single LHP developed based on control volume approach is presented. By using the control volume approach, the complex geometry of the LHP, and the components attached to it, are reduced to thermal lump capacitances, thus reducing the computation effort. This paper also describes the effort to correlate the proposed LHP model using thermal vacuum (TVAC) test data.

Loop Heat Pipe↗

A Simplified Model of VIPER Thermal Management System. Part I: Loop Heat Pipe

NASA’s Volatiles Investigating Polar Exploration Rover (VIPER) uses a combination of loop heat pipes (LHPs) and heaters as the primary hardware to modulate the temperature of the electronic boxes and payloads instruments. LHPs have been designed and instrumented to be capable of passive and active shutdown. The goal of the passive shutdown, using a thermal control valve (TCV) to limit the flow of the working fluid, is to preserve thermal energy and avoid unnecessary heat leaks during lunar night. The goal of the active shutdown on the other hand, is to forcefully increase the vehicle internal temperature to bakeout the payload instruments. Several numerical modeling tools are available to analyze LHP performance under the complex set of operational requirements and environments levied on the VIPER Thermal Management System (TMS). Yet, while those numerical modeling tools are capable of detailed performance analysis, they tend to take large computational resources and long computational time. In this paper, a model of a single LHP developed based on control volume approach is presented. By using the control volume approach, the complex geometry of the LHP, and the components attached to it, are reduced to thermal lump capacitances, thus reducing the computation effort. This paper also describes the effort to correlate the proposed LHP model using thermal vacuum (TVAC) test data.

Loop Heat Pipe↗

Comparison of Fixed and Variable Time Step Trajectory Integration Methods for Cislunar Trajectories

Due to the nonlinear nature of the Earth-Moon-Sun three-body problem and non-spherical gravity, CEV cislunar targeting algorithms will require many propagations in their search for a desired trajectory. For on-board targeting especially, the algorithm must have a simple, fast, and accurate propagator to calculate a trajectory with reasonable computation time, and still be robust enough to remain stable in the various flight regimes that the CEV will experience. This paper compares Cowell s method with a fourth-order Runge- Kutta integrator (RK4), Encke s method with a fourth-order Runge-Kutta- Nystr m integrator (RKN4), and a method known as Multi-Conic. Additionally, the study includes the Bond-Gottlieb 14-element method (BG14) and extends the investigation of Encke-Nystrom methods to integrators of higher order and with variable step size.

Weeks, ichael W.↗

A Special Purpose Architecture for Finite Element Analysis

The analysis of aerospace structures by the finite element method consumes considerable computer time. The cost of this resource and the designer's desire to have rapid feedback concerning such questions as the effect of a change in loading of the structure or in a parameter of some structural material led to the design of a special purpose parallel computing system for finite element analysis. As a special purpose computer, the architecture of this finite element computer is closely tied to computational aspects of the particular problem. Various aspects of an MIMD array of microprocessors are related to the requirements of the class of finite element analysis problems which it is intended to solve.

Jordan, H. F.↗

Time domain flutter analysis of cascades using a full-potential solver

A time domain approach is used to determine the dynamic aeroelastic stability of a cascade of blades. The structural model for each blade is a typical section with two degrees of freedom. The aerodynamic model is the unsteady, two-dimensional, full-potential flow through the cascade of airfoils. The unsteady equations of motion for the structure and the fluid are integrated simultaneously in time starting with the steady flowfield and a small initial disturbance applied to the airfoils. The motion of each blade is analyzed to determine the aeroelastic stability of the cascade. The effect of interblade phase angle is included in the analysis by allowing each blade to have an independent motion and considering a number of blade passages. Calculations are made using an airfoil section and structural parameters that are representative of a propfan. The results are compared with those from a separate frequency domain analysis. Good agreement between the results is observed. With the time domain approach, it is possible to consider nonlinear structural models and nonlinear force-displacement relations. The method allows a realistic simulation of the motion of the fluid and the cascade blades for a better physical understanding and it also has the potential for saving computational time when compared to the frequency domain approach for the flutter analysis of cascades.

Bakhle, Milind A.↗

Comparison of several methods for obtaining the time response of linear systems to either a unit impulse or arbitrary input from frequency-response data

Several methods of obtaining the time response of Linear systems to either a unit impulse or an arbitrary input from frequency-response data are described and compared. Comparisons indicate that all the methods give good accuracy when applied to a second-order system; the main difference is the required computing time. The methods generally classified as inverse Laplace transform methods were found to be most effective in determining the response to a unit impulse from frequency-response data of higher order systems. Some discussion and examples are given of the use of such methods as flight-data-analysis techniques in predicting loads and motions of a flexible aircraft on the basis of simple calculations when the aircraft frequency response is known.

Donegan, James J↗

A thin-shock-layer solution for nonequilibrium, inviscid hypersonic flows in earth, Martian, and Venusian atmospheres

An approximate inverse solution is presented for the nonequilibrium flow in the inviscid shock layer about a vehicle in hypersonic flight. The method is based upon a thin-shock-layer approximation and has the advantage of being applicable to both subsonic and supersonic regions of the shock layer. The relative simplicity of the method makes it ideally suited for programming on a digital computer with a significant reduction in storage capacity and computing time required by other more exact methods. Comparison of nonequilibrium solutions for an air mixture obtained by the present method is made with solutions obtained by two other methods. Additional cases are presented for entry of spherical nose cones into representative Venusian and Martian atmospheres. A digital computer program written in FORTRAN language is presented that permits an arbitrary gas mixture to be employed in the solution. The effects of vibration, dissociation, recombination, electronic excitation, and ionization are included in the program.

Grose, W. L.↗

A Monte Carlo investigation of thrust imbalance of solid rocket motor pairs

A technique is described for theoretical, statistical evaluation of the thrust imbalance of pairs of solid-propellant rocket motors (SRMs) firing in parallel. Sets of the significant variables, determined as a part of the research, are selected using a random sampling technique and the imbalance calculated for a large number of motor pairs. The performance model is upgraded to include the effects of statistical variations in the ovality and alignment of the motor case and mandrel. Effects of cross-correlations of variables are minimized by selecting for the most part completely independent input variables, over forty in number. The imbalance is evaluated in terms of six time - varying parameters as well as eleven single valued ones which themselves are subject to statistical analysis. A sample study of the thrust imbalance of 50 pairs of 146 in. dia. SRMs of the type to be used on the space shuttle is presented. The FORTRAN IV computer program of the analysis and complete instructions for its use are included. Performance computation time for one pair of SRMs is approximately 35 seconds on the IBM 370/155 using the FORTRAN H compiler.

Sforzini, R. H.↗

Biaxial thermo-mechanical fatigue

Stress-strain and durability information is often desirable for situations in which strain and temperature are changing simultaneously. To obtain such information, strain controlled uniaxial push-pull tests have typically been done. In order to control the mechanical strain, it is necessary in such tests to compute the mechanical strain from the total measured strain using measured temperature and the thermal expansion properties of the specimen. A system for conducting torsional thermomechanical tests is described which has the great advantage that the torsional strain is unaffected by the changing temperature and thus real time computations of quantities is not required for control of the test and the mechanical strain need not be determined from the subtraction of two measured qnantities as is the case in the uniaxial test. In addition to describing torsional thermomechanical tests, guidelines for software to be used in running biaxial thermomechanical tests will also be presented.

Jordan, Eric H.↗

Image processing enhancement of high-resolution TEM micrographs of nanometer-size metal particles

The high-resolution TEM detectability of lattice fringes from metal particles supported on substrates is impeded by the substrate itself. Single value decomposition (SVD) and Fourier filtering (FFT) methods were applied to standard high resolution micrographs to enhance lattice resolution from particles as well as from crystalline substrates. SVD produced good results for one direction of fringes, and it can be implemented as a real-time process. Fourier methods are independent of azimuthal directions and allow separation of particle lattice planes from those pertaining to the substrate, which makes it feasible to detect possible substrate distortions produced by the supported particle. This method, on the other hand, is more elaborate, requires more computer time than SVD and is, therefore, less likely to be used in real-time image processing applications.

Artal, P.↗

A solution to the Navier-Stokes equations based upon the Newton Kantorovich method

An implicit finite difference scheme based on the Newton-Kantorovich technique was developed for the numerical solution of the nonsteady, incompressible, two-dimensional Navier-Stokes equations in conservation-law form. The algorithm was second-order-time accurate, noniterative with regard to the nonlinear terms in the vorticity transport equation except at the earliest few time steps, and spatially factored. Numerical results were obtained with the technique for a circular cylinder at Reynolds number 15. Results indicate that the technique is in excellent agreement with other numerical techniques for all geometries and Reynolds numbers investigated, and indicates a potential for significant reduction in computation time over current iterative techniques.

Davis, J. E.↗

A Model for Predicting Thermoelectric Properties of Bi2Te3

A parameterized orthogonal tight-binding mathematical model of the quantum electronic structure of the bismuth telluride molecule has been devised for use in conjunction with a semiclassical transport model in predicting the thermoelectric properties of doped bismuth telluride. This model is expected to be useful in designing and analyzing Bi2Te3 thermoelectric devices, including ones that contain such nano - structures as quantum wells and wires. In addition, the understanding gained in the use of this model can be expected to lead to the development of better models that could be useful for developing other thermoelectric materials and devices having enhanced thermoelectric properties. Bi2Te3 is one of the best bulk thermoelectric materials and is widely used in commercial thermoelectric devices. Most prior theoretical studies of the thermoelectric properties of Bi2Te3 have involved either continuum models or ab-initio models. Continuum models are computationally very efficient, but do not account for atomic-level effects. Ab-initio models are atomistic by definition, but do not scale well in that computation times increase excessively with increasing numbers of atoms. The present tight-binding model bridges the gap between the well-scalable but non-atomistic continuum models and the atomistic but poorly scalable ab-initio models: The present tight-binding model is atomistic, yet also computationally efficient because of the reduced (relative to an ab-initio model) number of basis orbitals and flexible parameterization of the Hamiltonian.

Lee, Seungwon↗