Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Fluids”

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 1,369 records · Page 76

Coupled modeling of a directly heated tubular solar receiver for supercritical carbon dioxide Brayton cycle: Optical and thermal-fluid evaluation

In single phase performance and appealing thermo-physical properties supercritical carbon dioxide (s-CO 2 ) make a good heat transfer fluid candidate for concentrating solar power (CSP) technologies. The development of a solar receiver capable of delivering s-CO 2 at outlet temperatures ~973 K is required in order to merge CSP and s-CO 2 Brayton cycle technologies. A coupled optical and thermal-fluid modeling effort for a tubular receiver is undertaken to evaluate the direct tubular s-CO 2 receiver’s thermal performance when exposed to a concentrated solar power input of ~0.3–0.5 MW. Ray tracing, using SolTrace, is performed to determine the heat flux profiles on the receiver and computational fluid dynamics (CFD) determines the thermal performance of the receiver under the specified heating conditions. Moreover, an in-house MATLAB code is developed to couple SolTrace and ANSYS Fluent. CFD modeling is performed using ANSYS Fluent to predict the thermal performance of the receiver by evaluating radiation and convection heat loss mechanisms. Understanding the effects of variation in heliostat aiming strategy and flow configurations on the thermal performance of the receiver was achieved through parametric analyses. Finally, a receiver thermal efficiency ~85% was predicted and the surface temperatures were observed to be within the allowable limit for the materials under consideration.

54 ENVIRONMENTAL SCIENCES↗

Lagrange-Remap strategy for multi-material fluid-solid simulations using compressive limiters

In the present work, the Lagrange-Remap strategy proposed in [1] is extended to multi-material fluid-solid simulations. Both hypo-elastic and hyper-elastic material models are considered to describe the mechanical behavior of the solids. In practice, the deviatoric stress tensor (for hypo-elastic materials) and the left Cauchy-Green tensor (for isotropic hyper-elastic materials) are remapped, while the use of compressive limiters effectively reduces numerical diffusion during the remapping step. The simplicity of this diffuse interface approach is emphasized in the context of multi-material fluid-solid simulations. A series of Lagrange-Remap test cases, involving both solids and fluids, are conducted and compared with reference Lagrangian simulations, demonstrating the robustness and accuracy of the overall numerical strategy.

Compressive limiters↗

Fluid bed adsorption of carbon dioxide on immobilized polyethyenimine (PEI): kinetic analysis and breakthrough behavior

The adsorption of carbon dioxide (CO{sub 2}) by immobilized polyethylenimine (PEI) on mesoporous silica was investigated in a fluid bed. The tests were performed to determine breakthrough behavior with varying bed temperature, flow rates and feed concentrations. Experimental breakthrough curves were analyzed using a theoretical 1D model developed by Bohart and Adams. The results showed that Bohart-Adams model was suitable for the normal description of breakthrough curve for the temperature ranges of 40-90{degree}C. The maximum capacity increased with temperature up to 70{degree}C and then decreased. The adsorption rate constant exhibited a negative temperature dependence decreasing as the temperature increased. Parameters characteristic of a fluid bed adsorber were inferred from these breakthrough curves including the breakthrough time, saturation time, critical reactor length, and length of mass transfer zone LMTZ. These parameters can be used to design fluid bed adsorption system without resolving the mechanistic contributions of dispersion, mixing, and intraparticle diffusion.

Breakthrough↗

Quantification of the REE 3+ aqua ions and chloride species in aqueous fluids by in situ Raman spectroscopy using perturbations of the water band

Acidic NaCl-rich aqueous fluids play a crucial role in forming hydrothermal rare earth elements (REE) mineral deposits. Aqueous REE mobility is mostly controlled by the stabilities of REE 3+ and REE chloride species. Our current knowledge of REE speciation is based on solubility data, thermodynamic models and in situ spectroscopic measurements, sometimes coupled with molecular simulations. Here, in this study, we investigate Nd and Yb speciation in pH2 Cl-bearing solutions at 25 °C and 0.1 MPa with variable Cl/REE ratios using Raman Spectroscopy in solutions with 0.1 to 0.6 mol/kg NdCl 3 or YbCl 3 and 0.2 to 3.2 mol/kg NaCl. Due to the challenges in resolving the REE-Cl band, we developed a new method using the water vibrational mode and multivariate curve resolution (MCR) analysis. The Raman spectra for the vibrational band of water (2700 to 3900 cm –1 ) were collected at 25 °C and fitted by three Gaussian sub-peaks, then quantified using MCR analysis to de-convolute the water band into bulk H 2 O and the perturbations caused by of Cl – , REE 3+ , and REE chloride species. REE speciation based on the perturbations of the water band indicates that REE 3+ aqua ions dominate acidic solutions at 25 °C, but up to ~20 mol% YbCl 2+ forms at high YbCl 3 concentrations. The new method is promising for quantifying in situ speciation of the REE 3+ aqua ions and REE chloride species in aqueous fluids while providing information on the hydration of ions. This method improves our molecular level understanding of REE aqueous species stability and their role in REE mobilization during fluid-rock interaction.

