Integral-based constitutive equation for rubber at high strain rates

High-speed experiments were conducted to characterize the deformation and failure of Styrene Butadiene Rubber at impact rates. Dynamic tensile stress–strain curves of uniaxial strip specimens and force–extension curves of thin sheets were obtained from a Charpy tensile impact apparatus. Results from the uniaxial tension tests indicated that although the rubber became stiffer with increasing strain rates, the stress–strain curves remained virtually the same above 280 s⁻¹. Above this critical strain rate, strength, fracture strain and toughness decreased with increasing strain rates. When strain rates were below 180 s⁻¹, the initial modulus, tensile strength and breaking extension increased as the strain rate increased. Between strain rates of 180 and 280 s⁻¹, the initial modulus and tensile strength increased with increasing strain rates but the extension at break decreased with increasing strain rates. A hyper-viscoelastic constitutive relation of integral form was used to describe the rate-dependent material behavior of the rubber. Two characteristic relaxation times, 5 ms and 0.25 ms, were needed to fit the proposed constitutive equation to the data. The proposed constitutive equation was implemented in ABAQUS Explicit via a user-defined subroutine and used to predict the dynamic response of the rubber sheets in the experiments. Numerical predictions for the transient deformation and failure of the rubber sheet were within 10% of experimental results.     

Mechanics of the rate-dependent elastic–plastic deformation of glassy polymers from low to high strain rates

A combined experimental and analytical investigation has been performed to understand the mechanical behavior of two amorphous polymers—polycarbonate and poly(methyl methacrylate)—at strain rates ranging from 10⁻⁴ to 10⁴ s⁻¹. This range in strain rates was achieved in uniaxial tension and compression tests using a dynamic mechanical analyzer (DMA), a servo-hydraulic testing machine, and an aluminum split-Hopkinson pressure bar. DMA tension tests were used to characterize the viscoelastic behavior of these materials, with focus on the rate-dependent shift of material transition temperatures. Uniaxial compression tests on the servo-hydraulic machine (10⁻⁴ to 1 s⁻¹) and the split-Hopkinson pressure bar (10³ to 10⁴ s⁻¹) were used to characterize the rate-dependent yield and post-yield behavior. Both materials were observed to exhibit increased rate sensitivity of yield under the same strain rate/temperature conditions as the β-transition of the viscoelastic behavior. A physically based constitutive model for large strain deformation of thermoplastics was then extended to encompass high-rate conditions. The model accounts for the contributions of different molecular motions which become operational and important in different frequency regimes. The new features enable the model to not only capture the transition in the yield behavior, but also accurately predict the post-yield, large strain behavior over a wide range of temperatures and strain rates.     

Influence of temperature and strain rate on the mechanical behavior of three amorphous polymers: Characterization and modeling of the compressive yield stress

Uniaxial compression stress–strain tests were carried out on three commercial amorphous polymers: polycarbonate (PC), polymethylmethacrylate (PMMA), and polyamideimide (PAI). The experiments were conducted under a wide range of temperatures (−40 °C to 180 °C) and strain rates (0.0001 s⁻¹ up to 5000 s⁻¹). A modified split-Hopkinson pressure bar was used for high strain rate tests. Temperature and strain rate greatly influence the mechanical response of the three polymers. In particular, the yield stress is found to increase with decreasing temperature and with increasing strain rate. The experimental data for the compressive yield stress were modeled for a wide range of strain rates and temperatures according to a new formulation of the cooperative model based on a strain rate/temperature superposition principle. The modeling results of the cooperative model provide evidence on the secondary transition by linking the yield behavior to the energy associated to the β mechanical loss peak. The effect of hydrostatic pressure is also addressed from a modeling perspective.     

Thermo-mechanical response of nylon 101 under uniaxial and multi-axial loadings: Part I,Experimental results over wide ranges of temperatures and strain rates

