Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Fluid interfaces”

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 181 records · Page 10

Simulations of three-layer Richtmyer–Meshkov mixing in a shock tube

Here, the Richtmyer–Meshkov instability causes perturbations to grow after a shock traverses a fluid density interface. This increases the mixing rate between fluid from either side of the interface. We use the Flash Eulerian hydrodynamic code to investigate alterations when a thin third layer of intermediate density is placed along the interface, effectively creating two adjacent unstable interfaces. This is a common occurrence in engineering applications where a thin barrier initially separates two materials. We find that the width of the mixing layer is similar or slightly reduced; however, the total mass of mixed material can actually increase. The mixing layer becomes more compact and efficient. However, the normalized mixed mass decreases, meaning that finger entrainment becomes more important than in the simple two-layer case. The effect of adding the central layer appears to decrease when the Atwood number is decreased. The Flash results are also benchmarked against two-layer experimental data from a shock tube at the University of Arizona.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Orbital Fluid Transfer System

An automated fluid and power interface system needs to be developed for future space missions which require on orbit consumable replenishment. Current method of fluid transfer require manned vehicles and extravehicular activity. Currently the US does not have an automated capability for consumable transfer on-orbit. This technology would benefit both Space Station and long duration satellites. In order to provide this technology the Automated Fluid Interface System (AFIS) was developed. The AFIS project was an advanced development program aimed at developing a prototype satellite servicer for future space operations. This mechanism could transfer propellants, cryogens, fluids, gasses, electrical power, and communications from a tanker unit to the orbiting satellite. The development of this unit was a cooperative effort between Marshall Space Flight Center in Huntsville, Alabama, and Moog, Inc. in East Aurora, New York. An engineering model was built and underwent substantial development testing at Marshall Space Flight Center (MSFC). While the AFIS is not suitable for spaceflight, testing and evaluation of the AFIS provided significant experience which would be beneficial in building a flight unit. The lessons learned from testing the AFIS provided the foundation for the next generation fluid transfer mechanism, the Orbital Fluid Transfer System (OFTS). The OFTS project was a study contract with MSFC and Moog, Inc. The OFTS was designed for the International Space Station (ISS), but its flexible design could used for long duration satellite missions and other applications. The OFTS was designed to be used after docking. The primary function was to transfer bipropellants and high pressure gases. The other items addressed by this task included propellant storage, hardware integration, safety and control system issues. A new concept for high pressure couplings was also developed. The results of the AFIS testing provided an excellent basis for the OFTS design. The OFTS meet the servicing requirements for ISS and could also provide the automated fluid and power interface system needed for on orbit consumable resupply of spacecraft into the new century.

Johnston, A. S., (Nick)↗

Meshfree methods on manifolds for hydrodynamic flows on curved surfaces: A Generalized Moving Least-Squares (GMLS) approach

We utilize generalized moving least squares (GMLS) to develop meshfree techniques for discretizing hydrodynamic flow problems on manifolds. We use exterior calculus to formulate incompressible hydrodynamic equations in the Stokesian regime and handle the divergence-free constraints via a generalized vector potential. This provides less coordinate-centric descriptions and enables the development of efficient numerical methods and splitting schemes for the fourth-order governing equations in terms of a system of second-order elliptic operators. Using a Hodge decomposition, we develop methods for manifolds having spherical topology. We show the methods exhibit high-order convergence rates for solving hydrodynamic flows on curved surfaces. The methods also provide general high-order approximations for the metric, curvature, and other geometric quantities of the manifold and associated exterior calculus operators. The approaches also can be utilized to develop high-order solvers for other scalar-valued and vector-valued problems on manifolds.

97 MATHEMATICS AND COMPUTING↗

Enriched immersed boundary method (EIBM) for interface-coupled multi-physics and applications to convective conjugate heat transfer

