Search NASA⌕ Search

SEARCH · Search NASA

Results for “high 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 199 records · Page 11

Development and demonstration of a rolling ball viscometer for molten salts with near-minimum liquidus NaCl–KCl

The thermophysical characterization of molten chloride and fluoride salts is a key interest for the development of advanced molten salt nuclear reactors. Viscosity is one thermophysical property of interest; a high-accuracy understanding of this property is necessary for modeling the thermal energy transfer and mass flow rate for a given molten salt system. Only a few methods have consistently been used to measure the viscosity of molten chloride and fluoride salt systems, and large discrepancies have been observed in the literature between independent studies of nearly identical salt systems. As such, exploration of novel measurement methods or implementations may allow for further validation, improved measurement accuracy, and overall improved understanding about optimal approaches for measuring molten salt viscosity. Here, in this study, a novel implementation of the rolling ball viscometry method was used to measure the viscosity of a 44/56 mol % mixture of NaCl–KCl within a temperature range of 714–838 °C, with the intent to demonstrate the feasibility of this technique for molten salt viscosity measurement cross-validation. The purity of this salt was confirmed through laser-induced breakdown spectroscopy. Measured viscosities were approximately 1.3–1.0 cP within this temperature range, and the overall uncertainty was 10.8%, which was determined through Gaussian error propagation. The results agree reasonably well with measured viscosities of similar NaCl–KCl mixtures in the literature (within 2%–20% depending on the temperature and study), although the experimental error in some of these comparative studies may be significant.

47 OTHER INSTRUMENTATION↗

Renewable low viscosity polyester-polyols for biodegradable thermoplastic polyurethanes

In the transition to renewably sourced, biodegradable polymers, the preparation of low viscosity polyester-polyols has posed a challenge for renewable polyurethane (PU) development. Low viscosity polyols not only reduce the requirement for high process temperatures but also decrease manufacturing time. In our efforts to incorporate increasing ratios of bio-based monomers into renewable PUs, we mixed diacids such as even carbon sebacic acid and odd carbon azelaic acid along with a renewable diol. This provided library of 2000 g/mol molecular weight polyester-polyols, and structures were established by 1 H and 13 C NMR analysis. The prepared polyester-polyols offered lower viscosity and enable lower fabrication temperatures to make TPUs, and their structure and material metrics were evaluated. The formation of TPUs is ascertained from FTIR and NMR analysis. The final TPUs displayed good physical and mechanical properties. Further, these TPUs exhibited T g in the range of -56.5 to -39.7°C, corresponding to TPU soft block structure, and T m between 98.3 and 105.1°C originating from the hard segment. Prepared TPUs exhibit excellent biodegradation under compost environmental conditions. These TPUs showed up to 57% decrease in molecular weight by GPC analysis after 9 weeks of biodegradation, and respirometer analysis displayed up to 97% biodegradation over 120 days.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Impacts of phase state and water content on secondary organic aerosol formation and partitioning (Final Technical Report)

Secondary organic aerosols (SOA) account for a major fraction of particulate matter in the atmosphere, affecting climate, air quality, and public health. SOA formation and evolution are highly complex processes involving both chemical reactions and transport of molecules in air, at particle surfaces, and within particles. The state of the matter, or “phase state,” comprising SOA can vary from a liquid, over an amorphous semi-solid, to a glassy solid, depending on chemical composition, water content, relative humidity (RH), and temperature. The occurrence of glassy and amorphous semi-solid states can pose limitations on the rate of transport of molecules, affecting gas-particle interactions and challenging the treatment of SOA in atmospheric models. The objective of this project was to improve fundamental understanding of the interplay of the phase state of particles and water content on the evolution of SOA formation. We have developed a method to estimate the glass transition temperature (T g ) of organic compounds to predict the phase state and viscosity of SOA. This method has been applied to high resolution mass spectrometry data of various types of SOA, achieving a good agreement with viscosity measurements. In addition, the viscosity estimation method was implemented into a regional model CMAQ to simulate diel and seasonal variations of SOA phase state over the U.S. Applying kinetic multilayer modeling, we analyzed laboratory experiments and ARM field measurements (e.g., HI-SCALE) on kinetic limitations of amine uptake. We have also conducted numerical simulations to estimate equilibration timescales of SOA partitioning. Overall, we have evaluated the impacts of phase state and water content on SOA lifecycle that should contribute to reducing the uncertainty of SOA representation in regional climate and air quality predictions. The project has yielded 17 publications, in which DE-SC0018349 is explicitly acknowledged. Please see the list of publications as below, followed by summary of research activities.

