Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “FLUID FLOW”

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 883 records · Page 49

Machine learning to identify geologic factors associated with production in geothermal fields: A casestudy using 3D geologic data, Brady geothermal field, Nevada

In this paper, we present an analysis using unsupervised machine learning (ML) to identify the key geologic factors that contribute to the geothermal production in the Brady geothermal field. Brady is a hydrothermal system in northwestern Nevada that supports both electricity production and direct use of hydrothermal fluids. Transmissive fluid flow pathways are relatively rare in the subsurface but are critical components of hydrothermal systems like Brady and many other types of fluid flow systems in fractured rock. The ML method, non-negative matrix factorization with k-means clustering (NMFk), is applied to a library of fourteen 3D geologic characteristics hypothesized to control hydrothermal circulation in the Brady geothermal field. Our results indicate the macro-scale faults and a local step-over in the fault system preferentially occur along with production wells when compared to injection wells and non-productive wells. We infer that these are the key geologic characteristics that control the through-going hydrothermal transmission pathways at Brady. Our results demonstrate 1) the specific geologic controls on the Brady hydrothermal system and 2) the efficacy of pairing ML techniques with 3D geologic characterization to enhance the understanding of subsurface processes.

15 GEOTHERMAL ENERGY↗

The Role of Marangoni Convection for the FZ-Growth of Silicon

Fluctuations of the electrical resistivity due to inhomogeneous dopant distribution are still a serious problem for the industrial processing yield of doped silicon crystals. In the case of silicon floating-zone growth, the main sources of these inhomogeneities are time- dependent flows in the liquid phase during the growth process. Excluding radio frequency (RF) induced convection, buoyancy and thermocapillary (Marangoni) convection are the two natural reasons for fluid flow. Both originate from temperature/concentration gradients in the melt, buoyancy convection through thermal/concentrational volume expansion, and thermocapillary convection through the temperature/concentration dependence of the surface tension. To improve the properties of grown crystals, knowledge of the strength, the characteristic, and the relation of these two flow mechanisms is essential. By the use of microgravity, the effect and the strength of buoyancy (gravity dependent) and thermocapillary (gravity independent) convection can be separated and clarified. Applying magnetic fields, both convective modes can be influenced: fluid flow can either be damped (static magnetic fields) or overlaid by a regular flow regime (rotating magnetic fields). Two complementary approaches have been pursued: Silicon full zones (experiments on the German sounding rockets TEXUS 7, 12, 22, 29, and 36) with the maximum temperature at half of the zone height and silicon half zones (experiments on the Japanese sounding rockets TR-IA4 and 6) with the maximum temperature at the top of the melt. With the full zone arrangement, the intensity and the frequency of the dopant striations could be determined and the critical Marangoni number could be identified. The half zone configuration is suited to classify the flow pattern and to measure the amplitude and the frequency of temperature fluctuations in the melt by inserting thermocouples or temperature sensors into the melt. All experiments have been carried out in monoellipsoid mirror furnaces. Typical zone geometries are approx. 8 to 14 mm in diameter and height. The crystals grown under microgravity are compared to crystals grown in static axial magnetic fields (B<5 tesla) and in transversal rotating magnetic fields (B<7.5 mT / f=50 Hz). Experimental results are completed by 3D numerical simulations: the obtained temperature and concentration distribution in the melt confirm the damping effect of rotating magnetic fields and explain the change in the radial segregation under static magnetic fields.

Dold, P.↗

A Multiscale Approach to Simulate Non‐Isothermal Multiphase Flow in Deformable Porous Materials

