Evolution of anodic stress corrosion cracking in a coated material

In the present paper, we investigate the influence of corrosion driving forces and interfacial toughness for a coated material subjected to mechanical loading. If the protective coating is cracked, the substrate material may become exposed to a corrosive media. For a stress corrosion sensitive substrate material, this may lead to detrimental crack growth. A crack is assumed to grow by anodic dissolution, inherently leading to a blunt crack tip. The evolution of the crack surface is modelled as a moving boundary problem using an adaptive finite element method. The rate of dissolution along the crack surface in the substrate is assumed to be proportional to the chemical potential, which is function of the local surface energy density and elastic strain energy density. The surface energy tends to flatten the surface, whereas the strain energy due to stress concentration promotes material dissolution. The influence of the interface energy density parameter for the solid–fluid combination, interface corrosion resistance and stiffness ratios between coating and substrate is investigated. Three characteristic crack shapes are obtained; deepening and narrowing single cracks, branched cracks and sharp interface cracks. The crack shapes obtained by our simulations are similar to real sub-coating cracks reported in the literature

In search of patterns of land-use in Mexico City using logistic regression at the plot level

The study of big cities’ tendency to decentralisation is in the current agenda to understand the structure of Latin American cities. In general, centres and subcentres are related to specific functions. According to the theories of the movement economy and centrality as a process, the urban grid shapes land use distribution through movement and therefore is the main determinant of the location of ‘live centres’, a key component of centres. Activities related to ‘live centres’ include retail, catering and other movement dependent uses. However, the distribution of this kind of activity in cities like Mexico is not as spatially clear as it is in organically grown cities. In this paper we show that, nonetheless complex, there is a relationship between the location of ‘live centre’ uses and spatial configuration. We use multiple logistic regression to evaluate exactly how much influence each variable has on the outcome ‘shop’ given the presence of all the others. The results also suggest different spatial influences for different types of retail on different scales of centres

Nanovoid nucleation by vacancy aggregation and vacancy-cluster coarsening in high-purity metallic single crystals

A numerical model to estimate critical times required for nanovoid nucleation in high-purity aluminum single crystals subjected to shock loading is presented. We regard a nanovoid to be nucleated when it attains a size sufficient for subsequent growth by dislocation-mediated plasticity. Nucleation is assumed to proceed by means of diffusion-mediated vacancy aggregation and subsequent vacancy cluster coarsening. Nucleation times are computed by a combination of lattice kinetic Monte Carlo simulations and simple estimates of nanovoid cavitation pressures and vacancy concentrations. The domain of validity of the model is established by considering rate-limiting physical processes and theoretical strength limits. The computed nucleation times are compared to experiments suggesting that vacancy aggregation and cluster coarsening are feasible mechanisms of nanovoid nucleation in a specific subdomain of the pressure-strain rate-temperature space

Symplectic-energy-momentum preserving variational integrators

The purpose of this paper is to develop variational integrators for conservative mechanical systems that are symplectic and energy and momentum conserving. To do this, a space–time view of variational integrators is employed and time step adaptation is used to impose the constraint of conservation of energy. Criteria for the solvability of the time steps and some numerical examples are given

Dislocation subgrain structures and modeling the plastic hardening of metallic single crystals

A single crystal plasticity theory for insertion into finite element simulation is formulated using sequential laminates to model subgrain dislocation structures. It is known that local models do not adequately account for latent hardening, as latent hardening is not only a material property, but a nonlocal property (e.g. grain size and shape). The addition of the nonlocal energy from the formation of subgrain structure dislocation walls and the boundary layer misfits provide both latent and self-hardening of a crystal slip. Latent hardening occurs as the formation of new dislocation walls limits motion of new mobile dislocations, thus hardening future slip systems. Self-hardening is accomplished by an evolution of the subgrain structure length scale. The substructure length scale is computed by minimizing the nonlocal energy. The minimization of the nonlocal energy is a competition between the dislocation wall energy and the boundary layer energies. The nonlocal terms are also directly minimized within the subgrain model as they affect deformation response. The geometrical relationship between the dislocation walls and slip planes affecting the dislocation mean free path is taken into account, giving a first-order approximation to shape effects. A coplanar slip model is developed due to requirements while modeling the subgrain structure. This subgrain structure plasticity model is noteworthy as all material parameters are experimentally determined rather than fit. The model also has an inherit path dependence due to the formation of the subgrain structures. Validation is accomplished by comparison with single crystal tension test results

Frictional Collisions Off Sharp Objects

This work develops robust contact algorithms capable of dealing with multibody nonsmooth contact geometries for which neither normals nor gap functions can be defined. Such situations arise in the early stage of fragmentation when a number of angular fragments undergo complex collision sequences before eventually scattering. Such situations precludes the application of most contact algorithms proposed to date

Exponential Mixing for a Stochastic PDE Driven by Degenerate Noise

We study stochastic partial differential equations of the reaction-diffusion type. We show that, even if the forcing is very degenerate (i.e. has not full rank), one has exponential convergence towards the invariant measure. The convergence takes place in the topology induced by a weighted variation norm and uses a kind of (uniform) Doeblin condition.Comment: 10 pages, 1 figur

Arbitrary Dimensional Majorana Dualities and Network Architectures for Topological Matter

Motivated by the prospect of attaining Majorana modes at the ends of nanowires, we analyze interacting Majorana systems on general networks and lattices in an arbitrary number of dimensions, and derive various universal spin duals. Such general complex Majorana architectures (other than those of simple square or other crystalline arrangements) might be of empirical relevance. As these systems display low-dimensional symmetries, they are candidates for realizing topological quantum order. We prove that (a) these Majorana systems, (b) quantum Ising gauge theories, and (c) transverse-field Ising models with annealed bimodal disorder are all dual to one another on general graphs. As any Dirac fermion (including electronic) operator can be expressed as a linear combination of two Majorana fermion operators, our results further lead to dualities between interacting Dirac fermionic systems. The spin duals allow us to predict the feasibility of various standard transitions as well as spin-glass type behavior in {\it interacting} Majorana fermion or electronic systems. Several new systems that can be simulated by arrays of Majorana wires are further introduced and investigated: (1) the {\it XXZ honeycomb compass} model (intermediate between the classical Ising model on the honeycomb lattice and Kitaev's honeycomb model), (2) a checkerboard lattice realization of the model of Xu and Moore for superconducting $(p+ip)$ arrays, and a (3) compass type two-flavor Hubbard model with both pairing and hopping terms. By the use of dualities, we show that all of these systems lie in the 3D Ising universality class. We discuss how the existence of topological orders and bounds on autocorrelation times can be inferred by the use of symmetries and also propose to engineer {\it quantum simulators} out of these Majorana networks.Comment: v3,19 pages, 18 figures, submitted to Physical Review B. 11 new figures, new section on simulating the Hubbard model with nanowire systems, and two new appendice

Variational integrators, the Newmark scheme, and dissipative systems

Variational methods are a class of symplectic-momentum integrators for ODEs. Using these schemes, it is shown that the classical Newmark algorithm is structure preserving in a non-obvious way, thus explaining the observed numerical behavior. Modifications to variational methods to include forcing and dissipation are also proposed, extending the advantages of structure preserving integrators to non-conservative systems