SEARCH · Engineering Papers

Results for “Nonlinear projection”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

339 records · Page 19

All sources NASA NTRS DOE OSTI

NASA NTRS2023

Lessons Learned from Large-Scale Aerospace Structural Testing

Large-scale testing of aerospace structures is frequently the final step in a development project to validate the structural performance, and that step typically involves a large cost and time investment. To ensure that the testing provides the required data, avoiding errors that can result in an unsuccessful test and failure to meet objectives is critical. Five lessons learned are presented herein to provide insight to those conducting tests in order to help them avoid known pitfalls that may result in an unsuccessful test. Five large-scale tests are described, and include two composite wing tests, a composite hybrid-wing body center section test, a full-scale 27.5-ft diameter metallic barrel test, and an 8-ft diameter metallic barrel test. Problems identified during the testing and mitigation approaches to solve the problems are presented, then the lessons learned are identified and discussed.

Lessons Learned↗

NASA NTRS1999

Advanced Techniques for Assessment of Postural and Locomotor Ataxia, Spatial Orientation, and Gaze Stability

In addition to adapting to microgravity, major neurovestibular problems of space flight include postflight difficulties with standing, walking, turning corners, and other activities that require stable upright posture and gaze stability. These difficulties inhibit astronauts' ability to stand or escape from their vehicle during emergencies. The long-ter7n goal of the NSBRI is the development of countermeasures to ameliorate the effects of long duration space flight. These countermeasures must be tested with valid and reliable tools. This project aims to develop quantitative, parametric approaches for assessing gaze stability and spatial orientation during normal gait and when gait is perturbed. Two of this year's most important findings concern head fixation distance and ideal trajectory analysis. During a normal cycle of walking the head moves up and down linearly. A simultaneous angular pitching motion of the head keeps it aligned toward an imaginary point in space at a distance of about one meter in front of a subject and along the line of march. This distance is called the head fixation distance. Head fixation distance provides the fundamental framework necessary for understanding the functional significance of the vestibular reflexes that couple head motion to eye motion. This framework facilitates the intelligent design of counter-measures for the effects of exposure to microgravity upon the vestibular ocular reflexes. Ideal trajectory analysis is a simple candidate countermeasure based upon quantifying body sway during repeated up and down stair stepping. It provides one number that estimates the body sway deviation from an ideal sinusoidal body sway trajectory normalized on the subject's height. This concept has been developed with NSBRI funding in less than one year. These findings are explained in more detail below. Compared to assessments of the vestibuo-ocular reflex, analysis of vestibular effects on locomotor function is relatively less well developed and quantified. We are improving this situation by applying methodologies such as nonlinear orbital stability to quantify responses and by using multivariate statistical approaches to link together the responses across separate tests. In this way we can exploit the information available and increase the ability to discriminate between normal and pathological responses. Measures of stability and orientation are compared to measures such as dynamic visual acuity and with balance function tests. The responses of normal human subjects and of patients having well documented pathophysiologies are being characterized. When these studies are completed, we should have a clearer idea about normal and abnormal patterns of eye, head, and body movements during locomotion and their stability in a wide range of environments. We plan eventually to use this information to validate the efficacy of candidate neurovestibular and neuromuscular rehabilitative techniques. Some representative studies made during this year are summarized.

Wall, Conrad., III↗

NASA NTRS1995

Effects of target fragmentation on evaluation of LET spectra from space radiation in low-earth orbit (LEO) environment: impact on SEU predictions

Recent improvements in the radiation transport code HZETRN/BRYNTRN and galactic cosmic ray environmental model have provided an opportunity to investigate the effects of target fragmentation on estimates of single event upset (SEU) rates for spacecraft memory devices. Since target fragments are mostly of very low energy, an SEU prediction model has been derived in terms of particle energy rather than linear energy transfer (LET) to account for nonlinear relationship between range and energy. Predictions are made for SEU rates observed on two Shuttle flights, each at low and high inclination orbit. Corrections due to track structure effects are made for both high energy ions with track structure larger than device sensitive volume and for low energy ions with dense track where charge recombination is important. Results indicate contributions from target fragments are relatively important at large shield depths (or any thick structure material) and at low inclination orbit. Consequently, a more consistent set of predictions for upset rates observed in these two flights is reached when compared to an earlier analysis with CREME model. It is also observed that the errors produced by assuming linear relationship in range and energy in the earlier analysis have fortuitously canceled out the errors for not considering target fragmentation and track structure effects.

