Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “analytical 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 469 records · Page 26

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.↗

A new method to evaluate human-robot system performance

One of the key issues in space exploration is that of deciding what space tasks are best done with humans, with robots, or a suitable combination of each. In general, human and robot skills are complementary. Humans provide as yet unmatched capabilities to perceive, think, and act when faced with anomalies and unforeseen events, but there can be huge potential risks to human safety in getting these benefits. Robots provide complementary skills in being able to work in extremely risky environments, but their ability to perceive, think, and act by themselves is currently not error-free, although these capabilities are continually improving with the emergence of new technologies. Substantial past experience validates these generally qualitative notions. However, there is a need for more rigorously systematic evaluation of human and robot roles, in order to optimize the design and performance of human-robot system architectures using well-defined performance evaluation metrics. This article summarizes a new analytical method to conduct such quantitative evaluations. While the article focuses on evaluating human-robot systems, the method is generally applicable to a much broader class of systems whose performance needs to be evaluated.

Space Flight/instrumentation↗

Advanced Methods for Aircraft Engine Thrust and Noise Benefits: Nozzle-Inlet Flow Analysis

Three connected sub-projects were conducted under reported project. Partially, these sub-projects are directed to solving the problems conducted by the HU/FM&AL under two other NASA grants. The fundamental idea uniting these projects is to use untraditional 3D corrugated nozzle designs and additional methods for exhaust jet noise reduction without essential thrust lost and even with thrust augmentation. Such additional approaches are: (1) to add some solid, fluid, or gas mass at discrete locations to the main supersonic gas stream to minimize the negative influence of strong shock waves forming in propulsion systems; this mass addition may be accompanied by heat addition to the main stream as a result of the fuel combustion or by cooling of this stream as a result of the liquid mass evaporation and boiling; (2) to use porous or permeable nozzles and additional shells at the nozzle exit for preliminary cooling of exhaust hot jet and pressure compensation for non-design conditions (so-called continuous ejector with small mass flow rate; and (3) to propose and analyze new effective methods fuel injection into flow stream in air-breathing engines. Note that all these problems were formulated based on detailed descriptions of the main experimental facts observed at NASA Glenn Research Center. Basically, the HU/FM&AL Team has been involved in joint research with the purpose of finding theoretical explanations for experimental facts and the creation of the accurate numerical simulation technique and prediction theory for solutions for current problems in propulsion systems solved by NASA and Navy agencies. The research is focused on a wide regime of problems in the propulsion field as well as in experimental testing and theoretical and numerical simulation analysis for advanced aircraft and rocket engines. The F&AL Team uses analytical methods, numerical simulations, and possible experimental tests at the Hampton University campus. We will present some management activity and theoretical numerical simulation results obtained by the FM&AL Team in the reporting period in accordance with the schedule of the work.

Morgan, Morris H.↗

Aerospace Toxicology and Microbiology

Toxicology dates to the very earliest history of humanity with various poisons and venom being recognized as a method of hunting or waging war with the earliest documentation in the Evers papyrus (circa 1500 BCE). The Greeks identified specific poisons such as hemlock, a method of state execution, and the Greek word toxos (arrow) became the root of our modern science. The first scientific approach to the understanding of poisons and toxicology was the work during the late middle ages of Paracelsus. He formulated what were then revolutionary views that a specific toxic agent or "toxicon" caused specific dose-related effects. His principles have established the basis of modern pharmacology and toxicology. In 1700, Bernardo Ramazzini published the book De Morbis Artificum Diatriba (The Diseases of Workers) describing specific illnesses associated with certain labor, particularly metal workers exposed to mercury, lead, arsenic, and rock dust. Modern toxicology dates from development of the modern industrial chemical processes, the earliest involving an analytical method for arsenic by Marsh in 1836. Industrial organic chemicals were synthesized in the late 1800 s along with anesthetics and disinfectants. In 1908, Hamilton began the long study of occupational toxicology issues, and by WW I the scientific use of toxicants saw Haber creating war gases and defining time-dosage relationships that are used even today.

James, John T.↗

Inverse Compton Scattering in Mildly Relativistic Plasma

We investigated the effect of inverse Compton scattering in mildly relativistic static and moving plasmas with low optical depth using Monte Carlo simulations, and calculated the Sunyaev-Zel'dovich effect in the cosmic background radiation. Our semi-analytic method is based on a separation of photon diffusion in frequency and real space. We use Monte Carlo simulation to derive the intensity and frequency of the scattered photons for a monochromatic incoming radiation. The outgoing spectrum is determined by integrating over the spectrum of the incoming radiation using the intensity to determine the correct weight. This method makes it possible to study the emerging radiation as a function of frequency and direction. As a first application we have studied the effects of finite optical depth and gas infall on the Sunyaev-Zel'dovich effect (not possible with the extended Kompaneets equation) and discuss the parameter range in which the Boltzmann equation and its expansions can be used. For high temperature clusters (k(sub B)T(sub e) greater than or approximately equal to 15 keV) relativistic corrections based on a fifth order expansion of the extended Kompaneets equation seriously underestimate the Sunyaev-Zel'dovich effect at high frequencies. The contribution from plasma infall is less important for reasonable velocities. We give a convenient analytical expression for the dependence of the cross-over frequency on temperature, optical depth, and gas infall speed. Optical depth effects are often more important than relativistic corrections, and should be taken into account for high-precision work, but are smaller than the typical kinematic effect from cluster radial velocities.

Molnar, S. M.↗

Flow analysis and design optimization methods for nozzle afterbody of a hypersonic vehicle

This report summarizes the methods developed for the aerodynamic analysis and the shape optimization of the nozzle-afterbody section of a hypersonic vehicle. Initially, exhaust gases were assumed to be air. Internal-external flows around a single scramjet module were analyzed by solving the three dimensional Navier-Stokes equations. Then, exhaust gases were simulated by a cold mixture of Freon and Argon. Two different models were used to compute these multispecies flows as they mixed with the hypersonic airflow. Surface and off-surface properties were successfully compared with the experimental data. In the second phase of this project, the Aerodynamic Design Optimization with Sensitivity analysis (ADOS) was developed. Pre and post optimization sensitivity coefficients were derived and used in this quasi-analytical method. These coefficients were also used to predict inexpensively the flow field around a changed shape when the flow field of an unchanged shape was given. Starting with totally arbitrary initial afterbody shapes, independent computations were converged to the same optimum shape, which rendered the maximum axial thrust.

Baysal, Oktay↗

Experimental and Computational Investigation of Flow Boiling in Microgravity

This study explores use of Computational Fluid Dynamics (CFD) to predict near-saturated flow boiling of FC-72 in microgravity. The computational method employs transient analysis to predict detailed interfacial behavior and heat transfer characteristics along a rectangular channel heated along two opposite walls. Predicted results are validated against experimental temperature measurements and high-speed video images captured during a series of parabolic aircraft maneuvers for three sets of operating conditions which include variations of both mass velocity and wall heat flux. The computational method is based on the multi-phase volume of fluid (VOF) model, which is combined with appropriate phase change and turbulence models, and accounts for both shear-lift force on bubbles and conjugate heat transfer along the heating walls. A key advantage of the CFD method is ability to capture details that are very difficult to measure experimentally, including detailed spatial variations of bubble shape, void fraction, mixture fluid temperature, liquid velocity, and vapor velocity, results for which are presented for each of the three test cases. Different flow regimes predicted along the heated length exhibit a number of dominant mechanisms including bubble nucleation, bubble growth, coalescence, vapor blankets, interfacial waviness, and residual liquid sub-layer, all of which agree well with experiment. Vapor velocity is shown to increase appreciably along the heated length because of increased void fraction, while liquid velocity experiences large fluctuations. Non-equilibrium effects are accentuated with increasing mass velocity, contributing minor deviations of fluid temperature from simulations compared to those predicted by the analytical method. Predicted wall temperature is fairly uniform in the middle of the heated length but increases in the entrance region, due to sensible heat transfer in the subcooled liquid, and decreases toward the exit, mostly because of flow acceleration resulting from increased void fraction.

two-phase cooling; flow boiling; CFD; microgravity↗

Monte Carlo Simulation of a Knudsen Effusion Mass Spectrometer Sampling System

Knudsen flow is easily simulated with a Monte Carlo method. In this study we develop a Visual Basic for Excel (VBA) code to simulate the molecular beam from a vaporizing solid. The system at NASA Glenn uses the restricted collimation method of Chatillon and colleagues, which consists of two apertures between the effusion cell and the ionizer. The diameter of the first aperture is smaller than the diameter of the effusion cell orifice, so the ionizer effectively only sees inside the effusion cell. The code is able to calculate the transmission coefficient through the cell orifice, through the cell orifice and the first aperture, and through the cell orifice and first and second aperatures. Calculated transmission coefficients through the cell orifice are compared to tabulated values to validate the code. Then transmission coefficients are calculated through the cell orifice and both apertures to the ionizer. This allows the geometry (aperture spacing and diameters) of the sampling system the sampling system to be optimized. Calculated transmission factors are also compared to literature values calculated via an analytic method.

Radke, Michael J.↗

An Overview of the Null-Field Method. II: Convergence and Numerical Stability

In this paper we provide an analysis of the convergence and numerical stability of the null-field method with discrete sources. We show that (i) if the null-field scheme is numerically stable then we can decide whether or not convergence can be achieved; (ii) if the null-field scheme is numerically unstable then we cannot draw any conclusion about the convergence issue; and (iii) the numerical stability is closely related to the property of a tangential system of radiating discrete sources to form a Riesz basis. Our numerical analysis indicates that for prolate spheroids and localized vector spherical wave functions, the null-field scheme is numerically unstable (this system of vector functions does not form a Riesz basis), while for distributed vector spherical wave functions, the numerical instability is not so pronounced (this system of discrete sources almost possesses the property of being a Riesz basis). We also describe an analytical method for computing the surface integrals in the framework of the conventional null-field method with localized vector spherical wave functions which increases the stability of the numerical scheme.

Null-field method↗

Resin transfer molding of textile preforms for aircraft structural applications

The NASA LaRC is conducting and supporting research to develop cost-effective fabrication methods that are applicable to primary composite aircraft structures. One of the most promising fabrication methods that has evolved is resin transfer molding (RTM) of dry textile material forms. RTM has been used for many years for secondary structures, but has received increased emphasis because it is an excellent method for applying resin to damage-tolerant textile preforms at low cost. Textile preforms based on processes such as weaving, braiding, knitting, stitching, and combinations of these have been shown to offer significant improvements in damage tolerance compared to laminated tape composites. The use of low-cost resins combined with textile preforms could provide a major breakthrough in achieving cost-effective composite aircraft structures. RTM uses resin in its lowest cost form, and storage and spoilage costs are minimal. Near net shape textile preforms are expected to be cost-effective because automated machines can be used to produce the preforms, post-cure operations such as machining and fastening are minimized, and material scrap rate may be reduced in comparison with traditional prepreg molding. The purpose of this paper is to discuss experimental and analytical techniques that are under development at NASA Langley to aid the engineer in developing RTM processes for airframe structural elements. Included are experimental techniques to characterize preform and resin behavior and analytical methods that were developed to predict resin flow and cure kinetics.

Hasko, Gregory H.↗

Accretion disk coronae in high-luminosity systems

We present the results of self-consistent models of Compton-heated accretion disk coronae. The models are calculated using a new method for computing monochromatic radiative transfer n two dimensions. The method splits the radiation into direct and scattered components. The direct radiation is computed by calculating the optical depth along rays, while transfer of the scattered radiation is approximated by flux-limited diffusion. The resulting code agrees with more accurate treatments to within 50%, and is highly efficient, making it practical for use in large hydrodynamic simulations. The coronal models are used to confirm the results of earlier work, and to extend it to higher luminosities. In contrast to earlier work, which found the outer disks to be shadowed by the inner corona at high luminosities, we find our results to form an almost continuous extension of the models at lower luminosities. This is due to the presence of multiply scattered radiation, which acts to partially offset the loss of direct radiation from the central source. Although the analytic methods derived at lower luminosities cannot be used to derive the coronal structure for L/L(sub Edd) approx. greater than 0.1, the results of the models are amenable to semiempirical fits. We also discuss possible observational consequences of the results for coronal veiling and line fluorescence from the disk.

Murray, Stephen D.↗

SSME - Materials and Methods for Addressing High-Pressure Hydrogen Embrittlement

From the humid, corrosion-friendly atmosphere of KSC, to the extreme heat of ascent, to the cold vacuum of space, the Space Shuttle faced one hostile environment after another. One of those harsh environments the hydrogen environment existed within the shuttle itself. Liquid hydrogen was the fuel that powered the shuttle s complex, powerful, and reusable main engine. Hydrogen provided the high specific impulse the bang per pound of fuel needed to perform the shuttle s heavy lifting duties. Hydrogen, however, was also a potential threat to the very metal of the propulsion system that used it. The diffusion of hydrogen atoms into a metal can make it more brittle and prone to cracking a process called hydrogen embrittlement. This effect can reduce the toughness of carefully selected and prepared materials. A concern that exposure to hydrogen might encourage crack growth was present from the beginning of the Space Shuttle Program, but the rationale for using hydrogen was compelling. This paper outlines the material characterization, anomaly resolution, and path to understanding of hydrogen embrittlement on superalloys through the course of the SSME program. Specific examples of nickel alloy turbine housings and single crystal turbine blades are addressed. The evolution of fracture mechanics analytical methods is also addressed.

Matejczk, Daniel↗

Computation of iodine species concentrations in water

During an evaluation of the use of iodine as a water disinfectant and the development of methods for measuring various iodine species in water onboard Space Freedom, it became necessary to compute the concentration of the various species based on equilibrium principles alone. Of particular concern was the case when various amounts of iodine, iodide, strong acid, and strong base are added to water. Such solutions can be used to evaluate the performance of various monitoring methods being considered. The authors of this paper present an overview of aqueous iodine chemistry, a set of nonlinear equations which can be used to model the above case, and a computer program for solving this system of equations using the Newton-Raphson method. The program was validated by comparing results over a range of concentrations and pH values with those previously presented by Gottardi for a given pH. Use of this program indicated that there are multiple roots to many cases and selecting an appropriate initial guess is important. Comparison of program results with laboratory results for the case when only iodine is added to water indicates the program gives high pH values for the iodine concentrations normally used for water disinfection. Extending the model to include the effects of iodate formation results in the computer pH values being closer to those observed, but the model with iodate does not agree well for the case in which base is added in addition to iodine to raise the pH. Potential explanations include failure to obtain equilibrium conditions in the lab, inaccuracies in published values for the equilibrium constants, and inadequate model of iodine chemistry and/or the lack of adequate analytical methods for measuring the various iodine species in water.

Schultz, John R.↗

Improving the Fidelity of Capability & Resource Weighting in A Probalistic Risk Assessment Model for Spaceflight

INTRODUCTION NASA’s Informing Mission Planning via Analysis of Complex Tradespaces (IMPACT) tool uses Probabilistic Risk Assessment (PRA) to provide an evidence-based, data-driven estimate of how medical system capabilities affect mission outcomes. IMPACT maps condition incidence to available resources thereby facilitating the calculation of outcome metrics that allow the estimation of mission medical risk. Conditions can be nominally categorized as either treated or untreated depending on the availability of necessary diagnostic and therapeutic capabilities. This categorization enables IMPACT to estimate the effect of various medical system configurations on mission outcomes such as crew mortality, disability, crew member down-time, return to duty/recovery, and need for evacuation. IMPACT currently employs an equal weighting, “partial credit” approach to define treatment in which each of the capabilities associated with a given condition contributes an equal amount to management of the condition. This feature enables IMPACT to report values in between “fully untreated” and “fully treated” based on the proportion of capabilities available within the model. However, as is normal in medical/clinical practice, not all individual capabilities contribute equally to medical care. For example, the ability to provide intramuscular epinephrine during an anaphylactic episode contributes more likelihood of overall management success than does the administration of oral diphenhydramine. We hypothesize that weighting the relative contribution of each capability to each specific condition will improve outcome prediction and therefore will better provide mission planners with more nuanced and accurate options when designing space medical systems. METHODS Using a five-point Fibonacci scaling sequence (1, 2, 3, 5, 8) subject matter experts from NASA’s Exploration Medical Capabilities (ExMC) element assigned relative contribution weighting values to each identified capability within IMPACT. Since the relative importance of each capability varies depending on the specific condition, the resulting “partial” weighting was completed for more than 1,600 individual weighting assignments for 666 capabilities across 121 conditions. Each assignment required three-physician concurrence based on the overall importance of the capability to the diagnosis and management of the condition being considered and the difficulty with which it could be improvised by the crew. Once complete, 100,000 IMPACT simulations were run for a 6-month Lunar mission with a 30-day surface stay to evaluate the effect of this modification of the model. RESULTS Partial weighting significantly decreased predicted task time loss (TTL), evacuation, and loss of crew life without causing significant changes to the recommended medical system design. CONCLUSIONS The paucity of real-world referent data to support long-duration space missions of this type limits the ability to judge one predictive analytics method as superior to another. However, since the proposed method significantly reduces and optimizes outcome risks—without changing the medical system design—incorporating a partial weighting methodology is likely to provide a more accurate and operationally-relevant representation of medical risk without compromising IMPACTs ability to inform overarching medical system requirements.

Steller JG↗

Sonic Boom Prediction and Minimization of the Douglas Reference OPT5 Configuration

Conventional CFD methods and grids do not yield adequate resolution of the complex shock flow pattern generated by a real aircraft geometry. As a result, a unique grid topology and supersonic flow solver was developed at Northrop Grumman based on the characteristic behavior of supersonic wave patterns emanating from the aircraft. Using this approach, it was possible to compute flow fields with adequate resolution several body lengths below the aircraft. In this region, three-dimensional effects are diminished and conventional two-dimensional modified linear theory (MLT) can be applied to estimate ground pressure signatures or sonic booms. To accommodate real aircraft geometries and alleviate the burdensome grid generation task, an implicit marching multi-block, multi-grid finite-volume Euler code was developed as the basis for the sonic boom prediction methodology. The Thomas two-dimensional extrapolation method is built into the Euler code so that ground signatures can be obtained quickly and efficiently with minimum computational effort suitable to the aircraft design environment. The loudness levels of these signatures can then be determined using a NASA generated noise code. Since the Euler code is a three-dimensional flow field solver, the complete circumferential region below the aircraft is computed. The extrapolation of all this field data from a cylinder of constant radius leads to the definition of the entire boom corridor occurring directly below and off to the side of the aircraft's flight path yielding an estimate for the entire noise "annoyance" corridor in miles as well as its magnitude. An automated multidisciplinary sonic boom design optimization software system was developed during the latter part of HSR Phase 1. Using this system, it was found that sonic boom signatures could be reduced through optimization of a variety of geometric aircraft parameters. This system uses a gradient based nonlinear optimizer as the driver in conjunction with a computationally efficient Euler CFD solver (NIIM3DSB) for computing the three-dimensional near-field characteristics of the aircraft. The intent of the design system is to identify and optimize geometric design variables that have a beneficial impact on the ground sonic boom. The system uses a simple wave drag data format to specify the aircraft geometry. The geometry is internally enhanced and analytic methods are used to generate marching grids suitable for the multi-block Euler solver. The Thomas extrapolation method is integrated into this system, and hence, the aircraft's centerline ground sonic boom signature is also automatically computed for a specified cruise altitude and yields the parameters necessary to evaluate the design function. The entire design system has been automated since the gradient based optimization software requires many flow analyses in order to obtain the required sensitivity derivatives for each design variable in order to converge on an optimal solution. Hence, once the problem is defined which includes defining the objective function and geometric and aerodynamic constraints, the system will automatically regenerate the perturbed geometry, the necessary grids, the Euler solution, and finally the ground sonic boom signature at the request of the optimizer.

Siclari, Michael J.↗

The ground station high-power traveling-wave tube

The M4040 is a 2-kw CW traveling-wave tube developed specifically for communications applications. It was used as the power amplifier in the Project Telstar ground transmitter. Analytic methods employed in the design of the electron gun, the beam collector, the RF circuit and the output match are presented. Typical performance characteristics and drawings of the tube subassemblies serve to describe the electrical and physical features of the M4040 traveling-wave tube. Methods used to inhibit oscillation and the effects of ion drainage and collector depression on tube performance are also discussed.

TRAVELING WAVE TUBE↗

Digital Distributed Parameter Model for Analysis of Unsteady Flow in Liquid-Filled Lines

A distributed parameter model for computing unsteady flow conditions in liquid systems is presented. The model is adapted for use in conjunction with a digital computer. The analytical method employed is essentially a solution synthesized from the effects of incremental step pressure pulses. The pressure pulses are generated because of incremental flow-rate changes that originate in a fluid system from a variety of sources, including the mechanical motion of the system structure. The pressure pulses propagate throughout the system at sonic velocity and are partially transmitted and reflected at each discontinuity. The velocity change caused by each pressure pulse is obtained from the relation for the characteristic acoustic impedance. Pressure and velocity time histories at any point in the system are obtained by a timewise summation of the contributions of the incremental pressure pulses passing that point. The analysis is presented in a form general enough to be applied to a variety of liquid-filled fluid systems. To illustrate the application of the method to a specific system, the response of a long, straight propellant line to a sinusoidal inlet flow and pressure perturbation is computed. Both a constant-cross-section line and a tapered line are analyzed in the example.

UNSTEADY FLOW↗