Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “complex 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 469 records · Page 26

Computer simulation of reaction and transport of core-scale serpentinization of Fe-bearing olivine and evolution of geological hydrogen

We present reactive flow and transport simulations of core-scale serpentinization of forsteritic, Fe-bearing olivine under hydrothermal conditions. The model captures fluid flow, evolving porosity–permeability, and redox-controlled H 2 production in a cylindrical core over 55 days at 245 °C and 37 bar. Key processes include olivine dissolution, precipitation of serpentine and magnetite, and oxidation of Fe2+ to Fe 3+ . Here, results show strong coupling among flow velocity, alteration front propagation, and spatially heterogeneous H 2 generation. Elevated H2 concentrations align with inlet-localized magnetite precipitation, consistent with experiments. Secondary mineral formation reduces porosity and permeability, altering transport pathways. Reactive flow at the fluid–solid interface governs H 2 generation rates and distribution. Dissolved SiO2 promotes serpentine and talc formation by suppressing brucite, while bicarbonate extends reaction duration and moderates surface complexation effects.

Geologic hydrogen↗

Wyoming CarbonSAFE: Advancing Commercialization of Low-Carbon Energy Technologies in Fossil-Rich Wyoming

Wyoming’s Powder River Basin (PRB) is the most prolific coal producer and exporter in the United States and the State of Wyoming’s largest oil-producing basin. In addition to being a leading energy producer, the PRB is the site of research programs whose aim is to develop and integrate low-carbon technologies into existing fossil fuel energy industries. Much of the PRB’s low-carbon research is focused around Dry Fork Station (DFS), which is the newest coal-fired power station in the western US. Currently, DFS hosts five carbon capture projects, including pilot-scale capture projects that utilize different capture technologies and a full front-end engineering and design (FEED) study, a commercial CO 2 pipeline for nearby CO 2 -enhanced oil recovery industry, and is also co-located with the Wyoming Integrated Test Center, which is a host facility for unconventional carbon utilization research (i.e. flue gas carbon-to-products innovations). DFS is also the host site for the Wyoming CarbonSAFE project. Wyoming CarbonSAFE, funded by the Department of Energy/National Energy Technology Laboratory, is a multi-phased program (currently in Phase II) whose core objective is developing and validating storage sites within a complex capable of storing 50 million metric tonnes of CO 2 over a 30 year project period using carbon capture, utilization and storage (CCUS). This paper will provide an overview of the project to-date, showing that CarbonSAFE goals are achievable with respect to geologic, environmental, regulatory, CO 2 source and economic conditions, and why the State of Wyoming offers one of the more favorable environments to advance the commercialization studies. In 2019, the project team designed and completed a ~10,000 foot stratigraphic test well just south of DFS. From this well, the team collected over 600 feet of core, fluid samples from target injection zones, and a full petrophysical log suite. In addition, legacy 2D seismic lines were acquired and a 3D seismic survey was acquired in the fall of 2020. The objective of these field activities has been to identify and characterize target storage reservoirs and associated caprock. This complex has several reservoirs that could meet commercial injection goals, and over 4000 feet of associated caprock. Simulations of site performance suggest that stacked injection provides the most effective storage strategy, and would necessitate several sites within the greater complex to meet final project injection goals. Other work within this project includes the development of a business-case around proximal fields with CO 2 -enhanced oil recovery potential and existing tax credits, assessment of regulatory conditions, including pore space ownership and Class VI injection well permitting requirements, implementation of a robust public outreach Electronic copy available at: https://ssrn.com/abstract=3821220 GHGT-15 McLaughlin 2 program, and surface site characterization activities that have focused on environmental factors. In the PRB and at DFS, the State of Wyoming and other entities have made carbon management a priority of its future energy industry by providing a regulatory and business framework that is favorable to advancing these technologies. These endeavors will become more realistic with the successful implementation of Wyoming CarbonSAFE, and its ability to secure and validate commercial-scale CCUS at the center of Wyoming’s premier low-carbon research efforts.

McLaughlin, Fred↗

Mechanistic Study of Controlled Zinc Electrodeposition Behaviors Facilitated by Nanoscale Electrolyte Additives at the Electrode Interface