STS-51 Shuttle Project↗

NASA NTRSDate not supplied

Assessment of Coupling One-Dimension and Three-Dimension CFD in Conjugate Heat Transfer Modeling and Simulations for Cooling of Turbine Engine

The CHT (conjugate heat transfer) modeling and simulation work fits under a bigger area in the TTT (Transformational Tools and Technologies) project called “Cooling Tools and Technologies for Combustor and Turbine.” The existing validation data sets such as TTT GRC (Glenn Research Center) THX (Turbulent Heat FLUX) IV experiments will be used to assess the solid/fluid heat transfer development and validation. The experimental configuration was a subsonic convergent nozzle with a square exit blowing a heated exhaust over an effusion-cooled deck. The simulation solutions were obtained and compared to experimental measurements of the profiles of flow velocities and temperatures, as well as deck surface temperatures. Cooling air blowing ratios of 0, 1, and 2 were considered. The three-dimensional cooling tubes were modelled by the medial axis of the cylindrical geometries due to large length over width ratios. It represented a typical coupling of one-dimension and three-dimension simulations of conjugate heat transfer approaches. The surface temperatures on the solid surfaces and the profiles of velocities and temperatures near the surfaces compared to the experimental data were presented.

nonlinear sub-filter stresses↗

NASA NTRS1991

Control algorithm implementation for a redundant degree of freedom manipulator

This project's purpose is to develop and implement control algorithms for a kinematically redundant robotic manipulator. The manipulator is being developed concurrently by Odetics Inc., under internal research and development funding. This SBIR contract supports algorithm conception, development, and simulation, as well as software implementation and integration with the manipulator hardware. The Odetics Dexterous Manipulator is a lightweight, high strength, modular manipulator being developed for space and commercial applications. It has seven fully active degrees of freedom, is electrically powered, and is fully operational in 1 G. The manipulator consists of five self-contained modules. These modules join via simple quick-disconnect couplings and self-mating connectors which allow rapid assembly/disassembly for reconfiguration, transport, or servicing. Each joint incorporates a unique drive train design which provides zero backlash operation, is insensitive to wear, and is single fault tolerant to motor or servo amplifier failure. The sensing system is also designed to be single fault tolerant. Although the initial prototype is not space qualified, the design is well-suited to meeting space qualification requirements. The control algorithm design approach is to develop a hierarchical system with well defined access and interfaces at each level. The high level endpoint/configuration control algorithm transforms manipulator endpoint position/orientation commands to joint angle commands, providing task space motion. At the same time, the kinematic redundancy is resolved by controlling the configuration (pose) of the manipulator, using several different optimizing criteria. The center level of the hierarchy servos the joints to their commanded trajectories using both linear feedback and model-based nonlinear control techniques. The lowest control level uses sensed joint torque to close torque servo loops, with the goal of improving the manipulator dynamic behavior. The control algorithms are subjected to a dynamic simulation before implementation.

Cohan, Steve↗

NASA NTRS2004

Relaxation and Preconditioning for High Order Discontinuous Galerkin Methods with Applications to Aeroacoustics and High Speed Flows

