Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “dynamic loading”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,765 records · Page 98

Pt Particles on a Dynamic TiO 2 Support in Near-Ambient Conditions−Disentangling Size, Pressure, and Support Effects

Platinum particles on reducible oxides are known to form complex and highly dynamic catalyst systems at elevated pressures and temperatures, often adopting active structures that differ from those found at room temperature and under ultrahigh vacuum (UHV). Here, we study the oxidation and structural evolution of subnanometer Pt clusters and nanoparticles supported on rutile TiO 2 (110) across an oxygen pressure range from UHV to 0.1 mbar, using near-ambient pressure X-ray photoelectron spectroscopy (NAP-XPS), scanning tunneling microscopy (STM) under UHV and NAP conditions, and low-energy ion scattering (LEIS). Our results reveal distinct differences in oxidation behavior and thermal stability between Pt nanoparticles and clusters, which are further modulated by the support stoichiometry and oxygen pressure. Small Pt clusters become oxidized even at room temperature but are susceptible to accelerated sintering in 0.1 mbar O 2 at elevated temperatures. In contrast, well-crystallized Pt nanoparticles on near-stoichiometric TiO 2 show weaker oxidation. On a reduced, defective TiO 2 support, Pt instead quickly becomes deeply buried by new titania layers, which are formed during support reoxidation. This process appears to result primarily from interactions of the support with the gas phase, unlike the classical, selflimited encapsulation that is induced by the strong metal−support interaction (SMSI). Finally, we address the full complexity of real catalysts in a direct side-by-side comparison of the single-crystalline model system with a Pt-loaded TiO 2 powder catalyst (P25). We conclude that the stoichiometry of the model supports must be carefully chosen and controlled to accurately reproduce the expected state of powder supports during redox reactions.

metal nanoparticles↗

A Transported Livengood–Wu Integral Model for Knock Prediction in Computational Fluid Dynamics Simulation

This work describes the development of a transported Livengood–Wu (L–W) integral model for computational fluid dynamics (CFD) simulation to predict autoignition and engine knock tendency. The currently employed L–W integral model considers both single-stage and two-stage ignition processes, thus can be generally applied to different fuels such as paraffin, olefin, aromatics, and alcohol. The model implementation is first validated in simulations of homogeneous charge compression ignition (HCCI) combustion for three different fuels, showing good accuracy in prediction of autoignition timing for fuels with either single-stage or two-stage ignition characteristics. Then, the L–W integral model is coupled with G-equation model to indicate end-gas autoignition and knock tendency in CFD simulations of a direct-injection spark-ignition engine. This modeling approach is about 10 times more efficient than the ones that based on detailed chemistry calculation and pressure oscillation analysis. Two fuels with same Research Octane Number (RON) but different octane sensitivity are studied, namely, Co-Optima alkylate and Co-Optima E30. Feed-forward neural network model in conjunction with multivariable minimization technique is used to generate fuel surrogates with targets of matched RON, octane sensitivity, and ethanol content. The CFD model is validated against experimental data in terms of pressure traces and heat release rate for both fuels under a wide range of operating conditions. The knock tendency—indicated by the fuel energy contained in the autoignited region—of the two fuels at different load conditions correlates well with the experimental results and the fuel octane sensitivity, implying the current knock modeling approach can capture the octane sensitivity effect and can be applied to further investigation on composition of octane sensitivity.

33 ADVANCED PROPULSION SYSTEMS↗

Molecular Alignment of a Meta-Aramid on Carbon Nanotubes by In Situ Interfacial Polymerization

Molecularly organized nanocomposites of polymers and carbon nanotubes (CNTs) have great promise as high-performance materials; in particular, conformal deposition of polymers can control interfacial properties for mechanical load transfer, electrical or thermal transport, or electro/chemical transduction. However, controllability of polymer-CNT interaction remains a challenge with common processing methods that combine CNTs and polymers in melt or in solution, often leading to nonuniform polymer distribution and CNT aggregation. In this work, we demonstrate CNTs within net-shape sheets can be controllably coated with a conformal coating of meta-aramid by simultaneous capillary infiltration and interfacial polymerization. We determine that π-interaction between the polymer and CNTs results in chain alignment parallel to the CNT outer wall. Subsequent nucleation and growth of the precipitated aramid forms a smooth continuous layered sheath around the CNTs. These findings motivate future investigation of mechanical properties of the resulting composites, and adaptation of the in situ polymerization method to other substrates.

