Volcanic eruptions pose hazards to life and insfrastructure, and contribute to the resurfacing of earth and other planetary bodies. Lavas and magmas are multi-phase suspensions of silicate melts (liquids), solid crystals, and vapor bubbles, and solidify into glass and rock upon cooling. The interactions between phases place important controls on the dynamics and timescales of magma and lava transport and emplacement. The purpose of this thesis is to explore the role of multiphase interactions in controlling eruption dynamics and inform conceptual and numerical models for hazard prediction. In Chapters 1 and 2, centimeter to meter scale analog experiments are used to explore the multi-phase rheological properties and flow behaviors of bubble- and particle-bearing suspensions. Optical imaging of dam-break experiments presented in Chapter 1 expand existing experimental parameter ranges for lava analogs to higher bubble concentrations than existing datasets (up to 82% by volume bubbles and 37% by volume particles). I develop a constitutive relationship for threephase relative viscosity, and demonstrate that at the low strain-rate conditions relevant to many natural lava flows, accounting for the rheological effect of bubbles can result in the prediction of slower runout speeds. Chapter 2 expands upon the work of Chapter 1 using different analog materials observed using nuclear magnetic resonance imaging (MRI) phase-contrast velocimetry (PCV) to measure velocity in the flow interior of three-phase dam-break experiments. I find that for high-aspect ratio particles (sesame seeds), phase segregation into shear bands readily occurs, even at low particle fraction (20%) and results in strain localization. I suggest that the presence of shear bands can lead to faster flow runout than predicted using assumptions of bulk rheology. Chapter 3 analyzes thermal infrared (IR) time-lapse photography and videography of Hawaiian to Strombolian explosive activity during the 2021 eruption of Cumbre Vieja volcano, La Palma, Canary Islands, Spain. Images are analyzed to find time series of apparent plume radius, velocity, and apparent volume flux of high-temperature gas and lava. I compare with other measures of eruptive activity, including remote observations of plume height, SO₂ flux, effusive flux, tremor, and events at the volcano edifice including edifice collapses and the opening of new vents. I find correlations between tremor and explosive flux, but no correlation with SO2 flux or effusive flux, which I interpret as evidence of bubble segregation, highlighting the role of phase segregation and temporal variability in material properties in natural systems. Finally, in Chapter 4, I develop a novel finite element model to explore the interaction between a viscous flow with a solidified crust, and the effect of these interactions on lava flow and lava dome emplacement. I develop a model that couples a temperature-dependent viscous interior with an elastic shell flowing into air, water, or dense atmospheres. The model expands upon existing numerical simulations used in volcanology to have direct applications to lava flows and domes on the sea floor, which accounts for a large portion of the volcanism on Earth, and volcanism on other planetary bodies. Additionally, the formation of levees or solidified flow fronts that fracture and lead to a restart of flow. These lava flow breakouts pose a significant hazard, but there are currently no volcanological community codes capable of using a physics-based approach to predict the timing or location of breakouts. The model in Chapter 4 is the first to allow for assessment of the likelihood of failure at the scale of a flow lobe. Chapter 4 describes the model formulation and verification, and validation against centimeter-scale molten basalt experiments. The dissertation as a whole integrates work using a variety of methods including analog experiments, observations of natural eruptions, and numerical simulations to contribute to our understanding of the effects of multi-phase interactions on volcanic eruptions.

