Indiana University-Purdue University Indianapolis (IUPUI) / A wave rotor constant volume combustor (WRCVC) was designed and built as a collaborative work of Rolls Royce LibertyWorks, Indiana University-Purdue University at Indianapolis (IUPUI), and Purdue University, and ran experimental tests at Purdue's Zucrow Laboratories in 2009. Instrumentation of the WRCVC rig inlet flow included temperature and pressure transducers upstream of the venturi and at the fuel delivery plane. Other instrumentation included exhaust pressures and temperatures. In addition, ion sensors, dynamic pressure sensors, and accelerometers were used to instrument the rotating hardware. The rig hardware included inlet guide vanes directly in front of the rotating hardware, which together with concern for damage potential, prevented use of any pressure transducers at the entrance to the rotor. For this reason, a complete understanding of the conditions at the WRCVC inlet is unavailable, requiring simulations of the WRCVC to estimate the inlet pressure at a specific operating condition based on airflow. The operation of a WRCVC rig test is a sequence of events over a short time span. These events include introduction of the main air flow followed by time-sequenced delivery of fuel, lighting of the ignition source, and the combustion sequence. The fast changing conditions in the rig inlet hardware make necessary a time-dependent computation of the rig inlet section in order to simulate the overall rig operation. The chosen method for computing inlet section temperature and pressure was a time-dependent lumped volume model of the inlet section hardware, using a finite difference modified Euler predictor-corrector method for computing the continuity and energy equations. This is coupled with perfect gas prediction of venturi air and fuel flow rates, pressure drag losses at the fuel nozzles, pressure losses by mass addition of the fuel or nitrogen purge, friction losses at the inlet guide vanes, and a correlation of the non-dimensional flow characteristics of the WRCVC. The flow characteristics of the WRCVC are computed by varying the non-dimensional inlet stagnation pressure and the WRCVC's operational conditions, assuming constant rotational speed and inlet stagnation temperature. This thesis documents the creation of a computer simulation of the entire WRCVC rig, to understand the pressure losses in the inlet system and the dynamic coupling of the inlet section and the WRCVC, so that an accurate prediction of the WRCVC rotor inlet conditions can be computed. This includes the computational development of the WRCVC upstream rig dynamic model, the background behind supporting computations, and results for one test sequence. The computations provide a clear explanation of why the pressures at the rotor inlet differ so much from the upstream measured values. The pressure losses correlate very well with the computer predictions and the dynamic response tracks well with the estimation of measured airflow. A simple Fortran language computer program listing is included, which students can use to simulate charging or discharging of a container.

<p>The overall goal of the research presented in this thesis is to extend the physical understanding of the unsteady external aerodynamics associated with highly maneuverable delta-wing aircraft by using and developing novel, more efficient computational fluid dynamics (CFD) tools. More specific, the main purpose is to simulate and better understand the basic fluid phenomena, such as vortex breakdown, that limit the performance of delta-wing aircraft. The problem is approached by going from the most simple aircraft configuration - a pure delta wing - to more complex configurations. As the flow computations of delta wings at high angle of attack have a variety of unusual aspects that make accurate predictions challenging, best practices for the CFD codes used are developed and documented so as to raise their technology readiness level when applied to this class of flows.</p><p>Initially, emphasis is put on subsonic steady-state CFD simulations of stand-alone delta wings to keep the phenomenon of vortex breakdown as clean as possible. For half-span models it is established that the essential characteristics of vortex breakdown are captured by a structured CFD code. The influence of viscosity on vortex breakdown is studied and numerical results for the aerodynamic coefficients, the surface pressure distribution and breakdown locations are compared to experimental data where possible.</p><p>In a second step, structured grid generation issues, numerical aspects of the simulation of this nonlinear type of flow and the interaction of a forebody with a delta wing are explored.</p><p>Then, on an increasing level of complexity, time-accurate numerical studies are performed to resolve the unsteady flow field over half and full-span, stationary delta wings at high angle of attack. Both Euler and Detached Eddy Simulations (DES) are performed to predict the streamwise oscillations of the vortex breakdown location about some mean position, asymmetry in the breakdown location due to the interaction between the left and right vortices, as well as the rotation of the spiral structure downstream of breakdown in a time-accurate manner. The computed flow-field solutions are visualized and analyzed in a virtual-reality environment.</p><p>Ultimately, steady-state and time-dependent simulations of a full-scale fighter-type aircraft configuration in steady flight are performed using the advanced turbulence models and the detached-eddy simulation capability of an edge-based, unstructured flow solver. The computed results are compared to flight-test data.</p><p>The thesis also addresses algorithmic efficiency and presents a novel implicit-explicit algorithm, the Recursive Projection Method (RPM), for computations of both steady and unsteady flows. It is demonstrated that RPM can accelerate such computations by up to 2.5 times.</p>