77 NANOSCIENCE AND NANOTECHNOLOGY↗

Demonstration of Complete Recycling Processes of Reversible Epoxies Using Solar Energy Conversion

Reversible epoxies using the Diels–Alder chemistry enables recycling processes through depolymerizing the polymer at higher temperature and then repolymerizing upon cooling. Compared to conventional bulk heating, photothermal heating can save time and resource and, consequently, reduce costs to reach an elevated temperature for recycling processes of the reversible epoxies. In previous studies, self‐healing of cracks and reattachments of two broken pieces have been presented using a laser; however, recycling of a sample as a whole is not feasible by using such a point light source. Herein, complete recycling processes are demonstrated utilizing an area light source, i.e., sunlight. Reversible epoxies are incorporated with carbon black and refractory plasmonic titanium nitride nanoparticles (NPs). Under concentrated (10 times) sunlight, they can generate sufficient heat (≈140 °C) to completely liquefy, reprocess, and reshape the samples multiple times. Recycling processes are validated by evaluation of mechanical properties for each cycle. Using an integrated experimental and theoretical approach, photothermal performance is investigated in terms of the dispersion and loading of photothermal NPs in the matrix, as well as the sample thickness. In this study, an insight is provided into the design of polymer/photothermal nanomaterial composites which can be sustainably recycled using abundant solar energy.

14 SOLAR ENERGY↗

Advanced blade-shaped thermal energy storage device: Development and application

Thermal energy storage (TES) using phase change materials (PCMs) is a promising approach for capturing and reusing excess thermal energy, yet widespread adoption is limited by low thermal conductivity, bulky configurations, and inadequate scalability. Here, this study presents a modular, blade-shaped TES prototype designed to address these challenges. The device integrates a lightweight aluminum shell, an embedded serpentine coil for active or passive heat exchange, and a cost-effective corrugated metal mesh for enhanced PCM thermal conductivity. With thickness-to-length and thickness-to-width ratios of 0.03 and 0.08, respectively, the blade-shaped TES achieves a compact, modular form factor suitable for space-constrained applications. Experimental testing demonstrated the efficient charge and discharge behavior of blade-shaped TES, capturing PCM superheating, phase-change transitions, and subcooling dynamics, with charging and discharging efficiencies of 94.9% and 94.6%, respectively. Also, the system can potentially achieve higher energy density than that of conventional TES designs. When integrated into a household refrigerator during the study, three blade-shaped TES modules successfully shifted 100% of peak-time compressor operation to off-peak hours, reducing energy consumption while maintaining more stable compartment temperatures. The blade-shaped TES's thin geometry, modularity, and enhanced thermal performance support scalable deployment across residential, commercial, and industrial applications, providing a versatile, cost-effective solution for high-efficiency, demand-flexible thermal energy management.

Blade-shaped↗

Control co-design of 13 MW downwind two-bladed rotors to achieve 25% reduction in levelized cost of wind energy

Wind energy is recognized worldwide as cost-effective and environmentally friendly and is among the fastest-growing sources of electrical energy. To further decrease the cost of wind energy, wind turbines are being designed at ever larger scales, which is challenging due to greater structural loads and deflections. Large-scale systems such as modern wind turbines increasingly require a control co-design approach, whereby the system design and control design are performed in a more integrated fashion. We overview a two-bladed downwind morphing rotor concept that is expected to lower the cost of energy at wind turbine sizes beyond 13 megawatts (MW) compared with continued upscaling of traditional three-bladed upwind rotor designs. We describe an aero-structural-control co-design process that we have used in designing such extreme-scale wind turbines, and we discuss how we were able to achieve a 25% reduction in levelized cost of energy for our final turbine design compared to a conventional upwind three-bladed rotor design.

17 WIND ENERGY↗

A non-volatile thermal switch for building energy savings

Compared with traditional static insulation, a thermally switchable building envelope could reduce annual heating and cooling loads by intermittently coupling to the outside environment when beneficial. Here, we demonstrate a voltage-actuated, contact/non-contact thermal switch that meets the unique challenges of this application. The switch is non-volatile, consuming electricity only briefly while switching and none to hold steady state. The switch ratio is 12, the off state has a low effective thermal conductivity of 0.045 W m -1 K -1 , comparable to fiberglass insulation, and the performance is stable over 1,000 switching cycles. Numerical simulations using real-world climate data show that combining this thermal switch with a thermal storage layer in a building envelope can yield annual energy savings of 9%-55% (heating) and 17%-76% (air conditioning), depending on the climate zone. The greatest benefits are realized when the exterior temperature crosses well above and below the desired interior temperature within a single 24 h period.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

