Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “hydrodynamics 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 271 records · Page 15
All sourcesNASA NTRSDOE OSTI
NASA NTRS1991

Non-LTE analysis of the Ofpe/WN9 star HDE 269227 (R84)

The paper presents the results of a spectral analysis of the Ofpe/WN9 star HD 269227 (R84), which assumes a spherically expanding atmosphere to find solutions for equations of radiative transfer. The spectra of hydrogen and helium were predicted with a non-LTE model. Six stellar parameters were determined for R84. The shape of the velocity law is empirically found, since it can be probed from the terminal velocity of the wind. The six stellar parameters are further employed in a hydrodynamic model where stellar wind is assumed to be directed by radiation pressure, duplicating the mass-loss rate and the terminal wind velocity. The velocity laws found by computation and analysis are found to agree, supporting the theory of radiation-driven stellar wind. R84 is surmised to be a post-red supergiant which lost half of its initial mass, possibly during the red-supergiant phase. This mass loss is also suggested by its spectroscopic similarity to S Doradus.

Schmutz, Werner↗
DOE OSTI2022

Model comparisons for two-temperature plasma equations of state

When a plasma is generated in the laboratory, energy is often deposited preferentially into either the electrons or the ions, giving rise to a quasiequilibrium state in which the two species, electrons and ions, are well described by two effective temperatures, T e and T i . Accurate hydrodynamic modeling of such a two-temperature plasma requires an equation of state that captures the relevant many-body physics without assuming a strict local thermodynamic equilibrium. Several models have been proposed within the literature, which extend conventional statistical approaches, each employing a different combination of assumptions for modifying the equilibrium equations. Here, in this work, we compare the predictions for several models, presenting derivations of the internal energy and pressure for each microscopic model within a unified framework so that the assumptions of each model may be more easily compared to one another. We find that for sufficiently weak coupling, all models agree with one another. However, as the coupling strength is increased, the disagreement between the models becomes more pronounced. Moreover, the relative sizes of the corrections predicted by each model depend on which species has the higher temperature, T e >T i vs T e < T i .

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Embedding a critical point in a hadron to quark-gluon crossover equation of state

Lattice QCD simulations have shown unequivocally that the transition from hadrons to quarks and gluons is a crossover when the baryon chemical potential is zero or small. Many model calculations predict the existence of a critical point at a value of the chemical potential where current lattice simulations are unreliable. We show how to embed a critical point in a smooth background equation of state so as to yield the critical exponents and critical amplitude ratios expected of a transition in the same universality class as the liquid-gas phase transition and the three-dimensional Ising model. There are only two independent critical exponents; the relations α + 2β + γ = 2 and β(δ-1) = γ arise automatically, as does a relation between the two critical amplitudes. The resulting equation of state has parameters that may be inferred by hydrodynamic modeling of heavy-ion collisions in the Beam Energy Scan II at the BNL Relativistic Heavy Ion Collider or in experiments at other accelerators.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
NASA NTRS2006

An Initial Non-Equilibrium Porous-Media Model for CFD Simulation of Stirling Regenerators

The objective of this paper is to define empirical parameters (or closwre models) for an initial thermai non-equilibrium porous-media model for use in Computational Fluid Dynamics (CFD) codes for simulation of Stirling regenerators. The two CFD codes currently being used at Glenn Research Center (GRC) for Stirling engine modeling are Fluent and CFD-ACE. The porous-media models available in each of these codes are equilibrium models, which assmne that the solid matrix and the fluid are in thermal equilibrium at each spatial location within the porous medium. This is believed to be a poor assumption for the oscillating-flow environment within Stirling regenerators; Stirling 1-D regenerator models, used in Stirling design, we non-equilibrium regenerator models and suggest regenerator matrix and gas average temperatures can differ by several degrees at a given axial location end time during the cycle. A NASA regenerator research grant has been providing experimental and computational results to support definition of various empirical coefficients needed in defining a noa-equilibrium, macroscopic, porous-media model (i.e., to define "closure" relations). The grant effort is being led by Cleveland State University, with subcontractor assistance from the University of Minnesota, Gedeon Associates, and Sunpower, Inc. Friction-factor and heat-transfer correlations based on data taken with the NASAlSunpower oscillating-flow test rig also provide experimentally based correlations that are useful in defining parameters for the porous-media model; these correlations are documented in Gedeon Associates' Sage Stirling-Code Manuals. These sources of experimentally based information were used to define the following terms and parameters needed in the non-equilibrium porous-media model: hydrodynamic dispersion, permeability, inertial coefficient, fluid effective thermal conductivity (including themal dispersion and estimate of tortuosity effects}, and fluid-solid heat transfer coefficient. Solid effective thermal conductivity (including the effect of tortuosity) was also estimated. Determination of the porous-media model parameters was based on planned use in a CFD model of Infinia's Stirling Technology Demonstration Convertor (TDC), which uses a random-fiber regenerator matrix. The non-equilibrium porous-media model presented is considered to be an initial, or "draft," model for possible incorporation in commercial CFD codes, with the expectation that the empirical parameters will likely need to be updated once resulting Stirling CFD model regenerator and engine results have been analyzed. The emphasis of the paper is on use of available data to define empirical parameters (and closure models) needed in a thermal non-equilibrium porous-media model for Stirling regenerator simulation. Such a model has not yet been implemented by the authors or their associates. However, it is anticipated that a thermal non-equilibrium model such as that presented here, when iacorporated in the CFD codes, will improve our ability to accurately model Stirling regenerators with CFD relative to current thermal-equilibrium porous-media models.

