Dissertations / Theses on the topic 'Physical and numerical modeling of furnace'

Thesis: Ph. D., Massachusetts Institute of Technology, Department of Earth, Atmospheric, and Planetary Sciences, 2014.<br>This electronic version was submitted by the student author. The certified thesis is available in the Institute Archives and Special Collections.<br>Cataloged from student-submitted PDF version of thesis.<br>Includes bibliographical references (pages 141-146).<br>The ubiquity of cloud clusters and their role in modulating radiative cooling and the moisture distribution underlines the importance of understanding how and why tropical convection organizes. In this work, the fundamental mechanism underlying the self-aggregation of convection is explored using a cloud resolving model. The objective is to identify and quantify the interactions between the environment and the convection that allow the convection to spontaneously organize into a single cluster. Specifically, the System for Atmospheric Modeling is used to perform 3-d cloud system resolving simulations of radiative-convective equilibrium in a non-rotating framework, with interactive radiation and surface fluxes and fixed sea surface temperature. Self-aggregation only occurs at sea surface temperatures above a certain threshold. As the system evolves to an aggregated state, there are large changes to domain averaged quantities important to climate, such as radiative fluxes and moisture. Notably, self-aggregation begins as a dry patch that expands, eventually forcing all the convection into a single clump. Thus, when examining the initiation of self-aggregation, we focus on processes that can amplify this initial dry patch. Sensitivity tests suggest that wind-dependent surface fluxes and interactive longwave radiative fluxes are important for permitting self-aggregation. A novel method is introduced to quantify the magnitudes of the various feedbacks that control self-aggregation within the framework of the budget for the spatial variance of column - integrated frozen moist static energy. The absorption of shortwave radiation by atmospheric water vapor is found to be a key positive feedback in the evolution of aggregation. In addition, there is a positive wind speed - surface flux feedback whose role is to counteract a negative air-sea enthalpy disequilibrium - surface flux feedback. The longwave radiation - water vapor feedback transitions from positive to negative in the early and intermediate stages of aggregation. The long-wave radiation - cloud feedback is the dominant positive feedback that maintains the aggregated state once it develops. Importantly, the mechanisms that maintain the aggregated state are distinct from those that instigate the evolution of self-aggregation. These results and those of a companion study suggest that the temperature dependence of self-aggregation enters through the longwave feedback term.<br>by Allison A. Wing.<br>Ph. D.

Thesis (Ph. D.)--Joint Program in Physical Oceanography (Massachusetts Institute of Technology, Dept. of Earth, Atmospheric, and Planetary Sciences, and the Woods Hole Oceanographic Institution), 2003.<br>This electronic version was submitted by the student author. The certified thesis is available in the Institute Archives and Special Collections.<br>Includes bibliographical references (p. 190-197).<br>In this thesis, the effects of sea ice, downwelling favorable winds and barotropic background currents on shelf fronts are examined using numerical models. The models are configured with the characteristics of the East Greenland Current, north of the Denmark Strait, in mind. While the models are heavily idealized, basic physical parameters match the observed ocean. The first part of the thesis uses a three dimensional, primitive equation model to examine the behavior of a shelf front under steady, along shelf winds and barotropic currents. The wind stress generates shoreward surface Ekman transport and the barotropic current generates an offshore bottom Ekman transport. In both cases, the Ekman transport causes the creation of mixed layers and a relationship describing the mixed layer thickness is derived relating the cross shelf flux of density to the along shelf flux of density. When there is a barotropic inflow, the cross shelf Ekman mass transport is balanced by a return flow of mass in the interior. When there is a wind stress, in the present model configuration, the influence of the offshore boundary obscures the effect of the Ekman layers. The second part of the thesis focuses on the ice-ocean interaction using a simple, two layer, one dimensional toy model. The interaction of sea-ice, geostrophic currents, and wind are examined. In the presence of a current and the absence of wind, the ice is transported downstream with the current. In the presence of wind and the absence of a current, the net ice-ocean transport is perpendicular to the wind, as is expected for Ekman layer theory. The two layer system acts like a poorly resolved Ekman spiral: the ice has down wind and shoreward transport while the ocean has up wind and shoreward transport.<br>by Sarah L. Russell.<br>Ph.D.

<p> The generation of attosecond-duration (1 as = 10^&ndash;18 s) coherent light through the process of high-order harmonic generation has opened the perspective for probing fundamental processes, such as photoionization, on the natural timescale of electron dynamics in matter. One probing technique is the attosecond streaking method, in which the momentum of the photoelectron is measured as a function of the time delay between the ionizing, attosecond extreme ultraviolet (XUV) pulse, and a weak, femtosecond near-infrared (NIR) pulse which streaks the momentum of the photoelectron, known as a streaking trace. The observed trace contains time information about the photoionization process in the form of a time offset to the vector potential of the streaking field, known as the streaking time delay. Theoretical simulations show that for one-photon ionization this time delay is accumulated by the photoelectron in the continuum when propagating away from the parent ion, whereas for resonant two-photon ionization there exists an additional absorption delay which depends on the properties of the XUV pulse. In this thesis, we use both analytical techniques and numerical simulations to study the contributions of the total time delay observed in streaking, and further explore applications of the streaking time delay to gain insights into the electron dynamics. We first derive an analytical formula for the streaking time delay in one-photon ionization. The predictions based on the model formula, which can be performed within seconds of computation time, are in good agreement with those of computationally extensive numerical simulations. </p><p> We demonstrate that the analytical formula not only allows deeper insight into the nature of the time delay, but also offers the opportunity to effectively analyze other theoretical interpretations and potential effects, such as the effect of a chirp in the ionizing attosecond pulse on the time delay measurement. We then apply time-dependent perturbation theory to derive an analytical formula for the absorption delay in resonant two-photon ionization. We use the analytical formula to demonstrate how the absorption delay can be controlled further by the attosecond pulse duration and central frequency in case of an isolated resonance. Furthermore, we show how multiple resonances within the bandwidth of the ionizing pulses as well as the streaking field influence the absorption delay in model systems as well as simple atoms and molecules. We conclude by exploring the option to apply isolated elliptically polarized attosecond pulses to obtain sub-attosecond temporal information via the observation of photoelectron angular distributions as a function of the ellipticity of the pulse.</p><p>

A comprehensive methodology to physically and numerically model upsetting in cold heading was developed.<br>The physical model was a Drop Weight Test with a guided pocket-die set capable of approximating industrial cold heading conditions. The results show that the test is sensitive to the critical parameters for cold heading. These include surface quality, residual element level, nitrogen content, microstructure, decarburization, and specimen geometry. The test is capable of assessing the fracture behavior of cold heading materials.<br>One goal of the study was to reveal differences in fracture behavior with varying steel sources. Accordingly, the matrix of test materials consisted of grade 1038 steels from three different steel sources.<br>Material preparation and conditioning of test materials approximated industrial procedures for cold heading materials. These procedures included hot rolling, controlled rod cooling, descaling, straightening, lime coating and lubricating, and wire drawing. Spheroidization of test specimens was performed in an industrial batch furnace using an industrial heat treatment cycle.<br>A finite element program (FEM) enabled the simulation of upsetting in cold heading. The inputs required to model the cold heading process include flow stress behavior and friction conditions representative of cold heading. These inputs were obtained using the CANMET Cam Plastometer and the Friction Ring Test.<br>The Cockcroft and Latham fracture constants for an as-rolled and a spheroidize annealed 1038 material were computed by FEM modeling and the critical values were calibrated using the Drop Weight Test. The fracture criterion constant was found to be independent of strain path for upsetting in cold heading and thus is material-related.

Hydropower plants with spillways exposed to high pressures because of large discharges can result in a dam failure which ultimately results in devastating consequences. Therefore properly designed and constructed spillways are important. In this project there is an interest in examining the pressure distribution on a guided wall, measuring the water level inside the spillway and examine the water flow characteristics. The main purpose was to numerically model the flow behaviors at two upper pool levels +29.82 m and +30.92 m in order to achieve a comparison between a numerical model and a physical model. A previous experiment was carried out at Vattenfall R&amp;D in Älvkarleby where a physical model was built. The main purpose of that study was to find an optimal form of a guide wall for free weir flow. This study served as guidelines for what should be investigated in this project. The methodology of the project began with a development of an initial CAD model created by Vattenfall R&amp;D. The CFD model could be constructed into a volume model containing also an upstream reservoir. Next, a mesh that could numerically calculate the flow behavior was created in order to enable execution of the calculations. Furthermore, a big part of the method was to adjust the settings in Fluent so that the numerical model could recreate the flow behaviors of the physical model tests. The results indicate that the numerical model is an accurate replica of the physical model. The deviation comprising the mass flow rate came as close as 0.6 % from the physical model for the lower upper pool level. Highest calculated pressure was positioned at the bottom of the guided wall for both upper pool levels, which the results conducted from the physical model also indicated. Furthermore, results regarding the water level at the guided wall confirmed same statement Vattenfall's examinations pointed out, that the roof placed at the guided wall is needed in order to avoid major splashing. The present project has shown great results regarding the flow characteristics, pressure distribution, and the water level at the guided wall. Therefore, it is considered trustworthy enough to be used as a tool for Vattenfall R&amp;D in future studies comprising the targeted spillway.

<p> Partially water saturated condition in soils may change the cone penetration resistance comparing with that of dry or saturated conditions. This effect was investigated in this study using numerical finite element modeling and experimental centrifuge testing. The results showed suction in unsaturated soil significantly influenced the soil resistance to cone penetration. Two approaches were implemented to numerically consider the partially saturated soil condition; i.e. modifying simple constitutive models using an apparent cohesion strategy and implementing Barcelona Basic Model for unsaturated soils. Both successfully captured the cone resistance profiles inside a calibration chamber and also in free field. In addition, details of developing a miniature cone setup capable of for cone penetration inside geotechnical centrifuge was explained. Further, the use of Linear Softening Cohesive Model (LCFM) to predict the fracture growth in shale rocks during Brazilian Test was examined. The application and importance of considering two different compressive and tensile elastic modulus and soil anisotropy during the fracture modeling of shales were demonstrated.</p><p>

Scientific research in areas of physics, chemistry, and biology traditionally depends purely on experimental and theoretical methods. Recently numerical simulation is emerging as the third way of science discovery beyond the experimental and theoretical approaches. This work describes some general procedures in numerical computation, and presents several applications of numerical modeling in bioheat transfer and biomechanics, jet diffusion flame, and bio-molecular interactions of proteins in blood circulation. A three-dimensional (3D) multilayer model based on the skin physical structure is developed to investigate the transient thermal response of human skin subject to external heating. The temperature distribution of the skin is modeled by a bioheat transfer equation. Different from existing models, the current model includes water evaporation and diffusion, where the rate of water evaporation is determined based on the theory of laminar boundary layer. The time-dependent equation is discretized using the Crank-Nicolson scheme. The large sparse linear system resulted from discretizing the governing partial differential equation is solved by GMRES solver. The jet diffusion flame is simulated by fluid flow and chemical reaction. The second-order backward Euler scheme is applied for the time dependent Navier-Stokes equation. Central difference is used for diffusion terms to achieve better accuracy, and a monotonicity-preserving upwind difference is used for convective ones. The coupled nonlinear system is solved via the damped Newton's method. The Newton Jacobian matrix is formed numerically, and resulting linear system is ill-conditioned and is solved by Bi-CGSTAB with the Gauss-Seidel preconditioner. A novel convection-diffusion-reaction model is introduced to simulate fibroblast growth factor (FGF-2) binding to cell surface molecules of receptor and heparan sulfate proteoglycan and MAP kinase signaling under flow condition. The model includes three parts: the flow of media using compressible Navier-Stokes equation, the transport of FGF-2 using convection-diffusion transport equation, and the local binding and signaling by chemical kinetics. The whole model consists of a set of coupled nonlinear partial differential equations (PDEs) and a set of coupled nonlinear ordinary differential equations (ODEs). To solve the time-dependent PDE system we use second order implicit Euler method by finite volume discretization. The ODE system is stiff and is solved by an ODE solver VODE using backward differencing formulation (BDF). Findings from this study have implications with regard to regulation of heparin-binding growth factors in circulation.

