Higher-order analysis of crack tip fields in elastic power-law hardening materials

A HIGHER-ORDER asymptotic analysis of a stationary crack in an elastic power-law hardening material has been carried out for plane strain, Mode 1. The extent to which elasticity affects the near-tip fields is determined by the strain hardening exponent n. Five terms in the asymptotic series for the stresses have been derived for n = 3. However, only three amplitudes can be independently prescribed. These are K1, K2 and K5 corresponding to amplitudes of the first-, second- and fifth-order terms. Four terms in the asymptotic series have been obtained for n = 5, 7 and 10; in these cases, the independent amplitudes are K1, K2 and K4. It is found that appropriate choices of K2 and K4 can reproduce near-tip fields representative of a broad range of crack tip constraints in moderate and low hardening materials. Indeed, fields characterized by distinctly different stress triaxiality levels (established by finite element analysis) have been matched by the asymptotic series. The zone of dominance of the asymptotic series extends over distances of about 10 crack openings ahead of the crack tip encompassing length scales that are microstructurally significant. Furthermore, the higher-order terms collectively describe a spatially uniform hydrostatic stress field (of adjustable magnitude) ahead of the crack. Our results lend support to a suggestion that J and a measure of near-tip stress triaxiality can describe the full range of near-tip states.

General-domain compressible Navier-Stokes solvers exhibiting quasi-unconditional stability and high-order accuracy in space and time

This thesis presents a new class of solvers for the subsonic compressible Navier-Stokes equations in general two- and three-dimensional spatial domains. The proposed methodology incorporates: 1) A novel linear-cost implicit solver based on use of higher-order backward differentiation formulae (BDF) and the alternating direction implicit approach (ADI); 2) A fast explicit solver; 3) Dispersionless spectral spatial discretizations; and 4) A domain decomposition strategy that negotiates the interactions between the implicit and explicit domains. In particular, the implicit methodology is quasi-unconditionally stable (it does not suffer from CFL constraints for adequately resolved flows), and it can deliver orders of time accuracy between two and six in the presence of general boundary conditions. In fact this thesis presents, for the first time in the literature, high-order time-convergence curves for Navier-Stokes solvers based on the ADI strategy---previous ADI solvers for the Navier-Stokes equations have not demonstrated orders of temporal accuracy higher than one. An extended discussion is presented in this thesis which places on a solid theoretical basis the observed quasi-unconditional stability of the methods of orders two through six. The performance of the proposed solvers is favorable. For example, a two-dimensional rough-surface configuration including boundary layer effects at Reynolds number equal to one million and Mach number 0.85 (with a well-resolved boundary layer, run up to a sufficiently long time that single vortices travel the entire spatial extent of the domain, and with spatial mesh sizes near the wall of the order of one hundred-thousandth the length of the domain) was successfully tackled in a relatively short (approximately thirty-hour) single-core run; for such discretizations an explicit solver would require truly prohibitive computing times. As demonstrated via a variety of numerical experiments in two- and three-dimensions, further, the proposed multi-domain parallel implicit-explicit implementations exhibit high-order convergence in space and time, useful stability properties, limited dispersion, and high parallel efficiency.

Spatiotemporal instabilities in nonlinear Kerr media in the presence of arbitrary higher-order dispersions

Spatiotemporal instabilities in nonlinear Kerr media with arbitrary higher-order dispersions are studied by use of standard linear-stability analysis. A generic expression for instability growth rate that unifies and expands on previous results for temporal, spatial, and spatiotemporal instabilities is obtained. It is shown that all odd-order dispersions contribute nothing to instability, whereas all even-order dispersions not only affect the conventional instability regions but may also lead to the appearance of new instability regions. The role of fourth-order dispersion in spatiotemporal instabilities is studied exemplificatively to demonstrate the generic results. Numerical simulations confirm the obtained analytic results. (C) 2002 Optical Society of America.

Higher order vortex gyrotropic modes in circular ferromagnetic nanodots