58 GEOSCIENCES↗

A thermodynamically consistent discretization of 1D thermal-fluid models using their metriplectic 4-bracket structure

Thermodynamically consistent models in continuum physics, i.e. models which satisfy the first and second laws of thermodynamics, may be expressed using the metriplectic formalism. In this work, we leverage the structures underlying this modeling formalism to preserve thermodynamic consistency in discretizations of a fluid model. The procedure relies (1) on ensuring that the spatial semi-discretization retains certain symmetries and degeneracies of the Poisson and metriplectic 4-brackets, and (2) on the use of an appropriate energy conserving time-stepping method. Here, the minimally simple yet nontrivial example of a one-dimensional thermal-fluid model is treated. It is found that preservation of the requisite symmetries and degeneracies of the 4-bracket is relatively simple to ensure in Galerkin spatial discretizations, suggesting a path forward for thermodynamically consistent discretizations of more complex fluid models using more specialized Galerkin methods.

Hamiltonian structure↗

Extension of Landau-fluid closure to weakly collisional plasma regime

Both collisionless and collisional Landau fluid (LF) closure have been developed and implemented in two-fluid plasma simulation framework BOUT++, taking full advantage of the fast non-Fourier method (Dimits, 2014). The spectral range of good fit can be conveniently extended to larger domain by adding more Lorentzians. The scaling for upper limit of wavenumber resolved in collisionless LF operator and the fitting coefficients for the collisional LF operator, the later is extended to weakly collisional regime (ν$^{*}_{e}$ ~ 0.01), are presented in this paper. The flux limited expression q$^{FL}_{∥e}$ is reduced to free streaming expression q$^{FS}_{∥e}$ in collisionless limit and gives classical Spitzer–Härm result q$^{SH}_{∥e}$ in collisional limit as we expected. The collisional LF closure recovers the classical Spitzer–Härm result q$^{SH}_{∥e}$ in collisional limit and is same as the collisionless LF closure in collisionless limit. The test cases also obviously reflect the nonlocal effects from LF closures and responding to different boundary conditions in open surface region. Our report will make it more applicable and possible to include the Landau damping kinetic effect in the fluid simulation of tokamak plasma.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Heat Transfer and Fluid Transport of Supercritical CO 2 in Enhanced Geothermal System with Local Thermal Non-equilibrium Model

The heat transfer and fluid transport of supercritical CO 2 in enhanced geothermal system (EGS) is studied numerically with local thermal non-equilibrium model, which accounts for the temperature difference between solid matrix and fluid components in porous media and uses two energy equations to describe heat transfer in the solid matrix and in the fluid, respectively. As compared with the previous results of our research group, the effect of local thermal non-equilibrium mainly depends on the volumetric heat transfer coefficient ah, which has a significant effect on the production temperature at reservoir outlet and thermal breakthrough time. The uniformity of volumetric heat transfer coefficient ah has little influence on the thermal breakthrough time, but the temperature difference become more obvious with time after thermal breakthrough with this simulation model. The thermal breakthrough time reduces and the effect of local thermal non-equilibrium becomes significant with decreasing ah.

58 GEOSCIENCES↗

Integrated experimental studies of pore structure and fluid uptake in the Bossier Shale in eastern Texas, USA

