Search NASA⌕ Search

SEARCH · Search NASA

Results for “viscous fluid”

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 127 records · Page 7
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Fluid dynamic and thermal performance of a slotted cylinder at low Reynolds number

The fluid dynamic and thermal performance of a circular cylinder with a slot parallel to the flow is numerically investigated. The study utilized the semi-implicit finite volume multi-material algorithm MPM-ICE, a component of the Uintah framework. The normalized slot width s/D ranges from 0.1 - 0.3, introducing an additional heat transfer surface area between ~ 10 and ~ 50%, and a mass reduction between ~ 13 and ~ 38% in the cylinder. We assumed two-dimensional incompressible flow and simulated a Reynolds number Re D between 100 and 1000. The slotted cylinders are found to have a total drag force reduction up to ~ 45%, compared to a solid cylinder despite the additional viscous drag force in the slot. Convection heat transfer is enhanced up to ~ 70%. Further, the slotted cylinder performance index, defined as the ratio of the heat rate to the drag force, increases up to maximum of ~ 3, indicating better overall thermal fluid performance. An entropy analysis showed the best performance index occurs at the highest Re D . Correlations for drag coefficient and Nusselt number are proposed along with an entropy optimization method.

42 ENGINEERING↗
DOE OSTI2022

Quantifying transport and electrocatalytic reaction processes in a gastight rotating cylinder electrode reactor via integration of Computational Fluid Dynamics modeling and experiments

Understanding the complexity of the multiple processes of mass, momentum, charge, and heat transport, and how these affect reaction kinetics at the electrode/electrolyte interface is one of the major challenges in the field of energy and catalysis. The rapid and rational scale-up of electrocatalytic systems to industrial scales require a detailed understanding of nonlinear transport-reaction processes, accessible only through the building of multi-physics models that capture with high fidelity the complexity of real-world devices. The gastight rotating cylinder electrode (RCE) reactor is a promising lab-scale tool that can decouple transport from intrinsic kinetics to generate data for first-principle models useful in the design of industrial, electrochemical reactors. Computational Fluid Dynamics (CFD) studies have previously been used to investigate the bulk flow in RCE reactors for simple corrosion and electroplating processes. However, the quantification of changes in local concentration within the viscous layer where catalysis takes place requires capturing the correct flow conditions inside the hydrodynamic boundary layer near the surface of the electrode. Further, this requires simulations with spatial resolution in the nm and μm scale and temporal resolutions between ms and s scales that are similar to the timescales for reactions on the electrode surface. In this study, experimental electrocatalysis is combined with CFD modeling to elucidate and parameterize the hydrodynamics in a gastight RCE reactor. CFD simulations of the electrochemical ferricyanide reduction reaction under mass transport limited conditions are used to evaluate the validity of the CFD model parameters by comparing calculated dimensionless mass transport descriptors to dimensionless correlations obtained experimentally. Justifications for assumptions and details of the simulation methods used in this study are presented to provide a detailed understanding of the effect that each model parameter has on the ability to accurately simulate electrocatalysis in RCE systems. The simulation methodology reported here is a first step towards the development of multi-scale models for the study of transport dependent electrocatalytic processes, such as the electrochemical transformation of CO 2 to fuels and chemicals.

42 ENGINEERING↗
DOE OSTI2026

Fragile-to-strong transition in liquid As 2 S 3 under pressure: The effect of melt metallization

The well-known classification of glass-forming melts into fragile and strong liquids has several notable exceptions, including water, silica, and certain phase-change materials (PCMs). These exceptional fluid systems exhibit a fragile-to-strong transition (FST) upon cooling: a transformation from a high-temperature liquid with fast atomic dynamics, low viscosity, and low flow activation energy, to a viscous supercooled melt with high energy barriers near the glass transition temperature T g . This behavior is critically important for non-volatile memories, photonic tensor cores, reconfigurable metamaterials, and other devices, that use PCMs, enabling nanosecond-scale crystallization in the fragile regime and long data retention in the strong regime near or below T g . A significant structural transformation is expected between these two viscosity regimes, along with a semiconductor-metal (SC-M) transition upon heating, driven by high internal pressure and associated density increase. By applying high external pressure to the canonical low-conducting chalcogenide melt As 2 S 3 , we observed both the FST and the SC-M transition, occurring simultaneously within the same domain of the P, T−phase space. These findings suggest that the FST is not limited to a few exceptional liquids but is a common phenomenon, at least in systems that exhibit melt metallization within specific regions of their P, T−phase diagrams.

