Search NASA⌕ Search

SEARCH · Search NASA

Results for “low viscosity”

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 163 records · Page 9
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

A dynamic solvent chamber propagation estimation framework using RNN for warm solvent injection in heterogeneous reservoirs

Warm solvent injection (WSI), injecting low-temperature solvent into formations to reduce the viscosity of heavy oil, is a clean technology for heavy oil production through reducing greenhouse gas emissions and water usage. The success of WSI operation depends on the uniform development and propagation of solvent chambers in reservoirs. However, reservoir heterogeneity stemming from shale barriers plays a detrimental role in the conformance of solvent chamber development and oil production rate. In this work, we developed a novel recurrent neural network (RNN)-based framework with the capability of efficiently tracking and estimating the solvent chamber positions in heterogeneous reservoirs based on only production time-series data. The developed estimation model utilizes the “sequence-to-sequence" mapping methodology to correlate observed production time-series sequence and solvent chamber edge sequence via a long short-term memory (LSTM) algorithm. The trained RNN models exhibit high accuracy, evidenced by the predicted dynamic solvent chamber locations match the corresponding true locations from numerical simulation, with a high coefficient of determination (R 2 ) and a low mean squared error. Specifically, the achieved R 2 values exceed 0.98 on both the training and testing data. The developed RNN-based workflow was tested via several cases from both regularly- and irregularly-shaped shale barriers, and the results were promising. The predicted solvent chambers showed strong agreement with those obtained from numerical simulations. The major benefits of this workflow include reducing computational time and saving overall monitoring and tracking costs for conventional techniques. In conclusion, the present work would provide a good demonstration of the capability of practical integration of machine learning methods in solving engineering problems.

58 GEOSCIENCES↗
DOE OSTI2023

Viscosity anomaly of a metallic glass-forming liquid under high pressure

Viscosity, as a critical property closely associated with the glass-forming ability of a liquid, has been extensively studied with varying temperatures. However, its pressure dependence has not been well explored yet due to experimental difficulties. Here, we measured the viscosity of a metallic glass-forming liquid, Zr 46 Cu 37.6 Ag 8.4 Al 8 , at pressures up to 6.1 GPa above the melting points by falling-sphere viscometry with ultrafast synchrotron x-ray imaging. Overall, the viscosity increases with pressure, while surprisingly, there is an abrupt drop between 3.2 GPa and 3.7 GPa, indicating the possible existence of a pressure-induced liquid-to-liquid transition. Pressure could change the short- and medium-range orders in the multicomponent glass-forming liquid, as suggested by the different crystalline outcomes after cooling to room temperature at high- and low-pressure ranges. In conclusion, our work extended the viscosity investigation of metallic glass-forming liquids to the high-pressure regime, which will expand our understanding of liquid-liquid transitions and metallic glass formation.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Non-linear flow modes of identified particles in Pb-Pb collisions at $ \sqrt{s_{\mathrm{NN}}} $ = 5.02 TeV

The pT-differential non-linear flow modes, v 4,22 , v 5,32 , v 6,33 and v 6,222 for π ± , K ± , K$^0_S$ , $p + \bar{p}, Λ + \bar{Λ}$ and ϕ-meson have been measured for the first time at $\sqrt{s_{NN}}$ = 5.02 TeV in Pb-Pb collisions with the ALICE detector at the Large Hadron Collider. The results were obtained with a multi-particle technique, correlating the identified hadrons with reference charged particles from a different pseudorapidity region. These non-linear observables probe the contribution from the second and third order initial spatial anisotropy coefficients to higher flow harmonics. All the characteristic features observed in previous p T -differential anisotropic flow measurements for various particle species are also present in the non-linear flow modes, i.e. increase of magnitude with increasing centrality percentile, mass ordering at low p T and particle type grouping in the intermediate p T range. Hydrodynamical calculations (iEBE-VISHNU) that use different initial conditions and values of shear and bulk viscosity to entropy density ratios are confronted with the data at low transverse momenta. These calculations exhibit a better agreement with the anisotropic flow coefficients than the non-linear flow modes. These observations indicate that non-linear flow modes can provide additional discriminatory power in the study of initial conditions as well as new stringent constraints to hydrodynamical calculations.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2024

