Engineering topics

Tiwari, S. N.

Publications and source records attributed to Tiwari, S. N..

At least 19 records

NASA NTRS2004

Residual Strength Characterization of a Curved Integrally-Stiffened Panel

Over the years, Finite-element fracture simulation methodology has been very well established at NASA Langley to predict the residual strength of damaged aircraft structures. This methodology has been experimentally verified at NASA Langley for structures ranging from laboratory coupons up to full-scale built-up structural components with single and multiple-site damage cracking. The methodology uses the critical crack-tip-opening-angle (CTOA) fracture criterion to characterize the fracture behavior of the material. The CTOA fracture criterion assumes that stable crack growth occurs when the crack-tip angle reaches a constant critical value. The use of the CTOA criterion requires an elastic-plastic, finite-element analysis. The critical CTOA value is determined by simulating fracture behavior in laboratory specimens, such as a compact specimen, to obtain the angle that best fits the observed test behavior. The critical CTOA value appears to be independent of loading, crack length, and in-plane dimensions. However, it is a function of material thickness and local crack-front constraint. Modeling the local constraint requires either a three-dimensional analysis or a two-dimensional analysis with an approximation to account for the constraint effects. In recent times as the aircraft industry is leaning towards monolithic structures with the intension of reducing part count and manufacturing cost, there has been a consistent effort at NASA Langley to extend critical CTOA based numerical methodology in the analysis of integrally-stiffened panels. In this regard, a series of fracture tests were conducted on curved aluminum-alloy integrally-stiffened panels. These curved panels were subjected to uniaxial tension and pressure loading. During the test, applied load-crack extension, out-of-plane displacements and local deformations around the crack tip region were measured. Compact and middle-crack tension specimens were tested to determine the critical angle (psi(sub c) using three-dimensional code (ZIP3D) and the plane-strain core height (h(sub c)) using two-dimensional code (STAGS).

Seshadri, B. R.↗

NASA NTRS2004

Numerical Investigation of Dual-Mode Scramjet Combustor with Large Upstream Interaction

Dual-mode scramjet combustor configuration with significant upstream interaction is investigated numerically, The possibility of scaling the domain to accelerate the convergence and reduce the computational time is explored. The supersonic combustor configuration was selected to provide an understanding of key features of upstream interaction and to identify physical and numerical issues relating to modeling of dual-mode configurations. The numerical analysis was performed with vitiated air at freestream Math number of 2.5 using hydrogen as the sonic injectant. Results are presented for two-dimensional models and a three-dimensional jet-to-jet symmetric geometry. Comparisons are made with experimental results. Two-dimensional and three-dimensional results show substantial oblique shock train reaching upstream of the fuel injectors. Flow characteristics slow numerical convergence, while the upstream interaction slowly increases with further iterations. As the flow field develops, the symmetric assumption breaks down. A large separation zone develops and extends further upstream of the step. This asymmetric flow structure is not seen in the experimental data. Results obtained using a sub-scale domain (both two-dimensional and three-dimensional) qualitatively recover the flow physics obtained from full-scale simulations. All results show that numerical modeling using a scaled geometry provides good agreement with full-scale numerical results and experimental results for this configuration. This study supports the argument that numerical scaling is useful in simulating dual-mode scramjet combustor flowfields and could provide an excellent convergence acceleration technique for dual-mode simulations.

Mohieldin, T. O.↗

NASA NTRS2003

Residual Strength Analyses of Monolithic Structures

Finite-element fracture simulation methodology predicts the residual strength of damaged aircraft structures. The methodology uses the critical crack-tip-opening-angle (CTOA) fracture criterion to characterize the fracture behavior of the material. The CTOA fracture criterion assumes that stable crack growth occurs when the crack-tip angle reaches a constant critical value. The use of the CTOA criterion requires an elastic- plastic, finite-element analysis. The critical CTOA value is determined by simulating fracture behavior in laboratory specimens, such as a compact specimen, to obtain the angle that best fits the observed test behavior. The critical CTOA value appears to be independent of loading, crack length, and in-plane dimensions. However, it is a function of material thickness and local crack-front constraint. Modeling the local constraint requires either a three-dimensional analysis or a two-dimensional analysis with an approximation to account for the constraint effects. In recent times as the aircraft industry is leaning towards monolithic structures with the intention of reducing part count and manufacturing cost, there has been a consistent effort at NASA Langley to extend critical CTOA based numerical methodology in the analysis of integrally-stiffened panels.In this regard, a series of fracture tests were conducted on both flat and curved aluminum alloy integrally-stiffened panels. These flat panels were subjected to uniaxial tension and during the test, applied load-crack extension, out-of-plane displacements and local deformations around the crack tip region were measured. Compact and middle-crack tension specimens were tested to determine the critical angle (wc) using three-dimensional code (ZIP3D) and the plane-strain core height (hJ using two-dimensional code (STAGS). These values were then used in the STAGS analysis to predict the fracture behavior of the integrally-stiffened panels. The analyses modeled stable tearing, buckling, and crack branching at the integral stiffener using different values of critical CTOA for different material thicknesses and orientation. Comparisons were made between measured and predicted load-crack extension, out-of-plane displacements and local deformations around the crack tip region. Simultaneously, three-dimensional capabilities to model crack branching and to monitor stable crack growth of multiple cracks in a large thick integrally-stiffened flat panels were implemented in three-dimensional finite element code (ZIP3D) and tested by analyzing the integrally-stiffened panels tested at Alcoa. The residual strength of the panels predicted from STAGS and ZP3D code compared very well with experimental data. In recent times, STAGS software has been updated with new features and now one can have combinations of solid and shell elements in the residual strength analysis of integrally-stiffened panels.

Forth, Scott↗

NASA NTRS2002

Institute for Scientific and Educational Technology (ISET)-Education, Research and Training Programs in Engineering and Sciences

The 'Institute for Scientific and Educational Technology' has been established to provide a mechanism through which universities and other research organizations may cooperate with one another and with different government agencies and industrial organizations to further and promote research, education, and training programs in science, engineering, and related fields. This effort has been undertaken consistent with the national vision to 'promote excellence in America s educational system through enhancing and expanding scientific and technological competence.' The specific programs are directed in promoting and achieving excellence for individuals at all levels (elementary and secondary schools, undergraduate and graduate education, and postdoctoral and faculty research). The program is consistent with the existing activities of the Institute for Computational and Applied Mechanics (ICAM) and the American Society for Engineering Education (ASEE) at NASA Langley Research Center (LaRC). The efforts will be directed to embark on other research, education, and training activities in various fields of engineering, scientific, and educational technologies. The specific objectives of the present program may be outlined briefly as follows: 1) Cooperate in the various research, education, and technology programs of the Office of Education at LaRC. 2) Develop procedures for interactions between precollege, college, and graduate students, and between faculty and students at all levels. 3) Direct efforts to increase the participation by women and minorities in educational programs at all levels. 4) Enhance existing activities of ICAM and ASEE in education, research, and training of graduate students and faculty. 5) Invite distinguished scholars as appropriate and consistent with ISET goals to spend their summers and/or sabbaticals at NASA Langley andor ODU and interact with different researchers and graduate students. Perform research and administrative activities as needed to carry out the above mentioned activities. 6) The implementation of various activities of the ISET programs is carried out through cooperative efforts between Old Dominion University (ODU) and the Office of Education at LaRC. At present, major efforts are directed on the following ISET Programs: ICAM Programs, Academic Programs, Educational Research, Outreach Programs, Educational Technology and Cooperative Programs. These programs are described in the following sections.