first-principles molecular dynamics↗
DOE OSTI2021

Optimal surface-tension isotropy in the Rothman-Keller color-gradient lattice Boltzmann method for multiphase flow

The Rothman-Keller color-gradient (CG) lattice Boltzmann method is a popular method to simulate two-phase flow because of its ability to deal with fluids with large viscosity contrasts and a wide range of interfacial tensions. Here, two fluids are labeled red and blue, and the gradient in the color difference is used to compute the effect of interfacial tension. It is well known that finite-difference errors in the color-gradient calculation lead to anisotropy of interfacial tension and errors such as spurious currents. Here, we investigate the accuracy of the CG calculation for interfaces between fluids with several radii of curvature and find that the standard CG calculations lead to significant inaccuracy. Specifically, we observe significant anisotropy of the color gradient of order 7% for high curvature of an interface such as when a pinchout occurs. We derive a second order accurate color gradient and find that the diagonal nearest neighbors can be weighted differently than in the usual color-gradient calculation such that anisotropy is minimized to a fraction of a percent. The optimal weights that minimize anisotropy for the smallest radius of curvature interface are found to be w = (0.298, 0.284, 0.275) for diagonal nearest neighbors for the cases of the interface smoothing parameter β = (0.5, 0.7, 0.99), somewhat higher than the w = 0.25 value derived by Leclaire et al. [Leclaire, Reggio, and Trepanier, Computers and Fluids 48, 98 (2011)] based on obtaining isotropic errors to second order. We find that use of these optimal w values yields over a factor of 10 decrease in anisotropy and over a factor of 30 decrease in mean anisotropy relative to using the standard w = 1 value. And we find a factor of about 2 decrease in the anisotropic error and up to factor 15 decrease in mean anisotropic error relative to the choice of w = 0.25 for small radius of curvature interfaces. The improved CG calculations will allow the method to be more reliably applied to studies of phenomenology and pore scale processes such as viscous and capillary fingering, and droplet formation where surface-tension isotropy of narrow fingers and small droplets plays a crucial role in correctly capturing phenomenology. We present an example illustrating how different phenomena can be captured using the improved color-gradient method. Namely, we present simulations of a wetting fluid invading a fluid filled pipe where the viscosity ratio of fluids is unity in which droplets form at the transition to fingering using the improved CG calculations that are not captured using the standard CG calculations. We present an explanation of why this is so which relates to anisotropy of the surface tension, which inhibits the pinchouts needed to form droplets.

58 GEOSCIENCES↗
DOE OSTI2020

Probing the Partonic Degrees of Freedom in High-Multiplicity p–Pb collisions at √ s NN = 5.02 TeV

We investigate the role of partonic degrees of freedom in high-multiplicity p–Pb collisions at √ s NN = 5.02 TeV carried out at the Large Hadron Collider (LHC) by studying the production and collective flow of identified hadrons at intermediate p T via the coalescence of soft partons from the viscous hydrodynamics (VISH 2 + 1) and hard partons from the energy loss model, linear Boltzmann transport (LBT). We find that combining these intermediate p T hadrons with the low p T hadrons from the hydrodynamically expanding fluid and high p T hadrons from the fragmentation of quenched jets, the resulting hydro-dynamics-coalescence-fragmentation model provides a nice description of measured p T spectra and differential elliptic flow v 2 (p T ) of pions, kaons, and protons over the p T range from 0 to 6 GeV. We further demonstrate the necessity of including the quark coalescence contribution to reproduce the experimentally observed approximate number of constituent quark scaling of hadron v 2 at intermediate p T . Furthermore, our results thus indicate the importance of partonic degrees of freedom and also hint at the possible formation of quark-gluon plasma in high-multiplicity p–Pb collisions at the LHC.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2021

Activated penetrant dynamics in glass forming liquids: size effects, decoupling, slaving, collective elasticity and correlation with matrix compressibility

