From the nano scale to the macro scale, large models are used to simulate and predict the responses of dynamic systems. The construction and evaluation of such models, often in the form of finite element models, require tremendous computational resources and time. Due to this large computational endeavor, it is paramount to learn as much as possible from the models and their solutions. In this work, analyses and methods for efficiently deriving significant knowledge of damped systems from models and their solutions are presented. Of primary interest to this work is the analysis of damped structures. Damping, the means by which energy is dissipated, often adds an additional layer of complexity to finite element models and any subsequent analyses. This added complexity is due to the relative complexity of many damping models and their accompanying computational burden. Furthermore, on the micro and nano scale, a variety of damping mechanisms, each with their own unique set of physics, may be present. The research presented in this work is organized in two parts. The first part presents methods for deriving knowledge from models and their solutions. Here, the developed methods perform approximate yet highly efficient analysis on the matrices and solution vectors of finite element models. In this work, methods utilizing the Neumann series approximation are presented. These methods efficiently predict how the response of a structure depends on its damping or any other input model parameter. Additionally, a method for analyzing the spatial dependence of damping with the use of loss factor images is presented. Research presented in the second part derives knowledge solely from solutions of models. In this part, it is assumed that the matrices of the models are not available, and therefore analysis is restricted to the solution itself. Here, research is focused on the analyses of structures on the micro and nano scale. Specifically, micro and nano beams surrounded by a viscous compressible fluid are analyzed. The dynamic responses of the structure and the surrounding fluid are analyzed to determine the prominent damping mechanisms. Here, results from 2--Dimensional analytical models and 3--Dimensional finite element models are complemented by experimental measurements to analyze damping due to viscous dissipation and acoustic radiation.

All organs, including the brain, need a pathway to remove neurotoxic extracellular proteins. In the brain this is called the glymphatic system. The glymphatic system works by exchanging proteins from interstitial fluids to cerebrospinal fluids. The extracellular proteins are then removed through the cerebrospinal fluid drains. The glymphatic system is believed to be driven by arterial pulsatility, cerebrospinal fluid production and respiration. Cerebrospinal fluids enters the brain alongside arteries. In this project, we investigate if a simulated pulsatile flow in a common carotid artery can drive cerebrospinal fluid flow running along the artery, using computational simulations of a linearly elastic and fluid-structure multiphysical model in COMSOL. Our simulations show that a heartbeat pulse increases the arterial radius of the common carotid artery by 6 %. Experimental data, assessed using 4D magnetic resonance imaging of a living human, show an increase of 13 %. Moreover, our results indicate that arterial displacement itself is not able to drive cerebrospinal fluid flow. Instead, it seems to create a back and forth flow that possibly could help with the protein exchange between the cerebrospinal and interstitial fluids. Overall, the results indicate that the COMSOL Multiphysics linearly elastic model is not ideal for approximations of soft non-linearly elastic solids, such as soft polydimethylsiloxane and artery walls work for stiffer materials. The long term aim is to simulate a part of the glymphatic system and the present work is a starting point to reach this goal. As the simulations in this work are simplified there are more things to test in the future. For example, using the same geometries a non-linear elastic model could be tested. The pulsatile waveform or the geometry could be made more complex. Furthermore the model could be scaled down to represent a penetrating artery in the brain instead of the common carotid artery.

