Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Hydrodynamical simulations”

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 595 records · Page 33
All sourcesNASA NTRSDOE OSTI
NASA NTRS2001

Quasi-Equilibrium States in the Tropics Simulated by a Cloud-Resolving Model: Specific Features and Budget Analysis - Part 1

A series of long-term integrations using the two-dimensional Goddard Cumulus Ensemble (GCE) model were performed by altering imposed environmental components to produce various quasi-equilibrium thermodynamic states. Model results show that the genesis of a warm/wet quasi-equilibrium state is mainly due to either strong vertical wind shear (from nudging) or large surface fluxes (from strong surface winds), while a cold/dry quasi-equilibrium state is attributed to a remarkably weakened mixed-wind shear (from vertical mixing due to deep convection) along with weak surface winds. In general, latent heat flux and net large-scale temperature forcing, the two dominant physical processes, dominate in the beginning stage of the simulated convective systems, then considerably weaken in the final stage, which leads to quasi-equilibrium states. A higher thermodynamic regime is found to produce a larger rainfall amount, as convective clouds are the leading source of rainfall over stratiform clouds even though the former occupy much less area. Moreover, convective clouds are more likely to occur in the presence of strong surface winds (latent heat flux), while stratiform clouds (especially the well-organized type) are favored in conditions with strong wind shear (large-scale forcing). The convective systems, which consist of distinct cloud types due to the variation in horizontal winds, are also found to propagate differently. Accordingly, convective systems with mixed-wind shear generally propagate in the direction of shear, while the system with strong (multidirectional) wind shear propagates in a more complex way. Based on the results from the temperature (Q1) and moisture (Q2) budgets, cloud-scale eddies are found to act as a hydrodynamic 'vehicle' that cascades the heat and moisture vertically. Several other specific features such as atmospheric stability, CAPE, and mass fluxes are also investigated and found to be significantly different between diverse quasi-equilibrium states. Detailed comparisons between the various states are presented.

Shie, C.-L.↗
NASA NTRS2003

Foam on Tile Impact Modeling for the Space Shuttle Program

Following the breakup of the Space Shuttle Columbia during reentry a NASA-wide investigation team was formed to examine the probable damage inflicted on Orbiter Thermal Protection System (TPS) elements by impact of External Tank insulating foam projectiles. Our team was to apply rigorous, physics-based analysis techniques to help determine parameters of interest for an experimental test program, utilize validated codes to investigate the full range of impact scenarios, and use analysis derived models to predict aero-thermal-structural responses to entry conditions. We were to operate on a non-interference basis with the j Team, and were to supply significant findings to that team and to the Orbiter Vehicle Engineering Working Group, being responsive to any solicitations for support from these entities. The authors formed a working sub-group within the larger team to apply the Smooth Particle Hydrodynamics code SPHC to the damage estimation problem. Numerical models of the LI-900 TPS tiles and of the BX-250 foam were constructed and used as inputs into the code. Material properties needed to properly model the tiles and foam were obtained from other working sub-groups who performed tests on these items for this purpose. Two- and three- dimensional models of the tiles were constructed, including the glass outer layer, the densified lower layer of LI-900 insulation, the Nomex felt Strain Isolation Pad (SIP) mounting layer, and the underlying aluminum 2024 vehicle skin. A model for the BX-250 foam including porous compression, elastic rebound, and surface erosion was developed. Code results for the tile damage and foam behavior were extensively validated through comparison with the Southwest Research Institute (SwRI) foam-on-tile impact experiments carried out in 1999. These tests involved small projectiles striking individual tiles and small tile arrays. Following code and model validation we simulated impacts of larger ET foam projectiles on the TPS tile systems used on the wings of the orbiter. Tiles used on the Wing Acreage, the Main Landing Gear Door, and the Carrier Panels near the front edge of the wing were modeled. Foam impacts shot for the CAB investigation were modeled, as well as impacts at larger angles, including rapid rotation of the projectile, and with varying foam properties. General results suggest that foam impacts on tiles at about 500 mph could cause appreciable damage if the impact angle is greater than about 20 degrees. Some variations of the foam properties, such as increased brittleness or increased density could increase damage in some cases. Rapid (17 rps) rotation failed to increase the damage for the two cases considered. This does not rule out other cases in which the rotational energy might lead to an increase in tile damage, but suggests that in most cases rotation will not be an important factor. Similar models will be applied for other impacting materials, other velocities, and other geometries as part of the Return to Flight process.

Stellingwerf, R. F.↗
NASA NTRS1994

The heating of nova ejecta by radioactive decays of the beta-unstable nuclei