This project is about the investigation of the development of the discontinuous Galerkin finite element methods, for general geometry and triangulations, for solving convection dominated problems, with applications to aeroacoustics. Other related issues in high order WENO finite difference and finite volume methods have also been investigated. methods are two classes of high order, high resolution methods suitable for convection dominated simulations with possible discontinuous or sharp gradient solutions. In [18], we first review these two classes of methods, pointing out their similarities and differences in algorithm formulation, theoretical properties, implementation issues, applicability, and relative advantages. We then present some quantitative comparisons of the third order finite volume WENO methods and discontinuous Galerkin methods for a series of test problems to assess their relative merits in accuracy and CPU timing. In [3], we review the development of the Runge-Kutta discontinuous Galerkin (RKDG) methods for non-linear convection-dominated problems. These robust and accurate methods have made their way into the main stream of computational fluid dynamics and are quickly finding use in a wide variety of applications. They combine a special class of Runge-Kutta time discretizations, that allows the method to be non-linearly stable regardless of its accuracy, with a finite element space discretization by discontinuous approximations, that incorporates the ideas of numerical fluxes and slope limiters coined during the remarkable development of the high-resolution finite difference and finite volume schemes. The resulting RKDG methods are stable, high-order accurate, and highly parallelizable schemes that can easily handle complicated geometries and boundary conditions. We review the theoretical and algorithmic aspects of these methods and show several applications including nonlinear conservation laws, the compressible and incompressible Navier-Stokes equations, and Hamilton-Jacobi-like equations.

Shu, Chi-Wang↗

NASA NTRS1998

High Temperature Mechanical Characterization of Ceramic Matrix Composites

A high temperature mechanical characterization laboratory has been assembled at NASA Lewis Research Center. One contribution of this work is to test ceramic matrix composite specimens in tension in environmental extremes. Two high temperature tensile testing systems were assembled. The systems were assembled based on the performance and experience of other laboratories and meeting projected service conditions for the materials in question. The systems use frames with an electric actuator and a center screw. A PC based data acquisition and analysis system is used to collect and analyze the data. Mechanical extensometers are used to measure specimen strain. Thermocouples, placed near the specimen, are used to measure the specimen gage section temperature. The system for testing in air has a resistance element furnace with molybdenum disilicide elements and pneumatic grips with water cooling attached to hydraulic alignment devices. The system for testing in an inert gas has a graphite resistance element furnace in a chamber with rigidly mounted, water cooled, hydraulically actuated grips. Unidirectional SiC fiber reinforced reaction bonded Si3N4 and triaxially woven, two dimensional, SiC fiber reinforced enhanced SiC composites were tested in unidirectional tension. Theories for predicting the Young's modulus, modulus near the ultimate strength, first matrix cracking stress, and ultimate strength were applied and evaluated for suitability in predicting the mechanical behavior of SiC/RBSN and enhanced SiC/SiC composites. The SiC/RBSN composite exhibited pseudo tough behavior (increased area under the stress/strain curve) from 22 C to 1500 C. The rule of mixtures provides a good estimate of the Young's modulus of the SiC/RBSN composite using the constituent properties from room temperature to 1440 C for short term static tensile tests in air or nitrogen. The rule of mixtures significantly overestimates the secondary modulus near the ultimate strength. The ACK theory provides the best approximation of the first matrix cracking stress when residual stresses are ignored. The theory of Cao and Thouless, based on Weibull statistics, gave the best prediction for the composite ultimate strength. The enhanced SiC/SiC composite exhibited nonlinear stress/strain behavior from 24 C to 1370 C in air with increased ultimate strain when compared to monolithic SiC. The theory of Yang and Chou with the assumption of a frictional fiber/matrix interface provided the best estimate of the Young's modulus. The theory of Cao and Thouless gave the best estimate for the ultimate strength.

Gyekenyesi, John Z.↗

NASA NTRS2002

Vapor transport growth of organic solids in microgravity and unit gravity: some comparisons and results to date

Thin films of an organic nonlinear optical (NLO) material, N, N-dimethyl-p-(2,2-dicyanovinyl) aniline (DCVA), have been grown in space and on the ground by physical vapor transport in an effusive ampoule arrangement. The thin film growth technique developed on the ground is a direct result of information gleaned from experiments in microgravity. This paper covers the results of our experimental investigations for establishing "ideal" terrestrial conditions for deposition of a DCVA film. The active control during the deposition process was exercised by three deposition variables: the material source temperature, the background pressure external to the growth ampoule and the substrate temperature. Successful growth occurred when the difference in temperature between the source material and the copper substrate was 14 degrees C and the background nitrogen pressure was such that the transport was either diffusive or convective. A qualitative diffusion limited boundary was estimated to occur at a pressure of approximately 20 torr. We have probed the DCVA thin films with visible-near infrared reflection absorption spectroscopy, polarized Fourier transform infrared spectrometry, differential interference contrast optical microscopy, and stylus profilometry.

