Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “parametric analysis”

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 1,207 records · Page 67

Optimal Control within the Context of Multidisciplinary Design, Analysis, and Optimization

Multidisciplinary design, analysis and optimization involves modeling the interactions of complex systems across a variety of disciplines. The optimization of such systems can be a computationally expensive exercise with multiple levels of nested nonlinear solvers running under an optimizer.The application of optimal control in project development often involves performing trajectory optimization for fixed vehicle designs or parametric sweeps across some key vehicle properties.This information is then relayed to the subsystem design teams who update their designs and relay some bulk characteristics back to the trajectory optimization procedure.This iteration is then repeated until the design closes.However, with increasing interest in more tightly coupled systems, such as electric and hybrid-electric aircraft propulsion and boundary layer ingestion, this process is prone to ignore subtle coupling between vehicle subsystem designs and vehicle operation on a given mission.Integrating trajectory optimization into a tightly coupled multidisciplinary design procedure can be computationally prohibitive, depending on the complexity of the subsystem analyses and the optimal control technique applied.To address these issues a new optimal control software tool, Dymos, has been developed.Dymos is built upon NASA's OpenMDAO software and can leverage its capabilities to efficiently compute gradients for the optimization and optimize complex models in parallel on distributed memory systems.This report provides some explanation into the numerical methods employed in Dymos and provides several use cases that demonstrate its performance on traditional optimal control problems and improvements ino techniques have been used extensively in recent decades to solve a variety of optimal control problems, typically in the form of aerospace vehicle trajectory optimization.

pseudospectral↗

A Method for Calculating Strain Energy Release Rates in Preliminary Design of Composite Skin/Stringer Debonding Under Multi-Axial Loading

Three simple procedures were developed to determine strain energy release rates, G, in composite skin/stringer specimens for various combinations of unaxial and biaxial (in-plane/out-of-plane) loading conditions. These procedures may be used for parametric design studies in such a way that only a few finite element computations will be necessary for a study of many load combinations. The results were compared with mixed mode strain energy release rates calculated directly from nonlinear two-dimensional plane-strain finite element analyses using the virtual crack closure technique. The first procedure involved solving three unknown parameters needed to determine the energy release rates. Good agreement was obtained when the external loads were used in the expression derived. This superposition technique was only applicable if the structure exhibits a linear load/deflection behavior. Consequently, a second technique was derived which was applicable in the case of nonlinear load/deformation behavior. The technique involved calculating six unknown parameters from a set of six simultaneous linear equations with data from six nonlinear analyses to determine the energy release rates. This procedure was not time efficient, and hence, less appealing. A third procedure was developed to calculate mixed mode energy release rates as a function of delamination lengths. This procedure required only one nonlinear finite element analysis of the specimen with a single delamination length to obtain a reference solution for the energy release rates and the scale factors. The delamination was extended in three separate linear models of the local area in the vicinity of the delamination subjected to unit loads to obtain the distribution of G with delamination lengths. This set of sub-problems was Although additional modeling effort is required to create the sub- models, this local technique is efficient for parametric studies.

Krueger, Ronald↗

Icing Physics Studies Using the 3D SIDRM Test Article: Aerodynamic and Supercooled Liquid Icing Analysis

