A Multiscale-Domain Decomposition Method for High Reynolds Number Flows

The high Reynolds number flow contains a wide range of length and time scales, and the flow domain can be divided into several sub-domains with different characteristic scales. In some sub-domains, the viscosity dissipation scale can only be considered in a certain direction; in some sub-domains, the viscosity dissipation scales need to be considered in all directions; in some sub-domains, the viscosity dissipation scales are unnecessary to be considered at all. For laminar boundary layer region, the characteristic length scales in the streamwise and normal directions are L and L Re-1/ 2 , respectively. The characteristic length scale and the velocity scale in the outer region of the boundary layer are L and U, respectively. In the neighborhood region of the separated point, the length scale l<scale u<scale-domain decomposition method is developed for the high Reynolds number flow. First, the whole domain is decomposed to different sub-domains with the different characteristic scales. Then the different dominant equation of all sub-domains is defined according to the diffusion parabolized (DP) theory of viscous flow. Finally these different equations are solved simultaneously in whole computational region. For numerical tests of high Reynolds numerical flows, two-dimensional supersonic flows over rearward and frontward steps as well as an interaction flow between shock wave and boundary layer were solved numerically. The pressure distributions and local coefficients of skin friction on the wall are given. The numerical results obtained by the multiscale-domain decomposition algorithm are well agreement with those by NS equations. Comparing with the usual method of solving the Navier-Stokes equations in the whole flow, under the same numerical accuracy, the present multiscale domain decomposition method decreases CPU consuming about 20% and reflects the physical mechanism of practical flow more accurately.

Simulation of gas flow over a micro cylinder

Gas flow over a micro cylinder is simulated using both a compressible Navier-Stokes solver and a hybrid continuum /particle approach. The micro cylinder flow has low Reynolds number because of the small length scale and the low speed, which also indicates that the rarefied gas effect exists in the flow. A cylinder having a diameter of 20 microns is simulated under several flow conditions where the Reynolds number ranges from 2 to 50 and the Mach number varies from 0.1 to 0.8. It is found that the low Reynolds number flow can be compressible even when the Mach number is less than 0.3, and the drag coefficient of the cylinder increases when the Reynolds number decreases. The compressible effect will increase the pressure drag coefficient although the friction coefficient remains nearly unchanged. The rarefied gas effect will reduce both the friction and pressure drag coefficients, and the vortex in the flow may be shrunk or even disappear.

The short-timescale behavior of glacial ice

Glaciers are often assumed to deform only at slow (i.e., glacial) rates. However, with the advent of high rate geodetic observations of ice motion, many of the intricacies of glacial deformation on hourly and daily timescales have been observed and quantified. This thesis explores two such short timescale processes: the tidal perturbation of ice stream motion and the catastrophic drainage of supraglacial meltwater lakes. Our investigation into the transmission length-scale of a tidal load represents the first study to explore the daily tidal influence on ice stream motion using three-dimensional models. Our results demonstrate both that the implicit assumptions made in the standard two-dimensional flow-line models are inherently incorrect for many ice streams, and that the anomalously large spatial extent of the tidal influence seen on the motion of some glaciers cannot be explained, as previously thought, through the elastic or viscoelastic transmission of tidal loads through the bulk of the ice stream. We then discuss how the phase delay between a tidal forcing and the ice stream’s displacement response can be used to constrain in situ viscoelastic properties of glacial ice. Lastly, for the problem of supraglacial lake drainage, we present a methodology for implementing linear viscoelasticity into an existing model for lake drainage. Our work finds that viscoelasticity is a second-order effect when trying to model the deformation of ice in response to a meltwater lake draining to a glacier’s bed. The research in this thesis demonstrates that the first-order understanding of the short-timescale behavior of naturally occurring ice is incomplete, and works towards improving our fundamental understanding of ice behavior over the range of hours to days.

Optimal scaling in ductile fracture