A comprehensive study of the thermo-mechanical response of a thermoplastic polymer, nylon 101 is presented. Quasi-static and dynamic compression uniaxial and multi-axial experiments (stress states) were performed at a wide range of strain rates (10⁻⁵ to 5000 s⁻¹) and temperatures (−60 to 177 °C or −76 to 350 °F). The material is found to be non-linearly dependent on strain rate and temperature. The change in volume after plastic deformation is investigated and is found to be negligibly small. The relaxation and creep responses at room temperature are found to be dependent on strain rate and the stress–strain level at which these phenomena are initiated. Total deformation is decomposed into visco-elastic and visco-plastic components; these components have been determined at different levels of deformation. Results from non-proportional uniaxial to biaxial compression, and torsion experiments, are also reported for three different strain rates at room temperature. It is shown that nylon 101 has a response dependent on the hydrostatic pressure.     

Irreversible thermodynamics of triple junctions during the intergranular void motion under the electromigration forces

A rigorous reformulation of internal entropy production and the rate of entropy flow is developed for multi-component systems consisting of heterophases, interfaces and/or surfaces. The result is a well-posed moving boundary value problem describing the dynamics of curved interfaces and surfaces associated with voids and/or cracks that are intersected by grain boundaries. Extensive computer simulations are performed for void configuration evolution during intergranular motion. In particular we simulate evolution resulting from the action of capillary and electromigration forces in thin film metallic interconnects having a “bamboo” structure, characterized by grain boundaries aligned perpendicular to the free surface of the metallic film interconnects. Analysis of experimental data utilizing previously derived mean time to failure formulas gives consistent values for interface diffusion coefficients and enthalpies of voids. 3.0 × 10⁻⁶ exp(−0.62 eV/kT) m² s⁻¹ is the value obtained for voids that form in the interior of the aluminum interconnects without surface contamination. 6.5 × 10⁻⁶ exp(−0.84 eV/kT) m² s⁻¹ is obtained for those voids that nucleate either at triple junctions or at the grain boundary-technical surface intersections, where the chemical impurities may act as trap centers for hopping vacancies.     

The characteristics of plastic flow and a physically-based model for 3003 Al–Mn alloy upon a wide range of strain rates and temperatures

The uniaxial compressive responses of 3003 Al–Mn alloy upon strain rates ranging from 0.001/s to about 10⁴/s with initial temperatures from 77 K to 800 K were investigated. Instron servohydraulic testing machine and enhanced split Hopkinson bar facilities have been employed in such uniaxial compressive loading tests. The maximum true strain up to 80% has been achieved. The following observations have been obtained from the experimental results: 1) 3003 Al–Mn alloy presents remarkable ductility and plasticity at low temperatures and high strain rates; 2) its plastic flow stress strongly depends on the applied temperatures and strain rates; 3) the temperature history during deformation strongly affects the microstructure evolution within the material. Finally, paralleled with the systematic experimental investigations, a physically-based model was developed based on the mechanism of dislocation kinetics. The model predictions are compared with the experimental results, and a good agreement has been observed.     

Plastic anisotropy and the role of non-basal slip in magnesium alloy AZ31B

Mechanistic explanations for the plastic behavior of a wrought magnesium alloy are developed using a combination of experimental and simulation techniques. Parameters affecting the practical sheet formability, such as strain hardening rate, strain rate sensitivity, the degree of anisotropy, and the stresses and strains at fracture, are examined systematically by conducting tensile tests of variously oriented samples at a range of temperatures (room temperature to 250 °C) and strain rates (10⁻⁵–0.1 s⁻¹). Polycrystal plasticity simulations are used to model the observed anisotropy and texture evolution. Strong in-plane anisotropy observed at low temperatures is attributed to the initial texture and the greater than anticipated non-basal cross-slip of dislocations with 〈a〉 type Burgers vectors. The agreement between the measured and simulated anisotropy and texture is further validated by direct observations of the dislocation microstructures using transmission electron microscopy. The increase in the ductility with temperature is accompanied by a decrease in the flow stress, an increase in the strain rate sensitivity, and a decrease in the normal anisotropy. Polycrystal simulations indicate that an increased activity of non-basal, 〈c + a〉, dislocations provides a self-consistent explanation for the observed changes in the anisotropy with increasing temperature.     

Flow patterns and pressure drop of ionic liquid–water two-phase flows in microchannels

