Air velocities induced by underground vehicles in complex metro systems are obtained using artificial neural networks. Complex tunnel shaft-systems with any number of tunnels and shafts and with most of the practically possible geometries encountered in underground structures can be simulated with the proposed method. A single neural network, of type feed-forward back propagation, with a single hidden layer is trained for modelling a single tunnel segment. Train and tunnel parameters that have influence on the vehicle induced flow characteristics are used together to obtain non-dimensional input and target parameters. First input parameter is the major head loss coefficient of tunnel, (L/D)_Tunnel. Blockage ratio A_Train/A_Tunnel and train aspect ratio (D/L)_Train are selected to be non-dimensional input parameters to represent the system geometry. As the final input parameter, skin friction coefficient of the train, f_Train drag coefficient of the train, C_D / frontal area of the train, A_Train and lateral area of the train, A_Lateral are combined into a single overall drag coefficient based on the train frontal area. Non-dimensional V_Air/V_Train speed ratio is selected to be the target parameter. Using maximum air velocity predicted by the trained neural network together with non dimensional system parameters and time, an additional neural network is trained for predicting the deceleration of air in case of train stoppage within the tunnel system and departure of the train from the system. A simulation tool for predicting time dependent velocity profile of air in metro systems is developed with the trained neural networks.

A new version of the Deforming-Spatial-Domain/Stabilized Space--Time (DSD/SST) formulation of incompressible flows has been introduced recently to have additional subgrid-scale representation features. This is the space--time version of the residual-based variational multiscale (VMS) method. The new version is called DSD/SST-VMST (i.e. the version with the VMS turbulence model) and also Space--Time VMS (ST-VMS). The thesis starts with a brief overview of the ST-VMS method. It continues with a comprehensive set of test computations with 3D airfoil geometries at high Reynolds numbers and comparison with experimental data. The thesis shows that the test computations signal a promising future for the ST-VMS method.

A new version of the Deforming-Spatial-Domain/Stabilized Space--Time (DSD/SST) formulation of incompressible flows has been introduced recently to have additional subgrid-scale representation features. This is the space--time version of the residual-based variational multiscale (VMS) method. The new version is called DSD/SST-VMST (i.e. the version with the VMS turbulence model) and also Space--Time VMS (ST-VMS). The thesis starts with a brief overview of the ST-VMS method. It continues with a comprehensive set of test computations with 3D airfoil geometries at high Reynolds numbers and comparison with experimental data. The thesis shows that the test computations signal a promising future for the ST-VMS method.

The present thesis is devoted both to the study of strictly hyperbolic operators with low regularity coefficients and of the density-dependent incompressible Euler system. On the one hand, we show a priori estimates for a second order strictly hyperbolic operator whose highest order coefficients satisfy a log-Zygmund continuity condition in time and a log-Lipschitz continuity condition with respect to space. Such an estimate involves a time increasing loss of derivatives. Nevertheless, this is enough to recover well-posedness for the associated Cauchy problem in the space $H^infty$ (for suitably smooth second order coefficients).In a first time, we consider acomplete operator in space dimension $1$, whose first order coefficients were assumed Hölder continuous and that of order $0$only bounded. Then, we deal with the general case of any space dimension, focusing on a homogeneous second order operator: the step to higher dimension requires a really different approach. On the other hand, we consider the density-dependent incompressible Euler system. We show its well-posedness in endpoint Besov spaces embedded in the class of globally Lipschitz functions, producing also lower bounds for the lifespan of the solution in terms of initial data only. This having been done, we prove persistence of geometric structures, such as striated and conormal regularity, for solutions to this system. In contrast with the classical case of constant density, even in dimension $2$ the vorticity is not transported by the velocity field. Hence, a priori one can expect to get only local in time results. For the same reason, we also have to dismiss the vortex patch structure. Littlewood-Paley theory and paradifferential calculus allow us to handle these two different problems .A new version of paradifferential calculus, depending on a parameter $ggeq1$, is also needed in dealing with hyperbolic operators with nonregular coefficients. The general framework is that of Besov spaces, which includes in particular Sobolev and Hölder sets. Intermediate classes of functions, of logaritmic type, come into play as well