High-amplitude freestream turbulence and surface roughness elements can excite a laminar boundary-layer flow sufficiently enough to cause streamwise-oriented vortices to develop. These vortices resemble elongated streaks having alternate spanwise variations of the streamwise velocity. Following the transient growth phase, the fully developed vortex structures downstream undergo an inviscid secondary instability mechanism and, ultimately, transition to turbulence. This mechanism becomes much more complicated in high-speed boundary layer flows due to compressibility and thermal effects, which become more significant for higher Mach numbers. In this research, we formulate and test an optimal control algorithm to suppress the growth rate of the aforementioned streamwise vortex system. The derivation of the optimal control algorithm follows two stages. In the first stage, to optimize the computational cost of the analysis, the study develops an efficient numerical algorithm based on the nonlinear boundary region equations (NBREs), a reduced form of the compressible Navier-Stokes equations in a high-Reynolds-number asymptotic framework. The NBREs algorithm results agree well with direct numerical simulation (DNS) results. The numerical simulations are substantially less computationally costly than a full DNS and have a more rigorous theoretical foundation than parabolized stability equation (PSE) based models. The substantial reduction in computational time required to predict the full development of a G\"{o}rtler vortex system in high-speed flows allows investigation into feedback control in reasonable total computational time, which is the focus of the second part of the study. In the second stage, the method of Lagrange multipliers is utilized -- via an appropriate transformation of the original constrained optimization problem into an unconstrained form -- to obtain the adjoint compressible boundary-region equations (ACBREs) and corresponding optimality conditions, which constitute the basis of the optimal control approach. Numerical solutions for high-supersonic and hypersonic flows reveal a significant decrease in the kinetic energy and wall shear stress for all configurations considered. Streamwise velocity contour plots illustrate the qualitative effect of the optimal control iterations, demonstrating a significant decrease in the amplitude of the primary vortex instabilities.

<p> </p> <p>Interfaces separating two fluids are incredibly complex physical structures and are common throughout science, technology, and nature. Examples from daily life include the air-water interface separating a water drop that is dripping from a leaky faucet from the surrounding air and the interface of a soap bubble (which actually consists of two interfaces) separating the interior of the bubble from its exterior. Other common examples from nature include interfaces between falling rain drops and the surrounding air, and the mist that one encounters at beaches, waterfalls, and fountains where the spray droplets are separated from the surrounding air by an interface. Interfaces and manipulating them are key to technological applications such as thin film coating flows and diverse processes involving drop-by-drop processing such as ink-jet printing, drop-wise manufacturing, spray coating, DNA microarraying, and chemical separations, e.g. extraction. Aside from the coating flows example, the aforementioned situations are all examples of free surface flows that involve abrupt and catastrophic topological changes of interfaces that include physical processes such as breakup (also called pinch-off) as in drop breakup, rupture as in liquid-film or liquid-sheet rupture, and coalescence as in drop or bubble coalescence (similar phenomena also arise in sintering and/or fusion of ceramic, metallic, and polymer particles). These topological changes entail what are referred to as finite-time hydrodynamic singularities. For example, at the location(s) where a drop breaks, the thickness of the drop locally tends to zero while fluid pressure and velocity diverge (hence the reason for the word singularity). In addition to hydrodynamic singularities, the presence of surface-active agents or surfactants at fluid interfaces in free surface flows is another reason scientists have been attracted to the study of such problems.</p> <p>Adsorption onto and lowering of the surface tension of a fluid interface by surfactants are exploited in applications such as enhanced oil recovery, coating flows, lung surfactants, drop/jet breakup, and film/sheet rupture, with the latter two being among the prime motivators for this PhD thesis. However, surfactant concentration can be nonuniform at the interface because surfactant molecules can be transported along it by convection and diffusion and also due to normal dilatation and tangential stretching of the interface. Thus, aside from simply lowering surface tension, nonuniformity in surfactant concentration causes gradients in surface tension and gives rise to tangential interfacial (Marangoni) stresses. The latter brings about rich physics including tears of wine, interfacial turbulence in mass transfer, and droplet bouncing. In addition to lowering surface tension and the Marangoni effect, surfactants may also induce surface rheological or viscous effects as surfactant molecules deform against each other. The primary goal of this thesis is to advance the understanding of surface rheological effects in situations involving the breakup of surfactant-covered liquid threads (which also includes jets and drops) and liquid sheets. The fundamental understanding developed in this thesis is likely to prove indispensable in and/or assist the development of new technologies where surface rheological effects are central to the processes at hand, e.g. in controlling drop size distributions and avoiding undesirable satellite droplets and/or misting. An initially unexpected but highly rewarding outcome of the research has been the development of techniques for the measurement of surface viscosities, a task that has heretofore proven to be a formidable challenge to experimentalists.</p> <p>In this thesis, surface rheological effects in free surface flows are examined through both analytical and numerical solution of the incompressible Navier-Stokes equations subjected to the traction boundary condition augmented by the Boussinesq-Scriven constitutive equation to account for surface viscous effects. Rigorous and robust numerical algorithms based on the Galerkin finite element (GFEM) method are developed for predictions of surfactant transport, surface rheological effects and hydrodynamics in response to the motion of moving boundaries. The accuracy of computational predictions is verified by demonstrating that computed results accord well with scaling theories.</p>