Nanoparticle organic hybrid materials (NOHMs) are liquid-like materials composed of an inorganic core to which a polymeric canopy is ionically tethered. NOHMs have unique properties including negligible vapor pressure, high oxidative thermal stability, and the ability to bind to reactive species of interest due to the tunability of their polymeric canopy. This makes them promising multifunctional materials for a wide range of energy and environmental technologies, including electrolyte additives for electrochemical energy storage (e.g., flow batteries) and the electrochemical conversion of CO 2 to chemicals and fuels. Due to their unique transport behaviors in fluid systems, an understanding of the near-electrode surface behavior of NOHMs in electrolyte solutions and their effect on electrochemical reactions is still lacking. In this work, the complexation of zinc (Zn) by NOHMs with an ionically tethered polyetheramine canopy (HPE) (NOHM-I-HPE) was studied using attenuated total reflectance Fourier transform infrared and Carbon-13 nuclear magnetic resonance spectroscopy. Additionally, various electrochemical techniques were employed to discern the role of NOHM-I-HPE during zinc electrodeposition, and the results were compared to those of the electrochemical system containing untethered HPE polymers. Our findings confirmed that NOHM-I-HPE and HPE reversibly complex zinc in the aqueous electrolyte. NOHM-I-HPE and HPE were found to block some of the electrode active sites, reducing the overall current density during electrodeposition, while facilitating the formation of smooth zinc deposits, as revealed by surface imaging and diffraction techniques. Observed variations in the current density responses and the degree of passivation created by the NOHM-I-HPE and HPE adsorbed on the electrode surface revealed that their different packing behaviors at the electrode–electrolyte interface influence the zinc deposition mechanism. The presence of the nanoparticle and ordering offered by the NOHMs as well as the structured conformation of the polymeric canopy allowed the formation of void spaces and free volumes for enhanced transport behaviors. Furthermore, these findings provided insights into how structured electrolyte additives such as NOHMs can allow for advancements in electrolyte design for controlled deposition of metal species from energy-dense electrolytes or for other electrochemical reactions.

25 ENERGY STORAGE↗

Forecasting high-dimensional spatio-temporal systems from sparse measurements

This paper introduces a new neural network architecture designed to forecast high-dimensional spatio-temporal data using only sparse measurements. The architecture uses a two-stage end-to-end framework that combines neural ordinary differential equations (NODEs) with vision transformers. Initially, our approach models the underlying dynamics of complex systems within a low-dimensional space; and then it reconstructs the corresponding high-dimensional spatial fields. Many traditional methods involve decoding high-dimensional spatial fields before modeling the dynamics, while some other methods use an encoder to transition from high-dimensional observations to a latent space for dynamic modeling. In contrast, our approach directly uses sparse measurements to model the dynamics, bypassing the need for an encoder. This direct approach simplifies the modeling process, reduces computational complexity, and enhances the efficiency and scalability of the method for large datasets. We demonstrate the effectiveness of our framework through applications to various spatio-temporal systems, including fluid flows and global weather patterns. Although sparse measurements have limitations, our experiments reveal that they are sufficient to forecast system dynamics accurately over long time horizons. Our results also indicate that the performance of our proposed method remains robust across different sensor placement strategies, with further improvements as the number of sensors increases. This robustness underscores the flexibility of our architecture, particularly in real-world scenarios where sensor data is often sparse and unevenly distributed.

97 MATHEMATICS AND COMPUTING↗

Active Seismoelectric Monitoring of CO 2 Injection

Conference presentation at Carbon Capture, Utilization, and Storage (CCUS) Conference 2024, Houston, Texas, March 11–13, 2024. We investigated the use of active seismoelectric monitoring (ASE) to determine its effectiveness in detecting CO 2 injected deep underground and as an additional monitoring option for carbon storage projects. The EERC performed an ASE survey in February of 2023 at Pennel oil field near Baker, Montana, as part of the U.S. Department of Energy-sponsored Williston Basin CO 2 Field Laboratory in conjunction with operating partner Denbury Inc. The field site is part of the Cedar Creek Anticline geologic structure, with CO 2 planned for injection into a stacked storage complex (SSC) in the Interlake Formation containing both a residual oil zone and a conventional reservoir. The objective was to test the ASE method’s ability to characterize baseline fluid distribution in both the near surface and the SSC prior to CO 2 injection.