It is challenging to develop numerical methods that simultaneously maintain accuracy of resolving boundary conditions and mesh flexibility to handle the interface in interface-coupled multi-physics problems involving large property discontinuity. On the one hand, boundary-fitted methods possess high accuracy in capturing the interfacial phenomenon but involve complicated volumetric mesh generation, mesh-motion, and even re-meshing procedures. On the other hand, immersed boundary methods (IBM) provide mesh flexibility but sometimes suffer from inferior interface representations, leading to poor enforcement of boundary conditions. This paper presents an enriched immersed boundary method (EIBM) to overcome this challenge and demonstrates its efficacy in conjugate heat transfer, a representative example in interface-coupled multi-physics systems that have implications for many industrial processes. Further, the core technique of the method is to enhance IBM’s accuracy of the fluid–solid interface by enriching the degrees of freedom of the cut elements to enforce temperature and flux compatibilities and resolve all the physical unknowns on the background mesh to simplify the volumetric mesh generation. We implement the EIBM under the framework of a variational multiscale formulation for coupled Navier–Stokes and thermodynamics equations. The enriched DoFs enable better enforcement of temperature and flux compatibilities with large conductivity ratios across the fluid–solid interface. At the same time, the immersed nature of the proposed method still attains mesh flexibility. The EIBM’s accuracy is thoroughly evaluated through a set of examples, ranging from benchmark problems with analytical solutions to real-world cooling processes of a moving metallic structure with complex geometry.

42 ENGINEERING↗

Two-dimensional glass transition–like behavior of Janus particle–laden interface

Abstract Understanding the interactive behavior of Janus particles (JPs) is a growing field of research. The enhancement in binding energy, in comparison to homogenous particles, and the dual characteristic of JPs open up new possibilities for novel applications. In many such applications, interfacial materials become subjected to flows that produce dilational and shear stresses. Therefore, it is important to understand the impact that the Janus character brings to interfaces. In this work, we study the microstructure of two-dimensional (2D) JP monolayers formed at the air–water interface and examine the shear viscoelasticity with an interface rheometer that was adapted for in situ surface pressure control via a Langmuir trough. We extend concepts from bulk rheology to data obtained from interfacial rheology as a tool to understand and predict the monolayer’s viscoelastic behavior. Finally, by calculating the time relaxation spectrum from the measured 2D dynamic moduli, we conclude that a phenomenon similar to glass transition is taking place by analogy.

42 ENGINEERING↗

Acoustic responses of coupled fluid-structure system by acoustic-structural analogy

The use of an analogy between structural mechanics and acoustics makes it possible to solve fluid-structural interaction (FSI) problems using an existing structural analysis computer program. This method was implemented in MSC/NASTRAN program and the FSI analysis was performed using two dimensional coupled fluid beam model to assess and evaluate the adequacy of this approach. The coupled modal analysis of 3-D model is also briefly discussed. The normal mode, modal frequency response and transient response analysis of 2-D coupled fluid beam system is presented. The significant reduction of the acoustic pressure response at the fluid structure interface is observed as a result of fluid structure interaction.

Shin, Y. S.↗

Shear Stress Sensing using Elastomer Micropillar Arrays

The measurement of shear stress developed as a fluid moves around a solid body is difficult to measure. Stresses at the fluid-solid interface are very small and the nature of the fluid flow is easily disturbed by introducing sensor components to the interface. To address these challenges, an array of direct and indirect techniques have been investigated with various advantages and challenges. Hot wire sensors and other indirect sensors all protrude significantly into the fluid flow. Microelectromechanical systems (MEMS) devices, although facilitating very accurate measurements, are not durable, are prone to contamination, and are difficult to implement into existing model geometries. One promising approach is the use of engineered surfaces that interact with fluid flow in a detectable manner. To this end, standard lithographic techniques have been utilized to generate elastomeric micropillar arrays of various lengths and diameters. Micropillars of controlled length and width were generated in polydimethylsiloxane (PDMS) elastomer using a soft-lithography technique. The 3D mold for micropillar replication was fabricated using laser ablative micromachining and contact lithography. Micropillar dimensions and mechanical properties were characterized and compared to shear sensing requirements. The results of this characterization as well as shear stress detection techniques will be discussed.

Wohl, Christopher J.↗