A transported Livengood-Wu integral model for knock prediction in CFD simulation

This work describes the development of a transported Livengood-Wu (L-W) integral model for computational fluid dynamics (CFD) simulation to predict auto-ignition and engine knock tendency. The currently employed L-W integral model considers both single-stage and two-stage ignition processes, thus can be generally applied to different fuels such as paraffin, olefin, aromatics and alcohol. The model implementation is first validated in simulations of homogeneous charge compression ignition combustion for three different fuels, showing good accuracy in prediction of auto-ignition timing for fuels with either single-stage or two-stage ignition characteristics. Then, the L-W integral model is coupled with G-equation model to indicate end-gas auto-ignition and knock tendency in CFD simulations of a direct injection spark ignition engine. This modeling approach is about 10 times more efficient than the ones that based on detailed chemistry calculation and pressure oscillation analysis. Two fuels with same Research Octane Number (RON) but different octane sensitivity are studied, namely Co-Optima Alkylate and Co-Optima E30. Feed-forward neural network model in conjunction with multi-variable minimization technique is used to generate fuel surrogates with targets of matched RON, octane sensitivity and ethanol content. The CFD model is validated against experimental data in terms of pressure traces and heat release rate for both fuels under a wide range of operating conditions. The knock tendency indicated by the fuel energy contained in the auto-ignited region of the two fuels at different load conditions correlates well with the experimental results and the fuel octane sensitivity, implying the current knock modeling approach can capture the octane sensitivity effect and can be applied to further investigation on composition of octane sensitivity.

Yue, Zongyu↗

Comparison of aircraft measurements during GoAmazon2014/5 and ACRIDICON-CHUVA

The indirect effect of atmospheric aerosol particles on the Earth's radiation balance remains one of the most uncertain components affecting climate change throughout the industrial period. The large uncertainty is partly due to the incomplete understanding of aerosol–cloud interactions. One objective of the GoAmazon2014/5 and the ACRIDICON (Aerosol, Cloud, Precipitation, and Radiation Interactions and Dynamics of Convective Cloud Systems)-CHUVA (Cloud Processes of the Main Precipitation Systems in Brazil) projects was to understand the influence of emissions from the tropical megacity of Manaus (Brazil) on the surrounding atmospheric environment of the rainforest and to investigate its role in the life cycle of convective clouds. During one of the intensive observation periods (IOPs) in the dry season from 1 September to 10 October 2014, comprehensive measurements of trace gases and aerosol properties were carried out at several ground sites. In a coordinated way, the advanced suites of sophisticated in situ instruments were deployed aboard both the US Department of Energy Gulfstream-1 (G1) aircraft and the German High Altitude and Long-Range Research Aircraft (HALO) during three coordinated flights on 9 and 21 September and 1 October. Here, we report on the comparison of measurements collected by the two aircraft during these three flights. Such comparisons are challenging but essential for assessing the data quality from the individual platforms and quantifying their uncertainty sources. Similar instruments mounted on the G1 and HALO collected vertical profile measurements of aerosol particle number concentrations and size distribution, cloud condensation nuclei concentrations, ozone and carbon monoxide mixing ratios, cloud droplet size distributions, and downward solar irradiance. We find that the above measurements from the two aircraft agreed within the measurement uncertainties. The relative fraction of the aerosol chemical composition measured by instruments on HALO agreed with the corresponding G1 data, although the total mass loadings only have a good agreement at high altitudes. Furthermore, possible causes of the discrepancies between measurements on the G1 and HALO are examined in this paper. Based on these results, criteria for meaningful aircraft measurement comparisons are discussed.

47 OTHER INSTRUMENTATION↗

Hydropower Black Start: A Guidebook for Retrofitting Grid Dependent Hydropower