Aquesta tesi doctoral proposa formulacions de Galerkin discontinu (DG) d'alt ordre per fluxos viscosos incompressibles. Es desenvolupa un nou mètode de DG amb penalti interior (IPM-DG), que condueix a una forma feble simètrica i coerciva pel terme de difusió, i que permet assolir una aproximació espacial d'alt ordre. Aquest mètode s'aplica per resoldre les equacions de Stokes i Navier-Stokes. L'espai d'aproximació de la velocitat es descompon dins de cada element en una part solenoidal i una altra irrotacional, de manera que es pot dividir la forma dèbil IPM-DG en dos problemes desacoblats. El primer permet el càlcul de les velocitats i de les pressions híbrides, mentre que el segon calcula les pressions en l'interior dels elements. Aquest desacoblament permet una reducció important del número de graus de llibertat tant per velocitat com per pressió. S'introdueix també un paràmetre extra de penalti resultant en una formulació DG alternativa per calcular les velocitats solenoidales, on les pressions no apareixen. Les pressions es poden calcular com un post-procés de la solució de les velocitats. Es contemplen altres formulacions DG, com per exemple el mètode Compact Discontinuous Galerkin, i es comparen al mètode IPM-DG. Es proposen mètodes implícits de Runge-Kutta d'alt ordre per problemes transitoris incompressibles, permetent obtenir esquemes incondicionalment estables i amb alt ordre de precisió temporal. Les equacions de Navier-Stokes incompressibles transitòries s'interpreten com un sistema de Equacions Algebraiques Diferencials, és a dir, un sistema d'equacions diferencials ordinàries corresponent a la equació de conservació del moment, més les restriccions algebraiques corresponent a la condició d'incompressibilitat. Mitjançant exemples numèrics es mostra l'aplicabilitat de les metodologies proposades i es comparen la seva eficiència i precisió. / This PhD thesis proposes divergence-free Discontinuous Galerkin formulations providing high orders of accuracy for incompressible viscous flows. A new Interior Penalty Discontinuous Galerkin (IPM-DG) formulation is developed, leading to a symmetric and coercive bilinear weak form for the diffusion term, and achieving high-order spatial approximations. It is applied to the solution of the Stokes and Navier-Stokes equations. The velocity approximation space is decomposed in every element into a solenoidal part and an irrotational part. This allows to split the IPM weak form in two uncoupled problems. The first one solves for velocity and hybrid pressure, and the second one allows the evaluation of pressures in the interior of the elements. This results in an important reduction of the total number of degrees of freedom for both velocity and pressure. The introduction of an extra penalty parameter leads to an alternative DG formulation for the computation of solenoidal velocities with no presence of pressure terms. Pressure can then be computed as a post-process of the velocity solution. Other DG formulations, such as the Compact Discontinuous Galerkin method, are contemplated and compared to IPM-DG. High-order Implicit Runge-Kutta methods are then proposed to solve transient incompressible problems, allowing to obtain unconditionally stable schemes with high orders of accuracy in time. For this purpose, the unsteady incompressible Navier-Stokes equations are interpreted as a system of Differential Algebraic Equations, that is, a system of ordinary differential equations corresponding to the conservation of momentum equation, plus algebraic constraints corresponding to the incompressibility condition. Numerical examples demonstrate the applicability of the proposed methodologies and compare their efficiency and accuracy.

Time integration methods are necessary for the solution of transient flow problems. In recent years, interest in transient flow problems has increased, leading to a need for better understanding of the costs and benefits of various time integration schemes. The present work investigates two common time integration schemes, namely the Semi-Implicit Method for Pressure-Linked Equations (SIMPLE) and the Fractional Step (FS) method. Three two-dimensional, transient, incompressible flow problems are solved using a cell centered, finite volume code. The three test cases are laminar flow in a lid-driven skewed cavity, laminar flow over a square cylinder, and turbulent flow over a square cylinder. Turbulence is modeled using wall functions and the k - ε turbulence model with the modifications suggested by Kato and Launder. Solution efficiency as measured by the effort carried out by the flow equation solver and CPU time is examined. Accuracy of the results, generated using the SIMPLE and FS time integration schemes, is analyzed through a comparison of the results with existing experimental and/or numerical solutions. Both the SIMPLE and FS algorithms are shown to be capable of solving benchmark flow problems with reasonable accuracy. The two schemes differ slightly in their prediction of flow evolution over time, especially when simulating very slowly changing flows. As the time step size decreases, the SIMPLE algorithm computational cost (CPU time) per time step remains approximately constant, while the FS method experiences a reduction in cost per time step. Also, the SIMPLE algorithm is numerically stable for time steps approaching infinity, while the FS scheme suffers from numerical instability if the time step size is too large. As a result, the SIMPLE algorithm is recommended to be used for transient simulations with large time steps or steady state problems while the FS scheme is better suited for small time step solutions, although both time-stepping schemes are found to be most efficient when their time steps are at their maximum stable value.

