Em experimentação agrícola, é freqüente a necessidade de análise conjunta de grupos de experimentos. Em muitos casos, o pesquisador deseja generalizar resultados para condições gerais de regiões e/ou em avaliar o desempenho de vários genótipos (tratamentos) em diversos ambientes (locais e/ou ano). Quando um conjunto de experimentos é planejado para vários locais é necessário considerar o delineamento individual em cada local e a combinação total dos genótipos com os locais (interação genótipo × ambiente). Logo, os dados observados podem ser organizados em uma tabela de dupla entrada. Existem várias metodologias de análise e interpretação para a interação genótipo × ambiente proveniente de um grupo de cultivares testados em vários ambientes. Entre essas metodologias destaca-se os modelos AMMI ("additive main effects and multiplicative interaction model"), como o próprio nome diz é um método uni-multivariado, que engloba uma análise de variância para os efeitos principais, que são os efeitos dos genótipos (G) e os ambientes (E) e para efeitos multiplicativos (interação genótipo × ambiente), utiliza-se a decomposição em valor singular (DVS). Essa técnica multivariada baseia-se no uso dos autovalores e autovetores provenientes da matriz de interação genótipo × ambiente. Araújo e Dias (2005) verificaram o problema de superestimação e subestimação de autovalores estimados da maneira convencional. Para superar esses problemas de estimação de autovalores, Muirhead (1987) apresenta três métodos para corrigir autovalores estimados a partir das matrizes de covariâncias amostral e alerta que nem sempre essas correções mantêm a ordem decrescente de valores, assim é sugerido o uso de regressão isotônica para ordenar esses dados, mas propriamente um algoritmo apresentado por Lin e Pearlman (1985). Os resultados indicaram que: A regressão isotônica juntamente com o algoritmo foi necessária e se mostrou muito importante em todos conjuntos de dados; Houve uma redução no número de componentes significativos para serem retidos nos modelos, fazendo com que os modelos AMMI selecionados sejam mais parcimoniosos quando se utiliza qualquer um dos métodos de correção; O método 2 apresentou as menores taxa de correção da soma de quadrados da interação genótipo × ambiente e o método 3 apresentou a maiores taxa de correção; Em relação a medida RMSPDPRESS, os menores valores foram obtidos quando se utilizou o método de correção 2. Já o método de correção 3 apresentou os maiores valores para RMSPDPRESS; O método 2 também se mostrou melhor quando o interesse era verificar o ganho em número de repetições, sendo que este benefício esteve sempre próximo de 3 repetições. Já o método 3 é o que apresenta um menor ganho em número de repetições, em torno de duas repetições. / In agricultural research is common to analyse groups of experiments. In many cases, the researcher intends to generalize results to general conditions of areas and/or evaluate the responses of several genotypes (treatments) in several environments (places and/or years). When a group of experiments is planned for several places it is necessary to consider the of design in each place and the combinations of the genotypes with the places (the interaction of genotype × environment). The observed data can be organized in an array. There are several methods of analysis and interpretation for the genotype × environment interaction from a group of genotype tested in several environments. These methods include AMMI models ("additive main effect and multiplicative interaction models"). As the name says it is a uni-multivariate method, that includes an analysis of variance for the main effects (the effects of the genotypes (G) and environments (E)) and assumes multiplicative effects for the genotype × environment interaction, using a singular value decomposition (DVS). This method estimates the eigenvalues and eigenvectors deriving from the matrix of genotype × environment interaction. Araújo and Dias (2005) found an overestimation and underestimation problem with the eigenvalues in the conventional way. To correct these problems Muirhead (1987) presents three methods to correct the eigenvalues from covariance the matrix and noted that these do not always maintain the order of values. The author suggested the use of isotonic regression to correct the eigenvalues, using an algorithm presented by Lin and Pearlman (1985). The results indicated that: The isotonic regression with the algorithm is necessary and it showed very important in all groups of data; There was a reduction in the number of significant components to be kept in the models and the order that the AMMI model selected is more parsimonious when any of the correction methods is used; The method 2 has the smallest rate of correction to the sum of squares of the genotype × environment interaction and method 3 has the largest correction rate; The measure RMSPDPRESS was smallest when method of correction 2 was used. The method of correction 3 has the largest values for RMSPDPRESS; Method 2 was also better when the interest was to verify the gain in number of replicates, and this benefit was always close to 3 replicates. The method 3 gives the smaller gain in the number of replicates, of around two replicates.