In this report a discriminative and a generative image classifier, used for classification of images with handwritten digits from zero to nine, are compared. The aim of this project was to compare the accuracy of the two classifiers in absence and presence of perturbations to the images. This report describes the architectures and training of the classifiers using PyTorch. Images were perturbed in four ways for the comparison. The first perturbation was a model-specific attack that perturbed images to maximize likelihood of misclassification. The other three image perturbations changed pixels in a stochastic fashion. Furthermore, The influence of training using perturbed images on the robustness of the classifier, against image perturbations, was studied. The conclusions drawn in this report was that the accuracy of the two classifiers on unperturbed images was similar and the generative classifier was more robust against the model-specific attack. Also, the discriminative classifier was more robust against the stochastic noise and was significantly more robust against image perturbations when trained on perturbed images. / I den här rapporten jämförs en diskriminativ och en generativ bildklassificerare, som används för klassificering av bilder med handskrivna siffror från noll till nio. Syftet med detta projekt var att jämföra träffsäkerheten hos de två klassificerarna med och utan störningar i bilderna. Denna rapport beskriver arkitekturerna och träningen av klassificerarna med hjälp av PyTorch. Bilder förvrängdes på fyra sätt för jämförelsen. Den första bildförvrängningen var en modellspecifik attack som förvrängde bilder för att maximera sannolikheten för felklassificering. De andra tre bildförvrängningarna ändrade pixlar på ett stokastiskt sätt. Dessutom studerades inverkan av träning med störda bilder på klassificerarens robusthet mot bildstörningar. Slutsatserna som drogs i denna rapport är att träffsäkerheten hos de två klassificerarna på oförvrängda bilder var likartad och att den generativa klassificeraren var mer robust mot den modellspecifika attacken. Dessutom var den diskriminativa klassificeraren mer robust mot slumpmässiga bildförvrängningar och var betydligt mer robust mot bildstörningar när den tränades på förvrängda bilder. / Kandidatexjobb i elektroteknik 2022, KTH, Stockholm

This thesis focuses on evaluating the appropriate speed at which vehicles should travel under different traffic conditions on freeways and its impact on crash frequency. The common belief is that the lower speed results in fewer crashes as reduced speed provides drivers with more time to react effectively and avoid collisions. However, this perspective overlooks the interplay among traffic speed, average spacing between consecutive vehicles, and the distance available for stopping a vehicle. Hence, we propose a safety parameter termed ‘Safety Correlate' (SCORE), which is defined as the proportion of average spacing relative to the stopping distance. To determine the relationship between SCORE and crash frequency, data from 366 4-lane urban freeway segments located in Virginia was analyzed and a Random-effects Poisson Lognormal model was developed. The obtained result indicated that the safety parameter SCORE is negatively associated with the annual hourly crash frequency, implying that the lesser the average spacing as a proportion of the stopping distance while traffic flow remains constant, the more frequent will be the crashes. Additionally, this research presents an application of SCORE in setting variable speed limits under various traffic flows. Overall, the study results provide valuable insights by investigating SCORE to improve traffic safety. Also, this research would help practitioners and policymakers to incorporate safety aspects while setting variable speed limits on freeways.

