Search NASA⌕ Search

SEARCH · Search NASA

Results for “process simulation modeling”

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,117 records · Page 62
All sourcesNASA NTRSDOE OSTI
DOE OSTI2015

Voltage equilibration for reactive atomistic simulations of electrochemical processes

We introduce electrochemical dynamics with implicit degrees of freedom (EChemDID), a model to describe electrochemical driving force in reactive molecular dynamics simulations. The method describes the equilibration of external electrochemical potentials (voltage) within metallic structures and their effect on the self-consistent partial atomic charges used in reactive molecular dynamics. An additional variable assigned to each atom denotes the local potential in its vicinity and we use fictitious, but computationally convenient, dynamics to describe its equilibration within connected metallic structures on-the-fly during the molecular dynamics simulation. This local electrostatic potential is used to dynamically modify the atomic electronegativities used to compute partial atomic changes via charge equilibration. Validation tests show that the method provides an accurate description of the electric fields generated by the applied voltage and the driving force for electrochemical reactions. We demonstrate EChemDID via simulations of the operation of electrochemical metallization cells. The simulations predict the switching of the device between a high-resistance to a low-resistance state as a conductive metallic bridge is formed and resistive currents that can be compared with experimental measurements. In addition to applications in nanoelectronics, EChemDID could be useful to model electrochemical energy conversion devices.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2018

Model-Based Air Flow Path Optimization of Heat Pump Clothes Dryer

A heat pump clothes dryer (HPCD) uses a vapor compression system to dry clothes. The condenser heats air, which passes through the drum to evaporate moisture out of the clothes, and the evaporator condenses water out of the air stream. As a result, the HPCD can achieve 50% energy savings compared to a conventional electric resistance dryer. In this work we developed a physics-based, quasi-steady-state HPCD system model with detailed heat exchanger and compressor models. The system model can simulate the inherently transient HPCD drying process, to size components, and to reveal trends in key variables (e.g. compressor discharge temperature, power consumption, required drying time, etc.) The system model was calibrated using experimental data from a prototype HPCD. Air leakages, in and out, along the closed air circulation path of HPCD cause varied effects on the performance. Understanding the location, magnitude, and direction of air leakage of the heat pump clothes dryer is critical for accurately characterizing the performance and developing a high-performance design. The system model was used to reveal the impacts. In addition, model-based parametric optimizations were conducted to design the HPCD charge inventory, flow rate, air path and leak points for optimum performance.

Shen, Bo↗
DOE OSTI2026

Hybrid combustion modeling approach for turbulent jet ignition in natural-gas pre-chamber spark-ignition engines at high EGR

Here, this study presented a hybrid modeling approach for simulating turbulent jet ignition and combustion processes in a natural-gas pre-chamber spark-ignition engine operating under exhaust gas recirculation (EGR) diluted conditions. In-depth analyses of experimental data and simulation results from previous work [Chinnathambi et al., ICEF2021-67836; Kim et al., Fuel 409: 137815, 2026] revealed two key findings: (i) the magnitude of pressure difference between the pre-chamber and main chamber ($∆P_{PC-MC}$) was positively correlated with the combustion duration from the moment of $∆P_{PC-MC}=0$ to the point of 5% mass fraction burned, with larger $∆P_{PC-MC}$ associated with longer duration; and (ii) the turbulent combustion regime in the main chamber transitioned from the broken reaction zone to the corrugated flamelet regime, with the Karlovitz number exceeding 100 immediately after turbulent hot jets were ejected from nozzles, coinciding with observed local extinction events. To accurately simulate the entire combustion process, a hybrid approach was developed under Reynolds-Averaged Navier Stokes framework, combining the G-equation model for pre-chamber combustion with the multi-zone well-stirred reactor approach and a turbulence-chemistry interaction (TCI) submodel for main chamber combustion. The TCI submodel accounted for the attenuation of reaction rates due to turbulent strain and modeled local extinction by suppressing reaction rates under certain flow and flame conditions. When applied to three EGR rate conditions toward the dilution limit, the hybrid modeling approach accurately reproduced experimental data in terms of cylinder pressure, apparent heat release rate, and the observed positive correlation, including the delayed onset of main chamber combustion—a feature not captured by existing combustion models.