O som gerado por escoamentos, também conhecido como aeroacústica, tem se tornado cada vez mais importante em áreas industriais diversas desde aviação comercial até aparelhos eletrodomésticos, afetando diretamente os requisitos necessários para o desenvolvimento de novos produtos. Um caso particular é o ruído gerado por válvulas de compressores herméticos de refrigeração, sendo o compressor a principal fonte de ruído em refrigeradores domésticos. O presente trabalho tem por objetivo iniciar o desenvolvimento de uma ferramenta confiável de simulação capaz de auxiliar engenheiros na predição de problemas de aeroacústica, especialmente um que possa no futuro ser utilizado para estudar o ruído gerado pelo escoamento em válvulas de compressores herméticos. Para isso, foi desenvolvido um código para simulação numérica direta de aeroacústica. Utilizou-se processamento paralelo com decomposição de domínio para usar Simulação Numérica Direta em um tempo factível; esquemas de discretização espaciais e temporais de alta ordem para minimizar ao máximo os fenômenos de dissipação e dispersão do escoamento e das ondas acústicas e uma série de tratamentos no domínio como filtragem e estiramento da malha como também condições de contorno características com o intuito de obter uma solução adequada para estudo de aeroacústica. Assim, são apresentadas todas as etapas desenvolvidas no equacionamento, implementação e verificação. A verificação foi realizada segundo um processo matemático formal (Método das Soluções Manufaturadas) com o qual obteve-se que a ordem de precisão dos cálculos era a mesma da ordem formal dos esquemas de discretização utilizados para todas as variáveis. Também obteve-se a mesma concordância para análise do divergente da velocidade, verificando o código para simulação numérica direta de aeroacústica. Posteriormente, foram realizadas simulações de escoamentos compressíveis cisalhantes e seus resultados comparados com dados apresentados em literatura. Também foram calculadas as taxas de amplificação de perturbações e comparadas com a Teoria de Estabilidade Linear. Novamente, foram obtidos resultados satisfatórios nessas etapas, mostrando que a implementação do código DNS está verificada. / Sound generated by flow, known as aeroacoustics, is becoming more important in several industrial areas from commercial aircraft to household appliances, affecting directly the requirements to the development of new products. A particular case is the noise generated by valves of refrigeration hermetic compressors, being the compressor the main noise source in household refrigerators. This work has the goal of initiate the development of a reliable tool able to help engineers to predict aeroacoustics problems, specially one that can be used in the future to study the noise generated by the flow in valves of hermetic compressors. To do so, it was developed a numerical code to perform direct numerical simulation of aeroacoustics. It was used parallel processing with domain decomposition to use Direct Numerical Simulation in a feasible time; high order temporal and spatial discretization schemes to minimize the most the dispersion and dissipation phenomena of the flow field and of the acoustics waves and a series of treatments in the domain as filtering and mesh stretching as well as characteristics boundary conditions aiming a proper solution to study aeroacoustics. Thus, here all the steps developed in the formulation, implementation and verification are presented . The verification was done according to a formal mathematical procedure (Method of Manufactured Solutions) with which was found that the precision order of the calculations was the same of the formal order of the used discretization schemes for all variables. The same agreement was also obtained to the analysis of the divergence of the velocity, verifying the code to direct numerical simulation of aeroacoustics. Posteriorly, it were simulated shear compressible flows and the results were compared to literature data. Also, it were calculated the amplification rates of the disturbances and compared to Linear Stability Theory. Once more, it was obtained satisfactory results in these steps, showing that the implementation of the DNS code is verified.

