Smaller Deborah number inducing more serrated plastic flow of metallic glass

Spherical nano-indentations of Cu46Zr54 bulk metallic glass (BMG) model systems were performed using molecular dynamics (MD) computer simulations, focusing specifically on the physical origin of serrated plastic flow. The results demonstrate that there is a direct correlation between macroscopic flow serration and underlying irreversible rearrangement of atoms, which is strongly dependent on the loading (strain) rate and the temperature. The serrated plastic flow is, therefore, determined by the magnitude of such irreversible rearrangement that is inhomogeneous temporally. A dimensionless Deborah number is introduced to characterize the effects of strain rate and temperature on serrations. Our simulations are shown to compare favorably with the available experimental observations.

Towards understanding the mixing characteristics of turbulent buoyant flows

This work proposes a new simulation methodology in which variable density turbulent flows can be studied in the context of a mixing layer with or without the presence of gravity. Specifically, this methodology is developed to probe the nature of non-buoyantly-driven (i.e. isotropically-driven) or buoyantly-driven mixing deep inside a mixing layer. Numerical forcing methods are incorporated into both the velocity and scalar fields, which extends the length of time over which mixing physics can be studied. The simulation framework is designed to allow for independent variation of four non-dimensional parameters, including the Reynolds, Richardson, Atwood, and Schmidt numbers. Additionally, the governing equations are integrated in such a way to allow for the relative magnitude of buoyant energy production and non-buoyant energy production to be varied.

The computational requirements needed to implement the proposed configuration are presented. They are justified in terms of grid resolution, order of accuracy, and transport scheme. Canonical features of turbulent buoyant flows are reproduced as validation of the proposed methodology. These features include the recovery of isotropic Kolmogorov scales under buoyant and non-buoyant conditions, the recovery of anisotropic one-dimensional energy spectra under buoyant conditions, and the preservation of known statistical distributions in the scalar field, as found in other DNS studies.

This simulation methodology is used to perform a parametric study of turbulent buoyant flows to discern the effects of varying the Reynolds, Richardson, and Atwood numbers on the resulting state of mixing. The effects of the Reynolds and Atwood numbers are isolated by looking at two energy dissipation rate conditions under non-buoyant (variable density) and constant density conditions. The effects of Richardson number are isolated by varying the ratio of buoyant energy production to total energy production from zero (non-buoyant) to one (entirely buoyant) under constant Atwood number, Schmidt number, and energy dissipation rate conditions. It is found that the major differences between non-buoyant and buoyant turbulent flows are contained in the transfer spectrum and longitudinal structure functions, while all other metrics are largely similar (e.g. energy spectra, alignment characteristics of the strain-rate tensor). Also, despite the differences noted between fully buoyant and non-buoyant turbulent fields, the scalar field, in all cases, is unchanged by these. The mixing dynamics in the scalar field are found to be insensitive to the source of turbulent kinetic energy production (non-buoyant vs. buoyant).

Small scale turbulence in high Karlovitz number premixed flames

The purpose of this thesis is to characterize the behavior of the smallest turbulent scales in high Karlovitz number (Ka) premixed flames. These scales are particularly important in the two-way coupling between turbulence and chemistry and better understanding of these scales will support future modeling efforts using large eddy simulations (LES). The smallest turbulent scales are studied by considering the vorticity vector, ω, and its transport equation.

Due to the complexity of turbulent combustion introduced by the wide range of length and time scales, the two-dimensional vortex-flame interaction is first studied as a simplified test case. Numerical and analytical techniques are used to discern the dominate transport terms and their effects on vorticity based on the initial size and strength of the vortex. This description of the effects of the flame on a vortex provides a foundation for investigating vorticity in turbulent combustion.

