Numerical investigation of the propagation of shock waves in rigid porous materials: flow field behavior and parametric study

A numerical parametric study of the flow field which develops when a planar shock wave impinge on a rigid porous material is presented. This study complements an earlier study (Levy et al. 1996a) where the values of some dominating parameters were estimated and the dependence of the resulting flow field on these values was not checked. Received 22 April 1996 / Accepted 5 January 1997     

Shock tube study of the dynamical behavior of granular materials

The purpose of the present study is to clarify both the compression phenomenon and the gas filtration effect that take place inside a granular medium when it is dynamically loaded by a shock wave. In order to measure the pore pressure and the total stress at different locations along the granular medium, pressure transducers were placed along the side-wall and at the end-wall of the shock tube test section, which was filled with the granular material. In order to elucidate the gas filtration effect, the results of two experiments with identical granular media but with and without filtration were compared. The gas filtration was eliminated by means of a thin plastic film, which was placed at the front edge of the granular medium. Based on this comparison quantitative information on the gas filtration and its role in the stress formation inside granular media of different material and length was obtained. Furthermore, curves of the dynamic compression and the Young moduli of the granular medium for the range of the operating conditions were reconstructed.     

Distance effect on the behavior of an impinging swirling jet by PIV and flow visualizations

The present paper deals with the problem of an impinging swirling jet against a wall. The study concerned a detailed experimental investigation on the jet–wall interaction using PIV and flow visualizations over a range of operating conditions in which the distance of the ducted propeller from the wall was changed. The influence of the impingement distance and the swirl number (i.e., ratio between the axial fluxes of the swirl and the axial momentum) as well as the interaction between the jet deformation and the perturbation induced on the wall is discussed in this paper.     

Characterization of powder flow properties from micron to nanoscale using FT4 powder rheometer and PT-X powder tester

The rapid development of nanotechnology has led to a need to further understand the physical characteristics of nanoparticles. In this paper, the flow characteristics of micro-nano alumina particles with different particle sizes were characterized. The FT4 powder rheometer and the PT-X powder tester were used to measure the compression, friction, and dynamic properties of powders. Powder compressibility increased significantly as the particle size decreased from 27 μm to 30 nm. Pressure distribution in the silo was measured and predicted by Janssen's theory, with errors mostly less than 10%. The basic flow energy and the specific energy of the three powders were 4983, 1734, and 244 mJ, and 6.80, 11.70, and 6.70 mJ/g, respectively, indicating that there was no linear relationship between the change in flowability and particle size. The dynamic properties of the powders change from particle-dominated to agglomeration-dominated as the particle size decreases. The conclusion is supported by the results of field emission scanning electron microscopy.     

Experimental and numerical analyses on the effect of increasing inflow temperatures on the flow mixing behavior in a T-junction

Thermal degradation of piping induced by high cycle thermal fatigue (HCTF) is of significant importance as operating Nuclear Power Plants (NPP) become older and lifetime extension activities are initiated. In particular, HCTF incidents related to turbulent thermal mixing of fluids in a T-junction piping system are not well understood and could not be adequately monitored using common thermocouple instrumentation. To investigate this phenomenon, an experimental T-junction test facility was commissioned at the University of Stuttgart, known as the Fluid Structure Interaction (FSI) test facility. The paper presents the experimental investigation and the corresponding numerical validation using the large eddy simulation (LES) method to study T-junction flow mixing. Three experimental test cases are investigated with temperature differences (∆T) of 51.5 K (Case 1), 76 K (Case 2) and 97 K (Case 3) between the mixing fluids. A constant mass flow rate ratio (main/branch) of 4:1 is maintained in all the investigated cases. Flow mixing is observed to be incomplete in all the cases, resulting in a thermally stratified flow with an oscillating stratification layer downstream of the T-junction. Mean temperature and root mean square (RMS) temperature fluctuations predicted by LES in the mixing region are found to be in good agreement with measurement data, with the exception of few positions. Amplitudes of temperature fluctuations are observed to be higher near the stratification layer, ranging from 6.3–9.9% of ∆T. Power spectral density (PSD) analyses of temperature fluctuations indicate no dominant frequency (spectral peak) under prevailing flow conditions, an important factor in thermal fatigue analysis, and the energy of these fluctuations are mainly contained in the frequency range of 0.1–2 Hz for all the investigated cases. LES is performed using the CFD software ANSYS CFX 14.0.     

