Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Iterative limiters”

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 811 records · Page 45
All sourcesNASA NTRSDOE OSTI
NASA NTRS2020

INVESTIGATING THE FEASIBILITY AND STABILITY FOR MODELING ACOUSTIC WAVE SCATTERING USING A TIME-DOMAIN BOUNDARY INTEGRAL EQUATION WITH IMPEDANCE BOUNDARY CONDITION

Reducing aircraft noise is a major objective in the field of computational aeroacoustics. When designing next generation quiet and environmentally friendly aircraft, it is important to be able to accurately and efficiently predict the acoustic scattering by an aircraft body from a given noise source. Acoustic liners are an effective tool for aircraft noise reduction and are characterized by a frequency-dependent impedance. Converted into the time-domain using Fourier transforms, an impedance boundary condition can be used to simulate the acoustic wave scattering by geometric bodies treated with acoustic liners. This work considers using either an impedance or an admittance (inverse of impedance) boundary condition to allow for acoustic scattering problems to be modeled with geometries consisting of both un-lined and lined surfaces. Three acoustic liner models are discussed: the Extended Helmholtz Resonator Model, the Three-Parameter Impedance Model, and the Broadband Impedance Model. In both the Helmholtz and Three-Parameter models, liner impedance is specified at a given frequency, whereas the Broadband model allows for the investigation of multiple frequencies simultaneously. The impedance and admittance boundary conditions for acoustic liners are derived for each model and coupled with a time-domain boundary integral equation. The scattering solution is obtained iteratively using a boundary element method with constant spatial and third-order temporal basis functions. Time-domain boundary integral equations are unfortunately prone to numerical instabilities due to resonant frequencies resulting from non-trivial solutions in the interior domain. When reformulated with the Burton-Miller method, the instabilities are eliminated. Using a Burton-Miller reformulation, the stability of the boundary element method assuming a liner boundary condition is assessed using eigenvalue analysis. The stability of each liner model is discussed, and it is shown that the Three-Parameter and Broadband models are sufficient for modeling an acoustic liner on the surface of scattering bodies. The Helmholtz model demonstrates strict limitations for stability, whereas the Three-Parameter and Broadband models are stable for most cases. Also included in this work is an assessment of the spatial accuracy of the time-domain boundary element method with respect to the surface element basis functions, as well as a performance study of the numerical algorithm.

computational mathematics↗
DOE OSTI2021

Flattening the Radial Temperature Profile across the Transformational Challenge Reactor Core

The Transformational Challenge Reactor (TCR) program is demonstrating an agile development approach to advanced nuclear reactor design, which has traditionally utilized a linear design process. In leveraging artificial intelligence, additive manufacturing, advanced materials, and cutting-edge modeling and simulation, the TCR program aims to minimize the high cost and lengthy deployment timelines now standard in the nuclear industry. Within a relatively short period of time, a robust and mature advanced gas-cooled reactor was iteratively designed under the TCR program, using these cutting-edge technologies. The TCR is a 3 MWt gas-cooled microreactor fueled with uranium nitride (UN) tristructural isotropic (TRISO) fuel particles. Though manufactured via traditional means, these UN TRISO particles are loaded into additively manufactured silicon carbide (SiC) cans [4]. Once loaded with TRISO particles, the SiC cans are densified using a chemical vapor infiltration process. The additively manufactured SiC enables significantly more freedom in the design of the fuel form than could ever be achieved using traditionally manufactured SiC. The helium coolant, pressurized to 5 MPa, enters the core at 300°C and nominally exits it at 500°C. Typically, the most thermally limiting components in any reactor design are the fuel assemblies in the core center. To provide a wide thermal margin in these central fuel assemblies, the flow may be biased toward the center of the core to more effectively cool these fuel assemblies with more power deposition and flatten the core’s radial temperature distribution. An analytical fluid model of the TCR core was developed to explore methods for biasing the flow away from the cooler outer fuel assemblies and towards the hotter inner ones. Higher-fidelity models developed in STAR-CCM+ 2020.3.1, a computational fluid dynamics code, were then utilized to verify the analytical model’s findings.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
NASA NTRS2014

