The objective of this research is to assess the performance of two popularReynolds-averaged Navier-Stokes (RANS) models, standard k-E and k-w, andto suggest modifications to improve model predictions for high-speed flows. Numerical simulations of turbulent ow past a at plate are performed at M1 = 2:25; 7:2.The results from these two Mach number cases are compared with Direct NumericalSimulation (DNS) results from Pirozzoli et al. (2004) and experimental results fromHorstman & Owen (1975). The effect of the Boussinesq coefficient (Cu) and turbulenttransport coefficients (sigmak; sigmaE; sigma; sigma*) on the boundary layer ow is examined. Further,the performance of a new model with realizability-based correction to Cu and corresponding modifications to sigma; sigma* is examined. The modification to Cu is based oncontrolling the ratio of production to dissipation of kinetic energy (P/E=1). The firstchoice of P/E = 1 ensures that there is no accumulation of kinetic energy in stagnation or free-stream regions of the ow. The second choice of P/E= 1:6 holds underthe assumption of a homogeneous shear ow. It is observed that the new model'sperformance is similar to that of the existing RANS models, which is expected for asimple ow over a at plate. Finally, the role of turbulent Prandtl number (Prt) intemperature and density predictions is established. The results indicate that the k-wmodel's performance is better compared to that of the standard k-E model for highMach number flows. A modification to Cu must be accompanied with correspondingchanges to sigmak; sigmaE; sigma; sigma* for an accurate log-layer prediction. The results also indicate that a Prt variation is required across the boundary layer for improved temperatureand density predictions in high-speed flows.

The objective of this research is to assess the performance of two popularReynolds-averaged Navier-Stokes (RANS) models, standard k-E and k-w, andto suggest modifications to improve model predictions for high-speed flows. Numerical simulations of turbulent ow past a at plate are performed at M1 = 2:25; 7:2.The results from these two Mach number cases are compared with Direct NumericalSimulation (DNS) results from Pirozzoli et al. (2004) and experimental results fromHorstman & Owen (1975). The effect of the Boussinesq coefficient (Cu) and turbulenttransport coefficients (sigmak; sigmaE; sigma; sigma*) on the boundary layer ow is examined. Further,the performance of a new model with realizability-based correction to Cu and corresponding modifications to sigma; sigma* is examined. The modification to Cu is based oncontrolling the ratio of production to dissipation of kinetic energy (P/E=1). The firstchoice of P/E = 1 ensures that there is no accumulation of kinetic energy in stagnation or free-stream regions of the ow. The second choice of P/E= 1:6 holds underthe assumption of a homogeneous shear ow. It is observed that the new model'sperformance is similar to that of the existing RANS models, which is expected for asimple ow over a at plate. Finally, the role of turbulent Prandtl number (Prt) intemperature and density predictions is established. The results indicate that the k-wmodel's performance is better compared to that of the standard k-E model for highMach number flows. A modification to Cu must be accompanied with correspondingchanges to sigmak; sigmaE; sigma; sigma* for an accurate log-layer prediction. The results also indicate that a Prt variation is required across the boundary layer for improved temperatureand density predictions in high-speed flows.

Zugl.: Braunschweig, Techn. Univ., Diss., 2007

This thesis investigates three problems relating to thermoconvective stability in porous media. These are (i) the stability of an inclined boundary layer flow to vortex type instability, (ii) front propagation in the Darcy-B´enard problem and (iii) the onset of Prantdl-Darcy convection in a horizontal porous layer subject to a horizontal pressure gradient. The nonlinear, elliptic governing equations for the inclined boundary layer flow are discretised using finite differences and solved using an implicit, MultiGrid Full Approximation Scheme. In addition to the basic steady state three configurations are examined: (i) unforced disturbances, (ii) global forced disturbances, and (iii) leading edge forced disturbances. The unforced inclined boundary layer is shown to be convectively unstable to vortex-type instabilities. The forced vortex system is found to produce critical distances in good agreement with parabolic simulations. The speed of propagation and the pattern formed behind a propagating front in the Darcy-B´enard problem are examined using weakly nonlinear analysis and through numerical solution of the fully nonlinear governing equations for both two and three dimensional flows. The unifying theory of Ebert and van Saarloos (Ebert and van Saarloos (1998)) for pulled fronts is found to describe the behaviour well in two dimensions, but the situation in three dimensions is more complex with combinations of transverse and longitudinal rolls occurring. A linear perturbation analysis of the onset of Prandtl-Darcy convection in a horizontal porous layer subject to a horizontal pressure gradient indicates that the flow becomes more stable as the underlying flow increases, and that the wavelength of the most dangerous disturbances also increases with the strength of the underlying flow. Asymptotic analyses for small and large underlying flow and large Prandtl number are carried out and results compared to those of the linear perturbation analysis.