In-flight icing is an important safety issue and is a factor that affects aircraft design and performance. Newer regulations are driving a need for improvements in airframe and engine icing simulation capability. Experimental data is required for development of icing physics models and simulation validation. To that end, this paper presents the analysis of the supercooled liquid icing data subset from tests conducted in 2022 at the NASA Icing Research Tunnel that studied both supercooled water and ice-crystal icing. The test article that was utilized replicated 3D geometrical features of an inter-compressor duct and strut region of a turbofan engine. The surfaces of the Simulated Inter-compressor Duct Research Model (SIDRM) can be heated to simulate the warm surfaces of the turbofan inter-compressor duct. The test article is instrumented with pressure taps, heaters, heat flux gauges, and thermocouples, while a 3D laser scanner, cameras, and a scale to measure ice mass were utilized to characterize the icing behavior. The aim of these tests was to generate ice accretions on the SIDRM test article under well-characterized supercooled liquid icing and ice crystal icing conditions. This paper discusses measurements related to aerodynamic testing and supercooled liquid icing tests that were conducted. Aerodynamic measurements were analyzed and compared to computational simulations and were found to be in good agreement for the range of airspeeds (50 to 230 knots) and angles of attack (0 to 4°) tested. Various parametric sweeps were conducted during the supercooled liquid icing portion of the test entry (cloud median volumetric diameter ranged from 15 to 90 µm, total air temperature from -3 to -17 °C, angle of attack from 0 to 4°, and accretion time from 5 to 20 min). These sweeps were performed to measure that parameter’s impact on ice accretion size, location (icing extent), characteristics (such as glaze/rime ice and shedding behavior), and test article surface temperature. Analysis of the test data showed that clouds composed of larger drops, colder air temperatures, smaller angles of attack, and longer spray times were the primary parameters that resulted in accretions with greater ice mass. Test article angle of attack and cloud droplet size influenced the location of ice accretion as these two parameters directly impact collection efficiency. With respect to icing characteristics, total air temperature dictated icing type, and smaller cloud drop size along with warmer air temperatures resulted in greater amounts of ice shedding. Surface temperature increased during ice accretion from the release of latent (fusion) heat, where total air temperature and cloud drop size impacted the amount of surface temperature change. The icing measurements collected during the SIDRM tests will be used to develop and validate 3D computational engine icing tools, such as GlennICE, that predictively assesses the onset and growth of ice. One of the goals of the sponsoring NASA project is to develop simulation models and tools that can assist in the design and certification of engines for flight in icing conditions in a cost effective way.

supercooled liquid icing↗

System and Propagation Availability Analysis for NASA's Advanced Air Transportation Technologies

This report summarizes the research on the System and Propagation Availability Analysis for NASA's project on Advanced Air Transportation Technologies (AATT). The objectives of the project were to determine the communication systems requirements and architecture, and to investigate the effect of propagation on the transmission of space information. In this report, results from the first year investigation are presented and limitations are highlighted. To study the propagation links, an understanding of the total system architecture is necessary since the links form the major component of the overall architecture. This study was conducted by way of analysis, modeling and simulation on the system communication links. The overall goals was to develop an understanding of the space communication requirements relevant to the AATT project, and then analyze the links taking into consideration system availability under adverse atmospheric weather conditions. This project began with a preliminary study of the end-to-end system architecture by modeling a representative communication system in MATLAB SIMULINK. Based on the defining concepts, the possibility of computer modeling was determined. The investigations continue with the parametric studies of the communication system architecture. These studies were also carried out with SIMULINK modeling and simulation. After a series of modifications, two end-to-end communication links were identified as the most probable models for the communication architecture. Link budget calculations were then performed in MATHCAD and MATLAB for the identified communication scenarios. A remarkable outcome of this project is the development of a graphic user interface (GUI) program for the computation of the link budget parameters in real time. Using this program, one can interactively compute the link budget requirements after supplying a few necessary parameters. It provides a framework for the eventual automation of several computations required in many experimental NASA missions. For the first year of this project, most of the stated objectives were accomplished. We were able to identify probable communication systems architectures, model and analyze several communication links, perform numerous simulation on different system models, and then develop a program for the link budget analysis. However, most of the work is still unfinished. The effect of propagation on the transmission of information in the identified communication channels has not been performed. Propagation effects cannot be studied until the system under consideration is identified and characterized. To study the propagation links, an understanding of the total communications architecture is necessary. It is important to mention that the original project was intended for two years and the results presented here are only for the first year of research. It is prudent therefore that these efforts be continued in order to obtain a complete picture of the system and propagation availability requirements.

Ugweje, Okechukwu C.↗

Flood hazard studies in Central Texas using orbital and suborbital remote sensing machinery

Central Texas is subject to infrequent, unusually intense rainstorms which cause extremely rapid runoff from drainage basins developed on the deeply dissected limestone and marl bedrock of the Edwards Plateau. One approach to flood hazard evaluation in this area is a parametric model relating flood hydrograph characteristics to quantitative geomorphic properties of the drainage basins. The preliminary model uses multiple regression techniques to predict potential peak flood discharge from basin magnitude, drainage density, and ruggedness number. After mapping small catchment networks from remote sensing imagery, input data for the model are generated by network digitization and analysis by a computer assisted routine of watershed analysis. The study evaluated the network resolution capabilities of the following data formats: (1) large-scale (1:24,000) topographic maps, employing Strahler's "method of v's," (2) standard low altitude black and white aerial photography (1:13,000 and 1:20,000 scales), (3) NASA - generated aerial infrared photography at scales ranging from 1:48,000 to 1:123,000, and (4) Skylab Earth Resources Experiment Package S-190A and S-190B sensors (1:750,000 and 1:500,000 respectively).