computational fluid dynamics simulation↗
NASA NTRS1988

Numerical simulation of tropical cumulus clouds and their interaction with the subcloud layer

A two-dimensional numerical model, with an emphasis on turbulent processes in the boundary layer and in clouds, was developed for simulating an ensemble of cumulus clouds. The model was used to investigate the response of tropical cumulus clouds to imposed large-scale vertical advection by performing two simulations which were identical except that one was with upward large-scale vertical velocity (the 'disturbed' case) and the other was without any advection (the 'undisturbed' case). The results showed that deep cumulus clouds formed in the disturbed case, while only shallow clouds formed in the undisturbed case. The subcloud layer (SCL) in the disturbed case was strongly affected by cumulus circulations and rain evaporation, with both cumulus updrafts and downdrafts being important contributors to the SCL sensible and latent heat fluxes. Cumulus-scale circulation in the disturbed SCL determined where new clouds formed by creating convergence zones and moisture anomalies.

Krueger, Steven K.↗
DOE OSTI2019

Examination of Dendritic Growth During Solidification of Ternary Alloys via a Novel Quantitative 3D Cellular Automaton Model

A three-dimensional (3D) quantitative cellular automaton (CA) model was developed to simulate dendritic growth during solidification processing of ternary alloys. A detailed method was proposed to solve solute diffusion and calculate the solid fraction for ternary alloys during solidification. The present model has been shown to accurately predict dendrite morphologies and solute distributions of both single equiaxed dendrite and columnar dendrites during solidification. The model was also used to study the influence of the concentration of a third component (Mg in the Al-Si-Mg system) on dendritic growth. With increasing Mg concentration, the steady-state dendrite tip growth velocity was shown to decrease, resulting in a shorter primary dendrite length. The 3D CA simulation results agree well with the prediction of the LGK theoretical model. Multi-columnar dendrite growth was simulated with different cooling rates. Primary and secondary dendrite morphology was shown to have little variation at low cooling rates (1.0 to 2.0 K/s) and a constant undercooling. The average secondary dendrite arm spacing was shown to decrease with the increase of cooling rate within this range. The 3D CA simulation results are in good agreement with the experimental directional solidification experiments of a commercial A356 (Al-7 wt pct Si-0.3 wt pct Mg).

36 MATERIALS SCIENCE↗
DOE OSTI2023

Parameterization and Explicit Modeling of Cloud Microphysics: Approaches, Challenges, and Future Directions

Cloud microphysical processes occur at the smallest end of scales among cloud-related processes and thus must be parameterized not only in large-scale global circulation models (GCMs) but also in various higher-resolution limited-area models such as cloud-resolving models (CRMs) and large-eddy simulation (LES) models. Instead of giving a comprehensive review of existing microphysical parameterizations that have been developed over the years, this study concentrates purposely on several topics that we believe are understudied but hold great potential for further advancing bulk microphysics parameterizations: multi-moment bulk microphysics parameterizations and the role of the spectral shape of hydrometeor size distributions; discrete vs “continuous” representation of hydrometeor types; turbulence-microphysics interactions including turbulent entrainment-mixing processes and stochastic condensation; theoretical foundations for the mathematical expressions used to describe hydrometeor size distributions and hydrometeor morphology; and approaches for developing bulk microphysics parameterizations. Also presented are the spectral bin scheme and particle-based scheme (especially, super-droplet method) for representing explicit microphysics. Their advantages and disadvantages are elucidated for constructing cloud models with detailed microphysics that are essential to developing processes understanding and bulk microphysics parameterizations. Particle-resolved direct numerical simulation (DNS) models are described as an emerging technique to investigate turbulence-microphysics interactions at the most fundamental level by tracking individual particles and resolving the smallest turbulent eddies in turbulent clouds. Outstanding challenges and future research directions are explored as well.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS1984

Atmospheric Modeling And Sensor Simulation (AMASS) study