Self-organization in collisionless, high- β turbulence

The magnetohydrodynamic (MHD) equations, as a collisional fluid model that remains in local thermodynamic equilibrium (LTE), have long been used to describe turbulence in myriad space and astrophysical plasmas. Yet, the vast majority of these plasmas, from the solar wind to the intracluster medium (ICM) of galaxy clusters, are only weakly collisional at best, meaning that significant deviations from LTE are not only possible but common. Recent studies have demonstrated that the kinetic physics inherent to this weakly collisional regime can fundamentally transform the evolution of such plasmas across a wide range of scales. Here, we explore the consequences of pressure anisotropy and Larmor-scale instabilities for collisionless, β $\gg$ 1, turbulence, focusing on the role of a self-organizational effect known as ‘magneto-immutability’. We describe this self-organization analytically through a high-β, reduced ordering of the Chew–Goldberger–Low-MHD (CGL-MHD) equations, finding that it is a robust inertial-range effect that dynamically suppresses magnetic-field-strength fluctuations, anisotropic-pressure stresses and dissipation due to heat fluxes. As a result, the turbulent cascade of Alfvénic fluctuations continues below the putative viscous scale to form a robust, nearly conservative, MHD-like inertial range. These findings are confirmed numerically via Landau-fluid CGL-MHD turbulence simulations that employ a collisional closure to mimic the effects of microinstabilities. We find that microinstabilities occupy a small (~5%) volume-filling fraction of the plasma, even when the pressure anisotropy is driven strongly towards its instability thresholds. We discuss these results in the context of recent predictions for ion-vs-electron heating in low-luminosity accretion flows and observations implying suppressed viscosity in ICM turbulence.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Advanced Turbulence Models for Large-Scale Atmospheric Boundary Layer Flows

We present high-fidelity large-eddy-simulation (LES) modeling approaches for the turbulent atmospheric boundary layer (ABL) flows. Wind energy is a prime example of an application driven by ABL. Generation of electrical energy from farms of wind turbines at night in the stable ABL is a particularly interesting situation. In this report, we consider the well-known GEWEX (Global Energy and Water Cycle Experiment) Atmospheric Boundary Layer Study (GABLS) stably stratified benchmark LES case. We use a high-order spectral element code Nek5000/RS, which is supported under the DOE's Exascale Computing Project (ECP) Center for Efficient Exascale Discretizations (CEED) project, targeting application simulations on various acceleration-device based exascale computing platforms. In our earlier ANL report, we demonstrated our newly developed subgrid-scale (SGS) models based on high-pass filter (HPF), mean-field eddy viscosity (MFEV), and Smagorinsky (SMG) with no-slip and traction boundary conditions, provided with low-order statistics, convergence and turbulent structure analysis. In this report, we extend the range of our SGS modeling approaches in the context of the mean-field eddy viscosity (MFEV), to include the solution of an SGS turbulent kinetic energy equation (TKE). We demonstrate the model fidelity of Nek5000/RS in comparison to that of AMR-Wind, a block-structured second-order finite-volume code with adaptive-mesh-refinement capabilities, with which we studied scaling performance for both codes in comparison on DOE's leadership computing platforms.

17 WIND ENERGY↗
DOE OSTI2020

Enhancing Heavy-Oil-Recovery Efficiency by Combining Low-Salinity-Water and Polymer Flooding