Recent nucleosynthesis and hydrodynamic calculations of the consequences of accretion onto massive ONeMg white dwarf stars show that under certain circumstances significant amounts of the beta-unstable nuclei can be produced and ejected by the resulting explosion. We use these calculations as a guide in order to obtain the conditions under which the heating of the ejected material by the nonthermal electrons and positrons produced by the decays of the beta-unstable nuclei is sufficient to overcome the cooling from adiabatic expansion and lead to the production of X-ray-emitting coronal gas. These conditions are as follows: (1) a mass fraction for Na-22 of the order of 10(exp -3) or greater, (2) an expansion velocity in the range approximately 10(exp 2) - 10(exp 3) km/s, (3) a photospheric radius of approximately 10(exp 14) cm, (4) if the density distribution in the atmosphere satisfies a power law, then the exponent must be less than 3 for heating to overcome adiabatic cooling. Both the simulations of the outburst and the model atmosphere fits to the observed energy distributions, however, imply that the exponent is greater than or = 3 during the early phases of the outburst. Nevertheless, for a value of the exponent of 2, we predict the time when hot coronal gas can form during the expansion phases of the envelope.

Pistinner, Shlomi↗
NASA NTRS2023

Towards an Indian Land Data Assimilation System (ILDAS): A Coupled Hydrologic-Hydraulic System for Water Balance Assessments

Effective management of water resources requires reliable estimates of land surface states and fluxes, including water balance components. But most land surface models run in uncoupled mode and do not produce river discharge at catchment scales to be useful for water resources management applications. Such integrated systems are also rare over India where hydrometeorological extremes have wreaked havoc on the economy and people. So, an Indian Land Data Assimilation System (ILDAS) with a coupled land surface and a hydrodynamic model has been developed and driven by multiple meteorological forcings (0.1°, daily) to estimate land surface states, channel discharge, and floodplain inundation. ILDAS benefits from an integrated framework as well as the largest suite of observation records collected over India and has been used to produce a reanalysis product for 1981–2021 using four forcing datasets, namely, Modern-Era Retrospective Analysis for Research and Applications, Version 2 (MERRA-2), Climate Hazards Group InfraRed Precipitation with Station data (CHIRPS), ECMWF’s ERA-5, and Indian Meteorological Department (IMD) gridded precipitation. We assessed the uncertainty and bias in these precipitation datasets and validated all major components of the terrestrial water balance, i.e., surface runoff, soil moisture, terrestrial water storage anomalies, evapotranspiration, and streamflow, against a combination of satellite and in situ observation datasets. Our assessment shows that ILDAS can represent the hydrological processes reasonably well over the Indian landmass with IMD precipitation showing the best relative performance. Evaluation against ESA-CCI soil moisture shows that MERRA-2 based estimates outperform the others, whereas ERA-5 performs best in simulating evapotranspiration when evaluated against MODIS ET. Evaluations against observed records show that CHIRPS-based estimates have the highest performance in reconstructing surface runoff and streamflow. Once operational, this system will be useful for supporting transboundary water management decision making in the region.

Indian Land Data Assimilation System (ILDAS)↗
NASA NTRS1993

A hydrodynamic treatment of the tilted cold dark matter cosmological scenario

A standard hydrodynamic code coupled with a particle-mesh code is used to compute the evolution of a tilted cold dark matter (TCDM) model containing both baryonic matter and dark matter. Six baryonic species are followed, with allowance for both collisional and radiative ionization in every cell. The mean final Zel'dovich-Sunyaev y parameter is estimated to be (5.4 +/- 2.7) x 10 exp -7, below currently attainable observations, with an rms fluctuation of about (6.0 +/- 3.0) x 10 exp -7 on arcmin scales. The rate of galaxy formation peaks at a relatively late epoch (z is about 0.5). In the case of mass function, the smallest objects are stabilized against collapse by thermal energy: the mass-weighted mass spectrum peaks in the vicinity of 10 exp 9.1 solar masses, with a reasonable fit to the Schechter luminosity function if the baryon mass to blue light ratio is about 4. It is shown that a bias factor of 2 required for the model to be consistent with COBE DMR signals is probably a natural outcome in the present multiple component simulations.

Cen, Renyue↗
NASA NTRS2015

Turbulent Concentration of mm-Size Particles in the Protoplanetary Nebula: Scale-Dependent Cascades