In this thesis novel feedback attitude control algorithms and attitude estimation algorithms are developed for a three-axis stabilised spacecraft attitude control system. The spacecraft models considered include a rigid-body spacecraft equipped with (i) external control torque devices, and (ii) a redundant reaction wheel configuration. The attitude sensor suite comprises a three-axis magnetometer and three-axis rate gyroscope assembly. The quaternion parameters (also called Euler symmetric parameters), which globally avoid singularities but are subject to a unity-norm constraint, are selected as the primary attitude coordinates. There are four novel contributions presented in this thesis. The first novel contribution is the development of a robust control strategy for spacecraft attitude tracking maneuvers, in the presence of dynamic model uncertainty in the spacecraft inertia matrix, actuator magnitude constraints, bounded persistent external disturbances, and state estimation error. The novel component of this algorithm is the incorporation of state estimation error into the stability analysis. The proposed control law contains a parameter which is dynamically adjusted to ensure global asymptotic stability of the overall closedloop system, in the presence of these specific system non-idealities. A stability proof is presented which is based on Lyapunov's direct method, in conjunction with Barbalat's lemma. The control design approach also ensures minimum angular path maneuvers, since the attitude quaternion parameters are not unique. The second novel contribution is the development of a robust direct adaptive control strategy for spacecraft attitude tracking maneuvers, in the presence of dynamic model uncertainty in the spacecraft inertia matrix. The novel aspect of this algorithm is the incorporation of a composite parameter update strategy, which ensures global exponential convergence of the closed-loop system. A stability proof is presented which is based on Lyapunov's direct method, in conjunction with Barbalat's lemma. The exponential convergence results provided by this control strategy require persistently exciting reference trajectory commands. The control design approach also ensures minimum angular path maneuvers. The third novel contribution is the development of an optimal control strategy for spacecraft attitude maneuvers, based on a rigid body spacecraft model including a redundant reaction wheel assembly. The novel component of this strategy is the proposal of a performance index which represents the total electrical energy consumed by the reaction wheel over the maneuver interval. Pontraygin's minimum principle is applied to formulate the necessary conditions for optimality, in which the control torques are subject to timevarying magnitude constraints. The presence of singular sub-arcs in the statespace and their associated singular controls are investigated using Kelley's necessary condition. The two-point boundary-value problem (TPBVP) is formulated using Pontrayagin's minimum principle. The fourth novel contribution is an attitude estimation algorithm which estimates the spacecraft attitude parameters and sensor bias parameters from three-axis magnetometer and three-axis rate gyroscope measurement data. The novel aspect of this algorithm is the assumption that the state filtering probability density function (PDF) is Gaussian distributed. This Gaussian PDF assumption is also applied to the magnetometer measurement model. Propagation of the filtering PDF between sensor measurements is performed using the Fokker-Planck equation, and Bayes theorem incorporates measurement update information. The use of direction cosine matrix elements as the attitude coordinates avoids any singularity issues associated with the measurement update and estimation error covariance representation.

Submitted by Laison Junio da Silva Furlan (laisonfurlan@gmail.com) on 2018-09-11T21:45:47Z No. of bitstreams: 1 dissertacao_final.pdf: 2611301 bytes, checksum: b4c51f4e16b1f3e612c4d3a044c777c5 (MD5) / Approved for entry into archive by Claudia Adriana Spindola null (claudia@fct.unesp.br) on 2018-09-12T11:25:26Z (GMT) No. of bitstreams: 1 furlan_ljs_me_prud.pdf: 2445202 bytes, checksum: 1f1f2699158710f217b33ed602d0f51e (MD5) / Approved for entry into archive by Claudia Adriana Spindola null (claudia@fct.unesp.br) on 2018-09-12T11:36:28Z (GMT) No. of bitstreams: 1 furlan_ljs_me_prud.pdf: 2445202 bytes, checksum: 1f1f2699158710f217b33ed602d0f51e (MD5) / Made available in DSpace on 2018-09-12T11:36:28Z (GMT). No. of bitstreams: 1 furlan_ljs_me_prud.pdf: 2445202 bytes, checksum: 1f1f2699158710f217b33ed602d0f51e (MD5) Previous issue date: 2018-08-02 / Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES) / Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP) / O presente trabalho investiga a transição laminar-turbulenta devido a ondas de Tollmien-Schlichting para o escoamento de Poiseuille incompressível, bidimensional, de um fluido viscoelástico, utilizando a equação constitutiva Giesekus. A Teoria de Estabilidade Linear e a Simulação Numérica Direta são utilizadas para verificar a estabilidade de escoamentos de fluidos viscoelásticos a perturbações não estacionárias. Na análise LST a equação de Orr-Sommerfeld é modificada para um fluido viscoelástico e resolvida pelo método da estimativa (Shooting). Enquanto que, na formulação DNS, as equações de Navier-Stokes, juntamente com a equação constitutiva Giesekus, são resolvidas utilizando métodos de diferenças finitas compactas de alta ordem. Com o objetivo de avaliar as curvas neutras de estabilidade e as taxas de amplificação, diferentes simulações numéricas são realizadas variando-se os parâmetros adimensionais no modelo Giesekus e comparando com o fluido Newtoniano. As técnicas LST e DNS mostraram-se ferramentas eficientes na análise espacial da estabilidade de escoamentos viscoelásticos do tipo Giesekus, permitindo uma melhor compreensão da influência dos parâmetros adimensionais desses escoamentos e contribuindo com resultados originais na verificação da estabilidade de escoamentos viscoelásticos utilizando o fluido Giesekus. / The present work investigates the laminar-turbulent transition due to Tollmien-Schlichting waves for the incompressible two-dimensional Poiseuille flow of a viscoelastic fluid, using the Giesekus constitutive equation. Linear Stability Theory and Direct Numerical Simulation are used to verify the stability of viscoelastic fluid flows to unsteady disturbances. In the LST analysis, the Orr-Sommerfeld equation is modified to a viscoelastic fluid and solved by Shooting method. Whereas, in the DNS formulation, the Navier-Stokes equations with the Giesekus constitutive equation are solved using high-order compact finite difference methods. In order to evaluate the neutral stability curves and the amplification rates, different numerical simulations are performed by varying the dimensionless parameters in the Giesekus model and their results are compared with the Newtonian fluid. The LST and DNS techniques proved to be efficient tools to the spatial stability analysis of viscoelastic fluid flows of the Giesekus type, allowing a better comprehension of the dimensionless parameters influence of those flows, contributing with originals results to verification of the viscoelastics fluid flows stability using Giesekus fluid. / FAPESP: 2017/11068-6