<p>Severe floods caused by heavy autumn rains in year 2000 raised the question whether measures to reduce the damage from high water levels, in the area surrounding lake Glafsfjorden and along the river Byälven down to lake Vänern, are possible. One option is to reduce flow resistance along the river and thereby lower the maximum water level a given inflow would cause. Good knowledge of hydraulic and hydrological conditions is necessary in order to estimate the effect of such flow-reducing measures. In order to quantify such effect a 1-dimensional hydraulic flow model has been set up for the river Byälven using the software package MIKE 11. The model is more detailed, especially concerning topography and bathymetry, than earlier models used for studies of the river. Boundary conditions consist of measured inflows, the level of the lake Vänern and runoff calculated using the HBV-model. The model was calibrated for two different floods and a good fit to measured water levels was obtained for both these periods. Using the calibrated model critical sections, causing much flow resistance during high floods, were identified. With that knowledge different measures to reduce high water levels was adopted to the model both individually and combined with each other and the model was run with boundary conditions mainly from the flood in year 2000. The most radical measures simulated resulted in a lowering of the maximum water in the two largest reservoirs Glafsfjorden and Harefjorden with 78 and 97 cm respectively. A more modest combination of measures gave water levels 48 and 84 cm lower than a model run without changes. Some combinations of relatively small measures lowered the maximum water level by a few decimeters. The simulation results give good guidance to further investigations and decisions of actual changes. The model constitutes a useful tool when making flood maps of the area and if water level forecasts would be needed during future floods.</p> / <p>Allvarliga översvämningar i samband med höstregn år 2000 väckte frågan om det går att vidta åtgärder för att minska skadorna vid höga flöden i området runt Glafsfjorden och längs Byälvens sträckning ner till Vänern. Ett alternativ är att med åtgärder längs älven underlätta vattnets utflöde och på så sätt minska den högsta vattennivå ett givet flöde orsakar. God kunskap om hydrauliska och hydrologiska förhållanden behövs för att bedöma nyttan av olika åtgärdsalternativ. För att kunna avgöra effekterna av olika åtgärdsalternativ har en 1-dimensionell strömningsmodell satts upp för Byälven i programverktyget MIKE 11. Modellen är mer detaljerad, framförallt beträffande höjdinformationen, än tidigare modeller som använts för studier av Byälven varit. Randvillkor till modellen utgörs av registrerade inflöden, Vänerns vattenstånd och avrinning modellerad med HBV-modellen. Modellen har kalibrerats för två översvämningsperioder och god anpassning uppnåddes för de vattenstånd som dessa situationer representerar och med dess hjälp har sedan älven studerats och områden som bromsar flödet har kunnat identifieras. Med kunskap om vilka områden som begränsar flödet mest har ett antal olika åtgärder simulerats i modellen, både var för sig och kombinerade med varandra. Randvillkoren för översvämningen år 2000 behölls och förändringarna lades in i modellen. De extremaste åtgärderna som simulerats resulterade i minskningar av de högsta vattennivåerna i de två största vattenmagasinen Glafsfjorden och Harefjorden med 78 respektive 97 cm. Mer realistiska åtgärdspaket gav minskningar med i storleksordningen 48 respektive 84 cm och även relativt små ingrepp gav minskningar på några decimeter. Simuleringsresultaten ger god vägledning för vidare undersökningar av och beslut om konkreta åtgärder i Byälven. Dessutom utgör modellen ett bra verktyg för att ta fram översvämningskartor och för att prognostisera vattennivåer vid nya översvämningssituationer.</p>