Statistical Symbolic Execution with Informed Sampling

Symbolic execution techniques have been proposed recently for the probabilistic analysis of programs. These techniques seek to quantify the likelihood of reaching program events of interest, e.g., assert violations. They have many promising applications but have scalability issues due to high computational demand. To address this challenge, we propose a statistical symbolic execution technique that performs Monte Carlo sampling of the symbolic program paths and uses the obtained information for Bayesian estimation and hypothesis testing with respect to the probability of reaching the target events. To speed up the convergence of the statistical analysis, we propose Informed Sampling, an iterative symbolic execution that first explores the paths that have high statistical significance, prunes them from the state space and guides the execution towards less likely paths. The technique combines Bayesian estimation with a partial exact analysis for the pruned paths leading to provably improved convergence of the statistical analysis. We have implemented statistical symbolic execution with in- formed sampling in the Symbolic PathFinder tool. We show experimentally that the informed sampling obtains more precise results and converges faster than a purely statistical analysis and may also be more efficient than an exact symbolic analysis. When the latter does not terminate symbolic execution with informed sampling can give meaningful results under the same time and memory limits.

Reliability↗
DOE OSTI2025

Graph-based Reversible Evaluation and Tangents Library

GRETL is a C++ library for evaluation, re-evaluation and algorithmic differentiation of functional operations on an arbitrary computational graph with limited memory usage. Similar to popular machine learning frameworks in Python, like PyTorch and JAX, it tracks and stores both operations and output data as functions are evaluated. Once this composition of functions is built up, the entire chain of operations can be back propagated to compute sensitivities of the final result with respect to any number of inputs. In contrast to most machine learning applications, memory usage becomes the bottleneck for back propagation in many physics applications, especially for time-dependent PDEs. Dynamic check pointing becomes essential. An important distinguishing feature of GRETL is its ability to limit the maximum memory usage by automatically dynamic checkpointing the data output for each graph operation (see Wang, Moin, Iaccarino, 2009). During backpropagation, parts of the graph that are no longer in memory are automatically re-evaluated from upstream checkpointed states as needed for derivative sensitivity calculations (or more precisely, for vector-Jacobian products). GRETL is particularly beneficial for applications, such as coupled multi-physics, where deriving adjoint-based sensitivities and managing checkpoint memory across modules becomes onerous. Cases which can be readily handled by the GRETL library include: different time-integration algorithms per physics (e.g., coupled predictor-corrector algorithms, IMEX, etc.), sub-cycling, asynchronous integrators, state dependent timestep sizes, iterative solvers and coupling algorithms, controller algorithms, and more.

