Parallel Schwarz alternating methods for anisotropic diffusion of speckled medical images

This paper deals with ultrasound medical image processing, particularly to filter the noise while preserving the edges and structures of information. The mathematical processing consists in solving by a numerical way a nonlinear evolutive boundary value problem. Several numerical semi-implicit time marching schemes are considered and analyzed. At each time step, parallel synchronous or asynchronous Schwarz alternating methods are used to solve the linear system and its convergence is studied. Lastly, the results of sequential and parallel simulations are presented.      

Resonant Leading Term Geometric Optics Expansions with Boundary Layers for Quasilinear Hyperbolic Boundary Problems

We construct and justify leading order weakly nonlinear geometric optics expansions for nonlinear hyperbolic initial value problems, including the compressible Euler equations. The technique of simultaneous Picard iteration is employed to show approximate solutions tend to the exact solutions in the small wavelength limit. Recent work [2 Coulombel, J.-F., Gues, O., and Williams, M., 2011. Resonant leading order geometric optics expansions for quasilinear hyperbolic fixed and free boundary problems, Comm. Part. Diff. Eqs. 36 (2011), pp. 1797–1859.[Taylor &; Francis Online], [Web of Science ®] , [Google Scholar]] by Coulombel et al. studied the case of reflecting wave trains whose expansions involve only real phases. We treat generic boundary frequencies by incorporating into our expansions both real and nonreal phases. Nonreal phases introduce difficulties such as approximately solving complex transport equations and result in the addition of boundary layers with exponential decay. This also prevents us from doing an error analysis based on almost periodic profiles as in [2 Coulombel, J.-F., Gues, O., and Williams, M., 2011. Resonant leading order geometric optics expansions for quasilinear hyperbolic fixed and free boundary problems, Comm. Part. Diff. Eqs. 36 (2011), pp. 1797–1859.[Taylor &; Francis Online], [Web of Science ®] , [Google Scholar]].     

Nonlinear Schrödinger equation in a semi-strip: Evolution of the Weyl–Titchmarsh function and recovery of the initial condition and rectangular matrix solutions from the boundary conditions

Rectangular matrix solutions of the defocusing nonlinear Schrödinger equation (dNLS) are studied in quarter-plane and semi-strip. Evolution of the corresponding Weyl–Titchmarsh (Weyl) function is described in terms of the initial Weyl function and boundary conditions. In the next step, the initial Weyl function is recovered (for the quarter-plane case) from the long-time asymptotics of the wave function considered at the boundary. Thus, it is shown that the evolution of the Weyl function is uniquely defined by the boundary conditions. Moreover, a procedure to recover solutions of dNLS (uniquely defined by the boundary conditions) is given. In a somewhat different way, the same boundary value problem is also dealt with in a semi-strip (for the case of a quasi-analytic initial condition).     

Uniform stabilization to equilibrium of a nonlinear fluid–structure interaction model

We consider uniform stability to a nontrivial equilibrium of a nonlinear fluid–structure interaction (FSI) defined on a two or three dimensional bounded domain. Stabilization is achieved via boundary and/or interior feedback controls implemented on both the fluid and the structure. The interior damping on the fluid combining with the viscosity effect stabilizes the dynamics of fluid. However, this dissipation propagated from the fluid alone is not sufficient to drive uniformly to equilibrium the entire coupled system. Therefore, additional interior damping on the wave component or boundary porous like damping on the interface is considered. A geometric condition on the interface is needed if only boundary damping on the wave is active. The main technical difficulty is the mismatch of regularity of hyperbolic and parabolic component of the coupled system. This is overcome by considering special multipliers constructed from Stokes solvers. The uniform stabilization result obtained in this article is global for the fully coupled FSI model.     

Aeolian transport of sand

The airborne transport of particles on a granular surface by the saltation mechanism is studied through numerical simulation of particles dragged by turbulent air flow. We calculate the saturated flux q_s and show that its dependence on the wind strength u_* is consistent with several empirical relations obtained from experimental measurements. We propose and explain a new relation for fluxes close to the threshold velocity u_t, namely, q_s=a(u_*-u_t)^α with α≈2. We also obtain the distortion of the velocity profile of the wind due to the drag of the particles and find a novel dynamical scaling relation. We also obtain a new expression for the dependence of the height of the saltation layer as function of the strength of the wind.     

