Experimental analysis of coaxial jets: instability, flow and mixing characterization

The velocity and mixing field of two turbulent jets configurations have been experimentally characterized by means of cold- and hot-wire anemometry in order to investigate the effects of the initial conditions on the flow development. In particular, experiments have been focused on the effect of the separation wall between the two streams on the flow field. The results of the experiments have pointed out that the wake behind a thick wall separating wall has a strong influence on the flow field evolution. For instance, for nearly unitary velocity ratios, a clear vortex shedding from the wall is observable. This phenomenon enhances the mixing between the inner and outer shear layer. This enhancement in the fluctuating activity is a consequence of a local absolute instability of the flow which, for a small range of velocity ratios, behaves as an hydrodynamic oscillator with no sensibility to external perturbations. It has been suggested indeed that this absolute instability can be used as a passive method to control the flow evolution. Finally, acoustic excitation has been applied to the near field in order to verify whether or not the observed vortex shedding behind the separating wall is due to a global oscillating mode as predicted by the theory. A new scaling relationship has been also proposed to determine the preferred frequency for nearly unitary velocity ratios. The proposed law takes into account both the Reynolds number and the velocity ratio dependence of this frequency and, therefore, improves all the previously proposed relationships.

Studio dei metodi di valutazione della pericolosità associata ai "debris flow" in ambiente alpino

The present work concerns with the study of debris flows and, in particular, with the related hazard in the Alpine Environment. During the last years several methodologies have been developed to evaluate hazard associated to such a complex phenomenon, whose velocity, impacting force and inappropriate temporal prediction are responsible of the related high hazard level. This research focuses its attention on the depositional phase of debris flows through the application of a numerical model (DFlowz), and on hazard evaluation related to watersheds morphometric, morphological and geological characterization. The main aims are to test the validity of DFlowz simulations and assess sources of errors in order to understand how the empirical uncertainties influence the predictions; on the other side the research concerns with the possibility of performing hazard analysis starting from the identification of susceptible debris flow catchments and definition of their activity level. 25 well documented debris flow events have been back analyzed with the model DFlowz (Berti and Simoni, 2007): derived form the implementation of the empirical relations between event volume and planimetric and cross section inundated areas, the code allows to delineate areas affected by an event by taking into account information about volume, preferential flow path and digital elevation model (DEM) of fan area. The analysis uses an objective methodology for evaluating the accuracy of the prediction and involve the calibration of the model based on factors describing the uncertainty associated to the semi empirical relationships. The general assumptions on which the model is based have been verified although the predictive capabilities are influenced by the uncertainties of the empirical scaling relationships, which have to be necessarily taken into account and depend mostly on errors concerning deposited volume estimation. In addition, in order to test prediction capabilities of physical-based models, some events have been simulated through the use of RAMMS (RApid Mass MovementS). The model, which has been developed by the Swiss Federal Institute for Forest, Snow and Landscape Research (WSL) in Birmensdorf and the Swiss Federal Institute for Snow and Avalanche Research (SLF) takes into account a one-phase approach based on Voellmy rheology (Voellmy, 1955; Salm et al., 1990). The input file combines the total volume of the debris flow located in a release area with a mean depth. The model predicts the affected area, the maximum depth and the flow velocity in each cell of the input DTM. Relatively to hazard analysis related to watersheds characterization, the database collected by the Alto Adige Province represents an opportunity to examine debris-flow sediment dynamics at the regional scale and analyze lithologic controls. With the aim of advancing current understandings about debris flow, this study focuses on 82 events in order to characterize the topographic conditions associated with their initiation , transportation and deposition, seasonal patterns of occurrence and examine the role played by bedrock geology on sediment transfer.

Fiber beam-columns models with flexure-shear interaction for nonlinear analysis of reinforced concrete structures