Dams and power plants in watercourses pose a serious threat to diversity and aquatic species as they fragment habitats and create barriers. New legislation aims to facilitate the restoration of watercourses and in many places, various fish passages have been built. Evaluations of these are important in order to ensure that the desired function has been achieved. At Näkten regulation dam in Jämtland, Sweden, a combined technical and natural fishway has been built to enable migration for trout and grayling. The purpose of the study is to find out whether the technical solution has had suitable flows and gate settings for the fish that one wanted to enable migration for. To investigate this, unsteady flow analyzes were performedin HEC-RAS and then compared to field measurements as well as information from the dam owner. The results show that the gate during the study period has not been set according to recommended settings and that the flows at times have been too high to be suitable for migration. However, comparison with field measurements and previous calculations suggest that the modeled flows do not match reality but are too high in the model. The conclusion is therefore that further modeling and a more thorough comparison with field measurements is recommended to create reliable results and that one then possibly can adjust the settings on the gates. / Dammar och kraftverk i vattendrag utgör ett allvarligt hot mot diversiteten och vattenlevande arter då den fragmenterar livsmiljöer och skapar hinder. Nya lagstiftningar ska underlätta för återställning av vattendrag och på många platser har olika fiskpassager byggts. Utvärderingar av dessa är viktiga så att man kan säkerställa att önskad funktion är uppnådd. Vid Näktens regleringsdamm i Jämtlands län har ett kombinerat tekniskt oc hnaturlikt omlöp byggts för att åter möjliggöra vandring för öring och harr. Studiens syfte är att ta reda på om den tekniska lösningen har haft lämpliga flöden och luckinställningar för de fiskar som man velat möjliggöra vandring för. För att undersöka detta har flödesmodelleringar gjorts i programmet HEC-RAS och fältmätningar samt uppgifter från dammägaren har använts som referens. Resultaten visar att luckan under studieperioden inte varit inställd efter rekommenderade inställningar och att flödena ha varit för höga för att vara lämplig för vandring. Dock tyder jämförelse med fältmätningar och tidigare beräkningar att flödena inte stämmer överens med verkligheten utan blir för höga i modellen. Slutsatsen blir därför att ytterligare modellering och en grundligare jämförelse med fältmätningar rekommenderas för att skapa tillförlitliga resultat och att man därefter eventuellt kan justera inställningarna på luckorna. / <p>2023-06-02</p>

A swirl injector for a rocket engine being developed by \emph{AESIR} (Association of EngineeringStudents in Rocketry) was simulated with different geometric parameters. The swirl injector is usedto atomize the ethanol used as fuel and to create a spray that mixes well with the oxidizer withinthe combustion chamber. Inlet slot angle (90, 75, 60 and 45 degrees), swirl chamber length (15, 20and 25 mm) and outlet orifice diameter (3, 6 and 9 mm) were examined.Previous studies in swirl injectors show that CFD can be used to analyze the flow in such aninjector, furthermore theoretical models exist that can predict some of the general characteristicsof the flow. Previous studies have also simulated transient behavior and flow features effectingbreakup of fuel flowing through a swirl injector.A steady state simulation using Volume of Fluid (VOF) multiphase modeling and $k$-$\omega$ \emph{SST}turbulence modeling was used to simulate the swirl injector intended for the rocket engine. It wasfound that a wider outlet orifice would give a wider cone angle of spray. This is desirable in thecurrent rocket engine design as it will promote greater mixing of fuel and oxidizer higher up in thecombustion chamber. No large variances was observed when different inlet slot angles was simulated. Ashorter swirl chamber length reduced the amount of losses in energy due to viscous forces. The flowafter the outlet orifice was not simulated so the effect of turbulence kinetic energy and energylosses outside of the swirl injector have not been analyzed, previous studies have indicated thatturbulent kinetic energy does have an effect on the breakup and atomization of the fuel.It was concluded that using a wider outlet orifice of 9 mm gave the best results out of the differentgeometric parameters analyzed and the swirl chamber length should be a short as possible.