The two-phase flow of a hydrophobic ionic liquid and water was studied in capillaries made of three different materials (two types of Teflon, FEP and Tefzel, and glass) with sizes between 200 μm and 270 μm. The ionic liquid was 1-butyl-3-methylimidazolium bis{(trifluoromethyl)sulfonyl}amide, with density and viscosity of 1420 kg m⁻³ and 0.041 kg m⁻¹ s⁻¹, respectively. Flow patterns and pressure drop were measured for two inlet configurations (T- and Y-junction), for total flow rates of 0.065–214.9 cm³ h⁻¹ and ionic liquid volume fractions from 0.05 to 0.8. The continuous phase in the glass capillary depended on the fluid that initially filled the channel. When water was introduced first, it became the continuous phase with the ionic liquid forming plugs or a mixture of plugs and drops within it. In the Teflon microchannels, the order that fluids were introduced did not affect the results and the ionic liquid was always the continuous phase. The main patterns observed were annular, plug, and drop flow. Pressure drop in the Teflon microchannels at a constant ionic liquid flow rate, was found to increase as the ionic liquid volume fraction decreased, and was always higher than the single phase ionic liquid value at the same flow rate as in the two-phase mixture. However, in the glass microchannel during plug flow with water as the continuous phase, pressure drop for a constant ionic liquid flow rate was always lower than the single phase ionic liquid value. A modified plug flow pressure drop model using a correlation for film thickness derived for the current fluids pair showed very good agreement with the experimental data.     

On the interplay between material flaws and dynamic necking

In this paper we investigate the interplay between material defects and flow localization in elastoplastic bars subjected to dynamic tension. For that task, we have developed a 1D finite difference scheme within a large deformation framework in which the material is modelled using rate-dependent J₂ plasticity. A perturbation of the initial yield stress is introduced in each node of the finite difference mesh to model localized material flaws. Numerical computations are carried out within a wide spectrum of strain rates ranging from 500 s⁻¹ to 2500 s⁻¹. On the one hand, our calculations reveal the effect of the material defects in the necking process. On the other hand, our results show that the necking inception, instead of being a random type process, is the deterministic result of the interplay between the mechanical behaviour of the material and the boundary conditions. This conclusion agrees with the experimental evidence reported by Rittel et al. [1] and Rotbaum et al. [2].     

Numerical investigation of particle deposition inside aero-shielded solar cyclone reactor: A promising solution for reactor clogging

Solar cracking of methane is considered to be an attractive option due to its CO₂ free hydrogen production process. Carbon particle deposition on the reactor window, walls and exit is a major obstacle to achieve continuous operation of methane cracking solar reactors. As a solution to this problem a novel “aero-shielded solar cyclone reactor” was created. In this present study the prediction of particle deposition at various locations for the aero-shielded reactor is numerically investigated by a Lagrangian particle dispersion model. A detailed three dimensional computational fluid dynamic (CFD) analysis for carbon deposition at the reactor window, walls and exit is presented using a Discrete Phase Model (DPM). The flow field is based on a RNG k–ε model and species transport with methane as the main flow and argon/ hydrogen as window and wall screening fluid. Flow behavior and particle deposition have been observed with the variation of main flow rates from 10–20 L/min and with carbon particle mass flow rate of 7 × 10⁻⁶ and 1.75 × 10⁻⁵ kg/s. In this study the window and wall screening flow rates have been considered to be 1 L/min and 10 L/min by employing either argon or hydrogen. Also, to study the effect of particle size simulations have also been carried out (i) with a variation of particle diameter with a size distribution of 0.5–234 μm and (ii) by taking 40 μm mono sized particles which is the mean value for the considered size distribution. Results show that by appropriately selecting the above parameters, the concept of the aero-shielded reactor can be an attractive option to resolve the problem of carbon deposition at the window, walls and exit of the reactor.     

Ductility of selected metals under electromagnetic ring test loading conditions

