Techniques to accelerate convergence of stress-controlled molecular dynamics simulations of dislocation motion

Dislocation mobility —the relation between applied stress and dislocation velocity—is an important property to model the mechanical behavior of structural materials. These mobilities reflect the interaction between the dislocation core and the host lattice and, thus, atomistic resolution is required to capture its details. Because the mobility function is multiparametric, its computation is often highly demanding in terms of computational requirements. Optimizing how tractions are applied can be greatly advantageous in accelerating convergence and reducing the overall computational cost of the simulations. In this paper we perform molecular dynamics simulations of ½ 〈1 1 1〉 screw dislocation motion in tungsten using step and linear time functions for applying external stress. We find that linear functions over time scales of the order of 10–20 ps reduce fluctuations and speed up convergence to the steady-state velocity value by up to a factor of two.

Automatic motion compensation of free breathing acquired myocardial perfusion data by using independent component analysis

Images acquired during free breathing using first-pass gadolinium-enhanced myocardial perfusion magnetic resonance imaging (MRI) exhibit a quasiperiodic motion pattern that needs to be compensated for if a further automatic analysis of the perfusion is to be executed. In this work, we present a method to compensate this movement by combining independent component analysis (ICA) and image registration: First, we use ICA and a time?frequency analysis to identify the motion and separate it from the intensity change induced by the contrast agent. Then, synthetic reference images are created by recombining all the independent components but the one related to the motion. Therefore, the resulting image series does not exhibit motion and its images have intensities similar to those of their original counterparts. Motion compensation is then achieved by using a multi-pass image registration procedure. We tested our method on 39 image series acquired from 13 patients, covering the basal, mid and apical areas of the left heart ventricle and consisting of 58 perfusion images each. We validated our method by comparing manually tracked intensity profiles of the myocardial sections to automatically generated ones before and after registration of 13 patient data sets (39 distinct slices). We compared linear, non-linear, and combined ICA based registration approaches and previously published motion compensation schemes. Considering run-time and accuracy, a two-step ICA based motion compensation scheme that first optimizes a translation and then for non-linear transformation performed best and achieves registration of the whole series in 32 ± 12 s on a recent workstation. The proposed scheme improves the Pearsons correlation coefficient between manually and automatically obtained time?intensity curves from .84 ± .19 before registration to .96 ± .06 after registration

Seismic Hazard and Ground Motion Characterization at the Itoiz Dam (Northern Spain).

This paper presents a new hazard-consistent ground motion characterization of the Itoiz dam site, located in Northern Spain. Firstly, we propose a methodology with different approximation levels to the expected ground motion at the dam site. Secondly, we apply this methodology taking into account the particular characteristics of the site and of the dam. Hazard calculations were performed following the Probabilistic Seismic Hazard Assessment method using a logic tree, which accounts for different seismic source zonings and different ground-motion attenuation relationships. The study was done in terms of peak ground acceleration and several spectral accelerations of periods coinciding with the fundamental vibration periods of the dam. In order to estimate these ground motions we consider two different dam conditions: when the dam is empty (T = 0.1 s) and when it is filled with water to its maximum capacity (T = 0.22 s). Additionally, seismic hazard analysis is done for two return periods: 975 years, related to the project earthquake, and 4,975 years, identified with an extreme event. Soil conditions were also taken into account at the site of the dam. Through the proposed methodology we deal with different forms of characterizing ground motion at the study site. In a first step, we obtain the uniform hazard response spectra for the two return periods. In a second step, a disaggregation analysis is done in order to obtain the controlling earthquakes that can affect the dam. Subsequently, we characterize the ground motion at the dam site in terms of specific response spectra for target motions defined by the expected values SA (T) of T = 0.1 and 0.22 s for the return periods of 975 and 4,975 years, respectively. Finally, synthetic acceleration time histories for earthquake events matching the controlling parameters are generated using the discrete wave-number method and subsequently analyzed. Because of the short relative distances between the controlling earthquakes and the dam site we considered finite sources in these computations. We conclude that directivity effects should be taken into account as an important variable in this kind of studies for ground motion characteristics.

