This study focuses on the improvement of the axial buckling capacity of elliptical composite cylinders through the use of a circumferentially-varying lamination sequence. The concept of varying the lamination sequence around the circumference is considered as a viable approach for off-setting the disadvantages of having the cylinder radius of curvature vary with circumferential position, the source of the reduced buckling capacity when compared to a circular cylinder with the same circumference. Post-buckling collapse behavior and material failure characteristics are also of interest. Two approaches to implementing a circumferential variation of lamination are examined. For the first approach the lamination sequence is varied in a stepwise fashion around the circumference. Specifically, each quadrant of the cylinder circumference is divided into three equal-length regions denoted as the crown, middle, and side regions. Eight different cylinders designs, whereby each region is constructed of either a quasi-isotropic or an axially-stiff laminate of equal thickness, are studied. Results are compared to the baseline case of an elliptical cylinder constructed entirely of a quasi-isotropic laminate. Since the thickness of the quasi-isotropic and axially-stiff laminates are the same, all cylinders weight the same and thus comparisons are meaningful. Improvements upwards of 18% in axial buckling capacity can be achieved with one particular stepwise design. The second approach considers laminations that vary circumferentially in a continuous fashion to mitigate the effects of the continuously-varying radius of curvature. The methodology for determining how to tailor the lamination sequence circumferentially is based on the analytical predictions of a simple buckling analysis for simply-supported circular cylinders. With this approach, axial buckling load improvements upwards of 30% are realized. Of all the cylinders considered, very few do not exhibit material failure upon collapse in the post-buckled state. Of those that do not, there is little, if any, improvement in bucking capacity. Results for the pre-buckling, buckling, post-buckling, and material failure are obtained from the finite-element code ABAQUS using both static and dynamic analyses. Studies with the code demonstrate that the results obtained are converged. / Ph. D.

Investigation of momentum and heat transfer between the fluid and solid phase is critical to the study of fluid-particle systems. Dense suspensions are characterized by the solid fraction (ratio of solid volume to total volume), the particle Reynolds number, and the shape of the particle. The behavior of non-spherical particles deviates considerably from spherical particle shapes which have been studied extensively in the literature. Momentum transfer, to first-order, is driven by drag forces experienced by the particles in suspension, followed by lift and lateral forces, and also through the transmission of fluid torque to the particles. The subject of this thesis is a family of prolate ellipsoidal particle geometries of aspect ratios (AR) 2.5, 5.0 and 10.0 at nominal solid fractions (φ) between 0.1 and 0.3, and suspensions of cylinders of AR=0.25. The nominal particle Reynolds number (Re) is varied between 10 to 200, representative of fluidized beds. Fluid forces and heat transfer coefficients are obtained numerically by Particle Resolved Simulations (PRS) using the Immersed Boundary Method (IBM). The method enables the calculation of the interstitial flow and pressure field surrounding each particle in suspension leading to the direct integration of fluid forces acting on each particle in the suspension. A substantial outcome of the research is the development of a new drag force correlation for random suspensions of prolate ellipsoids over the full range of geometries and conditioned studied. In many practical applications, especially as the deviation from the spherical shape increases, particles are not oriented randomly to the flow direction, resulting in suspensions which have a mean preferential orientation. It is shown that the mean suspension drag varies linearly with the orientation parameter, which varies from -2.0 for particles oriented parallel to the flow direction to 1.0 for particles normal to the flow direction. This result is significant as it allows easy calculation of drag force for suspension with any preferential orientation. The heat transfer coefficient or Nusselt number is investigated for prolate ellipsoid suspensions. Significantly, two methods of calculating the heat transfer coefficient in the literature are reconciled and it is established that one asymptotes to the other. It is also established that unlike the drag force, at low Reynolds number the suspension mean heat transfer coefficient is very sensitive to the spatial distribution of particles or local-to-particle solid fractions. For the same mean solid fraction, suspensions dominated by particle clusters or high local solid fractions can exhibit Nusselt numbers which are lower than the minimum Nusselt number imposed by pure conduction on a single particle in isolation. This results from the dominant effect of thermal wakes at low Reynolds numbers. As the Reynolds number increases, the effect of particle clusters on heat transfer becomes less consequential. For the 0.25 aspect ratio cylinder, it was found that while existing correlations under predicted the drag forces, a sinusoidal function F_(d,θ)=F_(d,θ=0°)+(F_(d,θ=90°)-F_(d,θ=0°) )sin⁡(θ) captured the variation of normalized drag with respect to inclination angle over the range 10≤Re≤300 and 0≤φ≤0.3. Further the mean ensemble drag followed F_d=F_(d,θ=0°)+1/2(F_(d,θ=90°)-F_(d,θ=0°)). It was shown that lift forces were between 20% to 80% of drag forces and could not be neglected in models of fluid-particle interaction forces. Comparing the pitching fluid torque to collision torque during an elastic collision showed that as the particle equivalent diameter, density, and collision velocities decreased, fluid torque could be of the same order of magnitude as collisional torque and it too could not be neglected from models of particle transport in suspensions. / Doctor of Philosophy / Momentum and heat exchange between the fluids (air, water…) and suspensions of solid particles plays a critical role in power generation, chemical processing plants, pharmaceuticals, in the environment, and many other applications. One of the key components in momentum exchange are the forces felt by the particles in the suspension due to the flow of the fluid around them and the amount of heat the fluid can transfer to or from the particles. The fluid forces and heat transfer depend on many factors, chief among them being the properties of the fluid (density, viscosity, thermal properties) and the properties of the particles in the suspension (size, shape, density, thermal properties, concentration). This introduces a wide range of parameters that have the potential to affect the way the fluid and particles behave and move. Experimental measurements are very difficult and expensive to conduct in these systems and computational modeling can play a key role in characterization. For accuracy, computational models have to have the correct physical laws encoded in the software. The objective of this thesis is to use very high-fidelity computer models to characterize the forces and heat transfer under different conditions to develop general formulas or correlations which can then be used in less expensive computer models. Three basic particle shapes are considered in this study, a sphere, a disk like cylindrical particles, and particles of ellipsoidal shapes. More specifically, Particle Resolved Simulations of flow through suspensions of ellipsoids with aspect ratio of 2.5, 5, 10 and cylinders with aspect ratio of 0.25 are performed. The Reynolds number range covered is [10, 200] for ellipsoids and [10, 300] for cylinders with solid fraction range of [0.1, 0.3]. New fluid drag force correlations are proposed for the ellipsoid and cylinder suspensions, respectively, and heat transfer behavior is also investigated.