short duration↗

NASA NTRS2004

Fast Query-Optimized Kernel-Machine Classification

A recently developed algorithm performs kernel-machine classification via incremental approximate nearest support vectors. The algorithm implements support-vector machines (SVMs) at speeds 10 to 100 times those attainable by use of conventional SVM algorithms. The algorithm offers potential benefits for classification of images, recognition of speech, recognition of handwriting, and diverse other applications in which there are requirements to discern patterns in large sets of data. SVMs constitute a subset of kernel machines (KMs), which have become popular as models for machine learning and, more specifically, for automated classification of input data on the basis of labeled training data. While similar in many ways to k-nearest-neighbors (k-NN) models and artificial neural networks (ANNs), SVMs tend to be more accurate. Using representations that scale only linearly in the numbers of training examples, while exploring nonlinear (kernelized) feature spaces that are exponentially larger than the original input dimensionality, KMs elegantly and practically overcome the classic curse of dimensionality. However, the price that one must pay for the power of KMs is that query-time complexity scales linearly with the number of training examples, making KMs often orders of magnitude more computationally expensive than are ANNs, decision trees, and other popular machine learning alternatives. The present algorithm treats an SVM classifier as a special form of a k-NN. The algorithm is based partly on an empirical observation that one can often achieve the same classification as that of an exact KM by using only small fraction of the nearest support vectors (SVs) of a query. The exact KM output is a weighted sum over the kernel values between the query and the SVs. In this algorithm, the KM output is approximated with a k-NN classifier, the output of which is a weighted sum only over the kernel values involving k selected SVs. Before query time, there are gathered statistics about how misleading the output of the k-NN model can be, relative to the outputs of the exact KM for a representative set of examples, for each possible k from 1 to the total number of SVs. From these statistics, there are derived upper and lower thresholds for each step k. These thresholds identify output levels for which the particular variant of the k-NN model already leans so strongly positively or negatively that a reversal in sign is unlikely, given the weaker SV neighbors still remaining. At query time, the partial output of each query is incrementally updated, stopping as soon as it exceeds the predetermined statistical thresholds of the current step. For an easy query, stopping can occur as early as step k = 1. For more difficult queries, stopping might not occur until nearly all SVs are touched. A key empirical observation is that this approach can tolerate very approximate nearest-neighbor orderings. In experiments, SVs and queries were projected to a subspace comprising the top few principal- component dimensions and neighbor orderings were computed in that subspace. This approach ensured that the overhead of the nearest-neighbor computations was insignificant, relative to that of the exact KM computation.

Mazzoni, Dominic↗

NASA NTRS2010

Orion Abort Flight Test

The purpose of NASA's Constellation project is to create the new generation of spacecraft for human flight to the International Space Station in low-earth orbit, the lunar surface, as well as for use in future deep-space exploration. One portion of the Constellation program was the development of the Orion crew exploration vehicle (CEV) to be used in spaceflight. The Orion spacecraft consists of a crew module, service module, space adapter and launch abort system. The crew module was designed to hold as many as six crew members. The Orion crew exploration vehicle is similar in design to the Apollo space capsules, although larger and more massive. The Flight Test Office is the responsible flight test organization for the launch abort system on the Orion crew exploration vehicle. The Flight Test Office originally proposed six tests that would demonstrate the use of the launch abort system. These flight tests were to be performed at the White Sands Missile Range in New Mexico and were similar in nature to the Apollo Little Joe II tests performed in the 1960s. The first flight test of the launch abort system was a pad abort (PA-1), that took place on 6 May 2010 at the White Sands Missile Range in New Mexico. Primary flight test objectives were to demonstrate the capability of the launch abort system to propel the crew module a safe distance away from a launch vehicle during a pad abort, to demonstrate the stability and control characteristics of the vehicle, and to determine the performance of the motors contained within the launch abort system. The focus of the PA-1 flight test was engineering development and data acquisition, not certification. In this presentation, a high level overview of the PA-1 vehicle is given, along with an overview of the Mobile Operations Facility and information on the White Sands tracking sites for radar & optics. Several lessons learned are presented, including detailed information on the lessons learned in the development of wind placards for flight. PA-1 flight data is shown, as well as a comparison of PA-1 flight data to nonlinear simulation Monte Carlo data.