Experimental studies on ductility of selected metals differing mechanical properties under strain rates between 4 × 10³ and 2 × 10⁴ s^?1 are presented in this work. The electromagnetic expanding ring experiment was used as the primary tool for examining the ductility behaviour of metals. Through a use of the Phantom v12 digital high-speed camera and specialised TEMA Automotive software, rings expansion velocities were determined with satisfactory accuracy for all ring tests. In this paper, the experimental observations on cold-rolled copper Cu-ETP, aluminium alloy Al 7075, barrel steel and tungsten heavy alloy are reported. Ductility of studied materials was estimated by measuring changes in cross-sectional areas in the uniform strain portions of the recovered ring fragments. In a similar way the metals ductility was defined at the conventional tensile test condition. Moreover, results of analogue investigation for static and dynamic loading at the temperature of about ?40 °C were described. The experimental observations mainly revealed the different ductility behaviour of metals tested at applied dynamic loadings; Cu-ETP and barrel steel demonstrated an increase in ductility, whereas aluminium alloy Al 7075 and tungsten heavy alloy were characterised by lower ductility in comparison to static loading. These results appear to be partially in contrast with the observations reported recently by some other investigators.     

Effect of bubble size on void fraction fluctuations in dispersed bubble flows

It is known that bubble size affects seriously the average void fraction in bubbly flows where buoyant velocities vary considerably with bubble size. On the contrary, there is no systematic literature report about bubble size effects on the intensity and frequency of void fraction fluctuations around the average void fraction. This work aims to provide such information. An electrical impedance technique is employed along with non-intrusive ring electrodes to register void fraction fluctuations down to 10⁻⁵. Bubble size fluctuations are estimated from high resolution optical images. Experiments are conducted in co-current upward dispersed bubble flow inside a 21 mm tube with average bubble size between ∼50 and ∼700 μm. Water and blood simulant are used as test liquids with velocity from ∼3 to ∼30 cm s⁻¹. The above resemble conditions of Decompression Sickness (DCS) in the bloodstream of human vena cava. It is found that the intensity and frequency of void fraction fluctuations vary appreciably with bubble size at constant gas and liquid flow rates. Moreover, these variations are not random but scale with bubble size. As a first step to quantify this effect, an empirical expression is derived that relates average bubble size to the ratio standard deviation/average value of void fraction.     

Influences of ions and temperature on performance of carbon nano-particulates in supercapacitors with neutral aqueous electrolytes

A commercial product of carbon nano-particles, Cabot MONACH 1300 pigment black (CMPB), was studied for basic structural information and electrochemical performance in neutral aqueous electrolytes, aiming at applications in supercapacitors. As confirmed by SEM and HRTEM, the CMPB had a hierarchical structure, containing basic 10 nm nano-spheres which combined into ca. 50 nm agglomerates which further aggregated into larger particles of micrometres. The capacitance of this commercial material was found to increase with decreasing the size of hydrous cation (Li⁺ → Na⁺ → K⁺), instead of the cation crystal radius (K⁺ → Na⁺ → Li⁺) when coupled with the same anion (Cl⁻). In electrolytes with the same cation concentration (K⁺), changing the anion from the larger dianion (SO₄²⁻) to the smaller monoanion (Cl⁻) also increased the capacitance at high potential scan rates (>50 mV/s). Increasing electrolyte concentration produced expected effect, including raising the electrode capacitance, but lowering the equivalent series resistance (ESR), charge transfer resistance (CTR), and the diffusion resistance. At higher temperatures, the CMPB exhibited slightly higher capacitance, which does not agree with the Gouy–Chapman theory on electric double layer (EDL). A hypothesis is proposed to account for the capacitance increase with temperature as a result of the CMPB opening up some micro-pores for more ions to access in response to the temperature increase.     

Real-time measurements of PM2.5, PM10–2.5, and BC in an urban street canyon

A continuous dichotomous beta gauge monitor was used to characterize the hourly content of PM_2.5, PM_10–2.5, and Black Carbon (BC) over a 12-month period in an urban street canyon of Hong Kong. Hourly vehicle counts for nine vehicle classes and meteorological data were also recorded. The average weekly cycles of PM_2.5, PM_10–2.5, and BC suggested that all species are related to traffic, with high concentrations on workdays and low concentrations over the weekends. PM_2.5 exhibited two comparable concentrations at 10:00–11:00 (63.4 μg/m³) and 17:00–18:00 (65.0 μg/m³) local time (LT) during workdays, corresponding to the hours when the numbers of diesel-fueled and gasoline-fueled vehicles were at their maximum levels: 3179 and 2907 h⁻¹, respectively. BC is emitted mainly by diesel-fueled vehicles and this showed the highest concentration (31.2 μg/m³) during the midday period (10:00–11:00 LT) on workdays. A poor correlation was found between PM_2.5 concentration and wind speed (R = 0.51, P-value > 0.001). In contrast, the concentration of PM_10–2.5 was found to depend upon wind speed and it increased with obvious statistical significance as wind speed increased (R = 0.98, P-value < 0.0001).     

