Stress and fatigue life analyses of a five-piece rim and the proposed optimization with a two-piece rim

A five-piece rim and a two-piece bolt-connected rim were investigated to examine stress levels and fatigue lives on critical regions. The finite element models of the rim/tire assemblies were developed and validated through tire engineering data and previously validated modelling approaches. The rim/tire assemblies were simulated under two conditions, (1) application of a 23,100 kg static load followed by a 24.14 km/h travelling speed and an 82° wheel angle, and (2) application of a 26,900 kg static load followed by an 8.05 km/h travelling speed and an 82° wheel angle. The results revealed that travelling and steering speeds were the key factors in causing high stresses and bolt tension forces. Compared to the five-piece rim, the two-piece rim decreased the maximum stresses by over 30% for both loading conditions; consequently the fatigue lives were increased by over two orders of magnitude. The maximum bolt forces for the two-piece rim were estimated to be 195,680 N and 111,360 N separately.     

Squeeze flow of a power-law fluid between non-parallel plates

Squeeze flow in the gap between non-parallel circular plates of radius R is discussed. The test material is assumed to be a power-law fluid, with a no-slip boundary condition at the plates. If the mean separation between the plates is h, and the angle of inclination between the plates is ? ? h/R, the force on the plates is perturbed only at O(?²) and is increased by less than 10% if ? < 0.35h/R. A torque O(?) tends to return the plates to a parallel configuration.     

Experimental and computational studies of shock wave-to-bubbly water momentum transfer

Momentum transfer from shock waves (SWs) of various intensity (from 0.05 MPa to 0.5 MPa in amplitude) to water containing air bubbles 2.5 to 4 mm of mean diameter is studied both experimentally and by means of numerical simulation. Experiments are performed in a vertical shock tube of a 50 × 100 mm² rectangular cross section consisting of a 495-mm long high-pressure section (HPS), 495-mm long low-pressure section (LPS), and 990 mm long test section (TS) equipped with an air bubbler and filled with water. Experiments have shown that as the initial gas volume fraction in water increases from 0 to 0.3 the momentum imparted in bubbly water by SWs increases monotonically, gradually levelling off at an air volume fraction of about 0.30. The experimental data are confirmed by two-dimensional (2D) simulation of SW propagation in bubbly water in terms of the SW velocity versus the air content, pressure profiles, as well as liquid and gas velocity behind the shock front.     

CFD modelling and wind tunnel validation of airflow through plant canopies using 3D canopy architecture

The efficiency of pesticide application to agricultural fields and the resulting environmental contamination highly depend on atmospheric airflow. A computational fluid dynamics (CFD) modelling of airflow within plant canopies using 3D canopy architecture was developed to understand the effect of the canopy to airflow. The model average air velocity was validated using experimental results in a wind tunnel with two artificial model trees of 24 cm height. Mean air velocities and their root mean square (RMS) values were measured on a vertical plane upstream and downstream sides of the trees in the tunnel using 2D hotwire anemometer after imposing a uniform air velocity of 10 m s^?1 at the inlet. 3D virtual canopy geometries of the artificial trees were modelled and introduced into a computational fluid domain whereby airflow through the trees was simulated using Reynolds-Averaged Navier–Stokes (RANS) equations and k-ε turbulence model. There was good agreement of the average longitudinal velocity, U between the measurements and the simulation results with relative errors less than 2% for upstream and 8% for downstream sides of the trees. The accuracy of the model prediction for turbulence kinetic energy k and turbulence intensity I was acceptable within the tree height when using a roughness length (y₀ = 0.02 mm) for the surface roughness of the tree branches and by applying a source model in a porous sub-domain created around the trees. The approach was applied for full scale orchard trees in the atmospheric boundary layer (ABL) and was compared with previous approaches and works. The simulation in the ABL was made using two groups of full scale orchard trees; short (h = 3 m) with wider branching and long (h = 4 m) with narrow branching. This comparison showed good qualitative agreements on the vertical profiles of U with small local differences as expected due to the spatial disparities in tree architecture. This work was able to show airflow within and above the canopy in 3D in more details.     

CFD simulation of a gas–solid fluidized bed with two vertical jets