The initial accretion of primitive bodies (here, asteroids in particular) from freely-floating nebula particles remains problematic. Traditional growth-by-sticking models encounter a formidable "meter-size barrier" (or even a mm-to-cm-size barrier) in turbulent nebulae, making the preconditions for so-called "streaming instabilities" difficult to achieve even for so-called "lucky" particles. Even if growth by sticking could somehow breach the meter size barrier, turbulent nebulae present further obstacles through the 1-10km size range. On the other hand, nonturbulent nebulae form large asteroids too quickly to explain long spreads in formation times, or the dearth of melted asteroids. Theoretical understanding of nebula turbulence is itself in flux; recent models of MRI (magnetically-driven) turbulence favor low-or- no-turbulence environments, but purely hydrodynamic turbulence is making a comeback, with two recently discovered mechanisms generating robust turbulence which do not rely on magnetic fields at all. An important clue regarding planetesimal formation is an apparent 100km diameter peak in the pre-depletion, pre-erosion mass distribution of asteroids; scenarios leading directly from independent nebula particulates to large objects of this size, which avoid the problematic m-km size range, could be called "leapfrog" scenarios. The leapfrog scenario we have studied in detail involves formation of dense clumps of aerodynamically selected, typically mm-size particles in turbulence, which can under certain conditions shrink inexorably on 100-1000 orbit timescales and form 10-100km diameter sandpile planetesimals. There is evidence that at least the ordinary chondrite parent bodies were initially composed entirely of a homogeneous mix of such particles. Thus, while they are arcane, turbulent concentration models acting directly on chondrule size particles are worthy of deeper study. The typical sizes of planetesimals and the rate of their formation can be estimated using a statistical model with properties inferred from large numerical simulations of turbulence. Nebula turbulence is described by its Reynolds number Re = (L/eta)(exp 4/3), where L = H alpha(exp 1/2) is the largest eddy scale, H is the nebula gas vertical scale height, alpha the turbulent viscosity parameter, and eta is the Kolmogorov or smallest scale in turbulence (typically about 1km), with eddy turnover time t(sub eta). In the nebula, Re is far larger than any numerical simulation can handle, so some physical arguments are needed to extend the results of numerical simulations to nebula conditions. In this paper, we report new physics to be incorporated into our statistical models.

mm-size particles↗
NASA NTRS1993

Steady state and dynamical structure of a cosmic-ray-modified termination shock

A hydrodynamic model is developed for the structure of a cosmic-ray-modified termination shock. The model is based on the two-fluid equations of diffuse shock acceleration (Drury and Volk, 1981). Both the steady state structure of the shock and its interaction with outer heliospheric disturbances are considered. Under the assumption that the solar wind is decelerated by diffusing interstellar cosmic rates, it is shown that the natural state of the termination shock is a gradual deceleration and compression, followed by a discontinuous jump to a downstream state which is dominated by the pressure contribution of the cosmic rays. A representative model is calculated for the steady state which incorporates both interstellar cosmic ray mediation and diffusively accelerated anomalous ions through a proposed thermal leakage mechanism. The interaction of large-scale disturbances with the equilibrium termination shock model is shown to result in some unusual downstream structure, including transmitted shocks and cosmic-ray-modified contact discontinuities. The structure observed may be connected to the 2-kHz outer heliospheric radio emission (Cairns et al., 1992a, b). The time-dependent simulations also demonstrate that interaction with solar wind compressible turbulence (e.g., traveling interplanetary shocks, etc.) could induce the termination shock to continually fluctuate between cosmic-ray-dominated and gas-dynamic states. This fluctuation may represent a partial explanation of the galactic cosmic ray modulation effect and illustrates that the Pioneer and Voyager satellites will encounter an evolving shock whose structure and dynamic properties are strongly influence by the mediation of interstellar and anomalous cosmic rays.

Donohue, D. J.↗
NASA NTRS1999

Hydrodynamic Stability and Magnetic Reconnection in Disks and Stars