Tew, Roy↗
DOE OSTI2021

Slot-die-coating operability windows for polymer electrolyte membrane fuel cell cathode catalyst layers

Roll-to-roll (R2R) slot-die coating of polymer electrolyte membrane fuel cell (PEMFC) catalyst layers represents a scalable deposition method for producing 10–20 m 2 ·min -1 of catalyst-coated gas diffusion layers (GDLs). This high-throughput production technique will help lower the cost of PEMFC catalyst layers. The uniformity of the wet layer applied by slot die deposition is affected by process parameters such as substrate speed, vacuum pressure applied at the upstream meniscus, gap between the slot die lips and substrate, ink rheology, and other ink and substrate properties. The set of conditions for producing a defect-free coating with a dilute ink typically requires little to no upstream vacuum pressure, so suitable operating conditions can be found easily through trial and error and operator intuition. However, the higher viscosity of more concentrated inks dramatically shifts the range of settings that result in a homogeneous coating to higher vacuum levels, which are harder to find through hit or miss. A predictive model showing the range of operable conditions decreases material wastage inherent in experimentally searching for suitable parameters. In this study, the defect-free coating parameter window is explored experimentally and theoretically for two concentrations of PEFC cathode inks. Both a full capillary hydrodynamic model and a computationally cheaper viscocapillary model successfully predict the experimentally determined coating window within the experimental and model uncertainty limits for inks with 5.3 wt% and 12.0 wt% solids ink while maintaining the 0.1 mg Pt ·cm -2 Pt areal loading target. Here, this paper demonstrates a viable pathway for meeting the $30/kW net ultimate cost target of the United States Department of Energy (U.S. DOE) Hydrogen Fuel Cells Technologies Office (HFTO). The concentrated ink lowers the thermal energy and capital expenditure (CapEx) budget of the coating process by decreasing the amount of time, energy, and floorspace required for drying the coating.

25 ENERGY STORAGE↗
DOE OSTI2022

Embedding a Critical Point in a Hadron to Quark–Gluon Crossover Equation of State

It is shown how to embed a critical point in a smooth background equation of state so as to yield the critical exponents and critical amplitude ratios expected of a transition in the same universality class as the liquid–gas phase transition and the 3D Ising model. There are only two independent critical exponents; the relations α + 2β + $\gamma$ = 2 and β(δ - 1) = $\gamma$ arise automatically, as does a new relation between the two critical amplitudes. The resulting equation of state has parameters which may be inferred by hydrodynamic modeling of heavy-ion collisions in the Beam Energy Scan II at the Relativistic Heavy Ion Collider.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
NASA NTRS2003

Sounding Rocket Microgravity Experiments Elucidating Diffusive and Radiative Transport Effects on Flame Spread over Thermally-Thick Solids

