The presented dissertation thesis deals with a simulation of the human phonation in terms of latest theories. Phonation is considered here as a bi-directional fluid-structure-acoustic interaction, where the interaction between all three physical domains occurs due to the unsteady viscous compressible Navier-Stokes equations. There is a solid knowledge background in the first part of the thesis. It concerns the latest concepts in computational modeling of the human phonation, the most important and recent theories about the human voice production and some key aspects of the human anatomy, physiology and pathology. Also voice assessment is discussed. The second part of the thesis describes an in-depth analysis of a phonation simulation in a planar computational model. The basic concepts proceed from algorithms developed in the Institute of Solid Mechanics, Mechatronics and Biomechanics. Created models are able to reproduce sounds of all Czech vowels and the most common evaluated parameters very close to physiological ranges. The simulated pathology, Reinke's edema, is demonstrated in order to explore its influence on the vowel sound. The third part focuses on modeling of phonation in a spatial computational model. All Czech vowels are simulated also here and compared to the planar model and to actual measurement. The spatial model serves as the starting point to modeling of a longitudinal pretension incorporated in the vocal folds. In the last part of the thesis, a modeling of the phonation with vocal folds pretension is investigated. Although the models are tuned to a rather soft phonation, the results are in agreement with the relevant physiologic phenomena. While the spatial model is highly computationally expensive, a hybrid planar model with pretension is proposed. A special attention is paid to the analysis of self-sustained oscillation of the vocal folds. It is shown, the planar model cannot reproduce such kind of oscillation in the actual version, albeit time of oscillation was considerably extended. On the other hand, oscillation of the spatial vocal folds are stabilized without effects accompanying subduing of oscillation. It can be supposed that the spatial model is able to reproduce self-sustained oscillation as a basic principle present during the human phonation.

No description available.

Raman spectroscopy, laser-induced fluorescence, and sum-frequency generation (SFG) spectroscopy are used to investigate the behavior of compressible mobile phases and stationary phases under a variety of chromatographic conditions. Efforts to understand and optimize separations employing compressible mobile phases have been limited by a lack of understanding of the mobile phase density gradient. Mobile phase compressibility leads to gradients in linear velocity and solute retention and affects separation speed and efficiency, especially in packed columns. This work describes on-column density measurement of CO2, a common carrier fluid for SFC and SGC, in packed capillary columns using Raman microspectroscopy. On-column detection by laser-induced fluorescence is used to observe the effect of the mobile phase density gradient on separation speed and efficiency, and experimental efficiency is compared to a theoretical model. Additionally, SFG spectroscopy allows for probing the structure of model monomeric and polymeric C18 stationary phases under pressure; this provides a basis for correlating selectivity with pressure-induced structural changes in stationary phase materials. Together, this work provides a more complete understanding of the role of column pressure and fluid compressibility on the speed, efficiency, and selectivity of chemical separations.

<p>  </p> <p>Acoustic streaming is a phenomenon which occurs when acoustic excitations interact with a fluid (stationary or non-stationary). Exploitation of this phenomenon has the potential to open doors to new methods of flow control through the enhancement or diminishment of the present flow instabilities. A particular use of acoustic streaming shown by previous numerical studies is the enhancement of heat transfer in violation of the Reynold’s Analogy within a small range of Mach numbers and frequencies of periodic excitation. The focus of this thesis is to experimentally assess the usage of a Rossiter cavity in generating periodic acoustic excitations and its effects on the shear stress and heat transfer. </p> <p>In the present research, two large models are tested using a blow-down facility. The models are made of aluminum and Teflon and were developed to ensure optical access for infrared thermography. The geometries are tested at Mach number ranging from 0.373 to 0. 866. The target Mach number-frequency pair where significant heat transfer enhancement is a free stream Mach number at the cavity, Mc, of 0.75 and the frequency, fc, of 7.5 kHz. The cavity is tuned using the Rossiter equation with Rossiter constants k = 0.66 and y = 0.25. The heat transfer and skin friction enhancement are measured immediately upstream and downstream of the cavity and compared to the previous numerical studies.</p> <p>When testing the Teflon model with an ambient back pressure and 11 lb/s mass flow, a frequency of 7.8 kHz was generated by the cavity. For the aluminum model tested at a high vacuum and 3 lb/s mass flow, frequencies near 7, 10, and 20 kHz were generated by the cavity with 10 and 20 kHz appearing most often. High speed schlieren imaging was used to confirm the flow structures being generated in the flow. There was good agreement with the Rossiter modes at lower Mach numbers and moderate agreement at transonic Mach numbers. A correlation is presented which defines a band of Mach number-Reynolds number pairs which present with a discontinuous frequency behavior during operation of the wind tunnel. Measurable effects on both skin friction and heat transfer between tests with comparable operating conditions to a reference were observed and are presented.</p>