Hayes, Peggy Sue↗

NASA NTRS2003

Preliminary Planar Formation: Flight Dynamics Near Sun-Earth L2 Point

NASA's Goddard Space Flight Center is planning a series of missions in the vicinity of the Sun-Earth L2 libration point. Some of these projects will involve a distributed space system of telescope spacecraft acting together as a single telescope for high-resolution. The individual telescopes will be configured in a plane, surrounding a hub, where the telescope plane can be aimed toward various astronomical targets of interest. In preparation for these missions, it is necessary to develop an improved understanding of the dynamical behavior of objects in a planar configuration near L2. The classical circular restricted three body problem is taken as the basis for the analysis. At first order, the motion of such a telescope relative to the hub is described by a system of linear second order differential equations. These equations are identical to the circular restricted problem's linear equations describing the hub motion about L2. Therefore, the fundamental frequencies, both parallel to and normal to the ecliptic plane, are the same for the relative telescope motion as for the hub motion. To maintain the telescope plane for the duration necessary for the planned observations, a halo-type orbit of the telescopes about the hub is investigated. By using a halo orbit, the individual telescopes remain in approximately the same plane over the observation duration. For such an orbit, the fundamental periods parallel to and normal to the ecliptic plane are forced to be the same by careful selection of the initial conditions in order to adjust the higher order forces. The relative amplitudes of the resulting oscillations are associated with the orientation of the telescope plane relative to the ecliptic. As in the circular restricted problem, initial conditions for the linearized equations must be selected so as not to excite the convergent or divergent linear modes. In a higher order analysis, the telescope relative motion equations include the effects of the position of the hub relative to L2. In this paper, the differential equations are developed through second order in the distance of the hub from the libration point. A modified Lindstedt-Poincad perturbation method is employed to construct the solution of these differential equations through that same order of magnitude. In the course of the solution process, relationships are determined between the initial conditions of the telescopes, selected in order to avoid resonance excitation. As the differential equations include the hub position, it is necessary to simultaneously develop the solution for the hub. As has been done in past analyses of the circular restricted problem, the hub position is written in a power series formulation in terms of its distance from L2. Then, in order to be included in the telescope equations, the hub solution is cast in terms of the nonlinear frequency of the relative telescope motion. In the course of the analysis, it is determined that the hub should also maintain a halo orbit - about L2. Additionally, relationships are formed between the initial conditions of the telescopes and the hub. These relationships may be used to associate sets of initial conditions with particular orientations of the telescope plane. The accuracy of the analytical solution is verified through various simulations and comparison to numerical integration of the differential equations. The results of the simulations are presented, along with a graphical representation of the relationships between the initial conditions of the telescopes and hub.

Segerman, Alan M.↗

NASA NTRSDate not supplied

Aeroservoelastic Control Law Development for the Integrated Adaptive Wing Technology Maturation Wind-Tunnel Test