In this work, we employ the microscopic self-consistent cooperative hopping theory of penetrant activated dynamics in glass forming viscous liquids and colloidal suspensions to address new questions over a wide range of high matrix packing fractions and penetrant-to-matrix particle size ratios. The focus is on the mean activated relaxation time of smaller tracers in a hard sphere fluid of larger particle matrices. This quantity also determines the penetrant diffusion constant and connects directly with the structural relaxation time probed in an incoherent dynamic structure factor measurement. The timescale of the non-activated fast dissipative process is also studied and is predicted to follow power laws with the contact value of the penetrant–matrix pair correlation function and the penetrant–matrix size ratio. For long time penetrant relaxation, in the relatively lower packing fraction metastable regime the local cage barriers are dominant and matrix collective elasticity effects unimportant. As packing fraction and/or penetrant size grows, much higher barriers emerge and the collective elasticity associated with the correlated matrix dynamic displacement that facilitates penetrant hopping becomes important. This results in a non-monotonic variation with packing fraction of the degree of decoupling between the matrix and penetrant alpha relaxation times. The conditions required for penetrant hopping to become slaved to the matrix alpha process are determined, which depend mainly on the penetrant to matrix particle size ratio. By analyzing the absolute and relative importance of the cage and elastic barriers we establish a mechanistic understanding of the origin of the predicted exponential growth of the penetrant hopping time with size ratio predicted at very high packing fractions. A dynamics-thermodynamics power law connection between the penetrant activation barrier and the matrix dimensionless compressibility is established as a prediction of theory, with different scaling exponents depending on whether matrix collective elasticity effects are important. Quantitative comparisons with simulations of the penetrant relaxation time, diffusion constant, and transient localization length of tracers in dense colloidal suspensions and cold viscous liquids reveal good agreements. Multiple new predictions are made that are testable via future experiments and simulations. Extension of the theoretical approach to more complex systems of high experimental interest (nonspherical molecules, semiflexible polymers, crosslinked networks) interacting via variable hard or soft repulsions and/or short range attractions is possible, including under external deformation.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2025

Scanning Electrochemical Microscopy for Kinetic Investigations in Viscous Deep Eutectic Solvents: Identifying Practical Approach Curves and Deviations from Electron Transfer Models

Determining heterogeneous electrochemical electron transfer (ET) kinetics in electrolytes with a wide range of physical properties is of great interest for achieving high-performance redox flow batteries. Among such electrolytes, concentrated hydrogen-bonded electrolytes (CoHBEs), including deep eutectic solvents (DESs), have recently garnered significant attention. Unfortunately, traditional Tafel analysis using macroelectrodes often encounters issues with mass transfer limitations in CoHBEs with high viscosities, thereby restricting kinetic analysis to a narrow potential window. Here, in this work, we introduce a methodology for evaluating ET kinetics in viscous DES using the scanning electrochemical microscopy (SECM). We first determined practical solutions to SECM tip positioning in ethaline DES, which yield pseudopositive feedback responses. Lattice Boltzmann method (LBM) simulations helped us rationalize the impact of the fluid and concentration fields, as well as tip geometry, tip approach velocity v, and the solvent viscosity ηs, on the shape of the approach curves. In addition to successfully recreating approach curves over a variety of conditions, we found that approaching a conductor ensured a practical point where the normalized tip response (Ni T = 2) converged at L = 0.7 within ∼10% error regardless of tip velocity. With positioning capabilities at hand, we investigated the kinetics of Fe 3+ /Fe 2+ redox couple in aqueous and the ethaline media. The experimental kinetic results were interpreted using the Butler–Volmer (BV) and Marcus–Hush–Chidsey (MHC) models. For ethaline, a nonideal kinetic behavior was observed, potentially attributed to solvent dynamics within DESs or to the interplay of chloride anions in the charge transfer process.

electrodes↗
DOE OSTI2020

High pressure/high temperature multiphase simulations of dodecane injection to nitrogen: Application on ECN Spray-A