To increase our understanding of how ice sheets and glaciers interact with the climate system, numerical models have become an indispensable tool. However, the complexity of these systems and the natural limitation in computational power is reflected in the simplifications of the represented processes and the spatial and temporal resolution of the models. Whether the effect of these limitations is acceptable or not, can be assessed by theoretical considerations and by validating the output of the models against real world data. Equally important is to verify if the numerical implementation and computational method accurately represent the mathematical description of the processes intended to be simulated. This thesis concerns a set of numerical models used in the field of glaciology, how these are applied and how they relate to other study areas in the same field. The dynamical flow of glaciers, which can be described by a set of non-linear partial differential equations called the Full Stokes equations, is simulated using the finite element method. To reduce the computational cost of the method significantly, it is common to lower the order of the used elements. This results in a loss of stability of the method, but can be remedied by the use of stabilization methods. By numerically studying different stabilization methods and evaluating their suitability, this work contributes to constraining the values of stabilization parameters to be used in ice sheet simulations. Erroneous choices of parameters can lead to oscillations of surface velocities, which affects the long term behavior of the free-surface ice and as a result can have a negative impact on the accuracy of the simulated mass balance of ice sheets. The amount of basal sliding is an important component that affects the overall dynamics of the ice. A part of this thesis considers different implementations of the basal impenetrability condition that accompanies basal sliding, and shows that methods used in literature can lead to a difference in velocity of 1% to 5% between the considered methods. The subglacial hydrological system directly influences the glacier's ability to slide and therefore affects the velocity distribution of the ice. The topology and dominant mode of the hydrological system on the ice sheet scale is, however, ill constrained. A third contribution of this thesis is, using the theory of R-channels to implement a simple numerical model of subglacial water flow, to show the sensitivity of subglacial channels to transient processes and that this limits their possible extent. This insight adds to a cross-disciplinary discussion between the different sub-fields of theoretical, field and paleo-glaciology regarding the characteristics of ice sheet subglacial hydrological systems. In the study, we conclude by emphasizing areas of importance where the sub-fields have yet to unify: the spatial extent of channelized subglacial drainage, to what degree specific processes are connected to geomorphic activity and the differences in spatial and temporal scales. As a whole, the thesis emphasizes the importance of verification of numerical models but also acknowledges the natural limitations of these to represent complex systems. Focusing on keeping numerical ice sheet and glacier models as transparent as possible will benefit end users and facilitate accurate interpretations of the numerical output so it confidently can be used for scientific purposes.<br><p>At the time of the doctoral defense, the following papers were unpublished and had a status as follows: Paper 3: Manuscript. Paper 4: Manuscript.</p><br>Greenland Analogue Project

Les satellites sont vulnérables aux particules de forte énergie piégées dans les ceintures de Van Allen. Afin d’en assurer la protection, il est nécessaire de prédire avec précision la dynamique des électrons au sein de la magnétosphère. Dans un premier temps nous proposons une méthode originale de résolution analytique de l’équation de Fokker-Planck réduite qui modélise le transport et les pertes des électrons de la magnétosphère interne. La résolution repose sur une technique de décomposition spectrale. Si la solution analytique s’avère utile pour mettre en exergue certaines propriétés physiques des ceintures de radiation, elle est également pertinente pour valider le code numérique de résolution de l’équation de Fokker-Planck réduite, développé durant la thèse. Ce dernier nous amène à généraliser l’étude précédente en illustrant l’évolution des flux d’électrons pour diverses énergies et positions. Nous prouvons notamment que la structure des ceintures de radiation ainsi que leur temps d’évolution ne dépendent que de quelques facteurs bien choisis. Dans la perspective de reproduire un événement particulier de retour au calme après un orage magnétique, mesuré par les satellites de la NASA dédiés aux ceintures de radiation, nous sommes en mesure de simuler la précipitation des électrons dans l’atmosphère terrestre causée par les interactions avec les ondes électromagnétiques de la magnétosphère. L’utilisation de conditions bâties sur des données empiriques et spécifiques à la période en question nous permet de corroborer les flux observés. Enfin, l’influence du champ magnétique terrestre sur la dynamique des ceintures de radiation est étudiée sous divers aspects. Nous nous concentrons sur la ceinture externe pour comprendre comment les asymétries du champ magnétique, considérablement façonnées par l’activité solaire, affectent notre manière de concilier théorie et observations. Nous explorons également l’importance de certains processus diffusifs nouveaux et cachés, qui émergent à cause de l’irrégularité naturelle du champ magnétique au plus proche voisinage de la Terre<br>Satellites are vulnerable to high-energy particles trapped in the Van Allen belts. To ensure their protection, it is necessary to predict properly the electron dynamics in the magnetosphere. We first propose an original method to find the analytical solution of the reduced Fokker-Planck equation that models the transport and loss of electrons in the inner magnetosphere. The resolution relies on an eigenfunction expansion approach. If the analytical solution is proven to be useful at uncovering some of the physical properties of the radiation belts, it is also relevant to validate the numerical code that solves the reduced Fokker-Planck equation, which has been developed during the PhD. The latter code is used to generalize the previous study in illustrating the evolution of the electron fluxes for various energies and locations. We demonstrate that the structure of the radiation belts as well as their dynamical timescales only depend on a few well-chosen parameters. In the perspective of reproducing a specific storm-recovery event reported by the NASA Van Allen Probes, we are able to simulate the electron scattering in the Earth’s atmosphere due to the interaction with magnetospheric electromagnetic waves. The consideration of data-driven and event-specific conditions enables us to corroborate the observed fluxes. Finally, various influences of the Earth’s magnetic field on the dynamics of the radiation belts are investigated. We focus on the outer belt to see how the magnetic field asymmetries, which are strongly shaped by solar activity, affect the way of conciliating theory and observations. We also explore the importance of new hidden diffusive processes that emerge due to the natural irregularity of the magnetic field in the closest vicinity of the Earth

Thesis (Ph. D.)--Joint Program in Physical Oceanography (Massachusetts Institute of Technology, Dept. of Earth, Atmospheric, and Planetary Sciences; and the Woods Hole Oceanographic Institution), 1988.<br>Includes bibliographical references (leaves 163-172).<br>A freshwater plume often forms when a river or an estuary discharges water onto the continental shelf. Freshwater plumes are ubiquitous features of the coastal ocean and usually leave a striking signature in the coastal hydrography. The present study combines both hydrographic data and idealized numerical simulations to examine how ambient currents and winds influence the transport and mixing of plume waters. The first portion of the thesis considers the alongshore transport of freshwater using idealized numerical simulations. In the absence of any ambient current, the downstream coastal current only carries a fraction of the discharged fresh water; the remaining fraction recirculates in a continually growing "bulge" of fresh water in the vicinity of the river mouth. The fraction of fresh water transported in the coastal current is dependent on the source conditions at the river mouth. The presence of an ambient current augments the transport in the plume so that its freshwater transport matches the freshwater source. For any ambient current in the same direction as the geostrophic coastal current, the plume will evolve to a steady-state width. A key result is that an external forcing agent is required in order for the entire freshwater volume discharged by a river to be transported as a coastal current. The next section of the thesis addresses the wind-induced advection of a river plume, using hydrographic data collected in the western Gulf of Maine. The observations suggest that the plume's cross-shore structure varies markedly as a function of fluctuations in alongshore wind forcing. Consistent with Ekman dynamics, upwelling favorable winds spread the plume offshore, at times widening it to over 50 km in offshore extent, while downwelling favorable winds narrow the plume width to a few Rossby radii. Near-surface current meters show significant correlations between cross-shore currents and alongshore wind stress, consistent with Ekman theory. Estimates of the terms in the alongshore momentum equation calculated from moored current meter arrays also indicate an approximate Ekman balance within the plume. A significant correlation between alongshore currents and alongshore wind stress suggests that interfacial drag may be important. The final section of the thesis is an investigation of the advection and mixing of a surface-trapped river plume in the presence of an upwelling favorable wind stress, using a three-dimensional model in a simple, rectangular domain. Model simulations demonstrate that the plume thins and is advected offshore by the cross shore Ekman transport. The thinned plume is susceptible to significant mixing due to the vertically sheared horizontal currents. The first order plume response is explained by Ekman dynamics and a Richardson number mixing criterion.<br>by Derek Allen Fong.<br>Ph.D.

Le procédé de dépôt physique en phase vapeur assisté par plasma thermique (PS-PVD) consiste à évaporer le matériau sous forme de poudre à l’aide d’un jet de plasma d’arc soufflé pour produire des dépôts de structures variées obtenus par condensation de la vapeur et/ou dépôt des nano-agrégats. Dans le procédé de PS-PVD classique, l’intégralité du traitement du matériau est réalisée dans une enceinte sous faible pression, ce qui limite les phénomènes d’évaporation ou nécessite d’utiliser des torches de puissance importante. Dans ce travail, une extension du procédé de PS-PVD conventionnel à un procédé à deux enceintes est proposée puis explorée par voie de modélisation et de simulation numérique : la poudre est évaporée dans une enceinte haute pression (105 Pa) reliée par une tuyère de détente à une enceinte de dépôt basse pression (100 ou 1 000 Pa), permettant une évaporation énergétiquement plus efficace de poudre de Zircone Yttriée de granulométrie élevée, tout en utilisant des torches de puissance raisonnable. L’érosion et le colmatage de la tuyère de détente peuvent limiter la faisabilité d’un tel système. Aussi, par la mise en oeuvre de modèles numériques de mécaniquedes fluides et basé sur la théorie cinétique de la nucléation et de la croissance d’agrégats, on montre que, par l’ajustement des dimensions du système et des paramètres opératoires ces deux problèmes peuvent être contournés ou minimisés. En particulier, l’angle de divergence de la tuyère de détente est optimisé pour diminuer le risque de colmatage et obtenir le jet et le dépôt les plus uniformes possibles à l'aide des modèles susmentionnés, associés à un modèle DSMC (Monte-Carlo) du flux de gaz plasmagène raréfié. Pour une pression de 100 Pa, les résultats montrent que la barrière thermique serait formée par condensation de vapeur alors que pour 1 000 Pa, elle serait majoritairement formée par dépôt de nano-agrégats<br>Plasma Spray Physical Vapor Deposition (PS-PVD) aims to substantially evaporate material in powder form by means of a DC plasma jet to produce coatings with various microstructures built by vapor condensation and/or by deposition of nanoclusters. In the conventional PS-PVD process, all the material treatment takes place in a medium vacuum atmosphere, limiting the evaporation process or requiring very high-power torches. In the present work, an extension of conventional PS-PVD process as a two-chamber process is proposed and investigated by means of numerical modeling: the powder is vaporized in a high pressure chamber (105 Pa) connected to the low pressure (100 or 1,000 Pa) deposition chamber by an expansion nozzle, allowing more energetically efficient evaporation of coarse YSZ powders using relatively low power plasma torches. Expansion nozzle erosion and clogging can obstruct the feasibility of such a system. In the present work, through the use of computational fluid dynamics, kinetic nucleation theory and cluster growth equations it is shown through careful adjustment of system dimensions and operating parameters both problems can be avoided or minimized. Divergence angle of the expansion nozzle is optimized to decrease the clogging risk and to reach the most uniform coating and spray characteristics using the aforementioned approaches linked with a DSMC model of the rarefied plasma gas flow. Results show that for 100 Pa, the thermal barrier coating would be mainly built from vapor deposition unlike 1,000 Pa for which it is mainly built by cluster deposition