Summary Combining low-salinity-water (LSW) and polymer flooding was proposed to unlock the tremendous heavy-oil resources on the Alaska North Slope (ANS). The synergy of LSW and polymer flooding was demonstrated through coreflooding experiments at various conditions. The results indicate that the high-salinity polymer (HSP) (salinity = 27,500 ppm) requires nearly two-thirds more polymer than the low-salinity polymer (LSP) (salinity = 2,500 ppm) to achieve the target viscosity at the condition of this study. Additional oil was recovered from LSW flooding after extensive high-salinity-water (HSW) flooding [3 to 9% of original oil in place (OOIP)]. LSW flooding performed in secondary mode achieved higher recovery than that in tertiary mode. Also, the occurrence of water breakthrough can be delayed in the LSW flooding compared with the HSW flooding. Strikingly, after extensive LSW flooding and HSP flooding, incremental oil recovery (approximately 8% of OOIP) was still achieved by LSP flooding with the same viscosity as the HSP. The pH increase of the effluent during LSW/LSP flooding was significantly greater than that during HSW/HSP flooding, indicating the presence of the low-salinity effect (LSE). The residual-oil-saturation (Sor) reduction induced by the LSE in the area unswept during the LSW flooding (mainly smaller pores) would contribute to the increased oil recovery. LSP flooding performed directly after waterflooding recovered more incremental oil (approximately 10% of OOIP) compared with HSP flooding performed in the same scheme. Apart from the improved sweep efficiency by polymer, the low-salinity-induced Sor reduction also would contribute to the increased oil recovery by the LSP. A nearly 2-year pilot test in the Milne Point Field on the ANS has shown impressive success of the proposed hybrid enhanced-oil-recovery (EOR) process: water-cut reduction (70 to less than 15%), increasing oil rate, and no polymer breakthrough so far. This work has demonstrated the remarkable economical and technical benefits of combining LSW and polymer flooding in enhancing heavy-oil recovery.

Engineering↗
DOE OSTI2022

Error quantification of the Arrhenius blending rule for viscosity of hydrocarbon mixtures

Six hundred and seventy-five measurements of dynamic viscosity and density have been used to assess the prediction error of the Arrhenius blending rule for kinematic viscosity of hydrocarbon mixtures. Major trends within the data show that mixture complexity–binary to hundreds of components—and temperature are more important determinants of prediction error than differences in molecular size or hydrogen saturation between the components of the mixtures. Over the range evaluated, no correlation between prediction error and mole fractions was observed, suggesting the log of viscosity truly is linear in mole fraction, as indicated by the Arrhenius blending rule. Mixture complexity and temperature also impact molar volume and its prediction. However, a linear regression between the two model errors explains less than 20% of the observed variation, indicating that mixture viscosity and/or molar volume are not linear with respect to temperature and/or mixture complexity. Extensive discussion of the intermolecular forces and the geometric arrangement of molecules and vacancies in liquids, which ultimately determines its viscosity, is brought into context with the implicit approximations behind the Arrhenius blending rule. The complexity of this physics is not compatible with a simple algebraic correction to the model. However, sufficient data is now available to determine confidence intervals around the prediction of fuel viscosity based on its component mole fractions and viscosities. At -40°C, when all identified components are pure molecules the modeling error is 13.2% of the predicted (nominal) viscosity times the root mean square of the component mole fractions.

10 SYNTHETIC FUELS↗
DOE OSTI2021

Temperature evolution of the phonon dynamics in the Kitaev spin liquid

In this work, we present a study of the phonon dynamics in the honeycomb Kitaev spin model at finite temperatures. We show that the fractionalized spin excitations of the Kitaev spin liquid, the itinerant Majorana fermions, and static Z 2 fluxes have distinct effects on the phonon dynamics, which makes the phonon dynamics a promising tool for exploring the Kitaev spin liquid candidate materials. In particular, we focus on the signature of the fractionalized excitations in the thermodynamic behavior of the sound attenuation and the phonon Hall viscosity: the former describes the phonon decay into the fractionalized excitations, and the latter is the leading order time-reversal symmetry breaking effect on the acoustic phonon. We find that the angular dependence of the attenuation coefficient and its magnitude is modified by the thermal excitation of the Z 2 fluxes. The strength of this effect strongly depends on the relative magnitude of the sound velocity and the Fermi velocity characterizing the low-energy Majorana fermions. We also show that the Hall viscosity is strongly suppressed by the increase of the density of the Z 2 fluxes at finite temperatures. All these observations reflect the effects of the emergent disorder on the Majorana fermions introduced by the Z 2 fluxes. Our analysis is based on the complementary analytical calculations in the low-temperature zero-flux sector, and numerical calculations in the inhomogeneous flux sectors at intermediate and high temperatures with stratified Monte Carlo (strMC) method.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Improved Polypropylene Thermoformability through Polyethylene Layering