The present work investigates the complex phenomena associated with pressure/high temperature dodecane injection for the Engine Combustion Network (ECN) Spray-A case, employing more elaborate thermodynamic closures, to avoid well known deficiencies concerning density and speed of sound prediction using traditional cubic models. A tabulated thermodynamic approach is proposed here, based on log 10 (p)-T tables, providing very high accuracy across a large range of pressures, spanning from 0 to 2500 bar, with only a small number of interpolation points. The tabulation approach is directly extensible to any thermodynamic model, existing or to be developed in the future. Here NIST REFPROP properties are used, combined with PC-SAFT Vapor-Liquid-Equilibrium to identify the liquid in mixtures penetration, hence avoiding the use of an arbitrary threshold for mass fraction. Identified liquid and vapour penetration are compared against experimental data from the ECN database showing a good agreement, within approximately 3–8% for axial penetration of liquid, 2% for vapor axial penetration and within experimental uncertainty for radial distribution of mass fraction. Analysis of the vortex evolution indicates that driving mechanisms behind the jet break-up are vortex tilting/stretching, then baroclinic torque, leading to Rayleigh-Taylor instabilities, closely followed by vortex dilation and finally viscous effects.

42 ENGINEERING↗
DOE OSTI2023

Impact of Flow Configurations on Response Time and Data Quality in Real-Time, In-Line Fourier Transform Infrared (FTIR) Monitoring of Viscous Flows

The real-time, in-line monitoring of continuous flow concentrations is widely conducted via infrared (IR) spectrometry by using a flow cell connected to a reactive flow stream. For protective purposes, the IR sensor tip is typically offset from the flow. This offset can cause the formation of a stagnant boundary layer above the sensor, especially when dealing with high-viscosity fluids. As a result, the IR signal response time is often controlled by the slow diffusional exchange of fluid in the boundary layer, as confirmed via 2D computational fluid dynamics (CFD) simulations. We evaluated several flow configuration modifications in a typical IR flow cell in order to identify the changes to the flow dynamics that enable improved response times with minimal changes to the cell configuration: the use of (i) vertical flow, where the standard horizontal flow over the sensor is redirected to contact vertically with the sensor, (ii) a static mixer to create radial flow momentum above the IR sensor, and (iii) horizontal or vertical nozzles to direct the flow toward the IR sensor. The vertical flow configuration did not show any significant improvement over the standard horizontal flow configuration. However, the static mixer, horizontal nozzle, and vertical nozzle configurations all resulted in markedly improved response times and signal quality, albeit at the expense of a higher pressure drop across the flow cell. These results point toward straightforward, user-accessible modifications of in-line IR flow cells that result in significant improvements in signal stability and acquisition times.

42 ENGINEERING↗
DOE OSTI2022

A mass conservative, well balanced, tangency preserving and energy decaying method for the shallow water equations on a sphere

Here a fully discrete surface finite element method is proposed for solving the viscous shallow water equations in a bounded Lipschitz domain on the sphere based on a general triangular mesh. The method consists of a modified Crank–Nicolson method in time and a Galerkin surface finite element method in space for the fluid thickness H and the fluid velocity u. A finite element space tangential to the sphere at all finite element nodes is proposed to approximate the fluid velocity u. The proposed method has second-order accuracy in time and first-order accuracy in space, and preserves mass conservation, well balancedness, tangency of velocity to the sphere, and energy decay. Numerical experiments are presented to illustrate the accuracy of the proposed method and the preservation of the physical properties, including mass conservation, well balancedness, and energy decay. A numerical simulation of ocean mesoscale activity on a circular basin with a continental shelf is provided.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Pressure drop in a prototypical 3D magnetohydrodynamic flow across contraction of a fusion blanket manifold

Predicting 3D magnetohydrodynamic (MHD) pressure losses associated with liquid metal flows in complex geometry ducts is required to design breeding blankets for fusion reactors. Of such components, manifolds exhibit major pressure losses. Recent correlations [T. Rhodeset al. Magnetohydrodynamic pressure drop and flow balancing of liquid metal flow in a prototypic fusion blanket manifold. Phys. Fluids. 2018; 30: 057101.] demonstrate promise for predicting pressure drops in electrically insulated inlet manifolds. In the present work, the extend to which these correlations can be applied to outlet manifolds, which feature sudden contractions, is investigated in both viscous-electromagnetic (VE) and inertial-electromagnetic (IE) regimes for Reynolds numbers 50 < Re < 1500, Hartmann numbers 2500 < Ha < 5475, and expansion/contraction ratio 10. Numerical computations have shown that the MHD pressure drops in the eutectic lead-lithium (PbLi) flows in contractions are almost identical to those in the flows featuring expansion with slightly higher magnitudes in the contraction cases. The small discripancy in the MHD pressure drop between contractions and expansions (<8%) suggests that the earlier obtained correlations for the 3D MHD pressure drop in a duct flow with a sudden expansion can also be applied to flows with a sudden contraction.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Coupled confined phase behavior and transport of methane in slit nanopores

