Search NASA⌕ Search

SEARCH · Search NASA

Results for “Energy Code Field Study”

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 91 records · Page 5
All sourcesNASA NTRSDOE OSTI
NASA NTRS2005

Analysis of Bare-Tether Systems as a Thruster for MXER Studies

The concept of the electrodynamic propulsion has a number of attractive features and has been widely discussed for different applications. A number of system designs have been proposed and compared during the last ten years. In spite of this, the choice of the proper design, for a specific mission, is far from evident. Such characteristics of tether performance as system acceleration, efficiency, etc. should be calculated and compared. The code that calculates the current for bare and partly insulated tethers with circular (wire) and rectangle (tape) cross-sections is presented. It takes into account the corrections to the OML current due to the tether cross-section geometry and the magnetic field produced by the tether current. There are two options in this code: for current calculation with the prescribed energy supply and with the prescribed end-point potential. This permits to calculate the parameters characterizing tether performance. Results for the current calculated for tethers with different designs for the currently proposed Momentum exchange Electrodynamic Reboost (MXER) Tether System are presented.

Khazanov, G. V.↗
NASA NTRS2005

Analysis of Bare-Tether Systems as a Thruster for MXER Studies

The concept of the electrodynamic propulsion has a number of attractive features and has been widely discussed for different applications. A number of system designs have been proposed and compared during the last ten years. In spite of this, the choice of the proper design, for a specific mission, is far from evident. Such characteristics of tether performance as system acceleration, efficiency, etc. should be calculated and compared. The code that calculates the current for bare and partly insulated tethers with circular (wire) and rectangle (tape) cross-sections is presented. It takes into account the corrections to the OML current due to the tether cross-section geometry and the magnetic field produced by the tether current. There are two options in this code: for current calculation with the prescribed energy supply and with the prescribed end-point potential. This permits us to calculate the parameters characterizing tether performance. Results for the current calculated for tethers with different designs for the currently proposed Momentum exchange Electrodynamic Reboost (MXER) Tether System are presented.

Khazanov, G. V.↗
NASA NTRS1991

A study of physical mechanisms for filament eruption and coronal mass ejection via numerical simulation

It is well-known from both observation and theory that photospheric shear motion has played a key role in physical processes of the energy build-up and release for the solar flare. In order to further our understanding of the occurrence of solar flares it is necessary to investigate the triggering mechanism. One popular scenario for the onset of eruptive solar flares is that in response to photospheric shear motions the magnetic field evolves slowly through a series of magnetohydrodynamic-equilibria until a threshold is reached where magnetohydrodynamic (MHD) non-equilibrium sets in. Thus, a magnetic eruption occurs, causing the solar flare. To substantiate this claim we have employed our newly developed three-dimensional time-dependent MHD code with gravity to simulate the evolution of the coronal field. We use plasma beta = 0.1 to closely approximate the condition in the actual corona. Some preliminary results are presented.

Wu, S. T.↗
NASA NTRS1987

Computational chemistry

With the advent of supercomputers, modern computational chemistry algorithms and codes, a powerful tool was created to help fill NASA's continuing need for information on the properties of matter in hostile or unusual environments. Computational resources provided under the National Aerodynamics Simulator (NAS) program were a cornerstone for recent advancements in this field. Properties of gases, materials, and their interactions can be determined from solutions of the governing equations. In the case of gases, for example, radiative transition probabilites per particle, bond-dissociation energies, and rates of simple chemical reactions can be determined computationally as reliably as from experiment. The data are proving to be quite valuable in providing inputs to real-gas flow simulation codes used to compute aerothermodynamic loads on NASA's aeroassist orbital transfer vehicles and a host of problems related to the National Aerospace Plane Program. Although more approximate, similar solutions can be obtained for ensembles of atoms simulating small particles of materials with and without the presence of gases. Computational chemistry has application in studying catalysis, properties of polymers, all of interest to various NASA missions, including those previously mentioned. In addition to discussing these applications of computational chemistry within NASA, the governing equations and the need for supercomputers for their solution is outlined.