The aim of this study was to develop a model capable to capture the different contributions which characterize the nonlinear behaviour of reinforced concrete structures. In particular, especially for non slender structures, the contribution to the nonlinear deformation due to bending may be not sufficient to determine the structural response. Two different models characterized by a fibre beam-column element are here proposed. These models can reproduce the flexure-shear interaction in the nonlinear range, with the purpose to improve the analysis in shear-critical structures. The first element discussed is based on flexibility formulation which is associated with the Modified Compression Field Theory as material constitutive law. The other model described in this thesis is based on a three-field variational formulation which is associated with a 3D generalized plastic-damage model as constitutive relationship. The first model proposed in this thesis was developed trying to combine a fibre beamcolumn element based on the flexibility formulation with the MCFT theory as constitutive relationship. The flexibility formulation, in fact, seems to be particularly effective for analysis in the nonlinear field. Just the coupling between the fibre element to model the structure and the shear panel to model the individual fibres allows to describe the nonlinear response associated to flexure and shear, and especially their interaction in the nonlinear field. The model was implemented in an original matlab® computer code, for describing the response of generic structures. The simulations carried out allowed to verify the field of working of the model. Comparisons with available experimental results related to reinforced concrete shears wall were performed in order to validate the model. These results are characterized by the peculiarity of distinguishing the different contributions due to flexure and shear separately. The presented simulations were carried out, in particular, for monotonic loading. The model was tested also through numerical comparisons with other computer programs. Finally it was applied for performing a numerical study on the influence of the nonlinear shear response for non slender reinforced concrete (RC) members. Another approach to the problem has been studied during a period of research at the University of California Berkeley. The beam formulation follows the assumptions of the Timoshenko shear beam theory for the displacement field, and uses a three-field variational formulation in the derivation of the element response. A generalized plasticity model is implemented for structural steel and a 3D plastic-damage model is used for the simulation of concrete. The transverse normal stress is used to satisfy the transverse equilibrium equations of at each control section, this criterion is also used for the condensation of degrees of freedom from the 3D constitutive material to a beam element. In this thesis is presented the beam formulation and the constitutive relationships, different analysis and comparisons are still carrying out between the two model presented.

CFD modeling of two-phase boiling flows in the slug flow regime with an interface capturing technique

The objective of this thesis was to improve the commercial CFD software Ansys Fluent to obtain a tool able to perform accurate simulations of flow boiling in the slug flow regime. The achievement of a reliable numerical framework allows a better understanding of the bubble and flow dynamics induced by the evaporation and makes possible the prediction of the wall heat transfer trends. In order to save computational time, the flow is modeled with an axisymmetrical formulation. Vapor and liquid phases are treated as incompressible and in laminar flow. By means of a single fluid approach, the flow equations are written as for a single phase flow, but discontinuities at the interface and interfacial effects need to be accounted for and discretized properly. Ansys Fluent provides a Volume Of Fluid technique to advect the interface and to map the discontinuous fluid properties throughout the flow domain. The interfacial effects are dominant in the boiling slug flow and the accuracy of their estimation is fundamental for the reliability of the solver. Self-implemented functions, developed ad-hoc, are introduced within the numerical code to compute the surface tension force and the rates of mass and energy exchange at the interface related to the evaporation. Several validation benchmarks assess the better performances of the improved software. Various adiabatic configurations are simulated in order to test the capability of the numerical framework in modeling actual flows and the comparison with experimental results is very positive. The simulation of a single evaporating bubble underlines the dominant effect on the global heat transfer rate of the local transient heat convection in the liquid after the bubble transit. The simulation of multiple evaporating bubbles flowing in sequence shows that their mutual influence can strongly enhance the heat transfer coefficient, up to twice the single phase flow value.

Solving the cooling flow problem through mechanical AGN feedback

