Dislocation dynamical approach to force fluctuations in nanoindentation experiments

We develop an approach that combines the power of nonlinear dynamics with the evolution equations for the mobile and immobile dislocation densities and force to explain force fluctuations in nanoindentation experiments. The model includes nucleation, multiplication, and propagation thresholds for mobile dislocations, and other well known dislocation transformation mechanisms. The model predicts all the generic features of nanoindentation such as the Hertzian elastic branch followed by several force drops of decreasing magnitudes, and residual plasticity after unloading. The stress corresponding to the elastic force maximum is close to the yield stress of an ideal solid. The predicted values for all the quantities are close to those reported by experiments. Our model allows us to address the indentation-size effect including the ambiguity in defining the hardness in the force drop dominated regime. At large indentation depths, the hardness remains nearly constant with a marginal decreasing trend.

Modelling of transient heat conduction with diffuse interface methods

We present a survey on different numerical interpolation schemes used for two-phase transient heat conduction problems in the context of interface capturing phase-field methods. Examples are general transport problems in the context of diffuse interface methods with a non-equal heat conductivity in normal and tangential directions to the interface. We extend the tonsorial approach recently published by Nicoli M et al (2011 Phys. Rev. E 84 1-6) to the general three-dimensional (3D) transient evolution equations. Validations for one-dimensional, two-dimensional and 3D transient test cases are provided, and the results are in good agreement with analytical and numerical reference solutions.

REDUCED ORDER MODELLING OF COMBUSTION INSTABILITY IN A BACKWARD FACING STEP COMBUSTOR

This paper develops a fully coupled time domain Reduced Order Modelling (ROM) approach to model unsteady combustion dynamics in a backward facing step combustor The acoustic field equations are projected onto the canonical acoustic eigenmodes of the systems to obtain a coupled system of modal evolution equations. The heat release response of the flame is modelled using the G-equation approach. Vortical velocity fluctuations that arise due to shear layer rollup downstream of the step are modelled using a simplified 1D-advection equation whose phase speed is determined from a linear, local, temporal stability analysis of the shear layer just downstream of the step. The hydrodynamic stability analysis reveals a abrupt change in the value of disturbance phase speed from unity for Re < Re-crit to 0.5 for Re > Re-crit, where Remit for the present geometry was found to be approximate to 10425. The results for self-excited flame response show highly wrinkled flame shapes that are qualitatively similar to those seen in prior experiments of acoustically forced flames. The effect of constructive and destructive interference between the two contributions to flame surface wrinkling results in high amplitude wrinkles for the case when K-c -> 1.

Simulations of impinging droplets with surfactant-dependent dynamic contact angle

An arbitrary Lagrangian-Eulerian (ALE) finite element scheme for computations of soluble surfactant droplet impingement on a horizontal surface is presented. The numerical scheme solves the time-dependent Navier-Stokes equations for the fluid flow, scalar convection-diffusion equation for the surfactant transport in the bulk phase, and simultaneously, surface evolution equations for the surfactants on the free surface and on the liquid-solid interface. The effects of surfactants on the flow dynamics are included into the model through the surface tension and surfactant-dependent dynamic contact angle. In particular, the dynamic contact angle (theta(d)) of the droplet is defined as a function of the surfactant concentration at the contact line and the equilibrium contact angle (theta(0)(e)) of the clean surface using the nonlinear equation of state for surface tension. Further, the surface forces are included into the model as surface divergence of the surface stress tensor that allows to incorporate the Marangoni effects without calculating the surface gradient of the surfactant concentration on the free surface. In addition to a mesh convergence study and validation of the numerical results with experiments, the effects of adsorption and desorption surfactant coefficients on the flow dynamics in wetting, partially wetting and non-wetting droplets are studied in detail. It is observed that the effects of surfactants are more in wetting droplets than in the non-wetting droplets. Further, the presence of surfactants at the contact line reduces the equilibrium contact angle further when theta(0)(e) is less than 90 degrees, and increases it further when theta(0)(e) is greater than 90 degrees. Nevertheless, the presence of surfactants has no effect on the contact angle when theta(0)(e) = 90 degrees. The numerical study clearly demonstrates that the surfactant-dependent contact angle has to be considered, in addition to the Marangoni effect, in order to study the flow dynamics and the equilibrium states of surfactant droplet impingement accurately. The proposed numerical scheme guarantees the conservation of fluid mass and of the surfactant mass accurately. (C) 2015 Elsevier Inc. All rights reserved.