The use of turbocharged Diesel engines is nowadays a widespread practice in the automotive sector: heavy-duty vehicles like trucks or buses, in particular, are often equipped with turbocharged engines. An accurate study of the flow field developing inside both the main components of a turbocharger, i.e. compressor and turbine, is therefore necessary: the synergistic use of CFD simulations and experimental tests allows to fulfill this requirement. The aim of this thesis is to investigate the performance and the flow field that develops inside a centrifugal compressor for automotive turbochargers. The study is carried out by means of numerical simulations, both steady-state and transient, based on RANS models (Reynolds Averaged Navier-Stokes equations). The code utilized for the numerical simulations is Ansys CFX.   The first part of the work is an engineering attempt to develop a CFD method for predicting the performance of a centrifugal compressor which is based solely on steady-state RANS models. The results obtained are then compared with experimental observations. The study continues with an analysis of the sensitivity of the developed CFD method to different parameters: influence of both position and model used for the rotor-stator interfaces and the axial tip-clearance on the global performances is studied and quantified.   In the second part, a design optimization study based on the Design of Experiments (DoE) approach is performed. In detail, transient RANS simulations are used to identify which geometry of the recirculation cavity hollowed inside the compressor shroud (ported shroud design) allows to mitigate the backflow that appears at low mass-flow rates. Backflow can be observed when the operational point of the compressor is suddenly moved from design to surge conditions. On actual heavy-duty vehicles, these conditions may arise when a rapid gear shift is performed. / Användningen av turboladdade dieselmotorer ärr numera utbredd inom bilindustrin: i synnerhet tunga fordon som lastbilar eller bussar ärr ofta utrustade med turbo-laddade motorer. En utförlig förståelse av flödesfältet som utvecklas innuti båda huvudkomponenterna hos en turboladdare, dvs kompressor och turbin, är därför nödvändig: den synergistiska användningen av CFD-simuleringar och experimentel-la tester möjliggör att detta krav uppfylls. Syftet med denna avhandling är att undersöka prestanda och det flödesfält som utvecklas i en centrifugalkompressor för turboladdare. Studien utförs genom nu-meriska simuleringar, både steady state och transient, baserat på RANS-modeller (Reynolds Averaged Navier-Stokes-ekvationer). Koden som används för de numeriska simuleringarna är Ansys CFX.   Den första delen av arbetet ¨ar ett försöka att utveckla en CFD-metod för att förutsäga prestanda för en centrifugalkompressor med hjälp av steady-state RANS-modeller. De erhållna resultaten jämförs sedan med experimentella observationer. Studien fortsätter med en analys av känsligheten hos den utvecklade CFD-metoden till olika parametrar: Inflytande av både position och modell som används för rotor-statorgränssnitt samt axiellt spel mellan rotor och hus på de globala prestationerna studeras och kvantifieras.   I andra delen utförs en designoptimeringsstudie baserad på Design of Experiments (DoE). I detalj används tidsupplösta RANS-simuleringar för att identifiera vilken utformning av ported shroud som minskar backflöde i kompressorn under en snabb minskning av massflöde och varvtal och därmed ger bättre prestanda i transient surge. På tunga fordon kan dessa förhållanden uppstå under växling.