This doctoral thesis provides the first comprehensive study on the physical processes controlling hydrodynamics and transport in Lake Garda. The investigation is carried out in parallel on three different levels: data collection and analysis, three-dimensional numerical modeling and theoretical study. On the first level, data are collected by building up a network of research institutes and local administrations in the lake area. New data are acquired through traditional field campaigns (CTD, thermistor chains, satellite imagery), while a citizen-science approach, based on local knowledge harvesting, is successfully tested to gather qualitative data on surface circulation. On the second level, a three-dimensional modeling chain is set up, by coupling one-way a mesoscale atmospheric model to a hydrodynamic model. Both models are validated on multiple temporal and spatial scales, allowing to identify the main interactions between the weather forcing and the hydrodynamic response of the lake. Circulations in Lake Garda are found to be very sensitive to the thermal stratification, to the spatial distribution of the wind forcing and to the Earth’s rotation. Surface cyclonic gyre patterns develop in the lake as a residual outcome of alternating wind forcing of local breezes and differential acceleration induced by Earth’s rotation, whereas unidirectional currents flow under a nearly uniform and constant wind. Both model and observations evidences show that, under weak thermal stratification, Ekman transport activates a secondary circulations in the northern part of the lake, driving surface water to the deep layers and possibly preconditioning the lake for subsequent buoyancy-driven deep mixing events. On the third level, the relevance of the Coriolis term in the equations of motion for relatively narrow closed basins is analytically addressed. The classical Ekman problem is solved by including the presence of lateral boundaries and a new analytical solution is formulated. The validity of the new solution is proved by numerical tests of idealized domains of different size, geographical location and turbulent regime, and on Lake Garda as a real test case. The meaningful length scales are discussed, and the significance of Rossby radious as a reference horizontal scale is disproved for steady-state circulations driven by wind and planetary rotation.

Located at the northern edge of the Northern Gulf of Alaska (NGOA), Prince William Sound(PWS) is a small, two-strait semi-enclosed sea. The general ocean circulation pattern inside PWS is significantly affected by external forcings, for instance, the large-scale circulation in NGOA, atmospheric pressure and surface winds, surface heating/cooling, runoff, and tides. Motivated by multi-year experience with a well-validated, quasi-operational ocean circulation nowcast/forecast system for PWS (viz., Extended PWS Nowcast/Forecast System (EPWS/NFS)), the present study addresses some aspects of the PWS response to various external forcings, via numerical simulations. Based on the Princeton Ocean Model (POM), four numerical implementations have been examined, viz., PWS-POM, Extended PWS-POM (EPWS-POM), Idealized PWS-POM (IPWS-POM), and a 2-D tidal model. These implementations are used to simulate physical processes with various spatial and temporal scales in PWS. A series of numerical simulations are conducted, driven by various external forcings ranging from large scale and mesoscale circulation in NGOA represented by the Global Navy Coastal Ocean Model (NCOM), to atmospheric pressure observed by National Data Buoy Center (NDBC) buoys and mesoscale winds predicted by Regional Atmospheric Modeling System (RAMS), and to tides simulated by the 2-D tidal model. These simulations, along with analysis from a Helmholtz resonance model, demonstrate and help interpret some phenomena in PWS; for instance, barotropic Helmholtz resonance in coastal sea levels, and volume transports through the two PWS straits, and a dominant cyclonic gyre in the Central Sound in August and September. The simulation results are used to study a wide range of oceanic phenomena in PWS; e.g., two-layer/three-layer baroclinic transports through the straits, a "transition band" in the coherence pattern between volume transports through the two straits, mesoscale circulation in the Central Sound, the deep water circulation, and the annual tidal energy budget.

Pile-supported T-walls, which are concrete floodwalls that are shaped like an inverted "T" and supported by batter piles, are commonly used by the United States Army Corps of Engineers (USACE) to protect low-lying portions of New Orleans and other areas. The design of a T-wall in southern Louisiana is complex, as the structure needs to resist both 1) large settlements caused by fill placed beneath and beside the T-wall before T-wall construction or by fill placed beside the T-wall after T-wall construction, and 2) large lateral flood loads that are imposed during a hurricane. As a result of these loading conditions, large bending moments can develop in the batter piles and these moments need to be accounted for as part of the T wall design. The goal of this research is to develop a more complete understanding of the pile bending moments in T wall systems, specifically for cross sections where large settlements may occur. As a first step towards this goal, Rensselaer Polytechnic Institute (RPI) performed a series of eight centrifuge tests to investigate and physically model the effects of settlement-induced bending moments on pile-supported T-walls. The centrifuge tests were evaluated and interpreted, and in order to better capture uncertainty, upper and lower bounds were estimated for the interpreted data. The centrifuge results offered some valuable insights on their own, but were ultimately used as the basis for validating and calibrating corresponding numerical models. The numerical models were developed following a rigorous modeling approach and using rational and reasonable assumptions based on widely-accepted and well-justified procedures. The numerical model results were in good agreement with the centrifuge results without the need for significant calibration or modifications. This good agreement indicates that similar numerical models can be developed to reliably analyze actual T-wall cross sections. Detailed recommendations were developed for using numerical models to analyze pile-supported T walls, and an example problem is presented herein that illustrates the application of this approach. These same techniques were then used to perform a parametric study to analyze the combined and separate effects of flood loading for a wide range of different T-wall cross sections. The range was selected in collaboration with the USACE in order to reasonably cover cross sections and conditions that 1) are typically encountered in practice, and 2) were expected to generate both upper and lower bound pile bending moments. In total, 3,648 cross sections were analyzed, and 29,184 sets of analysis results were generated since each cross section was analyzed for eight different loading conditions. Summary results are provided to show the influence of the loading conditions and parameters on T-wall response, including the influence of flood loading, new fill symmetry, pile fixity, number of piles, subsurface profile, pile batter, pile type, levee slope, T-wall elevation, and the presence of existing levee fill. In addition, the key results for all of the analyses are provided in the appendices and in an electronic database. Based on the parametric study results, a simplified analysis procedure was developed that can be used to calculated maximum pile bending moments for T walls installed directly on foundation soils due to settlements. In this procedure, the loads from new fill placed during or after T-wall construction are distributed onto the pile, and the pile response is analyzed using traditional p-y curves and a beam on elastic foundation formulation. This procedure shows good agreement with the numerical model results for a range of conditions. To demonstrate the application of the procedure, the same example problem that is analyzed numerically is reanalyzed using the simplified analysis procedure. Due to the complexity of the problem, it was not possible to modify this procedure or develop a similar procedure for T-walls installed on top of new or existing levees. Overall, this research demonstrates that numerical models can be used to calculate the bending moments that can develop in pile-supported T-walls due to settlements and flood loading, provides valuable insights into the behavior of T-walls and the influence of various parameters on T-wall response, presents a large database of T-wall analysis results, and recommends a simplified analysis procedure that can be used in some cases to calculate pile bending moments due to settlements.<br>Ph. D.

This doctoral thesis provides the first comprehensive study on the physical processes controlling hydrodynamics and transport in Lake Garda. The investigation is carried out in parallel on three different levels: data collection and analysis, three-dimensional numerical modeling and theoretical study. On the first level, data are collected by building up a network of research institutes and local administrations in the lake area. New data are acquired through traditional field campaigns (CTD, thermistor chains, satellite imagery), while a citizen-science approach, based on local knowledge harvesting, is successfully tested to gather qualitative data on surface circulation. On the second level, a three-dimensional modeling chain is set up, by coupling one-way a mesoscale atmospheric model to a hydrodynamic model. Both models are validated on multiple temporal and spatial scales, allowing to identify the main interactions between the weather forcing and the hydrodynamic response of the lake. Circulations in Lake Garda are found to be very sensitive to the thermal stratification, to the spatial distribution of the wind forcing and to the Earth’s rotation. Surface cyclonic gyre patterns develop in the lake as a residual outcome of alternating wind forcing of local breezes and differential acceleration induced by Earth’s rotation, whereas unidirectional currents flow under a nearly uniform and constant wind. Both model and observations evidences show that, under weak thermal stratification, Ekman transport activates a secondary circulations in the northern part of the lake, driving surface water to the deep layers and possibly preconditioning the lake for subsequent buoyancy-driven deep mixing events. On the third level, the relevance of the Coriolis term in the equations of motion for relatively narrow closed basins is analytically addressed. The classical Ekman problem is solved by including the presence of lateral boundaries and a new analytical solution is formulated. The validity of the new solution is proved by numerical tests of idealized domains of different size, geographical location and turbulent regime, and on Lake Garda as a real test case. The meaningful length scales are discussed, and the significance of Rossby radious as a reference horizontal scale is disproved for steady-state circulations driven by wind and planetary rotation.

Questa tesi affronta una varietà di temi di ricerca, spaziando dalla interazione uomo-macchina alla modellizzazione fisica. Ciò che unisce queste ampie aree di interesse è l'idea di utilizzare simulazioni numeriche di fenomeni acustici basate sulla fisica, al fine di implementare interfacce uomo-macchina che offrano feedback sonoro coerente con l'interazione dell'utente. A questo proposito, negli ultimi anni sono nate numerose nuove discipline che vanno sotto il nome di -- per citarne alcune -- auditory display, sonificazione e sonic interaction design. In questa tesi vengono trattate la progettazione e la realizzazione di algoritmi audio efficienti per la sonificazione interattiva. A tale scopo si fa uso di tecniche di modellazione fisica di suoni ecologici (everyday sounds), ovvero suoni che non rientrano nelle famiglie del parlato e dei suoni musicali.<br>The work presented in this thesis ranges over a variety of research topics, spacing from human-computer interaction to physical-modeling. What combines such broad areas of interest is the idea of using physically-based computer simulations of acoustic phenomena in order to provide human-computer interfaces with sound feedback which is consistent with the user interaction. In this regard, recent years have seen the emergence of several new disciplines that go under the name of -- to cite a few -- auditory display, sonification and sonic interaction design. This thesis deals with the design and implementation of efficient sound algorithms for interactive sonification. To this end, the physical modeling of everyday sounds is taken into account, that is sounds not belonging to the families of speech and musical sounds.