Coupled thermal, hydraulic, and mechanical processes in porous materials play important roles in several energy and environmental technologies. The Darcy-Brinkman-Biot (DBB) framework has proven effective in modeling multiphase fluid flow in deformable porous solids across both pore and Darcy scales, including in systems where fractures coexist with a porous matrix. In this study, we extend the DBB framework, originally designed for isothermal conditions, to address non-isothermal problems by incorporating an energy conservation equation. The resulting solver, hybridBiotThermalInterFoam, enables simulations of coupled multiphase fluid flow, heat transfer, and solid deformation in hybrid-scale systems containing both solid-free regions and ductile porous domains. The new solver is validated through comparisons with analytical solutions and, also, against established heat transfer solvers chtMultiRegionFoam and compressibleInterFoam. Further, a series of 2D and 3D case studies, including two-phase heat transfer in solid-free, static, or deformable porous media, highlights the solver's capacity to simulate complex flow dynamics and heat transport in systems involving high mobility ratios, viscous fingering, and fracture propagation. Our results establish the feasibility of incorporating thermal effects in simulations of a wide variety of energy geotechnics and environmental applications, including enhanced hydrocarbon recovery, soil remediation, and enhanced geothermal energy systems.

04 OIL SHALES AND TAR SANDS↗

Machine learning to identify geologic factors associated with production in geothermal fields: a case-study using 3D geologic data, Brady geothermal field, Nevada

Abstract In this paper, we present an analysis using unsupervised machine learning (ML) to identify the key geologic factors that contribute to the geothermal production in Brady geothermal field. Brady is a hydrothermal system in northwestern Nevada that supports both electricity production and direct use of hydrothermal fluids. Transmissive fluid-flow pathways are relatively rare in the subsurface, but are critical components of hydrothermal systems like Brady and many other types of fluid-flow systems in fractured rock. Here, we analyze geologic data with ML methods to unravel the local geologic controls on these pathways. The ML method, non-negative matrix factorization with k -means clustering (NMF k ), is applied to a library of 14 3D geologic characteristics hypothesized to control hydrothermal circulation in the Brady geothermal field. Our results indicate that macro-scale faults and a local step-over in the fault system preferentially occur along production wells when compared to injection wells and non-productive wells. We infer that these are the key geologic characteristics that control the through-going hydrothermal transmission pathways at Brady. Our results demonstrate: (1) the specific geologic controls on the Brady hydrothermal system and (2) the efficacy of pairing ML techniques with 3D geologic characterization to enhance the understanding of subsurface processes.

58 GEOSCIENCES↗

Low frequency vibrational stirring - A new method for rapidly mixing solutions and melts during growth

A new method for rapidly stirring liquids has been developed which may be useful in crystal growth applications. The method involves the introduction of orthogonally coupled low frequency vibrations to fluids contained in cylindrical crucibles. A mathematical model has been developed to describe the influence of this vibrational stirring method on fluid flow. The optimum sloshing frequencies calculated for a specific system geometry were very close to those observed experimentally. A systematic study of the fluid flow induced by the vibrational motion was carried out at room temperature using aqueous solutions containing glycerine. Preliminary experiments were carried out to study the effect of vibrational stirring on interface shape and position during the Bridgman growth of CsCdCl3, a model system.

Liu, W.-S.↗

Liquid Oxygen/Liquid Methane Integrated Propulsion System Test Bed

In support of NASA?s Propulsion and Cryogenic Advanced Development (PCAD) project, a liquid oxygen (LO2)/liquid methane (LCH4) Integrated Propulsion System Test Bed (IPSTB) was designed and advanced to the Critical Design Review (CDR) stage at the Johnson Space Center. The IPSTB?s primary objectives are to study LO2/LCH4 propulsion system steady state and transient performance, operational characteristics and to validate fluid and thermal models of a LO2/LCH4 propulsion system for use in future flight design work. Two phase thermal and dynamic fluid flow models of the IPSTB were built to predict the system performance characteristics under a variety of operating modes and to aid in the overall system design work. While at ambient temperature and simulated altitude conditions at the White Sands Test Facility, the IPSTB and its approximately 600 channels of system instrumentation would be operated to perform a variety of integrated main engine and reaction control engine hot fire tests. The pressure, temperature, and flow rate data collected during this testing would then be used to validate the analytical models of the IPSTB?s thermal and dynamic fluid flow performance. An overview of the IPSTB design and analytical model development will be presented.

Flynn, Howard↗

Navier-Stokes calculations for unsteady three-dimensional vortical flows in unbounded domains

