SEARCH · Engineering Papers
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.
A cooling apparatus includes a container filled with a quantity of coolant fluid initially cooled to a solid phase, a cooling loop disposed between a heat load and the container, a pump for circulating a quantity of the same type of coolant fluid in a liquid phase through the cooling loop, and a pair of couplings for communicating the liquid phase coolant fluid into the container in a direct interface with the solid phase coolant fluid.
Despite essentially identical crystallography and equilibrium structuring of water, nanoscopic channels composed of hexagonal boron nitride and graphite exhibit an order-of-magnitude difference in fluid slip. We investigate this difference using molecular dynamics simulations, demonstrating that its origin is in the distinct chemistries of the two materials. In particular, the presence of polar bonds in hexagonal boron nitride, absent in graphite, leads to Coulombic interactions between the polar water molecules and the wall. We demonstrate that this interaction is manifested in a large typical lateral force experienced by a layer of oriented hydrogen atoms in the vicinity of the wall, leading to the enhanced friction in hexagonal boron nitride. The fluid adhesion to the wall is dominated by dispersive forces in both materials, leading to similar wettabilities. Finally, our results rationalize recent observations that the difference in frictional characteristics of graphite and hexagonal boron nitride cannot be explained on the basis of the minor differences in their wettabilities.
Nuclear Electric Propulsion (NEP) is an in-space propulsion technology capable of enabling opposition and conjunction class crewed Mars missions. NEP subsystems include the reactor for heat generation, a power conversion system (PCS), power management and distribution, electric propulsion system, and heat rejection system. Specific mass, or α (kg/kWe), is a key performance parameter (KPP) of the propulsion system which is directly scalable with the performance and mass predictions for each individual component. To inform technology maturation planning activities, full system and component level parametric modeling is ongoing to explore the design trade space and illustrate the effect of subsystem design choices on the system KPPs. In this study, scaling of high-assay, low-enriched uranium reactor designs is assessed through coupled reactor physics and thermal hydraulics analyses. Scaling analyses evaluate the impact of system performance parameters (power level, interface temperatures) on mass for direct gas cooled, pumped liquid metal, and passively-cooled heat pipe reactor concepts. Each concept requires specific geometries, fluids, and power conversion interface conditions (temperature, pressure, flow rate) to meet desired performance and mass. The reactor assembly includes the active core (fuel, moderator, cladding, working fluid), axial and radial neutron reflectors, control drums, structural support / pressure vessel, and external radiation shielding. Each of these components are parametrically sized based on performance parameters for a megawatt-class power cycle. Results of this scaling analysis increase NEP propulsion system modeling fidelity and ultimately aim to support concept down-selection along with related technology development planning. The reactor and shield α are a function of several PCS and heat rejection system design choices, and reactor scaling with these parameters must be considered to enable an informed decision on an optimal reactor geometry and working fluid combination.
Vapor-liquid interfaces drive many important phenomena in cryogenic fluid management, including heat transfer, evaporation, and capillary flow. Design of cryogenic fluid systems, such as propellant storage, requires accurate predictions of fluid behavior, including evaporation rates. Many models have been proposed for heat and mass transfer at vapor liquid interfaces, but the accuracy of these models in the context of cryogenic fluids has not been performed. We use molecular dynamics simulations, which allow for nanometer scale resolution of fluid phenomena, to evaluate the accuracy of a variety of vapor-liquid boundary conditions at evaporating and condensing interfaces. We find that an anisotropic temperature distribution is a critical ingredient for accurate prediction of intensive evaporation and condensation.
The International Space Station (ISS) program is nearing an assembly complete configuration with the addition of the final resource node module in early 2010. The Node 3 module will provide critical functionality in support of permanent long duration crews aboard ISS. The new module will permanently house the regenerative Environment Control and Life Support Systems (ECLSS) and will also provide important habitability functions such as waste management and exercise facilities. The ISS program has selected the Port side of the Node 1 "Unity" module as the permanent location for Node 3 which will necessitate architecture changes to provide the required interfaces. The USOS ECLSS fluid and ventilation systems, Internal Thermal Control Systems, and Avionics Systems require significant modifications in order to support Node 3 interfaces at the Node 1 Port location since it was not initially designed for that configuration. This paper outlines the design, development, certification, and implementation of these changes in support of ISS assembly complete.
Controlling interfaces of phase-separating fluid mixtures is key to the creation of diverse functional soft materials. Traditionally, this is accomplished with surface-modifying chemical agents. Using experiment and theory, we studied how mechanical activity shapes soft interfaces that separate an active and a passive fluid. Chaotic flows in the active fluid give rise to giant interfacial fluctuations and noninertial propagating active waves. At high activities, stresses disrupt interface continuity and drive droplet generation, producing an emulsion-like active state composed of finite-sized droplets. When in contact with a solid boundary, active interfaces exhibit nonequilibrium wetting transitions, in which the fluid climbs the wall against gravity. These results demonstrate the promise of mechanically driven interfaces for creating a new class of soft active matter.
After an elementary introduction to solidification, recent research that attempts to understand the interaction of fluid flow with a crystal-melt interface under gravitational influence using linear stability analysis is reviewed. Recent numerical calculations of fluid flow and interface morphology during alloy solidification are briefly discussed.
Large amplitude symmetric and asymmetric irrotational motion of inviscid incompressible fluid with liquid-vapor interface in accelerating cylindrical container of revolution
This paper documents a model validation effort for the prediction of cryogenic propellant tank selfpressurization physics using the Loci/STREAM sharp interface CFD model developed at NASA MSFC. Cryogenic fluid management (CFM) applications for long duration spaceflight missions require accurate modeling of heat transfer and thermodynamics. Therefore, detailed validation of those physics is important to ensure simulation accuracy. In the past, the MSFC Propulsion Fluid Dynamics Branch has validated and used the Loci/STREAM Volume of Fluid (VOF) module extensively for propellant slosh applications, which are dominated by fluid motion. For cases where the interface gas/liquid interface is essentially static, the Loci/STREAM sharp interface was developed for more rapid and efficient CFM analysis. In the present study, Loci/STREAM sharp interface simulations were compared with the NASA K-site self-pressurization experiment of a flight weight, partially full, LH2 tank in normal gravity. This study investigated several parametric effects and made a number of important observations on the performance of the Loci/STREAM sharp interface method: - First, a homogeneous model was derived to verify the CFD model. The homogeneous model assumes both gas and liquid phases are uniform and are at saturation temperature. The predicted pressure rise rates at 3.5 W/m2 and 2.0 W/m2 with a 49% fill level are all lower than the experimental measurement. This is expected due to the observed temperature stratification near the gas-liquid interface during testing. - The k-w and k-e turbulence models predict comparable pressurization rates, which are much lower than the experiment measurement. On the other hand, the pressure values and pressurization rates from both models are almost the same as those of the homogeneous model. The current turbulence model is too diffusive and is not capable of capturing the thermal stratification in the liquid and ullage. - When using the laminar flow model, the pressurization rate is found to be very close to the experimental value when a quasi-steady flow field is allowed to develop by continuously venting the tank prior to pressurization as was done during testing. On the other hand, without initial venting, the rate is not constant. Additionally, the development time to reach a quasi-steady flow field differed from the experiment due to the lack of turbulent mixing away from the gas-liquid interface that is expected in reality. - Using the laminar flow model and appropriate initial conditioning, the pressurization rate is within 15% of the experimental value for 49% fill level at 3.5 W/m 2 and is within 18% for 2.0 W/m 2 heat loading. - Using the laminar flow model and appropriate initial conditioning, the pressurization rate is within 18% of the experimental value for a 29% fill level and within 15% for an 83% fill level.
Many important engineering and scientific applications such as cement slurries, foams, crude oil, and granular avalanches involve the concept of yield stress. Therefore, modeling yield stress fluids in different flow configurations, including the accurate prediction of the yield surface, is important. In this paper, we present a computational model based on the finite element method to study the flow of yield stress fluids in a thin mold and compare the results with data from flow visualization experiments. Here, we use the level set method to describe the interface between the filling fluid and air. We use polypropylene glycol as a model Newtonian fluid and Carbopol for the model yield stress fluid, as the Carbopol solution demonstrates yielding without thixotropy. To describe the yielding and shear-thinning behavior, we use a generalized Newtonian constitutive equation with a Bingham–Carreau–Yasuda form. We compare the results obtained from the mold filling experiments with the results from the three-dimensional (3D) model and from a reduced-order Hele-Shaw (HS) model that is two-dimensional, including the effect of shear-thinning along the thin direction only approximately. We show that both the 3D and the HS model can capture the experimental meniscus shape reasonably well for all the fluids considered at three different flow rates. This indicates that the shape evolution is insensitive to the dimensionality of the model. However, the viscosity and yield surfaces predicted by the 3D and HS models are different. The HS model underestimates the high viscosity and unyielded regions compared to the estimation by the 3D model.
More accurate characterization and prediction of the in-situ distribution of fracturing fluid in fractured reservoirs are needed for enhancing well productivity. In this study, an implicit-sequentially coupled flow/geomechanics simulator incorporating an efficient discrete fracture model is developed to model fluid distribution and recovery performance of ultra-low permeability gas reservoirs. The finite-volume and finite-element methods are used for space discretization of the flow and geomechanics equations, respectively, while the backward Euler method is employed for time discretization. The flow and geomechanics equations are solved sequentially based on fixed-stress splitting. An efficient discrete-fracture model is used to explicitly model the fractured system. Flexible unstructured gridding is employed to model arbitrarily-oriented fractures. The interrelations among pore volume, permeability and geomechanical conditions are considered dynamically using two-way coupled flow and geomechanics computations. The geometry of fracture (networks) due to hydraulic fracturing has significant impacts on the fracturing fluid recovery efficiency and ensuing fluid distribution. Under the same injection volume, the fracturing fluid recovery is higher when the fracture geometry is planar. Fluid recovery is relatively lower whenever natural fractures are activated during fracturing treatments; flowback time is also shortened when complex fracture network with enlarged fracture interface is present. Fracturing fluid in hydraulic fractures may leak off into the natural fractures and subsequently imbibes into the surrounding matrix due to capillarity effects. The fracturing fluid recovery and in-situ fluid distribution are sensitive to the shut-in duration and fracture closure behavior. This study analyzes the coupled flow-geomechanical responses of fractured gas reservoirs during the post-fracturing periods. Understanding the fate of the fracturing fluid can provide insights on, to some extent, the stimulated fracture volume, size of the water invasion zone, and efficiency of the fracturing design. The simulation predictions can also provide more accurate initial reservoir conditions (e.g. distributions of different phases and pressure) for long-term well performance estimation.
The flow experiment consisted of an oil tunnel, 76 x 76 mm in cross-section, packed with lucite spheres. The index of refraction of the working fluid and the spheres were matched such that the physical spheres invisible to the eye and camera. By seeding the oil and illuminating the packed bed with planar laser light sheet, aligned in the direction of the bulk flow, the system fluid dynamics becomes visible and the 2-D projection was recorded at right angles to the bulk flow. The planar light sheet was traversed from one side of the tunnel to the other providing a simulated 3-D image of the entire flow field. The boundary interface between the working fluid and the sphere rendered the sphere black permitting visualization of the exact locations of the circular interfaces in both the axial and transverse directions with direct visualization of the complex interstitial spaces between the spheres within the bed. Flows were observed near the surfaces of a plane and set of spheres as well as minor circles that appear with great circles and not always uniformly ordered. In addition to visualizing a very complex flow field, it was observed that flow channeling in the direction of the bulk flow occurs between sets of adjacent spheres. Still photographs and video recordings illustrating the flow phenomena will be presented.
We propose a non-local model for surface tension obtained in the form of an integral of a molecular-force-like function with support 3.5e added to the Navier-Stokes momentum conservation equation. We demonstrate analytically and numerically that with the non-local model interfaces with a radius of curvature larger than the support length behave macroscopically and microscopically, otherwise. For static droplets, the pressure difference Pe,in - Pe,out satisfies the Young-Laplace law for droplet radius greater than 3.5e and other- wise deviates from the Young-Laplace law. The latter indicates that the surface tension in the proposed model decreases with decreasing radius of curvature, which agrees with molecular dynamics and experimental studies of nanodroplets. Using the non-local model we perform numerical simulations of droplets under dynamic conditions, including a rising droplet, a droplet in shear flow, and two colliding droplets in shear flow, and compare results with a standard Navier-Stokes model subject to the Young-Laplace boundary condition at the fluid-fluid interface implemented via the Conservative Level Set (CLS) method. We find good agreement with existing numerical methods and analytical results for a rising macroscopic droplet and a droplet in a shear flow. For colliding droplets in shear flow, the non- local model converges (with respect to the grid size) to the correct behavior, including sliding, coalescence, and merging and breaking of two droplets depending on the capillary number. In contrast, we find that the results of the CLS model are highly grid-size dependent.
When two fluid phases coexist near their mutual critical point, one fluid phase forms a layer which intrudes between the other fluid phase and any third phase that happens to be present. As the two fluid phases are taken away from their critical point, a phase transition often occurs such that the intruding layer vanishes. Recent evidence is presented that the transition from the complete wetting (intruding layer) configuration to the incomplete wetting (three-phase contact) configuration is a first-order phase transition. The thickness of the intruding layers has been measured for diverse systems. These data are not fully understood. The intruding layers are extreme examples of multilayer adsorption. Certain theories predict that a 'prewetting' transition from high adsorption to low adsorption is associated with the transition from complete wetting to incomplete wetting. Present search for this transition is inconclusive.
The Microgravity Transport Phenomena Experiment (MTPE) is a fluids experiment supported by the Fundamentals in Biotechnology program in association with the Human Exploration and Development of Space (BEDS) initiative. The MTP Experiment will investigate fluid transport phenomena both in ground based experiments and in the microgravity environment. Many fluid transport processes are affected by gravity. Osmotic flux kinetics in planar membrane systems have been shown to be influenced by gravimetric orientation, either through convective mixing caused by unstably stratified fluid layers, or through a stable fluid boundary layer structure that forms in association with the membrane. Coupled transport phenomena also show gravity related effects. Coefficients associated with coupled transport processes are defined in terms of a steady state condition. Buoyancy (gravity) driven convection interferes with the attainment of steady state, and the measurement of coupled processes. The MTP Experiment measures the kinetics of molecular migration that occurs in fluids, in response to the application of various driving potentials. Three separate driving potentials may be applied to the MTP Experiment fluids, either singly or in combination. The driving potentials include chemical potential, thermal potential, and electrical potential. Two separate fluid arrangements are used to study membrane mediated and bulk fluid transport phenomena. Transport processes of interest in membrane mediated systems include diffusion, osmosis, and streaming potential. Bulk fluid processes of interest include coupled phenomena such as the Soret Effect, Dufour Effect, Donnan Effect, and thermal diffusion potential. MTP Experiments are performed in the Microgravity Transport Apparatus (MTA), an instrument that has been developed specifically for precision measurement of transport processes. Experiment fluids are contained within the MTA fluid cells, designed to create a one dimensional flow geometry of constant cross sectional area, and to facilitate fluid filling and draining operations in microgravity. The fluid cells may be used singly for bulk solutions, or in a Stokes diaphragm configuration to investigate membrane mediated phenomena. Thermal and electrical driving potentials are applied to the experiment fluids through boundary plates located at the ends of the fluid cells. In the ground based instrument, two constant temperature baths circulate through reservoirs adjacent to the boundary plates, and establish the thermal environment within the fluid cells. The boundary plates also serve as electrodes for measurement and application of electrical potentials. The Fluid Manipulation System associated with the MTA is a computer controlled system that enables storage and transfer of experiment fluids during on orbit operations. The system is used to automatically initiate experiments and manipulate fluids by orchestrating pump and valve operations through scripted sequences. Unique technologies are incorporated in the MTA for measurement of fluid properties. Volumetric Flow Sensors have been developed for precision measurement of total fluid volume contained within the fluid cells over time. This data is most useful for measuring the kinetics of osmosis, where fluid is transported from one fluid cell to another through a semipermeable membrane. The MicroSensor Array has been designed to perform in situ measurement of several important fluid parameters, providing simultaneous measurement of solution composition at multiple locations within the experiment fluids. Micromachined sensors and interface electronics have been developed to measure temperature, electrical conductivity, pH, cation activity, and anion activity. The Profile Refractometer uses a laser optical system to directly image the fluid Index of Refraction profile that exists along the MTA fluid cell axis. A video system acquires images of the RI profile over time, and records the transport kinetics that occur upon application of chemical, thermal, or electrical driving potentials. Image processing algorithms have been developed to analyze the refractometer images on a pixel by pixel basis, calibrating and scaling the measured Index of Refraction profile to correlated solution properties of interest such as density, concentration, and temperature. Additional software has been developed to compile the processed images into a three dimensional matrix that contains fluid composition data as a function of experiment time and position in the fluid cell. These data are combined with data from the other sensor systems, and analyzed in the context of transport coefficients associated with the various transport phenomena. Analysis protocols have been developed to measure the transient kinetics, and steady state distribution of fluid components that occur in response to the applied driving potentials. The results are expressed in terms of effective transport coefficients. Experiments have been performed using a variety of solutes, and results generated are that are in agreement with published transport coefficient values.
Thermal buoyancy and surface tension driven convection is numerically investigated in a system with two, immiscible fluids. The geometry investigated, has an open cavity configuration with the lighter fluid situated on top of the heavier fluid, forming a stable layered system. The upper fluid meniscus and the interface are assumed to be flat and undeformable in the calculations. The governing equations and boundary/interface conditions are solved by a control volume based finite difference scheme for two pairs of immiscible fluids; the water-hexadecane system and a so called generic system. The steady state calculations show, that dramatically different flows are predicted when the interfacial tension effects are included or excluded from the system model. These differences are particularly appreciable in surface tension dominated flows, that are typical of microgravity situations. Complex flow patterns, with induced secondary flows are noticed in both the fluids. In general, the overall system heat transfer is found to increase with increases in the thermal buoyancy and surface tension effects, but the behavior of the system flow and thermal fields is not easy to characterize, when different combinations of these forces are considered.
The field-induced force density on a magnetic fluid is proportional to the magnetic susceptibility times the gradient of the magnetic field squared. The direction of the force is towards increasing magnetic field (positive gradient). Applying a magnetic field to a magnetic fluid will result in a force from all directions towards the location of peak field. Since the magnetic field is conservative and there are no magnetic monopoles, the net field-induced force on any fluid of constant susceptibly will be zero. The only manner to obtain a nonzero net field-induced force is to vary the susceptibility of the fluid. At the gas/liquid interface of liquid oxygen, the susceptibility varies drastically, and the exploitation of the resultant large net forces. An alternative method of varying the magnetic susceptibility is to vary the temperature of the fluid. The magnetic susceptibility of paramagnetic liquid oxygen obeys the Curie-Weiss law: it is inversely proportional to temperature. By applying a temperature gradient in the presence of a symmetric magnetic field, a nonzero net force results. Much of the theory of the so-called Magnetocaloric Effect has previously been developed for and applied to ferromagnetic fluids, or ferrofluids, but is readily applied to paramagnetic liquid oxygen.
Almost one hundred years ago Korteweg published a theory of how stresses could be induced in miscible fluids by concentration gradients, causing phenomena that would appear to be the same as with immiscible fluids. Miscible fluids could manifest a transient or effective interfacial tension (EIT). To this day, there has been no definitive experiment to confirm Korteweg's model but numerous fascinating and suggestive experiments have been reported. The goal of TIPMPS is to answer the question: Can concentration and temperature gradients in miscible materials induce stresses that cause convection? Many polymer processes involving miscible monomer and polymer systems could be affected by fluid flow and so this work could help understand miscible polymer processing, not only in microgravity, but also on earth. Demonstrating the existence of this phenomenon in miscible fluids will open up a new area of study for materials science. The science objectives of TIPMPS are: (1) Determine if convection can be induced by variation of the width of a miscible interface; (2) Determine if convection can be induced by variation of temperature along a miscible interface; (3) Determine if convection can be induced by variation of conversion along a miscible interface An interface between two miscible fluids can best be created via a spatially-selective photopolymerization of dodecyl acrylate with a photoinitiator, which allows the creation of precise and accurate concentration gradients between polymer and monomer. Optical techniques will be used to measure the refractive index variation caused by the resultant temperature and concentration fields. The viscosity of the polymer will be measured from the increase in the fluorescence of pyrene. Because the large concentration and temperature gradients cause buoyancy-driven convection that prevents the observation of the predicted flows, the experiment must be done in microgravity. In this report, we will consider our efforts to estimate the square gradient parameter, k, and our use of the estimates in modeling of the planned TIPMPS experiments. We developed a model consisting of the heat and diffusion equations with convective terms and of the Navier-Stokes equations with an additional volume force written in the form of the Korteweg stresses arising from nonlocal interaction in the fluid. The fluid's viscosity dependence on polymer conversion and temperature was taken from measurements of poly(dodecyl acrylate). Numerical modeling demonstrated that significant flows would arise for conditions corresponding to the planned experiments.