This thesis is divided into two projects. The first project investigates the dissolution of the Hanford salt cakes, the chemical properties of the effluent and the physical properties such as viscosity of the effluent, the porosity and the permeability of the salt cake bed as the dissolution proceeds. The chemical results are compared to predictions using a thermodynamic model. Physical properties are important because they govern the rate at which the Hanford tanks can be emptied thus facilitating the remediation process. Two simulants were investigated for the dissolution process. The chemical analysis matched with the model predictions for both the simulants. A typical gibbsite layer formation was observed in the chemically complex simulant and experiments were performed to remediate the layer. The second project of this thesis studied the remediation of sodium phosphate dodecahydrate plug using water and sodium carbonate solutions at varying concentrations. A flow loop previously used to study the sodium phosphate dodecahydrate plugging mechanisms, was used to form a plug followed by the addition of water and sodium carbonate solutions. Results indicate that there was a drastic decrease in time to remediate the plug when sodium carbonate solutions were used.

Yes / Simulation of the suspended sediment concentration (SSC) has great significance in predicting the sediment transport rate, vegetation growth and the river ecosystem in the vegetated open channel flows. The present study focuses on investigating the vertical SSC profile in the vegetated open channel flows. To this end, a model of the dispersive flux is proposed in which the dispersive coefficient is expressed as partitioned linear profile above or below the half height of vegetation. The double-averaging method, i.e. time-spatial average, is applied to improve the prediction accuracy of the vertical SSC profile in the vegetated open channel flows. The analytical solution of SSC in both the submerged and the emergent vegetated open channel flows is obtained by solving the vertical double-averaging sediment advection-diffusion equation. The morphological coefficient, a key factor of the dispersive coefficient, is obtained by fitting the existing experimental data. The analytically predicted SSC agrees well with the experimental measurements, indicating that the proposed model can be used to accurately predict the SSC in the vegetated open channel flows. Results show that the dispersive term can be ignored in the region without vegetation, while the dispersive term has significant effect on the vertical SSC profile within the region of vegetation. The present study demonstrates that the dispersive coefficient is closely related to the vegetation density, the vegetation structure and the stem Reynolds number, but has little relation to the flow depth. With a few exceptions, the absolute value of the dispersive coefficient decreases with the increase of the vegetation density and increases with the increase of the stem Reynolds number in the submerged vegetated open channel flows. / the Natural Science Foundation of China (Nos. 11872285 and 11672213), The UK Royal Society – International Exchanges Program (IES\R2\181122) and the Open Funding of State Key Laboratory of Water Resources and Hydropower Engineering Science (WRHES), Wuhan University (Project No: 2018HLG01)

Settling and consolidation of suspended clay particles are significant issue in many fields such as geotechnical engineering, coastal and hydraulic engineering, and environmental engineering. A comprehensive literature review was conducted on the settling, consolidation and erosion of mixed soil material (cohesive and non-cohesive). Soil beds formed by sedimentation process of loose particles will be either show a segregated or homogeneous in structure, depending on the depositional environment. These sediments initially undergo self-weight consolidation and may be eroded under high flow rate. A number of studies have recently investigated the characteristic of consolidated clay bed in stagnant water. Hence, consolidation parameters were determined using a well-known vertical settling column consolidation test setup. However, limited research studies are available for deposition and consolidation of a mixture of sediment (clay, silt and sand) under flow conditions which are more representative of what happens in nature. A long flume and pump were used to create different turbulent conditions and simulate the natural process... The results for deposition and consolidation of different mixtures under stagnant and turbulence conditions were analyzed and compered in term of compressibility, permeability as well as shear strength. The results of this experimental research program indicated that the flow rate, initial concentration, height of settling and composition of sediment are all important factors that could affect the final bed dry unit weight. Two non-intrusive techniques were applied for measuring the dry unit weight at settling and consolidation stages. Impact echo technique has never been applied to measure the dry unit weight of self-weight consolidation along the vertical stratification of cohesive and non-cohesive particles. Also, a novel conductance sensor has been developed to improve the efficiency of this technique. The limitations of using these techniques will be highlighted in this study.