Search NASA⌕ Search

SEARCH · Search NASA

Results for “COMPUTER MODELS”

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,189 records · Page 66

Maximized Information Gain of Next Generation Pulsed Power Using Optimized Design of Z-Machine Experiments

This project develops a Bayesian optimization approach to extracting insights from Z Machine experimental data to determine if and how these insights can be used to extrapolate to a larger facility. The primary goal is to address the scientific challenge of informing how confidently experimental conditions can be predicted on a next generation facility, the design of which requires the reliable extrapolation of current high energy density technologies to regimes yet unobserved, except by costly high-fidelity computational models. Maximizing the use of presently available data and understanding how it informs future endeavors is critically important to enable transformative pulsed power and the science of extreme conditions. We explore a Bayesian optimization approach to experimental design which combines information theory, experimental data, and computational modeling to explore how information gain can be maximized.

97 MATHEMATICS AND COMPUTING↗

The spiral-compact galaxy pair AM 2208-251: Computer simulations versus observations

The system AM2208-251 is a roughly edge-on spiral extending east-west with a smaller round compact E system about 60 arcsec east of the spiral nucleus along the major axis of the spiral. Bertola, Huchtmeier, and Zeilinger (1990) have presented optical spectroscopic as well as single dish 21 cm observations of this system. Their spectroscopic data show, via emission lines lambda lambda 3727-29A, a rising rotation curve near the nucleus. These spectroscopic observations may indicate a tidal interaction in the system. In order to learn more about such pairs, the authors simulated the interaction using the computer model developed by Miller (1976 a,b, 1978) and modified by the authors (Byrd 1986, 1987, 1988). To do the simulation they need an idea of the mutual orbits of the two galaxies. Their computer model is a two-dimensional polar N-body program. It consists of a self-gravitating disk of particles, within an inert axially symmetric stabilizing halo potential. The particles are distributed in a 24(radial) by 36(azimuthal) polar grid. Self consistent calculations can be done only within the grid area. The disk is modeled with a finite Mestel disk, where all the particles initially move in circular orbits with constant tangential velocities (Mestel 1963), resulting in a flat rotation curve. The gas particles in the spiral's disk, which make up 30 percent of its mass, collide in the following manner. The number of particles in each bin of the polar grid is counted every time step. If it is greater than a given critical density, all the particles in the bin collide, obtaining in the result the same velocities, equal to the average for the bin. This process produces clumps of gas particles-the star formation sites. The authors suppress the collision in the inner part of the disk (within the circle r = 6) to represent the hole seen in the gas in the nuclear bulge of spirals. They thus avoid spurious effects due to collisions in that region. They also varied the size of the collisional bins, which did not affect their conclusions.

Klaric, Mario↗

Metabolic Adaptation to Muscle Ischemia

Although all tissues in the body can adapt to varying physiological/pathological conditions, muscle is the most adaptable. To understand the significance of cellular events and their role in controlling metabolic adaptations in complex physiological systems, it is necessary to link cellular and system levels by means of mechanistic computational models. The main objective of this work is to improve understanding of the regulation of energy metabolism during skeletal/cardiac muscle ischemia by combining in vivo experiments and quantitative models of metabolism. Our main focus is to investigate factors affecting lactate metabolism (e.g., NADH/NAD) and the inter-regulation between carbohydrate and fatty acid metabolism during a reduction in regional blood flow. A mechanistic mathematical model of energy metabolism has been developed to link cellular metabolic processes and their control mechanisms to tissue (skeletal muscle) and organ (heart) physiological responses. We applied this model to simulate the relationship between tissue oxygenation, redox state, and lactate metabolism in skeletal muscle. The model was validated using human data from published occlusion studies. Currently, we are investigating the difference in the responses to sudden vs. gradual onset ischemia in swine by combining in vivo experimental studies with computational models of myocardial energy metabolism during normal and ischemic conditions.

Cabrera, Marco E.↗

Calculating Dose Release Factors to Obtain the Worst-Case Maximally Exposed Location by Sector