Baker, V. R.↗

Theoretical fan velocity distortions due to inlets and nozzles

Nonuniform velocity profiles imposed on the propulsion system fan can cause fan blade stresses and thrust losses. This paper presents a theoretical parametric study of the effects of inlets with 0 deg and 90 deg nozzle deflection on the velocity profile at a hypothetical fan. The parameters investigated are fan-to-nozzle spacing and inlet centerline offset. The interaction between the inlet and nozzle is also investigated. The study is made using a two-dimensional analysis

Hawk, J. D.↗

Interactive design and analysis of future large spacecraft concepts

Recent developments in interactive computer-aided design and analysis of future large spacecraft concepts are presented. The primary emphases are on automatic spacecraft modeling for lattice (truss-like) structures and simplified multidiscipline design and analysis modules. Capabilities and performance of the twenty some applications modules, the executive and data management software, and graphics display features are reviewed. Data generated for an earth-orbiting large diameter (750 meter) antenna satellite are used to illustrate current capabilities. A single user at an interactive terminal can create, design, analyze, and conduct parametric studies of earth-orbiting spacecraft with relative ease. The approach is particularly useful in the conceptual design phase of advanced space missions when a multiplicity of concepts must be evaluated in a cost-effective and timely manner

Garrett, L. B.↗

Autonomous frequency domain identification: Theory and experiment

The analysis, design, and on-orbit tuning of robust controllers require more information about the plant than simply a nominal estimate of the plant transfer function. Information is also required concerning the uncertainty in the nominal estimate, or more generally, the identification of a model set within which the true plant is known to lie. The identification methodology that was developed and experimentally demonstrated makes use of a simple but useful characterization of the model uncertainty based on the output error. This is a characterization of the additive uncertainty in the plant model, which has found considerable use in many robust control analysis and synthesis techniques. The identification process is initiated by a stochastic input u which is applied to the plant p giving rise to the output. Spectral estimation (h = P sub uy/P sub uu) is used as an estimate of p and the model order is estimated using the produce moment matrix (PMM) method. A parametric model unit direction vector p is then determined by curve fitting the spectral estimate to a rational transfer function. The additive uncertainty delta sub m = p - unit direction vector p is then estimated by the cross spectral estimate delta = P sub ue/P sub uu where e = y - unit direction vectory y is the output error, and unit direction vector y = unit direction vector pu is the computed output of the parametric model subjected to the actual input u. The experimental results demonstrate the curve fitting algorithm produces the reduced-order plant model which minimizes the additive uncertainty. The nominal transfer function estimate unit direction vector p and the estimate delta of the additive uncertainty delta sub m are subsequently available to be used for optimization of robust controller performance and stability.

Yam, Yeung↗

Evaluation of thermal loading on a methane injector at high pressure and temperature

Experimental and numerical analyses are conducted to determine the surface temperature on a methane fuel injector used to produce a high enthalpy test stream for a combustion-fed subscale wind tunnel facility. It was found that the ratio of the methane fuel injection velocity to the air stream velocity is a significant factor in the production of high injector surface temperatures which lead to rapid deterioration of the fuel injector structure. The numerical code utilized for the computational analysis was found to be representative of the experimentally measured data since the experimental trends were reproduced by the numerical simulation. The quantitative accuracy of the numerical predictions could not be assessed from the data gathered because of the difficulty of making a noninterfering injector surface temperature measurement. The numerical code can be used for parametric evaluation of combustor parameters and thus will serve as an important tool in the design of such fuel injector systems.

Harvin, Stephen F.↗

An analytical investigation of transient effects on rewetting of heated thin flat plates

The rewetting of a hot surface is a problem of prime importance in the microgravity application of heat pipe technology, where rewetting controls the time before operations can be re-established following depriming of a heat pipe. Rewetting is also important in the nuclear industry (in predicting behavior during loss-of-coolant accidents), as well as in the chemical and petrochemical industries. Recently Chan and Zhang have presented a closed-form solution for the determination of the rewetting speed of a liquid film flowing over a finite (but long) hot plate subject to uniform heating. Unfortunately, their physically unreasonable initial conditions preclude a meaningful analysis of start-up transient behavior. A new nondimensionalization and closed-form solution for an infinitely-long, uniformly-heated plate is presented. Realistic initial conditions (step change in temperature across the wetting front) and boundary conditions (no spatial temperature gradients infinitely far from the wetting front) are employed. The effects of parametric variation on the resulting simpler closed-form solution are presented and compared with the predictions of a 'quasi-steady' model. The time to reach steady-state rewetting is found to be a strong function of the initial dry-region plate temperature. For heated plates it is found that in most cases the effect of the transient response terms cannot be neglected, even for large times.