The NASA Advanced Air Transport Technology (AATT) Project seeks to “explore and develop technologies and concepts for improved energy efficiency and environmental compatibility for fixed wing subsonic transports” [1]. The IAWTM project is a subproject within AATT and is a joint Boeing/NASA effort to, in part, demonstrate the active control systems necessary to employ high-aspect ratio wings for improved aerodynamic efficiency. The aircraft configuration under study for the IAWTM project is a generic jet transport based on the NASA Common Research Model (CRM) whose wing aspect ratio has been increased to 13.5 (as opposed to nominal CRM aspect ratio of 9) and whose development is given in Ref. [2]. A semi-span model of the high-aspect-ratio CRM with 10 trailing edge control surfaces is currently being fabricated and will be delivered to the NASA Langley Transonic Dynamics Tunnel (TDT). There, several tests will be conducted to demonstrate maneuver load alleviation (MLA), gust load alleviation (GLA), active flutter suppression (AFS), and drag reduction, all at transonic conditions. A layout of the wind tunnel model is shown in Fig. 1. The wind-tunnel model is a spar/pod wing design with a relatively rigid fuselage attached to a mounting structure (seen in Fig. 1), which in turn is attached to a load balance behind the tunnel wall. Shims exist that can be attached to force the load to go through the fuselage rather than through the balance so as not to damage the balance in the event of flutter. The spar/pod design is similar to that used in the Boeing truss-braced wing aeroelastic model tested in the TDT in 2013 [3]. There is a central spar running the span of the wing to which the pods are attached that form the aerodynamic shape of the wing. The pods are constructed of ribs for attachment to the spar and skins to create the outer mold line. There are ten active control surfaces that will be used for the test objectives discussed above. Three of these (two outboard ailerons and one inboard aileron) are high-speed electrohydraulic actuators, and the remaining seven are low-speed mini-plane flaps, which have slower servo motors that will be used primarily for drag and maneuver load reductions. Control laws must be developed for each test (MLA, GLA, AFS, and drag reduction) and must be able to work in tandem when two or more control strategies are being tested simultaneously (e.g., when control laws to mitigate maneuver loads and gust loads operate simultaneously). The development of these control laws requires accurate mathematical representations of the aerodynamics, structural dynamics, and control surface dynamics of the system. Taken together, this mathematical representation—which is often in state-space form— is the ASE model. Three low-order ASE models of differing fidelities were developed for the IAWTM test and will be presented here. The differences in fidelity appeared in the aerodynamic modeling—one model assumed purely linear aerodynamics, while the other two models used different techniques to linearize the aerodynamics about a nonlinear mean. Furthermore, simple observer-based full-state feedback controllers for each test (MLA, GLA, AFS, and drag reduction) were designed from these low-order ASE models and implemented in a higher-order, nonlinear computational ASE simulation to assess the effectiveness of each controller. Observers were designed from each of the three lower-order ASE models, and each (linear) observer was then used as an estimate of the nonlinear dynamics of the system. Using the estimated states and the full-state feedback controller, control surface commands were fed back to the higher-order, nonlinear ASE system. The insight gained during this development will be applied to control law development for the physical IAWTM wind-tunnel test.

Josiah M Waite↗

NASA NTRS2024

Test and Analysis of the 8-foot Diameter Cylindrical Sandwich Composite Test Article CTA8.2B: As part of the NASA Engineering and Safety Center Shell Buckling Knockdown Factor Project

The buckling response of thin-walled cylindrical structures have long been shown to be sensitive to various imperfections. That is, imperfections related to load introduction, boundary conditions, material property variation, and geometric deviations can all contribute to experimentally obtained buckling loads being lower than the buckling loads predicted for perfect representations of the same thin-walled cylindrical structure. To account for the influence of geometry-based imperfections of a thin-walled cylindrical structure on buckling response, a design factor (also known as a buckling knockdown factor) is used during the design phase. Guidelines for buckling knockdown factors are readily available in NASA SP-8007, Buckling of Thin Walled Circular Cylinders, which was updated in 2020. Prior to this recent update, NASA SP-8007 had not been updated since 1969, and since that time, computational analysis, experimental testing, and manufacturing methods have significantly improved. The NASA Engineering and Safety Center (NESC) Shell Buckling Knockdown Factor Project (SBKF, NESC Assessment 07-010-E) had the goal of developing buckling design recommendations for select classes of metallic and composite launch-vehicle structures [3]. Specifically, the test and analysis results from the fourth SBKF composite test article SBKF-P3-CYL-CTA8.2B, which is referred to as CTA8.2B, are described in this paper. This test was the fourth in a series of four tests on sandwich composite cylinders that will be used to experimentally validate analysis methods, which in turn can be used to develop new analysis-based shell buckling design guidelines for typical sandwich composite launch vehicle cylindrical structures. Publications describing previous test validation with analysis are available in Refs. CTA8.2B was designed to occupy one corner of the design space with approximately equal bending stiffnesses in the axial and transverse directions and have a larger ratio of radius to effective shell thickness. Each test article had a unique design and was designed to occupy a specific region of the design space. CTA8.2B was an 8-foot diameter honeycomb-core sandwich cylinder that was fabricated and tested at the NASA Marshall Space Flight Center (MSFC). The primary objectives of this test were to interrogate the structural capability of a composite test article, and to verify the test article design and analysis approach for cylinders subjected to axial compression and combined axial and bending loads. First, descriptions of the test article design and test are given in Section 2.0, and modeling and analysis methods used in support of the test article design and testing activities are described briefly in Section 3.0. Then, select test results are presented and compared to predicted results in Section 4.0, and concluding remarks are presented in Section 5.0. Finally, drawings of CTA8.2B and references used to support the test and analysis are provided in the Appendix.