This work is concerned with the derivation of optimal scaling laws, in the sense of matching lower and upper bounds on the energy, for a solid undergoing ductile fracture. The specific problem considered concerns a material sample in the form of an infinite slab of finite thickness subjected to prescribed opening displacements on its two surfaces. The solid is assumed to obey deformation-theory of plasticity and, in order to further simplify the analysis, we assume isotropic rigid-plastic deformations with zero plastic spin. When hardening exponents are given values consistent with observation, the energy is found to exhibit sublinear growth. We regularize the energy through the addition of nonlocal energy terms of the strain-gradient plasticity type. This nonlocal regularization has the effect of introducing an intrinsic length scale into the energy. We also put forth a physical argument that identifies the intrinsic length and suggests a linear growth of the nonlocal energy. Under these assumptions, ductile fracture emerges as the net result of two competing effects: whereas the sublinear growth of the local energy promotes localization of deformation to failure planes, the nonlocal regularization stabilizes this process, thus resulting in an orderly progression towards failure and a well-defined specific fracture energy. The optimal scaling laws derived here show that ductile fracture results from localization of deformations to void sheets, and that it requires a well-defined energy per unit fracture area. In particular, fractal modes of fracture are ruled out under the assumptions of the analysis. The optimal scaling laws additionally show that ductile fracture is cohesive in nature, i.e., it obeys a well-defined relation between tractions and opening displacements. Finally, the scaling laws supply a link between micromechanical properties and macroscopic fracture properties. In particular, they reveal the relative roles that surface energy and microplasticity play as contributors to the specific fracture energy of the material. Next, we present an experimental assessment of the optimal scaling laws. We show that when the specific fracture energy is renormalized in a manner suggested by the optimal scaling laws, the data falls within the bounds predicted by the analysis and, moreover, they ostensibly collapse---with allowances made for experimental scatter---on a master curve dependent on the hardening exponent, but otherwise material independent.

Large-eddy simulations of fully developed turbulent channel and pipe flows with smooth and rough walls

Studies in turbulence often focus on two flow conditions, both of which occur frequently in real-world flows and are sought-after for their value in advancing turbulence theory. These are the high Reynolds number regime and the effect of wall surface roughness. In this dissertation, a Large-Eddy Simulation (LES) recreates both conditions over a wide range of Reynolds numbers Re_τ = O(10²)-O(10⁸) and accounts for roughness by locally modeling the statistical effects of near-wall anisotropic fine scales in a thin layer immediately above the rough surface. A subgrid, roughness-corrected wall model is introduced to dynamically transmit this modeled information from the wall to the outer LES, which uses a stretched-vortex subgrid-scale model operating in the bulk of the flow. Of primary interest is the Reynolds number and roughness dependence of these flows in terms of first and second order statistics. The LES is first applied to a fully turbulent uniformly-smooth/rough channel flow to capture the flow dynamics over smooth, transitionally rough and fully rough regimes. Results include a Moody-like diagram for the wall averaged friction factor, believed to be the first of its kind obtained from LES. Confirmation is found for experimentally observed logarithmic behavior in the normalized stream-wise turbulent intensities. Tight logarithmic collapse, scaled on the wall friction velocity, is found for smooth-wall flows when Re_τ ≥ O(10⁶) and in fully rough cases. Since the wall model operates locally and dynamically, the framework is used to investigate non-uniform roughness distribution cases in a channel, where the flow adjustments to sudden surface changes are investigated. Recovery of mean quantities and turbulent statistics after transitions are discussed qualitatively and quantitatively at various roughness and Reynolds number levels. The internal boundary layer, which is defined as the border between the flow affected by the new surface condition and the unaffected part, is computed, and a collapse of the profiles on a length scale containing the logarithm of friction Reynolds number is presented. Finally, we turn to the possibility of expanding the present framework to accommodate more general geometries. As a first step, the whole LES framework is modified for use in the curvilinear geometry of a fully-developed turbulent pipe flow, with implementation carried out in a spectral element solver capable of handling complex wall profiles. The friction factors have shown favorable agreement with the superpipe data, and the LES estimates of the Karman constant and additive constant of the log-law closely match values obtained from experiment.