Shear and elongational flow behavior of acrylic thickener solutions. Part II: effect of gel content

We have investigated the effect of crosslink density on shear and elongational flow properties of alkali-swellable acrylic thickener solutions using a mixing series of the two commercial thickeners Sterocoll FD and Sterocoll D as model system. Linear viscoelastic moduli show a smooth transition from weakly elastic to gel-like behavior. Steady shear data are very well described by a single mode Giesekus model at all mixing ratios. Extensional flow behavior has been characterized using the CaBER technique. Corresponding decay of filament diameter is also well fitted by the Giesekus model, except for the highest crosslink densities, when filament deformation is highly non-uniform, but the non-linearity parameter α, which is independent of the mixing ratio, is two orders of magnitude higher in shear compared to elongational flow. Shear relaxation times increase by orders of magnitude, but the characteristic elongational relaxation time decreases weakly, as gel content increases. Accordingly, variation of gel content is a valuable tool to adjust the low shear viscosity in a wide range while keeping extensional flow resistance essentially constant.     

The flow behavior of fiber suspensions in Newtonian fluids and polymer solutions.

This is the second part of a study examining the mechanical properties and capillary flow of fiber suspensions in Newtonian fluids and in polymer solutions. In part I results for the viscous and elastic properties of the fiber suspensions were presented and it was shown that the fiber suspensions exhibited normal stresses in Newtonian as well as in viscoelastic suspending media. It was thus expected that circulating secondary flows would occur near the entrance to a capillary. Four types of fillers (glass, carbon, nylon and vinylon fibers) suspended in glycerin, HEC solutions and Separan solutions were investigated. The entrance flow patterns were visualized and the pressure fluctuations measured. The visualization enabled the eddies occurring in the fiber suspensions in Newtonian fluids to be analysed and classified into two tpyes. The results from the flow visualization were correlated with the pressure fluctuations. Empirical equations for the tube length correction factor due to the elasticity were obtained.     

The flow behavior of fiber suspensions in Newtonian fluids and polymer solutions.

The main purpose of this study was to examine the viscous and elastic properties and capillary flow of fiber suspensions in Newtonian fluids as well as in polymer solutions. The fillers used were glass, carbon, nylon and vinylon fibers. Glycerin was used as a Newtonian suspending medium and HEC and Separan solutions as viscoelastic suspending media. The viscosity and the first normal-stress difference were measured using a coaxial cylindrical rotating viscometer and a parallel-plate rheogoniometer respectively. The influence of the concentration, aspect ratio, diameter and flexibility of the fibers on the viscous and elastic properties of the fiber suspensions was investigated. Empirical equations were obtained for the relative viscosity and first normal-stress difference for the fiber suspensions in glycerin. The capillary flow of these suspensions is discussed in part II.     

Correlation mechanism of friction behavior and topological properties of the contact network during powder compaction

The effect of friction behavior on the compacted density is significant, but the relationship between the topological properties of the contact network and friction behavior during powder compaction remains unclear. Based on the discrete element method (DEM), a DEM model for die compaction was established, and the Hertz contact model was modified into an elastoplastic contact model that was more suitable for metal-powder compaction. The evolution of the topological properties of the contact network and its mechanism during powder compaction was explored using the elastoplastic contact model. The results demonstrate that the friction behavior between the particles is closely related to the topological properties of the contact network. Side wall friction results in smaller clustering coefficient (CC) and excess contact (EC) in the lower region near the side wall. Corresponding to this phenomenon, the upper region near the side wall has more high-stress particles when the major principal stress threshold was considered, and the CC and EC are significantly higher than those in the other regions. This study provides a theoretical basis for improving powder compaction behavior.     