Within the Haynesville-Bossier Shale complex, the Bossier Shale has not been extensively studied by either industry and academia, despite it being an unconventional gas reservoir and a potential caprock for carbon storage in the underlaying Haynesville Shale. The lack of knowledge of the complex pore structures and fluid-rock interactions hinders the effective extraction of gas and the characterization of fluid reservoirs and sealing capacity. Integrated experimental studies of pore structure and fluid-rock interactions were conducted in seven Bossier Shale core samples collected in eastern Texas. Petrographic, geochemical, and petrophysical properties such as mineral composition, organic richness, thermal maturity, porosity, pore/pore throat diameter distribution, water-accessible pores, liquid water imbibition, and water vapor adsorption were characterized using complementary approaches of thin-section petrography, scanning electron microscopy, X-ray diffraction, total organic matter, pyrolysis, mercury intrusion porosimetry, nuclear magnetic resonance, (Ultra-) small angle X-rays scattering as well as small angle neutron scattering with deuterated liquids and contrast variation. Further, the results show that the thermally mature Bossier Shales are composed of mixed argillaceous mudstone, mixed mudstone, and mixed carbonate mudstone. The shale contains both organic and inorganic pores, with porosities of 3.24-9.37 %, pore-to-throat ratios of 1.65 to 19.4, and water-accessible pores accounting for 28.7-72.6 % of total pores. Approaches of liquid water imbibition and water vapor adsorption, with and without direct contact of water with shale samples, indicate that liquid water first enters the nano-sized pores under high capillary pressures, and water vapor adsorption is mainly controlled by both clay minerals and pores with diameters less than 10 nm. These findings contribute to a better understanding of pore structures and water-shale interactions and their controlling factors in the Bossier Shale.

58 GEOSCIENCES↗

An oscillation free shock-capturing method for compressible van der Waals supercritical fluid flows

Numerical solutions of the Euler equations using real gas equations of state (EOS) often exhibit serious inaccuracies. The focus here is the van der Waals EOS and its variants (often used in supercritical fluid computations). The problems are not related to a lack of convexity of the EOS since the EOS are considered in their domain of convexity at any mesh point and at any time. The difficulties appear as soon as a density discontinuity is present with the rest of the fluid in mechanical equilibrium and typically result in spurious pressure and velocity oscillations. This is reminiscent of well-known pressure oscillations occurring with ideal gas mixtures when a mass fraction discontinuity is present, which can be interpreted as a discontinuity in the EOS parameters. We are concerned with pressure oscillations that appear just for a single fluid each time a density discontinuity is present. As a result, the combination of density in a nonlinear fashion in the EOS with diffusion by the numerical method results in violation of mechanical equilibrium conditions which are not easy to eliminate, even under grid refinement.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Thermodynamical effects and high resolution methods for compressible fluid flows

One of the fundamental differences of compressible fluid flows from incompressible fluid flows is the involvement of thermodynamics. This difference should be manifested in the design of numerical schemes. Unfortunately, the role of entropy, expressing irreversibility, is often neglected even though the entropy inequality, as a conceptual derivative, is verified for some first order schemes. In this paper, we refine the GRP solver to illustrate how the thermodynamical variation is integrated into the design of high resolution methods for compressible fluid flows and demonstrate numerically the importance of thermodynamic effects in the resolution of strong waves. As a by-product, we show that the GRP solver works for generic equations of state, and is independent of technical arguments.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

A solution algorithm for fluid–particle flows across all flow regimes

Many fluid–particle flows occurring in nature and in technological applications exhibit large variations in the local particle volume fraction. For example, in circulating fluidized beds there are regions where the particles are closepacked as well as very dilute regions where particle–particle collisions are rare. Thus, in order to simulate such fluid–particle systems, it is necessary to design a flow solver that can accurately treat all flow regimes occurring simultaneously in the same flow domain. In this work, a solution algorithm is proposed for this purpose. The algorithm is based on splitting the free-transport flux solver dynamically and locally in the flow. In close-packed to moderately dense regions, a hydrodynamic solver is employed, while in dilute to very dilute regions a kinetic-based finite-volume solver is used in conjunction with quadrature-based moment methods. To illustrate the accuracy and robustness of the proposed solution algorithm, it is implemented in OpenFOAM for particle velocity moments up to second order, and applied to simulate gravity-driven, gas–particle flows exhibiting cluster-induced turbulence. By varying the average particle volume fraction in the flow domain, it is demonstrated that the flow solver can handle seamlessly all flow regimes present in fluid–particle flows.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Modeling of wave processes in blocky media with porous and fluid-saturated interlayers