Not all United States (US) hydropower plants were designed to provide black start, but they are increasingly needed to uphold resilience in the evolving electric grid. This guidance is designed to help understand the minimal retrofits required for grid dependent hydropower (GDH) plants behind the point of interconnection (POI). For distribution connected hydropower plants or those with dedicated cranking paths, such upgrades can be sufficient for the plant to provide black start. For others, more coordination with the transmission system operator will be needed. This guidebook answers a number of questions relevant to retrofitting hydropower plants with black start capabilities. For example, the guidebook answers: • How flexible do the wicket gate controls need to be? • Who needs to do hydro governor model validation, why, and how? • How robust and flexible do the excitation and AVR controls need to be? • What protection settings need to be adjusted? • What relay(s) will need to be bypassed or overridden and at what risk? • What is the electrical energy demand of the station load or auxiliary power systems? • What should the strategy to energize transformer(s) along cranking path to address inrush currents be? • How should the critical load restoration be sequenced? In addition to outlining the specifications that hydropower plants need to meet for each component to be able to perform black start, this guidebook provides a set of case studies for specific upgrades needed at actual plants. Between the case studies of plants that have already performed black start retrofits and the examples of how this guidebook can be applied to scope future retrofits, five key themes have been identified for retrofit needs. 1. Protection needs “black start” mode: hydropower plants that are not designed with black start capabilities will have protections that prevent them from interconnecting to a “dead bus.” These protections will need to be overridden in every retrofit case and a separate black start mode should be established so that operators can safely switch between black start and grid connected modes, minimizing the risk to the plant. 2. Wicket gates need modern controls: digital governors accelerate the parameter tuning process and gate position sensors improve controllability, so plants with mechanical governors should be upgraded. Furthermore, a black start and islanding mode should be established for controls to maximize plant performance. 3. Robust excitation support: the DC system or excitation generator needs to be reliable enough to form and sustain the rotor electromagnetic field. These systems are typically undersized in plants that were not designed for black start, so they will need to be upgraded. 4. Turbine-governor model validation and operator training: validation of a standard hydro governor model is needed to characterize the dynamic response (i.e., inertial and primary frequency response) of the GDH. This is required for control development and old hydropower plants often have outdated or incorrect models. Operator training is also typically required to ensure the hardware retrofits are utilized correctly during the black start process. 5. Transformer and cranking path energization: any upgradation and control adjustment in front of the POI will depend upon the existing interconnection. Coordination with the transmission or distribution operator may be required.

13 HYDRO ENERGY↗

Chemical Modeling of ATR Fuel in DOE Standard Canisters with Borated Stainless Steel Corrosion

Road-ready and final disposition packaging configurations for the ATR fuel dictates storage within helium backfilled sealed DOE standard canisters. The packaging configuration for the DOE canisters places 10 ATR elements within a Type 1a basket, and three baskets are loaded within each DOE canister. In order to provide a criticality control mechanism, these stainless baskets are proposed to be built with 1-2 weight percentage of boron. This report seeks to identify effects of corrosion of this basket material over time. The primary reaction found through surveying literature was the oxidation of the iron in steel by water to produce some hydrogen. This reaction was built into the previous chemical model of the DOE canister to assess additional hydrogen and pressure build up over time. Literature for non-borated stainless steel, and carbon steel has shown that irradiation increases the rate of corrosion, and it is assumed a similar effect would occur with borated stainless steel. Due to the wide range of results in literature, and the variation of the corrosion with boron content, a wide range of corrosion rates was tested in the model from 80 nm/yr up to 1800 nm/yr.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Novel data interpretation method for DIII-D divertor retarding field energy analyzer with 3-D particle-in-cell simulations

A novel data interpretation process that utilizes comprehensive particle-in-cell (PIC) simulations is developed for the new retarding field energy analyzer (RFEA) currently being constructed at DIII-D for the lower divertor using the Divertor Material Evaluation System. Furthermore, this probe is expected to survive a heat load of up to 100 MW/m 2 for up to 5 s and reliably measure the main ion temperature (T i ) on the divertor target ranging from 10 to 200 eV. These extreme conditions posed significant engineering limitations on the probe geometry, thus extensive validation work has been performed. The conventional fitting method for the RFEA I–V characteristics is based on a simplified 1-D model without considering the ion space charge inside the probe cavity and may not be sufficient for probes designed for the DIII-D divertor environment. In this article, a more realistic description of the particle propagation process within the RFEA cavity is achieved by including both 3-D geometric effects and ion space charge in the PIC simulations, and the capability to reconstruct the ion energy distribution functions is demonstrated with reasonable consistency.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

A GPU-based compressible combustion solver for applications exhibiting disparate space and time scales

