Search NASA⌕ Search

SEARCH · Search NASA

Results for “constitutive flow model”

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.

175 records · Page 10
All sourcesNASA NTRSDOE OSTI
NASA NTRS2024

Particle Impact Simulation and Ignition Prediction

An experimentally calibrated tool is needed to predict if a system is susceptible to failure by particle impact ignition (PI) based on use conditions, materials, and flow geometry. This tool will accelerate new components, evaluating existing hardware, and help disposition anomalies. - Conduct particle impact testing with in-situ diagnostics and complementary simulations on subset of key engineering materials (IN718, M400, 316L, 6061, Ti64, Zr) to develop a proof-of-concept predictive tool for assessing the risk of PI for idealized geometries (spherical particles) in realistic environments. - Assess particle/target interactions (coefficient of restitution, ignition, kindling) using instrumented particle impact rigs while systematically varying key parameters (materials, particle size, environment, target configuration). - Determine key field variables (temperature, strain, stress) in particle impacts using Multiphysics finite element and hydrocode simulations validated through comparison with experimental measurements and observations. - Synthesize experiments and simulations into constitutive models for PI that can be integrated with existing computational fluid dynamics (CFD) and Debris Transport Analysis (DTA) tools in future efforts.

Jonathan Tylka↗
NASA NTRS1987

Characteristics of cometary picked-up ions in a global model of Giacobini-Zinner

Energetic ions observed during the International Cometary Explorer (ICE) spacecraft flyby of comet Giacobini-Zinner provide information about both the constitution of comets and the plasma physical processes associated with their interaction with the solar wind. In this investigation the details of ion 'pickup,' in the limit where small-scale fluctuations in the plasma and magnetic field are neglected, are modeled by following the motion of a large number of initially cold, heavy (mass 18) ions in a global magnetohydrodynamic model of the local plasma and magnetic field. The results indicate how the background or macroscopic velocity and magnetic field structure of the comet can affect the average spatial and spectral characteristics of the observed cometary ions. These effects, which occur by virtue of forces associated with the compression and the curvature of the magnetic field in the presence of the stagnating plasma flow, can explain the double maxima in the time series of the energetic ion flux observed along the ICE trajectory.

Kimmel, C. D.↗
NASA NTRS1994

Notes on rotating turbulence

This work investigates the turbulent constitutive relation when turbulence is subjected to solid body rotation. Laws regarding spectra and asymptotic decay of rotating homogeneous turbulence were confirmed through large-eddy simulation (LES) computations. Rotating turbulent flows exist in many industrial, geophysical, and astrophysical applications. From Lagrangian analysis a relation between turbulent stress and strain in rotating homogeneous turbulence was inferred. This relation was used to derive the spectral energy flux and, ultimately, the energy spectrum form. If the rotation wavenumber k(sub Omega) lies in the inertial subrange, then for wavenumbers less than k(sub Omega) the turbulence motions are affected by rotation and the energy spectrum slope is modified. Energy decay laws inferred in other reports and the present results suggest a modification of the epsilon model equation and eddy viscosity in k-epsilon models.

Zeman, Otto↗
NASA NTRS1992

Preconditioning for stratospheric sudden warmings - Sensitivity studies with a numerical model

Results are presented of numerical experiments performed to determine what constitutes preconditioning for stratospheric sudden warmings. A number of 20-d runs were performed with a 3D primitive equation model. Runs varied in that they had different (zonally symmetric) initial fields and different lower boundary planetary wave fields, both based on NMC data for particular days. Initial and boundary data were used from the two minor and the sole major sudden warmings of January-February 1979 and some minor and the sole major sudden warmings of January-February 1989. To isolate the properties of the initial conditions from preconditioned cases and those with no warming, hybrid initial fields were created by combining data from two different days, one preconditioned and the other not. Runs with the hybrid initial conditions indicate that it is the lower stratospheric winds that are important in characterizing a preconditioned flow. It is suggested that it would be possible for a preconditioned state to exist for at least a short time without a warming occurring.

Smith, Anne K.↗
NASA NTRS1999

Fluid Physics in a Fluctuating Acceleration Environment