Subsequently, enstrophy, ω² = ω • ω, and its transport equation are investigated in premixed turbulent combustion. For this purpose, a series of direct numerical simulations (DNS) of premixed n-heptane/air flames are performed, the conditions of which span a wide range of unburnt Karlovitz numbers and turbulent Reynolds numbers. Theoretical scaling analysis along with the DNS results support that, at high Karlovitz number, enstrophy transport is controlled by the viscous dissipation and vortex stretching/production terms. As a result, vorticity scales throughout the flame with the inverse of the Kolmogorov time scale, τ_η, just as in homogeneous isotropic turbulence. As τ_η is only a function of the viscosity and dissipation rate, this supports the validity of Kolmogorov’s first similarity hypothesis for sufficiently high Ka numbers (Ka ≳ 100). These conclusions are in contrast to low Karlovitz number behavior, where dilatation and baroclinic torque have a significant impact on vorticity within the flame. Results are unaffected by the transport model, chemical model, turbulent Reynolds number, and lastly the physical configuration.

Next, the isotropy of vorticity is assessed. It is found that given a sufficiently large value of the Karlovitz number (Ka ≳ 100) the vorticity is isotropic. At lower Karlovitz numbers, anisotropy develops due to the effects of the flame on the vortex stretching/production term. In this case, the local dynamics of vorticity in the strain-rate tensor, S, eigenframe are altered by the flame. At sufficiently high Karlovitz numbers, the dynamics of vorticity in this eigenframe resemble that of homogeneous isotropic turbulence.

Combined, the results of this thesis support that both the magnitude and orientation of vorticity resemble the behavior of homogeneous isotropic turbulence, given a sufficiently high Karlovitz number (Ka ≳ 100). This supports the validity of Kolmogorov’s first similarity hypothesis and the hypothesis of local isotropy under these condition. However, dramatically different behavior is found at lower Karlovitz numbers. These conclusions provides/suggests directions for modeling high Karlovitz number premixed flames using LES. With more accurate models, the design of aircraft combustors and other combustion based devices may better mitigate the detrimental effects of combustion, from reducing CO₂ and soot production to increasing engine efficiency.

Modification of the biological intercept model to account for ontogenetic effects in laboratory-reared delta smelt (Hypomesus transpacificus)*

We investigated age, growth, and ontogenetic effects on the proportionality of otolith size to fish size in laboratory-reared delta smelt (Hypomesus transpacificus) from the San Francisco Bay estuary. Delta smelt larvae were reared from hatching in laboratory mesocosms for 100 days. Otolith increments from known-age fish were enumerated to validate that growth increments were deposited daily and to validate the age of fish at first ring formation. Delta smelt were found to lay down daily ring increments; however, the first increment did not form until six days after hatching. The relationship between otolith size and fish size was not biased by age or growth-rate effects but did exhibit an interruption in linear growth owing to an ontogenetic shift at the postflexon stage. To back-calculate the size-at-age of individual fish, we modified the biological intercept (BI) model to account for ontogenetic changes in the otolith-size−fish-size relationship and compared the results to the time-varying growth model, as well as the modified Fry model. We found the modified BI model estimated more accurately the size-at-age from hatching to 100 days after hatching. Before back-calculating size-at-age with existing models, we recommend a critical evaluation of the effects that age, growth, and ontogeny can have on the otolith-size−fish-size relations

Yielding of engineering polymers at strain rates of up to 500 s-1

A method is described for measuring the mechanical properties of polymers in compression at strain rates in the range approximately 300-500 s-1. A gravity-driven pendulum is used to load a specimen on the end of an instrumented Hopkinson output bar and the results are processed by a microcomputer. Stress-strain curves up to high strains are presented for polycarbonate, polyethersulphone and high density polyethylene over a range of temperatures. The value of yield stress, for all three polymers, was found to vary linearly with log (strain rate) at strain rates up to 500 s-1. © 1985.

Dynamic compressive response of composite square honeycombs