A computational fluid dynamics (CFD) model is used to investigate the hydrodynamics of a gas–solid fluidized bed with two vertical jets. Sand particles with a density of 2660 kg/m³ and a diameter of 5.0 × 10^?4 m are employed as the solid phase. Numerical computation is carried out in a 0.57 m × 1.00 m two-dimensional bed using a commercial CFD code, CFX 4.4, together with user-defined Fortran subroutines. The applicability of the CFD model is validated by predicting the bed pressure drop in a bubbling fluidized bed, and the jet detachment time and equivalent bubble diameter in a fluidized bed with a single jet. Subsequently, the model is used to explore the hydrodynamics of two vertical jets in a fluidized bed. The computational results reveal three flow patterns, isolated, merged and transitional jets, depending on the nozzle separation distance and jet gas velocity and influencing significantly the solid circulation pattern. The jet penetration depth is found to increase with increasing jet gas velocity, and can be predicted reasonably well by the correlations of Hong et al. (2003) for isolated jets and of Yang and Keairns (1979) for interacting jets.     

Conjugate turbulent heat transfer in a square cavity with a solar control coating deposited to a vertical semitransparent wall

This paper presents a numerical study of the conjugate heat transfer (natural convection, surface thermal radiation and conduction) in a square cavity with turbulent flow. The cavity has one vertical isothermal wall, two horizontal adiabatic walls and one vertical semitransparent wall with a selective coating applied to the inner side to control the solar radiation transmission. Later on the semitransparent wall is replaced with another one without the selective coating. The mathematical model for the turbulent flow in the cavity was solved using the finite volume method. The system had the following conditions: the uniform temperature in the isothermal wall was 21 °C, the external ambient temperature was fixed at 35 °C and on the semitransparent wall the direct normal solar irradiation of 750 W/m² was considered constant. The Rayleigh number was varied in the range of 10⁹ ? Ra ? 10¹² by changing the lengths of the cavity from 0.70 m to 6.98 m, respectively. The results show that, even though the air temperature of the cavity with the solar control film coating semitransparent wall (case A) is higher compared with the one without solar film coating (case B), the total amount of heat going through the cavity is lower compared to the one going through the cavity without solar control film. The total amount of energy transferred to the air in cavity for the case A was 41.98% less than for the case B. A set of correlations for the Nusselt number was obtained for both cases considering the conjugate heat transfer.     

Surface responses induced by point load or uniform traction moving steadily on an anisotropic half-plane

Surface responses induced by point load or uniform traction moving steadily with subsonic speed on an anisotropic half-plane boundary are investigated. It is found that the effects of the material constant on surface displacements are through matrices L^?1(v) and S(v)L^?1(v), while those on surface stress components are through matrices Ω(v) and Γ(v). Explicit expressions for the elements of these four matrices are expressed in terms of elastic stiffness for general anisotropic materials. The special cases of monoclinic materials with symmetry plane at x₁ = 0, x₂ = 0 and x₃ = 0, and the case for orthotropic materials are all deduced. Results for isotropic material may be recovered from present results. For monoclinic materials with a plane of symmetry at x₃ = 0, two of the elements of matrix Ω(v) are found to be independent of subsonic speed.     

Soil compaction and tires for harvesting and transporting sugarcane

Four tire types (A, block-shape tread; B, rib-shape tread; C, low-lug tread; D, high-lug tread) used to harvest and transport sugarcane were compared regarding the compaction induced to the soil. Tires were tested at three inflation pressures (207, 276, 345 kPa) and six loads ranging from 20 to 60 kN/tire. Track impressions were traced, and 576 areas were measured to find equations relating inflation pressure, load, contact surface and pressure. Contact surface increased with increasing load and decreasing inflation pressure; however, the contact pressure presented no defined pattern of variation, with tire types A and B generating lower contact pressure. The vertical stresses under the tires were measured and simulated with sensors and software developed at the Colombian Sugarcane Research Center (Cenicaña). Sensors were placed at 10, 30, 50 and 70 cm depth. Tire types A and B registered vertical stresses below 250 kPa at the surface. These two tires were better options to reduce soil compaction. The equations characterizing the tires were introduced into a program to simulate the vertical stress. Simulated and measured stresses were adjusted in an 87–92% range. Results indicate a good correlation between the tire equations, the vertical stress simulation and the vertical stress measurement.     

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).     

An investigation of a model of the flow pattern transition mechanism in relation to the identification of annular flow of R134a in a vertical tube using various void fraction models and flow regime maps