The proximity of the order of magnitude molecular mean free path and pore sizes in nanopores leads to remarkable interactions between molecules and walls. In such systems, the thermodynamic property and the transport behavior of fluids deviate from those at bulk conditions. Molecular dynamics (MD) simulation may be used to investigate the effects of confinement on fluid physics in nanosize pores. However, translating these subscale observations into a larger scale of interest remains a challenging task. In this work, we propose a modified extension of Peng-Robinson equation of state (PR-EOS) motivated and guided by MD simulation results. The shift of critical properties, i.e., pressure and temperature, are evaluated independently. A temperature dependent parameter is introduced to account for the existence of capillary pressure in the two-phase region. Here, this formulation is capable of capturing the shift in critical properties as well as density phase diagrams under various confinement scenarios. We incorporate the proposed EOS in the lattice Boltzmann method (LBM) to study coupled confined phase behavior and transport of methane in nano-size slit pores. Adsorptive strength is determined such that the resulting density ratio matches that from MD simulation. Results indicate that the use of the proposed EOS leads to smaller fluid velocities compared to PR-EOS. This effect is due to stronger interactions between fluid particles under confinement. Also, results show that transport characteristics are impacted by the pressure of the system: in systems with a relatively low pressure, transport seems to be dominated by Knudsen diffusion; however, at a higher pressure, the contribution from viscous flow increases as the pore widens while the influence of Knudsen and surface diffusion diminishes.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2023

Novel features of attractors and transseries in nonconformal Bjorken flows

In this work we investigate the impact of conformal symmetry breaking on hydrodynamization of a far-from-equilibrium fluid. We find a new kind of transseries solutions for the nonconformal hydrodynamic equations of a longitudinal boost invariant expanding plasma. The new transseries solutions unveil a rich physical structure which arises due to the interplay of different physical scales. In the perfect fluid case the nonconformal speed of sound slows down the cooling of the temperature due to the emergence of logarithmic corrections that depends on the mass of the particle. Further, these terms propagate into the perturbative and nonperturbative sectors of the transseries once viscous corrections are included. The logarithmic mass contributions increase the asymptotic value of the Knudsen number while decreasing the damping rate of the transient nonhydrodynamic modes and thus, yielding to an extremely slow hydrodynamization process where flow lines merge to their forward attractor at extremely late times. The early time free streaming expansion is modified and receives logarithmic mass corrections induced by the shear-bulk couplings. The global flow structure and numerical analyses carried out in our work demonstrate the existence of local attraction at early and late times for the shear viscous tensor and bulk viscous pressure.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2023

Helical Fins for Concentrated Solar Receivers: Design Optimization and Entropy Analysis

Concentrated solar power (CSP) with thermal energy storage (TES) has the potential to achieve grid parity. This can be realized by operating CSP systems at temperatures above 700 °C with high-efficiency sCO 2 power cycles. However, operating CSP systems at such temperatures poses several challenges, among which the design of solar receivers to accommodate increased thermal loads is critical. To this end, this work explores and optimizes various swirl-inducing internal fin designs for solar receiver tubes. These fin designs not only improve the thermal performance of receiver tubes but also levelize temperature unevenness caused by non-uniform thermal loading. In this work, the geometric parameters of the fin designs are optimized to maximize the Nusselt number with a constraint on the friction factor. This optimization, however, is computationally intensive, requiring hundreds of simulation calls to computational fluid dynamics (CFD) models. Here, to circumvent this problem, surrogate models are used to approximate the simulation outputs needed during the optimization. In addition, this study also examines the fin designs from an entropy generation perspective. To this end, the entropy contributions from thermal and viscous effects are quantitatively compared while varying the operational Reynolds number.

14 SOLAR ENERGY↗
DOE OSTI2020

A phenomenological rate of injection model for predicting fuel injection with application to mixture formation in light-duty diesel engines