The main objective of this study is to explore the complex fluid flow phenomena that result in the generation of a high frequency noise in counterbalance valves through an experimental and numerical investigation of the flow. Once the influence of the different components involved in noise generation is established, a secondary objective is the introduction of design modifications that eliminate the undesired effect without altering the operation envelope or the performance of the valve. A hydraulic test bench was used to carry out an experimental investigation of the noise generation process. A computer based data acquisition system was used to record pressure fluctuations, flowrates and hydraulic oil temperatures in a production valve under a variety of operational conditions. Extensive experimental measurements and numerical modeling lead to the hypothesis that noise generation is the result of an acoustic resonance triggered by shear layer instability at the valve inlet. The pressure gradients developed when the shear layer entrains the stagnant fluid in the valve main cavity cause the layer to become unstable and oscillate. The oscillation frequency will depend on a great number of factors such as valve geometry, pressure and velocity gradients and the density and viscosity of the fluid. It is postulated that the observed noise is generated when this frequency matches one of the resonant frequencies of the valve cavity. The proposed mechanism is theoretically poorly understood and well beyond simplified analysis, its accurate numerical simulation is computational very intensive requiring sophisticated CFD codes. The numerical investigation was carried out using STAR–CCM+, a commercially available CFD code featuring 3-D capabilities and sophisticated turbulence modeling. Streamline, pressure, velocity-vector and velocity-scalar plots were obtained for several valve configurations using steady and unsteady state flow simulations. An experimental and numerical analysis of an alternative valve geometry was carried out. Experimental results demonstrated a greatly reduced instability range. The numerical analysis of the unsteady behavior of the shear-layer streamlines for both valves yielded results that were compatible with the experimental work. The results of this investigation promise a great positive impact on the design of this type of hydraulic valves.

The study of the motion of manoeuvring aircraft has traditionally considered the aircraft to be rigid. This simplifying assumption has been shown to give quite accurate results for the flight dynamics of many aircraft types. As modern transport aircraft have developed however, there has been a marked increase in the size and weight of these aircraft. This trend is likely to continue with the development of future blended-wing-body and supersonic transport aircraft. This increase in size and weight has brought about a unique set of aeroelastic and handling quality issues. The aerodynamic forces and moments acting on an aeroplane have traditionally been represented using the aerodynamic derivative approach. It has been shown that this quasisteady aerodynamic model inadequately predicts the aircraft’s stability characteristics, and that the inclusion of unsteady aerodynamics “greatly improves the fidelity” of aircraft models. This thesis thus presents a novel numerical simulation of an aeroelastic aeroplane for real-time analysis. The model is built around the standard six degree-of-freedom equations of motion for a rigid aeroplane using the mean-axes system, and includes unsteady aerodynamics and structural dynamics. This is suitable for pilot-in-the-loop simulation, handling qualities and flight loads analysis, and control law development. The dynamics of the structure are modelled as a set of normal modes, and the equations of motion are realised in state-space form. The unsteady aerodynamic forces acting on the aeroplane are described by an indicial state-space model, including unsteady tailplane downwash and compressibility effects. An implementation of the model is presented in the MATLAB/ Simulink environment. The interaction between the flight control system, the aeroelastic system and the rigidbody motion of the aeroplane can result in degraded handling qualities, excessive actuator control, and fatigue problems. The introduction of load alleviation systems for the management of loads due to manoeuvres and gusts is also likely to result in the handling qualities of the aeroplane being degraded. This thesis presents a number of studies into the impact of structural dynamics, unsteady aerodynamics, and load alleviation on the handling qualities of a flexible civil transport aeroplane. The handling qualities of the aeroplane are assessed against a number of different handling qualities criteria and flying specifications, including the Neal-Smith, Bandwidth, and CAP criterion. It is shown that aeroelastic effects alter the longitudinal and lateral-directional characteristics of the aeroplane, resulting in degraded handling qualities. Manoeuvre and gust load alleviation are similarly found to degrade handling qualities, while active mode control is shown to offer the possibility of improved handling qualities.

Evaluating the performance of filter media for filtration applications is essential to assure design engineers and users that filter device will deliver promised performance for specific applications under the environmental stress. The study of particle loading characteristics of filter media in the laboratory setting is typically performed under the steady flow conditions, i.e., at the constant particle concentration and flow rate. In reality, filtration products are operated under the situations that the flow rate and mass concentration of particles are varied in time. The success of translating the laboratory data to estimate the performance of filter media in the fields is thus limited. It is necessary to investigate the performance of filter media under the real-world conditions, i.e., unsteady flow rate and mass concentration to bridge the gap. The overall goals of this research are (1) to study the performance of filter medium under unsteady conditions (i.e., the performance of respirator filter media under simulated breathing conditions); and (2) to investigate the issue of non-uniform particle deposition on HVAC filter panels. A new experimental setup was developed to accomplish the former goal. Numerical modeling tool, Computational fluid dynamics (CFD), was applied to achieve the latter objective……