The contribution of small leaks in a baghouse filter to dust emission in the PM2.5 range—A system approach

The contribution of leakage in a baghouse filter (defined as a short circuit between the upstream and downstream sides of the filter) to the emission of fine particles is quantified in comparison to other dust emission sources, and the influence of key operating variables on overall system response is analyzed. The study was conducted on a well-maintained pilot-scale filter unit (9 bags of 500 g/m² calendered polyester needle felt; total surface area 4.2 m²) operated in Δp-controlled mode over a range of pulsing intensities, with two types of test dust (one free-flowing and the other cohesive) at inlet concentrations of 10 and 30 g/m³. Leaks included single holes between 0.5 and 4 mm diameter, intentionally placed in either the plenum plate or one of the filter bags, as well as seamlines from bag confectioning. Emissions were separated by source into a transient contribution due to dust penetration through the filter bags after each cleaning pulse, and a continuous contribution from leaks. This separation was based on a novel method of data processing that relies on time-resolved concentration measurements with a specially calibrated optical particle counter. Tiny leaks on the order of 1 mm generated the same emission level as all the bags combined, and dominated continuous emissions. The equivalent leak cross section (leakage = media emission) was about 1 ppm of the total installed filter surface, independent of upstream dust concentration. Leakage through open seamlines amounted to 75% of media emissions in case of free-flowing test dust. Leakage was restricted to aerodynamic diameters less than ∼5 μm (roughly the PM_2.5 mass fraction). For comparison, time-averaged mass penetration through conventional needle-felt media ranged from about 10⁻⁵ to 10⁻⁶, depending on cohesiveness of the particle material and pulse cleaning intensity, giving emission levels between about 0.02 and 0.2 mg/m³ at the reference concentration of 10 g/m².     

Emission of organic carbon,elemental carbon and water-soluble ions from crop straw burning under flaming and smoldering conditions

Emissions from major agricultural residues were measured using a self-designed combustion system. Emission factors (EFs) of organic carbon (OC), elemental carbon (EC), and water-soluble ions (WSIs) (K⁺, NH₄⁺, Na⁺, Mg²⁺, Ca²⁺, Cl⁻, NO₃⁻, SO₄^2–) in smoke from wheat and rice straw were measured under flaming and smoldering conditions. The OC₁/TC (total carbon) was highest (45.8% flaming, 57.7% smoldering) among carbon fractions. The mean EFs for OC (EF_OC) and EC (EF_EC) were 9.2 ± 3.9 and 2.2 ± 0.7 g/kg for wheat straw and 6.4 ± 1.9 and 1.1 ± 0.3 g/kg for rice straw under flaming conditions, while they were 40.8 ± 5.6 and 5.8 ± 1.0 g/kg and 37.6 ± 6.3 and 5.0 ± 1.4 g/kg under smoldering conditions, respectively. Higher EC ratios were observed in particulate matter (PM) mass under flaming conditions. The OC and EC for the two combustion patterns were significantly correlated (p < 0.01, R = 0.95 for wheat straw; p < 0.01, R = 0.97 for rice straw), and a higher positive correlation between OC₃ and EC was observed under both combustion conditions. WSIs emitted from flaming smoke were dominated by Cl⁻ and K⁺, which contributed 3.4% and 2.4% of the PM mass for rice straw and 2.2% and 1.0% for wheat straw, respectively. The EFs of Cl⁻ and K⁺ were 0.73 ± 0.16 and 0.51 ± 0.14 g/kg for wheat straw and 0.25 ± 0.15 and 0.12 ± 0.05 g/kg for rice straw under flaming conditions, while they were 0.42 ± 0.28 and 0.12 ± 0.06 g/kg and 0.30 ± 0.27 and 0.05 ± 0.03 g/kg under smoldering conditions, respectively. Na⁺, Mg²⁺, and NH₄⁺ were vital components in PM, comprising from 0.8% (smoldering) to 3.1% (flaming) of the mass. Strong correlations of Cl⁻ with K⁺, NH₄⁺, and Na⁺ ions were observed in rice straw and the calculated diagnostic ratios of OC/EC, K⁺/Na⁺ and Cl⁻/Na⁺ could be useful to distinguishing crop straw burning from other sources of atmospheric pollution.     