Due to its low cost, stiffness, and recyclability, isotactic polypropylene (iPP) is an excellent candidate for packaging applications. However, iPP is notoriously difficult to thermoform due to its low melt strength. The addition of just 10 thin layers of high-molecular-weight, linear low-density polyethylene (LLDPE) into iPP sheets by coextrusion significantly increased extensional viscosity and reduced sag. Both LLDPE and iPP were metallocene-catalyzed with excellent adhesion as measured in our previous work. Here, we performed a series of hot tensile tests and sheet sag measurements to determine the properties of the iPP sheet and the multilayer sheet between 130 and 180 °C. To evaluate the thermoformability of these multilayer sheets, truncated conical cups were positive vacuum formed at different temperatures and heating times, and the crush strength was measured. Cups that released easily from the mold with good shape retention and a crush strength within 80% of the maximum value were used to define a temperature–time thermoformability window. We estimated the maximum stress that occurred during the thermoforming process to be 5 MPa. Layer thicknesses before and after thermoforming were used to estimate an average strain of 0.78. The thin LLDPE layers decreased the yield stress below 5 MPa. This enabled thermoforming at sheet temperatures as low as 150 °C. The immiscible LLDPE interfaces increased extensional viscosity, which decreased sag in the multilayer sheets compared to iPP. This broadened the thermoforming range to temperatures as high as 180 °C and allowed longer heating times. These highly thermoformable, layered sheets may be recycled as iPP since they contain only 8% of LLDPE.

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

Modeling Volatility-Based Aerosol Phase State Predictions in the Amazon Rainforest

Organic aerosol (OA) is a complex matrix of various constituents—fresh (primary organic aerosols—POA) and aged via oxidation (secondary organic aerosols—SOA), generated from biogenic, anthropogenic, and biomass burning sources. The viscosity of OA can be critical in influencing new particle formation, reactive uptake processes that impact evaporation-growth kinetics, and the lifetime of particles in the atmosphere. We report this work utilizes a well-defined relationship between volatility and viscosity for pure compounds, which we incorporated within the Weather Research and Forecasting Model coupled to chemistry (WRF-Chem) to simulate the phase state and viscosity of bulk OA during the dry-to-wet transition season (September–October) in the Amazon rainforest during 2014. Our simulations indicate spatial and temporal heterogeneity in aerosol phase state often not captured by global-scale models. We show the strong role of water associated with organic aerosol ($w_s$) as the dominant factor that can be used to quantitatively estimate OA viscosity. Analysis of WRF-Chem simulations across the entire atmospheric column indicates a strong inverse log-linear relationship between $w_s$ and OA viscosity with a correlation coefficient approaching 1, in the background and biomass burning-influenced conditions. At high altitudes where relative humidity (RH) and temperatures are low, our simulations indicate that OA exists in a semisolid-/solid-like phase state, consistent with previous studies. OA hygroscopicity is strongly correlated (ca. –0.8) with OA viscosity at RH ca. 30–50%, but this RH range is found mostly at low OA concentrations and the middle troposphere (ca. 6–10 km altitudes) in our simulated domain. OA hygroscopicity is uncorrelated with viscosity at higher-RH (near surface) and lower-RH (upper troposphere) regimes. At the urban site near surface, where day–night differences in RH are significant, RH is found to drive the phase state. At the background forested site near surface, where day–night RH differences are small, biomass burning-influenced OA is semisolid and a significant OA associated with background conditions is liquid-like. Simulations indicate a long tail of OA viscosity frequency distributions extending in the semisolid/solid regimes over background biogenic-influenced conditions due to the role of low-volatility OA components such as monoterpene oxidation products.

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