A fundamental gap in the current understanding of collapsed structures in the universe concerns the thermodynamical evolution of the ordinary, baryonic component. Unopposed radiative cooling of plasma would lead to the cooling catastrophe, a massive inflow of condensing gas toward the centre of galaxies, groups and clusters. The last generation of multiwavelength observations has radically changed our view on baryons, suggesting that the heating linked to the active galactic nucleus (AGN) may be the balancing counterpart of cooling. In this Thesis, I investigate the engine of the heating regulated by the central black hole. I argue that the mechanical feedback, based on massive subrelativistic outflows, is the key to solving the cooling flow problem, i.e. dramatically quenching the cooling rates for several billion years without destroying the cool-core structure. Using an upgraded version of the parallel 3D hydrodynamic code FLASH, I show that anisotropic AGN outflows can further reproduce fundamental observed features, such as buoyant bubbles, cocoon shocks, sonic ripples, metals dredge-up, and subsonic turbulence. The latter is an essential ingredient to drive nonlinear thermal instabilities, which cause cold gas condensation, a residual of the quenched cooling flow and, later, fuel for the AGN feedback engine. The self-regulated outflows are systematically tested on the scales of massive clusters, groups and isolated elliptical galaxies: in lighter less bound objects the feedback needs to be gentler and less efficient, in order to avoid drastic overheating. In this Thesis, I describe in depth the complex hydrodynamics, involving the coupling of the feedback energy to that of the surrounding hot medium. Finally, I present the merits and flaws of all the proposed models, with a critical eye toward observational concordance.

Analysis of optimal control problems for the incompressible MHD equations and implementation in a finite element multiphysics code

This thesis deals with the study of optimal control problems for the incompressible Magnetohydrodynamics (MHD) equations. Particular attention to these problems arises from several applications in science and engineering, such as fission nuclear reactors with liquid metal coolant and aluminum casting in metallurgy. In such applications it is of great interest to achieve the control on the fluid state variables through the action of the magnetic Lorentz force. In this thesis we investigate a class of boundary optimal control problems, in which the flow is controlled through the boundary conditions of the magnetic field. Due to their complexity, these problems present various challenges in the definition of an adequate solution approach, both from a theoretical and from a computational point of view. In this thesis we propose a new boundary control approach, based on lifting functions of the boundary conditions, which yields both theoretical and numerical advantages. With the introduction of lifting functions, boundary control problems can be formulated as extended distributed problems. We consider a systematic mathematical formulation of these problems in terms of the minimization of a cost functional constrained by the MHD equations. The existence of a solution to the flow equations and to the optimal control problem are shown. The Lagrange multiplier technique is used to derive an optimality system from which candidate solutions for the control problem can be obtained. In order to achieve the numerical solution of this system, a finite element approximation is considered for the discretization together with an appropriate gradient-type algorithm. A finite element object-oriented library has been developed to obtain a parallel and multigrid computational implementation of the optimality system based on a multiphysics approach. Numerical results of two- and three-dimensional computations show that a possible minimum for the control problem can be computed in a robust and accurate manner.

Model reduction of stochastic groundwater flow and transport processes

This work presents a comprehensive methodology for the reduction of analytical or numerical stochastic models characterized by uncertain input parameters or boundary conditions. The technique, based on the Polynomial Chaos Expansion (PCE) theory, represents a versatile solution to solve direct or inverse problems related to propagation of uncertainty. The potentiality of the methodology is assessed investigating different applicative contexts related to groundwater flow and transport scenarios, such as global sensitivity analysis, risk analysis and model calibration. This is achieved by implementing a numerical code, developed in the MATLAB environment, presented here in its main features and tested with literature examples. The procedure has been conceived under flexibility and efficiency criteria in order to ensure its adaptability to different fields of engineering; it has been applied to different case studies related to flow and transport in porous media. Each application is associated with innovative elements such as (i) new analytical formulations describing motion and displacement of non-Newtonian fluids in porous media, (ii) application of global sensitivity analysis to a high-complexity numerical model inspired by a real case of risk of radionuclide migration in the subsurface environment, and (iii) development of a novel sensitivity-based strategy for parameter calibration and experiment design in laboratory scale tracer transport.

Testing of subgrid scale (SGS) models for large-eddy simulation (LES) of turbulent channel flow

