Search NASA⌕ Search

SEARCH · Search NASA

Results for “computational methods”

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 379 records · Page 21

An Automated Method to Compute Orbital Re-Entry Trajectories with Heating Constraints

Determining how to properly manipulate the controls of a re-entering re-usable launch vehicle (RLV) so that it is able to safely return to Earth and land involves the solution of a two-point boundary value problem (TPBVP). This problem, which can be quite difficult, is traditionally solved on the ground prior to flight. If necessary, a nearly unlimited amount of time is available to find the "best" solution using a variety of trajectory design and optimization tools. The role of entry guidance during flight is to follow the pre-determined reference solution while correcting for any errors encountered along the way. This guidance method is both highly reliable and very efficient in terms of onboard computer resources. There is a growing interest in a style of entry guidance that places the responsibility of solving the TPBVP in the actual entry guidance flight software. Here there is very limited computer time. The powerful, but finicky, mathematical tools used by trajectory designers on the ground cannot in general be made to do the job. Nonconvergence or slow convergence can result in disaster. The challenges of designing such an algorithm are numerous and difficult. Yet the payoff (in the form of decreased operational costs and increased safety) can be substantial. This paper presents an algorithm that incorporates features of both types of guidance strategies. It takes an initial RLV orbital re-entry state and finds a trajectory that will safely transport the vehicle to a Terminal Area Energy Management (TAEM) region. During actual flight, the computed trajectory is used as the reference to be flown by a more traditional guidance method.

Zimmerman, Curtis↗

A Survey of Methods for Computing Best Estimates of Endoatmospheric and Exoatmospheric Trajectories

Beginning with the mathematical prediction of planetary orbits in the early seventeenth century up through the most recent developments in sensor fusion methods, many techniques have emerged that can be employed on the problem of endo and exoatmospheric trajectory estimation. Although early methods were ad hoc, the twentieth century saw the emergence of many systematic approaches to estimation theory that produced a wealth of useful techniques. The broad genesis of estimation theory has resulted in an equally broad array of mathematical principles, methods and vocabulary. Among the fundamental ideas and methods that are briefly touched on are batch and sequential processing, smoothing, estimation, and prediction, sensor fusion, sensor fusion architectures, data association, Bayesian and non Bayesian filtering, the family of Kalman filters, models of the dynamics of the phases of a rocket's flight, and asynchronous, delayed, and asequent data. Along the way, a few trajectory estimation issues are addressed and much of the vocabulary is defined.

Bernard, William P.↗

Application of a finite element method for computing grazing incidence wave structure in an impedance tube - Comparison with experiment

The acoustic performance of a liner specimen, in a grazing incidence impedance tube, is analyzed using a finite element method. The liner specimen was designed to be a locally reacting, two-degree-of-freedom type with the resistance and reactance provided by perforated facesheets and compartmented cavities. Measured and calculated wave structures are compared for both normal and grazing incidence from 0.3 to 1.2 kHz. A finite element algorithm was incorporated into an optimization loop in order to predict liner grazing incidence impedance from measured SWR and null position data. Results suggest that extended reaction effects may have been responsible for differences between normal and grazing incidence impedance estimates.

Lester, H. C.↗

Computational Electromagnetic Modeling of SansEC(Trade Mark) Sensors

This paper describes the preliminary effort to apply computational design tools to aid in the development of an electromagnetic SansEC resonant sensor composite materials damage detection system. The computational methods and models employed on this research problem will evolve in complexity over time and will lead to the development of new computational methods and experimental sensor systems that demonstrate the capability to detect, diagnose, and monitor the damage of composite materials and structures on aerospace vehicles.

Smith, Laura J.↗

A Method for Computing the Core Flow in Three-Dimensional Leading-Edge Vortices