High-speed chemically active flows pose significant computational challenges due to their disparate space and time scales, with stiff chemistry often dominating simulation time. While modern scientific computing programs achieve exascale performance by leveraging graphics processing units (GPUs), existing GPU-based compressible combustion solvers face critical limitations in memory management, load balancing, and handling the highly localized nature of chemical reactions. To this end, we present a high-performance compressible reacting flow solver built on the AMReX framework and optimized for multi-GPU settings. Here, our approach addresses three GPU performance bottlenecks: memory access patterns through column-major storage optimization, computational workload variability via a bulk-sparse integration strategy for chemical kinetics, and multi-GPU load distribution for adaptive mesh refinement applications. The solver adapts existing matrix-based chemical kinetics formulations to multi-grid contexts. Using representative combustion applications, including 2D and 3D detonations and a 3D jet-in-crossflow configuration, we demonstrate 1.4–5× performance improvements over initial implementations on an in-house cluster of NVIDIA H100 GPUs, and near-ideal weak scaling on the Frontier supercomputer (Oak Ridge Leadership Computing Facility) with up to 1024 AMD Instinct MI250X GPUs. Roofline analysis reveals substantial improvements in arithmetic intensity for both convection (∼ 10 ×) and chemistry (∼ 4 ×) routines, confirming efficient utilization of GPU memory bandwidth and computational resources.

42 ENGINEERING↗

RADAMS: Resilient and adaptive alert and attention management strategy against Informational Denial-of-Service (IDoS) attacks

Attacks exploiting human attentional vulnerability have posed severe threats to cybersecurity. In this work, we identify and formally define a new type of proactive attentional attacks called Informational Denial-of-Service (IDoS) attacks that generate a large volume of feint attacks to overload human operators and hide real attacks among feints. Here, we incorporate human factors (e.g., levels of expertise, stress, and efficiency) and empirical psychological results (e.g., the Yerkes-Dodson law and the sunk cost fallacy) to model the operators’ attention dynamics and their decision-making processes along with the real-time alert monitoring and inspection. To assist human operators in dismissing the feints and escalating the real attacks timely and accurately, we develop a Resilient and Adaptive Data-driven alert and Attention Management Strategy (RADAMS) that de-emphasizes alerts selectively based on the abstracted category labels of the alerts. RADAMS uses reinforcement learning to achieve a customized and transferable design for various human operators and evolving IDoS attacks. The integrated modeling and theoretical analysis lead to the Product Principle of Attention (PPoA), fundamental limits, and the tradeoff among crucial human and economic factors. Experimental results corroborate that the proposed strategy outperforms the default strategy and can reduce the IDoS risk by as much as 20%. Besides, the strategy is resilient to large variations of costs, attack frequencies, and human attention capacities. We have recognized interesting phenomena such as attentional risk equivalency, attacker’s dilemma, and the half-truth optimal attack strategy.

97 MATHEMATICS AND COMPUTING↗

Variability of MHD instabilities in benign termination of high-current runaway electron beams in the JET and DIII-D tokamaks

Benign termination, in which magnetohydrodynamic (MHD) instabilities deconfine runaway electrons (REs) following hydrogenic injections, is a promising strategy for mitigating dangerous RE loads after disruptions. Recent experiments on the Joint European Torus (JET) have explored this scenario at higher pre-disruptive plasma currents than are achievable on other devices, revealing challenges in obtaining benign terminations at I p ≥2.5 MA. This work analyzes the evolution of these high-current RE beams and their terminating MHD events using fast magnetic sensor measurements and EFIT equilibrium reconstructions for approximately 40 JET and 20 DIII-D tokamak discharges. On JET, unsuccessful non-benign terminations occur at low edge safety factor (q edge ≈ 2), and are preceded by intermittent, non-terminating MHD events at higher rational qedge. Trends in the internal inductance I i indicate more peaked RE current profiles in the high-I p non-benign population, which may hinder successful recombination through re-ionization of the companion plasma. In contrast, benign terminations on JET typically occur at higher q edge ≥3 and exhibit less peaked RE current profiles. DIII-D displays a broader range of terminating edge safety factors, again correlated with the measured values. Across both tokamaks, the RE current peaking is therefore found to determine which MHD instability boundary is encountered, a result confirmed by linear resistive MHD modeling with the CASTOR3D code. Measured growth rates are similar for benign and non-benign cases, indicating that ideal MHD timescales at low density after hydrogenic injection do not alone explain efficient RE deconfinement. Instead, non-benign cases are most readily characterized by their comparably lower overall MHD perturbation amplitudes δB. These observations suggest that the interplay between ideal and resistive dynamics governs the termination process, with implications for extrapolating benign RE termination to high-I p reactor scenarios.