The wave processes in blocky media are analyzed by applying different mathematical models, wherein the elastic blocks interact with each other via pliant interlayers with the complex mechanical properties. Four versions of constitutive equations are considered. In the first version, an elastic interaction between the blocks is simulated within the framework of linear elasticity theory, and the model of elastic–plastic interlayers is constructed to take into account the appearance of irreversible deformation of interlayers at short time intervals. In the second one, the effects of viscoelastic shear in the interblock interlayers are taken into the consideration using the Poynting–Thomson rheological scheme. In the third option, the model of an elastic porous material is used in the interlayers, where the pores collapse if an abrupt compressive stress is applied. In the fourth case, the model of a fluid-saturated material with open pores is examined based on Biot's equations. The collapse of pores is modeled by the generalized rheological approach, wherein the mechanical properties of a material are simulated using four rheological elements. Three of them are the traditional elastic, viscous and plastic elements, the fourth element is the so-called rigid contact, which is used to describe the behavior of materials with the different resistance to tension and compression. It was shown that the thermodynamically consistent model is provided, which means that the energy balance equation is fulfilled for an entire blocky structure, where the kinetic and potential energy of the system is the sum of the kinetic and potential energies of the blocks and interlayers. Under numerical implementation of the interlayers models, the dissipationless finite difference Ivanov's method was used. The splitting method by spatial variables in the combination with the Godunov gap decay scheme was applied in the blocks. As a result, robust and stable computational algorithms are built and tested. Using MPI technology, the parallel software was designed for the modeling of wave processes in 2D setting. The numerical results are presented, discussed and future studies are outlined. - Highlights: • We model the propagation of waves in a blocky medium with compliant interlayers. • Models of porous and fluid saturated interlayers are proposed using Biot's approach. • A parallel computational algorithm is developed for supercomputers of cluster type. • Problems of instant rotation of central and boundary blocks are solved numerically. • Wave motions are compared in the cases of elastic and fluid-saturated interlayers.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

A Cartesian-diffusion Langevin method for hybrid kinetic-fluid Coulomb scattering in particle-in-cell plasma simulations

A novel, drag-diffusion Langevin method of hybrid, kinetic-fluid Coulomb scattering in plasmas is presented. Unlike previous methods, the frictional drag is always applied in the simulation frame of reference. The velocity-space diffusion is performed in the stationary-fluid frame of reference when anisotropic, and in the laboratory frame when isotropic. While the general method is mass-ratio independent, we focus on interactions of kinetic-ions and fluid-electrons to show first-order modifications to the electron velocity distribution function that are an important correction for the accurate calculation of electric resistivity. Inclusion of sub-cycling and a limit to the maximum collision frequency is shown to allow for arbitrarily large timesteps without numerical failure. Here the Langevin method is compared with a grid-based binary method and found to require a much less restrictive timestep in cases of ion–electron slowing and temperature equilibration; this finding differs from previous work and is dependent on the mass ratio.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Characterization of the mud displacement in an enlarged wellbore: An integrated rock-fluid model

Cement-mud displacement plays a crucial role in the sealability of cement sheaths. Irregular geometric features of a wellbore due to washout can have a negative impact on mud and cement mobilization. An unstable interface between two fluids always leads to mud channeling, interfluid mixing, and cement contamination, degrading the cement quality. Many factors, such as mechanical and rheological properties of fluids, annulus geometry, flow pattern, and flow rate, significantly influence the displacement efficiency. This study investigates the characterization of the mud displacement in an irregular horizontal well using a 3D computational fluid dynamics (CFD) model. Mud is displaced in an enlarged wellbore by geopolymer and neat class G cement. The wellbore geometry is developed based on the caliper log data from an unconventional shale well in the Tuscaloosa Marine Shale (TMS) lithology. The effects of pump rate, density difference, and mud contamination are evaluated by numerical simulations. The results present those residual muds mainly exist in the upper annulus of the enlarged section. Geopolymer has a better sealing performance and can resist more water-based mud (WBM) contaminations than neat class G cement. The scenario with a low mud-cement density difference and high cement injection rate results in a high cement volume fraction, mitigating the gas migration.

58 GEOSCIENCES↗

A process to verify numerical models for seismic fluid-structure interaction in advanced reactor vessels