Fabrication, Characterization, And Deformation of 3D Structural Meta-Materials

Current technological advances in fabrication methods have provided pathways to creating architected structural meta-materials similar to those found in natural organisms that are structurally robust and lightweight, such as diatoms. Structural meta-materials are materials with mechanical properties that are determined by material properties at various length scales, which range from the material microstructure (nm) to the macro-scale architecture (μm – mm). It is now possible to exploit material size effect, which emerge at the nanometer length scale, as well as structural effects to tune the material properties and failure mechanisms of small-scale cellular solids, such as nanolattices. This work demonstrates the fabrication and mechanical properties of 3-dimensional hollow nanolattices in both tension and compression. Hollow gold nanolattices loaded in uniaxial compression demonstrate that strength and stiffness vary as a function of geometry and tube wall thickness. Structural effects were explored by increasing the unit cell angle from 30° to 60° while keeping all other parameters constant; material size effects were probed by varying the tube wall thickness, t, from 200nm to 635nm, at a constant relative density and grain size. In-situ uniaxial compression experiments reveal an order-of-magnitude increase in yield stress and modulus in nanolattices with greater lattice angles, and a 150% increase in the yield strength without a concomitant change in modulus in thicker-walled nanolattices for fixed lattice angles. These results imply that independent control of structural and material size effects enables tunability of mechanical properties of 3-dimensional architected meta-materials and highlight the importance of material, geometric, and microstructural effects in small-scale mechanics. This work also explores the flaw tolerance of 3D hollow-tube alumina kagome nanolattices with and without pre-fabricated notches, both in experiment and simulation. Experiments demonstrate that the hollow kagome nanolattices in uniaxial tension always fail at the same load when the ratio of notch length (a) to sample width (w) is no greater than 1/3, with no correlation between failure occurring at or away from the notch. For notches with (a/w) > 1/3, the samples fail at lower peak loads and this is attributed to the increased compliance as fewer unit cells span the un-notched region. Finite element simulations of the kagome tension samples show that the failure is governed by tensile loading for (a/w) < 1/3 but as (a/w) increases, bending begins to play a significant role in the failure. This work explores the flaw sensitivity of hollow alumina kagome nanolattices in tension, using experiments and simulations, and demonstrates that the discrete-continuum duality of architected structural meta-materials gives rise to their flaw insensitivity even when made entirely of intrinsically brittle materials.

Topology optimization of silicon anode structures for lithium-ion battery applications

This thesis presents a topology optimization methodology for the systematic design of optimal multifunctional silicon anode structures in lithium-ion batteries. In order to develop next generation high performance lithium-ion batteries, key design challenges relating to the silicon anode structure must be addressed, namely the lithiation-induced mechanical degradation and the low intrinsic electrical conductivity of silicon. As such, this work considers two design objectives of minimum compliance under design dependent volume expansion, and maximum electrical conduction through the structure, both of which are subject to a constraint on material volume. Density-based topology optimization methods are employed in conjunction with regularization techniques, a continuation scheme, and mathematical programming methods. The objectives are first considered individually, during which the iteration history, mesh independence, and influence of prescribed volume fraction and minimum length scale are investigated. The methodology is subsequently extended to a bi-objective formulation to simultaneously address both the compliance and conduction design criteria. A weighting method is used to derive the Pareto fronts, which demonstrate a clear trade-off between the competing design objectives. Furthermore, a systematic parameter study is undertaken to determine the influence of the prescribed volume fraction and minimum length scale on the optimal combined topologies. The developments presented in this work provide a foundation for the informed design and development of silicon anode structures for high performance lithium-ion batteries.

Double proximity effect in lateral S-NF-S junctions