A theory is presented for calculating the flow in the core of a separation-induced leading-edge vortex. The method is based on matching inner and outer representations of the vortex. The inner model of the vortex is based on the quasicylindrical Navier-Stokes equations; the flow is assumed to be steady, axially symmetric, and incompressible and in addition, gradients in the radial direction are assumed to be much larger then gradients in the axial direction. The outer model is based on the three-dimensional free-vortex-sheet theory, a higher-order panel method which solves the Prandtl-Glauert equation including nonlinear boundary conditions pertinent to the concentrated vorticity representation of the leading edge vortex. The resultant flow is evaluated a posteriori for evidence of incipient vortex breakdown and the critical helix angle concept, in conjunction with an adverse longitudinal pressure gradient, is found to correlate well with the occurrence of vortex breakdown at the trailing edge of delta, arrow, and diamond wings.

Luckring, J. M.↗

Computational aeroelasticity challenges and resources

In the past decade, there has been much activity in the development of computational methods for the analysis of unsteady transonic aerodynamics about airfoils and wings. Significant features are illustrated which must be addressed in the treatment of computational transonic unsteady aerodynamics. The flow regimes for an aircraft on a plot of lift coefficient vs. Mach number are indicated. The sequence of events occurring in air combat maneuvers are illustrated. And further features of transonic flutter are illustrated. Also illustrated are several types of aeroelastic response which were encountered and which offer challenges for computational methods. The four cases illustrate problem areas encountered near the boundaries of aircraft envelopes, as operating condition change from high speed, low angle conditions to lower speed, higher angle conditions.

Edwards, John W.↗

Some methods of computing platform transmitter terminal location estimates

A position estimation algorithm was developed to track a humpback whale tagged with an ARGOS platform after a transmitter deployment failure and the whale's diving behavior precluded standard methods. The algorithm is especially useful where a transmitter location program exists; it determines the classical keplarian elements from the ARGOS spacecraft position vectors included with the probationary file messages. A minimum of three distinct messages are required. Once the spacecraft orbit is determined, the whale is located using standard least squares regression techniques. Experience suggests that in instances where circumstances inherent in the experiment yield message data unsuitable for the standard ARGOS reduction, (message data may be too sparse, span an insufficient period, or include variable-length messages). System ARGOS can still provide much valuable location information if the user is willing to accept the increased location uncertainties.

Hoisington, C. M.↗

Mapping methods for computationally efficient and accurate structural reliability

The influence of mesh coarseness in the structural reliability is evaluated. The objectives are to describe the alternatives and to demonstrate their effectiveness. The results show that special mapping methods can be developed by using: (1) deterministic structural responses from a fine (convergent) finite element mesh; (2) probabilistic distributions of structural responses from a coarse finite element mesh; (3) the relationship between the probabilistic structural responses from the coarse and fine finite element meshes; and (4) probabilistic mapping. The structural responses from different finite element meshes are highly correlated.

Shiao, Michael C.↗

New techniques for the analysis of manual control systems

Studies are summarized on the application of advanced analytical and computational methods to the development of mathematical models of human controllers in multiaxis manual control systems. Specific accomplishments include the following: (1) The development of analytical and computer methods for the measurement of random parameters in linear models of human operators. (2) Discrete models of human operator behavior in a multiple display situation were developed. (3) Sensitivity techniques were developed which make possible the identification of unknown sampling intervals in linear systems. (4) The adaptive behavior of human operators following particular classes of vehicle failures was studied and a model structure proposed.

Bekey, G. A.↗

Testing relativistic gravity theories using radio tracking data from planetary orbiting spacecraft.

We present a thorough analysis of a computational method for determining the numerical values of the relativity and other related dynamical parameters using two-way Doppler and ranging data from planetary orbiting spacecraft. The computational method consists of two parts. From Doppler data we first determine the earth-planet components of the position of the orbiting spacecraft relative to the center of gravity of the planet to high accuracy; adding the observed spacecraft range yields a range value to the center of the planet. These constructed earth-planet range data, referred to as normal points, are then treated as raw data in a regression analysis combined with planetary radar delay and meridian circle measurements to solve for the significant solar system dynamical parameters. The major errors sources in the planetary orbiter process are enumerated and their individual effects on the overall accuracy of the normal point accuracies are presented. The accuracies of the parameter estimates as a function of time, data sampling, and a priori assumptions are illustrated.