Sub-Ambient Performance of Potassium Sarcosinate for Direct Air Capture Applications: CO 2 Flux and Viscosity Measurements

Absorption based direct air capture (DAC) technologies have garnered significant interest in recent years due to their scalability, competitive regeneration energy requirements, and low susceptibility to degradation. One of the key advantages of DAC lies in its flexible siting options and the potential to utilize low-value land. However, most of the research in this field has been focused on ambient climate zones (T > 20 °C), overlooking sub-ambient (–30 °C < T < 20 °C) regions, which comprise approximately 70 % of the Earth’s surface. To fully realize the potential of DAC, it is essential to understand how DAC solvents perform in these sub-ambient conditions before any large-scale deployment can be considered. Among DAC solvents of interest, potassium sarcosinate (K-SAR) has emerged as a promising candidate due to its high CO 2 capacity, fast uptake kinetics, compatibility with contactor packing materials, low volatility, good thermal and oxidative stability, and competitive regeneration energy requirements compared to current industry standards. This paper characterizes the CO 2 flux and viscosity of K-SAR at sub-ambient conditions and explores the potential of using ethylene glycol and triethylene glycol as additives to prevent solvent freezing in DAC applications. For 1 M K-SAR, the CO 2 flux ranges between 1.3 × 10 -5 and 8.0 × 10 -5 mol m –2 s –1 across a temperature range of –5 °C to 45 °C. Ethylene glycol is shown to effectively suppress the freezing point of K-SAR below –30 °C with volumetric loadings of the additive as low as 0.1. Here, a reaction model was developed to predict the CO 2 flux for 1 M K-SAR at different temperatures, demonstrating good agreement between experimental and theoretical fluxes.

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

NH 4 OH Looping with Membrane CO 2 Absorber and Distributed Stripper for Enhanced Algae Growth