Carbon fibre-epoxy composite square honeycombs, and the parent composite material, were tested in quasi-static compression at a strain rate of 10 -3 s -1 and in dynamic compression at strain rates of 10 3-10 4 s -1 using an instrumented Kolsky bar arrangement. Taken together, these tests provide an assessment of the potential of this composite topology for use as a lightweight sandwich core. The honeycombs had two relative densities, 0.12 and 0.24, and two material orientations, ±45° and 0/90° with respect to the prismatic, loading direction of the honeycomb. Honeycomb manufacture was by slotting, assembling and bonding together carbon fibre/epoxy woven plies of composite sheets of 2 × 2 twill weave construction. The peak value of wall stress in the honeycombs was about one third that of the parent material, for all strain rates. An elastic finite element analysis was used to trace the source of this knock-down in strength: a stress concentration exists at the root of the slots and leads to premature failure by microbuckling. Shock-wave effects were evident at impact velocities exceeding 50 ms -1 for the honeycomb of relative density 0.12. This was traced to stubbing of the buckled cell walls against the face of the Kolsky bar. © 2011 Elsevier Ltd. All rights reserved.

Study of the turbulence in the air-side and water-side boundary layers in experimental laboratory wind induced surface waves

This study detailed the structure of turbulence in the air-side and water-side boundary layers in wind-induced surface waves. Inside the air boundary layer, the kurtosis is always greater than 3 (the value for normal distribution) for both horizontal and vertical velocity fluctuations. The skewness for the horizontal velocity is negative, but the skewness for the vertical velocity is always positive. On the water side, the kurtosis is always greater than 3, and the skewness is slightly negative for the horizontal velocity and slightly positive for the vertical velocity. The statistics of the angle between the instantaneous vertical fluctuation and the instantaneous horizontal velocity in the air is similar to those obtained over solid walls. Measurements in water show a large variance, and the peak is biased towards negative angles. In the quadrant analysis, the contribution of quadrants Q2 and Q4 is dominant on both the air side and the water side. The non-dimensional relative contributions and the concentration match fairly well near the interface. Sweeps in the air side (belonging to quadrant Q4) act directly on the interface and exert pressure fluctuations, which, in addition to the tangential stress and form drag, lead to the growth of the waves. The water drops detached from the crest and accelerated by the wind can play a major role in transferring momentum and in enhancing the turbulence level in the water side.On the air side, the Reynolds stress tensor's principal axes are not collinear with the strain rate tensor, and show an angle α σ≈=-20°to-25°. On the water side, the angle is α σ≈=-40°to-45°. The ratio between the maximum and the minimum principal stresses is σ a/σ b=3to4 on the air side, and σ a/σ b=1.5to3 on the water side. In this respect, the air-side flow behaves like a classical boundary layer on a solid wall, while the water-side flow resembles a wake. The frequency of bursting on the water side increases significantly along the flow, which can be attributed to micro-breaking effects - expected to be more frequent at larger fetches. © 2012 Elsevier B.V.

Fracture mechanics of idealised bituminous mixes

© 2014 Taylor & Francis. The durability of asphalt pavements is strongly impaired by cracks, caused primarily by traffic loads and environmental effects. In this work, fracture behaviour of idealised asphalt mixes is investigated. Experiments on idealised asphalt mixes under pure-tension mode (mode I cracking) were performed and fracture parameters were evaluated. In these three-point bend fracture tests, the test variables were temperature and load rate. The test data were stored in an asphalt materials database and special-purpose tools were implemented to analyse and handle the laboratory data automatically. Fracture mechanism maps were constructed, showing the conditions associated with ductile, brittle and ductile-brittle transition regimes of behaviour. The mechanism maps show the failure response of the material in terms of the stress intensity factor, strain energy release rate and J-integral as a function of the temperature-compensated crack mouth opening strain rate. Fracture behaviour of asphalt mix specimens was simulated by cohesive zone model in conjunction with a novel material constitutive model for asphalt mixes. The finite element model agrees well with the experimental results and provides insights into fracture response of the notched asphalt mix beam specimens.

Effect of a low-temperature thin buffer layer on the strain accommodation of In0.25Ga0.75As grown on a GaAs(001) substrate