We investigate planar Josephson junctions where the intermediate spacer between the two superconductors is an hybrid structure made by a normal metal and a ferromagnet. The different behaviors of the S-N-S junctions with thicknesses of 50 nm in both Cu and Nb layers, and S-N/F-S junctions with 10 nm of Co, 50 nm of Cu and 50 nm of Nb are studied. In this way, we analyze the influence of the ferromagnetic exchange interaction on the proximity effect. A dramatic supression of the josephson critical current of the Nb-(Cu/Co)-Nb junctions is observed. We believe that the reason for this is due to the length scale of the superconducting correlations of the electrons and holes of the weak link is larger than the thickness of Cu/Co bilayer.

Prediction of sound emission from open turbulent premixed flames

The sound emission from open turbulent flames is dictated by the two-point spatial correlation of rate of change of fluctuating heat release rate and this correlation has not been investigated directly in the past studies. Turbulent premixed flame data from DNS and laser diagnostics are analyzed to study this correlation function and the two-point spatial correlation of the fluctuating heat release rate. This shows that the correlation functions have simple Gaussian forms whose integral length scale is related to the laminar flame thickness and amplitude depends on the spatial distribution of the time-mean rate of heat release. These results and RANS-CFD solution of open turbulent premixed flames are post-processed to obtain the far field SPL, which agrees well with measured values. © 2010 by the American Institute of Aeronautics and Astronautics, Inc.

Effects of rotating inlet distortion on multistage compressor stability

In multi-spool engines, rotating stall in an upstream compressor will impose a rotating distortion on the downstream compressor, thereby affecting its stability margin. In this paper experiments are described in which this effect was simulated by a rotating screen upstream of several multistage low-speed compressors. The measurements are complemented by, and compared with, a theoretical model of multistage compressor response to speed and direction of rotation of an inlet distortion. For co-rotating distortions (i.e., distortions rotating in the same direction as rotor rotation), experiments show that the compressors exhibited significant loss in stability margin and that they could be divided into two groups according to their response. The first group exhibited a single peak in stall margin degradation when the distortion speed corresponded to roughly 50% of rotor speed. The second group showed two peaks in stall margin degradation corresponding to distortion speeds of approximately 25-35% and 70-75% of rotor speed. These new results demonstrate that multistage compressors can have more than a single resonant response. Detailed measurements suggest that the two types of behavior are linked to differences between the stall inception processes observed for the two groups of compressors and that a direct connection thus exists between the observed forced response and the unsteady flow phenomena at stall onset. For counter-rotational distortions, all the compressors tested showed minimal loss of stability margin. The results imply that counter-rotation of the fan and core compressor, or LP and HP compressors, could be a worthwhile design choice. Calculations based on the two-dimensional theoretical model show excellent agreement for the compressors which had a single peak for stall margin degradation. We take this first-of-a-kind comparison as showing that the model, though simplified, captures the essential fluid dynamic features of the phenomena. Agreement is not good for compressors which had two peaks in the curve of stall margin shift versus distortion rotation speed. The discrepancy is attributed to the three-dimensional and short length scale nature of the stall inception process in these machines; this includes phenomena that have not yet been addressed in any model.

Compressor wake/leading-edge interactions at off design incidences

In this paper, the effects of wake/leading-edge interactions were studied at off-design conditions. Measurements were performed on the stator-blade suction surface at midspan. The leading-edge flow-field was investigated using hotwire micro-traverses, hotfilm surface shear-stress sensors and pressure micro-tappings. The trailing-edge flow-field was investigated using hotwire boundary-layer traverses. Unsteady CFD calculations were also performed to aid the interpretation of the results. At low flow coefficients, the time-averaged momentum thickness of the leading-edge boundary layer was found to rise as the flow coefficient was reduced. The time-resolved momentum-thickness rose due to the interaction of the incoming rotor wake. As the flow coefficient was reduced, the incoming wakes increased in pitch-wise extent, velocity deficit and turbulence intensity. This increased both the time-resolved rise in the momentum thickness and the turbulent spot production within the wake affected boundary-layer. Close to stall, a drop in the leading-edge momentum thickness was observed in-between wake events. This was associated with the formation of a leading-edge separation bubble in-between wake events. The wake interaction with the bubble gave rise to a shedding phenomenon, which produced large length scale disturbances in the surface shear stress. Copyright © 2008 by ASME.