Our program of research aims at developing a stochastic description of the residual acceleration field onboard spacecraft (g-jitter) to describe in quantitative detail its effect on fluid motion. Our main premise is that such a statistical description is necessary in those cases in which the characteristic time scales of the process under investigation are long compared with the correlation time of g-jitter. Although a clear separation between time scales makes this approach feasible, there remain several difficulties of practical nature: (i), g-jitter time series are not statistically stationary but rather show definite dependences on factors such as active or rest crew periods; (ii), it is very difficult to extract reliably the low frequency range of the power spectrum of the acceleration field. This range controls the magnitude of diffusive processes; and (iii), models used to date are Gaussian, but there is evidence that large amplitude disturbances occur much more frequently than a Gaussian distribution would predict. The lack of stationarity does not constitute a severe limitation in practice, since the intensity of the stochastic components changes very slowly during space missions (perhaps over times of the order of hours). A separate analysis of large amplitude disturbances has not been undertaken yet, but it does not seem difficult a priori to devise models that may describe this range better than a Gaussian distribution. The effect of low frequency components, on the other hand, is more difficult to ascertain, partly due to the difficulty associated with measuring them, and partly because they may be indistinguishable from slowly changing averages. This latter effect is further complicated by the lack of statistical stationarity of the time series. Recent work has focused on the effect of stochastic modulation on the onset of oscillatory instabilities as an example of resonant interaction between the driving acceleration and normal modes of the system, and on cavity flow as an example of how an oscillatory response under periodic driving becomes diffusive if the forcing is random instead. This paper describes three different topics that illustrate behavior that is peculiar to a stochastic acceleration field. In the first case, we show that g-jitter can induce effective attractive or repulsive forces between a pair of spherical particles that are suspended in an incompressible fluid of different density provided that the momentum diffusion length is larger than the interparticle separation (as in the case in most colloidal suspensions). Second, a stochastic modulation of the control parameter in the vicinity of a pitchfork or supercritical bifurcation is known not to affect the location of the threshold. We show, however, that resonance between the modulation and linearly stable modes close to onset can lead to a shift in threshold. Finally, we discuss the classical problem of vorticity diffusion away from a plane boundary that is being vibrated along its own plane. Periodic motion with zero average vorticity production results in an exponential decay of the vorticity away from the boundary. Random vibration, on the other hand, results in power law decay away from the boundary even if vorticity production averages to zero.

Drolet, Francois↗
NASA NTRS2000

Numerical Simulation of Receptivity for a Transition Experiment

The cost of fuel to overcome turbulence induced viscous drag on a commercial airplane constitutes a significant fraction of the operating cost of an airline. Achieving laminar flow and maintaining it over a large portion of the wing can significantly reduce the viscous drag, and hence the cost. Design of such laminar-flow-control wings and their practical operation requires the ability to accurately and reliably predict the transition from laminar to turbulent flow. The transition process begins with the conversion of environmental and surface disturbances into the instability waves of the flow by a process called receptivity. The goal of the current research project has been to improve the prediction of transition through a better understanding of the physics of receptivity. The initial objective of this work was to investigate the specific stability and receptivity characteristics of a particular experimental investigation of boundary layer receptivity at NASA Langley. Some simulation results using direct solutions of the linearized Navier-Stokes equations which modeled this experiment where presented in the 1999 APS DFD meeting. However, based on these initial investigations, it became clear that to cover the vast receptivity parameter space required for a practical transition prediction tool, more efficient methods would be required. Thus, the focus of this research was shifted from modeling this particular experiment to formulating and developing new techniques that could efficiently yet accurately predict receptivity for a wide range of disturbance conditions.

Collis, S. Scott↗
NASA NTRS1997

Remote Sensing of Surficial Process Responses to Extreme Meteorological Events