Arnold, J. O.↗
NASA NTRS1991

Definition of the unsteady vortex flow over a wing/body configuration

A problem of current interest in computational aerodynamics is the prediction of unsteady vortex flows over aircraft at high angles of attack. A six-month experimental effort was conducted at the John H. Harper Wind Tunnel to acquire qualitative and quantitative information on the unsteady vortex flow over a generic wing-body configuration at high angles of attack. A double-delta flat-plate wing with beveled edges was combined with a slender sharp-nosed body-of-revolution fuselage to form the generic configuration. This configuration produces a strong attached leading edge vortex on the wing, as well as sharply-peaked flow velocity spectra above the wing. While it thus produces flows with several well-defined features of current interest, the model was designed for efficiency of representation in computational codes. A moderate number of surface pressure ports and two unsteady pressure sensors were used to study the pressure distribution over the wing and body surface at high angles of attack; the unsteady pressure sensing did not succeed because of inadequate signal-to-noise ratio. A pulsed copper vapor laser sheet was used to visualize the vortex flow over the model, and vortex trajectories, burst locations, mutual induction of vortex systems from the forebody, strake, and wing, were quantified. Laser Doppler velocimetry was used to quantify all 3 components of the time-average velocity in 3 data planes perpendicular to the freestream direction. Statistics of the instantaneous velocity were used to study intermittency and fluctuation intensity. Hot-film anemometry was used to study the fluctuation energy content in the velocity field, and the spectra of these fluctuations. In addition, a successful attempt was made to measure velocity spectra, component by component, using laser velocimetry, and these were compared with spectra measured by hot-film anemometry at several locations.

Liou, S. G.↗
NASA NTRS2004

Magnetic-Nozzle Studies for Fusion Propulsion Applications: Gigawatt Plasma Source Operation and Magnetic Nozzle Analysis

