<div>The dynamic performance of the momentum and thermal boundary layer linked to the acoustic response dictate the efficiency of heat exchangers and the operational limits of fluid machinery. The specific time required by the boundary layer to establish or adapt to the free stream variations is vital to optimize flow control strategies as well as the thermal management of fluid systems. The proper understanding of the wall fluxes, separated flow regions and free stream response to transient conditions becomes the fulcrum of the further improvement of fluid machinery performance and endurance. Throughout this dissertation the establishment sequence and the main parameters dictating the acoustic response and the boundary layer settlement are quantified together with their implication on the wall fluxes and boundary layer detachment. </div><div><br></div><div>Unsteady Reynolds Average Navier Stokes evaluations, Large Eddy Simulations, Direct Numerical Simulations and wind tunnel experiments are exploited to analyze the transient behavior of attached and detached flow aerodynamics. The core of the research is built upon URANS simulations allowing the realization of multiple detailed parametric analyses. Thanks to its reduced computational cost, hundreds of transient flow evaluations are carried out, enabling the determination of the establishment sequence, the main flow features and relevant non-dimensional numbers. The URANS methodology is verified against experimental and analytic results on the flow conditions of the study. The Large Eddy Simulations and Direct Numerical Simulations allow further characterization of the near wall flow region behavior with much higher resolution while providing an additional source of verification for the coarser numerical tools. An experimental campaign on a novel full visual access linear wind tunnel explores the impact of mean flow sudden accelerations on the boundary layer detachment and reattachment phenomena over an ad-hoc wall mounted hump. The wind tunnel is designed based on the premises of: full visual access, spatial and temporal stability of total and static pressure together with the total temperature and fast flow settlement, minimizing the start-up phase duration of the wind tunnel. A wall mounted hump that mimics the behavior of the aft portion of a low pressure turbine is inserted in the wind tunnel guaranteeing a 2D flow separation phenomena. After steady state test article characterization series of sudden flow discharge experiments reveal the impact of mean flow transients on the boundary layer detachment inception. Finally, taking advantage of the knowledge on transient flow performance, optimum flow control mechanisms to abate boundary layer detachment are proposed. The recommended control approach effectively prevents the boundary layer separation while minimizing the energy requirement.</div>

The laminar flow in the small bronchial tubes is quite complex due to the presence of vortex-dominated, secondary flows. Contributing to this complexity are the geometrical characteristics of the bronchial tubes that include asymmetric and nonplanar branching. These secondary flow fields play a crucial role in particle deposition; however, the actual mechanisms that determine the particle distributions are not fully understood. The research reported here increases understanding of this phenomenon by studying flow structure/ particle interaction in the small bronchial tubes for steady and unsteady respiratory conditions. Specifically, the effects of simultaneous nonplanar and asymmetric branching were investigated. The nonplanar model was generated by applying a 90◦ out-of-plane rotation to the third-generation branches. Steady-state inspiratory flows for a Reynolds number of 1,000 and unsteady periodic flows with a 30-respiration-per-minute breathing frequency were simulated in three-generation, asymmetric, planar and nonplanar models. The asymmetry and nonplanarity produced asymmetric secondary flow patterns and unequal mass flow partitioning in the third-generation branches. Ten micron water droplet deposition in the nonplanar model was found to be significantly different from the planar model, demonstrating the impact of simultaneous nonplanar and asymmetric branching. The unsteady nature of the flow also affected particle deposition. Particles released at the same instantaneous inflow conditions during off-peak inhalation conditions, generated significantly different particle deposition patterns. The differences were attributed to the high temporal variations of the fluid velocities at these off-peak times and history effects in the flows. It was also observed that the initial particle velocities had a significant impact on particle deposition. The study of flow structure and particle interaction was facilitated by the development of a novel visualization technique that employs finite-time Lyapunov exponents (FTLE). This research provides a better understanding of the fluid dynamics driving the particle deposition in the bronchial tubes.