Phenomenology and theory of horizontally oscillated granular mixtures

We overview the physics of a granular mixture subject to horizontal oscillations, recently investigated via experiments and molecular dynamics simulations. First we discuss the rich phenomenology exhibited by this system, which encompasses both segregation and dynamical instabilities. Then we show that the phenomenology can be explained via an effective interaction approach, by which the driven, non-thermal, granular mixture in mapped into a monodispersed thermal system of particles interacting via an effective potential. After determining the effective interaction we discuss its microscopic origin and investigate how it induces the observed phenomenology. Finally, as much as in thermal fluids, from the effective interaction we derive a Cahn-Hilliard dynamics equation, which appears to capture the essential characteristics of the dynamics of the granular mixture.     

Free vibration of non-uniform nanobeams using Rayleigh–Ritz method

In this article, boundary characteristic orthogonal polynomials have been implemented in the Rayleigh–Ritz method to investigate free vibration of non-uniform Euler–Bernoulli nanobeams based on nonlocal elasticity theory. Non-uniform cross section of nanobeams has been considered by taking linear as well as quadratic variations of Young's modulus and density along the space coordinate. Detailed analysis has been reported for all the possible cases of such variations. The objective of the present study is to analyze the effects of nonlocal parameter, boundary condition, length-to-diameter ratio and non-uniform parameter on the frequency parameters. It is found that clamped nanobeams are having highest frequency parameters than other types of boundary conditions for a particular set of parameters. It is also observed that frequency parameters decrease with increase in scaling effect parameter. First four deflection shapes of non-uniform nanobeams have also been incorporated. In this analysis, some of the new results in terms of boundary conditions have also been included.     

One-dimensional migration of interstitial clusters in SUS316L and its model alloys at elevated temperatures

For self-interstitial atom (SIA) clusters in various concentrated alloys, one-dimensional (1D) migration is induced by electron irradiation around 300 K. But at elevated temperatures, the 1D migration frequency decreases to less than one-tenth of that around 300 K in iron-based bcc alloys. In this study, we examined mechanisms of 1D migration at elevated temperatures using in situ observation of SUS316L and its model alloys with high-voltage electron microscopy. First, for elevated temperatures, we examined the effects of annealing and short-term electron irradiation of SIA clusters on their subsequent 1D migration. In annealed SUS316L, 1D migration was suppressed and then recovered by prolonged irradiation at 300 K. In high-purity model alloy Fe-18Cr-13Ni, annealing or irradiation had no effect. Addition of carbon or oxygen to the model alloy suppressed 1D migration after annealing. Manganese and silicon did not suppress 1D migration after annealing but after short-term electron irradiation. The suppression was attributable to the pinning of SIA clusters by segregated solute elements, and the recovery was to the dissolution of the segregation by interatomic mixing under electron irradiation. Next, we examined 1D migration of SIA clusters in SUS316L under continuous electron irradiation at elevated temperatures. The 1D migration frequency at 673 K was proportional to the irradiation intensity. It was as high as half of that at 300 K. We proposed that 1D migration is controlled by the competition of two effects: induction of 1D migration by interatomic mixing and suppression by solute segregation.     

A regularised boundary element formulation for contactless SAR evaluations within homogeneous and inhomogeneous head phantoms

This work presents a Boundary Element Method (BEM) formulation for contactless electromagnetic field assessments. The new scheme is based on a regularised BEM approach that requires the use of electric measurements only. The regularisation is obtained by leveraging on an extension of Calderón techniques to rectangular systems leading to well-conditioned problems independent of the discretisation density. This enables the use of highly discretized Huygens surfaces that can be consequently placed very near to the radiating source. In addition, the new regularised scheme is hybridised with both surfacic homogeneous and volumetric inhomogeneous forward BEM solvers accelerated with fast matrix–vector multiplication schemes. This allows for rapid and effective dosimetric assessments and permits the use of inhomogeneous and realistic head phantoms. Numerical results corroborate the theory and confirms the practical effectiveness of all newly proposed formulations.