Jordan, J. F.↗

Computational predictive methods for fracture and fatigue

The damage-tolerant design philosophy as used by aircraft industries enables aircraft components and aircraft structures to operate safely with minor damage, small cracks, and flaws. Maintenance and inspection procedures insure that damages developed during service remain below design values. When damage is found, repairs or design modifications are implemented and flight is resumed. Design and redesign guidelines, such as military specifications MIL-A-83444, have successfully reduced the incidence of damage and cracks. However, fatigue cracks continue to appear in aircraft well before the design life has expired. The F16 airplane, for instance, developed small cracks in the engine mount, wing support, bulk heads, the fuselage upper skin, the fuel shelf joints, and along the upper wings. Some cracks were found after 600 hours of the 8000 hour design service life and design modifications were required. Tests on the F16 plane showed that the design loading conditions were close to the predicted loading conditions. Improvements to analytic methods for predicting fatigue crack growth adjacent to holes, when multiple damage sites are present, and in corrosive environments would result in more cost-effective designs, fewer repairs, and fewer redesigns. The overall objective of the research described in this paper is to develop, verify, and extend the computational efficiency of analysis procedures necessary for damage tolerant design. This paper describes an elastic/plastic fracture method and an associated fatigue analysis method for damage tolerant design. Both methods are unique in that material parameters such as fracture toughness, R-curve data, and fatigue constants are not required. The methods are implemented with a general-purpose finite element package. Several proof-of-concept examples are given. With further development, the methods could be extended for analysis of multi-site damage, creep-fatigue, and corrosion fatigue problems.

Cordes, J.↗

Parallelization of a Parabolized Navier-Stokes Solver with a Design Optimizer

The design of future supersonic aircraft, such as the High Speed Civil Transport (HSCT), will rely heavily on computational methods for aircraft design and the prediction of the complex aerodynamics encountered in flight. Parabolized Navier-Stokes (PNS) equation flow solvers are recognized as efficient and accurate computational tools for the solution of supersonic and hypersonic flow-fields, while design optimizers have the potential to be valuable tools within the overall design process. Presently, however, the execution of the flow solver in conjunction with a design optimizer presents a computationally intensive and formidable problem. To meet the challenges and increasing demand for multidisciplinary numerical tools which are faster, more robust and provide greater functionality, alternative strategies are explored to increase computational throughput by coupling a design optimizer and flow solver in a parallel processing environment. To address this problem the parallel processing of a PNS flow solver with a nonlinear constraint design optimizer is investigated as an alternative computational method.

Pallis, J. M.↗

Dynamic Shape Reconstruction of Three-Dimensional Frame Structures Using the Inverse Finite Element Method

A robust and efficient computational method for reconstructing the three-dimensional displacement field of truss, beam, and frame structures, using measured surface-strain data, is presented. Known as shape sensing , this inverse problem has important implications for real-time actuation and control of smart structures, and for monitoring of structural integrity. The present formulation, based on the inverse Finite Element Method (iFEM), uses a least-squares variational principle involving strain measures of Timoshenko theory for stretching, torsion, bending, and transverse shear. Two inverse-frame finite elements are derived using interdependent interpolations whose interior degrees-of-freedom are condensed out at the element level. In addition, relationships between the order of kinematic-element interpolations and the number of required strain gauges are established. As an example problem, a thin-walled, circular cross-section cantilevered beam subjected to harmonic excitations in the presence of structural damping is modeled using iFEM; where, to simulate strain-gauge values and to provide reference displacements, a high-fidelity MSC/NASTRAN shell finite element model is used. Examples of low and high-frequency dynamic motion are analyzed and the solution accuracy examined with respect to various levels of discretization and the number of strain gauges.

Gherlone, Marco↗