Finite-difference Navier-Stokes calculations for unsteady, three-dimensional, incompressible, viscous flows induced by initial vorticity distributions are presented and discussed in this paper. The initial vorticity distributions are assumed to be embedded in a flow field of infinite extent that is quiescent at infinity. These vorticity distributions are typical of vortex rings and other closed vortical tubes or structures. Such structures are important elements in fluid flows such as jets, atmospheric convection and the far-field wakes of aircraft; studies of their interaction may aid in an understanding of complex fluid flows. The calculations employ a method recently proposed by Ting to approximate the infinite-domain boundary value problem with a finite boundary computational domain, and this method is shown to yield accurate three-dimensional results for reasonable expenditures of computer time. Because of the efficiency of the boundary condition technique and the resulting Navier-Stokes code, a 16-bit minicomputer with virtual memory was capable of performing the calculations for the unsteady motion of two obliquely colliding vortex rings. The results of these calculations are presented in the paper.

Chamberlain, J. P.↗

Modeling Three-Dimensional Fluid-Driven Propagation of Multiple Fractures using TOUGH-FEMM

Here, a novel numerical simulation tool named TOUGH-FEMM is presented and applied to model three-dimensional (3D) hydraulic fracturing in porous rock. The fluid flow in both fractures and porous rock is modeled using TOUGH2, which is a well-established code for analysis of multiphase and multi-component fluid flow. Rock deformations associated with fracture propagation are modeled using finite element-meshfree method (FEMM). FEMM is an approach to simulate fracture propagation without remeshing, in which the fracture path does not need to be predetermined. Fracture mechanics with mixed-mode stress intensity factors are employed to detect fracture instability and determine the direction of fracture propagation. TOUGH-FEMM is confirmed for modeling fluid-driven fracture propagation in 3D through a number of simulation examples, including modeling of hydraulic fracturing laboratory experiments and by comparison to independent numerical simulation results for multiple interacting hydraulic fractures at ten to hundred meter scale.

58 GEOSCIENCES↗

An analytical thermohydraulic model for discretely fractured geothermal reservoirs

In discretely fractured reservoirs such as those found in Enhanced/Engineered Geothermal Systems (EGS), knowledge of the fracture network is important in understanding the thermal hydraulics, i.e., how the fluid flows and the resulting temporal evolution of the subsurface temperature. Here, the purpose of this study was to develop an analytical model of the fluid flow and heat transport in a discretely fractured network that can be used for a wide range of modeling applications and serve as an alternative analysis tool to more computationally intensive numerical codes. Given the connectivity and structure of a fracture network, the flow in the system was solved using a linear system of algebraic equations for the pressure at the nodes of the network. With the flow determined, the temperature in the fracture was solved by coupling convective heat transport in the fracture with one-dimensional heat conduction perpendicular to the fracture, employing the Green’s function derived solution for a single discrete fracture. The predicted temperatures along the fracture surfaces from the analytical solution were compared to numerical simulations using the TOUGH2 reservoir code. Through two case studies, we showed the capabilities of the analytical model and explored the effect of uncertainty in the fracture apertures and network structure on thermal performance. While both sources of uncertainty independently produce large variations in production temperature, uncertainty in the network structure, whenever present, had a predominant influence on thermal performance.

58 GEOSCIENCES↗

Harnessing surface-bound enzymatic reactions to organize microcapsules in solution

By developing new computational models, we examine how enzymatic reactions on an underlying surface can be harnessed to direct the motion and organization of reagent-laden microcapsules in a fluid-filled microchannel. In the presence of appropriate reagents, surface-bound enzymes can act as pumps, which drive large-scale fluid flows. When the reagents diffuse through the capsules’ porous shells, they can react with enzymatic sites on the bottom surface. The ensuing reaction generates fluid density variations, which result in fluid flows. These flows carry the suspended microcapsules and drive them to aggregate into “colonies” on and near the enzyme-covered sites. This aggregation continues until the reagent has been depleted and the convection stops. Here, we show that the shape of the assembled colonies can be tailored by patterning the distribution of enzymes on the surface. This fundamental physicochemical mechanism could have played a role in the self-organization of early biological cells (protocells) and can be used to regulate the autonomous motion and targeted delivery of microcarriers in microfluidic devices.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Multiscale Simulation of Moist Global Atmospheric Flows