Tiwari, S. N.↗

NASA NTRS2001

Entropy generation method to quantify thermal comfort

The present paper presents a thermodynamic approach to assess the quality of human-thermal environment interaction and quantify thermal comfort. The approach involves development of entropy generation term by applying second law of thermodynamics to the combined human-environment system. The entropy generation term combines both human thermal physiological responses and thermal environmental variables to provide an objective measure of thermal comfort. The original concepts and definitions form the basis for establishing the mathematical relationship between thermal comfort and entropy generation term. As a result of logic and deterministic approach, an Objective Thermal Comfort Index (OTCI) is defined and established as a function of entropy generation. In order to verify the entropy-based thermal comfort model, human thermal physiological responses due to changes in ambient conditions are simulated using a well established and validated human thermal model developed at the Institute of Environmental Research of Kansas State University (KSU). The finite element based KSU human thermal computer model is being utilized as a "Computational Environmental Chamber" to conduct series of simulations to examine the human thermal responses to different environmental conditions. The output from the simulation, which include human thermal responses and input data consisting of environmental conditions are fed into the thermal comfort model. Continuous monitoring of thermal comfort in comfortable and extreme environmental conditions is demonstrated. The Objective Thermal Comfort values obtained from the entropy-based model are validated against regression based Predicted Mean Vote (PMV) values. Using the corresponding air temperatures and vapor pressures that were used in the computer simulation in the regression equation generates the PMV values. The preliminary results indicate that the OTCI and PMV values correlate well under ideal conditions. However, an experimental study is needed in the future to fully establish the validity of the OTCI formula and the model. One of the practical applications of this index is that could it be integrated in thermal control systems to develop human-centered environmental control systems for potential use in aircraft, mass transit vehicles, intelligent building systems, and space vehicles.