A series of 6-minute microgravity combustion experiments of opposed flow flame spread over thermally-thick PMMA has been conducted to extend data previously reported at high opposed flows to almost two decades lower in flow. The effect of flow velocity on flame spread shows a square root power law dependence rather than the linear dependence predicted by thermal theory. The experiments demonstrate that opposed flow flame spread is viable to very low velocities and more robust than expected from the numerical model, which predicts that at very low velocities (less than 5 centimeters per second), flame spread rates fall off more rapidly as flow is reduced. It is hypothesized that the enhanced flame spread observed in the experiments may be due to three- dimensional hydrodynamic effects, which are not included in the zero-gravity, two-dimensional hydrodynamic model. The effect of external irradiation was found to be more complex that the model predicted over the 0-2 Watts per square centimeter range. In the experiments, the flame compensated for the increased irradiation by stabilizing farther from the surface. A surface energy balance reveals that the imposed flux was at least partially offset by a reduced conductive flux from the increased standoff distance, so that the effect on flame spread was weaker than anticipated.

Olson, Sandra L.↗
NASA NTRS2015

Multiscale and Multiphysics Modeling of Additive Manufacturing of Advanced Materials

The objective of this proposed project is to research and develop a prediction tool for advanced additive manufacturing (AAM) processes for advanced materials and develop experimental methods to provide fundamental properties and establish validation data. Aircraft structures and engines demand materials that are stronger, useable at much higher temperatures, provide less acoustic transmission, and enable more aeroelastic tailoring than those currently used. Significant improvements in properties can only be achieved by processing the materials under nonequilibrium conditions, such as AAM processes. AAM processes encompass a class of processes that use a focused heat source to create a melt pool on a substrate. Examples include Electron Beam Freeform Fabrication and Direct Metal Deposition. These types of additive processes enable fabrication of parts directly from CAD drawings. To achieve the desired material properties and geometries of the final structure, assessing the impact of process parameters and predicting optimized conditions with numerical modeling as an effective prediction tool is necessary. The targets for the processing are multiple and at different spatial scales, and the physical phenomena associated occur in multiphysics and multiscale. In this project, the research work has been developed to model AAM processes in a multiscale and multiphysics approach. A macroscale model was developed to investigate the residual stresses and distortion in AAM processes. A sequentially coupled, thermomechanical, finite element model was developed and validated experimentally. The results showed the temperature distribution, residual stress, and deformation within the formed deposits and substrates. A mesoscale model was developed to include heat transfer, phase change with mushy zone, incompressible free surface flow, solute redistribution, and surface tension. Because of excessive computing time needed, a parallel computing approach was also tested. In addition, after investigating various methods, a Smoothed Particle Hydrodynamics Model (SPH Model) was developed to model wire feeding process. Its computational efficiency and simple architecture makes it more robust and flexible than other models. More research on material properties may be needed to realistically model the AAM processes. A microscale model was developed to investigate heterogeneous nucleation, dendritic grain growth, epitaxial growth of columnar grains, columnar-to-equiaxed transition, grain transport in melt, and other properties. The orientations of the columnar grains were almost perpendicular to the laser motion's direction. Compared to the similar studies in the literature, the multiple grain morphology modeling result is in the same order of magnitude as optical morphologies in the experiment. Experimental work was conducted to validate different models. An infrared camera was incorporated as a process monitoring and validating tool to identify the solidus and mushy zones during deposition. The images were successfully processed to identify these regions. This research project has investigated multiscale and multiphysics of the complex AAM processes thus leading to advanced understanding of these processes. The project has also developed several modeling tools and experimental validation tools that will be very critical in the future of AAM process qualification and certification.

Liou, Frank↗
NASA NTRS1997

Simulation of Ultra-Small Electronic Devices: The Classical-Quantum Transition Region

Concern is increasing about how quantum effects will impact electronic device operation as down-scaling continues along the SIA Roadmap through 2010. This document describes part of a new semiconductor device modeling (SDM) program at NAS to investigate these concerns by utilizing advanced NAS and third-party numerical computation software to rapidly implement and investigate electronic device models including quantum effects. This SDM project will investigate quantum effects in devices in the classical-quantum transition region, especially sub-0.1 mm MOSFETs. Specific tasks planned for this project include the use of quantum corrections to the classical drift-diffusion and hydrodynamic models of electron transport, arid the use of nominally quantum models including significant scattering.

Biegel, Bryan A.↗
DOE OSTI2022

Explicit quantification of residence and flushing times in the Salish Sea using a sub-basin scale shoreline resolving model