Submitted by Ellen Silva Gervazoni null (emillen_@hotmail.com) on 2016-08-25T17:45:32Z No. of bitstreams: 1 Disserataçaofinal_EllenGervazoni.pdf: 1426674 bytes, checksum: 25ce3b189028e33485b050b9f7676a94 (MD5) / Approved for entry into archive by Juliano Benedito Ferreira (julianoferreira@reitoria.unesp.br) on 2016-08-29T19:17:40Z (GMT) No. of bitstreams: 1 gervazoni_es_me_prud.pdf: 1426674 bytes, checksum: 25ce3b189028e33485b050b9f7676a94 (MD5) / Made available in DSpace on 2016-08-29T19:17:40Z (GMT). No. of bitstreams: 1 gervazoni_es_me_prud.pdf: 1426674 bytes, checksum: 25ce3b189028e33485b050b9f7676a94 (MD5) Previous issue date: 2016-06-27 / Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES) / Diversos escoamentos de interesse prático são de fluidos viscoelásticos e muitas vezes é desejável saber se estes escoamentos propagam-se no estado laminar ou no turbulento. Embora a hidrodinâmica de fluidos viscoelásticos sejam fortemente afetadas pelo balanço entre forças inerciais e elásticas no escoamento, o efeito da elasticidade sobre a estabilidade de escoamentos inerciais não foi completamente estabelecida. No presente trabalho, estuda-se o que ocorre entre estes dois estados, na transição laminar-turbulenta. Especi- ficamente, é investigada a convecção de ondas de Tollmien-Schlichting para o escoamento incompressível de Poiseuille para um fluido viscoelástico, utilizando a equação constitutiva Oldroyd-B. Para isto, utiliza-se a Simulação Numérica Direta para verificar a estabilidade dos escoamentos de fluidos viscoelásticos a perturbações não estacionárias. Os resultados numéricos obtidos para escoamentos de fluidos viscoelásticos são comparados com os resultados de escoamentos de fluidos Newtonianos, que já estão bem documentados na comunidade científica. Além disso, uma equação de Orr-Sommerfeld modificada é deduzida para um escoamento viscoelástico utilizando a Teoria de Estabilidade Linear. / Several flows of practical interest are of viscoelastic fluids and it is often desirable to know if these flows are in a laminar or turbulent state. Although the hydrodynamics of viscoelastic fluids are strongly affected by the balance between inertia and elastic forces in the flow, the effect of elasticity on the stability of inertial flows has not been completely established. In this work is studied what happens between these two states, the laminar-turbulent transition. Specifically, it will be investigated the convection of Tollmien-Schlichting waves to incompressible Poiseuille flow of viscoelastic fluid, using the constitutive equation Oldroyd-B. For this, the analysis is carried out by means of Direct Numerical Simulation to verify the stability of the non-stationary disturbances viscoelastic fluids flows. The numerical results obtained for viscoelastic fluids flows are compared with the results of Newtonian fluids flows, which are already well documented in scientific community. In addition, an Orr-Sommerfeld modified equation is deducted for a viscoelastic flow using Linear Stability Theory.