Nonlinear description of transversal motion in a laminar boundary layer with streaks

The nonlinear streamwise growth of a spanwise periodic array of steady streaks in a flat plate boundary layer is numerically computed using the well known Reduced Navier- Stokes formulation. It is found that the flow configuration changes substantially when the amplitude of the streaks grows and the nonlinear effects come into play. The transversal motion (in the wall normal-spanwise plane), which is normally not considered, becomes non-negligible in the nonlinear regime, and it strongly distorts the streamwise velocity profiles, which end up being quite different from those predicted by the linear theory. We analyze in detail the resulting flow patterns for the nonlinearly saturated streaks, and compare them with available experimental results.

Technical evaluation of a third generation optical pose tracker for motion analysis and image-guided surgery

Laparoscopic instrument tracking systems are an essential component in image-guided interventions and offer new possibilities to improve and automate objective assessment methods of surgical skills. In this study we present our system design to apply a third generation optical pose tracker (Micron- Tracker®) to laparoscopic practice. A technical evaluation of this design is performed in order to analyze its accuracy in computing the laparoscopic instrument tip position. Results show a stable fluctuation error over the entire analyzed workspace. The relative position errors are 1.776±1.675 mm, 1.817±1.762 mm, 1.854±1.740 mm, 2.455±2.164 mm, 2.545±2.496 mm, 2.764±2.342 mm, 2.512±2.493 mm for distances of 50, 100, 150, 200, 250, 300, and 350 mm, respectively. The accumulated distance error increases with the measured distance. The instrument inclination covered by the system is high, from 90 to 7.5 degrees. The system reports a low positional accuracy for the instrument tip.

Control motion approach of a lower limb orthosis to reduce energy consumption

By analysing the dynamic principles of the human gait, an economic gait‐control analysis is performed, and passive elements are included to increase the energy efficiency in the motion control of active orthoses. Traditional orthoses use position patterns from the clinical gait analyses (CGAs) of healthy people, which are then de‐normalized and adjusted to each user. These orthoses maintain a very rigid gait, and their energy cosT is very high, reducing the autonomy of the user. First, to take advantage of the inherent dynamics of the legs, a state machine pattern with different gains in eachstate is applied to reduce the actuator energy consumption. Next, different passive elements, such as springs and brakes in the joints, are analysed to further reduce energy consumption. After an off‐line parameter optimization and a heuristic improvement with genetic algorithms, a reduction in energy consumption of 16.8% is obtained by applying a state machine control pattern, and a reduction of 18.9% is obtained by using passive elements. Finally, by combining both strategies, a more natural gait is obtained, and energy consumption is reduced by 24.6%compared with a pure CGA pattern.

Self-similar motion of laser fusion plasmas. Absorption in an unbounded plasma

The one-dimensional motion generated in a cold, infinite, uniform plasma of density na by the absorption, in a certain plane, of a linear pulse of energy per unit time and area = 4>0t/r, 0< t< r, is considered, the analysis allows for thermal conduction and viscosity of ions and electrons, their energy exchange, and an electron heat flux limiter The resulting motion is self-similar and governed by a single nondimensional parameter a«(n0 2T/0)2/3 Detailed asymptotic results are obtained for both a < l and a > l , the general behavior of the solution for arbitrary a is discussed The analysis can be extended to the case of a plasma initially occupying a half-space, and throws light on how to optimize the hydrodynamics of laser fusion plasmas Known approximate results corresponding to motion of a plasma submitted to constant irradiation (<()) are recovered in the present work under appropriate limiting processes

Self-similar motion of laser half-space plasmas. I. Deflagration regime