The purpose of this grant is to study parametric instability. The simplest example of parametric instability is a harmonic oscillator with a periodic modulation of the spring constant. If the modulation frequency is close to twice the natural frequency of the oscillator, the amplitude of oscillation tends to grow exponentially. The growth rate is proportional to the strength of the modulation, but it also depends upon the closeness to resonance of the two frequencies, and upon natural damping rate or "Q" of the oscillator. Parametric instabilities are very common in physics. A familiar example is a jogger's ponytail--normally a very strongly damped pendulum, it can be destabilized by the variation in effective gravity during the jogger's stride. Observation confirms that the period of the pendulum is half that of the jogger's vertical motion. In astrophysics, parametric instability may occur by external tidal forcing, or by interaction among eigenmodes. In the latter case, an energetic eigenmode may destabilize modes of half its frequency, provided some weak nonlinearity exists to couple them. Under a previous Astrophysical Theory grant (NAGW-2419), the PI discovered a parametric instability of tidally forced disks such as the accretion disks in cataclysmic variables and X ray binaries [2]. The destabilized modes are tightly-wound, incompressible, three-dimensional waves analogous to g-modes and r-modes in stars. Later work has confirmed our analysis [4]. It was hoped that these modes might provide a source of turbulence and angular momentum transport in accretion disks. However, a follow-up investigation of this instability by local numerical simulations, although confirming the analytically estimated growth rates, found negligible angular momentum flux [3]. Other work, partly supported by the ATP, now strongly indicates that the transport mechanism in such disks is magnetohydrodynamic turbulence [6]. Nevertheless, the parametric mechanism may truncate the outer edges of disks in close binaries [2], and it may be important in disks of very low ionization such as protostellar disks, or even cataclysmic-variable disks in quiescence where the MHD mechanism may be ineffective [5]. All analyses up to 1996 were done in a local approximation where the orbital frequency, shear rate, and tidal field were treated as constants. The locally computed growth rate turns out to depend strongly on radius, and it was unclear how to average these local rates to obtain the correct global rate. This is a critical issue for accretion disks in close binaries, because the local growth rate is comparable to the orbital frequency towards the outer edge of the disk but decreases rapidly inwards. Paper #1 examined this issue in a simplified global model where the destabilizing terms vary with position. We found that the global growth rate is essentially equal to the maximum local rate, provided that the latter is smoothed over a radial range equal to the distance that the destabilized wave propagates at its group speed in one growth time. Thus, in an accretion disk, waves would grow rapidly in the outer parts but would propagate both inwards and outwards at a maximum group speed of order the disk thickness divided by the orbital period.

Goodman, Jeremy↗
NASA NTRS2013

Electrical Capacitance Volume Tomography for the Packed Bed Reactor ISS Flight Experiment

Fixed packed bed reactors are compact, require minimum power and maintenance to operate, and are highly reliable. These features make this technology a highly desirable unit operation for long duration life support systems in space. NASA is developing an ISS experiment to address this technology with particular focus on water reclamation and air revitalization. Earlier research and development efforts funded by NASA have resulted in two hydrodynamic models which require validation with appropriate instrumentation in an extended microgravity environment. To validate these models, the instantaneous distribution of the gas and liquid phases must be measured.Electrical Capacitance Volume Tomography (ECVT) is a non-invasive imaging technology recently developed for multi-phase flow applications. It is based on distributing flexible capacitance plates on the peripheral of a flow column and collecting real-time measurements of inter-electrode capacitances. Capacitance measurements here are directly related to dielectric constant distribution, a physical property that is also related to material distribution in the imaging domain. Reconstruction algorithms are employed to map volume images of dielectric distribution in the imaging domain, which is in turn related to phase distribution. ECVT is suitable for imaging interacting materials of different dielectric constants, typical in multi-phase flow systems. ECVT is being used extensively for measuring flow variables in various gas-liquid and gas-solid flow systems. Recent application of ECVT include flows in risers and exit regions of circulating fluidized beds, gas-liquid and gas-solid bubble columns, trickle beds, and slurry bubble columns. ECVT is also used to validate flow models and CFD simulations. The technology is uniquely qualified for imaging phase concentrations in packed bed reactors for the ISS flight experiments as it exhibits favorable features of compact size, low profile sensors, high imaging speed, and flexibility to fit around columns of various shapes and sizes. ECVT is also safer than other commonly used imaging modalities as it operates in the range of low frequencies (1 MHz) and does not radiate radioactive energy. In this effort, ECVT is being used to image flow parameters in a packed bed reactor for an ISS flight experiment.

Tomography↗
NASA NTRS1999

Chain Dynamics in Magnetorheological Suspensions