Sub-grid scale (SGS) models are required in order to model the influence of the unresolved small scales on the resolved scales in large-eddy simulations (LES), the flow at the smallest scales of turbulence. In the following work two SGS models are presented and deeply analyzed in terms of accuracy through several LESs with different spatial resolutions, i.e. grid spacings. The first part of this thesis focuses on the basic theory of turbulence, the governing equations of fluid dynamics and their adaptation to LES. Furthermore, two important SGS models are presented: one is the Dynamic eddy-viscosity model (DEVM), developed by \cite{germano1991dynamic}, while the other is the Explicit Algebraic SGS model (EASSM), by \cite{marstorp2009explicit}. In addition, some details about the implementation of the EASSM in a Pseudo-Spectral Navier-Stokes code \cite{chevalier2007simson} are presented. The performance of the two aforementioned models will be investigated in the following chapters, by means of LES of a channel flow, with friction Reynolds numbers $Re_\tau=590$ up to $Re_\tau=5200$, with relatively coarse resolutions. Data from each simulation will be compared to baseline DNS data. Results have shown that, in contrast to the DEVM, the EASSM has promising potentials for flow predictions at high friction Reynolds numbers: the higher the friction Reynolds number is the better the EASSM will behave and the worse the performances of the DEVM will be. The better performance of the EASSM is contributed to the ability to capture flow anisotropy at the small scales through a correct formulation for the SGS stresses. Moreover, a considerable reduction in the required computational resources can be achieved using the EASSM compared to DEVM. Therefore, the EASSM combines accuracy and computational efficiency, implying that it has a clear potential for industrial CFD usage.

Finite element modeling of flow/compression-induced deformation of alginate scaffolds for bone tissue engineering

Trauma or degenerative diseases such as osteonecrosis may determine bone loss whose recover is promised by a "tissue engineering“ approach. This strategy involves the use of stem cells, grown onboard of adequate biocompatible/bioreabsorbable hosting templates (usually defined as scaffolds) and cultured in specific dynamic environments afforded by differentiation-inducing actuators (usually defined as bioreactors) to produce implantable tissue constructs. The purpose of this thesis is to evaluate, by finite element modeling of flow/compression-induced deformation, alginate scaffolds intended for bone tissue engineering. This work was conducted at the Biomechanics Laboratory of the Institute of Biomedical and Neural Engineering of the Reykjavik University of Iceland. In this respect, Comsol Multiphysics 5.1 simulations were carried out to approximate the loads over alginate 3D matrices under perfusion, compression and perfusion+compression, when varyingalginate pore size and flow/compression regimen. The results of the simulations show that the shear forces in the matrix of the scaffold increase coherently with the increase in flow and load, and decrease with the increase of the pore size. Flow and load rates suggested for proper osteogenic cell differentiation are reported.

Differential roles for endothelial ICAM-1, ICAM-2, and VCAM-1 in shear-resistant T cell arrest, polarization, and directed crawling on blood-brain barrier endothelium

Endothelial ICAM-1 and ICAM-2 were shown to be essential for T cell diapedesis across the blood-brain barrier (BBB) in vitro under static conditions. Crawling of T cells prior to diapedesis was only recently revealed to occur preferentially against the direction of blood flow on the endothelial surface of inflamed brain microvessels in vivo. Using live cell-imaging techniques, we prove that Th1 memory/effector T cells predominantly crawl against the direction of flow on the surface of BBB endothelium in vitro. Analysis of T cell interaction with wild-type, ICAM-1-deficient, ICAM-2-deficient, or ICAM-1 and ICAM-2 double-deficient primary mouse brain microvascular endothelial cells under physiological flow conditions allowed us to dissect the individual contributions of endothelial ICAM-1, ICAM-2, and VCAM-1 to shear-resistant T cell arrest, polarization, and crawling. Although T cell arrest was mediated by endothelial ICAM-1 and VCAM-1, T cell polarization and crawling were mediated by endothelial ICAM-1 and ICAM-2 but not by endothelial VCAM-1. Therefore, our data delineate a sequential involvement of endothelial ICAM-1 and VCAM-1 in mediating shear-resistant T cell arrest, followed by endothelial ICAM-1 and ICAM-2 in mediating T cell crawling to sites permissive for diapedesis across BBB endothelium.