Prestressed Concrete Cylinder Pipe (PCCP) is a large-diameter and high-pressure conduit for drinking water and wastewater transmission. Due to its large diameter, high pressure, and mode of breakdown, PCCP failures usually have catastrophic consequences. To mitigate failures, it is very important to assess the condition of the pipe and take proactive measures, such as repair, rehabilitation, or replacement. There are many challenges in assessing the condition of PCCP. PCCP has a complex structure with several layers of materials (e.g. mortar coating, prestressing wire, steel cylinder, and concrete core) working together under loading. This means that there are many factors that can cause pipe failure, and that failure mechanisms are complicated. Data collection could be difficult, and existing data are often unavailable or unreliable. Considerable research has been conducted by scholars and engineers in developing models to evaluate the condition of PCCP. There are mainly two types of models: statistical models, and numerical models using finite element method. Statistical models consider only a few factors, such as pipe age and failure rate, to predict the failure of PCCP. However, the failure of PCCP can be caused by many other factors including pipe material, and loading conditions. Models only considering a few factors are not robust enough for reliable results. The current numerical models assume that all broken wires are centrally distributed in the same location and broken wires have no prestress, thus all broken wires are completely removed from the model. These assumptions could be overly conservative when actual broken wires are distributed in different locations along the pipeline and broken wires have remaining prestress due to the bond between the wire and mortar coating. Therefore, a more comprehensive numerical model is needed to have a better understanding of the condition of PCCP. In this research, an extensive literature and practice review was conducted on PCCP failures to understand the critical factors that affect pipe condition. The available technologies commonly used to detect pipe defects were reviewed in order to better understand the accuracy and uncertainties of the collected data. Existing models were reviewed to better understand their limitations and to advance the research on condition analysis of PCCP using numerical models. Based on these comprehensive reviews, this dissertation proposed a numerical model to analyze the condition of PCCP for its long-term performance management. Detailed structural components such as concrete cores, prestressing wires, steel cylinder, and mortar coating were modelled. The interactions between different layers of pipe components were considered. An algorithm was proposed to account for the bond between the prestressing wire and mortar coating, which is a critical factor for the condition of PCCP with broken wires. A FORTRAN program was developed to assign linear stress distribution between the broken point and the full-prestress resuming point. The proposed numerical model was verified utilizing data from lab tests and forensic study. Lab test data helped to understand the functionality of the model and to verify the model parameters used in analyzing pipe components and the simulation of interactions between different layers. The forensic data helped to verify the model under actual field working conditions of the pipe. Validation of the proposed numerical model was conducted using a 66-inch Embedded Cylinder Pipe and two Lined Cylinder Pipes (42-inch and 48-inch, respectively) from a water utility. In the validation, field data were collected for model development. The simulation results were consistent with the field observation, which proved the validity and applicability of the proposed numerical model in practice. A series of sensitivity studies were conducted to investigate the impact of longitudinal and circumferential location on the structural integrity of the pipe. These investigations showed that considering the actual longitudinal and circumferential location of broken wires is very important to get accurate analysis of pipe condition, while assuming that all broken wires fail in one longitudinal location (assumptions by current numerical models for PCCP) will overestimate the actual damage to the pipe caused by broken wires. To consider the bedding condition, a critical factor for PCCP, the four most common bedding types found in practice were analyzed. Results show that poor bedding could lead to cracks in PCCP, which could cause corrosion in prestressing wires. Therefore, it is very important to account for bedding conditions in the PCCP analysis. The model presented in this dissertation is more comprehensive and robust compared with existing numerical models, and could provide a better understanding of the condition of PCCP. This is because the proposed model considers the contribution of remaining prestress in broken wires due to the bond between the wire and mortar coating. This model can consider the actual longitudinal and circumferential location of broken wires rather than centrally distribute them, and it can consider the actual bedding locations, and the interaction between different layers of materials. This model was calibrated using lab test data and forensic data, and was further validated using field data which showed consistence between simulation results and field observations. The proposed model does have limitations due to limited availability of data and assumptions. Material tests were not conducted to verify the material properties used in the model, which could cause accuracy issues in the results. A full-scale simulation of the interaction between prestressing wire and mortar coating was not considered because it could significantly increase the computation time. Lab tests were not conducted to verify the parameters used for the simulation of interaction between concrete core and steel cylinder which could lead to accuracy problems. Finally, it is acknowledged that the model was only validated in one water utility and validations in more geographically distributed utilities might further test the model's validity and robustness. Nonetheless, the comprehensiveness and robustness of this proposed model improved the analysis of the condition of PCCP. The findings and results of this research will provide guidance for better management of PCCP pipelines for water utilities, and provide reference for future research on numerical modeling of PCCP as well. / Ph. D.