Fuel injection rate laws are one of the most important pieces of information needed when modeling engine combustion with computational fluid dynamics. In this study, a simple phenomenological model of a common-rail injector was developed and calibrated for the Bosch CRI2.2 platform. The model requires three tunable parameter fits, making it relatively easy to calibrate and suitable for injector modeling when high-fidelity information about the internal injector’s geometry and electrical circuit details are not available. Each injection pulse is modeled as a sequence of up to four stages: an injection needle mechanical opening transient, a full-lift viscous flow inertial transient, a Bernoulli steady-state stage, and a needle descent transient. Parameters for each stage are obtained as polynomial fits from measured injection rate properties. The model enforces total injected mass, and the intermediate stages are only introduced if the injection pulse duration is long enough. Experimental rates of injection from two separate campaigns on the same injector were used to calibrate the model. The model was first validated against measured injection rate laws featuring pilot injections, short partially premixed combustion pulses, and conventional diesel combustion injection strategies. Then, it was employed as an input to engine computational fluid dynamics simulations, which were run to simulate experiments of mixture formation in an optically accessible light-duty diesel engine. It was found that, though simple, this model is capable of predicting both pilot and main injection pulse mass flow rates well: the simulations yielded accurate predictions of in-cylinder equivalence ratio distributions from injection strategies for both partially premixed combustion and pilot injections. Also, once calibrated, the model produced appropriate results for a wide range of injected mass and rail pressure values. Finally, it was observed that usage of such a relatively simple model can be a good choice when high-fidelity injection rate input and highly detailed information of the injector’s geometry and operation are not available, particularly as noticeable discrepancies can be present also among different experimental campaigns on similar hardware.

42 ENGINEERING↗
DOE OSTI2025

Development of fiber orientation in injection molding: Comparison of glass fiber, carbon fiber and their hybrid composites

Fiber orientation distribution (FOD) in injection-molded panels with respect to distance from the gate was analyzed using X-ray computed tomography (X-CT) for glass fiber (GF), carbon fiber (CF), and hybrid CF/GF (CGF) reinforced nylon 66. To understand the reason behind the FOD with different fiber types, computational fluid dynamics (CFD) and rheology were performed. Samples were extracted at three locations: near the gate, center, and opposite end. Thickness of the layers of typical skin-shell-core type FOD varied with fiber type and location. GF achieved flow direction alignment (in shell) earlier than viscous CF and CGF near the gate, whereas CF showed the highest flow-direction alignment at the center due to shear induced orientation. At the opposite end, GF experienced more backflow than others indicating faster mold filling owing to its lower viscosity. Hybrid CGF exhibited GF-dominated center and CF-dominated end region. The numerical model used to obtain FOD and rheological predictions for the CF and GF composites served to corroborate the trends observed in the experimental trials. The FOD responses across fiber types and location were reflected in their longitudinal and transverse properties. Only GF showed higher longitudinal modulus over transverse modulus near the gate attributed to rapid alignment, whereas CF and CGF exhibited opposite trend. However, fountain flow enhanced the longitudinal modulus over transverse modulus with the distance for all, particularly for CF. This study offers insights into mold filling behavior of different fibers which are critical in optimizing injection molding conditions for tailored final properties.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Design and Production for Maximum Structural Efficiency With Respect to Fiber Orientation With Increased Understanding of Hybrid Fiber Flow Behavior