Dans l’industrie nucléaire, la simulation de transitoires accidentels à l’échelle d’un réacteur devient une composante d’importance croissante de la démonstration de sûreté à destination des agences de surveillance nationales. Elle permet ainsi de limiter le recours à des expériences complexes et coûteuses tout en facilitant l’évaluation des stratégies de mitigation. Cependant, les modèles mis en jeu sont inévitablement volumineux et construits avec une finesse de modélisation rendant difficile la prise en compte de détails géométriques locaux pourtant susceptibles d’influencer significativement la solution globale. Dans ce travail de thèse, on propose ainsi des approches multi-modèles pour l’intégration de tels détails dans un modèle global sans modification du maillage initial (on parle aussi de zoom numérique). Des techniques sont proposées aussi bien pour les structures que pour les fluides, avec un souci de démontrer la précision et la stabilité de la solution multi-modèles couplée comparée à une solution de référence à une seule échelle. Ce travail intègre deux spécificités propres, à savoir son adéquation avec les contraintes de la dynamique rapide avec intégration temporelle explicite et l’objectif de traiter simultanément la superposition de modèles et l’interaction fluide-structure. / In nuclear industry, simulating accidental transient sequences at full reactor scale is becoming an increasingly important feature of the safety demonstration towards national agencies. It thus allows limiting the number of complex and costly experiments, while simplifying and accelerating the evaluation of mitigation strategies. However, the implemented numerical models are inevitably heavy to build and maintain, with a global modelling scale making it difficult to account for local geometric details yet able to significantly influence the physical solution. To provide an answer to these problems, this PhD work is dedicated to multi-model approaches designed to integrate such details into bigger models with no modification at the global level (techniques often designated as numerical zoom techniques). Some methods are proposed for both structures and fluids, with special care given to the accuracy and stability of the coupled multi-scale solution compared to a single-scale reference solution. This work handles two very specific topics, namely its compatibility with numerical features imposed by fast transient dynamics with explicit time integration, and the general objective of simultaneously dealing with superimposed models and fluid-structure interaction.

Wind storms pose dangerous threats to human lives and are an enormous drain on the economy. Their damage to buildings usually starts with the failure of structural components that are subjected to excessive wind loads. In this dissertation, we investigate the characteristics of extreme loads on low-rise structures through analysis of full-scale and numerical data. We also use numerical simulations to evaluate different approaches to control the separated flow over a surface-mounted prism with the objective of reducing extreme pressure coefficients or loads on its surface. In the first part, we use a probabilistic approach to characterize peak loads as measured on a subject house during Hurricane Ivan on 2004. Time series of pressure coefficients collected on the roof of that house are analyzed. Rather than using peak values, which could vary due to the stochastic nature of the data, a probabilistic analysis is used to determine the probability of non-exceedence of specific values of pressure coefficients and associated wind loads. The results show that the time series of the pressure coefficients follow a three-parameter Gamma distribution, while the peak pressure follows a two-parameter Gumbel distribution. The results of the analysis are contrasted with the design values. In the second part, we perform numerical simulations of the flow over a surface-mounted prism as a simplified example for the flow over a low-rise structure. A Direct Numerical Simulation (DNS) code is developed to solve the unsteady two-dimensional incompressible Navier-Stokes equations of the flow past the prism. The pressure coefficients are then computed on the prism surface in order to assess the wind loads. The code is written on a parallel platform using the Message Passing Interface (MPI) library. We use the simulations to study the effects of inflow disturbances on the extreme loads on structures. The sensitivities of peak loads on a surface mounted prism to variations in incident gust parameters are determined. Latin Hypercube Sampling (LHS) is applied to obtain different combinations of inflow parameters. A non-intrusive polynomial chaos expansion is then applied to determine the sensitivities. The results show that the gust enhances the destabilization of the separation shear layer, forces it to break down and moves it closer to the roof of the prism. As for the sensitivities, the results show that the extreme loads are most sensitive to the transverse amplitude of the disturbance. Because the separated flow over sharp edges is responsible for the extreme pressure peaks, we investigate the use of active and passive control strategies to reduce wind loads. The studied active flow control strategies include blowing, suction, and synthetic jets. We implement them by using different flux injections, different slot locations and different angles. Investigation of the possible peak pressure reduction for two Reynolds numbers is performed. For Re = 1000, a reduction by nearly 50% of the peak pressure is obtained. For Re = 10, 000, the highest achieved reduction is nearly 25%. For passive control, we mount a flexible membrane on the top of the prism. In a two-dimensional framework, the membrane equation is modeled by a forced string equation. This mechanical equation is coupled with the DNS solver and integrated in time using a fourth order Hamming predictor corrector scheme. The results show that this strategy is as efficient as the active control approach, in terms of reducing extreme loads, for Re = 10, 000. / Ph. D.