Magnetorheological (MR) suspensions are composed of colloidal particles which acquire dipole moments when subjected to an external magnetic field. At sufficient field strengths and concentrations, the dipolar particles rapidly aggregate to form long chains. Subsequent lateral cross-linking of the dipolar chains is responsible for a rapid liquid-to-solid-like rheological transition. The unique, magnetically-activated rheological properties of MR suspensions make them ideal for interfacing mechanical systems to electronic controls. Additionally, the ability to experimentally probe colloidal suspensions interacting through tunable anisotropic potentials is of fundamental interest. Our current experimental work has focused on understanding the fluctuations of dipolar chains. It has been proposed by Halsey and Toor (HT) that the strong Landau-Peierls thermal fluctuations of dipolar chains could be responsible for long-range attractions between chains. Such interactions will govern the long-time relaxation of MR suspensions. We have synthesized monodisperse neutrally buoyant MR suspensions by density matching stabilized ferrofluid emulsion droplets with D2O. This allows us to probe the dynamics of the dipolar chains using light scattering without gravitational, interfacial, and polydispersity effects to resolve the short-wavelength dynamics of the dipolar chains. We used diffusing wave spectroscopy to measure these dynamics. The particle displacements at short times that show an independence to the field strength, but at long times exhibit a constrained, sub-diffusive motion that slows as the dipole strength is increased. The experiments are in good qualitative agreement with Brownian dynamics simulations of dipolar chains. Although there have been several important and detailed studies of the structure and interactions in MR suspensions, there has not been conclusive evidence that supports or contradicts the HT model prediction that long-range interactions exist between fluctuating chains of dipolar particles. Resolving this issue would contribute greatly to the understanding of these interesting and important materials. We have begun to test the predictions of the HT model by both examining the dynamics of individual chains and by measuring the forces between dipolar chains directly to accurately and quantitatively assess the interactions that they experience. To do so, we employ optical trapping techniques and video-microscopy to manipulate and observe our samples on the microscopic level. With these techniques, it is possible to observe chains that are fluctuating freely in three-dimensions, independent of interfacial effects. More importantly, we are able to controllably observe the interactions of two chains at various separations to measure the force-distance profile. The techniques also allow us to study the mechanical properties of individual chains and chain clusters. Our work to this point has focused on reversibly-formed dipolar chains due to field induced dipoles where the combination of this chaining, the dipolar forces, and the hydrodynamic interactions that dictate the rheology of the suspensions. One can envision, however, many situations where optical, electronic, or rheological behavior may be optimized with magneto-responsive anisotropic particles. Chains of polarizable particles may have the best properties as they can coil and flex in the absence of a field and stiffen and orient when a field is applied. We have recently demonstrated a synthesis of stable, permanent paramagnetic chains by both covalently and physically linking paramagnetic colloidal particles. The method employed allows us to create monodisperse chains of controlled length. We observed the stability, field-alignment, and rigidity of this new class of materials. The chains may exhibit unique rheological properties in an applied magnetic field over isotropic suspensions of paramagnetic particles. They are also useful rheological models as bead-spring systems. These chains form the basis for our current experiments with optical traps.

Gast, A. P.↗
NASA NTRS1994

Differential rotation and turbulent convection: A new Reynolds stress model and comparison with solar data

In most hydrodynamic cases, the existence of a turbulent flow superimposed on a mean flow is caused by a shear instability in the latter. Boussinesq suggested the first model for the turbulent Reynolds stresses bar-(u(sub i)u(sub j)) in which the mean shear S(sub ij) is the cause (or source) of turbulence represented by the stress bar-(u(sub i)u(sub j)). In the case of solar differential rotation, exactly the reverse physical process occurs: turbulence (which must pre-exist) generates a mean flow which manifests itself in the form of differential rotation. Thus, the Boussinesq model is wholly inadequate because in the solar case, cause and effect are reversed. Since the Boussinesq model is inadequate, one needs an alternative model for the Reynolds stresses. We present a new dynamical model for the Reynolds stresses, convective fluxes, turbulent kinetic energy, and temperature fluctuations. The complete model requires the solution of 11 differential equations. We then introduce a set of simplifying assumptions which reduce the full dynamical model to a set of algebraic Reynolds stress models. We explicitly solve one of these models that entails only one differential equation. The overall agreement with the data is obtained with a model that is neither phenomenological nor one that requires a full numerical simulation, since it is algebraic in nature. The new model can play an important role in understanding the complex physics underlying the interplay between solar differential rotation and convection, as many physical processes can naturally be incorporated into the model.

Canuto, V. M.↗
NASA NTRS2002

Analysis of the Magneto-Hydrodynamic (MHD) Energy Bypass Engine for High-Speed Air-Breathing Propulsion

The performance of the MHD energy bypass air-breathing engine for high-speed propulsion is analyzed in this investigation. This engine is a specific type of the general class of inverse cycle engines. In this paper, the general relationship between engine performance (specific impulse and specific thrust) and the overall total pressure ratio through an engine (from inlet plane to exit plane) is first developed and illustrated. Engines with large total pressure decreases, regardless of cause or source, are seen to have exponentially decreasing performance. The ideal inverse cycle engine (of which the MHD engine is a sub-set) is then demonstrated to have a significant total pressure decrease across the engine; this total pressure decrease is cycle-driven, degrades rapidly with energy bypass ratio, and is independent of any irreversibility. The ideal MHD engine (inverse cycle engine with no irreversibility other than that inherent in the MHD work interaction processes) is next examined and is seen to have an additional large total pressure decrease due to MHD-generated irreversibility in the decelerator and the accelerator. This irreversibility mainly occurs in the deceleration process. Both inherent total pressure losses (inverse cycle and MHD irreversibility) result in a significant narrowing of the performance capability of the MHD bypass engine. The fundamental characteristics of MHD flow acceleration and flow deceleration from the standpoint of irreversibility and second-law constraints are next examined in order to clarify issues regarding flow losses and parameter selection in the MM modules. Severe constraints are seen to exist in the decelerator in terms of allowable deceleration Mach numbers and volumetric (length) required for meaningful energy bypass (work interaction). Considerable difficulties are also encountered and discussed due to thermal/work choking phenomena associated with the deceleration process. Lastly, full engine simulations utilizing inlet shock systems, finite-rate chemistry, wall cooling with thermally balanced engine (fuel heat sink), fuel injection and mixing, friction, etc. are shown and discussed for both the MHD engine and the conventional scramjet. The MHD bypass engine has significantly lower performance in all categories across the Mach number range (8 to 12.2). The lower performance is attributed to the combined effects of 1) additional irreversibility and cooling requirements associated with the MHD components and 2) the total pressure decrease associated with the inverse cycle itself.