Este trabalho representa a continuidade de estudos envolvendo a problemática dos escoamentos com superfície livre, contemplando a análise do fenômeno transiente em canais, a partir do modelo matemático unidimensional de Saint-Venant. Para tanto, é desenvolvido um modelo computacional em linguagem FORTRAN, capaz de avaliar o comportamento do escoamento não permanente. As equações hidrodinâmicas completas são discretizadas por um esquema completamente implícito de diferenças finitas e aplicadas no modelo computacional para a avaliação de dois casos. O modelo é previamente testado para um caso simples, cujos resultados são analisados viabilizando o modelo. No primeiro caso, o modelo é aplicado ao canal de alimentação da Usina Hidrelétrica Monjolinho em São Carlos-SP, para avaliar a necessidade de vertedouro quando se dá o fechamento brusco da turbina, e a ocorrência da entrada de ar na mesma quando da sua abertura repentina. No segundo caso, procurou-se avaliar o desenvolvimento do escoamento no Canal do Trabalhador, responsável pelo abastecimento da cidade de Fortaleza-CE. Com manobras de enchimento e esvaziamento do sistema, é possível determinar o tempo de antecedência de liga-desliga do sistema de recalque a partir das alturas dágua e velocidades de ocorrência, permitindo também a automação para as operações de controle. Em ambos os casos o modelo reproduziu resultados que ilustram com coerência os conceitos pré-estabelecidos, constituindo numa ferramenta útil para análise do fenômeno transiente nos escoamentos em condutos livres. / This work presents a computational model developed in FORTRAN language for the study of unsteady open-channel flows with the use of Saint-Venant one-dimensional equation. The discretization of hydrodynamic equations are presented in a completely implicit method of finite differences and applied in the model for the investigation of two cases, besides the one used previously to test the model. In the first case, the model is applied for a channel that supplies the Monjolinho hydroelectric plant in Sao Carlos SP, aiming to evaluate the need of a spillway when the turbine is closed and the flow abruptly stopped, as well as the occurrence of air entering the turbine when it is opened instantaneously. In the second case, the model simulates the development of the flow in the Trabalhador channel, responsible for the water supply in the city of Fortaleza - CE, in order to make possible the automation of operational control, based on data of flow velocity and water level. In both cases the model is presented as a useful tool for the analysis of unsteady open-channel flows, showing results and coherency with theory.

In this thesis, the seminal work of Stokes concerning the ow of a Navier-Stokesuid due to a suddenly accelerated or oscillating plate and the ow due to torsionaloscillations of an innitely long rod in a Navier-Stokes uid is extended to a uid withpressure dependent viscosity. The viscosity of many uids varies signicantly withpressure, a fact recognized by Stokes; and Barus, in fact, conducted experiments thatshowed that the variation of the viscosity with pressure was exponential. Given sucha tremendous variation, we study how this change in viscosity aects the nature of thesolution to these problems. We nd that the velocity eld, and hence the structureof the vorticity and the shear stress at the walls for uids with pressure dependentviscosity, are markedly dierent from those for the classical Navier-Stokes uid.