Various low-temperature (LT) ultra-thin buffer layers have been fabricated on the GaAs (001) substrate. The buffer layer is decoupled from the host substrate by introducing low-temperature defects. The 400 nm In0.25Ga0.75As films were grown on these substrates to test the 'compliant' effects of the buffer layers. Atomic force microscopy, photoluminescence, double crystal x-ray diffraction and transmission electron microscopy were used to estimate the quality of the ln(0.25)Ga(0.75)As layer. The measurements indicated that the misfit strains in the epilayer can be accommodated by the LT ultra-thin buffer layer. The strain accommodation effects of the LT defects have been discussed in detail.

Effects of annealing time and Si cap layer thickness on the Si/SiGe/Si heterostructures thermal stability

The effects of annealing time and Si cap layer thickness: on the thermal stability of the Si/SiGe/Si heterostructures deposited by disilane and solid-Ge molecule beam epitaxy were investigated. It is found that in the same strain state of the SiGe layers the annealing time decreases with increasing Si cap layer thickness. This effect is analyzed by a force-balance theory and an equation has been obtained to characterize the relation between the annealing time and the Si cap layer thickness. (C) 2001 Elsevier Science B.V. All rights reserved.

Effects of annealing time and Si cap layer thickness on the Si/SiGe/Si heterostructures thermal stability

The effects of annealing time and Si cap layer thickness: on the thermal stability of the Si/SiGe/Si heterostructures deposited by disilane and solid-Ge molecule beam epitaxy were investigated. It is found that in the same strain state of the SiGe layers the annealing time decreases with increasing Si cap layer thickness. This effect is analyzed by a force-balance theory and an equation has been obtained to characterize the relation between the annealing time and the Si cap layer thickness. (C) 2001 Elsevier Science B.V. All rights reserved.

Brittle-ductile transition of polymers and its percolation model

The brittle-ductile transition (BDT) of particle toughened polymers was extensively studied in terms of morphology, strain rate, and temperature. The calculation results showed that both the critical interparticle distance (IDc) and the brittle-ductile transition temperature (T-BD) of polymers were a function of strain rate. The IDc reduced nonlinearly with increasing strain rate, whereas T-BD increased considerably with increasing strain rate. The effects of temperature and plasticizer concentration on BDT were discussed using a percolation model. The results were in agreement with the experiments.

Stress corrosion cracking of AISI 321 stainless steel in acidic chloride solution

The stress corrosion cracking (SCC) of LambdaISI 321 stainless steel in acidic chloride solution was studied by slow strain rate (SSR) technique and fracture mechanics method. The fractured surface was characterized by cleavage fracture. In order to clarify the SCC mechanism, the effects of inhibitor KI on SCC behaviour were also included in this paper. A study showed that the inhibition effects of KI on SCC were mainly attributed to the anodic reaction of the corrosion process. The results of strain distribution in front of the crack tip of the fatigue pre-cracked plate specimens in air, in the blank solution (acidic chloride solution without inhibitor KI) and in the solution added with KI measured by speckle interferometry (SPI) support the unified mechanism of SCC and corrosion fatigue cracking (CFC).

Effects of sulfate-reducing bacteria on the corrosion behavior of carbon steel

The influences of the growing process of sulfate-reducing bacteria (SRB) in seawater system on the medium state and corrosion behavior of carbon steel were studied by detecting solution state parameters and using corrosion electrochemical methods. The growing process of SRB in the seawater shows the three stages of growing, death and residual phases. The solution state parameters of the concentration of sulfide, the pH value and the redox potential changed during the three stages of the SRB growing process. And the corrosion rate of D36 carbon steel was accelerated during the growing phase and stable during the death and residual phases. The results indicate that the medium state and the corrosion rate of the steel do not depend on the number of active SRB, but depend on the accumulation of the metabolism products of SRB. (c) 2007 Elsevier Ltd. All rights reserved.

A viscous froth model for dry foams in the Surface Evolver

Cox, S.J. (2005) A Viscous Froth Model for Dry Foams in the Surface Evolver. Coll. Surf. A 263:81-89. Sponsorship: Partial support from both the Ulysses France?Ireland Exchange Scheme and the European Space Agency, Contract 14308/00/NL/SH (AO-99-031) CCN 002 MAP Project AO-99-075.