Riggins, David W.↗
NASA NTRS1998

Development of Numerical Extended Hydrodynamics for Transition-Regime Non-Equilibrium Flows Encountered in Semiconductor Manufacturing Processes

Six months of funding was received for the proposed three year research program (funding for the period from March 1, 1997 to August 31, 1997). Although the official starting date for the project was March 1, 1997, no funding for the project was received until July 1997. In the funded research period, considerable progress was made on Phase I of the proposed research program. The initial research efforts concentrated on applying the 10-, 20-, and 35-moment Gaussian-based closures to a series of standard two-dimensional non-reacting single species test flow problems, such as the flat plate, couette, channel, and rearward facing step flows, and to some other two-dimensional flows having geometries similar to those encountered in chemical-vapor deposition (CVD) reactors. Eigensystem analyses for these systems for the case of two spatial dimensions was carried out and efficient formulations of approximate Riemann solvers have been formulated using these eigenstructures. Formulations to include rotational non-equilibrium effects into the moment closure models for the treatment of polyatomic gases were explored, as the original formulations of the closure models were developed strictly for gases composed of monatomic molecules. The development of a software library and computer code for solving relaxing hyperbolic systems in two spatial dimensions of the type arising from the closure models was also initiated. The software makes use of high-resolution upwind finite-volumes schemes, multi-stage point implicit time stepping, and automatic adaptive mesh refinement (AMR) to solve the governing conservation equations for the moment closures. The initial phase of the code development was completed and a numerical investigation of the solutions of the 10-moment closure model for the simple two-dimensional test cases mentioned above was initiated. Predictions of the 10-moment model were compared to available theoretical solutions and the results of direct-simulation Monte Carlo (DSMC) calculations. The first results of this study were presented at a meeting last year.

Groth, Clinton P. T.↗
NASA NTRS2005

Proceedings of the 2004 NASA/ONR Circulation Control Workshop, Part 2

This conference proceeding is comprised of papers that were presented at the NASA/ONR Circulation Control Workshop held 16-17 March 2004 at the Radisson-Hampton in Hampton, VA. Over two full days, 30 papers and 4 posters were presented with 110 scientists and engineers in attendance, representing 3 countries. As technological advances influence the efficiency and effectiveness of aerodynamic and hydrodynamic applications, designs, and operations, this workshop was intended to address the technologies, systems, challenges and successes specific to Coanda driven circulation control in aerodynamics and hydrodynamics. A major goal of this workshop was to determine the state-of-the-art in circulation control and to assess the future directions and applications for circulation control. The 2004 workshop addressed applications, experiments, computations, and theories related to circulation control, emphasizing fundamental physics, systems analysis, and applied research. The workshop consisted of single session oral presentations, posters, and written papers that are documented in this unclassified conference proceeding. The format of this written proceeding follows the agenda of the workshop. Each paper is followed with the presentation given at the workshop. the editors compiled brief summaries for each effort that is at the end of this proceeding. These summaries include the paper, oral presentation, and questions or comments that occurred during the workshop. The 2004 Circulation Control Workshop focused on applications including Naval vehicles (Surface and Underwater vehicles), Fixed Wing Aviation (general aviation, commercial, cargo, and business aircraft); V/STOL platforms (helicopters, military aircraft, tilt rotors); propulsion systems (propellers, jet engines, gas turbines), and ground vehicles (automotive, trucks, and other); wind turbines, and other nontraditional applications (e.g., vacuum cleaner, ceiling fan). As part of the CFD focus area of the 2004 CC Workshop, CFD practitioners were invited to compute a two-dimensional benchmark problem for which geometry, flow conditions, grids, and experimental data were available before the workshop. The purpose was to accumulate a database of simulations for a single problem using a range of CFD codes, turbulence models, and grid strategies so as to expand knowledge of model performance/requirements and guide simulation of practical CC configurations.