Tupek, MichaelR [Lawrence Livermore National Labor↗
DOE OSTI2026

Deep Gaussian process-based cost-aware batch Bayesian optimization for complex materials design campaigns

The accelerating pace and expanding scope of materials discovery demand optimization frameworks that efficiently navigate vast design spaces with complex response surfaces while judiciously allocating limited evaluation resources. We present a cost-aware, batch Bayesian optimization scheme powered by deep Gaussian process (DGP) surrogates and a heterotopic querying strategy. Our DGP surrogate, formed by stacking GP layers, models complex hierarchical relationships among high-dimensional compositional features and captures correlations across multiple target properties, propagating uncertainty through successive layers. We integrate evaluation cost into an upper-confidence-bound acquisition extension, which, together with heterotopic querying, proposes small batches of candidates in parallel, balancing exploration of under-characterized regions with exploitation of high-mean, low-variance predictions across correlated properties. Applied to refractory high-entropy alloys for high-temperature applications, our framework converges to optimal formulations in fewer iterations with cost-aware queries than conventional GP-based BO, highlighting the value of deep, uncertainty-aware, cost-sensitive strategies in materials campaigns.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Low Cost, High Volume, Carbon Fiber Precursor for Plasma Oxidation

Light weighting with carbon fiber is critical to improving the energy efficiency of wind turbines, airplanes, and automobiles, but the high cost of carbon fiber continues to limit market growth and applications. Half of this cost is in the raw materials. For the carbon fiber market to reach the next level of adoption, the cost of both the raw materials and conversion must be reduced. This project considered the combination of low-cost, textile-grade precursors combined with the increased efficiency of plasma oxidation to greatly reduce the cost of carbon fiber. The primary objective of this project was to make industrial-grade carbon fiber (550 ksi tensile strength, 35 Msi tensile modulus) from textile-grade PAN precursors, something that cannot be achieved with conventional oxidation and to date has not been achieved with any technology. To accomplish this, the project team iteratively optimized the process conditions in the plasma oxidation and carbonization stages to maximize the tensile properties of the resultant carbon fibers. This involved implementing experimental designs for both stages, whereby process conditions (e.g. temperature) were varied and correlated to fiber properties. Samples of oxidized PAN fibers (OPFs) were prepared using a small pilot-scale plasma oxidation oven at 4XTechnologies (4XT). The thermal, chemical, and physical properties of OPF samples were analyzed by the University of Tennessee (UTK) and Oak Ridge National Laboratory (ORNL). The OPF samples were then converted to carbon fiber by subsequent low- and high-temperature carbonization using a lab-scale furnace at ORNL. The resultant carbon fibers were tested for tensile strength also at ORNL, and the results were fed back into the experimental design. Of the hundreds of carbon fiber samples produced over the course of the project, the best sample had an average break strength of 479 ksi and modulus of 33 Msi. While the best tensile properties achieved fell short of the final project targets of 550 ksi for break strength and 35 Msi for modulus, the team demonstrated the potential to produce industrial-grade carbon fiber using textile-grade precursors. The best performing fiber samples had tensile properties that are 87% of the goal for tensile strength and 94% of the goal for tensile modulus. With additional optimization work, it is likely that the targets could be achieved, since only a fraction of the process parameter space was explored during this project. 4XT and its affiliate (4M) have already secured funding through private equity, and there is a CRADA (Cooperative Research and Development Agreement) project currently underway with the Carbon Fiber Technology Center (CFTF) at ORNL to continue this research. Additionally, while not an original goal of the project, the carbon fibers for most samples had diameters >8.5μm. This is considerably larger than commercially available carbon fiber products, which can impart unique properties to the composites (e.g. increased compressive strength) and decrease manufacturing costs even further. Moreover, these fibers were oxidized with plasma oxidation using residence times ranging <60mins, which would be considered fast for normal diameter fibers, but it is as much as 5x faster when factoring in the diameter of the fibers.

36 MATERIALS SCIENCE↗
NASA NTRSDate not supplied

Evaluating the Feasibility of Data Simulation for Radiation Research: Estimating RBE Values for HZE Radiation

As the drive for deep space exploration intensifies, a comprehensive understanding of the health effects of radiation exposure becomes paramount to the future of human space flight. However, epidemiological data for radiation exposure, particularly to high-energy (HZE) ions, is limited, partially due to the financial and logistical costs of radiation studies. As an alternative, this study aims to assess the viability of data simulation strategies to accurately model potential study parameters prior to utilizing laboratory conditions. This study estimates a relative biological effectiveness (RBE) factor based on the solid tumor data for outbred mice provided by Edmundson et al. 2020. Excess relative risk (ERR) models for HZE-Fe ions and gamma radiation were estimated using Poisson regression with Weibull models to represent the background solid tumor hazard without radiation. RBE values were calculated from the ratio of the heavy ion linear slope to the gamma linear slope. The parameters from these models were then applied to simulate iterations of 300 datasets across HZE-Fe doses of 0.05, 0.2, 0.4, and 0.75 Gy; gamma radiation doses of 0.75, 2, and 3 Gy, as well as an unirradiated control group. Sample size per dose varied from 100-500 mice across simulations (800-4000 mice total per dataset). 1500 datasets were generated total (300 each for sample sizes 800, 1600, 3200, and 4000). For each dataset, ERR per radiation type and RBE for HZE-Fe were calculated. The RBE from the Edmondson data was calculated to be 5.55. The RBEs from the simulated data converged around this value across the increasing sample sizes. For n = 4000, the mean RBE across the 300 datasets was 5.71 (95% CI: 5.47-5.95). These results suggest that RBEs calculated from simulated data are credible. Based on this exploration, data simulation is a viable method of testing radiation studies. It provides a method of testing study assumptions and refining research questions prior to embarking on costly laboratory experiments.

Sneha Mehta↗
DOE OSTI2024

Improving Cathode Testing with a High-Gradient Cryogenic Normal Conducting RF Photogun

Future electron accelerator applications such as X-ray free electron lasers and colliders are dependent on significantly increasing beam brightness. With the observation that linac beam manipulation’s best preservation of max brightness is at the cathode, we are incentivized to create an environment where we can study how to achieve the highest possible photogun brightness. In order to do so, we intend to extract beams from high-brightness photocathodes with the highest achievable accelerating gradients we can manage in a klystron-powered radiofrequency (RF) photogun. We utilize here cryogenic normal conducting cavities to achieve ultra-high gradients via limitation of breakdown rates (BDR). The low temperatures should also reduce cathode emittance by reducing the mean transverse energy (MTE) of electrons near the photoemission threshold. To this end, we have designed and produced a new CrYogenic Brightness-Optimized Radiofrequency Gun (CYBORG) for use in a new beamline at UCLA. We will introduce the enabling RF and photoemission physics as a primer for the new regime of high field low temperature cathodes we intend to enter. We further report the current status of the beamline commissioning, including the cooling of the photogun to 100 K, and producing 0.5 MW of RF feed power, which corresponds to cathode accelerating fields in the range of 80–90 MV/m. We further plan iterative improvements to both to 77 K and 1 MW corresponding to our ultimate goal >120 MV/m. Our discussion will include future beamline tests and the consideration of the initial realization of an ultra-high-gradient photoinjector concept.

Lawler, Gerard Emile (ORCID:0000000160113435)↗
DOE OSTI2022

CMB-S4: Forecasting Constraints on Primordial Gravitational Waves

CMB-S4—the next-generation ground-based cosmic microwave background (CMB) experiment—is set to significantly advance the sensitivity of CMB measurements and enhance our understanding of the origin and evolution of the universe. Among the science cases pursued with CMB-S4, the quest for detecting primordial gravitational waves is a central driver of the experimental design. This work details the development of a forecasting framework that includes a power-spectrum-based semianalytic projection tool, targeted explicitly toward optimizing constraints on the tensor-to-scalar ratio, r, in the presence of Galactic foregrounds and gravitational lensing of the CMB. This framework is unique in its direct use of information from the achieved performance of current Stage 2–3 CMB experiments to robustly forecast the science reach of upcoming CMB-polarization endeavors. The methodology allows for rapid iteration over experimental configurations and offers a flexible way to optimize the design of future experiments, given a desired scientific goal. To form a closed-loop process, we couple this semianalytic tool with map-based validation studies, which allow for the injection of additional complexity and verification of our forecasts with several independent analysis methods. We document multiple rounds of forecasts for CMB-S4 using this process and the resulting establishment of the current reference design of the primordial gravitational-wave component of the Stage-4 experiment, optimized to achieve our science goals of detecting primordial gravitational waves for r > 0.003 at greater than 5σ, or in the absence of a detection, of reaching an upper limit of r < 0.001 at 95% CL.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

FLORA Equilibrium and Stability Code Archive

FLORA solves, in a 2-D domain (radial and axial dimensions) with a specified azimuthal Fourier mode, for the linearized stability of a long, thin, axisymmetric plasma equilibrium in an applied magnetic field. Before the stability equation is solved, FLORA solves a set of simple equations for pressure balance that specify the equilibrium magnetic field and plasma pressure in the long-thin limit given a simplified description of the magnetic coils. It uses an initial-value method for the linear stability problem in which an equilibrium is given an initial perturbation to its magnetic field, and the temporal behavior of the perturbation is followed. The perturbation has been Fourier expanded in the azimuthal coordinate; each azimuthal mode must be examined separately. The complex partial differential equation of motion for the perturbed radial displacement of the field lines is solved as a coupled system of two real p.d.e.'s and the solution consists of two parts, the real part and the imaginary part. The system is solved by bringing the coupling terms in each equation to the right side and using an iterative technique.

Cohen, Bruce [Lawrence Livermore National Laborato↗
DOE OSTI2023

Model predictive current profile control in tokamaks by exploiting spatially moving electron cyclotron current drives

The plasma control systems in next-generation tokamaks like ITER will balance competing control objectives to achieve the desired level of performance in advanced scenarios while preventing magnetohydrodynamic instabilities and disruptions. During normal tokamak operation, the points of incidence of the electromagnetic waves generated by the electron cyclotron heating and current drives (EC H&CDs) are usually fixed in space. However, the points of incidence can be modified in real-time by changing the angles of the mirrors that reflect the EC H&CD waves. Altering the points of incidence, in turn, varies the ability of the plasma control system to regulate a plasma property. For instance, changing the EC H&CD wave incidence location may place the power demands necessary to achieve a particular plasma target within saturation limits. Therefore, using the EC H&CD deposition location, which is related to the EC H&CD mirror angle, as a supplementary controllable variable may facilitate access to a given target scenario. However, active scenario-control algorithms have not been designed so far to fully exploit this capability in real time. In this work, a model predictive controller that can handle actuation locations as control inputs is developed. In particular, the controller is designed to regulate both the auxiliary powers and the EC H&CD deposition locations in a pre-defined optimal sense to achieve the control objective of attaining and sustaining a target current profile. Further, the proposed controller is tested for a DIII-D tokamak scenario in nonlinear simulations using the Control Oriented Transport SIMulator (COTSIM).

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Uncertainty quantification for Bayesian active learning in rupture life prediction of ferritic steels

Abstract Three probabilistic methodologies are developed for predicting the long-term creep rupture life of 9–12 wt%Cr ferritic-martensitic steels using their chemical and processing parameters. The framework developed in this research strives to simultaneously make efficient inference along with associated risk, i.e., the uncertainty of estimation. The study highlights the limitations of applying probabilistic machine learning to model creep life and provides suggestions as to how this might be alleviated to make an efficient and accurate model with the evaluation of epistemic uncertainty of each prediction. Based on extensive experimentation, Gaussian Process Regression yielded more accurate inference ( $$Pearson\;correlation\;coefficent> 0.95$$ P e a r s o n c o r r e l a t i o n c o e f f i c e n t > 0.95 for the holdout test set) in addition to meaningful uncertainty estimate (i.e., coverage ranges from 94 to 98% for the test set) as compared to quantile regression and natural gradient boosting algorithm. Furthermore, the possibility of an active learning framework to iteratively explore the material space intelligently was demonstrated by simulating the experimental data collection process. This framework can be subsequently deployed to improve model performance or to explore new alloy domains with minimal experimental effort.

20 FOSSIL-FUELED POWER PLANTS↗
DOE OSTI2020

QCOR; A Language Extension Specification for the Heterogeneous Quantum-Classical Model of Computation

Quantum computing (QC) is an emerging computational paradigm that leverages the laws of quantum mechanics to perform elementary logic operations. Existing programming models for QC were designed with fault-tolerant hardware in mind, envisioning stand-alone applications. However, the susceptibility of near-term quantum computers to noise limits their stand-alone utility. To better leverage limited computational strengths of noisy quantum devices, hybrid algorithms have been suggested whereby quantum computers are used in tandem with their classical counterparts in a heterogeneous fashion. This modus operandi calls out for a programming model and a high-level programming language that natively and seamlessly supports heterogeneous quantum-classical hardware architectures in a single-source-code paradigm. Motivated by the lack of such a model, we introduce a language extension specification, called QCOR, which enables single-source quantum-classical programming. Programs written using the QCOR library–based language extensions can be compiled to produce functional hybrid binary executables. After defining QCOR’s programming model, memory model, and execution model, we discuss how QCOR enables variational, iterative, and feed-forward QC. Additionally, QCOR approaches quantum-classical computation in a hardware-agnostic heterogeneous fashion and strives to build on best practices of high-performance computing. The high level of abstraction in the language extension is intended to accelerate the adoption of QC by researchers familiar with classical high-performance computing.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1974

System design considerations for free-fall materials processing

The design constraints for orbiting materials processing systems are dominated by the limitations of the flight vehicle/crew and not by the processes themselves. Although weight, size and power consumption are all factors in the design of normal laboratory equipment, their importance is increased orders of magnitude when the equipment must be used in an orbital facility. As a result, equipment intended for space flight may have little resemblance to normal laboratory apparatus although the function to be performed may be identical. The same considerations influence the design of the experiment itself. The processing requirements must be carefully understood in terms of basic physical parameters rather than defined in terms of equipment operation. Preliminary experiments and analysis are much more vital to the design of a space experiment than they are on earth where iterative development is relatively easy. Examples of these various considerations are illustrated with examples from the M518 and MA-010 systems. While these are specific systems, the conclusions apply to the design of flight materials processing systems both present and future.

Seidensticker, R. G.↗
DOE OSTI2024

Plasma performance and operational space with an RMP-ELM suppressed edge

Abstract The operational space and global performance of plasmas with edge-localized modes (ELMs) suppressed by resonant magnetic perturbations (RMPs) are surveyed by comparing AUG, DIII-D, EAST, and KSTAR stationary operating points. RMP-ELM suppression is achieved over a range of plasma currents, toroidal fields, and RMP toroidal mode numbers. Consistent operational windows in edge safety factor are found across devices, while windows in plasma shaping parameters are distinct. Accessed pedestal parameters reveal a quantitatively similar pedestal-top density limit for RMP-ELM suppression in all devices of just over 3 × 10 19 m −3 . This is surprising given the wide variance of many engineering parameters and edge collisionalities, and poses a challenge to extrapolation of the regime. Wide ranges in input power, confinement time, and stored energy are observed, with the achieved triple product found to scale like the product of current, field, and radius. Observed energy confinement scaling with engineering parameters for RMP-ELM suppressed plasmas are presented and compared with expectations from established H and L-mode scalings, including treatment of uncertainty analysis. Different scaling exponents for individual engineering parameters are found as compared to the established scalings. However, extrapolation to next-step tokamaks ITER and SPARC find overall consistency within uncertainties with the established scalings, finding no obvious performance penalty when extrapolating from the assembled multi-device RMP-ELM suppressed database. Overall this work identifies common physics for RMP-ELM suppression and highlights the need to pursue this no-ELM regime at higher magnetic field and different plasma physical size.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Bioblendstocks to Optimize Mixing Controlled Compression Ignition (MCCI) Engines

In this project, a team of researchers from the University of Massachusetts Lowell, the University of Maine, and Mainstream Engineering developed an integrated process for the product of bioblendstocks to optimize mixing controlled compression ignition (MCCI) engines. The objective was to improve the energy density, sooting propensity, and cetane number of base diesel fuel while maintaining cold weather behavior. The process converts woody biomass (e.g. sawmill residues) into bio-oil through selective fast pyrolysis; the bio-oil is then selectively upgraded to form selectively oxygenated, minimally-branched hydrocarbons using non-noble metal catalysts in combination with metal-catalyzed hydrogenation. Advanced predictive models, in conjunction with existing property databases, and experimental testing are used to evaluate overall bioblendstock properties and their impact on base diesel fuel. An iterative, targeted upgrading approach was implemented to optimize the proposed bioblendstock’s properties. Assessment methodologies included techno-economic analysis, life-cycle assessment, property testing, and engine testing. Ultimately, the project team successfully produced a viable bioblendstock while identifying critical process points related to scale-up efforts. It was found that producing pyrolysis oils at 500 degrees C and with pine particle sizes of 1-2 mm led to bio-oil with a higher yield (of approximately 45 wt%) and rich amounts of aromatic alcohols. The resultant pyrolysis oil was then upgraded using a sequence of mild hydrotreating, followed by catalytic etherification and esterification, followed by another final mild hydrotreating to produce a blendstock containing saturated species with a limited, but non-zero, amount of oxygen. The aromatic alcohols produced by pyrolysis were especially helpful in this regard, as the resulting bicycloethers and derivatives exhibited high cetane numbers. While most bulk properties of the bioblendstock met or exceeded targeted thresholds, viscosity and cloud point notably fell outside the expected range; this could be addressed by blending limits and/or through the use of additives that are commonplace in current refinding practices. Identification of a bioblendstock that can be produced economically at scale while improving the performance and emissions characteristics of internal combustion engines positively affects the economy by boosting domestic fuel production and the environment by decreasing harmful emissions and increasing efficiency.

09 BIOMASS FUELS↗
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 NTRS2005

International Space Station Crew Restraint Design

With permanent human presence onboard the International Space Station (ISS), crews will be living and working in microgravity, dealing with the challenges of a weightless environment. In addition, the confined nature of the spacecraft environment results in ergonomic challenges such as limited visibility and access to the activity areas, as well as prolonged periods of unnatural postures. Without optimum restraints, crewmembers may be handicapped for performing some of the on-orbit tasks. Currently, many of the tasks on ISS are performed with the crew restrained merely by hooking their arms or toes around handrails to steady themselves. This is adequate for some tasks, but not all. There have been some reports of discomfort/calluses on the top of the toes. In addition, this type of restraint is simply insufficient for tasks that require a large degree of stability. Glovebox design is a good example of a confined workstation concept requiring stability for successful use. They are widely used in industry, university, and government laboratories, as well as in the space environment, and are known to cause postural limitations and visual restrictions. Although there are numerous guidelines pertaining to ventilation, seals, and glove attachment, most of the data have been gathered in a 1-g environment, or are from studies that were conducted prior to the early 1980 s. Little is known about how best to restrain a crewmember using a glovebox in microgravity. In 2004, The Usability Testing and Analysis Facility (UTAF) at the NASA Johnson Space Center completed development/evaluation of several design concepts for crew restraints to meet the various needs outlined above. Restraints were designed for general purpose use, for teleoperation (Robonaut) and for use with the Life Sciences Glovebox. All design efforts followed a human factors engineering design lifecycle, beginning with identification of requirements followed by an iterative prototype/test cycle. Anthropometric modeling was completed for the 5th percentile Asian female and the 95th percentile American male for all restraints. A series of three evaluations was performed onboard NASA's reduced gravity aircraft (KC-135). For all evaluations, participants performed representative tasks while being videotaped, and then completed a questionnaire following each flight day. The questionnaire included ratings scales and free format questions to assess topics such as comfort, stability provided, flexibility provided, etc. Results from the three flight evaluations are being used to develop the human factors design requirements for crew restraint concepts for 1) general purpose restraints, 2) teleoperation restraints and 3) glovebox restraints. The poster presentation will describe the detailed methodology used, results from each of the three evaluations, and the resulting human factors recommendations for the design of these restraints.

Whitmore, M.↗