Platt, J. A.↗

A 3D-CFD code for accurate prediction of fluid flows and fluid forces in seals

Current and future turbomachinery requires advanced seal configurations to control leakage, inhibit mixing of incompatible fluids and to control the rotodynamic response. In recognition of a deficiency in the existing predictive methodology for seals, a seven year effort was established in 1990 by NASA's Office of Aeronautics Exploration and Technology, under the Earth-to-Orbit Propulsion program, to develop validated Computational Fluid Dynamics (CFD) concepts, codes and analyses for seals. The effort will provide NASA and the U.S. Aerospace Industry with advanced CFD scientific codes and industrial codes for analyzing and designing turbomachinery seals. An advanced 3D CFD cylindrical seal code has been developed, incorporating state-of-the-art computational methodology for flow analysis in straight, tapered and stepped seals. Relevant computational features of the code include: stationary/rotating coordinates, cylindrical and general Body Fitted Coordinates (BFC) systems, high order differencing schemes, colocated variable arrangement, advanced turbulence models, incompressible/compressible flows, and moving grids. This paper presents the current status of code development, code demonstration for predicting rotordynamic coefficients, numerical parametric study of entrance loss coefficients for generic annular seals, and plans for code extensions to labyrinth, damping, and other seal configurations.

Athavale, M. M.↗

Analysis of woven and braided fabric reinforced composites

A general purpose micromechanics analysis that discretely models the yarn architecture within the textile repeating unit cell, was developed to predict overall, three dimensional, thermal and mechanical properties. This analytical technique was implemented in a user-friendly, personal computer-based, windows compatible code called Textile Composite Analysis for Design (TEXCAD). TEXCAD was used to analyze plain, 5-harness satin, and 8-harness satin weave composites along with 2-D braided and 2x2, 2-D triaxial braided composites. The calculated overall stiffnesses correlated well with available 3-D finite element results and test data for both the woven and the braided composites. Parametric studies were performed to investigate the effects of yarn size on the yarn crimp and the overall thermal and mechanical constants for plain weave composites. The effects of braid angle were investigated for the 2-D braided composites. Finally, the effects of fiber volume fraction on the yarn undulations and the thermal and mechanical properties of 2x2, 2-D triaxial braided composites were also investigated.

Naik, Rajiv A.↗

Hybrid Solid Oxide Fuel Cell/Gas Turbine System Design for High Altitude Long Endurance Aerospace Missions

A system level analysis, inclusive of mass, is carried out for a cryogenic hydrogen fueled hybrid solid oxide fuel cell and bottoming gas turbine (SOFC/GT) power system. The system is designed to provide primary or secondary electrical power for an unmanned aerial vehicle (UAV) over a high altitude, long endurance mission. The net power level and altitude are parametrically varied to examine their effect on total system mass. Some of the more important technology parameters, including turbomachinery efficiencies and the SOFC area specific resistance, are also studied for their effect on total system mass. Finally, two different solid oxide cell designs are compared to show the importance of the individual solid oxide cell design on the overall system. We show that for long mission durations of 10 days or more, the fuel mass savings resulting from the high efficiency of a SOFC/GT system more than offset the larger powerplant mass resulting from the low specific power of the SOFC/GT system. These missions therefore favor high efficiency, low power density systems, characteristics typical of fuel cell systems in general.

Himansu, Ananda↗

Navigation accuracy analysis for a Halley intercept mission

A ballistic intercept mission to Comet Halley is currently under consideration by the United States. This paper describes the navigation system and the navigation strategy which would be employed in such a mission, assuming a launch in the summer of 1985 and an arrival in March of 1986. Spacecraft comet-relative orbit determination accuracies are presented as functions of time for the baseline navigation strategy, along with the results of a number of parametric sensitivity studies involving parameters such as data frequencies, data accuracies, stochastic acceleration levels, cometary ephemeris uncertainties, maneuver execution errors, and encounter date.