The aim of this project was to update the SRS Dose Release Factors (DRF's) for the Maximally Exposed Individual (MEI) for 20 area release points onsite. To update, potential MEIs now include onsite locations and modeled stack heights were updated to match onsite operations. The DRF's and MEI locations are used to ensure DOE and EPA compliance is met for dose standards. The Environmental Protection Agency (EPA) requires the use of the computer model CAP88 to estimate doses for demonstrating compliance with 40 CFR 61, the National Emission Standard for Hazardous Air Pollutants (NESHAP). The location of the maximally exposed individual (MEI) by sector and distance is required for execution of CAP88. The worst-case sector is the sector and distance where an offsite individual would receive the highest dose due to a potential release from a designated release location at a facility. The worst-case sector analysis addresses potential ground-level and stack-height releases from 20 areas of Savannah River Site (SRS). Weather patterns were analyzed at each of these areas relative to sixteen compass points. The analysis makes use of recently updated distance data and meteorological data. The computer model, CAP88 is EPA's software system used by the Department of Energy (DOE) to calculate the receptor dose due to routine atmospheric releases via a gaussian plume at the Savannah River Site (SRS). The dose-release factors (DRFs) are used to estimate the offsite and onsite dose from facility-specific releases and represent the dose to a receptor exposed to 1 Ci of a specified radionuclide being released into the atmosphere. The DRFs are simply a ratio of the TED to the release rate and, therefore, are directly scalable. The DRF values listed in the Appendix can be applied to expected release values from each area to estimate the potential dose to an individual adult at each location. NESHAP DRFs have been supplied to Environmental Compliance and Area Closure Projects (EC and ACP) for many years and are periodically updated to include changes in Cap88. estimate the potential dose to an individual adult at each location DRF in mrem Ci{sup -1} are estimated by modeling the total effective dose (TED) in mrem yr{sup -1} assuming an annual atmospheric release rate of 1 Ci yr-1 of the associated radionuclides. TED = Effective Dose (ED) + Committed Effective Dose (CED). TED = DRF/RR. Approach and Methods: Determine possible MEI locations Onsite/Offsite including buildings and farms; Determine the distances from the release area to the possible MEI location and the direction from the release to the MEI location; Using the most recent meteorological data (2007-2011) and parameters, run CAP88 for each release area for each compass sector and previously calculated distances to MEI locations using a 1 Ci release rate for tritium; Use the calculated the air concentration values to determine the MEI for each sector from CAP88. The relationship between air concentration and received dose is directly proportional and can be used to determine which locations to calculate DRF's and a dose for an MEI; The DRFs are found for each compass sector from 20 SRS release areas using CAP88-PC Version 4.0; Each location is calculated with ground level releases evaluated for each area for a list of radionuclides, some with and some without progenies; Compare the DRFs of the potential onsite MEIs to offsite MEIs. Based on the 2018 MEI Report, the Site worst case sector MEI is located 593 mW of B-Area. This MEI location represents a point along the main road that travels through a tract of land leased to the SC National Guard. While this location is estimated to have the highest exposure, minimal radioactivity is released from B-Area. In 2017, SRS requested and was approved by the EPA to change the location for modeling site releases from the Center of Site to H-Area, and to use six stack heights. H-Area emission account for approximately 90% of all radioactive operations onsite. The MEI for H-Area was found to be 11,458 m NNE which locates the Savannah River Ecology Lab Conference Center. Comparisons for 2018 data and 2014 shows the new methods of calculation are more conservative yet still remain well below the DOE 100 mrem limit.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

A Scheduling Algorithm for Replicated Real-Time Tasks

We present an algorithm for scheduling real-time periodic tasks on a multiprocessor system under fault-tolerant requirement. Our approach incorporates both the redundancy and masking technique and the imprecise computation model. Since the tasks in hard real-time systems have stringent timing constraints, the redundancy and masking technique are more appropriate than the rollback techniques which usually require extra time for error recovery. The imprecise computation model provides flexible functionality by trading off the quality of the result produced by a task with the amount of processing time required to produce it. It therefore permits the performance of a real-time system to degrade gracefully. We evaluate the algorithm by stochastic analysis and Monte Carlo simulations. The results show that the algorithm is resilient under hardware failures.

Yu, Albert C.↗

Determination of Ballistic Limit for IM7/8552 Using LS-DYNA MAT 261

The goal of the NASA ACC High Energy Dynamic Impact Project is to determine the state of the art of dynamic fracture simulations for high velocity impact for composite fuselage shielding applications. Using a building block approach, several computational models considered under NASA ACC are being validated against test data, starting at unconfigured panels and progressing to configured panels under combined out-of-plane and in-plane loading due to ballistic impact. The computational models being evaluated in this project include MAT 162, MAT 213, MAT 261, SPG, and Peridynamics. In this paper, the simulation results using LS-Dyna Material MAT 261 are presented. In particular, a series of blind predictions for unconfigured panels were performed to determine the ballistic limit or V50 velocity. MAT 261 employs failure approach that is generally physically-based using fracture toughness criteria. The overall material model relies on typical ply-level stiffness properties, similar to MAT 162 and other composite continuum damage material models. The fracture toughness values are based on standard tests, and thus are not subject to extensive calibration. This approach is more efficient than performing extensive optimization studies for calibration of parameters. Also, this approach of relying on physical properties reduces the uncertainty of results, as questions concerning the quality and extent of the calibration studies is no longer relevant. However, it was found that carefully controlled coupon-level tests are needed to accurately obtain the required fracture toughness values. Additionally, it should be noted that there is one significant parameter in MAT 261 that does appear to require calibration, and that is the overall failure strain. This is the strain at which the element is deleted, and is not the same as the strain at which damage begins to accumulate. This failure strain is termed EFS (Effective Failure Strain), and is the maximum effective stain for element failure. Simulations have shown that this value will significantly affect impact response and failure. The paper presents the effect of this element failure strain parameter, along with possible uncertainties in fracture toughness values. With an adjusted appropriate value for EFS, it is seen that simulation results compare well with impact test data for predicted penetration velocity.

Byar, Alan↗

UQ Toolkit v. 3.0

The Uncertainty Quantification (UQ) Toolkit is a software library for the characterization and propagation of uncertainties in computational models. This library provides Bayesian inference tools to infer uncertain model parameters, Bayesian compressive sensing methods for discovering sparse representations of high-dimensional input-output response surfaces, methods for constructing Karhunen-Loève representations of stochastic processes, and global sensitivity analysis tools used to compute Sobol indices in order to characterize the importance of uncertain inputs or the interactions between them. The basis for many of these methods relies on representing random variables with Polynomial Chaos expansions (PCEs). This library implements several spectral basis function types (e.g. Hermite basis functions in terms of Gaussian random variables or Legendre basis functions in terms of uniform random variables) that can be used to represent random variables with PCEs. For the propagation of uncertainty, tools are provided to propagate PCEs that describe the input uncertainty through the computational model using either intrusive methods (Galerkin projection of equations onto basis functions) or non-intrusive methods (perform deterministic operation at sampled values of the random values and project the obtained results onto basis functions).

Sargsyan, Khachik↗

In-Depth Modeling and Simulation Analysis of Artemis Missions Using the Impact Probabilistic Risk Assessment Tool

BACKGROUND The Artemis campaign is a Moon exploration program with a series of six planned missions, five of which will be crewed. These five crewed missions will contain a single mission segment (space flight), or multiple mission segments involving space flight (Orion), lunar landing (LTV) and/or space habitat (Gateway). Each crewed segment faces the risk of unique medical conditions, necessitating medical sets/kits tailored to those specificities. To support and enable a data-driven and evidence-based decision-making process through out a mission’s life cycle, a software tool called IMPACT was developed. Using probabilistic risk assessment (PRA) methodologies, IMPACT (Informing Mission Planning via Analysis of Complex Tradespaces) is a novel tool built for analyzing the possibility of encountering complex medical risks during space flight, and for identifying the medical resources and capabilities needed to treat those potential at-risk medical conditions. IMPACT achieves this by performing hundreds of thousands of Monte Carlo simulations of missions to build aggregate pictures of medical risk. During an extended simulation modeling phase, IMPACT generated analytical results for medical risks, and the medical resources and capabilities to address those risks, for every segment of every crewed Artemis mission. This presentation will highlight the reliability, consistency and validity of IMPACT’s computational modeling techniques and will showcase the library of analytical outcomes generated for the Artemis missions. OVERVIEW During the early stages of IMPACT’s design, architecture and technical requirements collection, “scenarios” (use cases) - achievement goals required for acceptance testing, were identified by stakeholders. IMPACT successfully completed the scenario testing requirements and undertook an extensive operational run phase utilizing a wide range of input combinations with a goal of delivering a cohesive, trustworthy, reliable, vast, and diverse body of evidence. The intent of these modeling runs was to validate consistency in output, ensure solidity of executable operations and to streamline processes by identifying areas requiring efficiency improvements. Using the many missions of Artemis, IMPACT ran variations of operational runs to assess the output for acceptable, as well as unusual characteristics. This rigorous long-term “shakedown” analysis was implemented to help build a collective body of evidence to aid in securing a high level of confidence, reliability, and validity in the output, whether from the applicational components of IMPACT, or the entirety of the operational process. ANTICIPATED ANALYSIS AND CONCLUSION This presentation will discuss the various categories of input criteria; the comparisons in the application of these input criteria to various Artemis missions; the preparation and collection of the body of evidence, and reliability of the computational modeling techniques. This paper serves as an initial analytical overview of IMPACT’s probabilistic risk assessment (PRA) medical risk outputs covering Artemis missions and is not intended to be deemed the official medical response for the Artemis campaign.

Crew Composition↗

A hierarchal model for bacterial cell inactivation in solution by direct and indirect treatment using cold atmospheric plasmas

Cold atmospheric plasma devices have shown promise for a variety of plasma medical applications, including wound healing and bacterial inactivation often performed in liquids. In the latter application, plasma-produced reactive oxygen and nitrogen species (RONS) interact with and damage bacterial cells, though the exact mechanism by which cell damage occurs is unclear. Computational models can help elucidate relationships between plasma-produced RONS and cell killing by enabling direct comparison between dissimilar plasma devices and by examining the effects of changing operating parameters in these devices. In biological applications, computational models of plasma-liquid interactions would be most effective in design and optimization of plasma devices if there is a corresponding prediction of the biological outcome. In this work, we propose a hierarchal model for planktonic bacterial cell inactivation by plasma produced RONS in liquid. A previously developed reaction mechanism for plasma induced modification of cysteine was extended to provide a basis for cell killing by plasma-produced RONS. Results from the model are compared to literature values to provide proof of concept. Differences in time to bacterial inactivation as a function of plasma operating parameters including gas composition and plasma source configuration are discussed. Results indicate that optimizing gas-phase reactive nitrogen species production may be key in the design of plasma devices for disinfection.

59 BASIC BIOLOGICAL SCIENCES↗

Predicting Transport of 3,5,6-Trichloro-2-Pyridinol Into Saliva Using a Combination Experimental and Computational Approach

A combination experimental and computational approach was developed to predict chemical transport into saliva. A serous-acinar chemical transport assay was established to measure chemical transport with non-physiological (standard cell culture medium) and physiological (using surrogate plasma and saliva medium) conditions using 3,5,6-trichloro-2-pyridinol (TCPy) a metabolite of the pesticide chlorpyrifos. High levels of TCPy protein binding was observed in cell culture medium and rat plasma resulting in different TCPy transport behaviors in the two experimental conditions. In the non-physiological transport experiment, TCPy reached equilibrium at equivalent concentrations in apical and basolateral chambers. At higher TCPy doses, increased unbound TCPy was observed, and TCPy concentrations in apical and basolateral chambers reached equilibrium faster than lower doses, suggesting only unbound TCPy is able to cross the cellular monolayer. In the physiological experiment, TCPy transport was slower than non-physiological conditions, and equilibrium was achieved at different concentrations in apical and basolateral chambers at a comparable ratio (0.034) to what was previously measured in rats dosed with TCPy (saliva:blood ratio: 0.049). A cellular transport computational model was developed based on TCPy protein binding kinetics and accurately simulated all transport experiments using different permeability coefficients for the two experimental conditions (1.4 vs 0.4 cm/hr for non-physiological and physiological experiments, respectively). The computational model was integrated into a physiologically based pharmacokinetic (PBPK) model and accurately predicted TCPy concentrations in saliva of rats dosed with TCPy. Overall, this study demonstrates an approach to predict chemical transport in saliva potentially increasing the utility of salivary biomonitoring in the future.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Computational study of tungsten and depleted uranium photoneutron targets for a 20 MeV electron linear accelerator

Neutron production can be realized with a high energy electron linear accelerator by using Bremsstrahlung and photoneutron converters. In this study, Monte Carlo N-Particle Code (MCNP) was used to evaluate potential photonuclear target designs for a high energy electron linear accelerator for applications such as neutron radiography and neutron resonance spectroscopy. A computational model was developed to inform a target design that would yield a high number of neutrons. It consists of a 20 MeV electron beam incident on a Bremsstrahlung target and a photonuclear target to generate neutrons. This computational model showed that a thickness of 0.75 inches for both tungsten and depleted uranium yields the most neutrons from photoneutron reactions. Saturation in the total number of generated neutrons was observed at over 0.75-inch thickness for both evaluated materials. Depleted uranium yielded approximately twice the number of neutrons overall compared to tungsten. The highest neutron surface flux for Depleted Uranium was 1.06 × 10-4 neutrons/cm2/source electron, and for Tungsten it was 5.12 × 10-5 neutrons/cm2/source electron. The optimal target design for this study’s application would consist of a 0.75 inch-thick block of depleted uranium with the length, width, and/or diameter varying dependent on application.

43 PARTICLE ACCELERATORS↗

A visco-plastic constitutive model for accurate densification and shape predictions in powder metallurgy hot isostatic pressing

Powder metallurgy hot isostatic pressing (PM-HIP) is an advanced manufacturing process that produces near net shape parts with high material utilization and uniform microstructures. Despite being used frequently to produce small-scale components, the application of PM-HIP to large-scale components is limited due to inadequate understanding of its complex mechanisms that cause unpredictable post-HIP shape distortions. A computational model can provide necessary information about the intermediate and final stages of the HIP process that can help understand it better and make accurate predictions. Generally, two types of computational models are employed for PM-HIP of metal powders, namely, plastic and visco-plastic models. Between these, the plastic model is preferred due to its cheaper calibration approach requiring less experimental data. However, the plastic model sometimes produces incorrect predictions when slight variations of the HIP conditions are encountered in practical situations. Therefore, this work presents a visco-plastic model that addresses these limitations of the plastic model. A novel modified calibration approach is employed for the visco-plastic model that utilizes less experimental data than existing approaches. With the new approach, the data requirement is same for both plastic and visco-plastic models. This also enables a quantitative comparison of plastic and visco-plastic models, which have been only qualitatively compared in the past. When calibrated with the same experimental data, both the models are found to produce similar results. In conclusion, the calibrated visco-plastic model is applied to several complex geometries, and the predictions are found to be in good agreement with experimental observations.

Hot isostatic pressing↗

Metal Big Area Additive Manufacturing: Process Modeling and Validation

Metal Big Area Additive Manufacturing (mBAAM) is a new additive manufacturing (AM) technology for printing large-scale 3D objects. mBAAM is based on the gas metal arc welding process and uses a continuous feed of welding wire to manufacture an object. An electric arc forms between the wire and the substrate, which melts the wire and deposits a bead of molten metal along the predetermined path. In general, the welding process parameters and local conditions determine the shape of the deposited bead. The sequence of the bead deposition and the corresponding thermal history of the manufactured object determine the long range effects, such as thermal-induced distortions and residual stresses. Therefore, the resulting performance or final properties of the manufactured object are dependent on its geometry and the deposition path, in addition to depending on the basic welding process parameters. Physical testing is critical for gaining the necessary knowledge for quality prints, but traversing the process parameter space in order to develop an optimized build strategy for each new design is impractical by pure experimental means. Computational modeling and optimization may accelerate development of a build process strategy and saves time and resources. Because computational modeling provides these opportunities, we have developed a physics-based Finite Element Method (FEM) simulation framework and numerical models to support the mBAAM process s development and design. In this paper, we performed a sequentially coupled heat transfer and stress analysis for predicting the final deformation of a small rectangular structure printed using the mild steel welding wire. Using the new simulation technologies, material was progressively added into the FEM simulation as the arc weld traversed the build path. In the sequentially coupled heat transfer and stress analysis, the heat transfer was performed to calculate the temperature evolution, which was used in a stress analysis to evaluate the residual stresses and distortions. In this formulation, we assume that physics is directionally coupled, i.e. the effect of stress of the component on the temperatures is negligible. The experiment instrumentation (measurement types, sensor types, sensor locations, sensor placements, measurement intervals) and the measurements are presented. The temperatures and distortions from the simulations show good correlation with experimental measurements. Ongoing modeling work is also briefly discussed.

Additive Manufacturing↗

Radiolytic degradation of formic acid and formate in aqueous solution: modeling the final stages of organic mineralization under advanced oxidation process conditions

The successful use of advanced oxidation processes to treat aqueous solutions containing undesirable organic species requires the degradation of these species to lower molecular weight, lower hazard compounds. Safe application of this technology requires a thorough understanding of the mechanisms of degradation. These oxidative transformations are mainly initiated by the reactions of reactive oxygen species, particularly hydroxyl radicals. These react with organic molecules to generate carbon-centered radicals. In the presence of dissolved oxygen, the carbon-centered radicals are next converted to peroxyl radicals, which then decay to lower molecular weight species by multiple mechanistic pathways. Formic acid and its conjugate base formate are the last stable chemical species produced immediately before the complete mineralization of any organic molecule undergoing oxidative degradation in aqueous solution. Once understood, the radical-induced chemistry of formic acid/formate under these conditions has wide applicability in all advanced oxidation technologies. To develop this quantitative knowledge, here we have performed a series of 60 Co gamma irradiation studies on aqueous formic acid/formate over different pH and solution conditions. The measured species concentration changes, as a function of applied dose, are compared with the predictions of a kinetic computer model constructed from literature reactions and reported rate coefficients. The excellent agreement found between the results and modeling gives confidence in the mechanism presented here and provide the first complete computer model for the radiolytic degradation of formic acid in water.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

A computational study of multiple jet and wall interaction

In the fuel preburner of the Space Shuttle Main Engine, face plate, injector, and baffle erosion have been observed. The observed patterns of erosion suggest that flame attachment to the walls is a contributing factor. To better understand the physical phenomena involved, a portion of the preburner was modeled computationally. The simulated 'preburner' had three two-dimensional jets entering a cavity adjacent to a baffle. The computational model employed the Patankar Spalding algorithm with upwind differencing. The turbulence model was a standard k-epsilon model with wall functions. The effects of incoming boundary conditions on turbulent kinetic energy and dissipation, k and epsilon, was studied. The results indicate a very strong sensitivity to these boundary conditions over certain ranges of values.

Campbell, W.↗

Probing the Origin of Interfacial Carriers in SrTiO 3 –LaCrO 3 Superlattices

Emergent phenomena at complex oxide interfaces could provide the basis for a wide variety of next-generation devices, including photovoltaics and spintronics. To date, detailed characterization and computational modeling of interfacial defects, cation intermixing, and film stoichiometry have helped to explain many of the novel behaviors observed at a single heterojunction. Unfortunately, many of the techniques employed to characterize a single heterojunction are less effective for a superlattice made up of a repeating series of interfaces that induce collective interfacial phenomena throughout a film. However, these repeating interfaces present an untapped opportunity to introduce an additional degree of complexity, such as confined electric fields, that cannot be studied in a single heterojunction. In this work, we explore the properties of SrTiO3-LaCrO3 superlattices to understand the role of defects, including variations in cation stoichiometry of individual layers of the superlattice, intermixing across interfaces, and interfacial oxygen vacancies. Using x-ray photoelectron spectroscopy (XPS) and scanning transmission electron microscopy electron energy-loss spectroscopy (STEM-EELS), we quantify the stoichiometry of individual layers of the superlattice and determine the degree of intermixing in these materials. By comparing these results to both density functional theory (DFT) models and STEM-EELS measurements of the Ti and Cr valence in each layer of the superlattice, we correlate different types of defects with the associated materials properties of the superlattice. We show that a combination of ab initio modeling and complementary structural characterization methods can offer unique insight into structure-property relationships in many oxide superlattice systems.

36 MATERIALS SCIENCE↗

Phosphoric acid fuel cell power plant system performance model and computer program

A FORTRAN computer program was developed for analyzing the performance of phosphoric acid fuel cell power plant systems. Energy mass and electrochemical analysis in the reformer, the shaft converters, the heat exchangers, and the fuel cell stack were combined to develop a mathematical model for the power plant for both atmospheric and pressurized conditions, and for several commercial fuels.

Alkasab, K. A.↗

A computationally efficient modelling of laminar separation bubbles

The goal of this research is to accurately predict the characteristics of the laminar separation bubble and its effects on airfoil performance. To this end, a model of the bubble is under development and will be incorporated in the analysis section of the Eppler and Somers program. As a first step in this direction, an existing bubble model was inserted into the program. It was decided to address the problem of the short bubble before attempting the prediction of the long bubble. In the second place, an integral boundary-layer method is believed more desirable than a finite difference approach. While these two methods achieve similar prediction accuracy, finite-difference methods tend to involve significantly longer computer run times than the integral methods. Finally, as the boundary-layer analysis in the Eppler and Somers program employs the momentum and kinetic energy integral equations, a short-bubble model compatible with these equations is most preferable.

Maughmer, Mark D.↗