Comparative studies on the effect of duct height on heat transfer and flow behavior between co-angular and co-rotating type finned surface

This paper presents the comparative studies on the effect of duct height on heat transfer and flow behavior between co-angular and co-rotating type finned surface in duct. Experiments were performed to investigate the effect of duct height on heat transfer enhancement of a surface affixed with arrays (7 × 7) of short rectangular plate fins of a co-angular and a co-rotating type pattern in the duct. An infrared imaging system with the camera of TVS 8000 was used to measure the temperature distributions to calculate the local heat transfer coefficients of the representative fin regions. Pressure drop and heat transfer experiments were performed for both types of fin pattern varying the duct to fin height ratio (H _d/H _f) of 2.0–5.0. The friction factor calculated from the pressure drop shows that friction factor decreases with increasing the duct to fin height ratio (H _d/H _f) regardless of fin pattern and this is expected because the larger friction occurs for smaller duct to fin height ratios. Detailed heat transfer distribution gives a clear picture of heat transfer characteristics of the overall surface as well as the influence of the duct height. In addition, different flow behavior and flow structure developed by both patterns were visualized by the smoke flow visualization technique.     

The effect of entanglements on the rheological behavior of polybutadiene critical gels

We investigated the stress relaxation behavior of critical gels originating from six nearly monodisperse, highly entangled polybutadiene melts of different molecular weight from 18000 to 97 000 g/mole. The polymers were vulcanized by a hydrosilation reaction which takes place nearly exclusively at the pendant 1,2-vinyl sites distributed randomly along the polybutadiene chain. The BSW spectrum represents the relaxation of the initial uncrosslinked precursor. A characteristic parameter is the longest relaxation time of the precursor. Crosslinking increases this longest time even further. Surprisingly, the relaxation spectrum of the precursor is not altered much by the crosslinking except for an additional long time behavior. At the gel point (critical gel), this long time behavior is self-similar. It follows the typical power law as described by the Chambon-Winter gel equation, G(t) = St ^–n , in the terminal zone. The critical relaxation exponent was found to be close to n = 0.5 over a range of stoichiometric ratios and for all precursor molecular weights analyzed. A new scaling relationship was found between the gel stiffness, S, and the precursor molecular weight of the form: S M _w ^zn , where exponent z from the zero shear viscosity-molecular weight relationship, ₀ M _w ^z , is commonly found to be z = 3.3 – 3.6.     

The effect of temperature on the viscoelastic behavior of linear low-density polyethylene

Summary Observations are reported on linear low-density polyethylene in isothermal torsional oscillation and relaxation tests at various temperatures ranging from room temperature to 110 C. Constitutive equations are derived for the viscoelastic response of a semicrystalline polymer at small strains. The polymer is treated as an equivalent network of strands bridged by junctions (entanglements, physical cross-links on the surfaces of crystallites and lamellar blocks). The network is thought of as an ensemble of meso-regions with various potential energies for rearrangement of strands. Two types of meso-domains are introduced: active, where strands separate from temporary junctions as they are excited by thermal fluctuations, and passive, where detachment of strands is prevented by the surrounding macromolecules. The time-dependent behavior of the ensemble reflects separation of active strands from their junctions and merging of dangling strands with the network. Stress–strain relations are developed by using the laws of thermodynamics. The governing equations involve six material constants that are found by fitting the experimental data. The study focuses on the effects of (i) temperature, (ii) the deformation mode (torsion versus bending), and (iii) the loading program (oscillations versus relaxation) on the adjustable parameters.This work was partially supported by the West Virginia Research Challenge Grant Program     

The effect of particle shape and distribution on the macroscopic behavior of magnetoelastic composites

The magnetoelastic homogenization framework and the partial decoupling approximation proposed by Ponte Castañeda and Galipeau (2011) are used to estimate material properties for a class of magnetically susceptible elastomers. Specifically, we consider composites consisting of aligned, ellipsoidal magnetic particles distributed randomly with “ellipsoidal” symmetry under combined magnetic and mechanical loading. The model captures the coupling between the magnetic and mechanical fields, including the effects of magnetic saturation. The results help elucidate the effects of particle shape, distribution, and concentration on properties such as the magnetostriction, actuation stress, magnetic modulus, and magnetization behavior of a magnetorheological composite.     

