The aim of this study is to investigate the fluid-structural response to pulsatile Newtonian and non-Newtonian blood flow through an axisymmetric stenosed vessel using FLOTRAN and ANSYS. This is to provide a basic understanding of atherosclerosis. The flow was set to be laminar and follows a sinusoidal waveform. The solid model was set to have isotropic elastic properties. The Fluid-Structural Interaction (FSI) coupling was two-way and iterative. Rigid and Newtonian cases were investigated to provide an understanding on the effects of incorporating FSI into the model. The wall expansion was found to decrease the axial velocity and increase the recirculation effects of the flow. To validate the models and methods used, the results were compared with the study by Lee and Xu [2002] and Ohja et al [1989]. Close comparisons were achieved, suggesting the models used were valid. Two non-Newtonian models were investigated with FSI: Carreau and Power Law models. The Carreau model fluid behaviour was very close to the Newtonian model. The Power Law model produced significant difference in viscosity, velocity and wall shear stress distributions. Pressure distribution for all models was similar. In order to quantify the changes, Importance Factor (IG) was introduced to determine the overall non-Newtonian effects at two regions: the entire flow model and about the vessel wall. The Carreau model showed reasonable values of IG whereas the Power Law model showed excessive values. Transient and geometrical effects were found to affect the Importance Factor. The stress distributions for all models were found to be similar. Highest stress occurred at the shoulders of the stenosis where a stress concentration occurred due to sharp corners of the geometry and large bending moments. The highest stresses were in the axial direction. Notable circumferential stress was found at the ends of the vessel. Carreau model produced slightly higher stresses than the other models. Wall stresses were found to be primarily influenced by internal pressure, rather than wall shear stresses.