This project has been a multiyear effort to assess the feasibility of a key process inherent to virtually all fusion propulsion concepts: the expansion of a fusion-grade plasma through a diverging magnetic field. Current fusion energy research touches on this process only indirectly through studies of plasma divertors designed to remove the fusion products from a reactor. This project was aimed at directly addressing propulsion system issues, without the expense of constructing a fusion reactor. Instead, the program designed, constructed, and operated a facility suitable for simulating fusion reactor grade edge plasmas, and to examine their expansion in an expanding magnetic nozzle. The approach was to create and accelerate a dense (up to l0(exp 20)/m) plasma, stagnate it in a converging magnetic field to convert kinetic energy to thermal energy, and examine the subsequent expansion of the hot (100's eV) plasma in a subsequent magnetic nozzle. Throughout the project, there has been a parallel effort between theoretical and numerical design and modelling of the experiment and the experiment itself. In particular, the MACH2 code was used to design and predict the performance of the magnetoplasmadynamic (MPD) plasma accelerator, and to design and predict the design and expected behavior for the magnetic field coils that could be added later. Progress to date includes the theoretical accelerator design and construction, development of the power and vacuum systems to accommodate the powers and mass flow rates of interest to out research, operation of the accelerator and comparison to theoretical predictions, and computational analysis of future magnetic field coils and the expected performance of an integrated source-nozzle experiment.

Gilland, James H.↗
NASA NTRS1991

Summary of completed project

The primary goal of the work in this proposal has been the determination of the requirement for acceleration of particles, especially electrons, during solar flares. This is one of the important and yet unsolved problems in solar flares and other astronomical objects. This is done by detailed theoretical investigation of interaction of accelerated particles with ambient plasma including interaction with particles, waves, and large scale magnetic fields. Comparison with observed impulsive phase radiation such as X-rays, gamma-rays, microwaves, and type III radiation is used for constraining the models and for determination of the characteristics of the accelerated particles. Steady state and time-dependent numerical codes have been developed based on the Fokker-Planck method for study of the particle acceleration and transport during a flare. Using these, the correlation between characteristics of radiation at various wavelengths have been investigated and some constraints have been set on the model parameters. During the last three year period of this grant, dissertation research of three graduate students have been partially supported by this grant. Each of these theses has dealt with various observed properties of flare emission with particular emphasis on the interrelationship between radiation in at least two energy bands: X-rays and gamma-rays, X-rays and microwaves, and X-rays and type III radiation. This has allowed considerable tightening of constraints on models, and inclusion of wave particle interaction has allowed direct calculations of acceleration processes.

Petrosian, Vahe↗
NASA NTRS2014

A System of Conservative Regridding for Ice-Atmosphere Coupling in a General Circulation Model (GCM)

The method of elevation classes, in which the ice surface model is run at multiple elevations within each grid cell, has proven to be a useful way for a low-resolution atmosphere inside a general circulation model (GCM) to produce high-resolution downscaled surface mass balance fields for use in one-way studies coupling atmospheres and ice flow models. Past uses of elevation classes have failed to conserve mass and energy because the transformation used to regrid to the atmosphere was inconsistent with the transformation used to downscale to the ice model. This would cause problems for two-way coupling. A strategy that resolves this conservation issue has been designed and is presented here. The approach identifies three grids between which data must be regridded and five transformations between those grids required by a typical coupled atmosphere-ice flow model. This paper develops a theoretical framework for the problem and shows how each of these transformations may be achieved in a consistent, conservative manner. These transformations are implemented in Glint2, a library used to couple atmosphere models with ice models. Source code and documentation are available for download. Confounding real-world issues are discussed, including the use of projections for ice modeling, how to handle dynamically changing ice geometry, and modifications required for finite element ice models.

Ice-atmosphere↗
NASA NTRS2023

Automated Fluidics Device for Extraction and Quantification of miRNA Biomarkers From Blood

Radiation Assessment DuRing Exposure And long-Duration Spaceflight (RADREADS) demonstrates space-compatible point-of-care technology for quantitative biological monitoring of blood miRNA biomarkers in response to long-term low dose radiation exposure. This individualized monitoring approach will inform targeted treatment strategies to maximize medical resource utilization by accounting for individual susceptibility to radiation-related illnesses. As human spaceflight progresses beyond Earth’s magnetic shielding, radiation exposure poses a significant risk to astronaut health and safety. Extended operation in this environment comes with an increased risk of radiation exposure, leading to higher risks of radiation sickness, cancer, central nervous system effects, and degenerative diseases. While conventional physical dosimetry techniques capture radiation dose, individualistic susceptibility to radiation damage is varied. Multiple characteristics, including age, body weight, sex, genetics, and immune status, have been found to influence radiosensitivity (Liu et al. 2011, and Bouffler 2016). This differential response necessitates individualized monitoring and targeted treatment strategies to maximize medical resource utilization; however, a practical diagnostic platform for quantifying long-term, low dose radiation-induced tissue damage does not currently exist. MicroRNAs (miRNAs) are a class of small, non-coding RNAs that regulate gene expression by mediating the degradation of messenger RNA. The levels of particular miRNAs are influenced by biological processes such as inflammation and serve as biomarkers for a variety of conditions including cancer (Singh et al. 2017). MicroRNAs are found in various bodily fluids and are amenable to collection via liquid biopsies, providing a minimally invasive and easily quantifiable readout for a variety of radiosensitive reporters. A preliminary signature of 15 spaceflight sensitive miRNA has been identified in rodent and human studies, including miR-21-5p, miR-24-3p, miR-92a-3p, miR-17-5p, miR-16a-3p, miR-34a-3p, and miR-223-3p. These targets generally increased expression with radiation dose and linear energy transfer, though variation between individuals is not yet described. Current gaps in the field include a lack of understanding of longitudinal biological responses to long-term, low dose radiation exposure and the absence of space-compatible point-of-care technology for quantitative biological monitoring. In this body of work, we aim to develop an automated bleed-to-read system to process whole blood for the detection of miRNA biomarkers in order to monitor individualistic responses to radiation exposure. This will be achieved via separating serum (or plasma) from whole blood, followed by extraction, amplification, and quantification of the miRNA using a RT-qPCR reaction. Previously, the WetLab-2 hardware enabled execution of a RT-qPCR reaction aboard ISS; however, it is a manual system that requires crew manipulation and bulky components (Parra et al. 2017). To address these issues, automated fluid handling hardware was developed for each stage of sample preparation. Extraction of total RNA is achieved by sequentially pumping reagents through an off-the-shelf nucleic acid binding column (miRNeasy Serum/Plasma Advanced Kit, Qiagen). This approach eliminates several manual pipetting and centrifuging steps and limits the use of toxic chemicals commonly found in other sample processing techniques. The resulting elution will then be automatically dispensed for RT-qPCR analysis using a compact rotary qPCR (Mic qPCR Cycler, Bio Molecular Systems) that will improve spaceflight compatibility by removing bubbles from the detection region, another challenge highlighted by WetLab-2 (Parra et al. 2017). Efforts are also being made to simplify the RT-qPCR reaction to a 1-step air-dryable mix to improve long-term reagent stability at room temperature and reduce system complexity. By automating the RT-qPCR processes via microfluidic manipulation, RADREADS will reduce crewmember hands-on time and enable the personalized detection of radiation-induced tissue damage during long duration missions. Minimally invasive, longitudinal monitoring of individual’s response to radiation exposure will inform how the physiological system responds to long-term low dose space radiation and enables development of targeted countermeasures by the medical team. Ultimately, this portable technology will require minimal technical expertise and can also be used to monitor miRNA biomarkers associated with other diseases.

Tristen Head↗
NASA NTRS1996

Excitation of Plasma Waves in Aurora by Electron Beams

In this paper, we study numerically the excitation of plasma waves by electron beams, in the auroral region above 2000 km of altitude. We have solved the fully kinetic dispersion relation, using numerical method and found the real frequency and the growth rate of the plasma wave modes. We have examined the instability properties of low-frequency waves such as the Electromagnetic Ion Cyclotron (EMIC) wave as well as Lower-Hybrid (LH) wave in the range of high-frequency. In all cases, the source of free energy are electron beams propagating parallel to the geomagnetic field. We present some features of the growth rate modes, when the cold plasma parameters are changed, such as background electrons and ions species (H(+) and O(+)) temperature, density or the electron beam density and/or drift velocity. These results can be used in a test-particle simulation code, to investigate the ion acceleration and their implication in the auroral acceleration processes, by wave-particle interaction.

daSilva, C. E.↗
NASA NTRS1993

Hybrid simulations of thin current sheets

A one-dimensional, hybrid simulation code is used to study current sheets with a nonzero normal magnetic field B(sub z) and a dawn-to-dusk electric field E(sub y). Such configurations are dependent upon only two parameters: we use the normalized normal magnetic field B-normalized (sub z) = B(sub z)/(4(pi)(n(sub b)) (v(exp 2 sub T))(exp 1/2) and normalized electric field V-normalized (sub D) = (1/V(sub T)(cE(sub y)/B(sub z)), where V(sub T) is the thermal velocity of ions prior to their interaction with the current sheet and n(sub b) is the number density outside the current sheet (at the simulation boundary). A third parameter that is relevant to the motion of particles in current sheets is kappa(sub A), the value of kappa = (R(sub min)/rho(sub max))(exp 1/2) for particles of average energy. We find that if either B-normalized (sub z) is close to or greater than 1, or if kappa(sub A) is close to 1, a rotational mode develops in which the z = 0 current rotates with the ion sense about the normal magnetic field, while for small values of both B-normalized (sub z) or kappa(sub A), the configuration is quasi-steady. To achieve values of kappa(sub A) of the order of or larger than 1, we decrease the value of V-normalized (sub D) uniformly. We find that the magnetic field fluctuations and particle distribution functions are similar in many respects to what was observed in the day 240, 1986, Active Magnetospheric Particle Tracer Explorer (AMPTE)/CCE current disruption event, an event that appears to be located at the site of initiation of current disruption and related particle energization.

Burkhart, G. R.↗
NASA NTRS2016

eXTP: Enhanced X-Ray Timing and Polarimetry Mission

eXTP is a science mission designed to study the state of matter under extreme conditions of density, gravity and magnetism. Primary goals are the determination of the equation of state of matter at supra-nuclear density, the measurement of QED effects in highly magnetized star, and the study of accretion in the strong-field regime of gravity. Primary targets include isolated and binary neutron stars, strong magnetic field systems like magnetars, and stellar-mass and supermassive black holes. The mission carries a unique and unprecedented suite of state-of-the-art scientific instruments enabling for the first time ever the simultaneous spectral-timing-polarimetry studies of cosmic sources in the energy range from 0.5-30 keV (and beyond). Key elements of the payload are: the Spectroscopic Focusing Array (SFA) - a set of 11 X-ray optics for a total effective area of approx. 0.9 m(exp. 2) and 0.6 m(exp. 2) at 2 keV and 6 keV respectively, equipped with Silicon Drift Detectors offering less than 180 eV spectral resolution; the Large Area Detector (LAD) - a deployable set of 640 Silicon Drift Detectors, for a total effective area of approx. 3.4 m(exp. 2), between 6 and 10 keV, and spectral resolution better than 250 eV; the Polarimetry Focusing Array (PFA) - a set of 2 X-ray telescope, for a total effective area of 250 cm(exp. 2) at 2 keV, equipped with imaging gas pixel photoelectric polarimeters; the Wide Field Monitor (WFM) - a set of 3 coded mask wide field units, equipped with position-sensitive Silicon Drift Detectors, each covering a 90 degrees x 90 degrees field of view. The eXTP international consortium includes major institutions of the Chinese Academy of Sciences and Universities in China, as well as major institutions in several European countries and the United States. The predecessor of eXTP, the XTP mission concept, has been selected and funded as one of the so-called background missions in the Strategic Priority Space Science Program of the Chinese Academy of Sciences since 2011. The strong European participation has significantly enhanced the scientific capabilities of eXTP. The planned launch date of the mission is earlier than 2025.

X-ray astronomy↗
NASA NTRS1981

The small community solar thermal power experiment. Parabolic dish technology for industrial process heat application

Aspects of incorporating a thermal energy transport system (ETS) into a field of parabolic dish collectors for industrial process heat (IPH) applications were investigated. Specific objectives are to: (1) verify the mathematical optimization of pipe diameters and insulation thicknesses calculated by a computer code; (2) verify the cost model for pipe network costs using conventional pipe network construction; (3) develop a design and the associated production costs for incorporating risers and downcomers on a low cost concentrator (LCC); (4) investigate the cost reduction of using unconventional pipe construction technology. The pipe network design and costs for a particular IPH application, specifically solar thermally enhanced oil recovery (STEOR) are analyzed. The application involves the hybrid operation of a solar powered steam generator in conjunction with a steam generator using fossil fuels to generate STEOR steam for wells. It is concluded that the STEOR application provides a baseline pipe network geometry used for optimization studies of pipe diameter and insulation thickness, and for development of comparative cost data, and operating parameters for the design of riser/downcomer modifications to the low cost concentrator.

Polzien, R. E.↗
NASA NTRS1992

Arcing in space structures in low Earth orbit

This report describes results of an experimental and theoretical program to investigate arcing of structures containing dielectric and conducting materials when they are biased negatively with respect to a plasma. An argon ion source generated Ar(+) ions of directed energy 20 to 40 eV and density approximately 10(exp 7) cm(exp -3) that impinged upon samples containing a dielectric material on top of a negatively biased Kovar plate. Arcing events were studied for bias voltages between -300 and -1000V with respect to the ion beam. The samples were Dow Corning 93-500 adhesive on Kovar, fused silica cover slips bonded on Kovar, and silicon solar cells mounted on Kovar. Measurements of discharge current, Kovar plate voltage, and radiation from the arc versus time were carried out. Microsecond duration exposure images and optical spectra in the 0.24 to 0.40 micron band were also acquired during arcing events. Arcing events were found to be associated with exposed adhesive and means were found to eliminate arcing altogether. The charging of a silica cover plate and the fields around the plate were calculated using a particle-in-cell code. Models were developed to explain the ignition of the arc and the physical processes occurring during the discharge.

Upschulte, B. L.↗
NASA NTRSDate not supplied

What Processes are Defining the Ionospheric Conductivity and its Variability During Geomagnetic Disturbances?

Modeling of electrodynamic coupling between the magnetosphere, ionosphere, and upper atmosphere (MIA) depends on accurate specification of ionospheric conductance produced by auroral precipitation of high-energy electrons and ions. The precipitation of energetic electrons into the ionosphere is the result of a three-step process that relies on the proper selections of the simulation tools for the ionospheric conductivity studies, while observation can only measure the results of the three steps. In the region of diffuse aurora, the first step is the initiation of electron precipitation into both magnetically conjugate foot points from the Earth’s magnetosphere via wave-particle interactions. The second step is the multiple atmospheric backscatters (or reflections) of electrons at the two magnetic conjugate points, which produces secondary superthermal electron fluxes. The third step is namely the self-consistent electric and magnetic fields that influence magnetospheric particle transport and re-distribute precipitating electrons and ions through the ionospheric electrodynamics. These steps are especially important for revealing electron precipitation dynamics that carry most of the energy in the aurora, resulting also from ion precipitation production and the formation of ionospheric conductance during geomagnetic disturbances. We demonstrate all above results based on SuperThermal Electron Transport (STET), Superthermal Proton, Electron and Atomic Hydrogen tRansport in the Ionosphere and Thermosphere (SPEAH-RIT), and Comprehensive Inner Magnetosphere and Ionospere (CIMI) codes developed at NASA Goddard Space Flight Center.

George V. Khazanov↗
NASA NTRS2001

Coarse Grid Modeling of Turbine Film Cooling Flows Using Volumetric Source Terms

The recent trend in numerical modeling of turbine film cooling flows has been toward higher fidelity grids and more complex geometries. This trend has been enabled by the rapid increase in computing power available to researchers. However, the turbine design community requires fast turnaround time in its design computations, rendering these comprehensive simulations ineffective in the design cycle. The present study describes a methodology for implementing a volumetric source term distribution in a coarse grid calculation that can model the small-scale and three-dimensional effects present in turbine film cooling flows. This model could be implemented in turbine design codes or in multistage turbomachinery codes such as APNASA, where the computational grid size may be larger than the film hole size. Detailed computations of a single row of 35 deg round holes on a flat plate have been obtained for blowing ratios of 0.5, 0.8, and 1.0, and density ratios of 1.0 and 2.0 using a multiblock grid system to resolve the flows on both sides of the plate as well as inside the hole itself. These detailed flow fields were spatially averaged to generate a field of volumetric source terms for each conservative flow variable. Solutions were also obtained using three coarse grids having streamwise and spanwise grid spacings of 3d, 1d, and d/3. These coarse grid solutions used the integrated hole exit mass, momentum, energy, and turbulence quantities from the detailed solutions as volumetric source terms. It is shown that a uniform source term addition over a distance from the wall on the order of the hole diameter is able to predict adiabatic film effectiveness better than a near-wall source term model, while strictly enforcing correct values of integrated boundary layer quantities.

Heidmann, James D.↗
NASA NTRS2003

Monte Carlo Methods in Materials Science Based on FLUKA and ROOT

A comprehensive understanding of mitigation measures for space radiation protection necessarily involves the relevant fields of nuclear physics and particle transport modeling. One method of modeling the interaction of radiation traversing matter is Monte Carlo analysis, a subject that has been evolving since the very advent of nuclear reactors and particle accelerators in experimental physics. Countermeasures for radiation protection from neutrons near nuclear reactors, for example, were an early application and Monte Carlo methods were quickly adapted to this general field of investigation. The project discussed here is concerned with taking the latest tools and technology in Monte Carlo analysis and adapting them to space applications such as radiation shielding design for spacecraft, as well as investigating how next-generation Monte Carlos can complement the existing analytical methods currently used by NASA. We have chosen to employ the Monte Carlo program known as FLUKA (A legacy acronym based on the German for FLUctuating KAscade) used to simulate all of the particle transport, and the CERN developed graphical-interface object-oriented analysis software called ROOT. One aspect of space radiation analysis for which the Monte Carlo s are particularly suited is the study of secondary radiation produced as albedoes in the vicinity of the structural geometry involved. This broad goal of simulating space radiation transport through the relevant materials employing the FLUKA code necessarily requires the addition of the capability to simulate all heavy-ion interactions from 10 MeV/A up to the highest conceivable energies. For all energies above 3 GeV/A the Dual Parton Model (DPM) is currently used, although the possible improvement of the DPMJET event generator for energies 3-30 GeV/A is being considered. One of the major tasks still facing us is the provision for heavy ion interactions below 3 GeV/A. The ROOT interface is being developed in conjunction with the CERN ALICE (A Large Ion Collisions Experiment) software team through an adaptation of their existing AliROOT (ALICE Using ROOT) architecture. In order to check our progress against actual data, we have chosen to simulate the ATIC14 (Advanced Thin Ionization Calorimeter) cosmic-ray astrophysics balloon payload as well as neutron fluences in the Mir spacecraft. This paper contains a summary of status of this project, and a roadmap to its successful completion.

Pinsky, Lawrence↗
NASA NTRS1993

Numerical modeling of the interstellar medium in galactic disks

We have been developing detailed hydrodynamic models of the global interstellar medium in the hope of understanding the mass and volume occupied by various phases, as well as their structure and kinematics. In our model, the gas is modeled by one fluid while representative Pop 1 stars are modeled by a second fluid. The two fluids are coupled in that the gas forms into stars at a rate given by a Schmidt law while stellar mass loss returns matter into the gas phase (on a time scale of 100 Myr). Also, the stars heat the gas through stellar winds and the gas cools through optically thin radiation. The time behavior of these two fluids is studied in two spatial dimensions with the Eulerian finite difference numerical hydrodynamic code Zen. The two spatial dimensions are along the plane of a disk (x, total length of 2 kpc) and perpendicular to the disk (z, total height of +/- 15 kpc) and a galactic gravitational field in the z direction, typical of that at the solar circle, is imposed upon the simulation; self-gravity and rotation are absent. For the boundary conditions, outflow is permitted at the top and bottom of the grid (z = +/- 15 kpc) while periodic boundary conditions are imposed upon left and right sides of the grid. As initial conditions, we assumed a gaseous distribution like that seen for the H1 by earlier researchers, although the results are insensitive to the initial conditions. We have run simulations in which the heating due to stars, parameterized as a stellar wind velocity, a, is varied from low (a = 150 km/s), to intermediate (a = 300 km/s), to high (a = 600 km/s). Since the intermediate case is roughly equivalent to the Galactic energy injection rate from supernovae, this summary will concentrate on results from this simulation.

Rosen, A.↗