This thesis presents a GPU accelerated implementation of a high order splitting scheme with a spectral element discretization for the incompressible Navier Stokes (INS) equations. While others have implemented this scheme on clusters of processors using the Nek5000 code, to my knowledge this thesis is the first to explore its performance on the GPU. This work implements several of the Nek5000 algorithms using OpenCL kernels that efficiently utilize the GPU memory architecture, and achieve massively parallel on chip computations. These rapid computations have the potential to significantly enhance computational fluid dynamics (CFD) simulations that arise in areas such as weather modeling or aircraft design procedures. I present convergence results for several test cases including channel, shear, Kovasznay, and lid-driven cavity flow problems, which achieve the proven convergence results.

The objectives in this thesis are to present a model of a variable venturi in an exhaust gas recirculation (EGR) system located in a heavy duty diesel engine. A new legislation called EURO~4 will come into force in 2005 which affects truck development and it will require an On-Board Diagnostic system in the truck. If model based diagnostic systems are to be used, one of the advantages is that the system performance will increase if a model of a variable venturi is used. Three models with different complexity are compared in ten different experiments. The experiments are performed in a steady flow rig at different percentage of EGR gases and venturi areas. The model predicts the mass flow through the venturi. The results show that the first model with fewer simplifications performs better and has fewer errors than the other two models. The simplifications that differ between the models are initial velocity before the venturi and the assumption of incompressible flow. The model that shows the best result is not proposed by known literature in this area of knowledge and technology. This thesis shows that further studies and work on this model, the model with fewer simplifications, can be advantageous.

Time integration methods are necessary for the solution of transient flow problems. In recent years, interest in transient flow problems has increased, leading to a need for better understanding of the costs and benefits of various time integration schemes. The present work investigates two common time integration schemes, namely the Semi-Implicit Method for Pressure-Linked Equations (SIMPLE) and the Fractional Step (FS) method. Three two-dimensional, transient, incompressible flow problems are solved using a cell centered, finite volume code. The three test cases are laminar flow in a lid-driven skewed cavity, laminar flow over a square cylinder, and turbulent flow over a square cylinder. Turbulence is modeled using wall functions and the k - ε turbulence model with the modifications suggested by Kato and Launder. Solution efficiency as measured by the effort carried out by the flow equation solver and CPU time is examined. Accuracy of the results, generated using the SIMPLE and FS time integration schemes, is analyzed through a comparison of the results with existing experimental and/or numerical solutions. Both the SIMPLE and FS algorithms are shown to be capable of solving benchmark flow problems with reasonable accuracy. The two schemes differ slightly in their prediction of flow evolution over time, especially when simulating very slowly changing flows. As the time step size decreases, the SIMPLE algorithm computational cost (CPU time) per time step remains approximately constant, while the FS method experiences a reduction in cost per time step. Also, the SIMPLE algorithm is numerically stable for time steps approaching infinity, while the FS scheme suffers from numerical instability if the time step size is too large. As a result, the SIMPLE algorithm is recommended to be used for transient simulations with large time steps or steady state problems while the FS scheme is better suited for small time step solutions, although both time-stepping schemes are found to be most efficient when their time steps are at their maximum stable value.

<p>The objectives in this thesis are to present a model of a variable venturi in an exhaust gas recirculation (EGR) system located in a heavy duty diesel engine. A new legislation called EURO~4 will come into force in 2005 which affects truck development and it will require an On-Board Diagnostic system in the truck. If model based diagnostic systems are to be used, one of the advantages is that the system performance will increase if a model of a variable venturi is used. </p><p>Three models with different complexity are compared in ten different experiments. The experiments are performed in a steady flow rig at different percentage of EGR gases and venturi areas. The model predicts the mass flow through the venturi. The results show that the first model with fewer simplifications performs better and has fewer errors than the other two models. The simplifications that differ between the models are initial velocity before the venturi and the assumption of incompressible flow. </p><p>The model that shows the best result is not proposed by known literature in this area of knowledge and technology. This thesis shows that further studies and work on this model, the model with fewer simplifications, can be advantageous.</p>