04 OIL SHALES AND TAR SANDS↗

MFIX-Exa: Performance prediction of multiphase energy conversion devices

MFIX-Exa targets the high-fidelity CFD-DEM model in which particles are unresolved by the fluid grid (typically using a fluid mesh approximately twice the particle diameter) but discrete particle dynamics including collisions are resolved with a simple linear-spring dashpot model. The fluid field is solved with a modern low Mach number formation of a cell-centered, nodal-pressure approximate projection method using a Godunov scheme . Physics capabilities include: complex geometries via embedded boundaries (EBs), open and closed system ideal gas equations of state, species transport and heterogeneous chemi- cal reactions. MFIX-Exa is built on the AMReX software framework, the ECP Block-Structured Adaptive Mesh Refinement (AMR) Co-Design Center, which allows the code to be portable and performant.Further integration into the ECP ecosystem includes HYPRE linear solvers and Ascent for in situ visualization. Current project efforts focus on scaling up realistic simulations to the KPP of the challenge problem: NETs 50kW chemical looping reactor, shakedown and scaling on JSLE and OLCF TDS systems and reaching for any remaining performance improvements.

Musser, Jordan↗

A Membrane-Associated Light-Harvesting Model is Enabled by Functionalized Assemblies of Gene-Doubled TMV Proteins

Photosynthetic light harvesting requires efficient energy transfer within dynamic networks of light-harvesting complexes embedded within phospholipid membranes. Artificial light-harvesting models are valuable tools for understanding the structural features underpinning energy absorption and transfer within chromophore arrays. Here, a method for attaching a protein-based light-harvesting model to a planar, fluid supported lipid bilayer (SLB) is developed. The protein model consists of the tobacco mosaic viral capsid proteins that are gene-doubled to create a tandem dimer (dTMV). Assemblies of dTMV break the facial symmetry of the double disk to allow for differentiation between the disk faces. A single reactive lysine residue is incorporated into the dTMV assemblies for the site-selective attachment of chromophores for light absorption. Further, on the opposing dTMV face, a cysteine residue is incorporated for the bioconjugation of a peptide containing a polyhistidine tag for association with SLBs. The dual-modified dTMV complexes show significant association with SLBs and exhibit mobility on the bilayer. The techniques used herein offer a new method for protein-surface attachment and provide a platform for evaluating excited state energy transfer events in a dynamic, fully synthetic artificial light-harvesting system.

59 BASIC BIOLOGICAL SCIENCES↗

Unitary Qubit Lattice Algorithms for Plasma Physics

This final technical report summarizes research conducted under DOE Award DE-SC0021653 to develop unitary Quantum Lattice Algorithms for modeling electromagnetic wave propagation and scattering in complex media, including plasmas. The project developed and validated quantum-inspired formulations of Maxwell's equations that preserve unitary evolution and can be evaluated on classical high-performance computing systems while providing a foundation for future quantum-computing implementations. Major accomplishments include the development of two- and three-dimensional algorithms for electromagnetic scattering; scalable, distributed-memory implementations demonstrated on the Perlmutter supercomputer; formulations for nonlinear lossless fluid dynamics and cold, lossless, inhomogeneous magnetized plasmas; and an explicit quantum algorithm for a time-discretized Lorenz model. Simulations reproduced a range of characteristic wave phenomena, including transient effects that are not readily apparent in conventional frequency-domain studies, demonstrating the effectiveness of the proposed approach for modeling complex electromagnetic and plasma systems. The work establishes a unified theoretical and computational framework for quantum and quantum-inspired simulation and provides a foundation for future implementation on fault-tolerant quantum systems.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Droplet formation simulation using mixed finite elements