The capabilities of the atmospheric modeling and sensor simulation (AMASS) system were studied in order to enhance them. This system is used in processing atmospheric measurements which are utilized in the evaluation of sensor performance, conducting design-concept simulation studies, and also in the modeling of the physical and dynamical nature of atmospheric processes. The study tasks proposed in order to both enhance the AMASS system utilization and to integrate the AMASS system with other existing equipment to facilitate the analysis of data for modeling and image processing are enumerated. The following array processors were evaluated for anticipated effectiveness and/or improvements in throughput by attachment of the device to the P-e: (1) Floating Point Systems AP-120B; (2) Floating Point Systems 5000; (3) CSP, Inc. MAP-400; (4) Analogic AP500; (5) Numerix MARS-432; and (6) Star Technologies, Inc. ST-100.

Parker, K. G.↗
DOE OSTI2023

Simulation of radioactive plume transport in the atmosphere including dynamics of particle aggregation and breakup

Accurate prediction of the atmospheric transport of debris particles relies heavily on our knowledge of the size distribution of the particles within a debris cloud. Assuming a fixed particle size during simulations is not always viable since the size distribution of the debris can change during transport. Various microphysical processes, such as aggregation and breakup, influence debris particles and dictate any changes to the size distribution. To track those changes that can occur, a population balance model can be adopted and instituted within a model framework. Nonetheless, many of the models that simulate the transport of radioactive debris following a device-driven fission incident have historically neglected to consider these processes. As such, this work describes our effort to develop a modeling framework capable of simulating the transport and deposition of a radioactive plume generated from a fission incident with a dynamic population balance including particle aggregation and breakup. Herein the impact of aggregation and breakup, individually and collectively, on the particle size distribution is explored using the developed framework. When simulating aggregation, for example, six mechanisms, including Brownian coagulation, the convective enhancement to Brownian coagulation, van der Waals-viscous force correction for Brownian coagulation, gravitational collection, turbulent inertial motion, and turbulent shear, are considered. Brownian coagulation and its corrections have, as one would expect, a large impact on relatively small aggregates. Aggregates with a diameter that is less than or equal to 1.0 μm, for instance, comprise 50.6 vol % of all aggregates in the absence of aggregation and 31.2 vol % when Brownian coagulation and its corrections are considered. Gravitational collection and, to a much lesser extent, turbulent shear and turbulent inertial motion are, conversely, of great importance to relatively large aggregates (i.e., diameter greater than 3.0 μm). Additionally, the individual effects of atmospheric and particle parameters, such as wind speed and particle density, are examined. Of the parameters examined, turbulent energy dissipation and aggregate fractal dimension (i.e., aggregate shape with lower values representing more irregular particles) were of substantial importance since both terms directly impact aggregate stability and, by extension, the breakup rate. Large-scale transport and deposition simulations in a dry atmosphere are also presented and discussed as a proof of concept.

98 NUCLEAR DISARMAMENT, SAFEGUARDS, AND PHYSICAL P↗
DOE OSTI2022

Toward Better Understanding of Microphysical Processes and Resulting Precipitation Physics: A Merger Of Observations and Cloud Models (Final Report)

A large database of global disdrometer observations and a diverse set of simulations from the Regional Atmospheric Modeling System (RAMS) were compared using principal component analysis (PCA) in order to better understand warm and ice-based precipitation processes. The analysis demonstrated that six distinct precipitation groups (PGs) with common characteristics were revealed in both the observations and model simulations. These PGs were defined by similar co-variability of rain parameters. However, the model showed large concentrations of drops with a given size compared to the observations. Detailed investigations determined the parameterization of drop breakup was forcing the drops to an equilibrium size with a stronger constraint than was indicated by the observations. A series of sensitivity studies permutating the rain shape parameter in the RAMS 2-moment bulk scheme demonstrated how the assumed shape of the rain distribution influences the microphysics as well as precipitation characteristics. Results were compared to the same simulations performed with a bin microphysics scheme where the rain shape parameter is freely evolving. In the bin scheme, a wide range of shape parameters were found, with horizontal and vertical variability, in addition to being a function of storm lifetime. The results of this study highlight the limitations of a fixed assumed rain distribution in 2-moment microphysics schemes. Model simulations were used to probe the microphysical origins of the six PGs. Rain budgets from a supercell case showed that several of the groups hypothesized to be associated with strong ice microphysical processes had significant contributions from melting hail, but the complexity of ice and warm-rain processes in the supercell precluded definitive determination of microphysical origins for many of the PGs. A case study from the Mid-latitude Continental Convection and Clouds Experiment (MC3E) compared the model PGs to disdrometer observations, and demonstrated the spatial and temporal variability of the PGs which is not possible from disdrometer point measurements alone. The analyses performed during this project demonstrated how statistical analysis can provide a robust method for comparing model simulations and observations, which ultimately resulted identification of limitations in some current microphysics parameterizations as well as insights into precipitation variability which ultimately benefits both observations and modeling efforts.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