Jones, Gregory S.↗
NASA NTRS2002

Efficient Low Dissipative High Order Schemes for Multiscale MHD Flows

Accurate numerical simulations of complex multiscale compressible viscous flows, especially high speed turbulence combustion and acoustics, demand high order schemes with adaptive numerical dissipation controls. Standard high resolution shock-capturing methods are too dissipative to capture the small scales and/or long-time wave propagations without extreme grid refinements and small time steps. An integrated approach for the control of numerical dissipation in high order schemes for the compressible Euler and Navier-Stokes equations has been developed and verified by the authors and collaborators. These schemes are suitable for the problems in question. Basically, the scheme consists of sixth-order or higher non-dissipative spatial difference operators as the base scheme. To control the amount of numerical dissipation, multiresolution wavelets are used as sensors to adaptively limit the amount and to aid the selection and/or blending of the appropriate types of numerical dissipation to be used. Magnetohydrodynamics (MHD) waves play a key role in drag reduction in highly maneuverable high speed combat aircraft, in space weather forecasting, and in the understanding of the dynamics of the evolution of our solar system and the main sequence stars. Although there exist a few well-studied second and third-order high-resolution shock-capturing schemes for the MHD in the literature, these schemes are too diffusive and not practical for turbulence/combustion MHD flows. On the other hand, extension of higher than third-order high-resolution schemes to the MHD system of equations is not straightforward. Unlike the hydrodynamic equations, the inviscid MHD system is non-strictly hyperbolic with non-convex fluxes. The wave structures and shock types are different from their hydrodynamic counterparts. Many of the non-traditional hydrodynamic shocks are not fully understood. Consequently, reliable and highly accurate numerical schemes for multiscale MHD equations pose a great challenge to algorithm development. In addition, controlling the numerical error of the divergence free condition of the magnetic fields for high order methods has been a stumbling block. Lower order methods are not practical for the astrophysical problems in question. We propose to extend our hydrodynamics schemes to the MHD equations with several desired properties over commonly used MHD schemes.

Sjoegreen, Bjoern↗
NASA NTRS2008

How Very Massive Metal-Free Stars Start Cosmological Reionization

The initial conditions and relevant physics for the formation of the earliest galaxies are well specified in the concordance cosmology. Using ab initio cosmological Eulerian adaptive mesh refinement radiation hydrodynamical calculations, we discuss how very massive stars start the process of cosmological reionization. The models include nonequilibrium primordial gas chemistry and cooling processes and accurate radiation transport in the case B approximation using adaptively ray-traced photon packages, retaining the time derivative in the transport equation. Supernova feedback is modeled by thermal explosions triggered at parsec scales. All calculations resolve the local Jeans length by at least 16 grid cells at all times and as such cover a spatial dynamic range of approx.10(exp 6). These first sources of reionization are highly intermittent and anisotropic and first photoionize the small-scale voids surrounding the halos they form in, rather than the dense filaments they are embedded in. As the merging objects form larger, dwarf-sized galaxies, the escape fraction of UV radiation decreases and the H II regions only break out on some sides of the galaxies, making them even more anisotropic. In three cases, SN blast waves induce star formation in overdense regions that were formed earlier from ionization front instabilities. These stars form tens of parsecs away from the center of their parent DM halo. Approximately five ionizing photons are needed per sustained ionization when star formation in 10(exp 6) stellar Mass halos is dominant in the calculation. As the halos become larger than approx.10(exp 7) Stellar Mass, the ionizing photon escape fraction decreases, which in turn increases the number of photons per ionization to 15-50, in calculations with stellar feedback only. Radiative feedback decreases clumping factors by 25% when compared to simulations without star formation and increases the average temperature of ionized gas to values between 3000 and 10,000 K.

Wise, John H.↗
NASA NTRS2016

Liquid Oxygen/Liquid Methane Integrated Power and Propulsion