A prerequisite for aeroelastic stability prediction inturbomachines is the understanding of the fluctuatingaerodynamic forces acting on the blades. Unsteady transonicflows are complex because of mutual interactions betweentravelling pressure waves, outlet disturbances, shock motion,and fluctuating turbulent boundary layers. Complex phenomenaappear in the shock/boundary layer region and produce phaselags and high time harmonics, which can give a significantcontribution to the overall unsteady lift and torque, andtherefore affect flutter boundaries, cause large localstresses, or even severely damage the turbomachine. The present research work is concerned with theunderstanding of phenomena associated with travelling waves innon-uniform transonic flows and how they affect the unsteadypressure distribution on the surface as well as the far fieldradiated sound. In similitude with turbomachines potentialinteraction, the emphasis was put on the unsteady interactionof upstream propagating acoustic waves with an oscillatingshock in 2D and 3D nozzle flows. Both numerical andexperimental studies are carried out and compared with eachother. Results shows that the unsteady pressure distribution, bothon the bump surface and within the channel, results from thesuperposition of upstream and downstream propagating waves. Itis believed that outlet pressure perturbations propagateupstream in the nozzle, interact in the high subsonic flowregion according to the acoustic blockage theory, and arepartly reflected or absorbed by the oscillating shock,depending on the frequency of the perturbations and theintensity of the SBLI. Furthermore the shock motion amplitudeis found to be related to the mean flow gradients and localwave length of the perturbations rather than to the shockboundary layer interaction. The phase angle between incomingpressure perturbations and the shock motion increases with theperturbation frequency but also depends on the intensity of theSBLI. Additionally the phase angle "shift" observed underneaththe shock location linearly increases with the perturbationfrequency and the shock strength. Such phase shift is criticalregarding aeroelastic stability and might have a significantimpact on the phase angle of the overall aerodynamic forceacting on the blade and shift the aerodynamic damping fromstable to exciting. <b>Keywords:</b>Shock Boundary Layer Interaction, ShockMotion, Unsteady Flows, Nozzle Flows, Potential Interaction,Back Pressure Perturbations.

In this thesis, the seminal work of Stokes concerning the ow of a Navier-Stokesuid due to a suddenly accelerated or oscillating plate and the ow due to torsionaloscillations of an innitely long rod in a Navier-Stokes uid is extended to a uid withpressure dependent viscosity. The viscosity of many uids varies signicantly withpressure, a fact recognized by Stokes; and Barus, in fact, conducted experiments thatshowed that the variation of the viscosity with pressure was exponential. Given sucha tremendous variation, we study how this change in viscosity aects the nature of thesolution to these problems. We nd that the velocity eld, and hence the structureof the vorticity and the shear stress at the walls for uids with pressure dependentviscosity, are markedly dierent from those for the classical Navier-Stokes uid.

<p>A prerequisite for aeroelastic stability prediction inturbomachines is the understanding of the fluctuatingaerodynamic forces acting on the blades. Unsteady transonicflows are complex because of mutual interactions betweentravelling pressure waves, outlet disturbances, shock motion,and fluctuating turbulent boundary layers. Complex phenomenaappear in the shock/boundary layer region and produce phaselags and high time harmonics, which can give a significantcontribution to the overall unsteady lift and torque, andtherefore affect flutter boundaries, cause large localstresses, or even severely damage the turbomachine.</p><p>The present research work is concerned with theunderstanding of phenomena associated with travelling waves innon-uniform transonic flows and how they affect the unsteadypressure distribution on the surface as well as the far fieldradiated sound. In similitude with turbomachines potentialinteraction, the emphasis was put on the unsteady interactionof upstream propagating acoustic waves with an oscillatingshock in 2D and 3D nozzle flows. Both numerical andexperimental studies are carried out and compared with eachother.</p><p>Results shows that the unsteady pressure distribution, bothon the bump surface and within the channel, results from thesuperposition of upstream and downstream propagating waves. Itis believed that outlet pressure perturbations propagateupstream in the nozzle, interact in the high subsonic flowregion according to the acoustic blockage theory, and arepartly reflected or absorbed by the oscillating shock,depending on the frequency of the perturbations and theintensity of the SBLI. Furthermore the shock motion amplitudeis found to be related to the mean flow gradients and localwave length of the perturbations rather than to the shockboundary layer interaction. The phase angle between incomingpressure perturbations and the shock motion increases with theperturbation frequency but also depends on the intensity of theSBLI. Additionally the phase angle "shift" observed underneaththe shock location linearly increases with the perturbationfrequency and the shock strength. Such phase shift is criticalregarding aeroelastic stability and might have a significantimpact on the phase angle of the overall aerodynamic forceacting on the blade and shift the aerodynamic damping fromstable to exciting.</p><p><b>Keywords:</b>Shock Boundary Layer Interaction, ShockMotion, Unsteady Flows, Nozzle Flows, Potential Interaction,Back Pressure Perturbations.</p>