A two-dimensional (2-D), mathematical model is adopted to investigate the development of circulation patterns for compressible, laminar, and shear driven flow inside a rectangular cavity. The bottom of the cavity is free to move at a specified speed and the aspect ratio of the cavity is changed from 1.0 to 1.5. The vertical sides and the bottom of the cavity are assumed insulated. The cavity is filled with a compressible fluid with Prandtl number, Pr =1. The governing equations are solved numerically using the commercial Computational Fluid Dynamics (CFD) package ANSYS FLUENT 2015 and compared with the results for the primitive variables of the problem obtained using in house CFD code based on Coupled Modified Strongly Implicit Procedure (CMSIP). The simulations are carried out for the unsteady, lid driven cavity flow problem with moving boundary (bottom) for different Reynolds number, Mach numbers, bottom velocities and high initial pressure and temperature.

O escoamento hidrodinâmico é descrito por meio das equações da continuidade e da quantidade de movimento em conjunto com as condições iniciais e de contorno. Neste trabalho de pesquisa, o escoamento não permanente na saída de um canal retangular curto de declividade nula, resultante da abertura rápida de uma comporta foi experimentalmente investigado. Para determinar o coeficiente de pressão e o comprimento do trecho a montante da comporta com distribuição de pressão não hidrostática, usou-se um modelo unidimensional, descrito pelas equações de Saint Venant, para representar o fenômeno. Utilizou-se o método de Lax-Wendroff para obter as curvas características do fenômeno e compará-las com os resultados obtidos experimentalmente. Várias comparações foram feitas possibilitando a avaliação de valores para os parâmetros que definem a região com distribuição de pressão não hidrostática. Os resultados calculados para os valores de parâmetros constantes, estão de acordo com os resultados experimentais. Os resultados experimentais ajustam-se melhor quando os valores dos parâmetros variam no decorrer do tempo. / The hydrodynamic flow is described by means of the equations of the continuity and momentum together with the initial conditions, the boundary conditions. In this research work, the unsteady flow at the exit of a short rectangular channel of zero slope, resulting from the fast opening of a floodgate was investigated experimentally. For determining the pressure coefficient and the length of the upstream reach of the floodgate with non-hydrostatic pressure distribution, a one-dimensional model described by Saint Venant equations is used to represent the phenomenon. The Lax-Wendroff method was used to obtain the curves characteristics of the phenomenon and to compare them with the results obtained experimentally. Several comparisons were made to facilitate the evaluation of values for the parameters that define the region with non-hydrostatic pressure distribution. The results calculated for the values of constant parameters are in agreement with the experimental results. The experimental results are better adjusted when the values of the parameters vary in elapsing of the time.