Intussusceptive angiogenesis: pillars against the blood flow

Adaptation of vascular networks to functional demands needs vessel growth, vessel regression and vascular remodelling. Biomechanical forces resulting from blood flow play a key role in these processes. It is well-known that metabolic stimuli, mechanical forces and flow patterns can affect gene expression and remodelling of vascular networks in different ways. For instance, in the sprouting type of angiogenesis related to hypoxia, there is no blood flow in the rising capillary sprout. In contrast, it has been shown that an increase of wall shear stress initiates the splitting type of angiogenesis in skeletal muscle. Otherwise, during development, both sprouting and intussusception act in parallel in building the vascular network, although with differences in spatiotemporal distribution. Thereby, in addition to regulatory molecules, flow dynamics support the patterning and remodelling of the rising vascular tree. Herewith, we present an overview of angiogenic processes with respect to intussusceptive angiogenesis as related to local haemodynamics.

Expression of ADAMTS1 in endothelial cells is induced by shear stress and suppressed in sprouting capillaries

ADAMTS1 inhibits capillary sprouting, and since capillary sprouts do not experience the shear stress caused by blood flow, this study undertook to clarify the relationship between shear stress and ADAMTS1. It was found that endothelial cells exposed to shear stress displayed a strong upregulation of ADAMTS1, dependent upon both the magnitude and duration of their exposure. Investigation of the underlying pathways demonstrated involvement of phospholipase C, phosphoinositide 3-kinase, and nitric oxide. Forkhead box protein O1 was identified as a likely inhibitor of the system, as its knockdown was followed by a slight increase in ADAMTS1 expression. In silico prediction displayed a transcriptional binding site for Forkhead box protein O1 in the promotor region of the ADAMTS1 gene, as well as sites for nuclear factor 1, SP1, and AP-1. The anti-angiogenic effects of ADAMTS1 were attributed to its cleavage of thrombospondin 1 into a 70-kDa fragment, and a significant enhancement of this fragment was indeed demonstrated by immunoblotting shear stress-treated cells. Accordingly, scratch wound closure displayed a slowdown in conditioned medium from shear stress-treated endothelial cells, an effect that could be completely blocked by a knockdown of thrombospondin 1 and partially blocked by a knockdown of ADAMTS1. Non-perfused capillary sprouts in rat mesenteries stained negative for ADAMTS1, while vessels in the microcirculation that had already experienced blood flow yielded the opposite results. The shear stress-dependent expression of ADAMTS1 in vitro was therefore also demonstrated in vivo and thereby confirmed as a mechanism connecting blood flow with the regulation of angiogenesis.

Shear stress modulates the expression of thrombospondin-1 and CD36 in endothelial cells in vitro and during shear stress-induced angiogenesis in vivo

Binding of thrombospondin-1 (TSP-1) to the CD36 receptor inhibits angiogenesis and induces apoptosis in endothelial cells (EC). Conversely, matrix-bound TSP-1 supports vessel formation. In this study we analyzed the shear stress-dependent expression of TSP-1 and CD36 in endothelial cells in vitro and in vivo to reveal its putative role in the blood flow-induced remodelling of vascular networks. Shear stress was applied to EC using a cone-and-plate apparatus and gene expression was analyzed by RT-PCR, Northern and Western blot. Angiogenesis in skeletal muscles of prazosin-fed (50 mg/l drinking water; 4 d) mice was assessed by measuring capillary-to-fiber (C/F) ratios. Protein expression in whole muscle homogenates (WMH) or BS-1 lectin-enriched EC fractions (ECF) was analyzed by Western blot. Shear stress downregulated TSP-1 and CD36 expression in vitro in a force- and time-dependent manner sustained for at least 72 h and reversible by restoration of no-flow conditions. In vivo, shear stress-driven increase of C/F in prazosin-fed mice was associated with reduced expression of TSP-1 and CD36 in ECF, while TSP-1 expression in WMH was increased. Down-regulation of endothelial TSP-1/CD36 by shear stress suggests a mechanism for inhibition of apoptosis in perfused vessels and pruning in the absence of flow. The increase of extra-endothelial (e.g. matrix-bound) TSP-1 could support a splitting type of vessel growth.