Changes in the frequency and magnitude of extreme meteorological events are associated with changing environmental means. Such events are important in human affairs, and can also be investigated by orbital remote sensing. During the course of this project, we applied ERS-1, ERS-2, Radarsat, and an airborne sensor (AIRSAR-TOPSAR) to measure flood extents, flood water surface profiles, and flood depths. We established a World Wide Web site (the Dartmouth Flood Observatory) for publishing remote sensing-based maps of contemporary floods worldwide; this is also an online "active archive" that presently constitutes the only global compilation of extreme flood events. We prepared an article for EOS concerning SAR imaging of the Mississippi Valley flood; an article for the International Journal of Remote Sensing on measurement of a river flood wave using ERS-2, began work on an article (since completed and published) on the Flood Observatory for a Geoscience Information Society Proceedings volume, and presented lectures at several Geol. Soc. of America Natl. Meetings, an Assoc. of Amer. Geographers Natl. Meeting, and a Binghamton Geomorphology Symposium (all on SAR remote sensing of the Mississippi Valley flood). We expanded in-house modeling capabilities by installing the latest version of the Army Corps of Engineers RMA two-dimensional hydraulics software and BYU Engineering Graphics Lab's Surface Water Modeling System (finite elements based pre- and post-processors for RMA work) and also added watershed modeling software. We are presently comparing the results of the 2-d flow models with SAR image data. The grant also supported several important upgrades of pc-based remote sensing infrastructure at Dartmouth. During work on this grant, we collaborated with several workers at the U.S. Army Corps of Engineers, Remote Sensing/GIS laboratory (for flood inundation mapping and modeling; particularly of the Illinois River using the AIRSAR/TOPSAR/ERS-2 combined data), with Dr. Karen Prestegaard at the University of Maryland (geomorphological responses to the extreme 1993 flood along the Raccoon drainage in central Iowa), and with Mr Tim Scrom of the Albany National Weather Service River Forecast Center (initial planning for the use of Radarsat and ERS-2 for flood warning). The work thus initiated with this proposal is continuing.

Brakenridge, G. Robert↗
NASA NTRS2016

Application of a Full Reynolds Stress Model to High Lift Flows