Size effects in indentation of hydrated biological tissues

Fluid flow in biological tissues is important in both mechanical and biological contexts. Given the hierarchical nature of tissues, there are varying length scales at which time-dependent mechanical behavior due to fluid flow may be exhibited. Here, spherical nanoindentation and microindentation testings are used for the characterization of length scale effects in the mechanical response of hydrated tissues. Although elastic properties were consistent across length scales, there was a substantial difference between the time-dependent mechanical responses for large and small contact radii in the same tissue specimens. This difference was far more obvious when poroelastic analysis was used instead of viscoelastic analysis. Overall, indentation testing is a fast and robust technique for characterizing the hierarchical structure of biological materials from nanometer to micrometer length scales and is capable of making quantitative material property measurements to do with fluid flow. © 2011 Materials Research Society.

Patterning of crystalline organic materials by electro-hydrodynamic lithography.

The control of semi-crystalline polymers in thin films and in micrometer-sized patterns is attractive for (opto-)electronic applications. Electro-hydrodynamic lithography (EHL) enables the structure formation of organic crystalline materials on the micrometer length scale while at the same time exerting control over crystal orientation. This gives rise to well-defined micro-patterned arrays of uniaxially aligned polymer crystals. This study explores the interplay of EHL structure formation with crystal alignment and studies the mechanisms that give rise to crystal orientation in EHL-generated structures.

Spatial correlation of heat release rate and sound emission from turbulent premixed flames

The two-point spatial correlation of the rate of change of fluctuating heat release rate is central to the sound emission from open turbulent flames, and a few attempts have been made to address this correlation in recent studies. In this paper, the two-point correlation and its role in combustion noise are studied by analysing direct numerical simulation (DNS) data of statistically multi-dimensional turbulent premixed flames. The results suggest that this correlation function depends on the separation distance and direction but, not on the positions inside the flame brush. This correlation can be modelled using a combination of Hermite-Gaussian functions of zero and second order, i.e. functions of the form (1-Ax2)e-Bx2 for constants A and B, to include its possible negative values. The integral correlation volume obtained using this model is about 0.2δL3 with the length scale obtained from its cube root being about 0.6δ L, where δ L is the laminar flame thermal thickness. Both of the values are slightly larger than the values reported in an earlier study because of the anisotropy observed for the correlation. This model together with the turbulence-dependent parameter K, the ratio of the root-mean-square (RMS) value of the rate of change of reaction rate to the mean reaction rate, derived from the DNS data is applied to predict the far-field sound emitted from open flames. The calculated noise levels agree well with recently reported measurements and show a sensitivity to K values. © 2012 The Combustion Institute.

Fracture toughness of fibrous membranes

Random fibrous networks exist in both natural biological and engineering materials. While the nonlinear deformation of fibrous networks has been extensively studied, the understanding of their fracture behaviour is still incomplete. To study the fracture toughness of fibrous materials, the near-tip region is crucial because failure mechanisms such as fibril rupture occur in this region. The consideration of this region in fracture studies is, however, a difficult task because it involves microscopic mechanical responses at a small length scale. This paper extends our previous finite element analysis by incorporating the microscopic responses into a macroscopic domain by using a submodeling technique. The detailed study of microstructures at crack tips show a stochastic toughness of membranes due to the random nature of fibrous networks. Further, the sizes of crack tip region, which are sufficient to provide a reasonable prediction of fracture behaviour in a specific type of fibrous network, were presented. Future work includes improving the current linear assumption in the macroscopic models to become nonlinear.