[ES] La naturaleza aleatoria que caracteriza algunos fenómenos en sistemas físicos reales (e.g., ingeniería, biología, economía, finanzas, epidemiología y otros) nos ha planteado el desafío de un cambio de paradigma del modelado matemático y el análisis de sistemas dinámicos, y a tratar los fenómenos aleatorios como variables aleatorias o procesos estocásticos. Este enfoque novedoso ha traído como consecuencia nuevas especificidades que la teoría clásica del modelado y análisis de sistemas dinámicos deterministas no ha podido cubrir. Afortunadamente, maravillosas contribuciones, realizadas sobre todo en el último siglo, desde el campo de las matemáticas por científicos como Kolmogorov, Langevin, Lévy, Itô, Stratonovich, sólo por nombrar algunos; han abierto las puertas para un estudio bien fundamentado de la dinámica de sistemas físicos perturbados por ruido. En la presente tesis se discute el uso de ecuaciones diferenciales algebraicas estocásticas (EDAEs) para el modelado de sistemas multifísicos en red afectados por perturbaciones estocásticas, así como la evaluación de su estabilidad asintótica a través de exponentes de Lyapunov (ELs). El estudio está enfocado en EDAEs d-index-1 y su reformulación como ecuaciones diferenciales estocásticas ordinarias (EDEs). Fundamentados en la teoría ergódica, es factible analizar los ELs a través de sistemas dinámicos aleatorios (SDAs) generados por EDEs subyacentes. Una vez garantizada la existencia de ELs bien definidas, hemos procedido al uso de técnicas de simulación numérica para determinar los ELs numéricamente. Hemos implementado métodos numéricos basados en descomposición QR discreta y continua para el cómputo de la matriz de solución fundamental y su uso en el cálculo de los ELs. Las características numéricas y computacionales más relevantes de ambos métodos se ilustran mediante pruebas numéricas. Toda esta investigación sobre el modelado de sistemas con EDAEs y evaluación de su estabilidad a través de ELs calculados numéricamente, tiene una interesante aplicación en ingeniería. Esta es la evaluación de la estabilidad dinámica de sistemas eléctricos de potencia. En el presente trabajo de investigación, implementamos nuestros métodos numéricos basados en descomposición QR para el test de estabilidad dinámica en dos modelos de sistemas eléctricos de potencia de una-máquina bus-infinito (OMBI) afectados por diferentes perturbaciones ruidosas. El análisis en pequeña-señal evidencia el potencial de las técnicas propuestas en aplicaciones de ingeniería.<br>[CA] La naturalesa aleatòria que caracteritza alguns fenòmens en sistemes físics reals (e.g., enginyeria, biologia, economia, finances, epidemiologia i uns altres) ens ha plantejat el desafiament d'un canvi de paradigma del modelatge matemàtic i l'anàlisi de sistemes dinàmics, i a tractar els fenòmens aleatoris com a variables aleatòries o processos estocàstics. Aquest enfocament nou ha portat com a conseqüència noves especificitats que la teoria clàssica del modelatge i anàlisi de sistemes dinàmics deterministes no ha pogut cobrir. Afortunadament, meravelloses contribucions, realitzades sobretot en l'últim segle, des del camp de les matemàtiques per científics com Kolmogorov, Langevin, Lévy, Itô, Stratonovich, només per nomenar alguns; han obert les portes per a un estudi ben fonamentat de la dinàmica de sistemes físics pertorbats per soroll. En la present tesi es discuteix l'ús d'equacions diferencials algebraiques estocàstiques (EDAEs) per al modelatge de sistemes multifísicos en xarxa afectats per pertorbacions estocàstiques, així com l'avaluació de la seua estabilitat asimptòtica a través d'exponents de Lyapunov (ELs). L'estudi està enfocat en EDAEs d-index-1 i la seua reformulació com a equacions diferencials estocàstiques ordinàries (EDEs). Fonamentats en la teoria ergòdica, és factible analitzar els ELs a través de sistemes dinàmics aleatoris (SDAs) generats per EDEs subjacents. Una vegada garantida l'existència d'ELs ben definides, hem procedit a l'ús de tècniques de simulació numèrica per a determinar els ELs numèricament. Hem implementat mètodes numèrics basats en descomposició QR discreta i contínua per al còmput de la matriu de solució fonamental i el seu ús en el càlcul dels ELs. Les característiques numèriques i computacionals més rellevants de tots dos mètodes s'illustren mitjançant proves numèriques. Tota aquesta investigació sobre el modelatge de sistemes amb EDAEs i avaluació de la seua estabilitat a través d'ELs calculats numèricament, té una interessant aplicació en enginyeria. Aquesta és l'avaluació de l'estabilitat dinàmica de sistemes elèctrics de potència. En el present treball de recerca, implementem els nostres mètodes numèrics basats en descomposició QR per al test d'estabilitat dinàmica en dos models de sistemes elèctrics de potència d'una-màquina bus-infinit (OMBI) afectats per diferents pertorbacions sorolloses. L'anàlisi en xicotet-senyal evidencia el potencial de les tècniques proposades en aplicacions d'enginyeria.<br>[EN] The random nature that characterizes some phenomena in the real-world physical systems (e.g., engineering, biology, economics, finance, epidemiology, and others) has posed the challenge of changing the modeling and analysis paradigm and treat these phenomena as random variables or stochastic processes. Consequently, this novel approach has brought new specificities that the classical theory of modeling and analysis for deterministic dynamical systems cannot cover. Fortunately, stunning contributions made overall in the last century from the mathematics field by scientists such as Kolmogorov, Langevin, Lévy, Itô, Stratonovich, to name a few; have opened avenues for a well-founded study of the dynamics in physical systems perturbed by noise. In the present thesis, we discuss stochastic differential-algebraic equations (SDAEs) for modeling multi-physical network systems under stochastic disturbances, and their asymptotic stability assessment via Lyapunov exponents (LEs). We focus on d-index-1 SDAEs and their reformulation as ordinary stochastic differential equations (SDEs). Supported by the ergodic theory, it is feasible to analyze the LEs via the random dynamical system (RDSs) generated by the underlying SDEs. Once the existence of well-defined LEs is guaranteed, we proceed to the use of numerical simulation techniques to determine the LEs numerically. Discrete and continuous QR decomposition-based numerical methods are implemented to compute the fundamental solution matrix and use it in the computation of the LEs. Important numerical and computational features of both methods are illustrated through numerical tests. All this investigation concerning systems modeling through SDAEs and their stability assessment via computed LEs finds an appealing engineering application in the dynamic stability assessment of power systems. In this research work, we implement our QR-based numerical methods for testing the dynamic stability in two types of single-machine infinite-bus (SMIB) power system models perturbed by different noisy disturbances. The analysis in small-signal evidences the potential of the proposed techniques in engineering applications.<br>Mi agradecimiento al estado ecuatoriano que, a través del Programa de Becas para el Fortalecimiento y Desarrollo del Talento Humano en Ciencia y Tecnología 2012 de la Secretaría Nacional de Educación Superior, Ciencia y Tecnología (SENESCYT), han financiado mis estudios de doctorado.<br>González Zumba, JA. (2020). Dynamic Modeling and Stability Analysis of Stochastic Multi-Physical Systems Applied to Electric Power Systems [Tesis doctoral]. Universitat Politècnica de València. https://doi.org/10.4995/Thesis/10251/158558<br>TESIS

Cette étude porte sur les écoulements de gaz à faible vitesse dans des réseaux de tuyaux, en se concentrant sur un régime caractérisé par un faible nombre de Mach. Nous utilisons un modèle unidimensionnel obtenu en moyennant les équations de Navier-Stokes pour un fluide compressible sur la section d'une conduite. Notre approche utilise le développement asymptotique à faible nombre de Mach pour décrire les effets compressibles de manière plus précise, contrairement à l'approximation classique de Boussinesq qui en est un cas limite lorsque l'élévation de température est assez faible. Nous utilisons un schéma numérique fondé sur la méthode des caractéristiques et la méthode de projection pour traiter ce modèle à faible Mach. Nous présentons des résultats numériques pour une configuration appelée "thermosiphon". Cette configuration consiste en une boucle fermée constituée de deux tuyaux horizontaux adiabatiques et de deux tuyaux verticaux avec des températures de paroi prescrites, qui entraînent l'écoulement. L'algorithme développé permet de prendre en compte des distributions de Dirac, qui peuvent apparaître en terme source dans le modèle pour représenter les coins de la géométrie. La méthode proposée est également adaptée au modèle pour des conditions aux limites de type périodique en plus des conditions de type Dirichlet. Nous incorporons dans notre algorithme le traitement des distributions de Dirac en tant que dérivées du terme de gravité discontinu aux coins et des conditions périodiques. Afin d'établir une référence pour le problème du thermosiphon, nous fournissons une solution exacte mais semi-implicite d'un écoulement laminaire en régime permanent. Cette solution sert de référence pour valider la méthode proposée. Nous proposons également des lois qui régissent les jonctions entre plusieurs conduites et présentons des résultats numériques pour des configurations de conduites plus complexes. Nous nous intéressons en particulier à l'échelle à trois barreaux, une configuration fermée composée de tuyaux horizontaux adiabatiques et de tuyaux verticaux dont les parois sont soumises à des températures imposées afin d'induire un gradient de température qui alimente l'écoulement. qui corresponds à une extension du thermosiphon. À partir de celle-ci, plusieurs configurations générales peuvent être dérivées. Nous développons un algorithme pour garantir des conditions de transmission adéquates aux jonctions, en fournissant autant que possible des solutions stationnaires semi-analytiques pour valider nos résultats numériques. Cette étude contribue à une meilleure compréhension des principes qui régissent les écoulements de gaz à faible nombre de Mach en utilisant des techniques numériques avancées et en les comparant à des références établies<br>This project aims to develop a model for low Mach flow in pipelines and an industrial code implementing it. The model can describe low mach regimes while avoiding blunt approximations, improving over legacy approaches like Boussinesq; as a result, our program is more accurate. In order to build our model and program, we investigate gas flow at low velocities in a network of pipes. We consider a one-dimensional system of equations obtained by averaging the Navier-Stokes equations for a compressible fluid over the pipe section. In contrast to the classical Boussinesq approximation, we employ the Low Mach Expansion to describe asymptotically compressible effects, aiming for a more accurate solution capable of characterizing flows with significant temperature variations. We first apply the model we obtained thus far to a well-known configuration of pipes called the "thermosyphon.". This setup consists of a loop of two horizontal adiabatic pipes and two vertical pipes with prescribed wall temperatures, resulting in a temperature gradient that drives the flow. The application of the model to this configuration gives us an exact but semi-implicit solution under laminar and steady-state conditions. This solution serves as a benchmark against which we validate our numerical results. By comparing our computed values with the quasi-exact solution, we demonstrate the accuracy and reliability of our approach. To implement the low Mach averaged model, we use a numerical method based on the characteristics method and the projection technique. We incorporate in our algorithm the treatment of periodic conditions and Dirac distributions as derivatives of the discontinuous gravity term at the corners. To generalize the model to more complex configurations of pipes, we propose laws that govern the junctions between multiple pipes. We study the "three-rung ladder," a closed configuration consisting of horizontal adiabatic pipes and vertical ones with imposed wall temperatures to induce a temperature-driven flow. %From this setup, more complex pipeline configurations can be derived. To tackle the challenges junctions pose in this context, we implemented an algorithm in the program capable of ensuring proper transmission conditions. Whenever feasible, we provide quasi-exact solutions under laminar and steady-state conditions to validate our numerical results further. Overall, this study investigates further low Mach number gas flows, employing advanced numerical techniques and validating our findings against established benchmarks

In high frequency computational electromagnetics, rigorous numerical methods be come unrealistic tools due to computational demand increasing with the frequency. Instead approximations to the solutions of the Maxwell equations can be employed to evaluate th electromagnetic fields. In this thesis, we present the implementations of three high frequency approximat methods. The first two, namely the Geometrical Theory of Diffraction (GTD) and th Physical Optics (PO), are commonly used approximations. The third is a new invention that will be referred to as the Surface Current Extraction-Extrapolation (SCEE). Specifically, the GTD solver is a flexible and modular software package which use Non-Uniform Rational B-spline (NURBS) surfaces to model complex geometries. The PO solver is based on a triangular description of the surfaces and includes fas shadowing by ray tracing as well as contribution from edges to the scattered fields. GTD ray tracing was combined with the PO solver by a well thought-out software architecture Both implementations are now part of the GEMS software suite, the General ElectroMag netic Solvers, which incorporates state-of-the-art numerical methods. During validations both GTD and PO techniques turned out not to be accurate enough to meet the indus trial standards, thus creating the need for a new fast approximate method providing bette control of the approximations. In the SCEE approach, we construct high frequency approximate surface currents ex trapolated from rigourous Method of Moments (MoM) models at lower frequency. T do so, the low frequency currents are projected onto special basis vectors defined on th surface relative to the direction of the incident magnetic field. In such configuration, w observe that each component displays systematic spatial patterns evolving over frequenc in close correlation with the incident magnetic field, thus allowing us to formulate a fre quency model for each component. This new approach is fast, provides good control of th error and represents a platform for future development of high frequency approximations. As an application, we have used these tools to analyse the radar detectability of a new marine distress signaling device. The device, called "Rescue-Wing", works as an inflatabl radar reflector designed to provide a strong radar echo useful for detection and positionin during rescue operations of persons missing at sea.<br>QC 20101004

The focus of this thesis was to investigate which baseball metrics are most conducive to run creation and prevention. Stepwise regression and Liu estimation were used to formulate two models for the dependent variables and also used for cross validation. Finally, the predicted values were fed into the Pythagorean Expectation formula to predict a team’s most important goal: winning. Each model fit strongly and collinearity amongst offensive predictors was considered using variance inflation factors. Hits, walks, and home runs allowed, infield putouts, errors, defense-independent earned run average ratio, defensive efficiency ratio, saves, runners left on base, shutouts, and walks per nine innings were significant defensive predictors. Doubles, home runs, walks, batting average, and runners left on base were significant offensive regressors. Both models produced error rates below 3% for run prediction and together they did an excellent job of estimating a team’s per-season win ratio.