The Salish Sea, located in the Pacific Northwest region of North America has complex currents and circulation features distributed over numerous interconnected deep basins with islands. Increased risk of exposure to oil spills and untreated wastewater from maritime emergencies and treatment plant failures have led to a need for quantifying residence and flushing characteristics at a sub-basin scale using a shoreline resolving hydrodynamic model. An unstructured grid model of the Salish Sea was developed using FVCOM with a ≈75-100m shoreline resolution. In addition to 7 tides and 23 salinity and temperature monitoring stations, an extensive currents data set from 135 stations collected over a span of three years was used for skill assessment and validation. Explicit forward computations were then conducted to define and quantify residence and flushing times in various sub-basins of interest using (a) Lagrangian particles and (b) Numerical/virtual dye experiments. The results in most basins show expected seasonal variability with longer flushing time associated with summer lower tides and lower freshwater inflows. However, contrary to expectation, flushing time is significantly longer in wintertime in large fjord-like basins such as Hood Canal (≈138 days), likely due to increased stratification and reduced mixing. The flushing time for the Puget Sound region of the Salish Sea is ≈115 days, while Georgia Basin is 240 days when analyzed as stand-alone basins with zero background concentrations. When examined as part of the flushing of the entire Puget Sound filled with virtual dye, the compounded flushing times for embedded sub-basins can be significantly longer in order of magnitudes and largely dictated by the flushing time of Puget Sound. The computed residence and flushing time scales tabulated over 36 sub-basins provide an improved understanding of water renewal in the system, informing pollution management actions.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS1986

Observational discrimination between modes of shock propagation in interstellar clouds: Predictions of CH+ and SH+ column densities in diffuse clouds

Considerable effort in recent years has been devoted to the study of shocks in the diffuse interstellar medium. This work has been motivated partly by the observations of rotationally excited states of H2, and partly by the realization that species such as CH(+), OH and H2O might be formed preferentially in hot, post-shock gas. The problem of CH(+) and the difficulties encountered when trying to explain the high column densities, observed along lines of sight to certain hot stars, have been reviewed earlier. The importance of a transverse magnetic field on the structure of an interstellar shock was also demonstrated earlier. Transverse magnetic fields above a critical strength give rise to an acceleration zone or precursor, in which the parameters on the flow vary continuously. Chemical reactions, which change the degree of ionization of the gas, also modify the structure of the shock considerably. Recent work has shown that large column densities of CH(+) can be produced in magnetohydrodynamic shock models. Shock speeds U sub s approx. = 10 km/s and initial magnetic field strengths of a few micro G are sufficient to produce ion-neutral drift velocities which can drive the endothermic C(+)(H2,H)CH(+) reaction. It was also shown that single-fluid hydrodynamic models do not generate sufficiently large column densities of CH(+) unless unacceptably high shock velocities (u sub s approx. 20 km/s) are assumed in the models. Thus, the observed column densities of CH(+) provide a constraint on the mode of shock propagation in diffuse clouds. More precisely, they determine a lower limit to the ion-neutral drift velocity.

Flower, D. R.↗
NASA NTRSDate not supplied

Two-Phase Flow in Filters and Orifices: In preparation for Measurements on ISS Packed Bed Reactor Experiment-Water Recovery Series (PBRE-WRS)

Understanding the hydrodynamics of adiabatic two-phase flows in packed bed reactors offers numerous benefits. It enables the optimization of chemical reactions rates and products, crucial for pharmaceuticals and energy production industries. Additionally, this understanding aids in designing more efficient and compact reactors, reducing space and resource requirements. Insights gained from studying such flows in microgravity contribute to advancements of space technologies and the enhancement of our capabilities for undertaking future long duration safe, and sustainable space exploration missions. The Packed Bed Reactor Experiment-Water Recovery Series (PBRE-WRS) is a flight experiment planned for operation on the Microgravity Science Glove box (MSG) facility of the International Space Station (ISS). The objective of the experiment is to validate hydrodynamic models at a range of gas and liquid flow rates. The microgravity environment allows for measurement and observation of aspects of fluid dynamics that are unique when compared to observations made in full or partial gravity environments. The experiment consists of testing two-phase flow hydrodynamics in three different filters, four (4) orifices and one check valve test articles. The experiment utilizes the hardware of the previously flown PBRE and PBRE-2 with packed beds of different packings materials and sizes. The fluid system of the PBRE-WRS consists of a nitrogen gas delivery and water delivery subsystems. The gas delivery subsystem can accommodate two ranges of gas flow rates, just as for the water delivery system. The nitrogen gas loop is open whereas the water loop is closed. Gas is separated from the water using a gas-liquid phase separator which is located in the water delivery module. Preliminary results from ground testing show the pressure drop increasing at different rates with the liquid and gas flow rates in the Brine filter. In this work, a detailed system description is presented along with a summary of results from ground performance testing of different test sections in an end-to-end preflight testing campaign.