The University of Kentucky Center for Applied Energy (UK CAER) has devised a unique, integrated CO2 capture and utilization technology. CO2 from coal-fired power generation flue gas is first captured at half the operating cost of a typical aqueous CO2 capture system (CCS), distributed in an aqueous stream and then fixed by algae in bioreactors where the algae production is increased by 50% over that with a typical intermittent nutrient feeding system. Lower CCS operating cost is achieved by eliminating the flue gas pretreatment step for cooling and SO2 removal, eliminating steam extraction from the power generation steam cycle for solvent regeneration, and eliminating CO2 compression. Higher algae production is achieved by continuous, just-in-time nutrient feed to the bioreactors directly from a distributed solvent regenerator, which maintains the bioreactor pH for optimum growth. The process starts with a uniquely configured membrane absorber, where the flue gas is indirectly contacted with an ammonium hydroxide (NH4OH) solvent. Dissolved NH3 is attractive for both CO2 capture and as an algae nutrient. For CO2 capture it is inexpensive, has a low regeneration energy, is thermally- and oxidatively-stable and has a viscosity near that of water, which makes is easy to transport. Numerous studies have shown that the scrubbing capacity of NH3 is approximately 0.9-1.2 kg of CO2/kg of NH3, with a CO2 removal efficiency of ~99% and half the solvent regeneration energy than that of 30 wt% MEA[1, 2, 3]. NH3 is attractive as an algae nutrient due to its low cost. The rich NH4OH solvent is pumped to a set of distributed regenerators which are co-located with the algae bioreactors. Solvent pumping, transport and distribution reduces the balance of plant (BOP) cost compared to a typical aqueous CCS related to the flue gas duct and boost fan required to transport the flue gas. The energy required for the distributed solvent regeneration is supplied by solar-thermal panels eliminating the need for steam extraction from the power generation steam cycle. After solvent regeneration, the product stream contains both the CO2 captured from the flue gas and volatized NH3 from the solvent. This product stream is fed directly to the bioreactors, eliminating the need for compression of the CO2 stream. The relative amounts of CO2 and NH3 in the product stream are adjusted and controlled by a controlling the regeneration conditions (pressure and temperature). The continuous feed of the right ratio of nutrients overcomes the typical inhibition of algae growth resulting from frequent pH swings in the bioreactor due to unbalanced (intermittent) feeding systems for CO2 and N. Also, because the regenerators will operate at pressure and be located in close proximity to the bioreactors, there is no worry about pressure drop when sparging the gas into the algae. Sparging produces small bubbles which is beneficial for mass transfer efficiency. One known challenge when using an NH4OH solvent is high NH3 emission. Hydrophobic membranes are used for CO2 capture using an aqueous NH3 solution[4, 5] without the direct contact between flue gas and aqueous solution. Additionally, UK CAER CO2 capture and utilization process manages NH3 slip in three extra measures. First, NH3 slip is minimized by working with minimal species partial pressure, which is proportional to the concentration in the liquid. Hence, lowering the capture solvent concentration will lower the NH3 partial pressure. Second, UK CAER’s previous work has demonstrated that the addition of Zn2+ into NH3 solutions to chelate the NH3 can reduce NH3 volatility. Third, the configuration of the membrane CO2 absorber utilizes condensed water from the flue gas to continually wash the gas-side of the membrane to reduce fouling and recapture NH3 slip. Additional details about the UK CAER unique, integrated CO2 capture and utilization technology will be presented along with technology development plans. Diao, N., Q. Li, and Z. Fang. 2004. Heat transfer in ground heat exchangers with groundwater advection. International Journal of Thermal Sciences. 43: 1203-1211, He, Q., M. Chen, L. Meng, K. Liu, and W. Pan. 2004. Study on Carbon Dioxide Removal from Flue Gas by Absorption of Aqueous Ammonia. Western Kentucky University. Yeh, A.C., and H. Bai. 1999. Comparison of ammonia and monoethanolamine solvents to reduce CO2 greenhouse gas emissions. The Science of the Total Environment. 228: 121-133, Villeneuve, K., D. Roizard, J.C. Remigy, M. Iacono, and S. Rode. 2018. CO2 capture by aqueous ammonia with hollow fiber membrane contactors: Gas phase reactions and performance stability. Separation and Purification Technology, 199: 189-197, Toro Molina, C., and C. Bouallou. 2016. Carbon dioxide absorption by ammonia intensified with membrane contactors. Clean Techn Environ Policy 18, 2133–2146 (2016)

20 FOSSIL-FUELED POWER PLANTS↗
DOE OSTI2024

Molecular dynamics simulations of uranyl and plutonyl cations in a task-specific ionic liquid

Ionic liquids (ILs) are a unique class of solvents with potential applications in advanced separation technologies relevant to the nuclear industry. ILs are salts with low melting points and a wide range of tunable physical properties, such as viscosity, hydrophobiciy, conductivity, and liquidus range. ILs have negligible vapor pressure, are often non-flammable, and can have high thermal stability and a wide electrochemical window, making them attractive for use in separations processes relevant to the nuclear industry. Metal salts generally have a low solubility in ILs; however, by incorporating new functional groups into the IL cation or anion that promote complexation with the metal, the solubility can be greatly increased. One such task-specific ionic liquid (TSIL) is 1-carboxy-N, N, N-trimethylglycine bis(trifluoromethylsulfonyl)imide ([Hbet][Tf 2 N]). Water, which is detrimental for electrochemical separations, is a common impurity in ILs and can coordinate with actinyl cations, particularly in ILs containing only weakly coordinating components. Understanding the behavior of actinides in TSIL/water mixtures on a molecular level is vital for designing improved separations processes. Classical molecular dynamics simulations of uranyl(VI) and plutonyl(VI) in 1-ethyl-3-methylimidazolium bis(trifluoromethylsulfonyl)imide ([EMIM][Tf 2 N]) with deprotonated Hbet (betaine) and water have been performed to understand the coordination and dynamics of the actinyl cations. We find that betaine is a much stronger ligand than water and prefers to coordinate the metal in a bidentate manner. Potential of mean force simulations yield a relative free energy for betaine coordination of approximately -120 to -90 kJ/mol in mixtures with water. As the amount of betaine coordinated to the actinide increases, the diffusion coefficient of the actinyl cation decreases. Moreover, the betaine ligand is able to bridge between two metal centers, resulting in dimeric complexes with actinide–actinide distances of ~5 Å. Potential of mean force simulations show that these structures are stable, with relative free energies of up to -40 kJ/mol. The crystal structure for [(UO 2 ) 2 (bet) 6 (H 2 O) 2 ][Tf 2 N] 4 shows that the betaine bridges between two uranium atoms to form dimeric complexes similar to those found in our simulations.

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