Este trabalho representa a continuidade de estudos envolvendo a problemática dos escoamentos com superfície livre, contemplando a análise do fenômeno transiente em canais, a partir do modelo matemático unidimensional de Saint-Venant. Para tanto, é desenvolvido um modelo computacional em linguagem FORTRAN, capaz de avaliar o comportamento do escoamento não permanente. As equações hidrodinâmicas completas são discretizadas por um esquema completamente implícito de diferenças finitas e aplicadas no modelo computacional para a avaliação de dois casos. O modelo é previamente testado para um caso simples, cujos resultados são analisados viabilizando o modelo. No primeiro caso, o modelo é aplicado ao canal de alimentação da Usina Hidrelétrica Monjolinho em São Carlos-SP, para avaliar a necessidade de vertedouro quando se dá o fechamento brusco da turbina, e a ocorrência da entrada de ar na mesma quando da sua abertura repentina. No segundo caso, procurou-se avaliar o desenvolvimento do escoamento no Canal do Trabalhador, responsável pelo abastecimento da cidade de Fortaleza-CE. Com manobras de enchimento e esvaziamento do sistema, é possível determinar o tempo de antecedência de liga-desliga do sistema de recalque a partir das alturas dágua e velocidades de ocorrência, permitindo também a automação para as operações de controle. Em ambos os casos o modelo reproduziu resultados que ilustram com coerência os conceitos pré-estabelecidos, constituindo numa ferramenta útil para análise do fenômeno transiente nos escoamentos em condutos livres. / This work presents a computational model developed in FORTRAN language for the study of unsteady open-channel flows with the use of Saint-Venant one-dimensional equation. The discretization of hydrodynamic equations are presented in a completely implicit method of finite differences and applied in the model for the investigation of two cases, besides the one used previously to test the model. In the first case, the model is applied for a channel that supplies the Monjolinho hydroelectric plant in Sao Carlos SP, aiming to evaluate the need of a spillway when the turbine is closed and the flow abruptly stopped, as well as the occurrence of air entering the turbine when it is opened instantaneously. In the second case, the model simulates the development of the flow in the Trabalhador channel, responsible for the water supply in the city of Fortaleza - CE, in order to make possible the automation of operational control, based on data of flow velocity and water level. In both cases the model is presented as a useful tool for the analysis of unsteady open-channel flows, showing results and coherency with theory.

京都大学 / 0048 / 新制・課程博士 / 博士(工学) / 甲第22055号 / 工博第4636号 / 新制||工||1723(附属図書館) / 京都大学大学院工学研究科都市社会工学専攻 / (主査)教授 細田 尚, 教授 戸田 圭一, 准教授 音田 慎一郎 / 学位規則第4条第1項該当 / Doctor of Philosophy (Engineering) / Kyoto University / DFAM

This work deals with the design of a hybrid time integrator that couples spatially explicit and implicit time integrators. In order to cope with the industrial solver of Ariane Group called FLUSEPA, the explicit scheme of Heun and the implicit scheme of Crank-Nicolson are hybridized using the transition parameter : the whole technique is called AION time integration. The latter is studied into details with special focus on spectral behaviour and on its ability to keep the accuracy. It is shown that the hybrid technique has interesting dissipation and dispersion properties while maintaining precision and avoiding spurious waves. Moreover, this hybrid approach is validated on several academic test cases for both convective and diffusive fluxes. And as expected the method is more interesting in term of computational time than standard time integrators. For the extension of this hybrid approach to the temporal adaptive method implemented in FLUSEPA, it was necessary to improve some treatments in order to maintain conservation and acceptable spectral properties. Finally the hybrid time integration was also applied to a RANS/LES turbulent test case with interesting computational time while capturing the flow physics.