Packed Bed Reactor Experiment↗
DOE OSTI2022

Tracking the environmental impacts of ecological engineering on coastal wetlands with numerical modeling and remote sensing

Coastal wetlands are the most valuable ecosystems on the earth but facing severe degradation and losses owing to climate change and anthropogenic activities. Many ecological engineering projects (EEP) have been conducted to mitigate the degradation of coastal wetlands. However, the geomorphological impacts of EEP on coastal wetlands have not been well documented. In this study, a method employed a process-based hydrodynamic model and remote sensing (RS) was developed to evaluate the impacts of EEP on the geomorphological change of a prototype Ramsar site. Results demonstrated that RS has great potential in improving the quality of bathymetry data for the numerical model with a decrease of RMSE from 0.52 m to 0.3 m. It also showed good capacity in trend detection of geomorphological change spatially. Results showed the Chongming Dongtan wetland experienced erosion with an annual rate of -0.035 m/yr from 2013 to 2016 after the implementation of EEP. The deposition rate changed significantly in the area within 200 m of the EEP. It is found that the EEP modified the composition of vegetation, sediment transportation, as well as substrate stability, affecting the geomorphological change of coastal wetlands. It is suggested that the EEP with moderate anthropogenic disturbance is a direct and effective way to recover the coastal habitats for waterbirds. However, the modification of the coastal wetland ecosystem by EEP will lead to the potential vulnerability to global climate change. Therefore, how to mitigate the advantages and disadvantages of the EEP is needed to be further studied to find a more sustainable way for coastal management.

geomorphological change, coastal wetlands, Ecologi↗
NASA NTRS1978

A three-dimensional model of corotating streams in the solar wind. 1: Theoretical foundations

The theoretical and mathematical background pertinent to the study of steady, corotating solar wind structure in all three spatial dimensions (3-D) is discussed. The dynamical evolution of the plasma in interplanetary space (defined as the region beyond roughly 35 solar radii where the flow is supersonic) is approximately described by the nonlinear, single fluid, polytropic (magneto-) hydrodynamic equations. Efficient numerical techniques for solving this complex system of coupled, hyperbolic partial differential equations are outlined. The formulation is inviscid and nonmagnetic, but methods allow for the potential inclusion of both features with only modest modifications. One simple, highly idealized, hydrodynamic model stream is examined to illustrate the fundamental processes involved in the 3-D dynamics of stream evolution. Spatial variations in the rotational stream interaction mechanism were found to produce small nonradial flows on a global scale that lead to the transport of mass, energy, and momentum away from regions of relative compression and into regions of relative rarefaction.

Pizzo, V. J.↗
DOE OSTI2026

Assessing shellfish water exposure to fecal bacteria pollution in Salish Sea: three-dimensional modeling and implications for monitoring

Fecal bacteria (FB) contamination poses significant risks to shellfish safety and management in coastal and estuarine waters. Despite extensive pollution identification and correction efforts, FB contamination in shellfish-growing areas persists in the Salish Sea, highlighting the need to identify overlooked sources and better understand FB transport from riverine and shoreline inputs to shellfish beds. To address this, a high-resolution three-dimensional hydrodynamic model coupled with FB kinetics was developed and applied to a case study site in Salish Sea—Portage Bay—to simulate freshwater plume circulation, flushing dynamics, and bacterial transport. Daily FB loading from the major freshwater inflow—Nooksack River was generated by both linear interpolation and integrating a machine learning approach (XGBoost), trained on historical hydrological and meteorological data. The model successfully reproduced both the magnitude and seasonal variation of FB concentrations in Portage Bay for the year of 2021, demonstrating that simplified FB kinetics with first-order decay due to mortality was effective in this dynamic coastal environment with short flushing time. Model results identified the Nooksack River as the dominant far-field FB source, while scenario simulations showed that near-field coastal stormwater outfalls elevated local FB levels following rainfall, particularly under low-flow conditions. The XGBoost prediction provided comparable or superior accuracy to linear interpolation, particularly during periods of missing observational data, by capturing short-term variability and event-driven loading more effectively. Integrating data-driven riverine FB inputs with mechanistic coastal numerical modeling provides a robust framework for operational forecasting of shellfish bed exposure risk and supports adaptive monitoring and management of shellfish growing areas in the Salish Sea and similar coastal systems.