The one-dimensional self-similar motion of an initially cold, half-space plasma of electron density 0,produced by the (anomalous) absorption of a laser pulse of irradiation = (j>0f/T(0< (< T) at the critical density nc(«c/«0=eoKe)213, where k, m, are Boltzmann's constant and the ion mass, and Ke X (electron temperature)5'2 = heat conductivity. If a >e- 4 ' 3 , a deflagration wave separates an isentropic compression with a shock bounding the undisturbed plasma, and an isentropic expansion flow to the vacuum. The structures of these three regions are completely determined; in particular, the Chapman-Jouguet condition is proved and the density behind the deflagration is found. The deflagration-compression thickness ratio is large (small) for a^e- 5 ' 3(a>e- 5 ' 3 ) . The compression to expansion ratio for both energy and thickness is 0(e"2). For Z,- large, a deflagration exists even if a~e~413. Condition a>e~4'3 may be applied to pulses that are not linear.

Self-similar motion of laser half-space plasmas. II. Thermal waves and intermediate regimes

The one-dimensional self-similar motion of an initially cold, half-space plasma of electron density n,produced by the (anomalous) absorption of a laser pulse of irradiation

€~4'3, a qualitative discussion of how plasma behavior changes with a, is given.

The orbital-motion-limited regime of cylindrical Langmuir probes

An asymptotic analysis of electron collection at high bias Fp serves to determine the domain of validity of the orbital-motion-limited regime of cylindrical Langmuir probes, which is basic for the workings of conductive bare tethers. The radius of a wire collecting OML current in an unmagnetized plasma at rest cannot exceed a value, Rmax , which is found to exhibit a minimum as a function of Fp ; atFp values of interest, Rmax is already increasing and is larger than the electron Debye length lDe . The breakdown of the regime relates to conditions far fromthe probe, at electron energies comparable to the ion thermal energy, kTi ; Rmax is found to increase with Ti . It is also found that ~1! the maximumwidth of a thin tape, if used instead of a wire, is 4Rmax ; ~2! the electron thermal gyroradius must be larger than both R and lDe for magnetic effects to be negligible; and ~3! conditions applying to the tether case are such that trapped-orbit effects are negligible.

Cylindrical Langmuir probes beyond the orbital-motion-limited regime

The current I to a cylindrical probe at rest in an unmagnetized plasma, with probe bias highly positive, is determined. The way I lags behind the orbital-motion-limited OMLcurrent, 1 OML R, as the radius R exceeds the maximum radius for the OML regime to hold, is of interest for space-tether applications. The ratio I/I OML is roughly a decreasing function of R/lD R max /lDe , which is independent of bias, with lDe the electron Debye length and Rmax /l De roughly an increasing function of the temperature ratio, Ti /Te. The dependence of current on ion energy is used to discuss the effect of probe motion through the plasma, a case applying to tethers in low orbit.

Macroscopic motion and gravitation in thermodynamics

A series of examples rarely presented to students is discussed to illustrate a property of thermodynamic equilibrium: small parts of a fully isolated system move as if points of a rigid body, so as to minimize the macroscopic (kinetic) energy EM. Most examples lie in the fields of astronomy and astrophysics, EM then including the gravitational energy. The paradoxical behaviour of gravitation, in particular in the extreme case of black holes,is discussed.

Unattended motion detection system based on the mammalian retina

Sensing systems in living bodies offer a large variety of possible different configurations and philosophies able to be emulated in artificial sensing systems. Motion detection is one of the areas where different animals adopt different solutions and, in most of the cases, these solutions reflect a very sophisticated form. One of them, the mammalian visual system, presents several advantages with respect to the artificial ones. The main objective of this paper is to present a system, based on this biological structure, able to detect motion, its sense and its characteristics. The configuration adopted responds to the internal structure of the mammalian retina, where just five types of cells arranged in five layers are able to differentiate a large number of characteristics of the image impinging onto it. Its main advantage is that the detection of these properties is based purely on its hardware. A simple unit, based in a previous optical logic cell employed in optical computing, is the basis for emulating the different behaviors of the biological neurons. No software is present and, in this way, no possible interference from outside affects to the final behavior. This type of structure is able to work, once the internal configuration is implemented, without any further attention. Different possibilities are present in the architecture to be presented: detection of motion, of its direction and intensity. Moreover, some other characteristics, as symmetry may be obtained.