Droplet bioprinting of acellular and cell-laden structures at high-resolutions

Advances in digital light projection(DLP) based (bio) printers have made printing of intricate structures at high resolution possible using a wide range of photosensitive bioinks. A typical setup of a DLP bioprinter includes a vat or reservoir filled with liquid bioink, which presents challenges in terms of cost associated with bioink synthesis, high waste, and gravity-induced cell settling, contaminations, or variation in bioink viscosity during the printing process. Here, we report a vat-free, low-volume, waste-free droplet bioprinting method capable of rapidly printing 3D soft structures at high resolution using model bioinks and model cells. A multiphase many-body dissipative particle dynamics model was developed to simulate the dynamic process of droplet-based DLP printing and elucidate the roles of surface wettability and bioink viscosity. Process variables such as light intensity, photo-initiator concentration, and bioink formulations were optimized to print 3D soft structures (∼0.4–3 kPa) with a typical layer thickness of 50 µm, an XY resolution of 38 ± 1.5 μm and Z resolution of 237 ± 5.4 µm. To demonstrate its versatility, droplet bioprinting was used to print a range of acellular 3D structures such as a lattice cube, a Mayan pyramid, a heart-shaped structure, and a microfluidic chip with endothelialized channels. Droplet bioprinting, performed using model C3H/10T1/2 cells, exhibited high viability (90%) and cell spreading. Additionally, microfluidic devices with internal channel networks lined with endothelial cells showed robust monolayer formation while osteoblast-laden constructs showed mineral deposition upon osteogenic induction. Overall, droplet bioprinting could be a low-cost, no-waste, easy-to-use, method to make customized bioprinted constructs for a range of biomedical applications.

DLP↗
DOE OSTI2020

Examination of Solvent- Room Temperature Ionic Liquid Interactions using 2D NMR and Computational Chemistry

Rare earth metals (REMs) are the key functioning material found in microelectronics, advanced batteries, and high-efficiency electrical motors. Current methods for producing REMs require high temperatures (i.e. energy intensive) and produce an abundance of hazardous fluoride waste. Room temperature ionic liquids (RTILs) are an ideal alternative electrolyte because they allow for a low-temperature, energy-efficient electrochemical process. RTIL properties such as viscosity and conductivity determine how effective they can be in this process. Since viscosity is directly related to both intra and intermolecular forces, this study sets out to investigate two RTILs-Trihexyltetradecylphosphonium dicyanamide (Cyphos 105) and Trihexyltetradecylphosphonium bistrifimide (Cyphos 109) and their interactions. Results from 2-dimensional nuclear magnetic resonance (2D NMR) spectroscopy show that there is an interaction between the tetradecyl tail and the phosphorus atom of the Cyphos RTIL. In order to establish whether this is an intra or an inter-molecular interaction, DFT-based computational calculations were performed to determine the optimized geometry. Calculations were done both in the gas phase and with propylene carbonate as an augmenting solvent. Preliminary results show that there is minimal intramolecular interaction between the tetradecyl tail and the phosphorus atom of the RTILs, and that the interaction seen in the 2D NMR experiment is likely of an intermolecular nature instead.

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