Travaux portant sur le moussage physique par adjonction de CO2 supercritique de polyoléfines et thermoplastiques vulcanisés (TPV) que ce soit par procédé continu (extrusion) ou par procédé discontinu (réacteur).Grâce à la formulation, l’étude rhéologique et le moussage de polymères existants et inédits un critère rhéologique permettant de définir la capacité à mousser a été démontré et étendu au mélange complexe de polymère (TPV). Plus spécifiquement l’étude de la formulation des TPV sur le moussage a été investigué et des conclusions concernant la chimie de réticulation ainsi que l’ajout de charges ont pu être déterminées.La modélisation numérique et la mise en œuvre de l’extrusion moussage ont été réalisées et ont permis de démontrer la possibilité d’obtenir une mousse régulière et homogène ainsi que de prédire, en fonction des paramètres d’entrés, les caractéristiques finales de la mousse<br>Work on the physical foaming of vulcanised polyolefins and thermoplastics (TPV) by the addition of supercritical CO2, whether by continuous process (extrusion) or by batch process (reactor).Through the formulation, rheological study and foaming of existing and new polymers, a rheological criterion for defining the ability to foam was demonstrated and extended to complex polymer blends (as TPV). More specifically, the study of TPV formulation on foaming was investigated .Conclusions regarding cross-linking chemistry and addition of fillers were determined.Numerical modelling and implementation of extrusion foaming were carried out and demonstrated the possibility of obtaining a regular and homogeneous foam as well as predicting, depending on the input parameters, the final characteristics of the foam

Cette these propose une modelisation et une simulation d'un processus metallurgique de mise en forme dans un but de commande. Une description du procede en question est d'abord exposee. Puis nous analysons la commande et nous mettons en evidence l'interet de modeliser le procede. Notre volonte a ete d'etablir un modele physique simplifie : d'une part le modele doit contenir assez de physique pour faire de la prediction, d'autre part le modele doit etre exploitable et sa parametrisation compatible avec les donnees disponibles. Nous avons decompose le processus metallurgique de mise en forme en trois modeles. La conception de ces trois modeles est alors presentee, chaque modele faisant l'objet d'une etude particuliere. Les deux premiers modeles s'appuient essentiellement sur la theorie hydraulique, donc sur la theorie physique du processus, alors que le troisieme est davantage une representation simplifiee. Enfin cette modelisation est veritablement validee a partir de donnees de fabrication, dans un but d'utilisation en temps reel. Un logiciel de simulation a ete cree a partir de ce modele. La modelisation et la simulation ont permis d'obtenir une representation conceptuelle, ainsi qu'une meilleure connaissance du processus metallurgique de mise en forme.

In dieser Arbeit wird anfangs dargestellt, wie kostenintensiv und problematisch die Kanalreinigung ausfällt. Durch die demographisch bedingt rückläufigen Schmutzwassermengen und gleichzeitig prognostizierten langen Trockenperioden in den Sommermonaten muss allerdings mit einer verstärkten Akkumulation von Ablagerungen in den Kanalsystemen gerechnet werden. Zudem erweist es sich in Voraussicht auf den Klimawandel als notwendig, den gesamten Kanalquerschnitt für Starkregenereignisse von Ablagerungen freizuhalten; ansonsten ist die häufigere Entlastung des Abwassernetzes in den Vorfluter nach Füllung der vorhandenen Regenrückhaltebecken vorprogrammiert. Die Kanalreinigung wird somit in den nächsten Jahren an Bedeutung gewinnen. Vor allem bei der Reinigung von Großprofilen ist jedoch der Stand der Technik noch nicht ausgereift und die jetzigen Reinigungsmethoden erweisen sich wegen des niedrigen Reinigungserfolges sowie der starken Umweltbelastung durch die genutzten Geräte als ungeeignet. Durch die in dieser Arbeit vorgestellte und auf den Stand der Technik gebrachte Reinigungsmethode wird jedoch eine leistungsstarke Alternative geboten. In Zukunft wird es möglich sein, große Kanaldurchmesser fast vollkommen umweltneutral zu reinigen und nur mit Hilfe der Wasserkraft das gesammelte Sediment bis zum nachgeschalteten Sandfang zu mobilisieren. Mit den durchgeführten und hier dargestellten numerischen Versuchen wird in erster Linie eine geometrische Optimierung zur Erreichung der größten Mobilisierungskräfte im Hinblick auf einen ausgewählten Abschnitt der Dresdner Kanalisation vollzogen. In der nachfolgenden physikalischen Modellierung wird der Schwerpunkt auf die optimale Reinigungsgeschwindigkeit gelegt, um ein möglichst großes Volumen an Ablagerungen zu bewegen. Es muss allerdings präzisiert werden, dass während der Modellierung kein Grenzzustand erreicht werden konnte. Dieser wird schließlich durch die folgende analytische Beschreibung des Transportvorganges theoretisch bestimmt. Für das genutzte physikalische Modell wird analytisch ein Grenztransportkörper aus Modellsediment errechnet. Bei einer Übertragung auf den dazugehörigen vorhandenen Abschnitt im Dresdner Kanalnetz durch Ansetzen des Ähnlichkeitsgesetzes wird ein mobilisierbares Grenzvolumen überschlagen. Die hier dargestellten Ergebnisse der theoretischen Untersuchungen bzw. der physikalischen Modellierung wurden mit den Betreibern und den Herstellern des neuen Stauwagens diskutiert und hinsichtlich konstruktiver Schwerpunkte optimiert. Unter besonderer Berücksichtigung der Kostenoptimierung und der konstruktiven Realisierung wurde von der Firma "WKS Technik GmbH" ein Forschungsmuster zur Untersuchung in einer Pilottestphase entwickelt und gebaut. In der nachfolgenden dargestellten Testphase wurde das Ergebnis als Pilotprojekt im Dresdner Kanalnetz in mehreren Testläufen erprobt und untersucht. Nach jedem Testlauf wurden die aufgetretenen Probleme analysiert und durch Umbauten oder Veränderungen behoben. Regelmäßige Reinigungsfahrten sollten wegen der geringeren Ablagerungsmengen ein gleichmäßiges Fahren gewährleisten bzw. die Ablagerungen vollständig mobilisieren und bis zum Sandfang transportieren, um sie dort zu entnehmen. Während der Testphase wurde eine viel höhere Sedimentmenge gesammelt und transportiert als die, welche sich aus der analytischen Berechnung ergab. Dies folgte aus dem beim Naturversuch zugelassenen höheren Aufstauniveau hinter dem Spülwagen, welches ausnahmsweise zugelassen wurde, um auf eine zwischenzeitliche Sedimentabsaugung zu verzichten<br>The early pages of this research thesis demonstrate how costly and problematic sewer cleaning is. Due to changing demographics causing declining amounts of wastewater and, at the same time, to predicted long dry spells during the summer months, an increased accumulation of sediments in the sewer systems has to be expected. With respect to climate change it is necessary to keep the entire sewer cross section free of sediments in order to manage heavy rainfall events; otherwise a more frequent relief of the wastewater system into the receiving water is inevitable once the rain retention basins have been filled. Thus sewer cleaning will gain in importance within the next few years. Predominantly the cleaning of large-diameter profiles, however, has not reached a state-of-the-art technology, and the current cleaning methods prove themselves to be inapt due to their very limited success as well as to the environmental burden posed by the equipment they employ. The state-of-the-art cleaning method presented in this research paper constitutes a high-powered alternative. In the future, it will be possible to clean large sewer cross sections almost entirely in an environmentally neutral fashion and to mobilize the accumulated sediment all the way to the downstream sand trap by exclusively harnessing water power. The numerical experiments conducted for this paper primarily effected a geometric optimization to achieve the maximal mobilization forces in a selected section of Dresden's sewer system. In the subsequent physical modeling the emphasis is placed on the optimal cleaning speed in order to move a maximum volume of sediments. It has to be clarified, however, that no limit state could be reached during the modeling procedure. This limit state is determined theoretically only through the following analytical description of the transport process. For the physical model at hand, the maximum transportable sediment volume made of model sediment is calculated analytically. In an extrapolation for the associated existent section in the Dresden sewer system, a mobilizable limit volume is estimated applying the law of similarity. The results of the theoretical analysis and the physical modeling presented in this work were discussed with the operators and manufacturers of the new storage vehicle and were optimized with regard to constructional core aspects. In special consideration of the cost optimization and the constructional realization, the "WKS Technik GmbH" company developed and constructed a research prototype for a survey in a pilot test phase. During the subsequent test phase, the result was tried out and examined as a pilot project in several trial runs in Dresden's sewer system. After each trial run, the problems which had occurred were analyzed and then eliminated through constructional modifications or other alterations. Regular cleaning runs were intended to guarantee a steady movement of the cleaning vehicle due to a reduced amount of sediments, or to completely mobilize the sediments and transport them to the sand trap in order to remove them there. The sediment volume actually gathered and transported during the test phase was significantly larger than the one derived from the analytical calculation. This increase results from greater storage depths behind the flush car, which were permitted as an exception for the field test in order to go without any interim sediment extractions<br>La prima parte della presente tesi di dottorato mette in evidenza le problematicità legate alla pulizia dei canali delle acque reflue, sottolineando in particolare gli enormi costi generati dall 'alta frequenza di depositi e dalla grandezza delle reti di raccolta presenti sui territori. A causa della riduzione demografica e delle condizioni meteorologiche che si stanno alterando, per il prossimo futuro si prevedono lunghi periodi di siccità, interrotti da acquazzoni intensi e brevi. Questo porterà ad un accumulo di depositi nelle reti di raccolta delle acque miste soprattutto nei mesi estivi, che implicherà a sua volta la necessità di tenere le sezioni dei canali presenti libere da depositi per consentire un appropriato deflusso e ritenzione in presenza delle suddette intense precipitazioni. Se in questi casi la ritenzione della rete fognaria non sarebbe data, infatti, si andrebbe a scaricare acque sporche sempre in maggiore quantità e con maggiore frequenza nei corsi d 'acqua di superficie, inquinando in secondo luogo anche le falde acquifere. Da questo quadro, si evince come la pulizia delle canalizzazioni acquisterà in futuro sempre una maggiore importanza. Attualmente, per la pulizia di canali di grandi dimensioni (maggiori di 1000 mm) non esistono metodi efficaci ed ecologici per movimentare i depositi presenti. Il presente lavoro espone una soluzione competitiva ai suddetti problemi. Grazie al lavoro di ricerca e sviluppo svolto, infatti, è stato elaborato un metodo con cui sarà possibile in futuro pulire canali di grosse dimensioni pressoché senza l\'impiego di energia esterna, ma solamente grazie all 'utilizzo dell 'energia dell 'acqua reflua presente nei canali, mobilizzando i depositi fino al prossimo punto presente dove sarà possibile estrarre dalla rete fognaria il materiale raccolto. L 'analisi numerica eseguita e qui rappresentata porta innanzitutto ad un 'ottimizzazione geometrica ed idraulica che consente il raggiungimento delle maggiori forze di mobilizzazione possibili per un segmento scelto della rete fognaria della città di Dresda (Sassonia). Nelle prove in laboratorio correlate si è cercato di ottimizzare la velocità con la quale deve avvenire la pulizia del fondo del canale in modo da poter trasportare il più elevato volume di sedimenti possibile. Si precisa a questo riguardo che in laboratorio non è stato possibile raggiungere il limite volumetrico correlato al modello fisico utilizzato, perché i soli 24 m di lunghezza del banco di prova non hanno consentito il raggiungimento di tale stato limite. Questo limite è stato determinato successivamente attraverso un modello analitico sviluppato appositamente, il quale descrive la mobilizzazione dei sedimenti attraverso la tecnologia di pulizia utilizzata. Attraverso i fattori di conversione che regolano le leggi fisiche presenti tra modello di laboratorio e natura è stato possibile in seguito rapportare il limite volumetrico al segmento reale utilizzato per le prove numeriche. Usando i risultati ottenuti sia dalle analisi numeriche che dalle prove di laboratorio è stato sviluppato e poi realizzato un prototipo assieme al partner commerciale "WKS Technik GmbH", prestando particolare attenzione all 'ottimizzazione costruttiva e dei costi di produzione. Questo prototipo è stato testato per diverse volte in un tratto di circa quattro chilometri della canalizzazione di Dresda. Dopo ogni corsa sono stati analizzati i punti deboli del metodo di pulizia nonché del prototipo stesso e si è provveduto ad eliminarli. Gli ultimi test hanno confermato il funzionamento di tale prototipo offrendo inoltre una buona validazione del modello analitico sviluppato in precedenza