LandScan mosaic enables high-resolution gridded population estimates with explicit uncertainty

Gridded population datasets represent high-resolution distributions of human occupancy, enabling informed decision-making across a broad range of fields. These data products are valuable for assessing environmental risk, urban development, disaster preparedness and resource allocation—areas where accurate population estimates directly enhance policy effectiveness and optimize resource distribution. Despite the importance of gridded population datasets, traditional population modeling approaches often overlook inherent uncertainties in the estimation process. This limitation can create a false sense of certainty in population estimates, potentially leading to flawed decisions by those who rely on the data. To address this methodological gap, we introduce a probabilistic machine learning modeling framework, LandScan Mosaic, that explicitly incorporates uncertainty into the population modeling process. Our approach systematically quantifies uncertainty in three key modeling parameters of the LandScan HD gridded population dataset: building use types, floor counts, and occupancy rates. By employing Monte Carlo simulations, we propagate these uncertainties through the modeling process, yielding probability distributions of population counts in place of deterministic point estimates. We demonstrate the practical application of this framework in Iloilo City, Philippines, using structured decision-making techniques and our probabilistic estimates to identify and prioritize areas most affected by projected flooding, supporting targeted interventions that address both economic and social risks. In doing so, we propose a population-specific approach for incorporating confidence into structured decision making processes. Through a comparative analysis with conventional deterministic approaches and point estimate approaches, including LandScan HD and WorldPop, we evaluate how the incorporation of machine learning and uncertainty influences decision rankings. This research advances population distribution modeling by offering a robust, quantitative approach that explicitly accounts for uncertainty in the underlying data, along with guidance for how users can apply uncertainty in their decision-making.

Environmental sciences↗
DOE OSTI2018

Near Surface Nucleation and Particle Mediated Growth of Colloidal Au Nanocrystals

During non-classical growth of nanostructures via assembly of primary nuclei, nucleation and assembly are assumed to be distinct processes: nanoparticles nucleate randomly and aggregate to form extended structures through Brownian motion in the presence of long-range attractive interactions. Here we investigate the relationship between these two processes by using in situ AFM, in situ, ex situ and cryo TEM and UV-Vis spectroscopy to observe growth of colloidal gold and simulations to develop a mechanistic model of the process. Our results reveal an inexorable link between nucleation and assembly with nuclei forming almost exclusively within a ~ 1 nm interfacial region of existing particles. The new particles immediately close the gap either through a diffusive jump or via growth of a neck between the seed and new particle, generating aggregates exhibiting features commonly attributed to oriented attachment of independently nucleated particles. The results demonstrate that creation of initial particle interfaces leads to local environments that redirect growth towards non-classical processes.

Cheng, Yingwen↗
NASA NTRS1988

Appendices to the model description document for a computer program for the emulation/simulation of a space station environmental control and life support system

A Model Description Document for the Emulation Simulation Computer Model was already published. The model consisted of a detailed model (emulation) of a SAWD CO2 removal subsystem which operated with much less detailed (simulation) models of a cabin, crew, and condensing and sensible heat exchangers. The purpose was to explore the utility of such an emulation simulation combination in the design, development, and test of a piece of ARS hardware, SAWD. Extensions to this original effort are presented. The first extension is an update of the model to reflect changes in the SAWD control logic which resulted from test. Also, slight changes were also made to the SAWD model to permit restarting and to improve the iteration technique. The second extension is the development of simulation models for more pieces of air and water processing equipment. Models are presented for: EDC, Molecular Sieve, Bosch, Sabatier, a new condensing heat exchanger, SPE, SFWES, Catalytic Oxidizer, and multifiltration. The third extension is to create two system simulations using these models. The first system presented consists of one air and one water processing system. The second consists of a potential air revitalization system.