Current-to-pressure (I/P) converters are pneumatic devices which provide precise control of pressure in various industries – for example these devices are often used in valve positioner systems (typically found in the oil and gas industry) and tensioning systems (typically used in the packaging industry). With an increasing demand for such devices to operate in harsh environments all by delivering acceptable performance means that Current-to-pressure converters need to be carefully designed such that environmental factors have no or minimal effects on its performance. This work presents an investigation of the principles of operation of an existing I/P converter through mathematical modelling. A simulation model has been created and which allows prediction of performance of the I/P converter. This tool has been used to identify areas of poor performances through theoretical analysis and consequently led to optimisation of certain areas of the I/P converter through a design change to deliver improved performances, for instance the average percentage shift in gain at 1mA input signal (over the temperature range of -40°C to 85°C) on the new I/P converter is 2.13% compared to the average gain of 4.24% on the existing I/P converter, which represents an improvement of almost two fold. Experimental tests on prototypes have been carried out and tests results showed that improved linearity and temperature sensitivity can be expected from the new design.

Computational fluid dynamics (CFD) has become a significant tool routinely used in design and optimization in aerospace industry. Typical flows may be characterized by high-speed and compressible flow features and, in many cases, by massive flow separation and unsteadiness. Accurate and efficient numerical solution of time-dependent problems is hence required, and the efficiency of standard dual-time stepping methods used for unsteady flows in many CFD codes has been found inadequate for large-scale industrial problems. This has motivated the present work, in which major effort is made to replace the explicit relaxation methods with implicit time integration schemes. The CFD flow solver considered in this work is Edge, a node-based solver for unstructured grids based on a dual, edge-based formulation. Edge is the Swedish national CFD tool for computing compressible flow, used at the Swedish aircraft manufacturer SAAB, and developed at FOI, lately in collaboration with external national and international partners. The work is initially devoted to the implementation of an implicit Lower-Upper Symmetric Gauss-Seidel (LU-SGS) type of relaxation in Edge with the purpose to speed up the convergence to steady state. The convergence of LU-SGS has been firstly accelerated by basing the implicit operator on a flux splitting method of matrix dissipation type. An increase of the diagonal dominance of the system matrix was the principal motivation. Then the code has been optimized by means of performance tools Intel Vtune and CrayPAT, improving the run time. It was found that the ordering of the unknowns significantly influences the convergence. Thus, different ordering techniques have been investigated. Finding the optimal ordering method is a very hard problem and the results obtained are mostly illustrative. Finally, to improve convergence speed on the stretched computational grids used for boundary layers LU-SGS has been combined with the line-implicit method. / QC 20120626

