Start-up Flow Research Articles

Data-driven rheological characterization of stress buildup and relaxation in thermal greases

Thermal greases, often used as thermal interface materials, are complex paste-like mixtures composed of a base polymer in which dense metallic (or ceramic) filler particles are dispersed to improve the heat transfer properties of the material. They have complex rheological properties that impact the performance of the thermal interface material over its lifetime. We perform rheological experiments on thermal greases and observe both stress relaxation and stress buildup regimes. This time-dependent rheological behavior of such complex fluid-like materials is not captured by steady shear-thinning models often used to describe these materials. We find that thixo-elasto-visco-plastic (TEVP) and nonlinear-elasto-visco-plastic (NEVP) constitutive models characterize the observed stress relaxation and buildup regimes, respectively. Specifically, we use the models within a data-driven approach based on physics-informed neural networks (PINNs). PINNs are used to solve the inverse problem of determining the rheological model parameters from the dynamic response in experiments. These training data are generated by startup flow experiments at different (constant) shear rates using a shear rheometer. We validate the “learned” models by comparing their predicted shear stress evolution to experiments under shear rates not used in the training datasets. We further validate the learned TEVP model by solving a forward problem numerically to determine the shear stress evolution for an input step-strain profile. Meanwhile, the NEVP model is further validated by comparison to a steady Herschel–Bulkley fit of the material’s flow curve.

Molecular Processes Leading to Shear Banding in Entangled Polymeric Solutions.

The temporal and spatial evolution of shear banding during startup and steady-state shear flow was studied for solutions of entangled, linear, monodisperse polyethylene C3000H6002 dissolved in hexadecane and benzene solvents. A high-fidelity coarse-grained dissipative particle dynamics method was developed and evaluated based on previous NEMD simulations of similar solutions. The polymeric contribution to shear stress exhibited a monotonically increasing flow curve with a broad stress plateau at intermediate shear rates. For startup shear flow, transient shear banding was observed at applied shear rates within the steady-state shear stress plateau. Shear bands were generated at strain values where the first normal stress difference exhibited a maximum, with lifetimes persisting for up to several hundred strain units. During the lifetime of the shear bands, an inhomogeneous concentration distribution was evident within the system, with higher polymer concentration in the slow bands at low effective shear rate; i.e., γ˙<τR-1, and vice versa at high shear rate. At low values of applied shear rate, a reverse flow phenomenon was observed in the hexadecane solution, which resulted from elastic recoil of the molecules within the slow band. In all cases, the shear bands dissipated at high strains and the system attained steady-state behavior, with a uniform, linear velocity profile across the simulation cell and a homogeneous concentration.

On the thixotropy of mature fine tailings

The rheological behavior of mature fine tailings (MFTs) is investigated using transient and steady shear flow fields. The structure breakdown of intact MFT samples is examined by a startup flow experiment at various shear rates. The yield stress of MFTs is estimated by the steady shear stress values at low shear rates. Oscillatory shear (strain amplitude sweep) is also used to verify the obtained yield stress value. MFT samples exhibit thixotropy and a positive hysteresis loop at short shearing time intervals in increasing and decreasing stepwise shear rate tests. The observed hysteresis loops and thixotropy disappear by increasing the shearing time intervals, as the system reaches its equilibrium steady-state structure. The time-dependent rheological behavior of MFTs is quantified by a structural kinetics model through the dimensionless structure parameter, λ. [Toorman, “Modelling the thixotropic behaviour of dense cohesive sediment suspensions,” Rheol. Acta 36, 56–65 (1997).] The kinetic parameters are estimated based on steady-state stresses, elucidating the relative effects of shear rate and Brownian motion on buildup and breakdown of the structure. The flow behavior of MFTs predicted by the structural kinetics model is in agreement with the experimental data.

10 clinical tips for advancing patient safety when using syringe pump systems for microinfusion intravenous drug therapy.

10 clinical tips for advancing patient safety when using syringe pump systems for microinfusion intravenous drug therapy.

Hybrid bipolar VSC-DRU HVdc connection for offshore wind farms

This paper presents a hybrid bipolar high voltage direct current (HVdc) link for offshore wind farms (OWFs) connection. Voltage source converter (VSC) technology is used in one of the offshore poles, using a diode rectifier unit (DRU) in the other one. Therefore, economic and power-loss advantages of the DRU over the VSC are leveraged by changing one of the HVdc poles conventional technology. In the proposed control strategy, VSC reactive power is used to control OWF ac grid frequency, allowing the DRU commutation and active power transfer. VSC active power is used for OWF ac grid start-up and hybrid bipolar HVdc link power flow control. OWF ac grid voltage magnitude is controlled while the DRU pole is blocked, and it is automatically accommodated to the required value when the DRU is conducting. Moreover, VSC is also used as an active filter for the main DRU current harmonic components, thus removing the huge passive ac filters from the offshore platform. Poles balanced and unbalanced power flow control strategies are proposed for the hybrid bipolar HVdc, considering both bipolar with neutral return and rigid bipolar configurations. Simulation results validate the proposal which aims reducing cost and size of the OWF HVdc connection.