Nonlinear analysis↗

NASA NTRS2015

Use of Generalized Fluid System Simulation Program (GFSSP) for Teaching and Performing Senior Design Projects at the Educational Institutions

This paper describes the experience of the authors in using the Generalized Fluid System Simulation Program (GFSSP) in teaching Design of Thermal Systems class at University of Alabama in Huntsville. GFSSP is a finite volume based thermo-fluid system network analysis code, developed at NASA/Marshall Space Flight Center, and is extensively used in NASA, Department of Defense, and aerospace industries for propulsion system design, analysis, and performance evaluation. The educational version of GFSSP is freely available to all US higher education institutions. The main purpose of the paper is to illustrate the utilization of this user-friendly code for the thermal systems design and fluid engineering courses and to encourage the instructors to utilize the code for the class assignments as well as senior design projects. The need for a generalized computer program for thermofluid analysis in a flow network has been felt for a long time in aerospace industries. Designers of thermofluid systems often need to know pressures, temperatures, flow rates, concentrations, and heat transfer rates at different parts of a flow circuit for steady state or transient conditions. Such applications occur in propulsion systems for tank pressurization, internal flow analysis of rocket engine turbopumps, chilldown of cryogenic tanks and transfer lines, and many other applications of gas-liquid systems involving fluid transients and conjugate heat and mass transfer. Computer resource requirements to perform time-dependent, three-dimensional Navier-Stokes computational fluid dynamic (CFD) analysis of such systems are prohibitive and therefore are not practical. Available commercial codes are generally suitable for steady state, single-phase incompressible flow. Because of the proprietary nature of such codes, it is not possible to extend their capability to satisfy the above-mentioned needs. Therefore, the Generalized Fluid System Simulation Program (GFSSP1) has been developed at NASA Marshall Space Flight Center (MSFC) as a general fluid flow system solver capable of handling phase changes, compressibility, mixture thermodynamics and transient operations. It also includes the capability to model external body forces such as gravity and centrifugal effects in a complex flow network. The objectives of GFSSP development are: a) to develop a robust and efficient numerical algorithm to solve a system of equations describing a flow network containing phase changes, mixing, and rotation; and b) to implement the algorithm in a structured, easy-to-use computer program. The analysis of thermofluid dynamics in a complex network requires resolution of the system into fluid nodes and branches, and solid nodes and conductors as shown in Figure 1. Figure 1 shows a schematic and GFSSP flow circuit of a counter-flow heat exchanger. Hot nitrogen gas is flowing through a pipe, colder nitrogen is flowing counter to the hot stream in the annulus pipe and heat transfer occurs through metal tubes. The problem considered is to calculate flowrates and temperature distributions in both streams. GFSSP has a unique data structure, as shown in Figure 2, that allows constructing all possible arrangements of a flow network with no limit on the number of elements. The elements of a flow network are boundary nodes where pressure and temperature are specified, internal nodes where pressure and temperature are calculated, and branches where flowrates are calculated. For conjugate heat transfer problems, there are three additional elements: solid node, ambient node, and conductor. The solid and fluid nodes are connected with solid-fluid conductors. GFSSP solves the conservation equations of mass and energy, and equation of state in internal nodes to calculate pressure, temperature and resident mass. The momentum conservation equation is solved in branches to calculate flowrate. It also solves for energy conservation equations to calculate temperatures of solid nodes. The equations are coupled and nonlinear; therefore, they are solved by an iterative numerical scheme. GFSSP employs a unique numerical scheme known as simultaneous adjustment with successive substitution (SASS), which is a combination of Newton-Raphson and successive substitution methods. The mass and momentum conservation equations and the equation of state are solved by the Newton-Raphson method while the conservation of energy and species are solved by the successive substitution method. GFSSP is linked with two thermodynamic property programs, GASP2 and WASP3 and GASPAK4, that provide thermodynamic and thermophysical properties of selected fluids. Both programs cover a range of pressure and temperature that allows fluid properties to be evaluated for liquid, liquid-vapor (saturation), and vapor region. GASP and WASP provide properties of 12 fluids. GASPAK includes a library of 36 fluids. GFSSP has three major parts. The first part is the graphical user interface (GUI), visual thermofluid analyzer of systems and components (VTASC). VTASC allows users to create a flow circuit by a 'point and click' paradigm. It creates the GFSSP input file after the completion of the model building process. GFSSP's GUI provides the users a platform to build and run their models. It also allows post-processing of results. The network flow circuit is first built using three basic elements: boundary node, internal node, and branch.