Assessment of diffuse-interface methods for compressible multiphase fluid flows and elastic-plastic deformation in solids

This work describes three diffuse-interface methods for the simulation of immiscible, compressible multiphase fluid flows and elastic-plastic deformation in solids. The first method is the localized-artificial-diffusivity approach of Cook, Subramaniam et al., and Adler and Lele, in which artificial diffusion terms are added to the individual phase mass fraction transport equations and are coupled with the other conservation equations. The second method is the gradient-form approach that is based on the quasi-conservative method of Shukla et al., in which the diffusion and sharpening terms (together called regularization terms) are added to the individual phase volume fraction transport equations and are coupled with the other conservation equations. The third approach is the divergence-form approach that is based on the fully conservative method of Jain et al., in which the regularization terms are added to the individual phase volume fraction transport equations and are coupled with the other conservation equations. In the present study, all three diffuse-interface methods are used in conjunction with a four-equation, multicomponent mixture model, in which pressure and temperature equilibria are assumed among the various phases. The primary objective of this work is to compare these three methods in terms of their ability to: maintain constant interface thickness throughout the simulation; conserve mass, momentum, and energy; and maintain accurate interface shape for long-time integration. The second objective of this work is to consistently extend these methods to model interfaces between solid materials with strength. To assess and compare the methods, they are used to simulate a wide variety of problems, including (1) advection of an air bubble in water, (2) shock interaction with a helium bubble in air, (3) shock interaction and the collapse of an air bubble in water, and (4) Richtmyer–Meshkov instability of a copper–aluminum interface. The current work focuses on comparing these methods in the limit of relatively coarse grid resolution, which illustrates the true performance of these methods. In conclusion, this is because it is rarely practical to use hundreds of grid points to resolve a single bubble or drop in large-scale simulations of engineering interest.

97 MATHEMATICS AND COMPUTING↗

Assessment of the Fluid Dynamics Boundary Condition in Ablating or Blowing Flows

Improved models of ablative thermal protection systems have enabled the treatment of materials and fluid behavior in a coupled manner. This paper reports a new approach to modeling the interface between fluid and material, with attention to the conservation of species mass flux and energy on the fluid side of the interface. The general equation is presented and is shown to recover the traditional uncoupled fluid/materials response interface. Including the chemical reaction terms on the CFD side of the interface makes the heat flux exchange independent of the thermodynamic reference state and, therefore, a measurable quantity. Doing so allows the material response solver to take as input the surface heat flux rather than a film coefficient. Removing the film coefficient approximation enables more direct solution of vehicle thermal response but requires consistency in the wall state. The mixing of the shock layer and pyrolysis gas is then computed with finite rate chemistry within the fluid solver. The boundary conditions described have been implemented in the DPLR v4.05.1 code. Char removal is captured using finite rate chemistry in DPLR’s gas surface interaction module. Aspects of coupling these solutions to material response are discussed.

Ablation↗

Assessment of the Fluid Dynamics Boundary Condition in Ablating or Blowing Flows

Improved models of ablative thermal protection systems have enabled the treatment of materials and fluid behavior in a coupled manner. This paper reports a new approach to modeling the interface between fluid and material, with attention to the conservation of species mass flux and energy on the fluid side of the interface. The general equation is presented and is shown to recover the traditional uncoupled fluid/materials response interface. Including the chemical reaction terms on the CFD side of the interface makes the heat flux exchange independent of the thermodynamic reference state and, therefore, a measurable quantity. Doing so allows the material response solver to take as input the surface heat flux rather than a film coefficient. Removing the film coefficient approximation enables more direct solution of vehicle thermal response but requires consistency in the wall state. The mixing of the shock layer and pyrolysis gas is then computed with finite rate chemistry within the fluid solver. The boundary conditions described have been implemented in the DPLR v4.05.1 code. Char removal is captured using finite rate chemistry in DPLR’s gas surface interaction module. Aspects of coupling these solutions to material response are discussed.

Ablation↗

Role of pH and Eh in geothermal systems: Thermodynamic examples and impacts on scaling and corrosion