Research on Solidity of Horizontal-Axis Tidal Current Turbine

In this paper, the blade solidity of the tidal current turbine was investigated. Based on the blade element momentum theory, different design flow velocities were selected to design two blade types with different solidities. The geometric parameters of the blade were calculated using MATLAB programming, and the performance of two blades was compared in terms of the start-up flow rate, power generation and thrust by test experiment, which showed that the blade with higher solidity has better start-up performance and higher energy capture efficiency at low flow rates. The performance is better than that of the blade with low solidity, but due to the high solidity, the thrust is also high, which should be taken into account when installing the turbine.

The motion of a layer of yield-stress material on an oscillating plate

The motion of a finite layer of Bingham material on a solid plate that executes in-plane oscillations was reported previously by Balmforth et al. (J. Non-Newtonian Fluid Mech., vol. 158, issue 1–3, 2009, pp. 46–53). There, it was suggested that multiple yielded regions may arise within the material; this contrasts to start-up flow of the same material for which only one yielded region is generated. Here, we explore quantitatively the fluid physics of this oscillatory flow problem through analytical approximations and further numerical computation. Four new key topological properties concerning the generation and termination of the yielded regions are reported. It is shown that the existence of multiple yielded regions is equivalent to the layer never becoming entirely rigid during the periodic motion. For small inertia, the flow is approximately time-reversible with only a single yielded region generated at the plate. For large inertia, shear stress in the material decays rapidly as a function of distance from the plate. A thin zone of yielded material detaches periodically from the plate, and subsequently terminates within the layer. At high oscillation frequency, there can be any number $N$ of distinct rigid regions, satisfying $N= \lfloor 1- {\rm \pi}^{-1} \log B \rfloor$ where $B$ is the Bingham number. It is also shown that for $B&gt;0.5370$ , there are at most one yielded region and one rigid region throughout the motion. These theoretical results can be used as a basis for oscillatory rheometry, allowing for measurement of the yield stress using existing apparatus.

Design and fabrication of a differential-pressure optical fiber grating sensor for monitoring the flow rate of fluid.

A fiber Bragg grating (FBG) flow sensor is designed and fabricated, in which two FBGs are fixed on the front and other side of the metal diaphragm, and differential pressure is used to monitor the flow rate of fluid. The temperature sensitivity of these two FBGs is 0.030 and 0.029nm/°C, which is almost the same, suggesting that the influence of temperature on the flow measurement can be effectively eliminated. The static pressure sensitivity of these two FBGs can be up to 86.7nm/MPa and 68.6nm/MPa, respectively; accordingly, the static pressure sensitivity of the sensor overall is 155.3nm/MPa. Furthermore, the flow rate sensitivity is 0.00029L/s. This FBG flow sensor exhibits high sensitivity, high accuracy, and a low start-up flow rate. Furthermore, the cross effect between the temperature and strain on the sensing sensitivity is eliminated, which makes this FBG flow sensor suitable for real-time monitoring of the trace flow rate in oil and gas wells.

Microstructural evolution and reverse flow in shear-banding of entangled polymer melts.

The temporal and spatial evolution of shear banding under startup of shear flow was simulated for highly entangled, linear, monodisperse polyethylene melts of differing molecular weight, C750H1502, C1200H2402, and C3000H6002, using a high-fidelity coarse-grained dissipative particle dynamics method. It was determined that shear stress was dominated by segmental orientation of entangled strands at low shear rates, but at a critical shear rate below the reciprocal of the Rouse time, flow-induced disentanglement resulted in the onset of chain tumbling that reduced the average degree of orientation, leading to a regime of decreasing shear stress, with a commensurate onset of increasing average chain extension imposed by the strong flow kinematics that ultimately drove the steady-state shear stress higher. During startup of shear flow, shear band development began immediately after the maximum in the first normal stress difference, where distinct fast and slow bands formed. The slow bands consisted of relatively entangled and coiled molecules, whereas the fast bands consisted of more disentangled and extended chains that experienced quasiperiodic rotation/retraction cycles. The simulation results often exhibited a generation of temporary reverse flow, in which the local fluid velocity was temporarily opposite to that of the bulk flow direction, at the onset of the shear-banding phenomena; this effect was consistent with earlier experiments and theoretical results. The physical mechanism for the generation of reverse flow during shear-band formation was investigated and found to be related to the recoil of the molecules comprising the slow band. Overall, the phenomenon of shear banding appeared to arise due to flow-induced disentanglement from orientational ordering and segmental stretching that affected individual chains to different degrees, ultimately resulting in regions of relatively coiled and entangled chains that evolved into a slow band, whereas the locally disentangled chains, experiencing quasiperiodic stretch-rotation cycles, formed a fast band. The transitional period resulted in a kinematic instability that generated the temporary reverse-flow phenomenon.