In the present study, new experimental data are presented for literature on the prediction of film thickness and identification of flow regime during the co-current downward condensation in a vertical smooth copper tube having an inner diameter of 8.1 mm and a length of 500 mm. R134a and water are used as working fluids in the tube side and annular side of a double tube heat exchanger, respectively. Condensation experiments are done at mass fluxes of 300 and 515 kg m^?2 s^?1. The condensing temperatures are between 40 and 50 °C; heat fluxes are between 12.65 and 66.61 kW m^?2. The average experimental heat transfer coefficient of the refrigerant HFC-134a is calculated by applying an energy balance based on the energy transferred from the test section. A mathematical model by Barnea et al. based on the momentum balance of liquid and vapor phases is used to determine the condensation film thickness of R134a. The comparative film thickness values are determined indirectly using relevant measured data together with various void fraction models and correlations reported in the open literature. The effects of heat flux, mass flux, and condensation temperature on the film thickness and condensation heat transfer coefficient are also discussed for the laminar and turbulent flow conditions. There is a good agreement between the film thickness results obtained from the theoretical model and those obtained from six of 35 void fraction models in the high mass flux region of R134a. In spite of their different valid conditions, six well-known flow regime maps from the literature are found to be predictive for the annular flow conditions in the test tube in spite of their different operating conditions.     

Some C0-continuous mixed formulations for general dipolar linear gradient elasticity boundary value problems and the associated energy theorems

The goal of this work is a systematic presentation of some classes of mixed weak formulations, for general multi-dimensional dipolar gradient elasticity (fourth order) boundary value problems. The displacement field main variable is accompanied by the double stress tensor and the Cauchy stress tensor (case 1 or μ ? τ ? u formulation), the double stress tensor alone (case 2 or μ ? u formulation), the double stress, the Cauchy stress, the displacement second gradient and the standard strain field (case 3 or μ ? τ ? κ ? ε ? u formulation) and the displacement first gradient, along with the equilibrium stress (case 4 or u ? θ ? γ formulation). In all formulations, the respective essential conditions are built in the structure of the solution spaces. For cases 1, 2 and 4, one-dimensional analogues are presented for the purpose of numerical comparison. Moreover, the standard Galerkin formulation is depicted. It is noted that the standard Galerkin weak form demands C¹-continuous conforming basis functions. On the other hand, up to first order derivatives appear in the bilinear forms of the current mixed formulations. Hence, standard C⁰-continuous conforming basis functions may be employed in the finite element approximations. The main purpose of this work is to provide a reference base for future numerical applications of this type of mixed methods. In all cases, the associated quadratic energy functionals are formed for the purpose of completeness.     

Differences in tractor performance parameters between single-wheel 4WD and dual-wheel 2WD driving systems

Vertical wheel load and tire pressure are both easily managed parameters which play a significant role in tillage operations for limiting slip which involves energy loss. This aspect to a great extent affects the fuel consumption and the time required for soil tillage. The main focus of this experiment was to determine the effect on the wheels’ slip, the fuel consumption and the field performance of a tractor running in a single-wheel 4WD driving system and in a dual-wheel 2WD driving system, due to the variations in air pressure of the tires as well as in the ballast mass. With no additional mass, the lowest fuel consumption was reached by a tractor with the least air pressure in the tires and running in a dual-wheel 2WD driving system. It was determined that for a stubble cultivation with a medium-power (82.3 kW) tractor running in a dual-wheel 2WD driving system, the hourly fuel consumption was by 1.15 L h⁻¹ (or 7.3%), the fuel consumption per hectare by 0.35 L ha⁻¹ (or 7.9%) and the field performance by 0.05 ha h⁻¹ (or 1.25%) lower compared to a single-wheel 4WD driving system, when driving wheels’ slip for both modes was the same, i.e., at 8–12%.     

Drag reduction for liquid flow through micro-apertures