The global demand for sustainable and clean renewable energy sources is rapidly increasing to meet growing energy demands while curtailing the CO2 emissions from fossil fuels. Harnessing the natural heat or geothermal energy of the Earth's interior can simultaneously address both issues by offering a clean and renewable energy source. Enhanced geothermal systems require the movement and cycling of fluid through a porous rock media. Further, this process requires understanding the chemistry occurring at the fluid-rock and fluid-metal interfaces both in the deep subsurface and surface, where fluid is circulated for energy transfer or storage. In this review paper, we evaluate the current state of research on how two key variables, pH and Eh, control the fluid-rock reactions controlling mineral scaling and the steel corrosion reactions associated with geothermal systems. These variables are frequently used in geochemical modeling to predict chemical reactions in geothermal systems. Prediction of the fluid-rock behavior can be utilized to formulate new strategies for extending the lifespan and efficiency of heat extraction and utilization and maximizing the safety of geothermal systems.

15 GEOTHERMAL ENERGY↗

Improving the reliability of machine learned potentials for modeling inhomogeneous liquids

The atomic-scale response of inhomogeneous fluids at interfaces and surrounding solute particles plays a critical role in governing chemical, electrochemical, and biological processes. Classical molecular dynamics simulations have been applied extensively to simulate the response of fluids to inhomogeneities directly, but are limited by the accuracy of the underlying interatomic potentials. Here, we use neural network potentials (NNPs) trained to ab initio simulations to accurately predict the inhomogeneous responses of two distinct fluids: liquid water and molten NaCl. Although NNPs can be readily trained to model complex bulk systems across a range of state points, we show that to appropriately model a fluid's response at an interface, relevant inhomogeneous configurations must be included in the training data. In order to sufficiently sample appropriate configurations of such inhomogeneous fluids, we develop protocols based on molecular dynamics simulations in the presence of external potentials. We demonstrate that NNPs trained on inhomogeneous fluid configurations can more accurately predict several key properties of fluids—including the density response, surface tension and size-dependent cavitation free energies—for liquid water and molten NaCl, compared to both empirical interatomic potentials and NNPs that are not trained on such inhomogeneous configurations. This work therefore provides a first demonstration and framework to extract the response of inhomogeneous fluids from first principles for classical density-functional treatment of fluids free from empirical potentials.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Surface-Tension Driven Convection

The present experiments are directed to examine surface-tension gradients. The surface-tension induced convection under reduced gravity it's comparison with normal gravity convection with all other identical conditions are studied. The intensities of such flows and their penetration depth below the free surface are also determined. In this way, the dominant role of a free surface in a reduced-gravity environment is indicated. The interface between two fluid phases can influence the motion of fluids when either the interface has finite curvature or when the interfacial tension varies from point to point. In both cases, forces appear in the interfacial region that can affect or generate fluid motions.

Ostrach, S.↗

Aerothermal Shape Optimization of Actively-Cooled Battery Packs using Conjugate Heat Transfer

Thermal management for battery is important for electric aircraft because battery temperature is critically important to vehicle safety, and it also has direct impact on the efficiency of the battery system. Because ambient air is a readily available resource for aircraft, this paper considers an active cooling concept with forced convection of ambient air through the battery pack. Conjugate heat transfer analysis is used to solve the coupled aero-thermal problem, which consists of a finite-volume computational fluid dynamics solver for the fluid domain, and a conduction heat transfer solver for the solid domain. A mixed Neumann and Dirichlet boundary condition is developed for the fluid-solid interface, which allows the solid domain to completely submerge in the fluid domain. A gradient-based optimization method is adopted, and the discrete adjoint approach implemented in DAFoam is used to efficiently compute the gradients. The aero-thermal coupling for primal analysis and gradient computation is handled using the OpenMDAO-based MPhys framework. A constant heat source is prescribed for the battery cells, and the battery shape (design variable) is optimized to minimize cooling pump power and battery weight (composite objective function) while keeping the battery temperature below a threshold (constraint). The optimized design achieves a 44.6% and 1.5% reduction in the cooling pump power and battery weight, respectively, and the maximal temperature constraint is satisfied. This work has the potential to reduce battery-pack weight, improve performance, and reduce the weight of thermal management systems for electric vertical take-off and landing aircraft.