Models, Biological↗

NASA NTRS2001

Fabrication of Thermoplastic Composite Laminates Having Film Interleaves By Automated Fiber Placement

Experiments were carried out at the NASA Langley Research Center automated Fiber placement facility to determine an optimal process for the fabrication of composite materials having polymer film interleaves. A series of experiments was conducted to determine an optimal process for the composite prior to investigation of a process to fabricate laminates with polymer films. The results of the composite tests indicated that a well-consolidated, void-free laminate could be attained. Preliminary interleaf processing trials were then conducted to establish some broad guidelines for film processing. The primary finding of these initial studies was that a two-stage process was necessary in order to process these materials adequately. A screening experiment was then performed to determine the relative influence of the process variables on the quality of the film interface as determined by the wedge peel test method. Parameters that were found to be of minor influence on specimen quality were subsequently held at fixed values enabling a more rapid determination of an optimal process. Optimization studies were then performed by varying the remaining parameters at three film melt processing rates. The resulting peel data were fitted with quadratic response surfaces. Additional specimens were fabricated at levels of high peel strength as predicted by the regression models in an attempt to gage the accuracy of the predicted response and to assess the repeatability of the process. The overall results indicate that quality laminates having film interleaves can be successfully and repeatably fabricated by automated fiber placement.

Hulcher, A. B.↗

NASA NTRS1999

Thermodynamic Modeling and Analysis of Human Stress Response

A novel approach based on the second law of thermodynamics is developed to investigate the psychophysiology and quantify human stress level. Two types of stresses (thermal and mental) are examined. A Unified Stress Response Theory (USRT) is developed under the new proposed field of study called Engineering Psychophysiology. The USRT is used to investigate both thermal and mental stresses from a holistic (human body as a whole) and thermodynamic viewpoint. The original concepts and definitions are established as postulates which form the basis for thermodynamic approach to quantify human stress level. An Objective Thermal Stress Index (OTSI) is developed by applying the second law of thermodynamics to the human thermal system to quantify thermal stress or dis- comfort in the human body. The human thermal model based on finite element method is implemented. It is utilized as a "Computational Environmental Chamber" to conduct series of simulations to examine the human thermal stress responses under different environmental conditions. An innovative hybrid technique is developed to analyze human thermal behavior based on series of human-environment interaction simulations. Continuous monitoring of thermal stress is demonstrated with the help of OTSI. It is well established that the human thermal system obeys the second law of thermodynamics. Further, the OTSI is validated against the experimental data. Regarding mental stress, an Objective Mental Stress Index (OMSI) is developed by applying the Maxwell relations of thermodynamics to the combined thermal and cardiovascular system in the human body. The OMSI is utilized to demonstrate the technique of monitoring mental stress continuously and is validated with the help of series of experimental studies. Although the OMSI indicates the level of mental stress, it provides a strong thermodynamic and mathematical relationship between activities of thermal and cardiovascular systems of the human body.