Comparison of temperature adaptive calibration methods for laser induced fluorescence based fuel-in-oil instrument

Fuel dilution of engine oil (or fuel-in-oil, FiO) is an important issue as multiple and late-cycle fuel injection, integral to many combustion efficiency and emissions improvements (e.g., downsized boosted gasoline engines and catalyst thermal management) increases FiO rate. In addition to causing general wear and corrosion in engine due to decreased oil viscosity and pH buffering, FiO is also believed to cause destructive low speed pre-ignition (or super knock) in boosted SI engines. To understand the effects of engine operating conditions on the FiO rate, an optical diagnostic capable of measuring transient FiO on minute timescales was recently developed and demonstrated on a modified GM Ecotech engine system (Neupane et al., Applied Spectroscopy 2021). The measurement is based on adding a dye to the fuel and monitoring for its presence in oil via laser-induced fluorescence (LIF). The measured LIF signal is related to FiO concentration via pre-determined calibration factors using a multivariate classical least square (CLS) method.Since fluorescence quantum yield is a function of temperature, measured LIF intensity not only depends on FiO concentration but also oil temperature. To expand the applicability of the FiO diagnostic to transient oil-temperature conditions (e.g., cold start in practical engines), this study develops a method to account for oil-temperature variations. The effect of oil temperature (20°C - 95°C) on the LIF spectra of eight FiO samples ranging from ~0.8-15% was investigated. LIF intensities of the FiO samples decreased linearly with increasing temperature; the reductions being more significant at dye peaks. We develop a new calibration model (T-adaptive CLS) incorporating the temperature (T) effects on LIF intensity that enables simultaneous calculation of FiO and oil temperature. The improved FiO diagnostic with T-adaptive calibration is more robust, and applicable to varying oil-temperature conditions. For example, when strategies such as multiple/late fuel injections are applied to overcome cold-start instability due to use of low vaporization bio-based fuels such as ethanol, FiO rate is expected to be very high; the improved T-adaptive FiO diagnostic is hence relevant for engine- and fuel-system calibration and optimization. The diagnostic could also provide validation data for flow-field and spray interaction CFD models, further broadening the diagnostic’s utility for advancing engine technology and efficiency.

Neupane, Sneha↗
DOE OSTI2020

The Sensitivity of Joint Inversions of Seismic and Geodynamic Data to Mantle Viscosity

Seismic tomography has revealed the existence of large-scale velocity heterogeneities in the mantle. The interpretation of seismic velocity anomalies in terms of temperature and chemical composition is nonunique. We use geodynamic observations including gravity, plate motions, dynamic topography, and excess ellipticity of the core-mantle boundary combined with seismic observations to investigate the thermo-chemical structure of the mantle through joint inversions. An outstanding issue, however, is the physical connection between mantle density anomalies and the surface geodynamic observations, which requires knowledge of the mantle viscosity structure. Here we perform joint inversions assuming different viscosity profiles and examine the dependence of the results on the viscosity. We first assume that mantle heterogeneity is due to thermal variations, which places a constraint on the relation between seismic velocity and density, and we subsequently relax the constraint to allow for potential nonthermal effects. In all of our joint inversions, a nonthermal origin of density anomalies is required to explain the geodynamic data, though the amount varies with the assumed viscosity structure. A common observation is a high-density chemical signal in the center of the large low-shear-velocity provinces at the base of the mantle resulting in a near neutral or slightly dense overall buoyancy there. Using the derived density models and their corresponding viscosity profiles, we also calculate instantaneous mantle flow fields. The predicted flow fields derived from joint inversions are generally similar but are quite different from flow fields using density models derived from a posteriori scaling of pure seismic tomography models.

58 GEOSCIENCES↗