The objective of this research is to develop an efficient and accurate methodology to resolve flow non-linearity of fluid-structural interaction. To achieve this purpose, a numerical strategy to apply the detached-eddy simulation (DES) with a fully coupled fluid-structural interaction model is established for the first time. The following novel numerical algorithms are also created: a general sub-domain boundary mapping procedure for parallel computation to reduce wall clock simulation time, an efficient and low diffusion E-CUSP (LDE) scheme used as a Riemann solver to resolve discontinuities with minimal numerical dissipation, and an implicit high order accuracy weighted essentially non-oscillatory (WENO) scheme to capture shock waves. The Detached-Eddy Simulation is based on the model proposed by Spalart in 1997. Near solid walls within wall boundary layers, the Reynolds averaged Navier-Stokes (RANS) equations are solved. Outside of the wall boundary layers, the 3D filtered compressible Navier-Stokes equations are solved based on large eddy simulation(LES). The Spalart-Allmaras one equation turbulence model is solved to provide the Reynolds stresses in the RANS region and the subgrid scale stresses in the LES region. An improved 5th order finite differencing weighted essentially non-oscillatory (WENO) scheme with an optimized epsilon value is employed for the inviscid fluxes. The new LDE scheme used with the WENO scheme is able to capture crisp shock profiles and exact contact surfaces. A set of fully conservative 4th order finite central differencing schemes are used for the viscous terms. The 3D Navier-Stokes equations are discretized based on a conservative finite differencing scheme, which is implemented by shifting the solution points half grid interval in each direction on the computational domain. The solution points are hence located in the center of the grid cells in the computational domain (not physical domain). This makes it possible to use the same code structure as a 2nd order finite volume method. A finite differencing high order WENO scheme is used since a finite differencing WENO scheme is much more efficient than a finite volume WENO scheme. The unfactored line Gauss-Seidel relaxation iteration is employed for time marching. For the time accurate unsteady simulation, the temporal terms are discretized using the 2nd order accuracy backward differencing. A pseudo temporal term is introduced for the unsteady calculation following Jameson's method. Within each physical time step, the solution is iterated until converged based on pseudo time step. A general sub-domain boundary mapping procedure is developed for arbitrary topology multi-block structured grids with grid points matched on sub-domain boundaries. The interface of two adjacent blocks is uniquely defined according to each local mesh index system (MIS) which is specified independently. A pack/unpack procedure based on the definition of the interface is developed to exchange the data in a 1D array to minimize data communication. A secure send/receive procedure is employed to remove the possibility of blocked communication and achieve optimum parallel computation efficiency. Two terms, "Order" and "Orientation", are introduced as the logics defining the relationship of adjacent blocks. The domain partitioning treatment of the implicit matrices is to simply discard the corner matrices so that the implicit Gauss-Seidel iteration can be implemented within each subdomain. This general sub-domain boundary mapping procedure is demonstrated to have high scalability. Extensive numerical experiments are conducted to test the performance of the numerical algorithms. The LDE scheme is compared with the Roe scheme for their behavior with RANS simulation. Both the LDE and the Roe scheme can use high CFL numbers and achieve high convergence rates for the algebraic Baldwin-Lomax turbulence model. For the Spalart-Allmaras one equation turbulence model, the extra equation changes the Jacobian of the Roe scheme and weakens the diagonal dominance. It reduces the maximum CFL number permitted by the Roe scheme and hence decreases the convergence rate. The LDE scheme is only slightly affected by the extra equation and maintains high CFL number and convergence rate. The high stability and convergence rate using the Spalart-Allmaras one equation turbulence model is important since the DES uses the same transport equation for the turbulence stresses closure. The RANS simulation with the Spalart-Allmaras one equation turbulence model is the foundation for DES and is hence validated with other transonic flows including a 2D subsonic flat plate turbulent boundary layer, 2D transonic inlet-diffuser, 2D RAE2822 airfoil, 3D ONERA M6 wing, and a 3D transonic duct with shock boundary layer interaction. The predicted results agree very well with the experiments. The RANS code is then further used to study the slot size effect of a co-flow jet (CFJ) airfoil. The DES solver with fully coupled fluid-structural interaction methodology is validated with vortex induced vibration of a cylinder and a transonic forced pitching airfoil. For the cylinder, the laminar Navier-Stokes equations are solved due to the low Reynolds number. The 3D effects are observed in both stationary and oscillating cylinder simulation because of the flow separations behind the cylinder. For the transonic forced pitching airfoil DES computation, there is no flow separation in the flow field. The DES results agree well with the RANS results. These two cases indicate that the DES is more effective on predicting flow separation. The DES code is used to simulate the limited cycle oscillation of NLR7301 airfoil. For the cases computed in this research, the predicted LCO frequency, amplitudes, averaged lift and moment, all agree excellently with the experiment. The solutions appear to have bifurcation and are dependent on the initial perturbation. The developed methodology is able to capture the LCO with very small amplitudes measured in the experiment. This is attributed to the high order low diffusion schemes, fully coupled FSI model, and the turbulence model used. This research appears to be the first time that a numerical simulation of LCO matches the experiment. The DES code is also used to simulate the CFJ airfoil jet mixing at high angle of attack. In conclusion, the numerical strategy of the high order DES with fully coupled FSI model and parallel computing developed in this research is demonstrated to have high accuracy, robustness, and efficiency. Future work to further maturate the methodology is suggested.

Brush seals are contact seals which are efficient and reliable and can beused in any rotatory machinery. A theoretical model of a brush seal wassuggested and used for a simulation study. Comparison between the resultsof the experiment and simulation is used to verify the accuracy of model.Following the basic simulation steps, more simulation will be done to geta further analysis. The further analysis will be studied in bristles’ bucklingcharacteristics, deformation and stress. The buckling, stress anddeformation is related to the brush seal’s performance. This work givesthe different geometry of bristle’s effect to its buckling characteristics andthe study of stress and deformation caused by fluid flow across bristlesduring operation.