Droplet formation happens in finite time due to the surface tension force. The linear stability analysis is useful to estimate the size of a droplet but fails to approximate the shape of the droplet. This is due to a highly nonlinear flow description near the point where the first pinch-off happens. A one-dimensional axisymmetric mathematical model was first developed by Eggers and Dupont [“Drop formation in a one-dimensional approximation of the Navier–Stokes equation,” J. Fluid Mech. 262, 205–221 (1994)] using asymptotic analysis. This asymptotic approach to the Navier–Stokes equations leads to a universal scaling explaining the self-similar nature of the solution. Numerical models for the one-dimensional model were developed using the finite difference [Eggers and Dupont, “Drop formation in a one-dimensional approximation of the Navier–Stokes equation,” J. Fluid Mech. 262, 205–221 (1994)] and finite element method [Ambravaneswaran et al., “Drop formation from a capillary tube: Comparison of one-dimensional and two-dimensional analyses and occurrence of satellite drops,” Phys. Fluids 14, 2606–2621 (2002)]. The focus of this study is to provide a robust computational model for one-dimensional axisymmetric droplet formation using the Portable, Extensible Toolkit for Scientific Computation. The code is verified using the Method of Manufactured Solutions and validated using previous experimental studies done by Zhang and Basaran [“An experimental study of dynamics of drop formation,” Phys. Fluids 7, 1184–1203 (1995)]. The present model is used for simulating pendant drops of water, glycerol, and paraffin wax, with an aspiration of extending the application to simulate more complex pinch-off phenomena.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Extended dynamic mode decomposition for model reduction in fluid dynamics simulations

High computational cost and storage/memory requirements of fluid dynamics simulations constrain their usefulness as a predictive tool. Reduced-order models (ROMs) provide a viable solution to this challenge by extracting the key underlying dynamics of a complex system directly from data. We investigate the efficacy and robustness of an extended dynamic mode decomposition (xDMD) algorithm in constructing ROMs of three-dimensional cardiovascular computations. Focusing on the ROMs' accuracy in representation and interpolation, we relate these metrics to the truncation rank of singular value decomposition, which underpins xDMD and other approaches to ROM construction. Our key innovation is to relate the truncation rank to the singular values of the original flow problem. This result establishes a priori guidelines for the xDMD deployment and its likely success as a means of data compression and reconstruction of the system's dynamics from dominant spatiotemporal structures present in the data.

Mechanics↗

A hybrid machine-learning approach for analysis of methane hydrate formation dynamics in porous media with synchrotron CT imaging

Fast multi-phase processes in methane hydrate bearing samples pose a challenge for quantitative micro-computed tomography study and experiment steering due to complex tomographic data analysis involving time-consuming segmentation procedures. This is because of the sample's multi-scale structure, which changes over time, low contrast between solid and fluid materials, and the large amount of data acquired during dynamic processes. Here, a hybrid approach is proposed for the automatic segmentation of tomographic data from time-resolved imaging of methane gas-hydrate formation in sandy granular media, which includes a deep-learning 3D U-Net model. To prepare a training dataset for the 3D U-Net, a technique to automate data labeling based on sample-specific information about the mineral matrix immobility and occasional fluid movement in pores is proposed. Automatic segmentation allowed for studying properties of the hydrate growth in pores, as well as dynamic processes such as incremental flow and redistribution of pore brine. Results of the quantitative analysis showed that for typical gas-hydrate stability parameters (100 bar methane pressure, 7°C temperature) the rate of formation is slow (less than 1% per hour), after which the surface area of contact between brine and gas increases, resulting in faster formation (2.5% per hour). Hydrate growth reaches the saturation point after 11 h of the experiment. Finally, the efficacy of the proposed segmentation scheme in on-the-fly automatic data analysis and experiment steering with zooming to regions of interest is demonstrated.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

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↗

Engine Combustion System Optimization Using Computational Fluid Dynamics and Machine Learning: A Methodological Approach

Gasoline compression ignition (GCI) engines are considered an attractive alternative to traditional spark-ignition and diesel engines. Here, a Machine Learning-Grid Gradient Ascent (ML-GGA) approach was developed to optimize the performance of internal combustion engines. ML offers a pathway to transform complex physical processes that occur in a combustion engine into compact informational processes. The developed ML-GGA model was compared with a recently developed Machine Learning-Genetic Algorithm (ML-GA). Detailed investigations of optimization solver parameters and variable limit extension were performed in the present ML-GGA model to improve the accuracy and robustness of the optimization process. Detailed descriptions of the different procedures, optimization tools, and criteria that must be followed for a successful output are provided here. The developed ML-GGA approach was used to optimize the operating conditions (case 1) and the piston bowl design (case 2) of a heavy-duty diesel engine running on a gasoline fuel with a research octane number (RON) of 80. The ML-GGA approach yielded >2% improvements in the merit function, compared with the optimum obtained from a thorough computational fluid dynamics (CFD) guided system optimization. The predictions from the ML-GGA approach were validated with engine CFD simulations. This study demonstrates the potential of ML-GGA to significantly reduce the time needed for optimization problems, without loss in accuracy compared with traditional approaches.