Static and dynamic interaction mechanism of the retained soil-compressible geofoam buffer and yielding retaining structures requires further investigation. The present study, initiated on this motive, discusses the results of 1-g physical model tests and numerical analyses of cantilever retaining walls with and without deformable geofoam buffers between the wall and cohesionless granular backfill. 0.7m high walls with various wall thicknesses were utilized in the physical modeling. Dynamic tests were carried out by using a laminar container placed on a uni-axial shaking table. Influence of buffer thickness, geofoam type and wall flexibility as well as base excitation characteristics on the lateral earth pressures and flexural wall deflections were under concern. Outcomes of the analyses performed with FLAC-2D (v6.0) finite difference code were validated against the results of the physical model tests. It was observed that the arching effect induced in the retained soil by the lateral compression of the lower half of the geofoam buffer has a positive effect, as this zone is able to absorb a portion of the total unbalanced lateral force exerted by the backfill thus causing a reduction in the static and seismic lateral wall pressures. Relative thickness and stiffness of the geofoam buffer appear to be the most dominant factors affecting the reduction in earth thrust. Lateral earth pressure coefficients determined from physical model tests were compared with those calculated using methods available in the literature. Good agreement was observed between the predictions. Graphs were provided to estimate the static and dynamic lateral earth pressure coefficients for various combinations of wall stiffness and buffer characteristics. Analysis of a 6m high prototype cantilever wall subjected to an excitation recorded in August 17, 1999 Kocaeli earthquake by finite difference method exhibited the contribution of geofoam buffers on seismic performance of cantilever earth retaining walls. It was observed that the presence of an EPS geofoam inclusion provides a reduction of the permanent flexural wall deflections as well as total seismic thrust likely to be experienced by the wall during an earthquake.

Cette thèse de doctorat est dédiée au développement de modèles physiques pour l’étude des systèmes d’aspersion de gouttelettes d’eau en milieu réactif d’hydrogène-air pré-mélangée dans les centrales nucléaires. Des modèles d’ordre réduit sont développés pour décrire l’évaporation des gouttelettes d’eau dans la flamme, la dispersion des nuages de particules après le passage des ondes de choc et l’évolution de l’échelle caractéristiques de turbulence avec la présence d’un jet d’eau. Une nouvelle méthodologie est proposée pour évaluer les effets de l’évaporation par l’aspersion sur la propagation de la flamme d’hydrogène turbulente à l’intérieur d’un volume fermé et un modèle simple est développé pour la quantification de la décélération de la vitesse laminaire avec l’évaporation des gouttelettes à l’intérieur de la flamme. Également, un modèle analytique est proposé pour la prédiction de la dispersion de nuage de particule après le passage d’une onde de choc en s’appuyant sur le one-way formalisme avec une extension afin de prédire l’apparition d’un pic de densité du nombre de particules en utilisant le two-way formalisme. En ce qui concerne la modulation de la turbulence induite par les particules, un modèle simple est utilisé pour l’estimation des échelles intégrales de la turbulence induites par l’injection de nuage des particules. Ces modèles numériques développés peuvent être couplés pour être mis en œuvre dans les simulations numériques à grande échelle de l’effet du système d’aspersion sur les explosions accidentelles d’hydrogène dans les centrales nucléaires<br>This PhD dissertation is dedicated to develop simple models to investigate the effect of water spray system on the premixed hydrogen-air combustion in the nuclear power plants. Specific simple models are developed to describe the water droplet evaporation in the flame, particle cloud dispersion after the shock wave passage, and turbulence length scale evolution with the presence of a water spray. A methodology is proposed to evaluate the spray evaporation effects on the propagation of the turbulent hydrogen flame inside a closed volume and a simple model is developed for the quantification of the laminar velocity deceleration with the droplets evaporation inside the flame. An analytical model is proposed for the prediction of particle cloud dispersion after the shock passage in the one-way formalism and another analytical model is dedicated to describe the spray-shock interaction mechanism and predict the appearance of a particle number density peak using the two-way formalism. A review of the important criteria and physical modelings related to the particle-induced turbulence modulation is given and a mechanistic model is used for the estimation of the turbulent integral length scales induced by the injection of particle clouds. These developed numerical models can be coupled to implement in the large-scale numerical simulations of the spray system effects on the accidental hydrogen explosions in the nuclear power plants

This work deals with the design and calculation of membrane structure for shelter placed above selected ground plan. The hyperbolic paraboloid was chosen as most optimal shape. Subsequently, three different sub-studies were conducted. In the first study, influence of mesh density on a finite element calculation model was studied. Second study deal with the effect of the geometry changes of the computational model. In the third study, deformations and tensions of membrane structure in relation to material formation calculation model was studied. For the static solution of those selected variant, structural analysis software RFEM and RF-FORM-FINDING add-on module was used. Geometric nonlinear solution was applied in the calculation.

In recent years, Polymer Electrolyte Membrane Fuel Cells (PEMFC) have attracted much attention due to their potential as a clean power source for many applications, including automotive, portable and stationary devices. This resulted in a tremendous technological progress, such as the development of new membranes and electro-catalysts or the improvement of electrode structures. However, in order to compete within the most attractive markets, the PEMFC technologies did not reach all the required characteristics yet, in particular in terms of cost and durability.Because of the strong coupling between different physicochemical phenomena, the interpretation of experimental observations is difficult, and analysis through modeling becomes crucial to elucidate the degradation and failure mechanisms, andto help improving both PEMFC electrochemical performance and durability.The development of a theoretical tool is essential for industrials and the scientific community to evaluate the PEMFC degradation and to predict itsperformance and durability in function of the materials properties and in a diversity of operating conditions. This manuscript summarizes my scientific research efforts in this exciting topic during the last 9 years in France, including my invention of the MEMEPhys multiscale simulation package,developed on the basis of my childhood passion for the New Technologies for Energyin Argentina. My perspectives of adapting this approach to other electrochemical systems such as water electrolyzers and batteries are also discussed.

Дисертація на здобуття наукового ступеня кандидата технічних наук за спеціальністю 05.11.13 – прилади і методи контролю та визначення складу речовин. – Національний технічний університет "Харківський політехнічний інститут", м. Харків, 2019 р. У дисертації вирішено важливе наукове завдання розвитку методів теплового контролю технічного стану індукційних установок виробництва високоякісної мідної катанки на основі використання математичного моделювання електротеплових процесів, аналізу розподілу температури в багатошаровій термоізоляції установок з урахуванням її деградації та розподілу температури на їх поверхні, що має суттєве значення для розвитку кабельної промисловості. Розроблено новий метод неруйнівного контролю технічного стану індукційних установок для виготовлення мідної катанки енергетичного призначення, в якому крім вимірювання реактивного та активного опорів індуктора та різниці температури води у системі охолодження додатково визначається ступінь деградації багато-шарової термоізоляції шляхом порівняння розрахункового розподілу поточної температури в об'ємі розплаву та ізоляції на тривимірній математичній моделі установки з практичним вимірюванням температури локальних областей (у роботі їх було 72) поверхні корпусу печі. Це дає можливість точніше прогнозувати ресурс печі та завчасно виявляти режими її роботи близькі до аварійних. Розроблено нове технічне рішення з удосконалення багатошарової термоізоляції індукційних установок з метою зменшення їхніх енерговитрат і збільшення ресурсу. Воно полягає у застосуванні замість третього термоізоляційного шару з легковагої цегли трьох шарів: монолітного вогнетривкого бетону, легковагої цегли та вогнетривкого паперу, зі збереженням загальної товщини ізоляції. Впровадження результатів дисертації в індукційній установці UPCAST US20X-10 для виготовлення мідної катанки в ПАТ "ЗАВОД ПІВДЕНКАБЕЛЬ" дало позитивні результати у зниженні на 20°C температури на корпусі печі, скороченні на 5-15% питомої витрати електроенергії на 1 тонну катанки та поліпшенні її якості шляхом зменшення відсотку сторонніх домішок і забезпечення вмісту в ній міді 99,99%.<br>The thesis for a Candidate of Engineering Sciences degree by specialty 05.11.13 – devices and methods of control and determination of composition of substances. – National technical university "Kharkov polytechnic institute", Kharkiv, 2019. The thesis solves the important scientific task to improve methods of control of the technical condition of induction installations for the production of high-quality copper rod using the mathematical modeling of electrothermal processes in copper melt and multi-layer thermal insulation, with taking into account its degradation and temperature distribution on the installation surface, which is essential for the development of cable industry. A new method of non-destructive control of the technical condition of induction installations for the production of copper rod for power cables of energy value is developed, in which in addition to measuring the reactive and active resistances of the inductor and the reduction of the water temperature in the cooling system, the degree of degradation of multilayer thermal insulation is additionally determined by comparing the calculated distribution of the current temperature in the volume of melt and insulation on a three-dimensional mathematical model of the installation with practical measurement of temperature in local areas (in the work they were 72) of the surface of the furnace casing. This makes it possible to more accurately predict the furnace's profile and to detect the modes of its operation in a timely manner close to the emergency. A new technical solution for improving the multilayer thermoinsulation of induction installations has been developed in order to reduce their energy consumption and increase the resource. It consists in applying instead of the third thermal insulation layer of light-weight brick of three layers: monolithic refractory concrete, lightweight brick and refractory paper, with the preservation of the overall thickness of the insulation. The implementa-tion of the developed lining structure improvement at PJSC " YUZHCABLE WORKS" in the UPCAST US20X-10 induction installation for copper rod in the cable industry in Ukraine has yielded positive results in a 20°C decrease in the temperature of the furnace casing, a reduction of 5-15% of the specific cost Electricity for 1 ton of production and improvement of quality of copper rod, by reducing the percentage of foreign impurities and ensuring copper content in it 99,99%.

Дисертація на здобуття наукового ступеня кандидата технічних наук за спеціальністю 05.11.13 – прилади і методи контролю та визначення складу речовин. – Національний технічний університет "Харківський політехнічний інститут", м. Харків, 2019 р. У дисертації вирішено важливе наукове завдання розвитку методів теплового контролю технічного стану індукційних установок виробництва високоякісної мідної катанки на основі використання математичного моделювання електротеплових процесів, аналізу розподілу температури в багатошаровій термоізоляції установок з урахуванням її деградації та розподілу температури на їх поверхні, що має суттєве значення для розвитку кабельної промисловості. Розроблено новий метод неруйнівного контролю технічного стану індукційних установок для виготовлення мідної катанки енергетичного призначення, в якому крім вимірювання реактивного та активного опорів індуктора та різниці температури води у системі охолодження додатково визначається ступінь деградації багато-шарової термоізоляції шляхом порівняння розрахункового розподілу поточної температури в об'ємі розплаву та ізоляції на тривимірній математичній моделі установки з практичним вимірюванням температури локальних областей (у роботі їх було 72) поверхні корпусу печі. Це дає можливість точніше прогнозувати ресурс печі та завчасно виявляти режими її роботи близькі до аварійних. Розроблено нове технічне рішення з удосконалення багатошарової термоізоляції індукційних установок з метою зменшення їхніх енерговитрат і збільшення ресурсу. Воно полягає у застосуванні замість третього термоізоляційного шару з легковагої цегли трьох шарів: монолітного вогнетривкого бетону, легковагої цегли та вогнетривкого паперу, зі збереженням загальної товщини ізоляції. Впровадження результатів дисертації в індукційній установці UPCAST US20X-10 для виготовлення мідної катанки в ПАТ "ЗАВОД ПІВДЕНКАБЕЛЬ" дало позитивні результати у зниженні на 20°C температури на корпусі печі, скороченні на 5-15% питомої витрати електроенергії на 1 тонну катанки та поліпшенні її якості шляхом зменшення відсотку сторонніх домішок і забезпечення вмісту в ній міді 99,99%.<br>The thesis for a Candidate of Engineering Sciences degree by specialty 05.11.13 – devices and methods of control and determination of composition of substances. – National technical university "Kharkov polytechnic institute", Kharkiv, 2019. The thesis solves the important scientific task to improve methods of control of the technical condition of induction installations for the production of high-quality copper rod using the mathematical modeling of electrothermal processes in copper melt and multi-layer thermal insulation, with taking into account its degradation and temperature distribution on the installation surface, which is essential for the development of cable industry. A new method of non-destructive control of the technical condition of induction installations for the production of copper rod for power cables of energy value is developed, in which in addition to measuring the reactive and active resistances of the inductor and the reduction of the water temperature in the cooling system, the degree of degradation of multilayer thermal insulation is additionally determined by comparing the calculated distribution of the current temperature in the volume of melt and insulation on a three-dimensional mathematical model of the installation with practical measurement of temperature in local areas (in the work they were 72) of the surface of the furnace casing. This makes it possible to more accurately predict the furnace's profile and to detect the modes of its operation in a timely manner close to the emergency. A new technical solution for improving the multilayer thermoinsulation of induction installations has been developed in order to reduce their energy consumption and increase the resource. It consists in applying instead of the third thermal insulation layer of light-weight brick of three layers: monolithic refractory concrete, lightweight brick and refractory paper, with the preservation of the overall thickness of the insulation. The implementa-tion of the developed lining structure improvement at PJSC " YUZHCABLE WORKS" in the UPCAST US20X-10 induction installation for copper rod in the cable industry in Ukraine has yielded positive results in a 20°C decrease in the temperature of the furnace casing, a reduction of 5-15% of the specific cost Electricity for 1 ton of production and improvement of quality of copper rod, by reducing the percentage of foreign impurities and ensuring copper content in it 99,99%.