54 ENVIRONMENTAL SCIENCES↗

New scaling law for turbulent boundary layers with high surface mass transfer

A fundamental understanding of the phenomena occurring in the turbulent boundary layer in the presence of surface blowing is limited, and considerable disagreements persist even in describing primary quantities, such as the boundary layer profile. The theories based on the linear boundary layer equations show that the thickness of the sublayer increases in the presence of surface blowing; therefore, the viscous sublayer and law of the wall modify. In this study, direct numerical simulations (DNS) of turbulent boundary layers with uniform surface mass transfers are carried out in order to scale the velocity profile. Emphasis is placed on moderate to high mass transfer rates, which are relevant to the most common hybrid rockets configuration. DNS data are used to establish a functional law of the wall and a law of wake by means of the relation between the wall shear stress and surface mass transfer. Analysis of the mean kinetic energy budget shows that the magnitude of turbulent kinetic energy increases by surface mass transfer, and the production rate extends significantly in the inner layer as the injection rate increases. DNS data of various surface blowing are used to complete the closure of turbulence kinetic energy equation and develop an eddy viscosity model. The predicted turbulent kinetic energy and eddy viscosity agree with DNS data for moderate to high blowing rates.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Carbonic Anhydrase Enhanced UV-Crosslinked PEG-DA/PEO Extruded Hydrogel Flexible Filaments and Durable Grids for CO 2 Capture

In this study, poly (ethylene glycol) diacrylate/poly (ethylene oxide) (PEG-DA/PEO) interpenetrating polymer network hydrogels (IPNH) were extruded into 1D filaments and 2D grids. The suitability of this system for enzyme immobilization and CO 2 capture application was validated. IPNH chemical composition was verified spectroscopically using FTIR. The extruded filament had an average tensile strength of 6.5 MPa and elongation at break of 80%. IPNH filament can be twisted and bent and therefore is suitable for further processing using conventional textile fabrication methods. Initial activity recovery of the entrapped carbonic anhydrase (CA) calculated from esterase activity, showed a decrease with an increase in enzyme dose, while activity retention of high enzyme dose samples was over 87% after 150 days of repeated washing and testing. IPNH 2D grids that were assembled into spiral roll structured packings exhibited increased CO 2 capture efficiency with increasing enzyme dose. Long-term CO 2 capture performance of the CA immobilized IPNH structured packing was tested in a continuous solvent recirculation experiment for 1032 h, where 52% of the initial CO 2 capture performance and 34% of the enzyme contribution were retained. These results demonstrate the feasibility of using rapid UV-crosslinking to form enzyme-immobilized hydrogels by a geometrically-controllable extrusion process that uses analogous linear polymers for both viscosity enhancement and chain entanglement purposes, and achieves high activity retention and performance stability of the immobilized CA. Potential uses for this system extend to 3D printing inks and enzyme immobilization matrices for such diverse applications as biocatalytic reactors and biosensor fabrication.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

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↗

Methodologies for Design, Characterization and Testing of Electrolytes that Enable Extreme Fast Charging of Lithium-ion Cells

Selection, testing and validation of electrolyte candidates for Li-ion cells are discussed, based on a 10-minute target for extreme fast charge (XFC). A combination of modeling and laboratory measurements create a timely and synergistic approach to identifying candidate electrolyte formulations. Multi-solvent systems provide a balanced set of properties, wherein lower molecular-weight solvents offer reduced viscosity, increased species diffusivity, and mitigation of concentration polarization at high charge rates. Carefully selected formulations can exhibit peak conductivity and usable conductivity range of two to three times that of the baseline EC-EMC (3:7, wt.) + LiPF 6 . Candidates are also chosen based on stability and longevity within the cell environment. Lab testing coincides with property predictions from the Advanced Electrolyte Model (AEM) and a macro-scale cell model. Furthermore, cell testing utilized coin and pouch cells having NMC532 or NMC811 cathodes with graphite electrodes. Results indicate combinations of low-molecular weight solvents are key for fast-charge electrolytes as they extend the useful conductivity range to both low and higher salt concentrations, and possess higher self-diffusivities compared to conventional solvents. This reduces impacts from concentration polarization. The choice of electrolyte also influences the tendency for lithium metal deposition at the anode, as showcased by experimental and modeling results herein.

25 ENERGY STORAGE↗

The influence of laser power modulation on melt pool dynamics in laser powder bed fusion