Em experimentação agrícola, é freqüente a necessidade de análise conjunta de grupos de experimentos. Em muitos casos, o pesquisador deseja generalizar resultados para condições gerais de regiões e/ou em avaliar o desempenho de vários genótipos (tratamentos) em diversos ambientes (locais e/ou ano). Quando um conjunto de experimentos é planejado para vários locais é necessário considerar o delineamento individual em cada local e a combinação total dos genótipos com os locais (interação genótipo × ambiente). Logo, os dados observados podem ser organizados em uma tabela de dupla entrada. Existem várias metodologias de análise e interpretação para a interação genótipo × ambiente proveniente de um grupo de cultivares testados em vários ambientes. Entre essas metodologias destaca-se os modelos AMMI (“additive main effects and multiplicative interaction model”), como o próprio nome diz é um método uni-multivariado, que engloba uma análise de variância para os efeitos principais, que são os efeitos dos genótipos (G) e os ambientes (E) e para efeitos multiplicativos (interação genótipo × ambiente), utiliza-se a decomposição em valor singular (DVS). Essa técnica multivariada baseia-se no uso dos autovalores e autovetores provenientes da matriz de interação genótipo × ambiente. Araújo e Dias (2005) verificaram o problema de superestimação e subestimação de autovalores estimados da maneira convencional. Para superar esses problemas de estimação de autovalores, Muirhead (1987) apresenta três métodos para corrigir autovalores estimados a partir das matrizes de covariâncias amostral e alerta que nem sempre essas correções mantêm a ordem decrescente de valores, assim é sugerido o uso de regressão isotônica para ordenar esses dados, mas propriamente um algoritmo apresentado por Lin e Pearlman (1985). Os resultados indicaram que: A regressão isotônica juntamente com o algoritmo foi necessária e se mostrou muito importante em todos conjuntos de dados; Houve uma redução no número de componentes significativos para serem retidos nos modelos, fazendo com que os modelos AMMI selecionados sejam mais parcimoniosos quando se utiliza qualquer um dos métodos de correção; O método 2 apresentou as menores taxa de correção da soma de quadrados da interação genótipo × ambiente e o método 3 apresentou a maiores taxa de correção; Em relação a medida RMSPDPRESS, os menores valores foram obtidos quando se utilizou o método de correção 2. Já o método de correção 3 apresentou os maiores valores para RMSPDPRESS; O método 2 também se mostrou melhor quando o interesse era verificar o ganho em número de repetições, sendo que este benefício esteve sempre próximo de 3 repetições. Já o método 3 é o que apresenta um menor ganho em número de repetições, em torno de duas repetições. / In agricultural research is common to analyse groups of experiments. In many cases, the researcher intends to generalize results to general conditions of areas and/or evaluate the responses of several genotypes (treatments) in several environments (places and/or years). When a group of experiments is planned for several places it is necessary to consider the of design in each place and the combinations of the genotypes with the places (the interaction of genotype × environment). The observed data can be organized in an array. There are several methods of analysis and interpretation for the genotype × environment interaction from a group of genotype tested in several environments. These methods include AMMI models (“additive main effect and multiplicative interaction models”). As the name says it is a uni-multivariate method, that includes an analysis of variance for the main effects (the effects of the genotypes (G) and environments (E)) and assumes multiplicative effects for the genotype × environment interaction, using a singular value decomposition (DVS). This method estimates the eigenvalues and eigenvectors deriving from the matrix of genotype × environment interaction. Araújo and Dias (2005) found an overestimation and underestimation problem with the eigenvalues in the conventional way. To correct these problems Muirhead (1987) presents three methods to correct the eigenvalues from covariance the matrix and noted that these do not always maintain the order of values. The author suggested the use of isotonic regression to correct the eigenvalues, using an algorithm presented by Lin and Pearlman (1985). The results indicated that: The isotonic regression with the algorithm is necessary and it showed very important in all groups of data; There was a reduction in the number of significant components to be kept in the models and the order that the AMMI model selected is more parsimonious when any of the correction methods is used; The method 2 has the smallest rate of correction to the sum of squares of the genotype × environment interaction and method 3 has the largest correction rate; The measure RMSPDPRESS was smallest when method of correction 2 was used. The method of correction 3 has the largest values for RMSPDPRESS; Method 2 was also better when the interest was to verify the gain in number of replicates, and this benefit was always close to 3 replicates. The method 3 gives the smaller gain in the number of replicates, of around two replicates.