The overarching goal of this award was to include phase changes of the water substance and accompanying latent heating and precipitation processes into the all-scale nonhydrostatic atmospheric dynamics EUlerian/LAGrangian (EULAG) model. The model includes fluid flow solver that is based on either an unabbreviated set of the governing equations (i.e., compressible dynamics) or a simplified set of equations without sound waves (i.e., sound-proof, either anelastic or pseudo-incompressible). The latter set has been used in small-scale dynamics for decades, but its application to the all-scale dynamics (from small-scale to planetary) has never been studied in practical implementations. The highlight of the project is the development of the moist implicit compressible model that can be run by applying time steps, as long as the anelastic model is limited only by the computational stability of the fluid flow and not by the speed of sound waves that limit the stability of explicit compressible models. Applying various versions of the EULAG model within the same numerical framework allows for an unprecedented comparison of solutions obtained with various sets of the governing equations and straightforward evaluation of the impact of various physical parameterizations on the model solutions. The main outcomes of this study are reported in three papers, two published and one currently under review. These papers include comparisons between model solutions for idealized moist problems across the range of scales from small to planetary. These tests include: moist thermals rising in the stable-stratified environment (following Grabowski and Clark, J. Atmos. Sci. 1991) and in the moist-neutral environment (after Bryan and Fritsch, Mon. Wea. Rev. 2002), moist flows over a mesoscale topography (as in Grabowski and Smolarkiewicz, Mon. Wea. Rev. 2002), deep convection in a sheared environment (following Weisman and Klemp, Mon. Wea. Rev. 1982), moist extension of the baroclinic wave on the sphere of Jablonowski and Williamson (Q. J. R. Met. Soc. 2006), and moist extension of the Held-Suarez idealized climate benchmark (Held and Suarez, Bull. Amer. Met. Soc., 1994).

54 ENVIRONMENTAL SCIENCES↗

Fiber optic liquid mass flow sensor and method

A method and apparatus are provided for sensing the mass flow rate of a fluid flowing through a pipe. A light beam containing plural individual wavelengths is projected from one side of the pipe across the width of the pipe so as to pass through the fluid under test. Fiber optic couplers located at least two positions on the opposite side of the pipe are used to detect the light beam. A determination is then made of the relative strengths of the light beam for each wavelength at the at least two positions and based at least in part on these relative strengths, the mass flow rate of the fluid is determined.

Korman, Valentin↗

The influence of acceleration forces on dendritic growth and grain structure

The results of experiments on the tin-15 wt pct lead system are presented, showing the effects on microstructure of solidification in the presence of acceleration forces from 0.0001 to 5 g for three cooling rates. An increase in the acceleration level is shown to drive fluid flow and cause dendrite remelting, fragmentation, and macrosegregation. The cooling rate impacts the final structure through its control of dendrite arm spacings and permeability to fluid flow. At the low (0.0001 g) acceleration, dendrite arm spacings deviated from the predicted relationship to cooling rate. An explanation for this anomaly is given which considers the temperature and concentration gradients in the low-gravity environment.

Johnston, M. H.↗

System for creating at a site, remote from a sterile environment, a parenteral solution

The present invention relates to a container, system, and method for creating parenteral solutions at a site, remote from sterile environments. The system includes a flexible container that is empty except for a prepackaged amount of a solute that is housed in the interior of the container. The container includes at least one port and a sterilizing filter in communication with an interior of the port. The container is so constructed and arranged that a fluid flow path is created from the port through the filter and into the interior of the container. A sterile water source including means for establishing fluid flow from the sterile water source into the port is provided. Accordingly, sterile water can flow from the sterile water source through the filter into the container where it is mixed with the solute to create a parenteral solution that can then be infused into a patient. A method and container are also provided.

Scharf, Mike↗

Efficient simulation of press hardening process through integrated structural and CFD analyses