Magnetic vortex that consists of an in-plane curling magnetization configuration and a needle-like core region with out-of-plane magnetization is known to be the ground state of geometrically confined submicron soft magnetic elements. Here magnetodynamics of relatively thick (50-100 nm) circular Ni80Fe20 dots were probed by broadband ferromagnetic resonance in the absence of external magnetic field. Spin excitation modes related to the thickness dependent vortex core gyrotropic dynamics were detected experimentally in the gigahertz frequency range. Both analytical theory and micromagnetic simulations revealed that these exchange dominated modes are flexure oscillations of the vortex core string with n = 0,1,2 nodes along the dot thickness. The intensity of the mode with n = 1 depends significantly on both dot thickness and diameter and in some cases is higher than the one of the uniform mode with n = 0. This opens promising perspectives in the area of spin transfer torque oscillators.

Temporal difference models describe higher-order learning in humans.

The ability to use environmental stimuli to predict impending harm is critical for survival. Such predictions should be available as early as they are reliable. In pavlovian conditioning, chains of successively earlier predictors are studied in terms of higher-order relationships, and have inspired computational theories such as temporal difference learning. However, there is at present no adequate neurobiological account of how this learning occurs. Here, in a functional magnetic resonance imaging (fMRI) study of higher-order aversive conditioning, we describe a key computational strategy that humans use to learn predictions about pain. We show that neural activity in the ventral striatum and the anterior insula displays a marked correspondence to the signals for sequential learning predicted by temporal difference models. This result reveals a flexible aversive learning process ideally suited to the changing and uncertain nature of real-world environments. Taken with existing data on reward learning, our results suggest a critical role for the ventral striatum in integrating complex appetitive and aversive predictions to coordinate behaviour.

An investigation of higher-order multi-objective optimisation for 3D aerodynamic shape design

We investigate the performance of different variants of a suitably tailored Tabu Search optimisation algorithm on a higher-order design problem. We consider four objective func- tions to describe the performance of a compressor stator row, subject to a number of equality and inequality constraints. The same design problem has been previously in- vestigated through single-, bi- and three-objective optimisation studies. However, in this study we explore the capabilities of enhanced variants of our Multi-objective Tabu Search (MOTS) optimisation algorithm in the context of detailed 3D aerodynamic shape design. It is shown that with these enhancements to the local search of the MOTS algorithm we can achieve a rapid exploration of complicated design spaces, but there is a trade-off be- tween speed and the quality of the trade-off surface found. Rapidly explored design spaces reveal the extremes of the objective functions, but the compromise optimum areas are not very well explored. However, there are ways to adapt the behaviour of the optimiser and maintain both a very efficient rate of progress towards the global optimum Pareto front and a healthy number of design configurations lying on the trade-off surface and exploring the compromise optimum regions. These compromise solutions almost always represent the best qualitative balance between the objectives under consideration. Such enhancements to the effectiveness of design space exploration make engineering design optimisation with multiple objectives and robustness criteria ever more practicable and attractive for modern advanced engineering design. Finally, new research questions are addressed that highlight the trade-offs between intelligence in optimisation algorithms and acquisition of qualita- tive information through computational engineering design processes that reveal patterns and relations between design parameters and objective functions, but also speed versus optimum quality. © 2012 AIAA.

Aeroacoustic source measurement methods for characterizing the sound generated by ducted flow devices with higher-order modes.

The International Organization for Standardization (ISO) method 5136 is widely used in industry and academia to determine the sound power radiated into a duct by fans and other flow devices. The method involves placing the device at the center of a long cylindrical duct with anechoic terminations at each end to eliminate reflections. A single off-axis microphone is used on the inlet and outlet sides that can theoretically capture the plane-wave mode amplitudes but this does not provide enough information to fully account for higher-order modes. In this study, the "two-port" source model is formulated to include higher-order modes and applied for the first three modes. This requires six independent surface pressure measurements on each side or "port." The resulting experimental set-up is much shorter than the ISO rig and does not require anechoic terminations. An array of six external loudspeaker sources is used to characterize the passive part of the two-port model and the set-up provides a framework to account for transmission of higher-order modes through a fan. The relative importance of the higher-order modes has been considered and their effect on inaccuracies when using the ISO method to find source sound power has been analyzed.

Finite Element studies on indentation size effect using a higher order strain gradient theory