Boregowda, S. C.↗

NASA NTRS1999

Numerical Investigation of Hydrogen and Kerosene Combustion in Supersonic Air Streams

The effect of mixing schemes on the combustion of both gaseous hydrogen and liquid kerosene is investigated. Injecting pilot gaseous hydrogen parallel to the supersonic incoming air tends to maintain the stabilization of the main liquid kerosene, which is normally injected. Also the maximum kerosene equivalence ratio that can maintain stable flame can be increased by increasing the pilot energy level. The wedge flame holding contributes to an increased kerosene combustion efficiency by the generation of shock-jet interaction.

Taha, A. A.↗

NASA NTRS1998

Hybrid IR-Gas Heater for Automated Tow Placement

A hybrid infrared-hot gas heat source has been developed and tested for the NASA Langley Tow Placement Facility. The IR heat source provides supplemental heat to the nip region. The additional heat is intended to reduce the need for conduction heating by the compaction roller, which causes (he roller to stick to the panel surface. Initial bench scale testing was performed to identify the most effective means of focusing IR energy into the nip region. A compact lamp and reflector that placed the lamp as close to the nip point as possible was found to deliver the highest heat flux in the nip region. A prototype heater was installed on the NASA Langley tow placement robot. Panels placed with a 400 C (sticking) compaction roller gave DCB initiation toughness numbers comparable to those reported for autoclave processed panels but were found to have unexpectedly high void contents. Placement with compaction roller temperatures that prevented roller sticking resulted in mode I fracture toughness approximately 70% that reported for autoclave processed panels. The variability in strength among specimens placed with reduced roller temperature was found to be greatly reduced, which implies that use of supplemental nip point heat may improve the robustness of the tow placement process. Use of the IR heat source permitted placement with a compaction roller temperature that would have resulted in negligible interfacial strength with the hot gas torch alone. The roller temperature reductions eliminated the need for the robot operator to attend placement operations.

Grenoble, Ray W.↗

NASA NTRS1997

Activities and Accomplishments of ICAM

A brief historical background on establishing the Institute for Computational and Applied Mechanics (ICAM) is presented and basic goals and objectives are discussed. It is emphasized that the goal of the ICAM has been to develop and maintain a self-sustaining center of excellence in computational methods at Old Dominion University (ODU). Information is provided on funding sources and budget disposition, recent activities and accomplishments, list of graduate students supported on the program, and number of students who received graduate degrees (M.S. as well as Ph.D.). Information is also provided on research coordination with various scientists and engineers, and on different reports specifically written for ICAM. ICAM has been supported, in part, by NASA Langley Research Center through Grant NAG-1-363. This report constitutes the final report for ICAM for the period ending December 1996. The grant has been monitored by the University Affairs Officers at NASA Langley.

Tiwari, S. N.↗

NASA NTRS1997

Bidirectional Reflectance Functions for Application to Earth Radiation Budget Studies

Reflected solar radiative fluxes emerging for the Earth's top of the atmosphere are inferred from satellite broadband radiance measurements by applying bidirectional reflectance functions (BDRFs) to account for the anisotropy of the radiation field. BDRF's are dependent upon the viewing geometry (i.e. solar zenith angle, view zenith angle, and relative azimuth angle), the amount and type of cloud cover, the condition of the intervening atmosphere, and the reflectance characteristics of the underlying surface. A set of operational Earth Radiation Budget Experiment (ERBE) BDRFs is available which was developed from the Nimbus 7 ERB (Earth Radiation Budget) scanner data for a three-angle grid system, An improved set of bidirectional reflectance is required for mission planning and data analysis of future earth radiation budget instruments, such as the Clouds and Earth's Radiant Energy System (CERES), and for the enhancement of existing radiation budget data products. This study presents an analytic expression for BDRFs formulated by applying a fit to the ERBE operational model tabulations. A set of model coefficients applicable to any viewing condition is computed for an overcast and a clear sky scene over four geographical surface types: ocean, land, snow, and desert, and partly cloudy scenes over ocean and land. The models are smooth in terms of the directional angles and adhere to the principle of reciprocity, i.e., they are invariant with respect to the interchange of the incoming and outgoing directional angles. The analytic BDRFs and the radiance standard deviations are compared with the operational ERBE models and validated with ERBE data. The clear ocean model is validated with Dlhopolsky's clear ocean model. Dlhopolsky developed a BDRF of higher angular resolution for clear sky ocean from ERBE radiances. Additionally, the effectiveness of the models accounting for anisotropy for various viewing directions is tested with the ERBE along tract data. An area viewed from nadir and from the side give two different radiance measurements but should yield the same flux when converted by the BDRF. The analytic BDRFs are in very good qualitative agreement with the ERBE models. The overcast scenes exhibit constant retrieved albedo over viewing zenith angles for solar zenith angles less than 60 degrees. The clear ocean model does not produce constant retrieved albedo over viewing zenith angles but gives an improvement over the ERBE operational clear sky ocean BDRF.

Manalo-Smith, N.↗

NASA NTRS1997

Sensitivity Analysis and Optimization of Aerodynamic Configurations with Blend Surfaces

A novel (geometrical) parametrization procedure using solutions to a suitably chosen fourth order partial differential equation is used to define a class of airplane configurations. Inclusive in this definition are surface grids, volume grids, and grid sensitivity. The general airplane configuration has wing, fuselage, vertical tail and horizontal tail. The design variables are incorporated into the boundary conditions, and the solution is expressed as a Fourier series. The fuselage has circular cross section, and the radius is an algebraic function of four design parameters and an independent computational variable. Volume grids are obtained through an application of the Control Point Form method. A graphic interface software is developed which dynamically changes the surface of the airplane configuration with the change in input design variable. The software is made user friendly and is targeted towards the initial conceptual development of any aerodynamic configurations. Grid sensitivity with respect to surface design parameters and aerodynamic sensitivity coefficients based on potential flow is obtained using an Automatic Differentiation precompiler software tool ADIFOR. Aerodynamic shape optimization of the complete aircraft with twenty four design variables is performed. Unstructured and structured volume grids and Euler solutions are obtained with standard software to demonstrate the feasibility of the new surface definition.

Thomas, A. M.↗

NASA NTRS1997

Variational Methods in Design Optimization and Sensitivity Analysis for Two-Dimensional Euler Equations

Variational methods (VM) sensitivity analysis employed to derive the costate (adjoint) equations, the transversality conditions, and the functional sensitivity derivatives. In the derivation of the sensitivity equations, the variational methods use the generalized calculus of variations, in which the variable boundary is considered as the design function. The converged solution of the state equations together with the converged solution of the costate equations are integrated along the domain boundary to uniquely determine the functional sensitivity derivatives with respect to the design function. The application of the variational methods to aerodynamic shape optimization problems is demonstrated for internal flow problems at supersonic Mach number range. The study shows, that while maintaining the accuracy of the functional sensitivity derivatives within the reasonable range for engineering prediction purposes, the variational methods show a substantial gain in computational efficiency, i.e., computer time and memory, when compared with the finite difference sensitivity analysis.

Ibrahim, A. H.↗

NASA NTRS1996

Nonlinear stability of supersonic jets

The stability calculations made for a shock-free supersonic jet using the model based on parabolized stability equations are presented. In this analysis the large scale structures, which play a dominant role in the mixing as well as the noise radiated, are modeled as instability waves. This model takes into consideration non-parallel flow effects and also nonlinear interaction of the instability waves. The stability calculations have been performed for different frequencies and mode numbers over a range of jet operating temperatures. Comparisons are made, where appropriate, with the solutions to Rayleigh's equation (linear, inviscid analysis with the assumption of parallel flow). The comparison of the solutions obtained using the two approaches show very good agreement.

Tiwari, S. N.↗

NASA NTRS1996

Variational Methods in Sensitivity Analysis and Optimization for Aerodynamic Applications

Variational methods (VM) sensitivity analysis, which is the continuous alternative to the discrete sensitivity analysis, is employed to derive the costate (adjoint) equations, the transversality conditions, and the functional sensitivity derivatives. In the derivation of the sensitivity equations, the variational methods use the generalized calculus of variations, in which the variable boundary is considered as the design function. The converged solution of the state equations together with the converged solution of the costate equations are integrated along the domain boundary to uniquely determine the functional sensitivity derivatives with respect to the design function. The determination of the sensitivity derivatives of the performance index or functional entails the coupled solutions of the state and costate equations. As the stable and converged numerical solution of the costate equations with their boundary conditions are a priori unknown, numerical stability analysis is performed on both the state and costate equations. Thereafter, based on the amplification factors obtained by solving the generalized eigenvalue equations, the stability behavior of the costate equations is discussed and compared with the state (Euler) equations. The stability analysis of the costate equations suggests that the converged and stable solution of the costate equation is possible only if the computational domain of the costate equations is transformed to take into account the reverse flow nature of the costate equations. The application of the variational methods to aerodynamic shape optimization problems is demonstrated for internal flow problems at supersonic Mach number range. The study shows, that while maintaining the accuracy of the functional sensitivity derivatives within the reasonable range for engineering prediction purposes, the variational methods show a substantial gain in computational efficiency, i.e., computer time and memory, when compared with the finite difference sensitivity analysis.

Ibrahim, A. H.↗

NASA NTRS1996

Investigation of shock-induced combustion past blunt projectiles

A numerical study is conducted to simulate shock-induced combustion in premixed hydrogen-air mixtures at various free-stream conditions and parameters. Two-dimensional axisymmetric, reacting viscous flow over blunt projectiles is computed to study shock-induced combustion at Mach 5.11 and Mach 6.46 in hydrogen-air mixture. A seven-species, seven reactions finite rate hydrogen-air chemical reaction mechanism is used combined with a finite-difference, shock-fitting method to solve the complete set of Navier-Stokes and species conservation equations. The study has allowed an improved understanding of the physics of shock-induced combustion over blunt projectiles and the numerical results can now be explained more readily with one-dimensional wave-interaction model.

Ahuja, J. K.↗

NASA NTRS1995

Nonlinear Stability of Supersonic Jets

This paper presents stability calculations made for a shock-free supersonic jet using the model based on parabolized stability equations. In this analysis the large-scale structures, which play a dominant role in the mixing as well as the noise radiated, are modeled as instability waves. This model takes into consideration non-parallel flow effects and also nonlinear interaction of the instability waves. The stability calculations have been performed for different frequencies and mode numbers over a range of jet operating temperatures. Comparisons are made, where appropriate, with the solutions to Rayleigh's equation (linear, inviscid analysis with the assumption of parallel flow). The comparison of the solutions obtained using the two approaches show very good agreement.

Bhat, T. R. S.↗

NASA NTRS1994

Numerical simulation of shock-induced combustion past blunt bodies using shock-fitting technique

Two-dimensional axisymmetric, reacting viscous flow over blunt projectiles is computed to study shock-induced combustion at Mach 5.11 and Mach 6.46 in hydrogen-air mixture. A finite-difference, shock-fitting method is used to solve the complete set of Navier-Stokes and species conservation equations. In this approach, the bow shock represents a boundary of the computational domain and is treated as a discontinuity across which Rankine-Hugoniot conditions are applied. All interior details of the flow such as compression waves, reaction front, and the wall boundary layer are captured automatically in the solution. Since shock-fitting approach reduces the amount of artificial dissipation, all the intricate details of the flow are captured much more clearly than has been possible with the shock-capturing approach. This has allowed an improved understanding of the physics of shock-induced combustion over blunt projectiles and the numerical results can now be explained more readily with one-dimensional wave-interaction model than before.

Ahuja, J. K.↗