Majumdar, A. K.↗

NASA NTRS2017

The Cloud Feedback Model Intercomparison Project (CFMIP) contribution to CMIP6.

The primary objective of CFMIP is to inform future assessments of cloud feedbacks through improved understanding of cloud-climate feedback mechanisms and better evaluation of cloud processes and cloud feedbacks in climate models. However, the CFMIP approach is also increasingly being used to understand other aspects of climate change, and so a second objective has now been introduced, to improve understanding of circulation, regional-scale precipitation, and non-linear changes. CFMIP is supporting ongoing model inter-comparison activities by coordinating a hierarchy of targeted experiments for CMIP6, along with a set of cloud-related output diagnostics. CFMIP contributes primarily to addressing the CMIP6 questions 'How does the Earth system respond to forcing?' and 'What are the origins and consequences of systematic model biases?' and supports the activities of the WCRP Grand Challenge on Clouds, Circulation and Climate Sensitivity. A compact set of Tier 1 experiments is proposed for CMIP6 to address this question: (1) what are the physical mechanisms underlying the range of cloud feedbacks and cloud adjustments predicted by climate models, and which models have the most credible cloud feedbacks? Additional Tier 2 experiments are proposed to address the following questions. (2) Are cloud feedbacks consistent for climate cooling and warming, and if not, why? (3) How do cloud-radiative effects impact the structure, the strength and the variability of the general atmospheric circulation in present and future climates? (4) How do responses in the climate system due to changes in solar forcing differ from changes due to CO2, and is the response sensitive to the sign of the forcing? (5) To what extent is regional climate change per CO2 doubling state-dependent (non-linear), and why? (6) Are climate feedbacks during the 20th century different to those acting on long-term climate change and climate sensitivity? (7) How do regional climate responses (e.g. in precipitation) and their uncertainties in coupled models arise from the combination of different aspects of CO2 forcing and sea surface warming? CFMIP also proposes a number of additional model outputs in the CMIP DECK, CMIP6 Historical and CMIP6 CFMIP experiments, including COSP simulator outputs and process diagnostics to address the following questions. 1. How well do clouds and other relevant variables simulated by models agree with observations? 2. What physical processes and mechanisms are important for a credible simulation of clouds, cloud feedbacks and cloud adjustments in climate models? 3. Which models have the most credible representations of processes relevant to the simulation of clouds? 4. How do clouds and their changes interact with other elements of the climate system?

carbon dioxide↗