Wood, L. J.↗

Finite element analysis of the end notched flexure specimen for measuring Mode II fracture toughness

The paper presents a finite element analysis of the end-notched flexure (ENF) test specimen for Mode II interlaminar fracture testing of composite materials. Virtual crack closure and compliance techniques employed to calculate strain energy release rates from linear elastic two-dimensional analysis indicate that the ENF specimen is a pure Mode II fracture test within the constraints of small deflection theory. Furthermore, the ENF fracture specimen is shown to be relatively insensitive to process-induced cracks, offset from the laminate midplane. Frictional effects are investigated by including the contact problem in the finite element model. A parametric study investigating the influence of delamination length, span, thickness, and material properties assessed the accuracy of beam theory expressions for compliance and strain energy release rate, GII. Finite element results indicate that data reduction schemes based upon beam theory underestimate GII by approximately 20-40 percent for typical unidirectional graphite fiber composite test specimen geometries. Consequently, an improved data reduction scheme is proposed.

Gillespie, J. W., Jr.↗

Analysis of cracked laminates with holes using the boundary force method

The boundary force method is an indirect boundary element technique for the analysis of composite laminates with holes and cracks which uses the orthotropic elasticity solution for a concentrated horizontal and vertical force applied at a given point in a cracked infinite sheet. The orthotropic formulation of the method is verified by comparison with accepted solutions for a counter-crack specimen that is subjected to uniaxial tension. Parametric studies are conducted for two configurations for which no orthotropic solutions are available: cracks emanating from a circular hole, and a four-hole cracked specimen; both are loaded in uniaxial tension.

Tan, P. W.↗

Cavity Heating Experiments Supporting Shuttle Columbia Accident Investigation

The two-color thermographic phosphor method has been used to map the local heating augmentation of scaled idealized cavities at conditions simulating the windward surface of the Shuttle Orbiter Columbia during flight STS-107. Two experiments initiated in support of the Columbia Accident Investigation were conducted in the Langley 20-Inch Mach 6 Tunnel. Generally, the first test series evaluated open (length-to-depth less than 10) rectangular cavity geometries proposed as possible damage scenarios resulting from foam and ice impact during launch at several discrete locations on the vehicle windward surface, though some closed (length-to-depth greater than 13) geometries were briefly examined. The second test series was designed to parametrically evaluate heating augmentation in closed rectangular cavities. The tests were conducted under laminar cavity entry conditions over a range of local boundary layer edge-flow parameters typical of re-entry. Cavity design parameters were developed using laminar computational predictions, while the experimental boundary layer state conditions were inferred from the heating measurements. An analysis of the aeroheating caused by cavities allowed exclusion of non-breeching damage from the possible loss scenarios being considered during the investigation.

Everhart, Joel L.↗

A Parametric Model of Shoulder Articulation for Virtual Assessment of Space Suit Fit

Shoulder injury is one of the most severe risks that have the potential to impair crewmembers' performance and health in long duration space flight. Overall, 64% of crewmembers experience shoulder pain after extra-vehicular training in a space suit, and 14% of symptomatic crewmembers require surgical repair (Williams & Johnson, 2003). Suboptimal suit fit, in particular at the shoulder region, has been identified as one of the predominant risk factors. However, traditional suit fit assessments and laser scans represent only a single person's data, and thus may not be generalized across wide variations of body shapes and poses. The aim of this work is to develop a software tool based on a statistical analysis of a large dataset of crewmember body shapes. This tool can accurately predict the skin deformation and shape variations for any body size and shoulder pose for a target population, from which the geometry can be exported and evaluated against suit models in commercial CAD software. A preliminary software tool was developed by statistically analyzing 150 body shapes matched with body dimension ranges specified in the Human-Systems Integration Requirements of NASA ("baseline model"). Further, the baseline model was incorporated with shoulder joint articulation ("articulation model"), using additional subjects scanned in a variety of shoulder poses across a pre-specified range of motion. Scan data was cleaned and aligned using body landmarks. The skin deformation patterns were dimensionally reduced and the co-variation with shoulder angles was analyzed. A software tool is currently in development and will be presented in the final proceeding. This tool would allow suit engineers to parametrically generate body shapes in strategically targeted anthropometry dimensions and shoulder poses. This would also enable virtual fit assessments, with which the contact volume and clearance between the suit and body surface can be predictively quantified at reduced time and cost.

Kim, K. Han↗