MHD instabilities↗

Effect of void positioning on the detonation sensitivity of a heterogeneous energetic material

We show although it is well-established that voids profoundly influence the initiation and reaction behaviors of heterogeneous energetic materials such as polymer-bonded explosives (PBX) and propellants, there has been little study of how void location in different constituents in the microstructures of such materials affect the macroscale behavior. Here, we use three-dimensional (3D) mesoscale simulations to study how void placement within the reactive grains versus the polymer binder influences the shock-to-detonation transition (SDT) in a polymer-bonded explosive. The material studied here has a microstructure comprised of 75% PETN (pentaerythritol tetranitrate) grains and 25% HTPB (hydroxyl-terminated polybutadiene) polymer binder by volume. Porosities up to 10% in the form of spherical voids distributed in both the grains and polymer are considered. An Arrhenius reactive burn relation is used to model the chemical kinetics of the PETN grains under shock loading, thereby resolving the heterogeneous detonation behavior of the PBX. The influence of void location on the shock initiation sensitivity of the material is quantitatively ranked by comparing the predicted run distance to detonation (RDD) for each sample. The analysis includes inherent quantification of uncertainties arising from the stochastic variations in the microstructure morphologies and void distributions by using statistically equivalent microstructure sample sets (SEMSS), leading to probabilistic formulations for the RDD as a function of shock pressure. The calculations reveal that the location of voids in the composite microstructure significantly affects the RDD. Specifically, voids exclusively within the grains cause the PBX to be more sensitive (having shorter RDD) than voids in the polymer binder. Unique probabilistic relationships are derived to map the probability of observing RDD for each void location material case, allowing for prediction of initiation behavior anywhere in the shock pressure – RDD space. These findings agree with trends reported in the literature.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Multiscale Electrothermal Design of a Modular Multilevel Converter for Medium-Voltage Grid-Tied Applications

As a key feature of modular multilevel converters (MMCs), a large number of semiconductor devices are employed in the converter and distributed over a stack of submodules. Each submodule has a strict temperature limit that imposes constraints on the operating range of the converter. Different loading conditions/losses in the submodules lead to unavoidable temperature variations inside the MMC, which consequently affect the system-level performance and reliability. This paper is focused on the electrothermal analysis and design of a medium-voltage silicon carbide (SiC)-based MMC system, from submodule power semiconductors to the overall MMC system integration of multiple submodules, for grid-connected applications. Loss calculations are performed to estimate the cooling requirements and aid thermal design at different levels. The performance of a forced-air cooling approach is analyzed within a numerical modeling framework. Maximum temperature of the SiC power modules is predicted using numerical tools. Thermal design of the MMC cabinet with various arrangements of air inlet(s) and outlet(s) is investigated and compared from a cooling performance perspective. A fully resolved model of the arrangement that yields optimal results is developed to accurately predict the temperature profile of the essential components in the submodules.

computational fluid dynamics↗

A constitutive model for sheared dense fiber suspensions

Here we propose a constitutive model to predict the viscosity of fiber suspensions, which undergoes shear thinning, at various volume fractions, aspect ratios, and shear stresses/rates. We calibrate the model using the data from direct numerical simulation and prove the accuracy by predicting experimental measurements from the literature. We use a friction coefficient decreasing with the normal load between the fibers to quantitatively reproduce the experimentally observed shear thinning in fiber suspensions. In this model, the effective normal contact force, which is directly proportional to the bulk shear stress, determines the effective friction coefficient. A rise in the shear stress reduces the effective friction coefficient in the suspension. As a result, the jamming volume fraction increases with the shear stress, resulting in a shear thinning in the suspension viscosity. Moreover, we extend the model to quantify the effects of fiber volume fraction and aspect ratio in the suspension. We calibrate this model using the data from numerical simulations for the rate-controlled shear flow. Once calibrated, we show that the model can be used to predict the relative viscosity for different volume fractions, shear stresses, and aspect ratios. The model predictions are in excellent agreement with the available experimental measurements from the literature. The findings of this study can potentially be used to tune the fiber size and volume fraction for designing the suspension rheology in various applications.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