CALCULATION OF THE HEAT TRANSFER COEFFICIENT IN THE OUTER BODY FOR A ROTATIONAL DETONATION ENGINE

Unsteady heat transfer characterization on the combustion surfaces of rotational detonation engines (RDE) is not well understood. It is generally thought that the complex nature of the unsteady, reacting, compressible fluid flow inside the combustion annulus of the RDE causes the convective heat transfer coefficient to be significantly higher than it is in other applications. Empirical models that have been used to analyze this strictly apply to steady flows where dimensionless groups can be employed. Because of the complex flows in RDEs, it is not likely that they will lend themselves to explanation using simple dimensionless parameters. In this work a numerical study is performed where the inside wall surface temperature distribution in the RDE outer body is systematically changed over a given range that would be characteristic of the startup flows inside an RDE. For each case, temperature distributions inside the outer containment wall of the RDE were calculated and compared with experimental data. The closest match was then used to directly calculate the convective heat transfer coefficient on the inside surface of the RDE. The calculation results show that the local heat transfer coefficient C&lt;sub&gt;H&lt;/sub&gt; is lower than what has been commonly used for these flows. The results show that C&lt;sub&gt;H&lt;/sub&gt; is within the range of normal forced convection of a gas.

High-speed videomicroscopy of sheared carbonyl iron suspensions

The postyield rheological regime is investigated in sheared magnetic field-responsive composites (i.e. carbonyl iron based magnetorheological fluids). When subjected to uniaxial DC fields, high-speed videomicroscopy techniques and dedicated image analysis tools demonstrate that dispersed magnetic microparticles self-assemble to form concentric layered patterns above a particular shear rate (). This critical shear rate for layer formation is dictated by a critical Mason number 1 that is associated to the destruction of the last doublet in the chain-like aggregates. The number of layers, mean width, percentage of occupation and mean period are found to be very weakly dependent on the shear rate in start-up shearing flow tests. Experimental data for the mean period are in good agreement with an energy minimization theory.

Influence of Cellulose Nanofibers on the Behavior of Pickering Emulsions. Part 1. Microscopy and Startup Flow Test

The dispersibility of flexible polymer chains present at the emulsion's interface between the dispersed and continuous phase has obvious effects on rheology and dielectric properties of the whole emulsion. Cellulose nanofiber (CNF)-based Pickering emulsions are good systems to research these properties with respect to their microscopic phase structure, dielectric, and rheological properties by using CNF as a water-dispersible Pickering emulsifier, liquid paraffin as an oil phase, and didodecyldimethylammonium bromide (DDAB) as a cationic auxiliary surfactant. The CNF and DDAB contents were systematically varied while the water-to-paraffin oil ratio was kept constant to discern the influence of the Pickering emulsifiers. Polarized optical microscopic images reveal that the droplets tend to shrink at higher CNF content but grow bigger when increasing the DDAB content, which is proved by fluorescence analysis of the CNF dispersibility with varying DDAB content. The dielectric damping exhibits a minimum, whose value decreases with increasing DDAB and CNF content. Increasing the DDAB content promotes the solubilization of CNF in the aqueous phase, which will increase the overall viscosity and yield points. Similarly, a higher CNF content leads to a higher viscosity and yield point, but at high DDAB contents, the viscosity function exhibits an S-shape at intermediate CNF contents. To evaluate the results further, they were compared with CNF dispersions (without oil phase), which showed a surfactant effect slightly on maximum stress but strongly on yield stress τy, indicating that DDAB can promote the formation of a CNF network rather than the viscosity of the whole system. This paper provides information on how a systematical variation of the composition influences morphology and physico-chemical interactions as detected by broadband dielectric spectroscopy and rheological behavior.

Leveraging spreadsheet analysis tool for electrically actuated start-up flow of non-Newtonian fluid in small-scale systems