General framework for acoustic emission during plastic deformation

Despite the long history, so far there is no general theoretical framework for calculating the acoustic emission spectrum accompanying any plastic deformation. We set up a discrete wave equation with plastic strain rate as a source term and include the Rayleigh-dissipation function to represent dissipation accompanying acoustic emission. We devise a method of bridging the widely separated time scales of plastic deformation and elastic degrees of freedom. While this equation is applicable to any type of plastic deformation, it should be supplemented by evolution equations for the dislocation microstructure for calculating the plastic strain rate. The efficacy of the framework is illustrated by considering three distinct cases of plastic deformation. The first one is the acoustic emission during a typical continuous yield exhibiting a smooth stress-strain curve. We first construct an appropriate set of evolution equations for two types of dislocation densities and then show that the shape of the model stress-strain curve and accompanying acoustic emission spectrum match very well with experimental results. The second and the third are the more complex cases of the Portevin-Le Chatelier bands and the Luders band. These two cases are dealt with in the context of the Ananthakrishna model since the model predicts the three types of the Portevin-Le Chatelier bands and also Luders-like bands. Our results show that for the type-C bands where the serration amplitude is large, the acoustic emission spectrum consists of well-separated bursts of acoustic emission. At higher strain rates of hopping type-B bands, the burst-type acoustic emission spectrum tends to overlap, forming a nearly continuous background with some sharp acoustic emission bursts. The latter can be identified with the nucleation of new bands. The acoustic emission spectrum associated with the continuously propagating type-A band is continuous. These predictions are consistent with experimental results. More importantly, our study shows that the low-amplitude continuous acoustic emission spectrum seen in both the type-B and type-A band regimes is directly correlated to small-amplitude serrations induced by propagating bands. The acoustic emission spectrum of the Luders-like band matches with recent experiments as well. In all of these cases, acoustic emission signals are burstlike, reflecting the intermittent character of dislocation-mediated plastic flow.

Bimodality and re-entrant behaviour in the hierarchical self-assembly of polymeric nanoparticles