33 ADVANCED PROPULSION SYSTEMS↗

High‐Speed Embedded Ink Writing of Anatomic‐Size Organ Constructs

Embedded ink writing (EIW) is an emerging 3D printing technique that fabricates complex 3D structures from various biomaterial inks but is limited to a printing speed of ∼10 mm s −1 due to suboptimal rheological properties of particulate-dominated yield-stress fluids when used as liquid baths. In this work, a particle-hydrogel interactive system to design advanced baths with enhanced yield stress and extended thixotropic response time for realizing high-speed EIW is developed. In this system, the interactions between particle additive and three representative polymeric hydrogels enable the resulting nanocomposites to demonstrate different rheological behaviors. Accordingly, the interaction models for the nanocomposites are established, which are subsequently validated by macroscale rheological measurements and advanced microstructure characterization techniques. Filament formation mechanisms in the particle-hydrogel interactive baths are comprehensively investigated at high printing speeds. To demonstrate the effectiveness of the proposed high-speed EIW method, an anatomic-size human kidney construct is successfully printed at 110 mm s −1 , which only takes ∼4 h. This work breaks the printing speed barrier in current EIW and propels the maximum printing speed by at least 10 times, providing an efficient and promising solution for organ reconstruction in the future.

36 MATERIALS SCIENCE↗

Association Between Injection and Microseismicity in Geothermal Fields With Multiple Wells: Data-Driven Modeling of Rotokawa, New Zealand, and Húsmúli, Iceland

Understanding injection-induced microseismicity in geothermal systems can provide insight into reservoir connectedness. However, fault and reservoir complexity are difficult to represent in simple analytical models, which makes it difficult to discern clear relationships from incidental associations. Here, we have used data-driven models to study how fluid injection and microseismicity are related in the Rotokawa (New Zealand) and Hellisheiði (Iceland) geothermal fields. We tested two classes of model: (a) lagged linear regression of seismicity rate as a function of well injection rates; and (b) systematic extraction of injection time series features that are then evaluated for associations with the seismicity. These models allowed us to determine which wells had the greatest correlation with microseismicity and to explain this association in a reservoir context. Finally, exploring different data types and transformations, we were unable to establish a link between rapid changes in injection rate and seismicity spikes, as suggested by some theoretical models.

15 GEOTHERMAL ENERGY↗

Optimization and Evaluation of Energy Savings for Connected and Autonomous Off-Road Vehicles

Off-road vehicles, such as wheel loaders, excavators, and harvesters, are extensively utilized across a wide range of industries, including construction, agriculture, and mining. These machines have become indispensable in supporting the day-to-day operational needs of a nation, playing a critical role in various sectors' infrastructure and productivity. However, despite their utility, off-road vehicles are significant consumers of fossil fuels, resulting in substantial emissions that contribute to environmental degradation. This highlights the pressing need for research and technological advancements aimed at improving their energy efficiency and reducing their carbon footprint. There are, however, two primary challenges that must be addressed to achieve these goals. First, off-road vehicles typically perform both driving and working tasks simultaneously, which introduces a high level of complexity into their overall dynamic systems. Analysis the interactions between these functions is challenging. Second, research into off-road vehicles is inherently interdisciplinary, demanding expertise across several domains such as fluid power systems, vehicle dynamics, control theory, optimization techniques, and real-world implementation. Recognizing these challenges, we proposed the project titled "Optimization and Evaluation of Energy Savings for Connected and Autonomous Off-Road Vehicles" as a comprehensive solution to enhance fuel efficiency while simultaneously improving productivity. This project specifically focuses on autonomous off-road vehicles, with particular attention to wheel loaders, and seeks to develop novel methods to optimize energy consumption without sacrificing operational performance. The project integrates real-time control algorithms, vehicle dynamics modeling, and co-optimization of powertrain system and vehicle system to achieve these goals. Our optimization strategy dynamically co-optimizes critical parameters at both the powertrain and vehicle levels, including vehicle speed, working tool movements, powertrain dynamics, and engine operations in real-time. To streamline this optimization process, we developed a vehicle model that captures the key dynamics while significantly enhancing computational efficiency. This allows the system to intelligently minimize fuel consumption, all while maintaining or even improving productivity through real-time calculations during various off-road operations. To validate the effectiveness of this energy optimization method, we introduced a state-of-the-art Hardware-in-the-Loop (HIL) testbed. This reconfigurable testbed seamlessly integrates the actual engine with virtual models of the wheel loader's subsystems, allowing for accurate emulation of real-world operational loads and environments. By simulating these conditions, the HIL testbed enables us to evaluate the wheel loader’s performance under diverse working scenarios, ensuring the developed solution is applicable in real-world operations. This testbed proved to be instrumental in validating the optimization algorithms and demonstrating the system's practical effectiveness. During the evaluation and testing phase, we employed the HIL testbed to rigorously assess the energy savings and productivity improvements generated by the optimized system. The results were highly encouraging, revealing that the automated wheel loader achieved over 30% fuel savings compared to traditional, human-operated cycles, with comparable or even enhanced levels of productivity. The insights gained from this HIL-based testing provided critical validation of our approach and highlighted the potential for deploying these optimized autonomous technologies in real-world off-road vehicles.