Ce mémoire présente l’étude du comportement d’un flotteur bord à quai. Ce flotteur est constituéd’une bouée de section rectangulaire située à proximité d’une digue verticale. Un modèle analytique linéaire basé sur la théorie de l’écoulement potentiel est réalisé en profondeur d’eau intermédiaire. Un modèle numérique résolvant les équations de Navier-Stockes est développé. Des essais expérimentaux sont également mis au point pour des flotteurs de différentes dimensions et servent de références aux deux autres modèles. La comparaison de ces trois modèles permet de déduire le domaine de validité du modèle analytique, ainsi que la présence d’effets non-linéaires agissant sur le comportement du flotteur.Dans un premier temps, l’impact des dimensions du système sur le comportement du flotteur est quantifié. Les effets non-linéaires sont quantifiés en fonction des grandeurs du système puis une correction du modèle analytique est proposée. Par la suite, la correction du modèle est appliquée au comportement du système incluant un récupérateur d’énergie.Dans un second temps, l’effet du flotteur sur les franchissements de digue est étudié. Les franchissements de digue sont exprimés grâce à la hauteur de houle entre la digue et le flotteur (exprimée par le modèle analytique) et à la formule de Van Der Meer (2002). Ces franchissements sont également estimés de manière expérimentale. Il est montré que lorsque la longueur d’onde de la houle est importante et que le dégagement entre la digue et le flotteur est faible, les franchissements de digue sont réduits par la présence du flotteur<br>The study of a quayside wave energy converter is proposed. A simplified model composed ofa rectangular buoy oscillating in heave motion and installed close to a vertical dike is developedusing three different models. An analytical model based on potential flow theory is developed inintermediate water depth. A numerical model resolving the Navier-Stokes equation is developed.Experimental tests are performed for different floater dimensions and are used as references for the two other models. The comparison between the three models presents the validity domain of the analytical model. The existence of non-linear effects acting on the floater behaviour is shown.At first, the impact of the system dimensions on the floater behaviour are highlighted and quantified. According to these results a correction of the analytical model is proposed. Moreover, the corrected analytical model is applied to the wave energy converter behaviour. Secondly, the overtoppings are expressed using the wave height between the dike and the floater (provided by analytical model) and the Van Der Meer formula (2002). These overtoppings are also estimated using experimental tests. The comparison of these two models shows the presence of strong load losses when the floater is close to the dike. In those cases the analytical model needs a correction taking into account these load losses. It is also shown that when at the same time the wave length is high and the toe clearance between the dike and the floater is small, the dike overtopping are reduced by the presence of the floater

En raison de l'épuisement des combustibles fossiles et en raison de son coût croissant, les déchets ont été utilisés comme combustibles alternatifs dans les fours rotatifs à ciment pendant plusieurs années. Afin de répondre aux exigences de protection de l'environnement et de qualité du produit final, il est nécessaire de comprendre et de quantifier les différents processus intervenant dans le four. Dans le cadre de cette thèse, nous nous sommes intéressés à la substitution partielle des déchets (pneus usés et Refuse Derived Fuel « RDF ») et biomasse (résidus d’agriculture) dans les fours tournants de production de ciment à un taux de 50 % d’apport énergétique. Le présent projet de thèse a été financé par la République Démocratique Congo, porté par l’Université Kongo (UK), dans le cadre du renfoncement des capacités des universités congolaises. Les travaux de recherche ont été réalisés au sein du Laboratoire de Thermique, Energétique et Procédés (LaTEP) de l’Ecole Nationale Supérieure en Génie de Technologies Industrielles (ENSGTI)/Université de Pau et des Pays de l’Adour (France), dans le but d’étudier et de comprendre les phénomènes thermochimiques se déroulant dans les fours tournants de production de ciments, dans le cas de substitution partielle des combustibles fossiles par des combustibles alternatifs. Pour cela, un modèle mathématique à deux sous-modèles a été établi : - Le modèle monodimensionnel stationnaire du lit de solides et de l’enveloppe du four, caractérisé dans sa première partie par le mélange de la farine de ciment, de déchets et/ou biomasse, et de gaz, et alors que dans sa deuxième partie par les échanges thermiques. Le travail de modélisation est effectué à l'aide du langage de programmation Fortran 90 dans lequel sont implémentées les équations de conservations de la masse, des espèces et de l’énergie, l’équation de transport de la charge, ainsi que les équations algébriques. - Le modèle tridimensionnel instationnaire, repose sur l’utilisation de la mécanique des fluides numérique, et l’outil commercial Fluent a été retenu pour représenter les processus se déroulant au-dessus de la surface libre du lit.Le modèle complet est le couplage entre les deux sous modèles, étant donné les échanges qui s’y effectuent.Deux grands groupes de simulations ont été réalisés dans cette thèse : le cas monodimensionnel et le cas global (couplage Modèle 1D-Modèle 3D). Dans les deux situations, les simulations ont été effectuées dans le cas de fonctionnement au fuel lourd seul, puis dans celui du fonctionnement avec substitution partielle du combustible fossile par les combustibles alternatifs. Pour le modèle monodimensionnel, les simulations ont été effectuées en imposant les flux de chaleur dans la surface libre du lit et dans la paroi découverte. Alors que pour le modèle global, les simulations ont été effectuées en considérant les densités de flux de chaleur associées aux phénomènes physico-chimiques (combustion du fuel lourd, rayonnement, turbulence) se déroulant dans le volume de gaz situé au-dessus de la surface libre du lit.Les résultats ainsi obtenus nous ont permis de conclure sur l’importance de mener en amont une étude expérimentale dans un four tournant à l’échelle pilote afin de pouvoir disposer des ordres de grandeur des différents flux de chaleur au sein du four tournant. Les résultats obtenus nous ont également permis de nous rendre compte de la nécessité de disposer de la géométrie exacte (emplacement) du brûleur dans le four tournant<br>Because of the depletion of fossil fuels and because of its increasing cost, waste has been used as alternative fuels in cement rotary kilns for several years. In order to fulfill the requirements of environmental protection and quality of the final product, it is necessary to understand and quantify the different processes occurring in the kiln.As part of this thesis, we focused on the partial substitution of waste (used tires and Refuse Derived Fuel "RDF") and biomass (agricultural residues) in rotary kilns producing cement at a rate of 50 % of energy supply. This thesis project was funded by the Democratic Republic of Congo, led by Kongo University (UK), as part of the capacity building of Congolese universities. The research work was carried out within the Laboratory of Thermal, Energy and Processes (LaTEP) of the National School of Engineering in Industrial Technologies (ENSGTI) /University of Pau, with the aim of to study and understand the thermochemical phenomena taking place in rotary kilns for the production of cement, in the case of partial substitution of fossil fuels by alternative fuels. For this, a mathematical model with two sub-models has been established:- The stationary one-dimensional model of the bed of solids and furnace casing, characterized in its first part by the mixing of cement flour, waste and / or biomass, and gas, and while in its second part by thermal exchanges. The modeling work is performed using the Fortran 90 programming language in which the conservation equations for mass, species and energy, the load transport equation, and the algebraic equations are implemented.- The unsteady three-dimensional model relies on the use of digital fluid mechanics, and the Fluent commercial tool has been retained to represent the processes taking place above the free surface of the bed.The complete model is the coupling between the two sub models, because of the exchanges that take place there.Two large groups of simulations have been realized in this thesis: the one-dimensional case and the global case (Model 1D-3D model coupling). In both situations, the simulations were carried out in the case of operation with heavy fuel oil alone, then in the case of operation with partial substitution of the fossil fuel by alternative fuels. For the one-dimensional model, the simulations were performed by imposing the heat flux in the free surface of the bed and in the exposed wall. Whereas for the global model, the simulations were performed considering the heat flux densities associated with physico-chemical phenomena (combustion of heavy fuel oil, radiation, turbulence) taking place in the volume of gas located above the free surface of the bed.The results thus obtained allowed us to conclude on the importance of carrying out an experimental study in a pilot scale furnace in order to be able to have orders of magnitude of the different heat flows within the rotating furnace. The results obtained also allowed us to realize the necessity to model the exact geometry (location) of the burner in the rotary kiln

La diversité des régimes des instruments de la famille des flûtes a été mise en évidence à de nombreuses reprises : régimes statiques, périodiques, ou non périodiques. Cependant, de nombreux aspects de la dynamique de ces instruments demeurent mal compris. Pour les musiciens comme pour les facteurs d'instruments, les transitions entre régimes revêtent une importance particulière : d'une part elles correspondent à des changements de notes, et d'autre part la production d'un régime donné est conditionnée par les paramètres de facture (liés à la fabrication de l'instrument), et de contrôle (ajustés en permanence par l'instrumentiste). On s'attache dans ce document à caractériser les transitions entre régimes dans les flûtes, en lien avec des problématiques de facture et de jeu. Différentes approches sont mises en place. Des approches expérimentales d'une part, avec des mesures sur musicien et sur bouche artificielle. Par ailleurs, un modèle physique de l'instrument - un système dynamique à retard de type neutre - est étudié, par intégration temporelle d'une part, mais également par collocation orthogonale et continuation, donnant ainsi accès aux diagrammes de bifurcations.Croiser les résultats de ces différentes approches permet de mieux appréhender différents phénomènes : hystérésis associée aux changements de régime, ou mécanisme d'apparition des régimes non périodiques. L'influence de paramètres de facture et de contrôle est également étudiée : le rôle majeur de la géométrie interne du canal des flûtes à bec est mis en évidence, et l'influence de la dynamique de la pression dans la bouche du musicien sur les seuils de changement de régimes est caractérisée<br>Various studies have highlighted the diversity of regimes in flute-like instruments : static, periodic or non periodic regimes. However, some aspects of their dynamics remain poorly understood. Both for flute players and makers, transitions between regimes are particularly important : on the one hand, they correspond to a change of the note played, and on the other hand, production of a given regime is determined by parameters related to making and to playing of the instrument. In this document, we are interested in caracteristics of regime change in flute-like instruments, in relation with making and playing issues.Different approches are considered. First, experimental methods, with measurement on both musician and an artificial mouth. On the other hand, a physical model of the instrument - a system of delay differential equations of neutral type - is studied, through time-domain integration, and using orthogonal collocation coupled to numerical continuation. This last approach provides access to bifurcation diagrams.Considering results of these different methods, it becomes possible to better understand different experimental phenomena, such as regime change and associated hysteresis, or production mechanisms of non periodic regimes. Influence of different parameters is further studied : the crucial importance of the channel geometry in recorders is highlighted, and the influence of the mouth pressure dynamics on regime change thresholds is analysed