Invasive aquatic species such as Hydrilla verticillata (hydrilla) have become a pervasive and nearly ineradicable part of the waterways of the American south. Hydrilla is an aggressive colonizer; grows rapidly and rapidly blocks flow areas, which greatly reduces the capacity of water supply canals. Hydrilla grows up through the water column and is present throughout flow zones that are typically assumed to be free flowing and without resistance, other than that transmitted via the mechanics of a Newtonian fluid. Hydrilla is highly flexible and its morphology in the flow field is dependent on many parameters, including flow, growth stage, cross-section geometry and substrate. Traditional methods of calculating canal flow capacities assume that resistance to flow originates at the boundary of the channel. These methods typically attempt to account for vegetation by increasing resistance coefficients, which are associated with the boundary of the canal. A combination of field studies and experimentation in three separate laboratory channels was used to characterize the behavior of hydrilla and its impacts on open-channel flow. This work developed relationships for energy losses of flow within the vegetation, as well as velocity gradients within the vegetation and through the vegetation water interface to the open water. The information developed in this investigation was used to develop a model of the cross-section of flow with vegetation growing in the center of the channel. The model is based on the Prandtlvon Kármán universal-velocity-distribution law; and uses modifications to the method of calculating the hydraulic radius, to account for the increased frictional elements and reduced flow areas in the canal cross-section. A simple function was developed to estimate the remaining flow capacity in a canal as a function of the remaining unblocked area. The Prandtl-von Kármán universal-velocity-distribution law, together with modifications to the method for calculating the hydraulic radius, can improve estimates of the flow in channels impacted by submerged aquatic vegetation. The effects of a broad range of parameters can thus be represented by a relatively simple function, which was developed in this project.

A numerical study is performed to model the 2-D axisymmetric turbulent flow of supercritical water flowing upward in vertical pipes with constant wall heat fluxes, using ANSYS CFX v14.5. This study was aimed to use CFD in analyzing supercritical flow instability and heat transfer. Two types of flow instabilities are analyzed and results are compared with 1-D non-linear code solutions. Also, conditions for approximating the thresholds of instabilities based on steady-state results are assessed. It is determined that the results of instability thresholds obtained using the k-ɛ and the SST models are similar. Also the results of CFD and 1-D codes are different mainly due to the difference in the pressure drop predictions. Moreover, approximating the flow instability threshold by the conditions proposed holds true for a CFD solution. Results also indicate that Prt does not have a noticeable effect on the instability threshold for the cases examined.

Zugl.: Aachen, Techn. Hochsch., Diss., 2007.

Zugl.: Aachen, Techn. Hochsch., Diss., 2007

In dieser Arbeit entwickeln und untersuchen wir ein numerisches Schema für die eindimensionale Wellengleichung mit Hysterese für unterschiedliche Arten von Randbedingungen. Diese Gleichung ist ein Modell für die Longitudinal- oder Torsionsschwingungen eines homogenen Stabes unter dem Einfluß einer uniaxialen äußeren Kraftdichte, wobei wir ein elastoplastisches Materialgesetz annehmen. Hysterese-Operatoren sind ratenunabhängige Volterra-Operatoren, die Zeitfunktionen in Zeitfunktionen abbilden. Mit ihnen lassen sich Gedächtniseffekte modellieren, wie sie zum Beispiel in der Elastoplastizität oder im Ferromagnetismus auftauchen. Zunächst führen wir Hysterese-Operatoren allgemein ein und analysieren dann eine spezielle Klasse von Hysterese-Operatoren, die Prandtl-Ishlinskii-Operatoren. Wir untersuchen ihre Gedächtnisstruktur und erklären, wie sich die Operatoren numerisch auswerten lassen. Dazu stellen wir zwei verschiedene Approximationsansätze vor. Wir führen aus, wie sich die approximierenden Operatoren implementieren lassen und leiten lineare und quadratische Fehlerabschätzungen her. Zur numerischen Lösung des gekoppelten Systems aus der Wellengleichung mit einem Hysterese-Operator führen wir ein implizites Differenzenschema mit Gedächtnis ein. Für eine Klasse von Hysterese-Operatoren zeigen wir die Existenz und Eindeutigkeit der Lösung des numerischen Schemas, beweisen mit Hilfe von Kompaktheitsschlüssen und einem Monotonieargument die Konvergenz des Verfahrens und leiten eine Fehlerabschätzung der Ordnung 1/2 her. Wir diskutieren, wie das vorgestellte Verfahren auf die Prandtl-Ishlinskii-Operatoren angewendet werden kann. / In this thesis we develop and investigate a numerical scheme for the one-dimensional wave equation with hysteresis for different kinds of boundary conditions. This equation can be regarded as a model for the longitudinal or torsional oscillations of a homogeneous bar under the influence of an uniaxial external force density assuming an elastoplastic material law. Hysteresis operators are rate-independent Volterra operators mapping time functions to time functions. This kind of operator can be used to model memory effects as they appear in elastoplasticity or ferromagnetism, for example. We first give an introduction to the general concept of hysteresis operators before we analyze a special class of hysteresis operators called Prandtl-Ishlinskii operators. We investigate their memory structure and explain how the operators can be evaluated numerically. To that end we present two different kinds of approximation schemes. We point out how the approximating operators can be implemented and we derive linear and quadratic error estimates. For the numerical solution of the coupled system of the wave equation with a hysteresis operator we introduce an implicit difference scheme with memory. For a class of hysteresis operators we show the existence and uniqueness of the numerical solution. We prove the convergence of the scheme by compactness and monotonicity arguments. We derive an error estimate of order 1/2. We discuss the application of the method presented to Prandtl-Ishlinskii operators.

Ziel der Arbeit ist es, auf der Basis von Messungen ein rheologisches Materialmodell für technische Gummiwerkstoffe zu erstellen, welches deren Eigenschaften nachbildet, insbesondere vorhandene komplexe Zusammenhänge zwischen Relaxation, Erholung, Versuchsgeschwindigkeit und Belastungsamplitude. Dabei wird sich auf die Simulation von großen einfachen, aber beliebigen Scherverformungen beschränkt, woraus ein skalarwertiges Modell resultiert. Anwendung finden generalisierte Maxwell-Elemente und generalisierte kontinuierliche Prandtl-Elemente. Verschiedene Modellvarianten werden diskutiert. Es wird ein Berechnungsprogramm unter MATLAB erstellt.