Successful lower extremity angioplasty improves brachial artery flow-mediated dilation in patients with peripheral arterial disease

INTRODUCTION: Peripheral arterial disease (PAD) is associated with systemic impaired flow-mediated dilation (FMD) and increased risk for cardiovascular events. Decreased FMD may be caused by a decrease in arterial shear stress due to claudication and inflammation due to muscle ischemia and reperfusion. We assumed that endovascular revascularization of lower limb arterial obstructions ameliorates FMD and lowers inflammation through improvement of peripheral perfusion. METHODS: The study was a prospective, open, randomized, controlled, single-center follow-up evaluation assessing the effect of endovascular revascularization on brachial artery reactivity (FMD) measured by ultrasound, white blood cell (WBC) count, high-sensitive C-reactive protein (hs-CRP), and fibrinogen. We investigated 33 patients (23 men) with chronic and stable PAD (Rutherford 2 to 3) due to femoropopliteal obstruction. Variables were assessed at baseline and after 4 weeks in 17 patients (group A) who underwent endovascular revascularization and best medical treatment, and in 16 patients (group B) who received best medical treatment only. RESULTS: FMD did not differ between group A and B (4.96% +/- 1.86% vs 4.60% +/- 2.95%; P = .87) at baseline. It significantly improved after revascularization in group A (6.44% +/- 2.88%; P = .02) compared with group B at 4 weeks of follow-up (4.53% +/- 3.17%; P = .92), where it remained unchanged. The baseline ankle-brachial index (ABI) was similar for group A and B (0.63 +/- 0.15 vs 0.66 +/- 0.10; P = .36). At 4 weeks of follow-up, ABI was significantly increased in group A (1.05 +/- 0.15; P = .0004) but remained unchanged in group B (0.62 +/- 0.1). WBC counts of the two groups were comparable at baseline (group A: 7.6 +/- 2.26 x 10(6)/mL and group B: 7.8 +/- 2.02 x 10(6)/mL, P = .81). In group A, the leukocyte count significantly decreased after angioplasty from 7.6 +/- 2.26 to 6.89 +/- 1.35 x 10(6)/mL (P = .03). For group B, WBC count did not differ significantly compared with baseline (7.76 +/- 2.64 x 10(6)/mL; P = .94). No effects were observed on hs-CRP or fibrinogen from endovascular therapy. CONCLUSION: Endovascular revascularization with reestablishment of peripheral arterial perfusion improves FMD and reduces WBC count in patients with claudication. Revascularization may therefore have clinical implications beyond relief of symptoms, for example, reducing oxidative stress caused by repeated muscle ischemia or increased shear stress due to improved ambulatory activity.

Shear band evolution in Zr/Hf-based bulk metallic glasses under static and dynamic indentations