Press hardened steel parts are being increasingly used in automotive structures for their higher strength to meet safety standards while reducing vehicle weight to improve fuel consumption. However, manufacturing of sheet metal parts by press hardening process to achieve desired properties is extremely challenging as it involves complex interaction of plastic deformation, metallurgical change, thermal distribution, and fluid flow. Numerical simulation is critical for successful design of the process and to understand the interaction among the numerous process parameters to control the press hardening process in order to consistently achieve desired part properties. Until now there has been no integrated commercial software solution that can efficiently model the complete process from forming of the blank, heat transfer between the blank and tool, microstructure evolution in the blank, heat loss from tool to the fluid that flows through water channels in the tools. In this study, a numerical solution based on Altair HyperWorks® product suite involving RADIOSS®, a non-linear finite element based structural analysis solver and AcuSolve®, an incompressible fluid flow solver based on Galerkin Least Square Finite Element Method have been utilized to develop an efficient solution for complete press hardening process design and analysis. RADIOSS is used to handle the plastic deformation, heat transfer between the blank and tool, and microstructure evolution in the blank during cooling. While AcuSolve is used to efficiently model heat loss from tool to the fluid that flows through water channels in the tools. The approach is demonstrated through some case studies.

36 MATERIALS SCIENCE↗

Spatial characterization of heterogeneous nanopore surfaces from XCT scans of Niobrara shale

The multiphase fluid flow through natural and geo-architected nanopore structures is fundamentally controlled by the geometry of such nanopore structures, pore wettability and interfacial tension. Furthermore, multiphase fluid flow through the nanopore structures is compromised when the pore surfaces are heterogeneous, having surface regions of different composition, polarity, and wetting characteristics. In this regard, compositional analysis of the pore channel surface has been a big challenge in the past, since 3D analysis is required. Currently, due to advances in nano X-ray CT, with a voxel resolution of 20 nm, it is possible to identify nanopore network structures and quantify their surface composition in 3D. In this paper, the development of procedures for detailed analysis to quantify the composition of the pore channel surface regions, including the corresponding area, is described using nanoCT images. Image analysis procedures, including preprocessing of the 3D tomographic images, correction for the partial volume effect (PVE), and the use of a new algorithm for determination of pore channel surface composition and area are reported on. Videos of the nanopore structure have been provided as supplementary information to further describe the heterogeneity of the nanopore surfaces.

36 MATERIALS SCIENCE↗

Hot foil transducer skin friction sensor

The device utilizes foil transducers with only one edge exposed to the fluid flow. The surfaces are polished producing a foil transducer that does not generate turbulence while sufficiently thick to carry the required electrical current for high temperature fluid flow. The assembly utilizes a precut layered metal sandwich with attached electrodes eliminating a need for welding and individual sensor calibration.

Vranas, T.↗

Direct Numerical Simulation of the Flow through a Randomly Packed Pebble Bed

The proposition for molten salt and high-temperature gas-cooled reactors has increased the focus on the dynamics and physics in randomly packed pebble beds. Research is being conducted on the validity of these designs as a possible contestant for the fourth-generation nuclear power systems. A detailed understanding of the coolant flow behavior is required in order to ensure proper cooling of the reactor core during normal and accident conditions. In order to increase the understanding of the flow through these complex geometries and enhance the accuracy of lower-fidelity modeling, high-fidelity approaches such as direct numerical simulation (DNS) can be utilized. Nek5000, a spectral-element computational fluid dynamics (CFD) code, was used to develop DNS fluid flow data. The flow domain consisted of 147 pebbles enclosed by a bounding wall. In the work presented, the Reynolds numbers ranged from 430 to 1050 based on the pebble diameter and inlet velocity. Characteristics of the flow domain such as volume averaged porosity, axial porosity, and radial porosity were studied and compared with correlations available in the literature. Friction factors from the DNS results for all Reynolds numbers were compared with correlations in the literature. The first- and second-order statistics show good agreement with the available experimental data. Turbulence length scales were analyzed in the flow. Reynolds stress anisotropy was characterized by utilizing invariant analysis. Overall, the results of the analysis in this study provide deeper understanding of the flow behavior and the effect of the wall in packed beds.

Yildiz, Mustafa Alper↗