Two rotor configurations with different non-uniform inflows were studied: a rotor ingesting the wake of an upstream cylinder and a rotor ingesting a thick axially symmetric boundary layer from an upstream centerbody. In both cases, the undisturbed inflow was measured without the rotor present in order to characterize the inflow, in particular to calculate the unsteady upwash velocity distribution at the location of the rotor. In addition, detailed acoustic measurements were completed using a 251-channel large-area microphone array. In all, over 400 conditions covering different advance ratios, angles of yaw, and inflow conditions were measured. Measurements of the sound show that the source has a complex directivity, different from that of a streamwise aligned dipole, due to the inhomogeneous unsteady upwash distribution. In addition, observers at different far field locations will perceive sources from different locations on the rotor disk. The directivity is a function of both the rotor geometry and turbulent inflow. A simplified model of the sound source was developed using these inputs and accurately predicts trends observed in the far field noise. For the cylinder wake ingestion case, on-blade measurements of the flow field show that the wake is drawn to the center of the rotor disk with increasing thrust. This is particularly noticeable if the wake does not strike the center of the rotor disk. The effects of this flow distortion on the far field directivity are well predicted by the model. The effects of yaw to rotate the produced sound field can be inferred from this model as well. A novel beamforming procedure was used to isolate sources across the face of the rotor for the cylinder wake ingestion case for an upstream observer position. This method may be used to isolate different sound sources on a rotor if multiple sources are present or if different regions of the rotor disk need to be isolated. The directivity of a rotor ingesting an axially symmetric boundary layer is far less complex than the ingestion of a two-dimensional cylinder wake, but measurements still show the perceived source location shift with observer location. Overall, the proposed noise modeling technique is an efficient method to predict the directivity of turbulence ingestion noise for inhomogeneous inflows. This can enable quick absolute noise predictions at all far field locations using only a single point measurement or far field noise prediction to establish absolute levels. / Doctor of Philosophy / In many engineering applications, rotors interact with turbulence. Aircraft and ships with rear mounted propellers can have upstream appendages or discontinuities that generate turbulence that travels downstream and is drawn into the propeller. Wind turbines interact with turbulence in the atmosphere and with turbulent wakes from other turbines. Interaction of a rotor with turbulence results in unsteady loading on the rotor blades that can radiate as sound, causing unwanted community noise or vehicle detection. As such, prediction and reduction of noise due to turbulence ingestion is highly desirable and remains an active area of research. Turbulence ingestion noise is well understood from first principles and can be successfully predicted provided an accurate description of the turbulent inflow and unsteady aerodynamic response of the rotor blades. Much work has focused on homogenous, isotropic turbulence ingestion noise, however, in practical applications, the rotor inflow is often non-uniform, anisotropic, and can change dramatically with the thrusting condition of the rotor. Research efforts to develop noise predictions considering these more complex, but practical inflows have focused on the inflow modeling and measurement and have relied on a small subset of sound measurements for validation. The present study seeks to provide new physical insight into inhomogeneous, anisotropic turbulence ingestion noise through wind tunnel experiments. In particular, two rotor configurations with different practical non-uniform inflows are studied: a rotor ingesting the wake of an upstream cylinder and a rotor ingesting a thick axially symmetric boundary layer from an upstream center body. In both cases, the undisturbed inflow was measured without the rotor present in order to characterize the inflow, and detailed acoustic measurements were completed using a 251-channel large-area microphone array. In all, over 400 rotor operating conditions were measured. The acoustic directivity in each case is examined in detail as a function of rotor operating condition. A simplified directivity model is developed and validated with measurements. Ultimately, the directivity model can provide a good engineering approximation of the full directivity with reduced computational time or can be used to extrapolate measured results to positions in the far field where placement of sensors is not possible. The results can also be used to guide the analysis and interpretation of single point or microphone array measurements in the acoustic far field of a rotor.