Yanosy, James L.↗
NASA NTRS1988

Appendices to the user's manual for a computer program for the emulation/simulation of a space station environmental control and life support system

A user's Manual for the Emulation Simulation Computer Model was published previously. The model consisted of a detailed model (emulation) of a SAWD CO2 removal subsystem which operated with much less detailed (simulation) models of a cabin, crew, and condensing and sensible heat exchangers. The purpose was to explore the utility of such an emulation/simulation combination in the design, development, and test of a piece of ARS hardware - SAWD. Extensions to this original effort are presented. The first extension is an update of the model to reflect changes in the SAWD control logic which resulted from the test. In addition, slight changes were also made to the SAWD model to permit restarting and to improve the iteration technique. The second extension is the development of simulation models for more pieces of air and water processing equipment. Models are presented for: EDC, Molecular Sieve, Bosch, Sabatier, a new condensing heat exchanger, SPE, SFWES, Catalytic Oxidizer, and multifiltration. The third extension is to create two system simulations using these models. The first system presented consists of one air and one water processing system, the second a potential Space Station air revitalization system.

Yanosy, James L.↗
NASA NTRS2008

A Lunar Surface Operations Simulator

The Lunar Surface Operations Simulator (LSOS) is being developed to support planning and design of space missions to return astronauts to the moon. Vehicles, habitats, dynamic and physical processes and related environment systems are modeled and simulated in LSOS to assist in the visualization and design optimization of systems for lunar surface operations. A parametric analysis tool and a data browser were also implemented to provide an intuitive interface to run multiple simulations and review their results. The simulator and parametric analysis capability are described in this paper.

Nayar, H.↗
DOE OSTI2021

Defense Waste Processing Facility Nitric-Glycolic Flowsheet Chemical Process Cell Chemistry: Part 2

The conversions of nitrite to nitrate, the destruction of glycolate, and the conversion of glycolate to formate and oxalate were modeled for the Nitric-Glycolic flowsheet using data from Chemical Process Cell (CPC) simulant runs conducted by Savannah River National Laboratory (SRNL) from 2011 to 2016. The goal of this work was to develop empirical correlation models to predict these values from measurable variables from the chemical process so that these quantities could be predicted a-priori from the sludge or simulant composition and measurable processing variables. The need for these predictions arises from the need to predict the REDuction/OXidation (REDOX) state of the glass from the Defense Waste Processing Facility (DWPF) melter. This report summarizes the work on these correlations based on the aforementioned data. Previous work on these correlations was documented in a technical report covering data from 2011-2015. This current report supersedes this previous report. Further refinement of the models as additional data are collected is recommended. The glass REDOX depends on the concentrations of nitrate and manganese (oxidants), and of glycolate, formate, oxalate, carbon, and antifoam (reductants) in the melter feed. The waste sludge contains nitrite, nitrate, manganese (Mn), and oxalate. Virtually all of the nitrite is converted to nitrate or NO+NO 2 +N 2 O gases in the CPC. The portion of the nitrite converted to nitrate increases the amount of nitrate in the sludge. The amount of glycolate in the final melter feed depends on the amount of the glycolic acid feed that is destroyed. Similarly, the amounts of formate and oxalate formed during the decomposition of glycolic acid are required. The material balance on carbon was found to not close in most cases. Generally, there was less carbon at the end of testing compared to the inputs. The most uncertain product variable was glycolate, so material balances were performed where the glycolate concentration was adjusted, usually upward, to close the balance. Correlation versus the original, as-measured, data was generally poor, but correlation against the material balance adjusted values was greatly improved. It was also shown that the correlation of the measured REDOX versus the predicted REDOX was much better when the material balance adjusted glycolate values were used. Three data series were primarily used during the regressions of the data; these series were 1) Sludge Batch 9 NG flowsheet simulant runs NG51-62 (SB9-NG); 2) Scaled Runs + Bounding Hydrogen Runs (SR+BH); and 3) Runs GN43-50 and 57 (43-50,57). The glycolate destruction was found to correlate with acid stoichiometry (AS), percent reducing acid (PRA), and for some data series, headspace to simulant volume ratio (HSV), mercury (Hg), and nitrate. Although glycolate destruction for pairs of data series (e.g., [SB9-NG] and [SR+BH]) were found to depend on HSV, the combination of all three data series was not found to have significant dependence on this variable. The best model for glycolate destruction depended on AS, nitrate, and Hg. This model predicted the product glycolate compositions of the data to within 92-106%. The conversion of glycolate to formate was high when noble metals and Hg were not present, with values up to 100%. When noble metals and Hg were present, this conversion ranged from zero to 7%, and was dependent on AS. Lower AS gave higher conversions to formate. The conversion to oxalate was found to depend on the AS and the initial concentration of nitrite. An alternative fit versus AS and the form of ruthenium (Ru) used is a possible alternative. This fit was somewhat less statistically significant. This second model predicts that more oxalate is formed when Ru-nitrosyl nitrate is used rather than Ru chloride. The conversion of glycolate to oxalate ranged from zero to 6%. The conversion of nitrite to nitrate depended primarily on AS and PRA, with HSV and Hg being significant when these variables were varied. For multiple series of data, nitrite was also needed to SRNL-STI-2017-00172 5HYLVLRQ viL distinguish between data series, and the effect of HSV became insignificant. The best model for nitrite to nitrate conversion depended on AS, PRA, nitrite, and Hg. The 95% confidence intervals on the predicted values of glycolate destruction, glycolate to oxalate conversion, and nitrite to nitrate conversion were used to determine the uncertainty in the predicted REDOX when starting with only the composition of the sludge, AS, and PRA. Using the 95% confidences on an individual value (that is the confidence in getting a particular value for one single test as opposed to what the mean would be for multiple tests), the uncertainty in the predicted REDOX was calculated. The uncertainty in the actual product composition glycolate, oxalate, formate, and nitrate concentrations translated to an uncertainty in the REDOX value of ±0.1,which is approximately the uncertainty claimed in the REDOX model itself.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2020