thermal management↗

Aerothermal Shape Optimization of Actively-Cooled Battery Packs Using Conjugate Heat Transfer

Thermal management for battery is important for electric aircraft because battery temperature is critically important to vehicle safety, and it also has direct impact on the efficiency of the battery system. Because ambient air is a readily available resource for aircraft, this paper considers an active cooling concept with forced convection of ambient air through the battery pack. Conjugate heat transfer analysis is used to solve the coupled aero-thermal problem, which consists of a finite-volume computational fluid dynamics solver for the fluid domain, and a conduction heat transfer solver for the solid domain. A mixed Neumann and Dirichlet boundary condition is developed for the fluid-solid interface, which allows the solid domain to completely submerge in the fluid domain. A gradient-based optimization method is adopted, and the discrete adjoint approach implemented in DAFoam is used to efficiently compute the gradients. The aero-thermal coupling for primal analysis and gradient computation is handled using the OpenMDAO-based MPhys framework. A constant heat source is prescribed for the battery cells, and the battery shape (design variable) is optimized to minimize cooling pump power and battery weight (composite objective function) while keeping the battery temperature below a threshold (constraint). The optimized design achieves a 44.6% and 1.5% reduction in the cooling pump power and battery weight, respectively, and the maximal temperature constraint is satisfied. This work has the potential to reduce battery-pack weight, improve performance, and reduce the weight of thermal management systems for electric vertical take-off and landing aircraft.

heat transfer↗

High accuracy computation of fluid-structure interaction in transonic cascades

A coupling strategy for simulating fluid-structure interaction phenomena is formulated and applied to the prediction of flutter in transonic cascades. The flow is governed by the Euler equations and discretized using a finite volume flux-splitting scheme. The structure is modeled using an isoparametric finite element formulation. The coupling strategy successfully reconciles these two formulations at the fluid-structure interface by enforcing both kinematic and kinetic boundary conditions. In particular, the conservation laws applicable to the combined fluid-structure system are preserved across the interface. Since the primary mechanism driving aeroelastic phenomena involves energy exchange occurring at the interface, this highly accurate coupling mechanism is believed to improve the predictive capability of the scheme. The coupled equations are advanced simultaneously in time using an implicit time integration method. Results obtained using the coupling method are presented for cascade geometries operating in transonic flow.

Boschitsch, Alexander↗

Investigation of Conjugate Heat Transfer in Turbine Blades and Vanes

We report on work carried out to develop a 3-D coupled Finite Volume/BEM-based temperature forward/flux back (TFFB) coupling algorithm to solve the conjugate heat transfer (CHT) which arises naturally in analysis of systems exposed to a convective environment. Here, heat conduction within a structure is coupled to heat transfer to the external fluid which is convecting heat into or out of the solid structure. There are two basic approaches to solving coupled fluid structural systems. The first is a direct coupling where the solution of the different fields is solved simultaneously in one large set of equations. The second approach is a loose coupling strategy where each set of field equations is solved to provide boundary conditions for the other. The equations are solved in turn until an iterated convergence criterion is met at the fluid-solid interface. The loose coupling strategy is particularly attractive when coupling auxiliary field equations to computational fluid dynamics codes. We adopt the latter method in which the BEM is used to solve heat conduction inside a structure which is exposed to a convective field which in turn is resolved by solving the NASA Glenn compressible Navier-Stokes finite volume code Glenn-HT. The BEM code features constant and bi-linear discontinuous elements and an ILU-preconditioned GMRES iterative solver for the resulting non-symmetric algebraic set arising in the conduction solution. Interface of flux and temperature is enforced at the solid/fluid interface, and a radial-basis function scheme is used to interpolated information between the CFD and BEM surface grids. Additionally, relaxation is implemented in passing the fluxes from the conduction solution to the fluid solution. Results from a simple test example are reported.

Kassab, A. J.↗