Indiana University-Purdue University Indianapolis (IUPUI) / Glaucoma is a multifactorial ocular disease progressively leading to irreversible blindness. There is clear evidence of correlations between alterations in ocular hemodynamics and glaucoma; however, the mechanisms giving rise to these correlations are still elusive. The objective of this thesis is to develop mathematical models and methods to help elucidate these mechanisms. First, we develop a mathematical model that describes the deformation of ocular structures and ocular blood flow using a reduced-order fluid-structure interaction model. This model is used to investigate the relevance of mechanical and vascular factors in glaucoma. As a first step in expanding this model to higher dimensions, we propose a novel energy-based technique for coupling partial and ordinary differential equations in blood flow, using operator splitting. Next, we combine clinical data and model predictions to propose possible explanations for the increase in venous oxygen saturation in advanced glaucoma patients. We develop a computer-aided manipulation process of color Doppler images to extract novel waveform parameters to distinguish between healthy and glaucomatous individuals. The results obtained in this work suggest that: 1) the increase in resistance of the retinal microcirculation contributes to the influence of intraocular pressure on retinal hemodynamics; 2) the influence of cerebrospinal fluid pressure on retinal hemodynamics is mediated by associated changes in blood pressure; 3) the increase in venous oxygen saturation levels observed among advanced glaucoma patients depends on the value of the patients’ intraocular pressure; 4) the normalized distance between the ascending and descending limb of the ophthalmic artery velocity profile is significantly higher in glaucoma patients than in healthy individuals.

京都大学 / 0048 / 新制・課程博士 / 博士(工学) / 甲第21350号 / 工博第4509号 / 新制||工||1702(附属図書館) / 京都大学大学院工学研究科社会基盤工学専攻 / (主査)教授 後藤 仁志, 教授 KIM Chul-Woo, 准教授 KHAYYER,Abbas / 学位規則第4条第1項該当 / Doctor of Philosophy (Engineering) / Kyoto University / DGAM

No description available.

Sweden is generally considered to be a low seismicity area, but for structures within nuclear power facilities, the safety level demand with respect to seismic events are high and thus, these structures are required to be earthquake-resistant. The seismic hazard is here primarily considered to be associated with near-field earthquakes. The nuclear power plants are further founded on hard rock and the expected ground motions are dominated by high frequencies. The design earthquake considered for the nuclear facilities has an annual probability of 10-5 events, that is, the probability of occurrence is once per 100 000 years. The focus of the study is the seismic response of large concrete structures for the nuclear power industry, with regard not only to the structure itself but also to non-structural components attached to the primary structure, and with emphasis on Swedish conditions. The aim of this licentiate thesis is to summarize and demonstrate some important aspects when the seismic load is dominated by high frequencies. Additionally, an overview of laws, regulations, codes, standards, and guidelines important for seismic analysis and design of nuclear power structures is provided. The thesis includes two case studies investigating the effect of seismic high-frequency content loads. The first study investigates the influence of gaps in the piping supports on the response of a steel piping system subjected to a seismic load dominated by high amplitudes at high frequencies. The gaps are found in the joints of the strut supports or are gaps between the rigid box supports and the pipe. The piping system is assessed to be susceptible to high-frequency loads and is located within the reactor containment building of a nuclear power plant. The stress response of the pipe and the acceleration response of the valves are evaluated. The second study investigates the effect of fluid-structure interaction (FSI) on the response of an elevated rectangular water-containing concrete pool subjected to a seismic load with dominating low and high frequencies, respectively. The pool is located within the reactor containment building of a boiling water reactor at a nuclear power plant. The hydrodynamic pressure distribution is evaluated together with the stress distribution in the walls of the tank. From the two case studies, it is evident that the response due to a seismic load dominated by high frequencies and low frequencies, respectively, is different. Although the seismic high-frequency load may be considered non-damaging for the structure, the effect may not be negligible for non-structural components attached to the primary structure. Including geometrical non-linear effects such as gaps may however reduce the response. It was shown that the stress response for most of the pipe elements in the first case study was reduced due to the gaps. It may also be that the inclusion of fluid-structure interaction effects changes the dynamic properties of a structural system so that it responds significantly in the high frequency range, thus making it more vulnerable to seismic loads dominated by high frequencies. In the second case study, it was shown that even for a seismic load with small amplitudes and short duration, but with dominating high-frequency content, as the Swedish 10-5 design earthquake, the increase of the dynamic response as fluid-structure interaction is accounted for is significant. / <p>QC 20150519</p>

No description available.

No description available.