Numerical methods for solving partial differential equations are commonplace in the engineering community and their popularity can be attributed to the rapid performance improvement of modern workstations and desktop computers. The ubiquity of computer technology has allowed all areas of engineering to have access to detailed thermal, stress, and fluid flow analysis packages capable of performing complex studies of current and future designs. The rapid pace of computer development, however, has begun to outstrip efforts to reduce analysis overhead. As such, most commercially available software packages are now limited by the human effort required to prepare, develop, and initialize the necessary computational models. Primarily due to the mesh-based analysis methods utilized in these software packages, the dependence on model preparation greatly limits the accessibility of these analysis tools. In response, the so-called meshless or mesh-free methods have seen considerable interest as they promise to greatly reduce the necessary human interaction during model setup. However, despite the success of these methods in areas demanding high degrees of model adaptability (such as crack growth, multi-phase flow, and solid friction), meshless methods have yet to gain notoriety as a viable alternative to more traditional solution approaches in general solution domains. Although this may be due (at least in part) to the relative youth of the techniques, another potential cause is the lack of focus on developing robust methodologies. The failure to approach development from a practical perspective has prevented researchers from obtaining commercially relevant meshless methodologies which reach the full potential of the approach. The primary goal of this research is to present a novel meshless approach called MIMS (Model Integrated Meshless Solver) which establishes the method as a generalized solution technique capable of competing with more traditional PDE methodologies (such as the finite element and finite volume methods). This was accomplished by developing a robust meshless technique as well as a comprehensive model generation procedure. By closely integrating the model generation process into the overall solution methodology, the presented techniques are able to fully exploit the strengths of the meshless approach to achieve levels of automation, stability, and accuracy currently unseen in the area of engineering analysis. Specifically, MIMS implements a blended meshless solution approach which utilizes a variety of shape functions to obtain a stable and accurate iteration process. This solution approach is then integrated with a newly developed, highly adaptive model generation process which employs a quaternary triangular surface discretization for the boundary, a binary-subdivision discretization for the interior, and a unique shadow layer discretization for near-boundary regions. Together, these discretization techniques are able to achieve directionally independent, automatic refinement of the underlying model, allowing the method to generate accurate solutions without need for intermediate human involvement. In addition, by coupling the model generation with the solution process, the presented method is able to address the issue of ill-constructed geometric input (small features, poorly formed faces, etc.) to provide an intuitive, yet powerful approach to solving modern engineering analysis problems.

We analyze a number of systems of evolution equations that arise in the study of physical chemistry. First we discuss the well-posedness of a system of mixing compressible barotropic multicomponent flows. We discuss the regularity of these variational solutions, their existence and uniqueness, and we analyze the emergence of a novel type of entropy that is derived for the system of equations. Next we present a numerical scheme, in the form of a discontinuous Galerkin (DG) finite element method, to model this compressible barotropic multifluid. We find that the DG method provides stable and accurate solutions to our system, and that further, these solutions are energy consistent; which is to say that they satisfy the classical entropy of the system in addition to an additional integral inequality. We discuss the initial-boundary problem and the existence of weak entropy at the boundaries. Next we extend these results to include more complicated transport properties (i.e. mass diffusion), where exotic acoustic and chemical inlets are explicitly shown. We continue by developing a mixed method discontinuous Galerkin finite element method to model quantum hydrodynamic fluids, which emerge in the study of chemical and molecular dynamics. These solutions are solved in the conservation form, or Eulerian frame, and show a notable scale invariance which makes them particularly attractive for high dimensional calculations. Finally we implement a wide class of chemical reactors using an adapted discontinuous Galerkin finite element scheme, where reaction terms are analytically integrated locally in time. We show that these solutions, both in stationary and in flow reactors, show remarkable stability, accuracy and consistency. / text