Bulk metallic glasses (BMGs) exhibit superior mechanical properties as compared with other conventional materials and have been proposed for numerous engineering and technological applications. Zr/Hf-based BMGs or tungsten reinforced BMG composites are considered as a potential replacement for depleted uranium armor-piercing projectiles because of their ability to form localized shear bands during impact, which has been known to be the dominant plastic deformation mechanism in BMGs. However, in conventional tensile, compressive and bending tests, limited ductility has been observed because of fracture initiation immediately following the shear band formation. To fully investigate shear band characteristics, indentation tests that can confine the deformation in a limited region have been pursued. In this thesis, a detailed investigation of thermal stability and mechanical deformation behavior of Zr/Hf-based BMGs is conducted. First, systematic studies had been implemented to understand the influence of relative compositions of Zr and Hf on thermal stability and mechanical property evolution. Second, shear band evolution under indentations were investigated experimentally and theoretically. Three kinds of indentation studies were conducted on BMGs in the current study. (a) Nano-indentation to determine the mechanical properties as a function of Hf/Zr content. (b) Static Vickers indentation on bonded split specimens to investigate the shear band evolution characteristics beneath the indention. (c) Dynamic Vickers indentation on bonded split specimens to investigate the influence of strain rate. It was found in the present work that gradually replacing Zr by Hf remarkably increases the density and improves the mechanical properties. However, a slight decrease in glass forming ability with increasing Hf content has also been identified through thermodynamic analysis although all the materials in the current study were still found to be amorphous. Many indentation studies have revealed only a few shear bands surrounding the indent on the top surface of the specimen. This small number of shear bands cannot account for the large plastic deformation beneath the indentations. Therefore, a bonded interface technique has been used to observe the slip-steps due to shear band evolution. Vickers indentations were performed along the interface of the bonded split specimen at increasing loads. At small indentation loads, the plastic deformation was primarily accommodated by semi-circular primary shear bands surrounding the indentation. At higher loads, secondary and tertiary shear bands were formed inside this plastic zone. A modified expanding cavity model was then used to predict the plastic zone size characterized by the shear bands and to identify the stress components responsible for the evolution of the various types of shear bands. The applicability of various hardness—yield-strength ( H −σγ ) relationships currently available in the literature for bulk metallic glasses (BMGs) is also investigated. Experimental data generated on ZrHf-based BMGs in the current study and those available elsewhere on other BMG compositions were used to validate the models. A modified expanding-cavity model, employed in earlier work, was extended to propose a new H −σγ relationship. Unlike previous models, the proposed model takes into account not only the indenter geometry and the material properties, but also the pressure sensitivity index of the BMGs. The influence of various model parameters is systematically analyzed. It is shown that there is a good correlation between the model predictions and the experimental data for a wide range of BMG compositions. Under dynamic Vickers indentation, a decrease in indentation hardness at high loading rate was observed compared to static indentation hardness. It was observed that at equivalent loads, dynamic indentations produced more severe deformation features on the loading surface than static indentations. Different from static indentation, two sets of widely spaced semi-circular shear bands with two different curvatures were observed. The observed shear band pattern and the strain rate softening in indentation hardness were rationalized based on the variations in the normal stress on the slip plane, the strain rate of shear and the temperature rise associated with the indentation deformation. Finally, a coupled thermo-mechanical model is proposed that utilizes a momentum diffusion mechanism for the growth and evolution of the final spacing of shear bands. The influence of strain rate, confinement pressure and critical shear displacement on the shear band spacing, temperature rise within the shear band, and the associated variation in flow stress have been captured and analyzed. Consistent with the known pressure sensitive behavior of BMGs, the current model clearly captures the influence of the normal stress in the formation of shear bands. The normal stress not only reduces the time to reach critical shear displacement but also causes a significant temperature rise during the shear band formation. Based on this observation, the variation of shear band spacing in a typical dynamic indentation test has been rationalized. The temperature rise within a shear band can be in excess of 2000K at high strain rate and high confinement pressure conditions. The associated drop in viscosity and flow stress may explain the observed decrease in fracture strength and indentation hardness. The above investigations provide valuable insight into the deformation behavior of BMGs under static and dynamic loading conditions. The shear band patterns observed in the above indentation studies can be helpful to understand and model the deformation features under complex loading scenarios such as the interaction of a penetrator with armor. Future work encompasses (1) extending and modifying the coupled thermo-mechanical model to account for the temperature rise in quasistatic deformation; and (2) expanding this model to account for the microstructural variation-crystallization and free volume migration associated with the deformation.