Dans le processus de conception des véhicules, la vibration et le bruit sont des sujets d’étude très importants. En effet, les vibrations sont susceptibles d’affecter le comportement dynamique des structures et le bruit dégrade le confort acoustique des passagers. L’objectif principal de la thèse est de proposer un ensemble de méthodes pour l’optimisation du comportement dynamique des systèmes complexes afin de réduire les vibrations des structures et le bruit dans l’habitacle. À cet effet, on s’intéresse à des stratégies de contrôle des interfaces, comme le collage de couches viscoélastiques sur les zones les plus déformées, ou l’introduction de dispositifs frottants calibrés pour ajouter de l’amortissement à certaines fréquences de résonance. Les structures assemblées résultantes sont étudiées numériquement par une méthode de synthèse modale généralisée. La méthode de synthèse modale proposée contient plusieurs niveaux de condensation. Le premier concerne les degrés de libertés (DDL) internes de chaque sous-structure. La deuxième condensation s’effectue sur les modes de branches, de sorte à réduire le nombre de DDL aux interfaces entre les sous-structures. Pour les systèmes couplés fluide/structure, une troisième condensation portant sur les DDL du fluide est proposée. Suite à ces condensations, la dimension du système est fortement réduite. Cette méthode permet alors d’obtenir aussi bien la réponse forcée de la structure que les fluctuations du champ de pression dans le fluide. Les chemins de transmission acoustiques et vibratoires peuvent également être déduits des contributions modales intermédiaires. On montre que ces paramètres modaux peuvent être utilisés comme fonctions objectif pour une démarche d’optimisation des interfaces. Le front de Pareto des conceptions optimales est obtenu avec un algorithme génétique multi-objectif élitiste, appliqué à une approximation par krigeage de la fonction objectif. Cette approche modale est étendue à l’étude de systèmes non-linéaires. L’hypothèse fondamentale est que les modes non-linéaires sont faiblement couplés. Les paramètres modaux non-linéaires (fréquences propres, amortissements...), dépendent des amplitudes modales. L’idée est alors de calculer des modes normaux non-linéaires en fonction de leur amplitude et de superposer leurs réponses pour obtenir celle de la structure. La méthode est appliquée à des systèmes incorporant des non-linéarités de type Duffing et de frottement sec. Le cas particulier du frottement sec est considéré à travers un modèle de Masing généralisé. Deux approches modales sont développées : l’une basée sur les modes complexes, et l’autre basée sur les modes réels. L’utilisation de modes complexes ou réels dans la synthèse modale conduit à des termes d’amortissement par frottement différents. On montre que la synthèse modale non-linéaire combinée au modèle de Masing généralisé aboutit à une méthode numérique simple, rapide et efficace pour décrire le comportement non-linéaire de structures soumise à du frottement sec. / Noise and vibration are important topics in the automotive industry for several reasons, including passenger comfort and structural integrity. The main objective of this thesis is to propose a series of appropriate methods to optimize structural system characteristics, so that the vibration and noise can be reduced. To achieve this goal, interface control strategies are employed, including bonding viscoelastic layers onto the most heavily deformed zones and introducing frictional damping devices calibrated on certain resonance frequencies. Such built-up structural systems are numerically investigated via a generalized modal synthesis approach that incorporates several groups of modes. The employed modal synthesis approach consists of several levels of condensation. The first one is on the internal degrees of freedoms (DOFs) of each substructure, and the second condensation is on the branch modes so as to reduce the boundary DOFs among substructures. For coupled fluid-structural systems, a third condensation on the fluid DOFs is suggested. With these condensation techniques, the system dimension can be significantly reduced. The method allows us to obtain the forced response of the structures as well as the pressure variation of the fluids. Additionally, modal parameters characterizing vibration and noise transmission paths can be deduced as mid-stage results. We show that these modal parameters can be used as optimization objective during the interface configuration design. The Pareto front of the optimal design is achieved by employing Kriging approximations followed with an elitist multi-objective genetic algorithm. Another advantage of the modal approach is that a modal overview on the system characteristics is provided by analyzing the natural frequencies, modal damping ratios and the aforementioned modal parameters. The modal synthesis approach is further extended to study nonlinear systems. The basic assumption is that the nonlinear modes are weakly coupled. Nonlinear modal parameters, such as modal frequency and modal damping ratio, contain the essential nonlinear information and depend on modal amplitude. The main idea is to compute nonlinear normal modes according to their modal amplitude and superimpose the response of several nonlinear modes to obtain the overall forced response. The method is applied to systems involving Duffing and dry friction nonlinearities. In the case of dry friction, a generalized Masing model is considered to capture the dry friction nature. Both complex modes and real modes are used in the modal synthesis, leading to different frictional damping terms. We show that the nonlinear modal synthesis combined with the generalized Masing model yields a simple, fast and efficient numerical method to describe nonlinear performance of structures with dry friction.

No description available.