The proposed paper will cover ongoing work at the National Aeronautics and Space Administration (NASA) Johnson Space Center (JSC) on integrated power and propulsion for advanced human exploration. Specifically, it will present findings of the integrated design, testing, and operational challenges of a liquid oxygen / liquid methane (LOx/LCH4) propulsion brassboard and Solid Oxide Fuel Cell (SOFC) system. Human-Mars architectures point to an oxygen-methane economy utilizing common commodities, scavenged from the planetary atmosphere and soil via In-Situ Resource Utilization (ISRU), and common commodities across sub-systems. Due to the enormous mass gear-ratio required for human exploration beyond low-earth orbit, (for every 1 kg of payload landed on Mars, 226 kg will be required on Earth) increasing commonality between spacecraft subsystems such as power and propulsion can result in tremendous launch mass and volume savings. Historically, propulsion and fuel cell power subsystems have had little interaction outside of the generation (fuel cell) and consumption (propulsion) of electrical power. This was largely due to a mismatch in preferred commodities (hypergolics for propulsion; oxygen & hydrogen for fuel cells). Although this stove-piped approach benefits from simplicity in the design process, it means each subsystem has its own tanks, pressurization system, fluid feed system, etc. increasing overall spacecraft mass and volume. A liquid oxygen / liquid methane commodities architecture across propulsion and power subsystems would enable the use of common tankage and associated pressurization and commodity delivery hardware for both. Furthermore, a spacecraft utilizing integrated power and propulsion could use propellant residuals - propellant which could not be expelled from the tank near depletion due to hydrodynamic considerations caused by large flow demands of a rocket engine - to generate power after all propulsive maneuvers are complete thus utilizing previously wasted mass. Such is the case for human and robotic planetary landers. Although many potential benefits through integrated power & propulsion exist, integrated operations have yet to be successfully demonstrated and many challenges have already been identified the most obvious of which is the large temperature gradient. SOFC chemistry is exothermic with operating temperatures in excess of 1,000 K; however, any shared commodities will be undoubtedly stored at cryogenic temperatures (90-112 K) for mass efficiency reasons. Spacecraft packaging will drive these two subsystems in close proximity thus heat leak into the commodity tankage must be minimized and/or mitigated. Furthermore, commodities must be gasified prior to consumption by the SOFC. Excess heat generated by the SOFC could be used to perform this phase change; however, this has yet to be demonstrated. A further identified challenge is the ability of the SOFC to handle the sudden power spikes created by the propulsion system. A power accumulator (battery) will likely be necessary to handle these sudden demands while the SOFC thermally adjusts. JSC's current SOFC test system consists of a 1 kW fuel cell designed by Delphi. The fuel cell is currently undergoing characterization testing at the NASA JSC Energy Systems Test Area (ESTA) after which a Steam Methane Reformer (SMR) will be integrated and the combined system tested in closed-loop. The propulsion brassboard is approximately the size of what could be flown on a sounding rocket. It consists of one 100 lbf thrust "main" engine developed for NASA by Aerojet and two 10 lbf thrusters to simulate a reaction control system developed at NASA JSC. This system is also under development and initial testing at ESTA. After initial testing, combined testing will occur which will provide data on the fuel cell's ability to sufficiently handle the power spikes created by the propulsion system. These two systems will also be modeled using General-Use Nodal Network Solver (GUNNS) software. Once anchored with test data, this model will be used to extrapolate onto other firing profiles and used to size the power accumulator.

Banker, Brian↗
NASA NTRS1993

Large eddy simulation and direct numerical simulation of high speed turbulent reacting flows

The objective of this research is to make use of Large Eddy Simulation (LES) and Direct Numerical Simulation (DNS) for the computational analyses of high speed reacting flows. Our efforts in the first phase of this research conducted within the past three years have been directed in several issues pertaining to intricate physics of turbulent reacting flows. In our previous 5 semi-annual reports submitted to NASA LaRC, as well as several technical papers in archival journals, the results of our investigations have been fully described. In this progress report which is different in format as compared to our previous documents, we focus only on the issue of LES. The reason for doing so is that LES is the primary issue of interest to our Technical Monitor and that our other findings were needed to support the activities conducted under this prime issue. The outcomes of our related investigations, nevertheless, are included in the appendices accompanying this report. The relevance of the materials in these appendices are, therefore, discussed only briefly within the body of the report. Here, results are presented of a priori and a posterior analyses for validity assessments of assumed Probability Density Function (PDF) methods as potential subgrid scale (SGS) closures for LES of turbulent reacting flows. Simple non-premixed reacting systems involving an isothermal reaction of the type A + B yields Products under both chemical equilibrium and non-equilibrium conditions are considered. A priori analyses are conducted of a homogeneous box flow, and a spatially developing planar mixing layer to investigate the performance of the Pearson Family of PDF's as SGS models. A posteriori analyses are conducted of the mixing layer using a hybrid one-equation Smagorinsky/PDF SGS closure. The Smagorinsky closure augmented by the solution of the subgrid turbulent kinetic energy (TKE) equation is employed to account for hydrodynamic fluctuations, and the PDF is employed for modeling the effects of scalar fluctuations. The implementation of the model requires the knowledge of the local values of the first two SGS moments. These are provided by additional modeled transport equations. In both a priori and a posteriori analyses, the predicted results are appraised by comparison with subgrid averaged results generated by DNS. Based on these results, the paths to be followed in future investigations are identified.

Adumitroaie, V.↗