While the majority of laser powder bed fusion (LPBF) metal additive manufacturing uses a continuous wave (CW) laser heat source, some commercial applications of LPBF additive manufacturing instead involve the modulation of the laser power on tens-of-microsecond timescales as an adjustable process variable. This article reports the use of in situ, high speed x-ray and optical imaging to probe melt pool fluid flow, defect formation, and nearby powder motion during LPBF with both modulated and CW laser heat sources. We observe melt pool dynamics unique to modulated laser melting even at very high duty cycles that are related to fluctuations in vapor depression depth, complex pore formation mechanisms, and changes to denudation physics when compared to CW melting. These behaviors are present in Ti–6Al–4V, 316L stainless steel, and AL1100 alloys but vary slightly as a function of material, indicating a substantial dependence on the viscosity and surface tension of the liquid metal. While high duty cycles produce weld tracks of comparable quality to CW melting, lower duty cycles introduce substantial defect concentrations. At intermediate duty cycles, careful control of modulation parameters can repeatably and precisely yield one pore per laser pulse, suggesting a method for intentionally inserting engineered porosity at specific sites during an LPBF build.

3D printing↗

Effect of Thermal Expansion Coefficient, Viscosity and Melting Range in Simulation of PCM Embedded Heat Exchangers With and Without Fins

Phase change material heat exchangers (PCMHX) have an important role in integrating renewable energy systems. PCMHX can offer high storage density at various temperatures of interest and can be used for grid load shifting purposes. Numerical models enable engineers to estimate PCMHX performance for different design parameters and operating conditions. Modeling phase change phenomena is challenging due to the complex time-dependent nature of the process. The accuracy of models is highly sensitive to PCM thermo-physical properties. Thermal expansion coefficient (β), viscosity (μ) and melting temperature range (MR) of a PCM are important properties, especially when natural convection is not negligible. In PCMHX modeling, using less than accurate values for these properties can have significant impact on the simulation outcomes. These properties and discussions thereof are not readily available in the literature. This paper presents a brief review of the literature and a numerical study investigating the model sensitivity to the above-mentioned properties for PCMHX with and without fins. CFD is used to evaluate the charging (melting) phenomena. The study quantifies the impact of uncertainty in these properties on the melting rate and temperature profile. Constant wall temperature was considered as heat source with no heat loss to ambient. Here, the results show that β and μ has significant effect on the melting rate and evolution of the melting front. For a non-finned domain when comparing results for different published values of β and μ, the deviation in melting time can be up to 12.9% and 57.6% respectively. For high wall temperatures, change in melting range did not impact melting time. But when the wall temperature is reduced, up to 9.8% deviation in melting time is observed.

CFD↗

The ultralow viscosity of volatile-rich kimberlite magma: Implications for the water content of primitive kimberlite melts

The eruption of deeply sourced kimberlite magma offers the fastest route to bring deep-seated volatiles back to the Earth’s surface. However, the viscosity of kimberlite magma, a factor governing its migration and eruption dynamics within Earth, remains poorly constrained. We conducted synchrotron in situ falling sphere viscometry experiments to examine kimberlite magma with different volatile contents (0 to 5 wt % H 2 O and 2 to 8 wt % CO 2 ) under high pressure-temperature conditions. The results reveal that the viscosity of volatile-rich kimberlite magma is ~1 to 2 orders lower than that of mid-ocean ridge basalt (MORB) and comparable to the ultramobile pure carbonate melt. Using the measured viscosity values, we simulated the ascent and eruption process of kimberlite magma. We found that a minimum content of ~0.5 wt % water in the primitive magma is necessary to allow the ultrafast eruption process of kimberlite, thereby enabling the preservation of diamonds and high-pressure mineral inclusions transported by the magma.

58 GEOSCIENCES↗

Unlocking Larger Scales and Aspect Ratios in 3D Printed Glass: Coupling Active Mixing and UV Curing for Advanced Printability and Crack Resistance