Salish Sea↗
NASA NTRSDate not supplied

Two-Phase Flow in Filters and Orifices: ISS Packed Bed Reactor Experiment Water Recovery Series (PBRE-WRS)

Understanding the hydrodynamics of adiabatic two-phase flows in packed bed reactors offers numerous benefits. It enables the optimization of chemical reactions rates and products which results in more efficient and compact reactors, thus reducing space and resource requirements, crucial for long duration space missions. The Packed Bed Reactor Experiment-Water Recovery Series (PBRE-WRS) is a flight experiment planned for operation on the Microgravity Science Glove box (MSG) facility of the International Space Station (ISS). The objective of the experiment is to validate hydrodynamic models at a range of gas and liquid flow rates, where these models are used in the design of processes for life support systems in microgravity. The microgravity environment allows for measurement and observation of aspects of fluid dynamics that are unique when compared to observations made in full or partial gravity environments. The experiment consists of testing two-phase flow hydrodynamics in three different filters, four (4) orifices and one check valve test articles. The experiment utilizes the hardware of the previously flown PBRE and PBRE-2 with packed beds of different packings materials and sizes. The fluid system of the PBRE-WRS consists of a nitrogen gas delivery and water delivery subsystems. The gas delivery subsystem can accommodate two ranges of gas flow rates, just as for the water delivery system. The nitrogen gas loop is open whereas the water loop is closed. Gas is separated from the water using a gas-liquid phase separator which is located in the water delivery module. Preliminary results from ground testing show the pressure drop across the filters increasing at different rates with the liquid and gas flow rates. In this work, a detailed system description is presented along with a summary of results from ground performance testing of different test sections in an end-to-end preflight testing campaign.

ISS↗
NASA NTRS1989

Gamma rays, X-rays, and optical light from the cobalt and the neutron star in SN 1987A

Recent developments in modeling the X-ray and gamma-ray emission from SN 1987A are discussed by taking into account both the decaying cobalt and the buried neutron star. The light curve and the spectra evolution of X-rays and gamma-rays are well modeled up to day of about 300 if mixing of Co-56 into hydrogen-rich envelope is assumed. However, the 16-28 keV flux observed by Ginga declines very slowly, whereas the spherical mixing model predicts that the flux should have decreased by a large factor at t greater than 300d. It is shown that this problem can be solved if the photoelectric absorption of X-rays is effectively reduced as a result of the formation of chemically inhomogeneous clumps. Based on the adopted hydrodynamical model and the abundance distribution, predictions are offered for future optical, X-ray, and gamma-ray light curves by taking into account other radioactive sources and various types of the central source, e.g., a buried neutron star accreting the reinfalling material or an isolated pulsar.

Kumagai, Shiomi↗
NASA NTRS1978

A three-dimensional model of corotating streams in the solar wind. I - Theoretical foundations

The paper is concerned with the development of the theoretical and mathematical background pertinent to the study of steady, corotating solar wind structure in all three spatial dimensions. The dynamical evolution of the plasma in interplanetary space (defined as the region beyond roughly 35 Rs where the flow is supersonic) is approximately described by the nonlinear, single-fluid, polytropic magnetohydrodynamic or hydrodynamic equations. Efficient numerical techniques are outlined for solving this complex system of coupled, hyperbolic partial differential equations. The present formulation is inviscid and nonmagnetic, but the methods used allow for the potential inclusion of both features with only modest modifications. A simple, highly idealized hydrodynamic model stream is examined to illustrate the fundamental processes involved in the three-dimensional dynamics of stream evolution. It is found that spatial variations in the rotational stream interaction mechanism produce small nonradial flows on a global scale that lead to the transport of mass, energy, and momentum away from regions of relative compression and into regions of relative rarefaction. Comparison with simpler models demonstrates the essential nonlinear, multidimensional nature of the interplanetary dynamics.

Pizzo, V.↗