The analysis of Gemini 7 EEG data by three different computer oriented methods

The problem in assessment of electroencephalograms has been the subject element in visual evaluations of such records as well as the need for extensive training (2 or 3 years) before suitable candidate is capable to perform this task. There have been variety of computer type analyses of EEG during the past ten years or more, and we are still searching for a reliable, inexpensive, online data processing procedure for this task. We are aware our Soviet colleagues are routinely using computer assessment of EEGs, using relatively inexpensive hardware. In our report we endeavor to show how four different assessment procedures (visual assessment, zero crossing technique, smoothing and peak counting technique, and Weibull Statistic) may be utilized, and indicate the advantages and disadvantages peculiar to each. We wish to point out that in the Gemini 7 EEG analog records supplied to the Henry Ford Hospital, Channel II (mid-line central to mid-line occipital electrodes) could not be assessed reliably after approximately 26 hours of flight. Channel I (left central to left occipital electrodes) provided the more assessable EEG up to approximately 54 hours of flight. These channels were compared during the first 26 hours of flight, and no significant difference was found by our various assessment procedures. We, therefore, have limited this report to the findings on the analyses of Channel I during its approximate 54 hours of flight.

Proctor, L. D.↗

Modifications Of Discrete Ordinate Method For Computations With High Scattering Anisotropy: Comparative Analysis

A numerical accuracy analysis of the radiative transfer equation (RTE) solution based on separation of the diffuse light field into anisotropic and smooth parts is presented. The analysis uses three different algorithms based on the discrete ordinate method (DOM). Two methods, DOMAS and DOM2+, that do not use the truncation of the phase function, are compared against the TMS-method. DOMAS and DOM2+ use the Small-Angle Modification of RTE and the single scattering term, respectively, as an anisotropic part. The TMS method uses Delta-M method for truncation of the phase function along with the single scattering correction. For reference, a standard discrete ordinate method, DOM, is also included in analysis. The obtained results for cases with high scattering anisotropy show that at low number of streams (16, 32) only DOMAS provides an accurate solution in the aureole area. Outside of the aureole, the convergence and accuracy of DOMAS, and TMS is found to be approximately similar: DOMAS was found more accurate in cases with coarse aerosol and liquid water cloud models, except low optical depth, while the TMS showed better results in case of ice cloud.

Korkin, Sergey V.↗

Transonic Navier-Stokes calculations about a 65 deg delta wing

A computational study has been conducted in which the CFL3D Navier-Stokes solver coupled with an algebraic and a one-equation nonequilibrium turbulence model has been used to predict the flow over a 65 degree delta wing at transonic conditions for Reynolds numbers ranging from 6 x 10(exp 6) to 120 x 10(exp 6) based on mean aerodynamic chord. Solutions obtained indicated that the computational method when used with the one-equation turbulence model predicts results that compare well with experiment for attached flow conditions. Comparisons with experimental pressure at separated conditions show that the computational method, even though primary flow-field features are predicted well, does not predict secondary flow features.

Londenberg, W. Kelly↗

Computational simulation methods for composite fracture mechanics

Structural integrity, durability, and damage tolerance of advanced composites are assessed by studying damage initiation at various scales (micro, macro, and global) and accumulation and growth leading to global failure, quantitatively and qualitatively. In addition, various fracture toughness parameters associated with a typical damage and its growth must be determined. Computational structural analysis codes to aid the composite design engineer in performing these tasks were developed. CODSTRAN (COmposite Durability STRuctural ANalysis) is used to qualitatively and quantitatively assess the progressive damage occurring in composite structures due to mechanical and environmental loads. Next, methods are covered that are currently being developed and used at Lewis to predict interlaminar fracture toughness and related parameters of fiber composites given a prescribed damage. The general purpose finite element code MSC/NASTRAN was used to simulate the interlaminar fracture and the associated individual as well as mixed-mode strain energy release rates in fiber composites.

Murthy, Pappu L. N.↗