In this article, we demonstrate the solution methodology of start-up electrokinetic flow of non-Newtonian fluids in a microfluidic channel having square cross-section using Spreadsheet analysis tool. In order to incorporate the rheology of the non-Newtonian fluids, we take into consideration the Ostwald-de Waele power law model. By making a comprehensive discussion on the implementation details of the discretized form of the transport equations in Spreadsheet analysis tool, and establishing the analytical solution for a special case of the start-up flow, we compare the results both during initial transience as well as in case of steady-state scenario. Also, to substantiate the efficacy of the proposed spreadsheet analysis in addressing the detailed flow physics of rheological fluids, we verify the results for several cases with the corresponding numerical results. It is found that the solution obtained from the Spreadsheet analysis is in good agreement with the numerical results—a finding supporting spreadsheet analysis's suitability for capturing the fine details of microscale flows. We strongly believe that our analysis study will open up a new research scope in simulating microscale transport process of non-Newtonian fluids in the framework of cost-effective and non-time consuming manner.

Peeling of linearly elastic sheets using complex fluids at low Reynolds numbers

We investigate the transient, fluid–structure interaction (FSI) of a non-Newtonian fluid peeling two Hookean sheets at low Reynolds numbers ( R e ≪ 1 ). Two different non-Newtonian fluids are considered – a simplified Phan–Thien–Tanner (sPTT) model fluid, and an inelastic fluid with shear-thinning viscosity (generalized Newtonian fluid). In the limit of small gap between the sheets, we invoke a lubrication approximation and numerically solve for the gap height between the two sheets during the start-up of a pressure-controlled flow. What we observe is that for an impulse pressure applied to the sheet inlet, the peeling front moves diffusively ( x f ∼ t 1 / 2 ) towards the end of the sheet when the fluid is Newtonian. However, when one examines a complex fluid with shear-thinning, the propagation front moves sub-diffusively in time ( x f ∼ t α , α < 0 . 5 ), but ultimately reaches the end faster due to an order of magnitude larger prefactor for the propagation speed. We provide scaling analyses and similarity solutions to delineate several regimes of peeling based on the sheet elasticity, flow Weissenberg number (for sPTT fluid), and shear-thinning exponent (for generalized Newtonian fluids). To conclude, this study aims to afford to the experimentalist a system of knowledge to a priori delineate the peeling characteristics of a certain class of complex fluids. • Modeled fluid–structure interaction of a complex fluid peeling two Hookean sheets. • Examined two different fluid constitutive rheological models. • Developed simulations and similarity solutions for the peeling front during start-up. • Quantified peeling time effects due to shear thinning. • Peeling time for a shear-thinning fluid is smaller than a Newtonian fluid.

Womersley's solution for the measurement of volume flow rates in transient laminar flow tubes

The characterization of transient flows within the Reynolds number range of Re = 10–100 and the Womersley number range α = 0.8–10 is required for the ongoing development of green monopropellant thrusters. However, at the ml/min scale flow rates of interest, these measurements are outside the capabilities of current commercial flow meters. It is proposed here that transient flows under the required dynamic conditions can be characterized via Womersley's solution for transient flow in a rigid tube. This solution method requires only the measured transient pressure gradient within a controlled laminar flow section and can be accomplished using existing commercial pressure measurement hardware. Experiments were performed where flow similarity was maintained with the ultimate thruster characterization application, but the radius of the flow passage was increased so that flows could be simultaneously characterized by both the proposed solution method and a commercial ultrasonic flow meter. It was shown that across the range of interest, applying Womersley's solution to a measured pressure gradient was an effective method of transient flow characterization. Additionally, it was shown that non-periodic flows can be characterized except for the initial flow startup transient with solution convergence times approximating the analytical solution to starting flow in a pipe. While these results were expected due to an experimental design matching the assumptions required for Womersley's analytic solution, this work demonstrates that this method is practically feasible as novel instrumentation enabling previously unobtainable measurements.

Topological Transition in Multicyclic Chains with Structural Symmetry Inducing Stress-Overshoot Phenomena in Multicyclic/Linear Blends under Biaxial Elongational Flow

Blends of multicyclic and linear polymers under biaxial elongational flow were analyzed using coarse-grained molecular dynamics simulations. The multicyclic/linear blends displayed overshoot in the normal stress difference at the start-up of biaxial elongational flow. This overshoot was steeper for multicyclic/linear blends than for our previously reported monocyclic/linear blends [T. Murashima, K. Hagita, and T. Kawakatsu, Macromolecules, 2021, 54, 7210]. Investigation of the origin of the overshoot in the multicyclic/linear blends revealed a different mechanism than that previously observed in our monocyclic/linear blends. Specifically, a "topological transition" mechanism comprising a morphological change from the open- to closed-ring state was observed in the multicyclic chains, but not in the monocyclic chains. This topological transition drastically changes the stress of the rings. Although the topological transition was also observed in asymmetric-multicyclic/linear blends, no stress overshoot appeared, owing to the asymmetry in the multicyclic chains. Therefore, we hypothesized that the structural symmetry in multicyclic chains is indispensable for overshoot behavior to occur. We determined that the topological transition in multicyclic chains with structural symmetry induces stress-overshoot behavior in multicyclic/linear blends under biaxial elongational flow.