The effect of die flow on the dynamics of isothermal melt spinning

The effect of die flow variables on the stability of isothermal melt spinning has been studied, both theoretically and experimentally. A die flow analysis provides the boundary conditions for a differential treatment of the spinline, both as a steady flow problem and as a linear stability problem. From the latter, one can predict the onset of draw resonance as a function of draw ratio, certain rheological parameters, and the stresses in the die. The experimental materials were two commercial polypropylenes and the apparatus consisted of a short (1.5–6 cm) isothermal spinning chamber; the agreement with theory was quite satisfactory. In most cases, high shear rates in the die (and subsequent high die swells) decrease the spinline stability but the magnitude of this interaction depends on many variables. In general, there is a high propensity for draw resonance (or ductile breakage) when the spinline is operating under conditions of severe thinning in a rheological sense.     

The void-size effect on plastic flow localization in the Gurson model

Recent studies have shown that the size of microvoids has a significant effect on the void growth rate. The purpose of this paper is to explore whether the void size effect can influence the plastic flow localization in ductile materials. We have used the extended Gurson‘s dilatational plasticity theory, which accounts for the void size effect, to study the plastic flow localization in porous solids with long cylindrical voids. The localization model of Rice is adopted, in which the material inside the band may display a different response from that outside the band at the incipient plastic flow localization. The present study shows that it has little effect on the shear band angle.     

The behavior of naturally fractured reservoirs including fluid flow in matrix blocks

An alternative analytical solution to the system of equations for fluid flow through a double-porosity medium with a boundary condition of an equipotential surface is given in this paper. The problem is reduced to solving an integral equation. The solution is straight-forward, and involves only ordinary Bessel functions. Numerical results show that D, the ratio of matrix system permeability to fracture system permeability, has a strong effect on the two semilog straight lines characteristic of the pressure response in a double-porosity medium. As D increases from zero (the Warren-Root model) to one, the first semilog straight line moves closer to the second. This is similar to the effect of increasing ω, the ratio of storage capacity of the fracture system to total storage capacity, in the Warren-Root model which neglects the flow within matrix blocks.     

The effect of orifices on the liquid distribution in annular two-phase flow

The effect of using orifices to disrupt the water film in air-water annular two-phase flow has been studied experimentally in a vertical tube by measuring the wall film flowrate at various distances upstream and downstream of several different sizes of orifice. The orifices cause a temporary reduction in the downstream water film flowrate, which returns to its equilibrium value further downstream. The experimental results have been used, together with those of other investigators, to compare the effects of orifices to those of swirl tapes, and further to compare the processes of entrainment and deposition within annular two-phase flow.     

Shear and elongational flow behavior of acrylic thickener solutions

We investigate the effect of hydrophobic aggregation in alkali-swellable acrylic thickener solutions on shear and extensional flow properties at technically relevant polymer concentrations using the commercial thickener Sterocoll FD as model system. Apparent molecular weight of aggregates in water is M _w  ≈ 10⁸ g/mol and decreases by more than an order of magnitude in ethanol. Zero shear viscosity η ₀ is low and shear thinning is weak compared to the high molecular weight of the aggregates. Linear viscoelastic relaxation is described by the Zimm theory up to frequencies of 10⁴ rad/s, demonstrating that no entanglements are present in these solutions. This is further supported by the concentration dependence of η ₀ and is attributed to strong association within the aggregates. Extensional flow behavior is characterized using the capillary break-up extensional rheometry technique including high-speed imaging. Solutions with ϕ ≥ 1% undergo uniform deformation and show pronounced strain hardening up to large Hencky strains. Elongational relaxation times are more than one order of magnitude lower than the longest shear relaxation times, suggesting that aggregates cannot withstand strong flows and do not contribute to the elongational viscosity.