Pressure drops in the flow through micro-orifices and capillaries were measured for silicone oils, aqueous solutions of polyethylene glycol (PEG), and surfactant aqueous solutions. The diameter of micro-orifices ranged from 5 μm to 400 μm. The corresponding length/diameter ratio was from 4 to 0.05 and capillary diameters were 105 μm and 450 μm. The following results were obtained: silicone oils of 10^?6 m²/s and 10^?5 m²/s in kinematic viscosity generated a reduction of pressure drop (RPD), that is, drag reduction, similar to the RPD of water and a glycerol/water mixture reported in the previous paper by the present authors. When RPD occurred, the pressure drop (PD) of silicone oils of 10^?6 m²/s and 10^?5 m²/s had nearly the same magnitude. Namely, the difference in viscosity did not influence RPD. A 10³ ppm aqueous solution of PEG20000 provided almost the same PD as that of PEG8000 for the 400 μm to 15 μm orifices, but a greater PD than that of PEG8000 for the 10 μm to 5 μm orifices. A non-ionic surfactant and a cationic surfactant were highly effective in RPD compared with anionic surfactants: the non-ionic and cationic surfactant solutions had PD one order of magnitude lower than that of water under some flow conditions in the concentration range from 1 ppm to 10⁴ ppm, but the anionic surfactant solutions did not generate RPD except in the case of the smallest orifice of 5 μm in diameter. The PD of the non-ionic surfactant solution showed a steep rise at a Reynolds number (Re_t) for 400 μm to 15 μm orifices. The Re_t provides the relationship Re_t = K/D, where D is the orifice diameter, and K is a constant of 2 × 10^?2 m for the 100–20 μm orifices irrespective of liquid concentration. Capillary flow experiment revealed that the PEG, non-ionic and cationic surfactant solutions generated RPD also in a laminar flow through the capillary of 105 μm in diameter, but not in the flow through the capillary of 450 μm in diameter. In order to clarify the cause of RPD, an additional experiment was carried out by changing the orifice material from metal to acrylic resin. The result gave a different appearance of RPD, suggesting that RPD is related to an interfacial phenomenon between the liquid and wall. The large RPDs found in the present experiment are very interesting from both academic and practical viewpoints.     

Chimney effect due to different vertical position of an isothermal horizontal cylinder confined between two adiabatic walls

The variation of natural convection heat transfer from an isothermal horizontal cylinder confined between two adiabatic walls of constant height is investigated by Mach-Zehnder interferometry technique. This paper focuses on the chimney effect due to the vertical position changes of cylinder (Y) located between two walls with a constant distance of W measuring 1.5 cylinder diameter. The cylinder’s local and average Nusselt numbers are determined for ratio of vertical position to its diameter ranging from Y/D = (0 to 10), and the Rayleigh number ranging from 3.5 × 10³ to 1.4 × 10⁴. There is an optimum distance between the walls in which the Nusselt number is maximum. Results are indicated with a single correlation which gives the average Nusselt number as a function of the ratio of vertical position to cylinder diameter and the Rayleigh number. The experimental data shows that there is an optimum vertical position for the cylinder at which the Nusselt number has a maximum value at each Rayleigh number. This optimal vertical position is derived from the correlation and is presented by an equation. The value of the optimum vertical position increases as the Rayleigh number increases.     

Plasma emission spectroscopy of solids irradiated by intense XUV pulses from a free electron laser

The FLASH XUV-free electron laser has been used to irradiate solid samples at intensities of the order 10¹⁶ W cm^?2 at a wavelength of 13.5 nm. The subsequent time integrated XUV emission was observed with a grating spectrometer. The electron temperature inferred from plasma line ratios was in the range 5–8 eV with electron density in the range 10²¹–10²² cm^?3. These results are consistent with the saturation of absorption through bleaching of the L-edge by intense photo-absorption reported in an earlier publication.     

Flow boiling of halogenated refrigerants at high saturation temperature in a horizontal smooth tube

Several correlations are available in the open literature for computing the heat transfer coefficient during flow boiling inside plain channels. With respect to halogenated refrigerants, these correlations are usually compared to data taken in a limited range of evaporation temperature and reduced pressure. More recently, the adoption of new refrigerants, such as high pressure HFCs and carbon dioxide, requires to largely extend the pressure range of application of such correlations. Besides, the design of evaporators for some heat pumping applications, where temperatures are set at higher values as compared to usual evaporating temperatures in air-conditioning equipment, requires proper validation of the computing methods.The present paper aims at comparing four well-known predicting models to a new database collected during flow boiling of HCFC (R22) and HFC refrigerants (R134a, R125 and R410A) in a horizontal 8 mm internal diameter tube. This database is characterized by saturation temperature ranging between 25 °C and 45 °C, reduced pressure spanning between 0.19 and 0.53. Mass velocity ranges between 200 and 600 kg m^?2 s^?1 and heat flux between 9 and 53 kW m^?2.Evaporating heat transfer coefficients of halogenated refrigerants at such high temperatures have not been reported in the open literature so far. The discussion of the results will enlighten some similarities with experimental trends presented in the literature for evaporating carbon dioxide.Two models tested here show good prediction capabilities of the present experimental data, but not for all the data sets in the same way. For the purpose of practical use, a simple modification of the correlation by Gungor and Winterton [1] is proposed, showing that this is able to catch the experimental trends of the present database with good agreement.     

LES of turbulent flow in a concentric annulus with rotating outer wall