A recently developed second-moment Reynolds stress model was applied to two challenging high-lift flows: (1) transonic flow over the ONERA M6 wing, and (2) subsonic flow over the DLR-F11 wing-body configuration from the second AIAA High Lift Prediction Workshop. In this study, the Reynolds stress model results were contrasted with those obtained from one- and two{equation turbulence models, and were found to be competitive in terms of the prediction of shock location and separation. For an ONERA M6 case, results from multiple codes, grids, and models were compared, with the Reynolds stress model tending to yield a slightly smaller shock-induced separation bubble near the wing tip than the simpler models, but all models were fairly close to the limited experimental surface pressure data. For a series of high-lift DLR{F11 cases, the range of results was more limited, but there was indication that the Reynolds stress model yielded less-separated results than the one-equation model near maximum lift. These less-separated results were similar to results from the one-equation model with a quadratic constitutive relation. Additional computations need to be performed before a more definitive assessment of the Reynolds stress model can be made.

Lee-Rausch, E. M.↗
NASA NTRS1999

The Isothermal Dendritic Growth Experiment (IDGE)

The Isothermal Dendritic Growth Experiment (IDGE) constituted a series of three NASA-supported microgravity experiments, all of which flew aboard the space shuttle, Columbia. This experimental space flight series was designed and operated to grow and record dendrite solidification in the absence of gravity-induced convective heat transfer, and thereby produce a wealth of benchmark-quality data for testing solidification scaling laws. The data and analysis performed on the dendritic growth speed and tip size in Succinontrie (SCN) demonstrates that although the theory yields predictions that are reasonably in agreement with experiment, there are significant discrepancies. However, some of these discrepancies can be explained by accurately describing the diffusion of heat. The key finding involves recognition that the actual three-dimensional shape of dendrites includes time-dependent side-branching and a tip region that is not a paraboloid of revolution. Thus, the role of heat transfer in dendritic growth is validated, with the caveat that a more realistic model of the dendrite then a paraboloid is needed to account for heat flow in an experimentally observed dendrite. We are currently conducting additional analysis to further confirm and demonstrate these conclusions. The data and analyses for the growth selection physics remain much less definitive. From the first flight, the data indicated that the selection parameter, sigma*, is not exactly a constant, but exhibits a slight dependence on the supercooling. Additional data from the second flight are being examined to investigate the selection of a unique dendrite speed, tip size and shape. The IDGE flight series is now complete. We are currently completing analyses and moving towards final data archiving. It is gratifying to see that the IDGE published results and archived data sets are being used actively by other scientists and engineers. In addition, we are also pleased to report that the techniques and IDGE hardware system that the authors developed with NASA, are being currently employed on both designated flight experiments, like EDSE, and on flight definition experiments, like TDSE.

Glicksman, Martin E.↗
NASA NTRS2000

Convective Melting of Particles in Flow Under Microgravity Conditions

Study of melting of dispersed or packed solid particles in a fluid under gravity and microgravity conditions provides benchmark information for many engineering applications such as material processing, environmental assessment and protection and space fire protection. During such processes, packed or dispersed solid particles are interacting with fluid flow at above-melting temperatures. By unmasking the buoyancy effects in coupled flow, phase change and heat transport phenomena, a better understanding of melting rate in non-thermal equilibrium, convective conditions can be studied. A series of flight experiments were conducted onboard the NASA KC-135 Microgravity Research Airplane. The Particle Melting in Plow (PMF) module was designed to allow flow through the initially packed ice particles at controlled temperature and velocity. To achieve this, a close-loop flow system was designed. Video images were taken to record the visualization of the melting process, from which a time variation of packed particle thickness distribution at different times can be obtained by the image analysis method. The fluid temperature distribution within the melting zone is measured by thermocouples. An infrared camera was mounted from the top of the test section to record the ice-water thermal images at a given location. The results from thermal images yield local temperature variation between melting solid and liquid and local Stephan number. Typical results for a number of cases are presented. The mathematical model, describing mass, energy, and momentum balance equations for the liquid and solid phases, is presented. It is found that melting rate is influenced mainly by the ratio of Reynolds number (based on the initial particle diameter) to the Froud number, and Stephan number. At the absence of gravity, Froud number approaches zero, Reynolds number and Stephan number become dominant factors governing the melting rate. The numerically determined results are compared with the experimental ones. It is found that the discrepancy between the predicted and measured melting rate is largely due to the inaccuracy in the constitutive equations for effective thermal physical properties, such as effective thermal conductivity and diffusivity, and transport properties, such as particle interaction coefficient, and local heat transfer coefficient of particles.

Stewart, C.↗
NASA NTRS2003

Minimizing Segregation during the Controlled Directional Solidification of Dendric Alloys

Gravity-driven convection induced in the liquid by density gradients of temperature or composition disrupts uniform dendritic growth during controlled directional solidification and promotes severe macrosegregation. The solute-rich region about the dendrite tip appears to play a pivotal role in channel initiation. Allen and Hunt referred to this region as an "initial transient" or dynamic region constituting steep concentration gradients. Experimental investigation also point to the role the tip region plays in developing microstructure. Hellawell and co-workers showed that flow-through dendritic channels could be effectively disrupted, and segregation minimized, during the gradient freezing of bulk castings by rotating the melt through a slight angle with respect to Earth's gravity vector. Adapting this principle to controlled directional solidification, it has been shown" that segregation in dendritic alloys can be minimized, and properties improved, by processing the sample near horizontal in conjunction with a slow axial rotation of the crucible. It is postulated that the observed microstructural uniformity arises by maintaining the developing solute field about the dendrite tip. Solute rejected during vertical directional solidification will rise or sink parallel to the primary dendrite arms during axial rotation setting the stage for accumulation, instabilities, and segregation. In contrast, during horizontal growth, the rejected solute will sink or rise perpendicular to the primary dendrite. Now, in the presence of a slight axial rotation, solute that was initially sinking (or rising) will find itself above (or below) its parent dendrite, i.e., still about the tip region. The following is intended to experimentally demonstrate the viability of this concept in coordination with a model that gives predictive insight regarding solute distribution about growing dendrites. Alloys based on the lead-tin eutectic system were used in this study. The system is well characterized, the constituent metals are available in a very pure form, and the thermophysical properties are well known. During solidification of hypoeutectic alloys, e.g., 55 wt pct Pb, the primary dendrites reject the less dense tin, and for the hypereutectic alloys, e.g., 75 wt pct Sn, the primary dendrites reject denser lead. Alloys were prepared by melting appropriate amounts of lead and tin in a glass crucible after which the homogeneous liquid was sucked directly into 5-mm i.d. glass tubes. The sample tube, containing approximately 30 cm of alloy, was then mechanically driven into the directional solidification furnace assembly and positioned such that approx. 20 cm of the sample was remelted. Subsequently, directional solidification was initiated by withdrawing the sample through a water-cooled jacket at a constant growth velocity of 2 ,microns/s. After 5 to 6 cm of growth, the sample was quickly removed from the furnace and quenched in a water bath to preserve the solid-liquid interface. Samples were directionally solidified vertically upward, nearly horizontally, and some in conjunction with an applied axial rotation of the crucible. Temperature gradients at the solid-liquid interface were measured with an in-siru K-type thermocouple. Solidified samples were cut perpendicular and parallel to the growth direction and conventionally prepared for microscopic examination.

Grugel, Richard N.↗
NASA NTRS2001

Combustion of Metals in Carbon Dioxide and Reduced-Gravity Environments

Ongoing exploration and future mission2001110444 s to Mars have given impetus to research on the use of natural resources of the planet. Since carbon dioxide (CO2) constitutes approximately 95% of the Mars atmosphere and since it reacts directly and vigorously with several metals, this investigation focuses on metal-CO2 reactions as a possible combination for rocket-propellant production and energy generation. Magnesium (Mg) has been initially selected as the metal fuel owing to its low ignition temperature and high specific impulse and burning rate in CO2. Our studies in this field started with low gravity (g) combustion tests of Mg in O2, CO2, and CO. Reduced gravity provided a clear picture of the burning phenomena by eliminating the intrusive buoyant flows in high-temperature metal reactions and by removing the destructive effect of gravity on the shape of molten metal samples. Suspended cylindrical metal samples of 2, 3, and 4-mm in diameter and length were radiatively ignited in low-g to generate free-floating samples exhibiting a spherically symmetric flame with increasing metal-oxide accumulation in an outer shell. For the Mg-CO2 combination, burning times twice as long as in normal-g and five times longer than in Mg-O2 flames were observed, revealing a diffusion-controlled reaction. The burning time is proportional to the square of the sample diameter. In tests conducted with pure CO, combustion was not possible without constant heating of the sample due to the formation of a thick carbon-containing coating around the Mg sample generated by surface reactions. The following work presents two new studies that attempt to explain some of the low-g experimental observations. First, a simplified one-dimensional, quasi-steady numerical model is developed to obtain temperature, species concentrations, and burning rates of the spherically symmetric diffusion flame around the Mg sample burning in O2 and CO2. Second, a Planar Laser Induced Fluorescence (PLIF) technique is implemented to provide spatially resolved measurements of magnesium oxide (MgO) in the reaction zone of Mg samples burning in O2 and CO2. These experiments reveal fundamental differences between the two combustion systems.

Branch, M. C.↗
NASA NTRS1999

A Study of Nucleate Boiling with Forced Convection in Microgravity

The ultimate objective of basic studies of flow boiling in microgravity is to improve the understanding of the processes involved, as manifested by the ability to predict its behavior. This is not yet the case for boiling heat transfer even in earth gravity, despite the considerable research activity over the past 30 years. The elements that constitute the nucleate boiling process - nucleation, growth, motion, and collapse of the vapor bubbles (if the bulk liquid is subcooled) - are common to both pool and flow boiling. It is well known that the imposition of bulk liquid motion affects the vapor bubble behavior relative to pool boiling, but does not appear to significantly influence the heat transfer. Indeed, it has been recommended in the past that empirical correlations or experimental data of pool boiling be used for design purposes with forced convection nucleate boiling. It is anticipated that such will most certainly not be possible for boiling in microgravity, based on observations made with pool boiling in microgravity. In earth gravity buoyancy will act to remove the vapor bubbles from the vicinity of the heater surface regardless of how much the imposed bulk velocity is reduced, depending, of course, on the geometry of the system. Vapor bubbles have been observed to dramatically increase in size in pool boiling in microgravity, and the heat flux at which dryout took place was reduced considerably below what is generally termed the critical heat flux (CHF) in earth gravity, depending on the bulk liquid subcooling. However, at heat flux levels below dryout, the nucleate pool boiling process was enhanced considerably over that in earth gravity, in spite of the large vapor bubbles formed in microgravity and perhaps as a consequence. These large vapor bubbles tended to remain in the vicinity of the heater surface, and the enhanced heat transfer appeared to be associated with the presence of what variously has been referred to as a liquid microlayer between the bubble and the heater surface. The enhancement of the boiling process with low velocities in earth gravity for those orientations producing the formation of a liquid macrolayer described above, accompanied by "sliding" vapor bubbles, has been demonstrated. The enhancement was presented as a function of orientation, subcooling, and heated length, while a criterion for the heat transfer for mixed natural/forced convection nucleate boiling was given previously. A major unknown in the prediction and application of flow boiling heat transfer in microgravity is the upper limit of the heat flux for the onset of dryout (or critical heat flux - CHF), for given conditions of fluid-heater surfaces, including geometry, system pressure and bulk liquid subcooling. It is clearly understood that the behavior in microgravity will be no different than on earth with sufficiently high flow velocities, and would require no space experimentation. However, the boundary at which this takes place is still an unknown. Previous results of CHF measurements were presented for low velocity flow boiling at various orientations in earth gravity as a function of flow velocity and bulk liquid subcooling, along with preliminary measurements of bubble residence times on a flat heater surface. This showed promise as a parameter to be used in modeling the CHF, both in earth gravity and in microgravity. The objective of the work here is to draw attention to and show results of current modeling efforts for the CHF, with low velocities in earth gravity at different orientations and subcoolings. Many geometrical possibilities for a heater surface exist in flowing boiling, with boiling on the inner and outer surfaces of tubes perhaps being the most common. If the vapor bubble residence time on and departure size from the heater surface bear a relationship to the CHF, as results to be given indicate, it is important that visualization of and access to vapor bubble growth be conveniently available for research purposes. In addition, it is desirable to reduce the number of variables as much as possible in a fundamental study. These considerations dictated the use of a flat heater surface, which is rectangular in shape, 1.91 cm by 3.81 cm (0.75 x 1.5 inches), consisting either of a 400 Angstrom thick semi-transparent gold film sputtered on a quartz substrate which serves simultaneously as a heater and a resistance thermometer, or a copper substrate of the same size. The heater substrate is a disc which can be rotated so that the heated length in the flow direction can be changed from 1.91 to 3.81 cm (0.75 to 1.5 inches). The fluid is R-113, and the velocities can be varied between 0.5 cm/s and 60 cm/s. For a sufficiently low velocity the CHF can be modeled reasonably well at various orientations by the correlation for pool boiling corrected for the influence of bulk liquid subcooling, multiplied by the square root of q, the angle relative to horizontal. This arises from equating buoyancy and drag forces in the inverted positions where the vapor bubbles are held against the heater surface as they slide. A distortion of the measurements relative to pool boiling occurs as the flow velocity increases. In modeling this effect at different levels of subcooling it appeared appropriate to estimate the volumetric rate of vapor generation, using measurements of bubble frequency (or residence time), void fraction and average bubble boundary layer thickness. These were determined with the use of a platinum hot wire probe 0.025 mm in diameter by 1.3 mm long, applying a constant current to distinguish between contact with liquid or vapor. Two-dimensional spatial variations are obtained with a special mechanism to resolve displacements in increments of 0.025 mm. From such measurements it was determined that the fraction of the surface heat transfer resulting in evaporation varies inversely with the subcooling correction factor for the CHF. The measured inverse bubble residence time is normalized relative to that predicted for an infinite horizontal flat plate at the CHF, and is correlated well with the CHF normalized relative to that for pool boiling, for various orientation angles and subcooling levels. This correspondence is then combined with a normalizing factor for the energy flux leaving the heater surface at the CHF and the computed bubble radius at departure, determined from the balance between the outward velocity of the interface due to evaporation and the buoyance induced velocity of the center of mass of the bubble. The product of the CHF and the corresponding residence time was determined to be a constant for all orientations at a given bulk flow velocity and liquid subcooling, and must be determined empirically for each velocity and subcooling at present. It then becomes possible to predict the CHF for the different orientations, velocities, and subcoolings. These are compared with normalized measurements of the CHF for velocities ranging from 4 cm/s to 55 cm/s, subcoolings from 2.8 to 22.2 K, over orientations angles of 360 degrees.

Merte, Herman, Jr.↗