33 ADVANCED PROPULSION SYSTEMS↗

TCCSP Site Characterization and Feasibility Report

The Tulare County Carbon Storage Project (TCCSP) focused on establishing the geological and commercial viability of a 50 million metric ton (MMT) carbon storage complex located in Tulare County, California. The geologic scope of work included collecting new and existing subsurface data, such as two-dimensional (2D) seismic data, advanced geophysical logs, core and fluid samples, and specialized core-analysis datasets, to support the integration and interpretation of findings from earlier tasks. These efforts contributed to the comprehensive feasibility study for commercial carbon capture and storage (CCS) development at the TCCSP site and in the surrounding region. The TCCSP project team evaluated the feasibility of the Lower Monterey Group (Santa Margarita Sand Member or its equivalent), Temblor Formation sand members (Olcese, Jewett, Vedder), and Lodo-Martinez Zone (Domengine Sandstone, Walker Formation, Lodo Formation, and the Martinez Sand) to serve as commercial CCS reservoirs along with potential for the Upper Monterey Group (Reef Ridge, Antelope Members) and the Kreyenhagen Formation to vertically seal underlying reservoirs from shallow above-zones.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

Insights into mixing of non-isothermal multi-polymer melts for complex plastics recycling

Catalytic recycling or upcycling of plastics is often limited not by catalyst performance, but by transport, arising from highly viscous, non-Newtonian polymer melts. In this work, we develop a reactor-scale framework that integrates rheological measurements, constitutive modeling, computational fluid dynamics (CFD), and experiments to quantify mixing, heat transfer, and dispersion in surrogate hydrocarbon melts representing mixed plastics systems. Temperature- and shear rate-dependent viscosity of low-density polyethylene (LDPE) and high-density polyethylene (HDPE) is measured to create two surrogate polymers (PLD and PHD) that capture the dominant shear-thinning flow behavior while neglecting strong elastic effects, enabling tractable simulation of non-isothermal, polymer-melt mixing using a Carreau-Arrhenius generalized Newtonian framework. Three-dimensional CFD simulations are employed to evaluate impeller performance in PLD using mixing time, cavern volume, thermal uniformity, and interfacial area for regimes in which viscoelastic effects are not dominant. We show that magnetic stir bars commonly used in lab-scale studies produce large thermal gradients (~60 °C) and poor mixing, even under idealized power delivery and polymer flow conditions. In contrast, close-clearance anchor impellers achieve near-isothermal operation, reduce mixing times by up to 5×, and provide >90% active circulation volume. We further demonstrate that, at low pseudo-Deborah number (De*), motor power requirements can be predicted directly from shear rate-dependent rheology using the Carreau-Arrhenius framework, enabling rational selection of operating conditions. Extension to surrogate immiscible multi-polymer systems based on PLD and PHD shows that interfacial area is highly sensitive to operating conditions and impeller design, with coaxial anchor-turbine configurations enhancing dispersion by up to 4 × .

Close-clearance impellers↗