Les émissions polluantes à l’échappement des véhicules automobiles sont l'une des principales sources de pollution de l'air dans le monde d'aujourd'hui. Par conséquent, la législation a évolué afin de limiter ces émissions. L'un des aspects clés pour répondre consiste à bien maîtriser les échanges gazeux au sein du moteur à combustion interne. Cette amélioration est possible par l'optimisation de répartiteurs d'admission d'air. Dans ces répartiteurs d'admission d'air, la maitrise de l’écoulement de type tumble est une piste de progrès. Des volets sont installés à la sortie du répartiteur afin d'améliorer le rapport de tumble et donc le mélange air-carburant (VTS-Variable Tumble System). Une autre caractéristique de l'écoulement à l'intérieur des répartiteurs est l'effet des écoulements pulsés qui engendrent des fluctuations de pression assez importante. Par conséquent, le but de cette étude consiste à simuler le flux d'air pulsé à l'intérieur des répartiteurs d'admission et à identifier l'effet des pulsations de pression sur les composants actifs tels que les volets. Le travail de simulation dans la présente thèse a été effectué à partir du code open source CFD OpenFOAM. Dans un premier temps, l'effet des pulsations de pression est simulé à l'intérieur d'un tube d'acier et une méthodologie de simulation est développée. Les résultats de la simulation sont validés à partir de résultats expérimentaux obtenus sur un dispositif spécifique, le banc dynamique. Ensuite, des simulations ont été effectuées sur le répartiteur d'admission principal avec des volets. Tout d’abord, les simulations sont effectuées en régime permanent avec cinq positions d'ouverture différentes du clapet. Les forces et les moments agissant sur le volet en régime permanent sont obtenus et analysés. Puis, des simulations en régime transitoire avec des effets de pulsation de pression sont effectuées. Les résultats de la simulation instationnaire sont comparés aux résultats expérimentaux en termes de fluctuations de pression relative. Les effets des pulsations de pression sur les forces aérodynamiques et les moments agissant sur les volets sont analysés et commentés. / The exhaust emissions from automobiles are one of the major sources of air pollution in today’s world. Thence,research and development is the key feature of the modern automotive industries to meet strict emission legislation. One of the key aspects to meet these requirements is to improve the gas exchange process within internal combustion engines. It is possible by the design optimization of the air intake manifolds for internal combustion engines. One of such advancement in air intake manifolds is variable tumble systems (VTS). In VTS system, tumble flaps are installed at the exit of the manifold runner in order to improve tumble ratio and hence air-fuel mixing. Another feature of the flow inside the intake manifolds is pressure pulsation effect. Therefore, the aim of the Ph.D. work is to simulate the pulsating air flow inside the air intake manifolds and to identify the effect of the pressure pulsations on the active components like tumble flaps. The simulation work in the present thesis has been carried out on open source CFD code OpenFOAM. In a first step, the effect of pressure pulsations is simulated inside a steel tube and a simulation methodology is developed. The results of the simulation are validated on a specific experimental device, the dynamic flow bench. Then,simulations have been carried out on the main intake manifold with tumble flaps. Firstly, the simulations are performed with five different opening positions of the tumble flap in a steady state configuration. The forces and moments acting on the flap in steady state are obtained and analyzed. Then, unsteady simulations with pressure pulsation effects are performed. The results of obtained from unsteady simulation are compared with the experimental results in terms of relative pressure fluctuations. The effect of the pressure pulsation on the aerodynamic forces and moments acting on the tumble flaps are analyzed and explained.

abstract: Flooding is a critical issue around the world, and the absence of comprehension of watershed hydrologic reaction results in lack of lead-time for flood forecasting and expensive harm to property and life. It happens when water flows due to extreme rainfall storm, dam breach or snowmelt exceeds the capacity of river system reservoirs and channels. The objective of this research was to develop a methodology for determining a time series operation for releases through control gates of river-reservoir systems during flooding events in a real-time using one- and/or two-dimensional modeling of flows through river-reservoir systems. The optimization-simulation methodology interfaces several simulation-software coupled together with an optimization model solved by genetic algorithm coded in MATLAB. These software include the U.S. Army Corps of Engineers HEC-RAS linked the genetic algorithm in MATLAB to come up with an optimization-simulation model for time series gate openings to control downstream elevations. The model involves using the one- and two-dimensional ability in HEC-RAS to perform hydrodynamic routing with high-resolution raster Digital Elevation Models. Also, the model uses both real-time gridded- and gaged-rainfall data in addition to a model for forecasting future rainfall-data. This new model has been developed to manage reservoir release schedules before, during, and after an extraordinary rainfall event that could cause extreme flooding. Further to observe and control downstream water surface elevations to avoid exceedance of threshold of flood levels in target cells in the downstream area of study, and to minimize the damage and direct effects in both the up and downstream. The application of the complete optimization-simulation model was applied to a portion of the Cumberland River System in Nashville, Tennessee for the flooding event of May 2010. The objective of this application is to demonstrate the applicability of the model for minimizing flood damages for an actual flood event in real-time on an actual river basin. The purpose of the application in a real-time framework would be to minimize the flood damages at Nashville, Tennessee by keeping the flood stages under the 100-year flood stage. This application also compared the three unsteady flow simulation scenarios: one-dimensional, two-dimensional and combined one- and two-dimensional unsteady flow. / Dissertation/Thesis / Doctoral Dissertation Civil, Environmental and Sustainable Engineering 2019