Numerical study on heat transfer enhancement by viscoelastic fluid pulsating laminar flow in rectangular microchannel heat sinks

Efficient heat transfer technique at microscale is a long-term demand for development of microelectronic device. Simple-structured microchannel heat sinks (MHSs) have the advantages of easier fabrication processes and profound applicability, which are however limited by the laminar flow with lower heat transfer efficiency. Against this backdrop, a new type of heat transfer enhancement (HTE) techniques by viscoelastic fluid pulsating laminar flow (VPL flow) is proposed and evaluated from different aspects in this research. First of all, the effects of VPL flow on the thermal performance in a rectangular microchannel heat sink (RMHS) are numerically studied with the focus on the parametric effects on HTE. It indicates that the overall heat transfer performance can be significantly augmented by VPL flow and possesses an optimal enhancement under cases of exerted pulsation frequencies closed to the characteristic frequency of viscoelastic fluid flow. The characteristic frequency is determined from the Poiseuille start-up flow of the viscoelastic fluid flow and can be summarized by a power-law function with the relationship of Reynold number (Re) and Elasticity number (El). Moreover, the optimal HTE of RMHS increases with external pulsation amplitude and El but decreases with Re, which can also be summed up with them in a power-law function. Then, by analyzing of flow performance, the reversal flow and the resonance effects by VPL flow are regarded as the reasons for the HTE, which results from a balance between relaxation time of viscoelastic microstructures deformation and external pulsation period. Finally, the law of HTE by VPL flow is extended and verified in a carefully designed leaf-type microchannel heat sink (LMHS), suggesting a promising micro-scale heat transfer enhancement technology by VPL flow which however still needs further investigations on the optimal operating conditions in the future.

The starting mechanism study on rainfall-induced sloop loose source in strong earthquake area

The 5.12 Wenchuan earthquake in 2008 led to a large number of loose deposits on the slope on both sides of the gully in the formation area. Under the triggering action of heavy rainfall, the previously unrecognized slope loose source in formation areas easily failed and started earlier than the overall start of the gully debris flow, which has practical significance for realizing the early warning of the gully debris flow starting. The special loose source conditions and postearthquake geological environment change the hydraulic mechanism of traditional debris flow start-up, which brings great difficulties to the monitoring and early warning of postearthquake debris flows. Therefore, based on hydraulics, on the premise of building a hydraulic model of groundwater level change of loose slope accumulation body, the characteristics and rules of groundwater level change are analyzed with the aid of hydraulic seepage theories, and the action characteristics of hydrodynamic pressure and hydrostatic pressure on the accumulation body are quantitatively studied. Based on field exploration data and average slope gradient, we divided the accumulation body in a reasonable manner; completed the sliding force, anti-sliding force and residual-sliding force mechanical expression and establishment; and performed a slope stability analysis. The research results show that the slope starting mode is divided into the thrust-type landslide and the retrogressive landslide mode of the subsection disintegration. The hydraulic deposit mechanism analytical results can also reverse the critical rainfall condition of slope failure. Finally, taking the Yindongzi gully in Dujiangyan as an example, the hydraulic mechanism of typical slope deposits in the formation area was analyzed and combined with the historical debris flow events in the study area to verify the physical simulation test, which was consistent with the field investigation results after the disaster.

Transient flow patterns of start-up flow in round microcavities

Transient flow patterns of start-up flow in round microcavities

Start-up flow in a pipe of a double distributed-order Maxwell fluid

The distributed-order model is equivalent to weighting the fractional-order model on a given interval. In this paper, a double time-distributed order Maxwell model is developed and formulated to study the characteristics of a start-up pipe flow. The numerical solution of the distributed-order governing equation is obtained by the finite difference method. In order to test the convergence of numerical solutions, the truncation error of the proposed numerical method be calculated. The characteristics of the flow and the effects of fractional parameters and weight functions on the velocity distribution are discussed in detail. By comparing the velocity distribution diagrams, it can be inferred that the distributed-order model controlled by the weight function can describe the memory and elastic properties of fluid more flexibly.