In this work, a Finite Element implementation of a higher order strain gradient theory (due to Fleck and Hutchinson, 2001) has been used within the framework of large deformation elasto-viscoplasticity to study the indentation of metals with indenters of various geometries. Of particular interest is the indentation size effect (ISE) commonly observed in experiments where the hardness of a range of materials is found to be significantly higher at small depths of indentation but reduce to a lower, constant value at larger depths. That the ISE can be explained by strain gradient plasticity is well known but this work aims to qualitatively compare a gamut of experimental observations on this effect with predictions from a higher order strain gradient theory. Results indicate that many of the experimental observations are qualitatively borne out by our simulations. However, areas exist where conflicting experimental results make assessment of numerical predictions difficult. © 2012 Elsevier Ltd. All rights reserved.

Reduction in microcystin concentrations in large and shallow lakes: Water and sediment-interface contributions

Blooms of cyanobacteria, or blue-greens, are known to produce chemicals, such as microcystins, which can be toxic to aquatic and terrestrial organisms. Although previous studies have examined the fate of microcystins in freshwater lakes, primary elimination pathways and factors affecting degradation and loss have not been fully explained. The goal of the present study was to explore sources of algal toxins and investigate the distribution and biodegradation of microcystins in water and sediment through laboratory and field analyses. Water and sediment samples were collected monthly from several locations in Lake Taihu from February 2005 to January 2006. Samples were analyzed for the presence of microcystin. Water and sediment were also used in laboratory studies to determine microcystin degradation rates by spiking environmental samples with known concentrations of the chemical and observing concentration changes over time. Some water samples were found to efficiently degrade microcystins. Microcystin concentrations dropped faster in water collected immediately above lake sediment (overlying water). Degradation in sediments was higher than in water. Based on spatial distribution analyses of microcystin in Lake Taihu, higher concentrations (relative to water concentrations) of the chemical were found in lake sediments. These data suggest that sediments play a critical role in microcystin degradation in aquatic systems. The relatively low levels of microcystins found in the environment are most likely due to bacterial biodegradation. Sediments play a crucial role as a source (to the water column) of bio-degrading bacteria and as a carbon-rich environment for bacteria to proliferate and metabolize microcystin and other biogenic toxins produced by cyanobacteria. These, and other, data provide important information that may be applied to management strategies for improvement of water quality in lakes, reservoirs and other water bodies. (C) 2007 Elsevier Ltd. All rights reserved.

Polarity dependence of hexagonal inclusions and cubic twins in GaN/GaAs(001) epilayers measured by conventional X-ray pole figure and grazing incident diffraction pole figure

The distribution of mixed phases and its dependence on the polarity of cubic GaN epilayers are investigated by conventional X-ray pole figure and grazing incident diffraction (GID) pole figure. The hexagonal inclusions and cubic twins can be classified into two portions: one is formed with strict crystalline orientations, the other with crystalline misorientations. The former can be measured by conventional pole figures which reveal that the density of lamellate hexagonal grains and cubic twins located on (1 1 1)(Ga) and ((1) over bar (1) over bar1)(Ga) along [1 (1) over bar 0] direction are higher than those on ((1) over bar 1 1), and (1 (1) over bar 1)(N) along [110] direction. However, the low signals from tiny mixed phases with crystalline misorientations, detected by GID pole figures, distribute in a larger phi region near the [1 1 0] and [(1) over bar (1) over bar 0] directions with much weaker intensity, and in a smaller phi region near the [1 (1) over bar 0] and [(1) over bar 1 0] directions with slightly stronger intensity. (C) 2001 Elsevier Science B.V. All rights reserved.

Higher-order Boussinesq-type equations for interfacial waves in a two-fluid system

Interfacial waves propagating along the interface between a three-dimensional two-fluid system with a rigid upper boundary and an uneven bottom are considered. There is a light fluid layer overlying a heavier one in the system, and a small density difference exists between the two layers. A set of higher-order Boussinesq-type equations in terms of the depth-averaged velocities accounting for stronger nonlinearity are derived. When the small parameter measuring frequency dispersion keeping up to lower-order and full nonlinearity are considered, the equations include the Choi and Camassa's results (1999). The enhanced equations in terms of the depth-averaged velocities are obtained by applying the enhancement technique introduced by Madsen et al. (1991) and Schaffer and Madsen (1995a). It is noted that the equations derived from the present study include, as special cases, those obtained by Madsen and Schaffer (1998). By comparison with the dispersion relation of the linear Stokes waves, we found that the dispersion relation is more improved than Choi and Camassa's (1999) results, and the applicable scope of water depth is deeper.