Seismic design and qualification of a liquid-filled advanced nuclear reactor will have to account for fluid-structure interaction (FSI). Interaction between the tank, internal components, and contained liquid will rely on analysis of numerical models that must be verified and validated. Here this study demonstrates a verification process for models of a base-supported cylindrical tank by comparing numerical predictions and analytical solutions. The numerical models are consistent with the assumptions made to derive analytical solutions, namely, either a rigid or a linear elastic tank, ideal fluid, and small-amplitude, unidirectional, horizontal inputs. One software platform is used to illustrate the process. Seismic FSI analysis is performed using the Arbitrary Lagrangian-Eulerian (ALE) and Incompressible Computational Fluid Dynamics (ICFD) solvers in LS-DYNA. Reported responses are those used for design, including hydrodynamic pressures on the tank wall, shear forces and moments at the tank base, and wave heights of the contained liquid. The accuracy of the numerical results is discussed. The numerical models are verified for calculating the pressures on the tank wall and reactions at its base. Accurate simulation of wave action is challenging for both solvers. Recommendations for modeling, code development, and steps for verification are provided. Although focused on reactor vessels and one software platform, the verification process described herein is broadly applicable to liquid-filled vessels and other finite element codes.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Local and global Λ polarization in a vortical fluid

We compute the fermion spin distribution in the vortical fluid created in off-central high energy heavy-ion collisions. We employ the event-by-event (3+1)D viscous hydrodynamic model. The spin polarization density is proportional to the local fluid vorticity in quantum kinetic theory. As a result of strong collectivity, the spatial distribution of the local vorticity on the freeze-out hyper-surface strongly correlates to the rapidity and azimuthal angle distribution of fermion spins. We investigate the sensitivity of the local polarization to the initial fluid velocity in the hydrodynamic model and compute the global polarization of Λ hyperons by the AMPT model. The energy dependence of the global polarization agrees with the STAR data.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Device-scale computational fluid dynamics modeling of carbon dioxide absorption using encapsulated sorbents

Micro-encapsulated carbon dioxide sorbent (MECS) is considered to be a promising material for carbon capture system because of high capture rate, low fabrication cost, and efficient energy consumption. To accelerate the development of MECS technology for industrial deployment, device-scale computational fluid dynamics (CFD) simulations are performed to study hydrodynamics and CO2 absorption behavior in a conceptual fluidized bed reactor. The two-fluid model solver, built in the Multiphase Flow with Interphase eXchanges (MFIX), is adopted to solve the multiphase flow hydrodynamics by representing both fluid and solids as continuous phases. Filtered models are implemented to capture the sub-grid effects that normally cannot be directly considered because MECS particles are too small to be explicitly resolved at such a large scale. Using established CFD tools, this study aims to characterize MECS behavior and predict its carbon capture performance in a conceptual device-scale absorber. The MECS particle size and gas flow rate have significant impact on CO2 capture efficiency. The CO2 capture efficiency decreases with increasing gas flow rate. Smaller particle size lead to better gas-solid homogeneity, as a results CO2 absorption efficiency increases.

Wang, Chao (ORCID:0000000279406614)↗

From 2D to 3D in fluid turbulence: unexpected critical transitions

How do the laws of physics change with changes in spatial dimension? Maybe not at all in some cases, but in important cases, the changes are dramatic. Fluid turbulence – the fluctuating, intermittent and many-degree-of-freedom state of a highly forced fluid – determines the transport of heat, mass and momentum and is ubiquitous in nature, where turbulence is found on spatial scales from microns to millions of kilometres (turbulence in stars) and beyond (galactic events such as supernovae). When the turbulent degrees of freedom are suppressed in one spatial dimension, the resulting turbulent state in two dimensions (2D) is remarkably changed compared with the turbulence in three dimensions (3D) – energy flows to small scales in 3D but towards large scales in 2D. Although this result has been known since the 1960s due to the pioneering work of Kraichnan, Batchelor and Leith, how one transitions between 3D and 2D turbulence has remained remarkably unexplored. For real physical systems, this is a highly significant question with important implications about transport in geophysical systems that determine weather on short time scales and climate on longer scales. Is the transition from 3D to 2D smooth or are there sharp transitions that signal a threshold of the dominance of one type of turbulence over another? Finally, recent results by Benavides & Alexakis (J. Fluid Mech., vol. 822 (2017), pp. 364–385) suggest that the latter may be the case – a surprising and provocative discovery.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