The desire to achieve higher heat transfer augmentation for turbine blades is fueled by the increased power output and efficiency that is achievable with high turbine inlet temperatures. The use of internal cooling channels fitted with pin fin arrays serves as one method of accomplishing this goal. Consequently, the addition of pin fin arrays comes at the expense of increased pressure drop. Therefore the pin fin geometry must be judiciously chosen to achieve the required heat transfer rate while minimizing the associated pressure drop. This project culminates in the measurement of both pin fin force and array pressure drop as they related to changes in the array geometry. Specifically, the effects of Reynolds number, spanwise pin spacing, streamwise pin spacing, pin aspect ratio, and flow incidence angle. Direct two-component force measurement is achieved with a cantilever beam force sensor that uses highly sensitive piezoresistive strain gauges, relating the strain at the base of the beam to the applied force. With proper characterization, forces as small as one-tenth the weight of a paper clip are successfully measured. Additionally, array pressure drop measurements are achieved using static pressure taps. Experiments were conducted over a range of Reynolds numbers between 7,500 and 35,000. Changes in the spanwise pin spacing were shown to substantially alter the pin fin drag and array pressure drop, while changes in the streamwise pin spacing were less influential. The experimental results also showed a dramatic reduction in the pin fin drag and array pressure drop for an inline flow incidence angle. Finally, changes in the pin aspect ratio were shown to have little effect on the array pressure drop. / Master of Science

An optimal solution of vortex-induced vibrations of structures would be a time-domain numerical simulation that simultaneously solves the fluid flow and structural response. Yet, the requirements in terms of computing power remains a major obstacle for implementing such a simulation. On the other hand, lower- or reduced-order models provide an alternative for determining structural response to forcing by fluid flow. The objective of this thesis is to provide a consistent approach for the development of reduced-order models for the lift and drag on oscillating cylinders and the identification of their parameters. Amplitudes and phases of higher-order spectral moments of the lift and drag coefficients data are combined with approximate solutions of the representative models to determine their parameters. The results show that the amplitude and phase of the trispectrum could be used to model the lift on the oscillating cylinder under different excitation conditions. Moreover, the amplitude and phase of the cross-bispectrum could be used to establish the lift-drag relation for oscillating cylinders. A forced van der Pol equation is used to represent the lift on a transversely oscillating cylinder, and a parametrically excited van der Pol equation is used to model the lift coefficient on an inline oscillating cylinder. All cases of excitations lead to close values for the damping and nonlinear parameters in the van der Pol equation. Consequently, and as shown in this thesis, different excitation cases could be used to identify the parameters in the governing equations. Moreover, the results show that the drag coefficient could be derived from the lift coefficient through a square relation that takes into account the effects of the forced motions. / Ph. D.