Discontinuous fiber-reinforced thermoplastic composites have gained considerable attention in automotive, aerospace, and other industries, due to their high-rate of production combined with their ability to attain complex and intricate shapes. Among other high-rate thermoplastic manufacturing processes, injection-molding is one of the most common manufacturing methods due to fast production and high surface finishing of complex geometries. Fiber orientation in discontinuous fiber composites plays a pivotal role in determining the mechanical, electrical, and thermomechanical properties, underscoring the necessity to comprehend fiber orientation in injection molded parts. Among different fiber types, glass and carbon fibers are most common in the composite industries. The recent trend of hybrid composites comprising both glass fiber (GF) and carbon fiber (CF) is also gaining importance in the automotive industry. Hybrid fiber options allow designers to optimize the balance between glass and carbon fibers by leveraging the high durability and low cost of GF while the strength and lightweight properties of CF. Consequently, comparing the fiber oriented distribution (FOD) of injection molded composites containing GF, CF, and a hybrid of GF/CF is critical to investigating the local mechanical properties of intricate structures for high-end applications. In Phase I of this project, FOD in injection-molded panels with respect to distance from the gate was analyzed using X-ray computed tomography (X-CT) for GF, CF, and hybrid CF/GF (CGF) reinforced nylon 66. To understand the reason behind the FOD with different fiber types, computational fluid dynamics (CFD) and rheology were performed. Samples were extracted at three locations: near the gate, center, and opposite end. Thickness of the layers of typical skin-shell-core type FOD varies with fiber type and location. GF achieved flow direction alignment (in shell) earlier than viscous CF and CGF near the gate, whereas CF showed the highest flow-direction alignment at the center due to shear induced orientation. At the opposite end, GF experienced more backflow than others indicating faster mold filling owing to its lower viscosity. Hybrid CGF exhibited GF-dominated center and CF-dominated end region. The numerical model used to obtain FOD and rheological predictions for the CF and GF composites served to corroborate the trends observed in the experimental trials. The FOD responses across fiber types and location were reflected in their longitudinal and transverse properties. Only GF showed higher longitudinal modulus over transverse modulus near the gate attributed to rapid alignment, whereas CF and CGF exhibited opposite trend. However, fountain flow enhanced the longitudinal modulus over transverse modulus with the distance for all, particularly for CF. This study offers insights into mold filling behavior of different fibers which are critical in optimizing injection molding conditions for tailored final properties.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Implementation of a High-Fidelity Interface Resolving Method in Nek5000

The development and utilization of computational fluid dynamics (CFD) models for large, high-temperature electric melters in the Waste Treatment and Immobilization Plant (WTP) in eastern Washington State have proven to be valuable for various purposes. These models allow for a better understanding of the physio-chemical processes occurring within the melter vessels and can contribute to improving operational efficiency, throughput, and addressing operational issues related to vitrification. The CFD models employed for these melter vessels incorporate multiphase fluid flow and heat transfer simulations in different regions, including the plenum, cold cap, and molten glass regions. As the tank waste and glass formers are introduced into the melter, a reacting batch layer known as the cold cap forms on top of the molten glass. To enhance the melt rate, forced convection bubblers located at the bottom of the melters generate convection currents that help homogenize the molten glass and provide heat to the cold cap. As the bubbles rise through the highly viscous glass, they adopt a spherical-cap shape [1]. Meanwhile, the conversion of the batch to glass generates significant amounts of gases (such as water vapor, carbon dioxide, sulfur dioxide, and NOx) due to thermal decomposition [2]. These gases become trapped between the cold cap and molten glass, forming a foam layer [3]. For modeling multiphase flow in CFD and heat transfer simulations of waste glass melters at different scales, efforts are underway to augment the capabilities of the Nek5000 [4] and NekRS [5] open-source codes [6]. Nek5000/NekRS is a scalable and efficient spectral element code that has been successfully applied to a wide range of fluid dynamics problems. By leveraging the Nek5000/NekRS software, it becomes possible to model the melter systems more affordably and with lower computational requirements compared to currently utilized commercial CFD software. The specific objectives of this ongoing effort include: 1. Implementation of a level set method in Nek5000/NekRS: The level set method is a numerical technique commonly employed in CFD simulations to track and represent the interface between different phases or materials accurately. By incorporating this method into Nek5000/NekRS, the ability to simulate multiphase flows in waste glass melters at a high level of fidelity can be achieved. 2. Demonstration of capability for air bubbling through molten glass: As part of the development process, a specific case of air bubbling through molten glass will be simulated using the augmented Nek5000/NekRS code. This demonstration aims to showcase the ability of the software to accurately capture and analyze the complex phenomena involved in the multiphase flow within waste glass melters. By achieving these objectives, the improved Nek5000/NekRS code will offer a powerful computational tool for simulating and analyzing waste glass melter systems, enabling better understanding, optimization, and troubleshooting of these vitrification processes. The ability to accurately model and simulate multiphase flows has broad relevance across many industries and scientific domains, and the improved functionality can contribute to advancements in various fields beyond waste glass melter simulations.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