Boundary-layer transition from a laminar to a turbulent regime is a critical driver in the design of high-speed vehicles. The aerothermodynamic loads associated with transitional or fully turbulent hypersonic boundary layers are several times higher than those associated with laminar flow. The presence of isolated roughness elements on the surface of a body can accelerate the growth of incoming disturbances and introduce additional instability mechanisms in the flow field, eventually leading to a premature occurrence of transition. This dissertation studies the instabilities induced by three-dimensional discrete roughness elements located inside a high-speed boundary layer developing on a flat plate. Two-dimensional local linear stability theory (2D-LST) is employed to identify the instabilities evolving in the three-dimensional flow field that characterizes the wake induced by the roughness elements and to investigate their evolution downstream. A formulation of the disturbance energy evolution equation available for base flows depending on a single spatial direction is generalized for the first time to base flows featuring two inhomogeneous directions and perturbations depending on three spatial directions. This generalization allows to obtain a decomposition of the temporal growth rate of 2D-LST instabilities into the different contributions that lead to the production and dissipation of the total disturbance energy. This novel extension of the formulation provides an additional layer of information for understanding the energy exchange mechanisms between a three-dimensional base flow and the perturbations resulting from 2D-LST. Stability computations for a calorically perfect gas illustrate that the wake induced by the roughness elements supports the growth of different sinuous and varicose instabilities which coexist together with the Mack-mode perturbations that evolve in the flat-plate boundary layer, and which become modulated by the roughness-element wake. A single pair of sinuous and varicose disturbances is found to dominate the wake instability in the vicinity of the obstacles. The application of the newly developed decomposition of the temporal growth rate reveals that the roughness-induced wake modes extract most of their potential energy from the transport of entropy fluctuations across the base-flow temperature gradients and most of their kinetic energy from the work of the disturbance Reynolds stresses against the base-flow velocity gradients. Further downstream, the growth rate of the wake instabilities is found to be influenced by the presence of Mack-mode disturbances developing on the flat plate. Strong evidence is observed of a continuous synchronization mechanism between the wake instabilities and the Mack-mode perturbations. This phenomenon leads to an enhancement of the amplification rate of the wake modes far downstream of the roughness element, ultimately increasing the associated integrated amplification factors for some of the investigated conditions. The effects of vibrational molecular excitation and chemical non-equilibrium on the instabilities induced by a roughness element are studied for the case of a high-temperature boundary layer developing on a sharp wedge configuration. For this purpose, a 2D-LST solver for chemical non-equilibrium flows is developed for the first time, featuring a fully consistent implementation of the thermal and transport models employed for the base flow and the perturbation fields. This is achieved thanks to the automatic derivation and implementation tool (ADIT) available within the von Karman Institute extensible stability and transition analysis (VESTA) tool-kit, which enables an automatic derivation and implementation of the 2D-LST governing equations for different thermodynamic flow assumptions and models. The stability computations for this configuration show that sinuous and varicose disturbances also dominate the wake instability in the presence of vibrational molecular energy mode excitation and chemical reactions. The resulting base-flow cooling associated with the modeling of such high-temperature phenomena is found to have opposite stabilizing and destabilizing effects on the streamwise evolution of the sinuous and varicose instabilities. The modeling of vibrational excitation and chemical non-equilibrium acting exclusively on the perturbations is found to have a stabilizing influence in all cases. / Doctorat en Sciences de l'ingénieur et technologie / info:eu-repo/semantics/nonPublished