The desire for smaller, faster, and more efficient products places a strain on thermal management in components ranging from gas turbine blades to computers. Heat exchangers that utilize internal cooling flows have shown promise in both of these industries. Although pin fins are often placed in the cooling channels to augment heat transfer, their addition comes at the expense of increased pressure drop. Consequently, the pin fin geometry must be judiciously chosen to achieve the desired heat transfer rate while minimizing the pressure drop and accompanying pumping requirements. This project culminates in the construction of a new test facility and the development of a unique force measurement methodology. Direct force measurement is achieved with a cantilever beam force sensor that uses sensitive piezoresistive strain gauges to simultaneously measure aerodynamic lift and drag forces on a pin fin. After eliminating the detrimental environmental influences, forces as small as one-tenth the weight of a paper clip are successfully measured. Although the drag of an infinitely long cylinder in uniform cross flow is well documented, the literature does not discuss the aerodynamic forces on a cylinder with an aspect ratio of unity in channel flow. Measured results indicate that the drag coefficient of a cylindrical pin in a single row array is greater than the drag coefficient of an infinite cylinder in cross flow. This phenomenon is believed to be caused by an augmentation of viscous drag on the pin fin induced by the increased viscous effects inherent in channel flow. / Master of Science

京都大学 / 新制・課程博士 / 博士(工学) / 甲第25236号 / 工博第5195号 / 京都大学大学院工学研究科社会基盤工学専攻 / (主査)教授 八木 知己, 教授 KIM Chul-Woo, 教授 高橋 良和 / 学位規則第4条第1項該当 / Doctor of Philosophy (Engineering) / Kyoto University / DFAM

Yes / The near-trapping effects on wave-induced dynamic seabed response and liquefaction close to a multi-cylinder foundation in storm wave conditions are examined. Momentary liquefaction near multi-cylinder structures is simulated using an integrated wave-structure-seabed interaction model. The proposed model is firstly validated for the case of interaction of wave and a four-cylinder structure, with a good agreement with available experimental measurements. The validated model is then applied to investigate the seabed response around a four-cylinder structure at 0° and 45° incident angles. The comparison of liquefaction potential around individual cylinders in an array shows that downstream cylinder is well protected from liquefaction by upstream cylinders. For a range of incident wave parameters, the comparison with the results for a single pile shows the amplification of pressure within the seabed induced by progressive wave. This phenomenon is similar to the near-trapping phenomenon of free surface elevation within a cylinder array. / Energy Technology Partnership (ETP), Wood Group / Full-text of this article will be released for public view at the end of the publisher embargo on 10 Oct 2021.

O presente estudo descreve as alterações metalúrgicas e topográficas do cilindro de blocos de motor de combustão interna fabricados em ferro fundido cinzento, após a realização de testes de durabilidade em dinamômetro. O motor testado tem tecnologia Flex-Fuel, e pode ser operado com qualquer proporção de mistura gasolina comum/etanol. Um motor foi testado com combustível etanol e o outro com gasolina comum, buscando-se assim representar duas condições extremas de trabalho em termos de combustível. A pesquisa constituiu-se da revisão bibliográfica e da parte prática que envolveu: a realização de testes de durabilidade de motor em dinamômetro; a caracterização metalúrgica e topográfica do cilindro após teste com ambos os combustíveis, utilizando técnicas de microscopia eletrônica de varredura, metalografia por microscopia óptica, perfilometria e quantificação de parâmetros de rugosidade, entre outros; a análise crítica, discussão dos resultados e apresentação da conclusão. / The present study approaches the metallurgical and topographical alterations on the cylinder of an internal combustion engine block made of gray cast iron, after durability dynamometer tests. The tested engine has flexible fuel technology (Flex-Fuel), and is capable to work with both gasoline/ethanol fuels, in any mixture proportion. One engine has been tested with ethanol and another one tested with gasoline, and so representing two extreme conditions on which the engine may work in terms of fuel. The research has been developed basically through the conceptual approach by the bibliographic review and the experimental steps that involved: general engine durability test at dynamometer bench; metallurgical and topographic characterization analysis at laboratory, after the test with both fuels, applying Scanning Electronic Microscopy SEM, optical microscopy metallography, profilometer and roughness parameters quantification, and so forth; critical analysis on the results, discussion and final conclusion.