Image Collection Simulation Using High-Resolution Atmospheric Modeling

A new method is described for simulating the passive remote sensing image collection of ground targets that includes effects from atmospheric physics and dynamics at fine spatial and temporal scales. The innovation in this research is the process of combining a high-resolution weather model with image collection simulation to attempt to account for heterogeneous and high-resolution atmospheric effects on image products. The atmosphere was modeled on a 3D voxel grid by a Large-Eddy Simulation (LES) driven by forcing data constrained by local ground-based and air-based observations. The spatial scale of the atmospheric model (10–100 m) came closer than conventional weather forecast scales (10–100 km) to approaching the scale of typical commercial multispectral imagery (2 m). This approach was demonstrated through a ground truth experiment conducted at the Department of Energy Atmospheric Radiation Measurement Southern Great Plains site. In this experiment, calibrated targets (colored spectral tarps) were placed on the ground, and the scene was imaged with WorldView-3 multispectral imagery at a resolution enabling the tarps to be visible in at least 9–12 image pixels. The image collection was simulated with Digital Imaging and Remote Sensing Image Generation (DIRSIG) software, using the 3D atmosphere from the LES model to generate a high-resolution cloud mask. The high-resolution atmospheric model-predicted cloud coverage was usually within 23% of the measured cloud cover. The simulated image products were comparable to the WorldView-3 satellite imagery in terms of the variations of cloud distributions and spectral properties of the ground targets in clear-sky regions, suggesting the potential utility of the proposed modeling framework in improving simulation capabilities, as well as testing and improving the operation of image collection processes.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS1994

A simulation of the intracluster medium with feedback from cluster galaxies