Fully-developed turbulent flow in a concentric annulus, r₁/r₂ = 0.5, Re_h = 12,500, with the outer wall rotating at a range of rotation rates N = U_θ,wall/U_b from 0.5 up to 4 is studied by large-eddy simulations. The focus is on the effects of moderate to very high rotation rates on the mean flow, turbulence statistics and eddy structure. For N up to ∼2, an increase in the rotation rate dampens progressively the turbulence near the rotating outer wall, while affecting only mildly the inner-wall region. At higher rotation rates this trend is reversed: for N = 2.8 close to the inner wall turbulence is dramatically reduced while the outer wall region remains turbulent with discernible helical vortices as the dominant turbulent structure. The turbulence parameters and eddy structures differ significantly for N = 2 and 2.8. This switch is attributed to the centrifuged turbulence (generated near the inner wall) prevailing over the axial inertial force as well as over the counteracting laminarizing effects of the rotating outer wall. At still higher rotation, N = 4, the flow gets laminarized but with distinct spiralling vortices akin to the Taylor–Couette rolls found between the two counter-rotating cylinders without axial flow, which is the limiting case when N approaches to infinity. The ratio of the centrifugal to axial inertial forces, Ta/Re² ∝ N² (where Ta is the Taylor number) is considered as a possible criterion for defining the conditions for the above regime change.     

Instantaneous flow field above the free end of finite-height cylinders and prisms

The flow above the free ends of surface-mounted finite-height circular cylinders and square prisms was studied experimentally using particle image velocimetry (PIV). Cylinders and prisms with aspect ratios of AR = 9, 7, 5, and 3 were tested at a Reynolds number of Re = 4.2 × 10⁴. The bodies were mounted normal to a ground plane and were partially immersed in a turbulent zero-pressure-gradient boundary layer, where the boundary layer thickness relative to the body width was δ/D = 1.6. PIV measurements were made above the free ends of the bodies in a vertical plane aligned with the flow centreline. The present PIV results provide insight into the effects of aspect ratio and body shape on the instantaneous flow field. The recirculation zone under the separated shear layer is larger for the square prism of AR = 3 compared to the more slender prism of AR = 9. Also, for a square prism with low aspect ratio (AR = 3), the influence of the reverse flow over the free end surface becomes more significant compared to that for a higher aspect ratio (AR = 9). For the circular cylinder, a cross-stream vortex forms within the recirculation zone. As the aspect ratio of the cylinder decreases, the reattachment point of the separated flow on the free end surface moves closer to the trailing edge. For both the square prism and circular cylinder cases, the instantaneous velocity vector field and associated in-plane vorticity field revealed small-scale structures mostly generated by the separated shear layer.     

Topological sensitivity to diagonal member flips of two-layered statically determinate trusses under worst loading

The paper deals with the topological sensitivity of free, unsupported, statically determinate plane trusses whose horizontal and vertical members form two horizontal layers of square cells and two or more vertical layers. The topology of a truss is decomposed into a form vector – the placement of cells containing diagonal members – and a binary vector describing the slopes of the diagonals. The construction of complete form and slope spaces is provided for any number of vertical layers. Using exhaustive search, forms with minimum and maximum sensitivity to slope change are found for trusses with 2 × 2 through 2 × 8 layers under worst static load condition, represented by the lowest eigenvalue of the least-squares equilibrium matrix. Typical features of the least and most sensitive forms and associated loads and internal forces are shown. Changes of absolute and relative topological sensitivities with increasing number of vertical layers are discussed.     

Development of a dynamometric vehicle to assess the drawbar performance of high-powered agricultural tractors

Agricultural tractors are machines originally designed to mechanize agricultural tasks, especially tillage and pulling. A large part of research activities have been interested in optimizing tractor efficiency, in particular in terms of emissions and energy. In this frame, the OECD Tractor Code 2 sets out a drawbar test in specific controlled conditions with the aim of evaluating the power of the tractor available at the drawbar. The principal measurement chain relies on dynamometric vehicles (DV) that are instrumented vehicles specifically engineered to develop horizontal force at the drawbar of agricultural tractors. The CREA Laboratory of Treviglio, Italy, engineered a new dynamometric vehicle to test tractors with up to 200 kW at the drawbar (245 kW at the engine flywheel) and a maximum of 118 kN drawbar force. The chosen basis is a FIAT 6605 N truck (TM 69 6 × 6) which has been transformed into a hydrostatic vehicle driven by a hydraulic system and an auxiliary gearbox. The maximum drawbar force was verified up to 122 kN. The drawbar power verification (200 kW) was successfully carried. The final verification confirmed that the project is valid for the investigation and optimization of the parameters regarding the traction efficiency of agricultural tractors.