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.
Fluid mechanics experiments to investigate methods for reducing mixing between confined coaxial flows in cylindrical chambers for application to open-cycle gaseous-core nuclear rockets
A model formulated in terms of both conservation and kinematic equations for phases and interfaces in two-fluid-phase flow in a porous medium system is summarized. Macroscale kinematic equations are derived as extensions of averaging theorems and do not rely on conservation principles. Models based on both conservation and kinematic equations can describe multiphase flow with varying fidelity.When only phase-based equations are considered, a model similar in form to the traditional model for two-fluid-phase flow results. When interface conservation and kinematic equations are also included, a novel formulation results that naturally includes evolution equations that express dynamic changes in fluid saturations, pressures, the capillary pressure, and the fluid-fluid interfacial area density in a two-fluid-system. This dynamic equation set is unique to this work, and the importance of the modeled physics is shown through both microfluidic experiments and high-resolution lattice Boltzmann simulations. The validation work shows that the relaxation of interface distribution and shape toward an equilibrium state is a slow process relativeto the time scale typically allowed for a system to approach an apparent equilibrium state based upon observations of fluid saturations and external pressure measurements. Consequently, most pressure-saturation data intended to denote an equilibrium state are likely a sampling from a dynamic system under-going changes of interfacial curvatures that are not typically monitored. The results confirm the importance of kinematic analysis in combination with conservation equations for faithful modeling of system physics.
Wetting films can develop in the corners of angular pores under strong wetting conditions. Modeling the dynamics of corner film remains elusive using direct numerical simulations because of the significant scale difference between main meniscus and corner film flow. In this paper, the modified interacting capillary bundle model (ICB), developed in our previous work to describe accurately corner film dynamics in a single square tube, is incorporated into a single-pressure dynamic pore network model (DPNM) to simulate imbibition in strongly wetting porous media with corner film flow. The traditional pore network is decomposed into several layers of interacting subpore networks where the 0th layer of subpore network simulates the main meniscus flow and higher layers the corner film flow. The fluid flow between different layers is captured by interlayer throats. In addition, the snap-off mechanism caused by the thickening of wetting corner film is considered. The accuracy of the developed model is validated for four cases: spontaneous imbibition in a single square tube, wetting fluid redistribution through corner films under a capillary pressure difference, snap off in a narrow throat connecting two large pores, and imbibition dynamics in a real microfluidic porous geometry. The validated model is then used to simulate both spontaneous and controlled imbibition in a pore network with random pore size distribution. Finally, the interaction between corner film and main meniscus flow in porous media is analyzed from a pore-scale perspective.
An apparatus for energy extraction from fluid flow including an assembly including a plenum. The assembly further includes an aperture extending from an exterior surface to the plenum to allow flow therethrough. The apparatus further includes a channel including an inlet and an outlet in fluid communication with the plenum. The apparatus yet further includes an energy extraction device. The assembly is configured to create a pressure differential between the plenum and the inlet of the channel. The pressure differential causes fluid flow from the inlet of the channel to the plenum. The energy extraction device is configured to extract energy from the fluid flow. The apparatus additionally includes a control system configured to modify the pressure differential to control the fluid flow between the inlet of the channel and the plenum based on a characteristic of an exterior environment.
A heat exchanger includes a shell housing a plurality of tubes and defining an exhaust fluid flow path within a first volume enclosed by the shell. The outer surfaces of the plurality of tubes are in fluid communication with the exhaust fluid flow path. The heat exchanger includes a cap attached to a first end of the shell and defining a second volume. A header is configured to separate the first volume from the second volume, flex with thermal expansion, and define tube inlet and outlet positions. The tube inlets and outlets are in fluid communication with a source fluid flow path, and each tube is substantially U-shaped and defines a flow path of the source fluid within the exhaust fluid flow path. The heat exchanger includes at least one longitudinal flow baffle within the shell configured to reduce an amount of exhaust fluid that may bypass the tubes.
A preliminary application of the underlying principles of the investigator's general system theory to the description and analyses of the fluid flow system is presented. An attempt is made to establish practical models, or elements of the general fluid flow system from the point of view of the general system theory fundamental principles. Results obtained are applied to a simple experimental fluid flow system, as test case, with particular emphasis on the understanding of fluid flow instability, transition and turbulence.
The hybrid optimization software suite (HOSS) is a general purpose fully 2D/3D parallel combined finite discrete element (FDEM) code that can be used to simulate problems involving fracture and fragmentation processes, large deformation and large rotations, discrete particle systems, solid-fluid interaction in hydro fracture problems, etc. HOSS uses a hybrid approach that combines finite-element and discrete-element methods with a novel computational fluid dynamics solver. The finite-element method is often used to analyze a material or object and how it responds to stress. The discrete-element method analyzes stresses and displacements in a volume containing a large number of particles, such as grains of sand. Fluid dynamics analyzes the fluid flow inside, around, or through solid domains. With these processes combined, HOSS represents a paradigm shift when it comes to generating accurate simulations of material deformations and failure. HOSS can perform the following: Resolves problems that consist of millions of deforming fracturing, and interacting solid particles; Requires no coupling—it naturally integrates with all regimes of fluid flow; Resolves all regimes of fluid flow, such as Stokes, low/high Reynolds, subsonic/transonic/supersonic/hypersonic, compressible/incompressible, viscid/inviscid, Newtonian/Non-Newtonian and turbulent/laminar flow; Combines all flow regimes in the same solver; Numerically stable for all flow regimes; Naturally combines different flow regimes in the same problem; Possesses material library for fluid behavior; Naturally matches fluid solver to solid solver; Includes non-inertial Eulerian formulation for fluids; Uses highly efficient parallel computing based on the Virtual Parallel Machine, thus enabling easy porting between different parallel architectures, including possible future architectures; Designed specifically for multi-physics problems in research and industry in virtual experimentation format, thereby complementing both theoretical-computational and experimental tools; Can be extended to address millions of deforming solid particles immersed in fluid flowing at all regimes; Can be extended to address biological applications, such as modeling bone fracture or red blood cells in blood plasma
The flow of fluids in point bars is affected by the existence of heterogeneities like shale drapes that are found on the surfaces of inclined heterolithic stratifications. In fact, these shale drapes can act as fluid flow baffles; therefore, developing a framework for modeling point bars and their associated heterogeneities is vital. In this study, a stochastic process-based modeling approach is presented for capturing the main point bar heterogeneities: accretion surfaces (i.e., the aerial heterogeneity) and inclined heterolithic stratifications (i.e., the vertical heterogeneity). The former was modeled using a sine-generation function and the latter, with a sigmoidal function after which they were combined into a 3D point bar model. To ensure proper modeling of petrophysical properties, we developed a more representative gridding scheme which generates curvilinear grids representative of the point bar geometry. This grid was then transformed into a rectilinear grid to allow for geostatistical simulation after which all petrophysical properties were mapped back into the original curvilinear grid. An essential element of this modeling approach is the stochastic representation of shale drapes at the interface between successive accretion surfaces. The workflow was tested using a real field dataset for the Cranfield field, Mississippi. The constructed model was then subjected to a flow simulation study mimicking a CO 2 storage scenario. Various sensitivities were simulated to evaluate the effect of heterogeneities on CO 2 flow within the point-bar. Results demonstrate the importance of representing point-bar related heterogeneity and the spatial distribution of shale drapes on CO 2 plume migration and storage.
Massive concrete structures such as the containments of nuclear power plant must maintain their tightness at any circumstances to prevent the escape of radioactive fission products into the environment. In the event of an accident like a Loss of Coolant Accident (LOCA), the concrete wall is submitted to both hydric and mechanical loadings. A new experimental device reproducing these extreme conditions (water vapor transfer, 140 °C and 5 bars) is developed in the GeM Laboratory to determine the effect of the saturation degree, the mechanical loading and the flowing fluid type on the concrete transfer properties. The experimental tests show that the previous parameters significantly affect the concrete permeability and the gas leakage rate. Their evolution as a function of the mechanical loading is characterized by two phases that are directly related to concrete microstructure and crack development.
Understanding the flow of fluid, either liquid or gas, through and around solid bodies has challenged man since the dawn of scientific inquiry. Many of the great minds of science and mathematics have progressively built up a hierarchy of fluid models since fluid flow impacts our lives in so many fundamental ways - from the early days of flow of water in viaducts to cities, to today's flight of planes. This chapter is concerned with the computational modeling of turbulent flow around aerodynamic bodies such as planes and wind turbines. In this case, viscous effects near the solid bodies create very thin boundary layers that yield highly anisotropic (gradients normal to the surface may be 10 6 larger than gradients along the surface) solutions to the governing non-linear partial differential equations (PDE); the Navier-Stokes equations. Furthermore, turbulent flows develop extremely broad ranges of length and time scales motivating the use of discretization methods capable of employing adaptivity and implicit time integration. The combination of these features - non-linear, anisotropy, adaptivity, and implicit - dramatically raise the complexity of the discretization, posing large challenges to efficient scalable parallel implementation. However, through careful design, the more complex algorithms can provide great reductions in computational cost relative to simpler methods (e.g., Cartesian grids with explicit time integration) that are easier to mate efficiently to hardware. In this chapter, we not only describe our approach but we also address the fact that while complex algorithms may never be as efficient flop-for-flop as simple methods, in the important measure of science-per-core-hour, they can still win big by making complex features like adaptivity and implicit methods as efficient and scalable as possible.
Technique applies to liquid or gas transfer through flexible ducting and proves valuable in high velocity fluid flow cases. Fluid mechanism responsible for free bellows vibrational excitation also causes flexible hose oscillation. Static pressure stress influences flexible ducting fatigue life and is considered separately.
A self-consistent derivation of the conservation laws is given for flows of a fluid-solid mixture. A unified analytical framework for obtaining constitutive relations is provided. This analysis uses a control volume/control surface approach that is widely used in fluid mechanics. All terms in the governing equations and the constitutive relations are written in terms of the mass-weighted averages except solid concentration. It is believed that the mass-weighted average is the natural bridge between micromechanics and constitutive relations. The derived momentum equations contain terms that differ from all existing models except that of Prosperetti and Jones (1984). However, their assumptions are not needed here. Special attention is given to the solid phase pressure. The physical basis of the previously assumed form for this pressure (Givler 1987) becomes clear. A number of related phenomena are also discussed. These include the anti-diffusion and anisotropic normal stresses. The energy equations are also different from existing models.
Sculptor(R) is a commercially available software tool, based on an Arbitrary Shape Design (ASD), which allows the user to perform shape optimization for computational fluid dynamics (CFD) design. The developed software tool provides important advances in the state-of-the-art of automatic CFD shape deformations and optimization software. CFD is an analysis tool that is used by engineering designers to help gain a greater understanding of the fluid flow phenomena involved in the components being designed. The next step in the engineering design process is to then modify, the design to improve the components' performance. This step has traditionally been performed manually via trial and error. Two major problems that have, in the past, hindered the development of an automated CFD shape optimization are (1) inadequate shape parameterization algorithms, and (2) inadequate algorithms for CFD grid modification. The ASD that has been developed as part of the Sculptor(R) software tool is a major advancement in solving these two issues. First, the ASD allows the CFD designer to freely create his own shape parameters, thereby eliminating the restriction of only being able to use the CAD model parameters. Then, the software performs a smooth volumetric deformation, which eliminates the extremely costly process of having to remesh the grid for every shape change (which is how this process had previously been achieved). Sculptor(R) can be used to optimize shapes for aerodynamic and structural design of spacecraft, aircraft, watercraft, ducts, and other objects that affect and are affected by flows of fluids and heat. Sculptor(R) makes it possible to perform, in real time, a design change that would manually take hours or days if remeshing were needed.
One phenomenon that contributes to nuclear fuel and clad failure mechanisms is the departure from nucleate boiling (DNB). Accident scenarios such as reactivity-initiated accidents (RIAs) provide conditions where the frequency of boiling crisis is much more likely to occur and on much shorter time scales than loss of coolant (LOCA) or loss of flow (LOFA) accidents. Here, boiling crisis in the form of DNB is attributed to transient-induced heat flux values that exceed the critical heat flux (CHF). The resulting effects, such as bubble crowding and eventual dryout, can affect flow boiling heat transfer and compromise the structural integrity of nuclear reactor materials. The focus of the work proposed here lies in developing a rigorous process using out-of-pile testing, under simulated RIA conditions, to determine the transient critical heat flux of cladding materials, without contributions from aforementioned radiation-induced effects that have been observed during in-pile experiments. This study provides a road-map for developing and performing experiments that simulate the power profile and energy deposition aspects of transient events, specific to RIAs, on the laboratory scale. The goal of our project is to observe and measure DNB in an annular flow channel around a single heater rod that is operating and nominal thermal hydraulic conditions of a PWR. It is also an aim of these experiments to determine how to quantify the gap between SS and transient CHF cases, based on state-dependent properties of the cladding material type that is being interrogated. Traditionally, CHF is determined as a combination of four factors: hydraulic diameter, system pressure, mass flow, fluid quality. The authors propose that an additional measureable quantity, which is dependent on the structural properties of the cladding material itself, is partly reasonable for the difference between measured and model-developed CHF values. Using a carefully metered approach will allow for several key concepts to be realized: compensation for transient-based effects to current models for CHF; adaptation of current models for different cladding material types, including Gen-IV concept materials; and implementation of a correction factor to account for both power transient-based and radiation-induced effects on the CHF that have not yet been accounted for.
Ti6Al4V is a widely used, multi-component alloy in additive manufacturing, during which the fluid flow in the molten pool significantly affects the solidified dendrites. To predict and further control the microstructure, modeling and simulating the microstructure evolution play a critical role. In this study, a newly developed, efficient, quantitative multi-component phase-field (PF) model is coupled with a lattice Boltzmann (LB) model to simulate Ti6Al4V solidified dendrite evolution under fluid flow. The accuracy and convergence behavior of the model is validated by the Gibbs–Thomson relation at the dendrite tip. Single and multiple two-dimensional (2D) equiaxed dendrite evolution cases under forced flow were simulated. Results show that the dendrite pattern is influenced remarkably by the fluid flow. Underlying mechanisms of the asymmetrical evolution are revealed by discussing the interaction among the flow, composition distribution and dendrite morphology, quantitatively. The dendrite kinetics are also derived, which ascertains the relationship between tip velocity and undercooling and inlet velocity and is the foundation for larger-scale simulation. We believe that the coupled quantitative multi-component PF–LB framework employed in this study helps in investigating the solidified dendrite morphology evolution in a deep and quantitate manner.
An electrohydrodynamic (EHD) pump increases refrigerant flow rate and the resulting pressure in a vapor compression based cooling system for permitting reduced compressor sizes and power demands. The EHD pump disposes electrodes in a liquid path of the refrigerant flow, and increases fluid flow and resulting pressure by an induced liquid flow between a pair of asymmetric electrodes. Voltage applied to these electrodes results in a conduction pumping mechanism associated with heterocharge layers in the vicinity of the electrodes based on disassociation of a neutral electrolyte species in the refrigerant fluid and recombination of the generated ions. The induced flow draws the liquid due to a net fluid flow toward one of the electrodes based on the asymmetry of the electrode pair. Electrodes are disposed on an inner surface of a refrigerant vessel, in communication with an annular liquid film that forms around the inner circumference in two-phase fluid systems.
Mineral dissolution rates measured in natural environments are much slower than those measured in laboratory settings. This project tested the hypothesis that the way fluid flows through rocks in natural systems creates areas where mineral dissolution is fast and areas where mineral dissolution is slow. This hypothesis was tested with a combination of laboratory experiments and numerical simulation. We demonstrated a separation of fluid flow pathways and rates of mineral dissolution in laboratory experiments for the first time using an experimental approach where we created synthetic rocks that have different ratios of connected and dead-end pathways for fluid flow. In laboratory experiments with higher proportions of dead-end pathways, the mineral dissolution rates were slower. We found that where fluid flows through connected pathways the continuous refreshing of fluid at the mineral surface creates conditions where dissolution is fast. In contrast, where fluid either flows slowly through poorly connected pathways or is stagnant in dead-end pathways, mineral dissolution is slow. The results from this project suggest that the overall slowing of rates of mineral dissolution is important when the proportion of dead-end pathways is greater than ~40%. This project informs our understanding of the way that fluids react with rocks in carbon dioxide sequestration and enhanced geothermal projects where fluids are purposefully injected into rocks for energy applications.
Measurements of drag and wake structure in magneto-fluid dynamic flow about a sphere