In this work, a DNS code was developed to investigate problems on transition in subsonic compressible boundary layer on a flat plate. Code validation tests were performed for linear and nonlinear stages of transition, on incompressible and compressible regimes. The focus of the present work is to investigate natural transition in subsonic boundary layers modeled by wave packets; and perform a preliminary study of transition induced by white noise. Three main problems were considered, namely, a DNS simulation and analysis of the ex- periment (MEDEIROS; GASTER, 1999b) of wave packet evolution on incompressible boundary layer, the influence of compressibility on wave packet evolution at subsonic Mach numbers and finally, a preliminary study of the evolution of a white noise perturbation in the boundary layer at Mach 0.2 and Mach 0.9. Comparisons between numerical and experimental results show remarkably good agreement in the linear and nonlinear stages, in both, spatial and Fourier spaces. A numerical simulation of this experiment and the analysis carried out is not available in the literature for wave packets in the incompressible boundary layer. The nonlinear modal analysis performed established the existence of tuned fundamental and subharmonic resonance of H-type and K-type in the packet. Influence of compressibility in the wave packet evolution was here investigated in boundary layers at Mach 0.7 and Mach 0.9. There are no works reported in the literature on wave packets in compressible subsonic boundary layer. In the linear regime, the oblique modes were the most unstable for Mach > 0.7, as expected by the results of the literature. In the nonlinear regime, strong streaks were observed, associated with low frequency modes that eventually decay downstream. An isolated wave packet at Mach 0.9 showed nonlinear amplification only in the subharmonic band, which may be associated to H-type or detuned resonance. However this packet has a relatively stable character. On the other hand, at Mach 0.9 spanwise interaction of wave packet pairs were more unstable than the isolated case, because stable modes for the isolated packet evolution becomes unstable in the wave packet interaction. This scenario evidenced the presence of oblique transition. Finally, the nonlinear evolution of the same white noise disturbance at Mach 0.2 and Mach 0.9 were observed to be completely different. In the incompressible boundary layer localized lambda vortex structures were observed, that could be associated to the local presence of H-type and/or K-type resonance. In the compressible regime, longitudinal vortex structures distributed across the entire domain seemed to be linked to oblique transition. In the white noise evolution, compressibility seems to have a stronger effect than in the wave packet evolution. In the conditions considered, the wave packet interaction appear to be a better representation of white noise compressible transition scenario. / No presente trabalho, um código DNS (Direct Numerical Simulation) foi desenvolvido para abordar problemas de transição para turbulência em camada limite subsônica compressível em uma placa plana. Foram realizados testes de validação de código , nos regimes linear e não linear do processo de transição, nos regimes incompressível e compressível. O foco do presente trabalho é estudar transição natural modelada por meio de pacotes de onda em camada limite compressível subsônica, e realizar uma análise preliminar da transição induzida por ruído branco. Três assuntos principais foram considerados: uma simulação DNS e uma análise comparativa com o experimento (MEDEIROS; GASTER, 1999b) sobre a evolução de um pacote de ondas em camada limite incompressível, a influência da compressibilidade na evolução de pacotes de ondas no regime subsônico, e por último, um estudo preliminar da transição induzida por ruído branco em Mach 0.2 e Mach 0.9. As comparações realizadas entre a solução numérica e os dados experimentais mostram uma boa concordância, nos regimes linear e não linear, tanto no espaço físico quanto no espaço de Fourier. A simulação numérica deste experimento e a análise realizada neste trabalho, não são encontradas na literatura para o regime incompressível. A análise modal não linear aplicada aos resultados, permitiu identificar a presença das ressonâncias tipo H e tipo K no pacote de ondas. A influência da compressibilidade na evolução dos pacotes de onda foi estudada em Mach 0.7 e Mach 0.9. Na literatura não há trabalhos sobre pacotes de ondas no regime sub- sônico. No regime linear da transição, os modos oblíquos resultam ser os mais instáveis para Mach > 0.7, como era de esperar, de acordo com os resultados da literatura. No regime não linear, foram observadas estrias de moderada amplitude, associadas com modos de baixa frequência que acabam decaindo. O pacote de ondas em Mach 0.9 apresentou amplificação não linear somente na banda subharmônica, que pode ser associada com transição tipo H ou ressonância dessintonizada. No entanto, o comportamento geral neste regime é estabilizante. Por sua vez, a interação entre pacotes de ondas em Mach 0.9 mostrou um comportamento desestabilizante, pois a interação acaba gerando amplificação não linear em modos que decaem no pacote isolado. Os modos amplificados sugerem a presença do mecanismo de transição oblíqua. Finalmente, a evolução da mesma perturbação constituída por ruído branco em Mach 0.2 e Mach 0.9, resultaram ser completamente diferentes. Na camada limite incompressível foram observados vórtices tipo lambda, que poderiam ser gerados pela presença localizada das ressonâncias tipo H e/ou tipo K. No regime compressível foram observados vórtices distribuidos em todo o domínio, o que sugere a presença da transição oblíqua. Na transição gerada por ruído branco a compressibilidade teve uma influência maior que no pacote de ondas. Nas condições estudadas, a interação entre pacotes de ondas parece ser uma melhor representação do ruído branco no regime compressível.