Large-scale bed features are often encountered in coastal waters, and include sandbanks and spoil heaps. The morphodynamic development of such features involves complicated nonlinear interactions between the flow hydrodynamics, sediment transport, and bed profile. Numerical modelling of the morphodynamic evolution and migration of large-scale bed features is necessary in order to understand their long-term behaviour in response to changing environmental conditions. This thesis describes detailed measurements of the morphodynamics of sand mounds in unidirectional and oscillatory (tidal) flows, undertaken at the U.K. Coastal Research Facility (UKCRF). High quality data were collected, including water velocities, water levels and overhead images. The parameters tested are: three types of mound shape (circular and elliptical in plan shape, and Gaussian, cosine and triangular in cross-section); underlying fixed or mobile bed conditions; and initial crest height (submerged, surface-touching and surface-piercing). Peak flow velocities are about 0.5 m/s, the sand median grain size is 0.454 mm, and transport occurring mostly as bedload. When analysing the data, the bed contours are determined by digitising the shoreline at different water levels. From these plots, the volume, height, and centroid position of the mound are calculated. A large-scale fit method, based on a Gaussian function has been used to separate small-scale ripples from the large-scale bed structure during the evolution of an isolated sand mound or spoil heap. The bed profile after the ripples are removed is comparable to typical predictions by shallow-flow numerical solvers. The UKCRF experiments investigated the morphodynamic response of a bed mound to hydrodynamic forcing: shape changes, migration rates, volume decay and sediment transport rates. The measured migration rate and decay of a submerged sand mound in the UKCRF are found to be in satisfactory agreement with results from various theoretical models, such as the analytical solution derived by De Vriend. Numerical predictions of mound evolution by a commercial code, PISCES, are also presented for a fully submerged sand mound; the bed evolution is reasonably similar to that observed in the UKCRF. The data provided as a result of the research reported in this thesis provide insight into the behaviour of sand mounds in steady and unsteady flows at laboratory scale, and should also be useful for benchmark (validation) purposes to numerical modellers of large-scale morphodynamics.

This thesis details the development, verification and validation of an unsteady unstructured high order (≥ 3) h/p Discontinuous Galerkin - Fourier solver for the incompressible Navier-Stokes equations on static and rotating meshes in two and three dimensions. This general purpose solver is used to provide insight into cross-flow (wind or tidal) turbine physical phenomena. Simulation of this type of turbine for renewable energy generation needs to account for the rotational motion of the blades with respect to the fixed environment. This rotational motion implies azimuthal changes in blade aero/hydro-dynamics that result in complex flow phenomena such as stalled flows, vortex shedding and blade-vortex interactions. Simulation of these flow features necessitates the use of a high order code exhibiting low numerical errors. This thesis presents the development of such a high order solver, which has been conceived and implemented from scratch by the author during his doctoral work. To account for the relative mesh motion, the incompressible Navier-Stokes equations are written in arbitrary Lagrangian-Eulerian form and a non-conformal Discontinuous Galerkin (DG) formulation (i.e. Symmetric Interior Penalty Galerkin) is used for spatial discretisation. The DG method, together with a novel sliding mesh technique, allows direct linking of rotating and static meshes through the numerical fluxes. This technique shows spectral accuracy and no degradation of temporal convergence rates if rotational motion is applied to a region of the mesh. In addition, analytical mappings are introduced to account for curved external boundaries representing circular shapes and NACA foils. To simulate 3D flows, the 2D DG solver is parallelised and extended using Fourier series. This extension allows for laminar and turbulent regimes to be simulated through Direct Numerical Simulation and Large Eddy Simulation (LES) type approaches. Two LES methodologies are proposed. Various 2D and 3D cases are presented for laminar and turbulent regimes. Among others, solutions for: Stokes flows, the Taylor vortex problem, flows around square and circular cylinders, flows around static and rotating NACA foils and flows through rotating cross-flow turbines, are presented.

The diploma thesis deals with experimental determination of hydrodynamic load on the overflowed bridge deck model. In the first part of the thesis the author describes the analysis of the problem together with the basic physical laws and principles that are used or assumed in the measurement itself. In the second part the author describes the measurements in the laboratory of the Faculty of Civil Engineering and its gradual processing. At the end of the work are described the results of measurement, their comparison with numerical modeling, other authors and their possible use in practice.

博士<br>國立臺灣大學<br>機械工程學研究所<br>96<br>In the present study, a three-dimensional simulation is performed for the turbulent reactive flow and radiactive heat transfer in the walking-beam type slab reheating furnace by STAR-CD software. The study employs the high-Reynolds-number k-ε turbulence model based on Favre-averaged governing equations. The pre-assumed PDF model associated with the fast chemistry assumption and a single diffusivity is used to account for turbulent combustion. The absorption coefficient of the gases mixture is calculated by WSGGM (weighted-sum-of-gray-gases model). The discrete ordinates radiation model is adopted to calculate the radiactive heat transfer. The geometric model takes care of all components of the furnace, including the burners, the walking beam system with skid buttons, the slab, the dam, the down-take, etc. The part 1 of the study, the surface temperature of the slab is prescribed via experimental measurements and the furnace wall is assumed adiabatic. The turbulent reactive flow is thus simulated. The part 2 of the study, the temperature distributions of the slab and the gas mixture are obtained through a coupled calculation. The slab is modeled as a laminar flow having a very high viscosity and thus moving at a nearly constant speed. No radiation is concerned within the slab. The temperature distributions of the slab and the gas mixture are obtained through a coupled calculation. And the part 3 of the study, a solution for improving the skid mark by varying either the distance between two static beams or the height of the skid buttons is targeted. To obtain a steady solution, the walking beams are assumed fixed in the furnace. The simulation results agree with the measurements very well. The difference between the predicted heating efficiency and the measured one of the furnace is only 6%. The prediction errors at six temperature-monitored points are all under 10%, except the one in the lower heating zone which appears to be 13%. The measured surface temperatures of the slab were used as boundary conditions and the flow field of the reheating furnace can be obtained. However, the measured temperature will affect the result of calculations. When coupled calculations are executed in the part 2 of the study, the surface temperature can be obtained simultaneously. Except for the upper surface temperature, all the other simulation results agree with the measurements very well. Most of all, the influence of the walking beam system on the skid marks is thoroughly explored. The simulation results show that the radiative shielding by the static beams is the main cause of the skid mark. The heat loss through the skid button to the cooling system worsens the skid mark. A parametric study then shows by shifting the static beams inward in the soaking zone, the skid mark gets even worse due to an enhanced radiative shielding on the lower surface of the slab in between the two static beams; the skid mark on the opposite can be improved by shifting the static beams outward in the soaking zone. Another parametric study shows the skid mark can also be improved by increasing the height of the skid button.

A thesis submitted in ful lment of the requirements for the degree of Master of Science in the School of Computer Science and Applied Mathematics, 2018<br>A linear and a non-linear partial di erential equation is used as a phenomenological model to describe the synthesis of a vibrating string. The models are numerically simulated using nite di erence method and an exact solution of the linear partial di erential equation is derived using a functional transformation method. The functional transformation method is based on the transformation to the frequency domain in time and in space by applying the Laplace transformation and the Finite Sine transformation. The string is excited through the initial conditions, initial displacement and an initial velocity, of the system. Another excitation mode that is considered is the force density, where the initial conditions are set to zero. The force density is added as a Source term to the partial di erential equation and the initial excitation enters the system by the use of this term. To analyse the model di erent observation positions are considered, the observation points are the positions of the string where the vibrations are observed. The exact solution is used to verify the Finite di erence approximation and a spectrogram analysis of the di erent excitation points is used to asses the methods. The nite di erence method is then used to analysis the vibrations of the string. The nonlinear model include the e ects of tension modulation, which is de ned by the Kirchho -Carrier equation. The results are analysed using a quantitative analysis and a Fourier analysis.<br>XL2019

碩士<br>國立臺灣大學<br>地質學研究所<br>87<br>Geological modeling is a new way in geological studies. By the help of physical modeling and numerical modeling, we can understand the possible models of the development of the kinematics and dynamics evolution of structure. These studies try to estimate the possibility of sandbox simulation experiments in physical modeling; and try to analyze the effect of boundary conditions and main controlling factors in physical modeling. Furthermore, by using physical modeling and three-dimensional numerical modeling, we try to simulate the consequence of geological settings, oblique convergence and basement high in the inter-relationship between the development, evolution and structure in the western Taiwan. We hope that will help us understand more about the regional tectonics of Taiwan. By the simulation of sandbox behavior, there are some conclusions: most models develop conjugate thrust faults and pup-up structures; in most cases, the same set of conjugate thrust faults developed earlier than the backthrust faults. The modeling result is highly affected by the friction coefficient of backstop; as the friction coefficient reduced, the strain would concentrate on the top part of the backstop and increase along the ramp of the backstop. The greater the angle of inclination of the backstop, the easier the forward propagation of the strain. Reduction in the angle of inclination results in the concentration of strain near the toe front of the backstop, and a set of pup-up structures could be found at the toe of the backstop. As the friction coefficient of the base reduced, the forward propagation of the strain becomes more easily. Thought the competence of materials would affect the displacement of the propagation, it is not as important as the others factors. This research aimed to realize the effect of the basement high of western Taiwan during oblique convergence. Results: Backstop mainly affects the earlier stage of the oblique convergence; the shape of the basement and the early stage structures formed by the indent backstop controls the shape and development of later stage structures. As a consequence, forward propagation of the deformation front occurred in the central western Taiwan, because of the resistance at the front of the basement high, the deformation front propagated slowly at the front of the basement high. There is low resistance in the southern Taiwan and the sediments there are thickened, the deformation front moved forward. The above reasons may explain the S-shaped deformation front at the western Taiwan. Most thrust faults are forward-breaking sequence and it develops later in the western region. The formation of Tatu terrace, Bakua terrace, west Hengchun terrace may be the result of thrust faults; Pingtung valley, Taichung basin, Hengchun plain may be piggy-back basins. And these results are supported by the preliminary results of these experiments. Sandbox physical modeling is very useful in the research of the development and evolution of geological structures, and it showed great potential in this research. From the analysis of the parameters of the three-dimensional numerical modeling, we suggest that the interface friction angle of the basement high should be smaller than 15°. The existence of the basement high in the numerical modeling results in the change in the direction of stress field and strain field. Based on the nature of faults, we suggest that the faults near basement ramp tend to develop thrust faults and strike-slip faults would develop at later stage. Since the strain energy concentrated at the front of the West-Foothill near Chiayi and Tainan, these areas may affect by earthquakes.

<div> <p>The electric arc furnace (EAF) is one of the most widely-used steelmaking process recycling the scrap as the raw material. Electric arc furnace steelmaking is a material-dependent and energy-consuming process. The electricity is utilized as the main heat source to melt the scrap by the electric arc discharging from the graphite electrodes. Since the plasma arc weld has the very similar heat transfer mechanism and melting phenomenon comparing with the scrap melting in EAF process, the model development for the free-burning arc, the arc-metal heat transfer, and the metal melting will refer to the research of plasma arc weld field and start with a relatively small scale. In this study, a comprehensive computational dynamics (CFD) model was developed to predict the heat transfer from the electric arc to the metal anode, which is composed of the DC electric arc model, solidification and melting model, and the arc-metal heat transfer model. The validation of the CFD models has been conducted utilizing the experimental data and simulation results in other literature. The commercial software, ANSYS FLUENT^®, was employed with the User-Defined Function scripts and the User-Defined Scalar to model the magnetic field, comprehensive flow field, and high temperature field. Furthermore, the parametric studies for the process of the anode melted by the plasma arc were performed to investigate the effects of the arc current and the initial anode temperature on the anode melting. The results reveal that the value of the arc current has a positive correlation with the arc temperature and velocity but has negative correlation with the penetration time of the anode. Meanwhile, with the increase of anode initial temperature, the metal penetrate time will decrease.</p> </div> <br>