An experimental and numerical study of the localization behavior of tantalum and stainless steel

In general, the shear localization process involves initiation and growth. Initiation is expected to be a stochastic process in material space where anisotropy in the elastic–plastic behavior of single crystals and inter-crystalline interactions serve to form natural perturbations to the material’s local stability. A hat-shaped sample geometry was used to study shear localization growth. It is an axi-symmetric sample with an upper “hat” portion and a lower “brim” portion with the shear zone located between the hat and brim. The shear zone length is 870–890 μm with deformation imposed through a split-Hopkinson pressure bar system at maximum top-to-bottom velocity in the range of 8–25 m/s. We present experimental results of the deformation response of tantalum and 316L stainless steel samples. The tantalum samples did not form shear bands but the stainless steel sample formed a late stage shear band. We have also modeled these experiments using both conductive and adiabatic continuum models. An anisotropic elasto-viscoplastic constitutive model with damage evolution was used within the finite element code EPIC. A Mie-Gruneisen equation of state and the rate and temperature sensitive MTS flow stress model together with a Gurson flow surface were employed. The models performed well in predicting the experimental data. The numerical results for tantalum suggested a maximum equivalent strain rate on the order of 7 × 10⁴ s⁻¹ in the gage section for an imposed top surface displacement rate of 17.5 m/s. The models also suggested that for an initial temperature of 298 K a temperature in the neighborhood of 900 K was reached within the shear section. The numerical results for stainless steel suggest that melting temperature was reached throughout the shear band shortly after peak load. Due to sample geometry, the stress state in the shear zone was not pure shear; a significant normal stress relative to the shear zone basis line was developed.     

Prediction of symmetry during intermittent and annular horizontal two-phase flows

An optical method is developed for a 2.95 mm inner diameter circular mini-channel to estimate liquid film thicknesses. The greyscale pictures obtained with a high-speed camera are processed to determine the liquid-vapour interface positions for annular and intermittent flows. The experiments are thus performed for a large range of flow conditions. The saturation temperatures tested ranged from 20°C to 100°C in steps of 10°C and the mass velocities are 50, 100, 200, 300 and 400 kg m⁻² s⁻¹. A new parameter ranging from 0 to 1, the symmetry, is defined to account for the level of non uniformity of liquid distribution around the tube perimeter. New experimental data are presented (270 data points), that cover a range of symmetry parameter from 0.35 to 1.00. These data are merged with those available in the literature (406 data points with symmetry parameter from 0.71 to 1.00). A sensitivity analysis of the dimensionless numbers of major influence is run and a new correlation is proposed, that enables to predict over 90% of the data points in an error bandwidth of 10%. This correlation is proposed as criterion for asymmetry.     

Influence of sudden contractions on in situ volume fractions for oil–water flows in horizontal pipes

Oil–water two-phase flow experiments were conducted in horizontal ducts made of Plexiglas® to determine the in situ oil fraction (holdup) by means of the closing valves technique, using mineral oil (viscosity: 0.838 Pa s at 20 °C; density: 890 kg m⁻³) and tap water. The ducts present sudden contractions from 50 mm to 40 mm i.d. and from 50 mm to 30 mm i.d., with contraction ratios of 0.64 and 0.36, respectively. About 200–320 tests were performed by varying the flow rates of the phases. Flow patterns were investigated for both the up- and downstream pipe in order to assess whether relevant variations of the flow patterns across the sudden contraction take place. Data were then compared with predictions of a specific correlation for oil–water flow and some correlations for gas–water flow. A drift-flux model was also applied to determine the distribution parameter.