<p>Accurate and reliable prediction of laminar-turbulent boundary layer transition at hypersonic velocities is important for the development of a variety of practical high-speed flight systems currently under development. Boundary layer transition can cause up to an order of magnitude increase in skin friction and heat flux on a flight vehicle, meaning that understanding boundary layer behavior is critical to the design of weight-efficient thermal protection systems. Despite the importance of the topic, significant gaps remain in the community's current understanding of boundary layer transition and control. </p> <p>One of the biggest areas of concern in the field of high-speed boundary layer transition is the effect of facility noise on wind tunnel measurements. Conventional hypersonic wind tunnels are contaminated by freestream fluctuations which can be as much as two orders of magnitude higher than free-flight atmospheric conditions. These disturbances are typically produced by turbulent boundary layers on the tunnel walls; they are acoustic in nature and consist of pressure waves which radiate into the test section. This facility noise plays a leading role in high-speed transition phenomena in conventional hypersonic tunnels.</p> <p><br></p> <p>The current work studies the effects of facility noise on hypersonic transition using both linear stability theory and direct numerical simulation. A model for the freestream disturbance environment of the von Karman Facility's Tunnel B based on experimental measurements of the disturbance spectra present in the tunnel is created and used to study a past experiment performed in the same wind tunnel using a sharp cone and hollow cylinder. The results show that while linear stability theory accurately captures the behavior of second-mode instability growth, it fails to predict the growth of low-frequency instabilities recorded in the experiments. The stability theory analysis also suggests that very fine scale variation in nose tip geometry can play an outsize role in the development of boundary layer instabilities significantly farther downstream.</p> <p><br></p> <p>The direct numerical simulation demonstrates that, using an artificial body forcing term to implement the constructed tunnel noise model, the experimental effects of facility noise can be adequately captured with a sufficiently dense computational grid. For the conical geometry used in the experiments, calculations of surface heat flux indicate good experimental agreement with in prediction of transition location, and total temperature spectra extracted from the flow compare favorably with the experimental data as well. Visualizations of the flowfield confirm the onset of turbulence as a result of the freestream forcing. The computations also suggest that nonlinear interactions may be present in the turbulent breakdown region, leading to the production of streamwise streaks along the cone's surface. Transition on the hollow cylinder was not achieved due to suspected resolution issues, so detailed physical comparison of the two cases was not possible.</p> <p><br></p> <p>Overall, the results of this work suggest that direct numerical simulation is a capable tool for studying the effects of facility noise on hypersonic transition for simple geometries, albeit one which can be difficult to practically realize considering the required computational cost. Computational results indicate that two phenomena may play a role in the development of boundary layer instabilities for a sharp cone --- the fine-scale shape of the tip, which may change the behavior of the entropy layer near the nose; and the interactions between low- and high-frequency waveforms, which seems to cause nonlinear breakdown in line with current experimental understanding.</p>

A hybrid system is a dynamical system that exhibits both continuous and discrete dynamic behavior. Hybrid systems arise in a wide variety of important applications in diverse areas, ranging from biology to computer science to air traffic dynamics. The interaction of continuous- and discrete-time dynamics in a hybrid system often leads to very rich dynamical behavior and phenomena that are not encountered in purely continuous- or discrete-time systems. Investigating the dynamical behavior of hybrid systems is of great theoretical and practical importance. The objectives of this thesis are to develop the qualitative theory of nonlinear hybrid systems with impulses, time-delay, switching modes, and stochastic disturbances, to develop algorithms and perform analysis for hybrid systems with an emphasis on stability and control, and to apply the theory and methods to real-world application problems. Switched nonlinear systems are formulated as a family of nonlinear differential equations, called subsystems, together with a switching signal that selects the continuous dynamics among the subsystems. Uniform stability is studied emphasizing the situation where both stable and unstable subsystems are present. Uniformity of stability refers to both the initial time and a family of switching signals. Stabilization of nonlinear systems via state-dependent switching signal is investigated. Based on assumptions on a convex linear combination of the nonlinear vector fields, a generalized minimal rule is proposed to generate stabilizing switching signals that are well-defined and do not exhibit chattering or Zeno behavior. Impulsive switched systems are hybrid systems exhibiting both impulse and switching effects, and are mathematically formulated as a switched nonlinear system coupled with a sequence of nonlinear difference equations that act on the switched system at discrete times. Impulsive switching signals integrate both impulsive and switching laws that specify when and how impulses and switching occur. Invariance principles can be used to investigate asymptotic stability in the absence of a strict Lyapunov function. An invariance principle is established for impulsive switched systems under weak dwell-time signals. Applications of this invariance principle provide several asymptotic stability criteria. Input-to-state stability notions are formulated in terms of two different measures, which not only unify various stability notions under the stability theory in two measures, but also bridge this theory with the existent input/output theories for nonlinear systems. Input-to-state stability results are obtained for impulsive switched systems under generalized dwell-time signals. Hybrid time-delay systems are hybrid systems with dependence on the past states of the systems. Switched delay systems and impulsive switched systems are special classes of hybrid time-delay systems. Both invariance property and input-to-state stability are extended to cover hybrid time-delay systems. Stochastic hybrid systems are hybrid systems subject to random disturbances, and are formulated using stochastic differential equations. Focused on stochastic hybrid systems with time-delay, a fundamental theory regarding existence and uniqueness of solutions is established. Stabilization schemes for stochastic delay systems using state-dependent switching and stabilizing impulses are proposed, both emphasizing the situation where all the subsystems are unstable. Concerning general stochastic hybrid systems with time-delay, the Razumikhin technique and multiple Lyapunov functions are combined to obtain several Razumikhin-type theorems on both moment and almost sure stability of stochastic hybrid systems with time-delay. Consensus problems in networked multi-agent systems and global convergence of artificial neural networks are related to qualitative studies of hybrid systems in the sense that dynamic switching, impulsive effects, communication time-delays, and random disturbances are ubiquitous in networked systems. Consensus protocols are proposed for reaching consensus among networked agents despite switching network topologies, communication time-delays, and measurement noises. Focused on neural networks with discontinuous neuron activation functions and mixed time-delays, sufficient conditions for existence and uniqueness of equilibrium and global convergence and stability are derived using both linear matrix inequalities and M-matrix type conditions. Numerical examples and simulations are presented throughout this thesis to illustrate the theoretical results.