We detail method and report first results from a three-dimensional hydrodynamical and N-body simulation of the formation and evolution of a Coma-sized cluster of galaxies, with the intent of studying the history of the hot, X-ray emitting intracluster medium. Cluster gas, galaxies, and dark matter are included in the model. The galaxies and dark matter fell gravitational forces; the cluster gas also undergoes hydrodynamical effects such as shock heating and PdV work. For the first time in three dimensions, we include modeling of ejection of processed gas from the simulated galaxies by winds, including heating and heavy element enrichment. For comparison, we employ a `pure infall' simulation using the same initial conditions but with no galaxies or winds. We employ an extreme ejection history for galactic feedback in order to define the boundary of likely models. As expected, feedback raises the entropy of the intracluster gas, preventing it from collapsing to densities as high as those attained in the infall model. The effect is more pronounced in subclusters formed at high redshift. The cluster with feedback is always less X-ray luminous, but experiences more rapid luminosity evolution, than the pure infall cluster. Even employing an extreme ejection model, the final gas temperature is only approximately 15% larger than in the infall model. The radial temperature profile is very nearly isothermal within 1.5 Mpc. The cluster galaxies in the feedback model have a velocity dispersion approximately 15% lower than the dark matter. This results in the true ratio of specific energies in galaxies to gas being less than one, beta(sub spec) approximately 0.7. The infall model predicts beta(sub spec) approximately 1.2. Large excursions in these values occur over time, following the complex dynamical history of the cluster. The morphology of the X-ray emission is little affected by feedback. The emission profiles of both clusters are well described by the standard beta-model with beta(sub fit) approximately equal to 0.7 - 0.9. X-ray mass estimates based on the assumptions of hydrostatic equilibrium and the applicability of the beta-model are quite accurate in both cases. A strong, radial iron abundance gradient is present, which develops as a consequence of the steepening of the galaxy density profile over time. Spectroscopic observations using nonimaging detectors with wide (approximately 45 min) fields of view dramatically smear the gradient. Observations with arcminute resolution, made available with the ASCA satellite, would readily resolve the gradient.

Metzler, Christopher A.↗
NASA NTRS1993

Models of the Cartwheel ring galaxy: Spokes and starbursts

Recent observations of this famous ring galaxy, including optical and near-infrared CCD surface photometry, and VLA radio continuum and 21 cm line mapping (Higdon 1992b, in prep.), have inspired a renewed modeling effort. Toomre's (1978, in The Large-scale Structure of the Universe, eds. Longair and Einasto) series of restricted three-body simulations demonstrated how the multiple rings could be produced in a nearly head-on galaxy collision. New models with a halo-dominated potential based on the 21 cm rotation curve are able to reproduce such details as the spacing between rings, ring widths, offset of the nucleus, and several kinematical features, thus providing strong support for the collisional theory. The new observations have shown there are little or no old stars in Cartwheel; it may consist almost entirely of gas and stars produced as a result of compression in the ring wave. To model this process Smooth Particle Hydrodynamics (SPH) simulations of the Cartwheel disk have been performed. Fixed gravitational potentials were used to represent the Cartwheel and a roughly 30 percent mass collision partner. The interaction dynamics was treated as in the usual restricted three-body approximation, and the effects of local self-gravity between disk particles were calculated. We are particularly interested in testing the theory that enhanced star formation in waves is the result of gravitational instability in the compressed region (see e.g. Kennicutt 1989, ApJ 344, 685). The gas surface density in a number of simulations was initialized to a value slightly below the threshold for local gravitational instability throughout most of the disk. The first ring wave produces relatively modest compressions (a factor of order a few), triggering instability in a narrow range of wavelengths. Self-gravity in the disk is calculated over a comparable range of scales. Simulations were run with isothermal, adiabatic, and adiabatic with radiative cooling characterized by a relatively short timescale. The isothermal approximation is good except in the vicinity of the strong second (inner) ring, and several snapshots from one case are shown in the figure below. Flocculent spiral segments are present before the collision, and these are compressed into dense knots in the ring wave. These knots are likely to be sites of vigorous star formation. In the strong rarefaction behind the outer ring most of the knots are radially stretched and sheared, giving rise to spoke-like features. A few dense knots are evidently very tightly bound, because they retain their coherence and are stretched relatively little through the rarefaction. This is in accord with evidence for continuing star formation in some spokes (Marcum, Appleton and Higdon 1992). The number and spacing of spokes is a direct function of the scale of the gravitational instability in the disk. Thus, the gravitational instability theory, together with the hypothesis that massive stars are only formed in dense knots of gas, can account for most of the distinct morphology of the Cartwheel.

Struck-Marcell, Curtis↗