Recent developments in additive manufacturing (AM) of glass via silica-filled inks have facilitated fabrication of previously unattainable geometries and compositions. However, the maximum processable size of 15 mm limits the use of these prints in applications such as optics. A key limitation lies in the trade-off between material printability and green strength: increasing silica content in the feedstock improves crack resistance and reduces shrinkage but results in dramatic changes in viscoelastic properties that hinder flowability. Here, this paper presents a novel approach that offers expanded versatility in processable size, feedstock formulation, and printing. Described here is a direct ink writing (DIW) system coupled with an active high-shear micromixer and UV light source, capable of simultaneously printing multiple inks with a wide range of rheological properties. Choice of silica sourc, solvent, UV-curable binder, and dispersant is used to tune the ink rheology and improve printability and mechanical properties. Imparting high shear with the micromixer while UV-curing the extrudate allows for increased ink viscosities and reduced nozzle diameters, enabling printing finer feature sizes. With these advances, thin-walled high-aspect ratio structures and a crack-free glass disk measuring 44 mm in diameter are demonstrated, an increase of 3× in the greatest dimension compared to current state-of-the-art.

36 MATERIALS SCIENCE↗

Techno-economic assessment and life cycle assessment of three potential pathways for biomass liquefaction

Biomass liquefaction can enhance corn-stover slurry transportation by reducing yield stress and overall viscosity, easing biomass handling, especially for slurries at high solids loading (≥30 % (w/v). Implementing a liquefaction process at the front-end of the biorefinery could reduce downtime caused by equipment fouling. Furthermore, liquefaction reduces process complexity by transforming high fiber mass into a pumpable slurry. Three scenarios for liquefaction of corn stover pellets: enzyme process, enzyme mimetic process, and a combined process of enzyme and enzyme mimetic were assessed. The ideal process will have the best performance-to-cost ratio, with minimal impact on the environment. Techno-economic and life-cycle analyses were used to understand the economic and environmental effects of the three liquefaction scenarios and to determine if one or more can achieve cost parity. Here, results show that the enzyme process had the best performance in terms of cost and emissions (gCO 2 -equivalent/kg slurry)compared to the other scenarios.

09 BIOMASS FUELS↗

Structural distortion and electron redistribution in dual-emitting gold nanoclusters

Deciphering the complicated excited-state process is critical for the development of luminescent materials with controllable emissions in different applications. Here we report the emergence of a photo-induced structural distortion accompanied by an electron redistribution in a series of gold nanoclusters. Such unexpected slow process of excited-state transformation results in near-infrared dual emission with extended photoluminescent lifetime. We demonstrate that this dual emission exhibits highly sensitive and ratiometric response to solvent polarity, viscosity, temperature and pressure. Thus, a versatile luminescent nano-sensor for multiple environmental parameters is developed based on this strategy. Furthermore, we fully unravel the atomic-scale structural origin of this unexpected excited-state transformation, and demonstrate control over the transition dynamics by tailoring the bi-tetrahedral core structures of gold nanoclusters. Overall, this work provides a substantial advance in the excited-state physical chemistry of luminescent nanoclusters and a general strategy for the rational design of next-generation nano-probes, sensors and switches.

77 NANOSCIENCE AND NANOTECHNOLOGY↗

Edge turbulence evolution and intermittency development near the density limit on the HL-2A tokamak

The development of intermittent non-Gaussian processes is researched in the edge turbulence of ohmically heated HL-2A discharges approaching the density limit. As the density limit is approached, the E × B shear flow at the last closed flux surface (LCFS) weakens, a strong positive skewness develops in the scrape-off layer (SOL), and negative skewness develops inside the LCFS of turbulent density fluctuations. A conditional averaging analysis confirms more frequent increased amplitude positive (negative) going density fluctuation activity in the SOL (inside the LCFS) as the density limit is approached. The measured turbulent stress across the edge, LCFS, and SOL region is decomposed into diffusive and residual stress components, and the nonlinear exchange of kinetic energy between the turbulence and the low-frequency E × B shear flow is determined. Residual stress acts to amplify the E × B flow at the LCFS, while the diffusive stress acts to dissipate the flow just inside this region, at the interface between the core plasma and the LCFS. The relative strength of the flow drive associated with the residual stress weakens as the density limit is approached, while the turbulent viscosity associated with the diffusive stress increases at high density. The adiabatic parameter, k z 2 v th 2 / ω ν e , drops significantly to about 0.5 in the SOL when the density limit is approached, indicating a transition from the adiabatic regime to the hydrodynamic regime due to increased collisionality. Such changes enhance the particle transport through the nonadiabatic electron response and hence should end in a stronger edge cooling at fixed edge plasma heat flux.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Simulating Atmospheric Boundary Layer Turbulence with Nek5000/RS

We present large-eddy-simulation (LES) modeling approaches for the simulation of atmospheric boundary layer turbulence that are of direct relevance to wind energy production. In this report, we study a GABLS benchmark problem using 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 [1, 2]. We demonstrate 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. The model fidelity and scaling performance of Nek5000/RS on DOE’s leadership computing platforms in comparison to those of AMR-Wind, a block-structured second-order finite-volume code with adaptive-mesh-refinement capabilities, are discussed in [3].

17 WIND ENERGY↗

Volatility Basis Set Distributions and Viscosity of Organic Aerosol Mixtures: Insights from Chemical Characterization Using Temperature-Programmed Desorption–Direct Analysis in Real-Time High-Resolution Mass Spectrometry

Quantitative assessment of gas-particle partitioning of individual components within complex atmospheric organic aerosols (OA) mixtures is critical for predicting and comprehending the formation and evolution of OA particles in the atmosphere. This investigation leverages previously documented data obtained through a temperature programmed desorption - direct analysis in real time – high resolution mass spectrometry (TPD-DART-HRMS) platform. This methodology facilitates the bottom-up construction of volatility basis set (VBS) distributions for constituents found in three biogenic secondary organic aerosol (SOA) mixtures produced through the ozonolysis of a-pinene, limonene, and ocimene. The apparent enthalpies (ΔH*, kJ mol -1 ) and saturated vapor mass concentrations (C T *, µg∙m -3 ) of individual SOA components, determined as a function of temperature (T, K), facilitated an assessment of changes in VBS distributions and gas-particle partitioning with respect to T and atmospheric total organic mass loadings (tOM, µg∙m -3 ). Further, the VBS distributions reveal distinct differences in volatilities among monomers, dimers, and trimers, enabling their categorization into separate volatility bins. At the ambient temperature of T = 298 K, only monomers efficiently partition between gas and particle phases across a broad range of atmospherically relevant total organic mass loadings (tOM) values of 1–100 µg∙m -3 . Partitioning of dimers and trimers becomes notable only at T > 360 K and T > 420 K, respectively. The viscosity of SOA mixtures is assessed using a bottom-up calculation approach, incorporating the input of elemental formulas, ΔH*, C T *, and particle-phase mass fractions of the SOA components. Through this approach, we are able to accurately estimate the variations in SOA viscosity that result from the evaporation of its components. These variations are, in turn, influenced by atmospherically relevant changes in tOM and T. Comparison of the calculated SOA viscosity and diffusivity values with literature reported experimental results shows close agreement, thereby validating the employed calculation approach. These findings underscore the significant potential for TPD-DART-HRMS measurements in enabling the untargeted analysis of organic molecules within OA mixtures. This approach facilitates quantitative assessment of their gas-particle partitioning and allows for the estimation of their viscosity and condensed-phase diffusion, thereby contributing valuable insights to atmospheric models.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Surface-Functionalized Cellulose Nanocrystals as Nanofillers for Crosslinking Processes: Implications for Thermosetting Resins

Understanding the response of fillers in the epoxy resin crosslinking process and characterizing polymer–filler dynamics are the key features that guide the engineering of new thermosetting resin composites. Here, in this work, X-ray photon correlation spectroscopy (XPCS) is used as a thermal analysis tool to track the microscopic changes occurring during the cure of a functionalized cellulose nanocrystal (mCNC)–epoxy composite. In contrast, traditional methods such as differential scanning calorimetry (DSC) and curing rheology are used to understand the kinetics and properties on macroscopic length scales. Of interest is the influence of the mCNC on the curing kinetics and properties of the thermosetting resin. Two levels of modification (increasing hydrophobicity) were chosen to observe the effect of functionalization. Before cure, the highly functionalized CNC (mCNC3) shows a 44% increase in complex viscosity (η*), while the less functionalized CNC (mCNC2) shows a η* value similar to that of the neat resin. As the cure cycle progresses, results from DSC, rheology, and XPCS further show the enhancement in dispersion for mCNC3. The results show a clear difference in the maximum drift velocity, maximum heat flow, and complex viscosity during the ramp to the isothermal cure temperature (T cure ). Such results suggest that mCNC3 contains a well-dispersed network of particles due to the higher level of functionalization. During T cure , a transition in elastic modulus (G') occurs only for the highly functionalized CNC particle system. We believe that heat-induced aggregation occurs, and the crosslinked resin ultimately dominates the macroscopic properties of the final cured system for all samples. The results from the three techniques are in good agreement and showcase XPCS as a beneficial experimental tool for characterizing the microscopic dynamics of particulate-filled thermosetting resins. Hence, we envision this to be a fundamental curing study for the design of thermosetting resin composites.

36 MATERIALS SCIENCE↗

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↗