A hybrid system is a dynamical system that exhibits both continuous and discrete dynamic behavior. Hybrid systems arise in a wide variety of important applications in diverse areas, ranging from biology to computer science to air traffic dynamics. The interaction of continuous- and discrete-time dynamics in a hybrid system often leads to very rich dynamical behavior and phenomena that are not encountered in purely continuous- or discrete-time systems. Investigating the dynamical behavior of hybrid systems is of great theoretical and practical importance. The objectives of this thesis are to develop the qualitative theory of nonlinear hybrid systems with impulses, time-delay, switching modes, and stochastic disturbances, to develop algorithms and perform analysis for hybrid systems with an emphasis on stability and control, and to apply the theory and methods to real-world application problems. Switched nonlinear systems are formulated as a family of nonlinear differential equations, called subsystems, together with a switching signal that selects the continuous dynamics among the subsystems. Uniform stability is studied emphasizing the situation where both stable and unstable subsystems are present. Uniformity of stability refers to both the initial time and a family of switching signals. Stabilization of nonlinear systems via state-dependent switching signal is investigated. Based on assumptions on a convex linear combination of the nonlinear vector fields, a generalized minimal rule is proposed to generate stabilizing switching signals that are well-defined and do not exhibit chattering or Zeno behavior. Impulsive switched systems are hybrid systems exhibiting both impulse and switching effects, and are mathematically formulated as a switched nonlinear system coupled with a sequence of nonlinear difference equations that act on the switched system at discrete times. Impulsive switching signals integrate both impulsive and switching laws that specify when and how impulses and switching occur. Invariance principles can be used to investigate asymptotic stability in the absence of a strict Lyapunov function. An invariance principle is established for impulsive switched systems under weak dwell-time signals. Applications of this invariance principle provide several asymptotic stability criteria. Input-to-state stability notions are formulated in terms of two different measures, which not only unify various stability notions under the stability theory in two measures, but also bridge this theory with the existent input/output theories for nonlinear systems. Input-to-state stability results are obtained for impulsive switched systems under generalized dwell-time signals. Hybrid time-delay systems are hybrid systems with dependence on the past states of the systems. Switched delay systems and impulsive switched systems are special classes of hybrid time-delay systems. Both invariance property and input-to-state stability are extended to cover hybrid time-delay systems. Stochastic hybrid systems are hybrid systems subject to random disturbances, and are formulated using stochastic differential equations. Focused on stochastic hybrid systems with time-delay, a fundamental theory regarding existence and uniqueness of solutions is established. Stabilization schemes for stochastic delay systems using state-dependent switching and stabilizing impulses are proposed, both emphasizing the situation where all the subsystems are unstable. Concerning general stochastic hybrid systems with time-delay, the Razumikhin technique and multiple Lyapunov functions are combined to obtain several Razumikhin-type theorems on both moment and almost sure stability of stochastic hybrid systems with time-delay. Consensus problems in networked multi-agent systems and global convergence of artificial neural networks are related to qualitative studies of hybrid systems in the sense that dynamic switching, impulsive effects, communication time-delays, and random disturbances are ubiquitous in networked systems. Consensus protocols are proposed for reaching consensus among networked agents despite switching network topologies, communication time-delays, and measurement noises. Focused on neural networks with discontinuous neuron activation functions and mixed time-delays, sufficient conditions for existence and uniqueness of equilibrium and global convergence and stability are derived using both linear matrix inequalities and M-matrix type conditions. Numerical examples and simulations are presented throughout this thesis to illustrate the theoretical results.