We show that a film of a suspension of polymer grafted nanoparticles on a liquid substrate can be employed to create two-dimensional nanostructures with a remarkable variation in the pattern length scales. The presented experiments also reveal the emergence of concentration-dependent bimodal patterns as well as re-entrant behaviour that involves length scales due to dewetting and compositional instabilities. The experimental observations are explained through a gradient dynamics model consisting of coupled evolution equations for the height of the suspension film and the concentration of polymer. Using a Flory-Huggins free energy functional for the polymer solution, we show in a linear stability analysis that the thin film undergoes dewetting and/or compositional instabilities depending on the concentration of the polymer in the solution. We argue that the formation via `hierarchical self-assembly' of various functional nanostructures observed in different systems can be explained as resulting from such an interplay of instabilities.

EXACT INTERNAL CONTROLLABILITY FOR THE WAVE EQUATION IN A DOMAIN WITH OSCILLATING BOUNDARY WITH NEUMANN BOUNDARY CONDITION

In this paper, we study the exact controllability of a second order linear evolution equation in a domain with highly oscillating boundary with homogeneous Neumann boundary condition on the oscillating part of boundary. Our aim is to obtain the exact controllability for the homogenized equation. The limit problem with Neumann condition on the oscillating boundary is different and hence we need to study the exact controllability of this new type of problem. In the process of homogenization, we also study the asymptotic analysis of evolution equation in two setups, namely solution by standard weak formulation and solution by transposition method.

Dislocation dynamics of strain relaxation in epitaxial layers

Many experimental observations have clearly shown that dislocation interaction plays a crucial role in the kinetics of strain relaxation in epitaxial thin films. A set of evolution equations are presented in this article. The key feature of the equations

High Order Multi-Moment Constrained Finite Volume Method. Part I: Basic Formulation

A new high-order finite volume method based on local reconstruction is presented in this paper. The method, so-called the multi-moment constrained finite volume (MCV) method, uses the point values defined within single cell at equally spaced points as the model variables (or unknowns). The time evolution equations used to update the unknowns are derived from a set of constraint conditions imposed on multi kinds of moments, i.e. the cell-averaged value and the point-wise value of the state variable and its derivatives. The finite volume constraint on the cell-average guarantees the numerical conservativeness of the method. Most constraint conditions are imposed on the cell boundaries, where the numerical flux and its derivatives are solved as general Riemann problems. A multi-moment constrained Lagrange interpolation reconstruction for the demanded order of accuracy is constructed over single cell and converts the evolution equations of the moments to those of the unknowns. The presented method provides a general framework to construct efficient schemes of high orders. The basic formulations for hyperbolic conservation laws in 1- and 2D structured grids are detailed with the numerical results of widely used benchmark tests. (C) 2009 Elsevier Inc. All rights reserved.

Invariants, Lie-Bäcklund operators and Bäcklund transformations

This thesis is mainly concerned with the application of groups of transformations to differential equations and in particular with the connection between the group structure of a given equation and the existence of exact solutions and conservation laws. In this respect the Lie-Bäcklund groups of tangent transformations, particular cases of which are the Lie tangent and the Lie point groups, are extensively used.

In Chapter I we first review the classical results of Lie, Bäcklund and Bianchi as well as the more recent ones due mainly to Ovsjannikov. We then concentrate on the Lie-Bäcklund groups (or more precisely on the corresponding Lie-Bäcklund operators), as introduced by Ibragimov and Anderson, and prove some lemmas about them which are useful for the following chapters. Finally we introduce the concept of a conditionally admissible operator (as opposed to an admissible one) and show how this can be used to generate exact solutions.

In Chapter II we establish the group nature of all separable solutions and conserved quantities in classical mechanics by analyzing the group structure of the Hamilton-Jacobi equation. It is shown that consideration of only Lie point groups is insufficient. For this purpose a special type of Lie-Bäcklund groups, those equivalent to Lie tangent groups, is used. It is also shown how these generalized groups induce Lie point groups on Hamilton's equations. The generalization of the above results to any first order equation, where the dependent variable does not appear explicitly, is obvious. In the second part of this chapter we investigate admissible operators (or equivalently constants of motion) of the Hamilton-Jacobi equation with polynornial dependence on the momenta. The form of the most general constant of motion linear, quadratic and cubic in the momenta is explicitly found. Emphasis is given to the quadratic case, where the particular case of a fixed (say zero) energy state is also considered; it is shown that in the latter case additional symmetries may appear. Finally, some potentials of physical interest admitting higher symmetries are considered. These include potentials due to two centers and limiting cases thereof. The most general two-center potential admitting a quadratic constant of motion is obtained, as well as the corresponding invariant. Also some new cubic invariants are found.

In Chapter III we first establish the group nature of all separable solutions of any linear, homogeneous equation. We then concentrate on the Schrodinger equation and look for an algorithm which generates a quantum invariant from a classical one. The problem of an isomorphism between functions in classical observables and quantum observables is studied concretely and constructively. For functions at most quadratic in the momenta an isomorphism is possible which agrees with Weyl' s transform and which takes invariants into invariants. It is not possible to extend the isomorphism indefinitely. The requirement that an invariant goes into an invariant may necessitate variants of Weyl' s transform. This is illustrated for the case of cubic invariants. Finally, the case of a specific value of energy is considered; in this case Weyl's transform does not yield an isomorphism even for the quadratic case. However, for this case a correspondence mapping a classical invariant to a quantum orie is explicitly found.

Chapters IV and V are concerned with the general group structure of evolution equations. In Chapter IV we establish a one to one correspondence between admissible Lie-Bäcklund operators of evolution equations (derivable from a variational principle) and conservation laws of these equations. This correspondence takes the form of a simple algorithm.

In Chapter V we first establish the group nature of all Bäcklund transformations (BT) by proving that any solution generated by a BT is invariant under the action of some conditionally admissible operator. We then use an algorithm based on invariance criteria to rederive many known BT and to derive some new ones. Finally, we propose a generalization of BT which, among other advantages, clarifies the connection between the wave-train solution and a BT in the sense that, a BT may be thought of as a variation of parameters of some. special case of the wave-train solution (usually the solitary wave one). Some open problems are indicated.

Most of the material of Chapters II and III is contained in [I], [II], [III] and [IV] and the first part of Chapter V in [V].

Métodos espectrais aplicados à relatividade numérica: determinação dos dados iniciais

Neste trabalho aplicamos métodos espectrais para a determinação da configuração inicial de três espaços-tempos contendo buracos negros. Para isto apresentamos primeiro a foliação do espaço-tempo em hipersuperfícies tridimensionais espaciais parametrizadas pela função temporal t. Este processo é chamado de decomposição 3+1 [2] [5]. O resultado deste processo são dois conjuntos de equações classificadas em equações de vínculo e evolução [4]. As equações de vínculo podem ser divididas em vínculos Hamiltoniano e dos momentos. Para a obtenção dos dados iniciais dos problemas estudados aqui, apenas a equação de vínculo Hamiltoniano será resolvida numericamente, pois as equações de vínculo dos momentos possuem solução analítica nestes casos. Uma pequena descrição dos métodos espectrais é apresentada, destacando-se os método de Galerkin, método pseudoespectral ou de colocação e método de Tau, que são empregados na resolução das equações de vínculo Hamiltoniano dos problemas estudados. Verificamos que os resultados obtidos neste trabalho superam aqueles produzidos por Kidder e Finn [15], devido a uma escolha diferente das funções de base, que aqui satisfazem uma das condições de contorno.

Modeling and optimizing of high-concentration erbium-doped fiber amplifiers with consideration of ion-clusters

Starting from the modeling of isolated ions and ion-clusters, a closed form rate and power evolution equations for high-concentration erbium-doped fiber amplifiers are constructed. Based on the equations, the effects of the fraction of ion-clusters in total ions and the number of ions per cluster on the performance of high-concentration erbium-doped fiber amplifiers are analyzed numerically. The results show that the presence of the ion-clusters deteriorates amplifier performance, such as the signal power, signal gain, the threshold pump power for zero gain, saturated signal gain, and the maximum gain efficiency, etc. The optimum fiber length or other parameters should be modified with the ion-clusters being taken into account for the amplifiers to achieve a better performance. (c) 2007 Elsevier B.V. All rights reserved.

A global shallow water model using high order multi-moment constrained finite volume method and icosahedral grid

A novel accurate numerical model for shallow water equations on sphere have been developed by implementing the high order multi-moment constrained finite volume (MCV) method on the icosahedral geodesic grid. High order reconstructions are conducted cell-wisely by making use of the point values as the unknowns distributed within each triangular cell element. The time evolution equations to update the unknowns are derived from a set of constrained conditions for two types of moments, i.e. the point values on the cell boundary edges and the cell-integrated average. The numerical conservation is rigorously guaranteed. in the present model, all unknowns or computational variables are point values and no numerical quadrature is involved, which particularly benefits the computational accuracy and efficiency in handling the spherical geometry, such as coordinate transformation and curved surface. Numerical formulations of third and fourth order accuracy are presented in detail. The proposed numerical model has been validated by widely used benchmark tests and competitive results are obtained. The present numerical framework provides a promising and practical base for further development of atmospheric and oceanic general circulation models. (C) 2009 Elsevier Inc. All rights reserved.

油田高含水期剩余油分布规律研究

After systemic investigation of the techniques，route lines and mechanisms about the remaining oil，the dynamic migration and congregation behavior of the remaining oil are discussed on base of interaction between flowing and enriching of water and oil.After the micro-scope modeling of the fluid flow in porous media and the changes in petrol-physical properties of the flowing system, the characters of fluid fields and the dynamic distribution of oil are discussed, among which the preference-flowing is focused on. Based on the preference-flowing in porous media, the concept of the preference-flowing channels is developed. According to above, heterogeneous distribution of water and oil in the field and dynamic mechanism of remaining oil are all obvious. media can be divided into three kinds, directional, stochastic, arbitrary porous media. The main research results are as following: 1. Treating the characteristic parameters such as permeability, porosity and wettability as regional parameter, the fluid field with high water-cut has been established by geostatistical method, among which the difference of flowing pores and the changes of its petrol-physical properties during flooding are studied. 2. The flow process of water and oil are recurrent in physical simulation experiments, in which the mechanisms and phenomena are caught and analyzed. Fluid flow mechanics in porous media with preference-flowing channels have been studied. 3. The mutual coupling between water and oil is induced and the mathematical evolution equations including this interaction were built. . 4. Through coupling effect between flowing water and oil, the dynamic migration and congregation behavior of remaining oil depend upon this coupling. 5. Coupling between water and oil act as driving force and trapping force for the remaining oil. The coupling model of thesis has been verified by simplified the numerical model and compared results with Ng35 oil reservoir in Gudao oil field, it has important theoretical and application values for improving precision of remaining oil and production performance prediction, and is a new method for studying the mechanics of remaining oil in channeled porous media has been established. Key words：flow field，high water-cut，coupling，dominant flow in porous media，remaining oil

Asymptotics of Perturbations to the Wave Equation

M. Hieber, I. Wood: Asymptotics of perturbations to the wave equation. In: Evolution Equations, Lecture Notes in Pure and Appl. Math., 234, Marcel Dekker, (2003), 243-252.