Dissipation mechanism in 3D magnetic reconnection
Fujimoto, Keizo
2011-11-15
Dissipation processes responsible for fast magnetic reconnection in collisionless plasmas are investigated using 3D electromagnetic particle-in-cell simulations. The present study revisits the two simulation runs performed in the previous study (Fujimoto, Phys. Plasmas 16, 042103 (2009)); one with small system size in the current density direction, and the other with larger system size. In the case with small system size, the reconnection processes are almost the same as those in 2D reconnection, while in the other case a kink mode evolves along the current density and deforms the current sheet structure drastically. Although fast reconnection is achieved in both the cases, the dissipation mechanism is very different between them. In the case without kink mode, the electrons transit the electron diffusion region without thermalization, so that the magnetic dissipation is supported by the inertia resistivity alone. On the other hand, in the kinked current sheet, the electrons are not only accelerated in bulk, but they are also partly scattered and thermalized by the kink mode, which results in the anomalous resistivity in addition to the inertia resistivity. It is demonstrated that in 3D reconnection the thickness of the electron current sheet becomes larger than the local electron inertia length, consistent with the theoretical prediction in Fujimoto and Sydora (Phys. Plasmas 16, 112309 (2009)).
NASA Astrophysics Data System (ADS)
Chai, Kil-Byoung; Zhai, Xiang; Bellan, Paul
2015-11-01
The Caltech astrophysical jet experiment provides a highly resolved demonstration of the interaction between single fluid and 2-fluid scales and possibly kinetic scales as well. The jet evolves through the following sequence: (i) a current-carrying MHD-driven plasma jet self-forms, (ii) the jet undergoes a kink instability, (iii) the kink provides the environment for a secondary, Rayleigh-Taylor (RT) instability, (iv) the RT instability erodes the current channel radius to a scale smaller than ion skin depth to cause fast magnetic reconnection, (v) the reconnection emits broadband obliquely-propagating, right-hand circularly polarized whistler waves, and (vi) the reconnection energizes electrons and ions. The observation of the whistler waves confirms that the reconnection is in the Hall MHD regime (i.e., 2-fluid). An energetic extreme ultraviolet burst is observed at the location of reconnection indicating strong, localized electron heating. Spectroscopic measurement shows simultaneous fast ion heating. The analysis shows that electrons are plausibly heated by Ohmic dissipation, and that ions are plausibly heated stochastically, i.e., the guiding center approximation fails, a kinetic effect.
Local Diagnosis of Reconnection in 3D
NASA Astrophysics Data System (ADS)
Scudder, J. D.; Karimabadi, H.; Daughton, W. S.; Roytershteyn, V.
2014-12-01
We demonstrate (I,II) an approach to find reconnection sites in 3D where there is no flux function for guidance, and where local observational signatures for the ``violation of frozen flux'' are under developed, if not non-existent. We use 2D and 3D PIC simulations of asymmetric guide field reconnection to test our observational hierarchy of single spacecraft kinetic diagnostics - all possible with present state of the art instrumentation. The proliferation of turbulent, electron inertial scale layers in the realistic 3D case demonstrates that electron demagnetization, while necessary, is not sufficient to identify reconnection sites. An excellent local, observable, single spacecraft proxy is demonstrated for the size of the theoretical frozen flux violation. Since even frozen flux violations need not imply reconnection is at hand, a new calibrated dimensionless method is used to determine the importance of such violations. This measure is available in 2D and 3D to help differentiate reconnection layers from weaker frozen flux violating layers. We discuss the possibility that this technique can be implemented on MMS. A technique to highlight flow geometries conducive to reconnection in 3D simulations is also suggested, that may also be implementable with the MMS flotilla. We use local analysis with multiple necessary, but theoretically independent electron kinetic conditions to help reduce the probability of misidentification of any given layer as a reconnection site. Since these local conditions are all necessary for the site, but none is known to be sufficient, the multiple tests help to greatly reduce false positive identifications. The selectivity of the results of this approach using PIC simulations of 3D asymmetric guide field reconnection will be shown using varying numbers of simultaneous conditions. Scudder, J.D., H. Karimabadi, W. Daughton and V. Roytershteyn I, II, submitted Phys. Plasma., 2014
Particle Acceleration in 3D Magnetic Reconnection
NASA Astrophysics Data System (ADS)
Dahlin, J.; Drake, J. F.; Swisdak, M.
2015-12-01
Magnetic reconnection is an important driver of energetic particles in phenomena such as magnetospheric storms and solar flares. Using kinetic particle-in-cell (PIC) simulations, we show that the stochastic magnetic field structure which develops during 3D reconnection plays a vital role in particle acceleration and transport. In a 2D system, electrons are trapped in magnetic islands which limits their energy gain. In a 3D system, however, the stochastic magnetic field enables the energetic electrons to access volume-filling acceleration regions and therefore gain energy much more efficiently than in the 2D system. We also examine the relative roles of two important acceleration drivers: parallel electric fields and a Fermi mechanism associated with reflection of charged particles from contracting field lines. We find that parallel electric fields are most important for accelerating low energy particles, whereas Fermi reflection dominates energetic particle production. We also find that proton energization is reduced in the 3D system.
Impulsive reconnection: 3D onset and stagnation in turbulent paradigms
Sears, Jason A; Intrator, Thomas P; Weber, Tom; Lapenta, Giovanni; Lazarian, Alexander
2010-12-14
Reconnection processes are ubiquitous in solar coronal loops, the earth's magnetotail, galactic jets, and laboratory configurations such as spheromaks and Z pinches. It is believed that reconnection dynamics are often closely linked to turbulence. In these phenomena, the bursty onset of reconnection is partly determined by a balance of macroscopic MHD forces. In a turbulent paradigm, it is reasonable to suppose that there exist many individual reconnection sites, each X-line being finite in axial extent and thus intrinsically three-dimensional (3D) in structure. The balance between MHD forces and flux pile-up continuously shifts as mutually tangled flux ropes merge or bounce. The spatial scale and thus the rate of reconnection are therefore intimately related to the turbulence statistics both in space and in time. We study intermittent 3D reconnection along spatially localized X-lines between two or more flux ropes. The threshold of MHD instability which in this case is the kink threshold is varied by modifying the line-tying boundary conditions. For fast inflow speed of approaching ropes, there is merging and magnetic reconnection which is a well known and expected consequence of the 2D coalescence instability. On the other hand, for slower inflow speed the flux ropes bounce. The threshold appears to be the Sweet Parker speed v{sub A}/S{sup 1/2}, where v{sub A} is the Alfven speed and S is the Lundquist number. Computations by collaborators at University of Wisconsin, Madison, Katholieke Universiteit Leuven, and LANL complement the experiment.
Does fast magnetic reconnection exist?
NASA Technical Reports Server (NTRS)
Priest, E. R.; Forbes, T. G.
1992-01-01
The main features of the Priest-Forbes (1986) and Priest-Lee (1990) models of magnetic reconnection in astrophysical plasmas are discussed, and the Priest-Lee model is generalized to include inflow pressure gradients and thus different regimes of reconnection. It is shown that different scaling results can be obtained depending on the boundary conditions. These results are compared to the ones observed in the numerical experiments of Biskamp (1986) and Lee and Fu (1986). It is concluded that numerical experiments with suitably designed boundary conditions are likely to exhibit fast reconnection, and that such reconnection is a common process in astrophysical and space plasmas.
Laboratory Study of Magnetic Reconnection in 3D Geometry Relevant to Magnetopause and Magnetotail
NASA Astrophysics Data System (ADS)
Ren, Y.; Lu, Q.; Ji, H.; Mao, A.; Wang, X.; E, P.; Wang, Z.; Xiao, Q.; Ding, W.; Zheng, J.
2015-12-01
Laboratory Study of Magnetic Reconnection in 3D Geometry Relevant to Magnetopause and Magnetotail Y. Ren1,2, Quaming Lu3, Hantao Ji1,2, Aohua Mao1, Xiaogang Wang1, Peng E1, Zhibin Wang1, Qingmei Xiao1, Weixing Ding4, Jinxing Zheng51 Harbin Institute of Technology, Harbin, China2 Princeton Plasma Physics Laboratory, Princeton University, Princeton, NJ 08543 3University of Science and Technology of China, Hefei, China 4University of California at Los Angeles, Los Angeles, CA, 90095 5ASIPP, Hefei, China A new magnetic reconnection experiment, Harbin reconnection eXperiment (HRX), is currently being designed as a key part of Space Plasma Environment Research Facility (SPERF) at Harbin Institute of Technology in Harbin, China. HRX aims to provide a unique experimental platform for studying reconnections in 3D geometry relevant to magnetopause and magnetotail to address: the role of electron and ion-scale dynamics in the current sheet; particle and energy transfer from magnetosheath to magnetosphere; particle energization/heating mechanisms during magnetic reconnection; 3D effects in fast reconnection, e.g. the role of 3D magnetic null point. HRX employs a unique set of coils to generate the required 3D magnetic geometry and provides a wide range of plasma parameters. Here, important motivating scientific problems are reviewed and the physics design of HRX is presented, including plasma parameters determined from Vlasov scaling law, reconnection scenarios explored using vacuum magnetic field calculations and numerical simulations of HRX using hybrid and MHD codes. Plasma diagnostics plan and engineering design of important coils will also be briefly presented.
Mechanisms for fast flare reconnection
NASA Technical Reports Server (NTRS)
Vanhoven, G.; Deeds, D.; Tachi, T.
1988-01-01
Normal collisional-resistivity mechanisms of magnetic reconnection have the drawback that they are too slow to explain the fast rise of solar flares. Two methods are examined which are proposed for the speed-up of the magnetic tearing instability: the anomalous enhancement of resistivity by the injection of MHD turbulence and the increase of Coulomb resistivity by radiative cooling. The results are described for nonlinear numerical simulations of these processes which show that the first does not provide the claimed effects, while the second yields impressive rates of reconnection, but low saturated energy outputs.
Wave modes facilitating fast magnetic reconnection
NASA Astrophysics Data System (ADS)
Singh, N.
2011-12-01
Whistler and kinetic Alfven waves are often invoked to explain fast magnetic reconnection in collsionless plasmas. But how these wave modes facilitate the reconnection has remained unclear. An important unanswered question deals with the meaning of the wave frequency in the context of magnetic reconnection. New measurement on a fast explosive reconnection event in the Versatile Toroidal Facility (VTF) at MIT provides an interesting example of the meaning of the wave mode and the associated frequency directly related to the time scale of the impulsive reconnection. We examine the measurements in VTF in view of the whistler wave mode, showing that the explosive growth in the reconnection is related to the thinning of the current sheet to a few electron skin depths. We further demonstrate that the fastest measured time scale (~ 3 microseconds) and the largest normalized reconnection rate (~0.35) agree with those predicted from the whistler mode dispersion relation.
3D outflow jets originating from turbulence in the reconnection current layer
NASA Astrophysics Data System (ADS)
Fujimoto, Keizo
2016-07-01
Satellite observations in the Earth's magnetosphere and in solar flares have suggested that the reconnection outflow jets are fully three dimensional, consisting of a series of narrow channels. The jet structure is important in evaluating the energy and flux transport in the reconnection process. Previous theoretical models based on fluid simulations have relied on patchy reconnection where reconnection takes place predominantly in patchy portions of the current layer. The problem of the previous models is that the gross reconnection rate is much smaller than that in the 2D reconnection case. The present study shows a large-scale 3D PIC simulation revealing that the 3D outflow jets are generated through the 3D flux ropes formed in the turbulent electron current layer around the x-line. Reconnection proceeds almost uniformly along the x-line, so that the gross reconnection rate is comparable to that in the 2D reconnection case. The flux ropes and resultant outflow channels have a typical current-aligned scale provided by the wavelength of an electron shear mode that is much larger than the typical kinetic scales. It is found that the structure of the 3D outflow jets obtained in the simulation is consistent with the bursty bulk flow observed in the Earth's magnetotail.
Fast Reconnection of Weak Magnetic Fields
NASA Technical Reports Server (NTRS)
Zweibel, Ellen G.
1998-01-01
Fast magnetic reconnection refers to annihilation or topological rearrangement of magnetic fields on a timescale that is independent (or nearly independent) of the plasma resistivity. The resistivity of astrophysical plasmas is so low that reconnection is of little practical interest unless it is fast. Yet, the theory of fast magnetic reconnection is on uncertain ground, as models must avoid the tendency of magnetic fields to pile up at the reconnection layer, slowing down the flow. In this paper it is shown that these problems can be avoided to some extent if the flow is three dimensional. On the other hand, it is shown that in the limited but important case of incompressible stagnation point flows, every flow will amplify most magnetic fields. Although examples of fast magnetic reconnection abound, a weak, disordered magnetic field embedded in stagnation point flow will in general be amplified, and should eventually modify the flow. These results support recent arguments against the operation of turbulent resistivity in highly conducting fluids.
NASA Astrophysics Data System (ADS)
Shimizu, Tohru; Torii, Hiroyuki; Kondoh, Koji
2016-05-01
The 3D instability of spontaneous fast magnetic reconnection process is studied with magnetohydrodynamic simulations, where 2D model of the spontaneous fast magnetic reconnection process is destabilized in three dimensions. In this 3D instability, the spontaneous fast magnetic reconnection process is intermittently and randomly caused in 3D. In this paper, as a typical event study, a single 3D fast magnetic reconnection process often observed in the 3D instability is studied in detail. As a remarkable feature, it is reported that, when the 3D fast magnetic reconnection process starts, plasma inflows along the magnetic neutral line are observed, which are driven by plasma static pressure gradient along the neutral line. The plasma inflow speed reaches about 15 in the upstream field region. The unmagnetized inflow tends to prevent the 3D reconnection process; nevertheless, the 3D reconnection process is intermittently maintained. Such high-speed plasma inflows along the neutral line may be observed as dawn-dusk flows in space satellite observations of magnetotail's bursty bulk flows.
Fast magnetic reconnection with large guide fields
Stanier, A.; Simakov, Andrei N.; Chacón, L.; Daughton, W.
2015-01-09
Here, we demonstrate using two-fluid simulations that low-βmagnetic reconnection remains fast, regardless of the presence of fast dispersive waves, which have been previously suggested to play a critical role. In order to understand these results, a discrete model is constructed that offers scaling relationships for the reconnection rate and dissipation region (DR) thickness in terms of the upstream magnetic field and DR length. We verify these scalings numerically and show how the DR self-adjusts to process magnetic flux at the same rate that it is supplied to a larger region where two-fluid effects become important. Ultimately, the rate is independent of the DR physics and is in good agreement with kinetic results.
Fast magnetic reconnection with large guide fields
Stanier, A.; Simakov, Andrei N.; Chacón, L.; Daughton, W.
2015-01-09
Here, we demonstrate using two-fluid simulations that low-βmagnetic reconnection remains fast, regardless of the presence of fast dispersive waves, which have been previously suggested to play a critical role. In order to understand these results, a discrete model is constructed that offers scaling relationships for the reconnection rate and dissipation region (DR) thickness in terms of the upstream magnetic field and DR length. We verify these scalings numerically and show how the DR self-adjusts to process magnetic flux at the same rate that it is supplied to a larger region where two-fluid effects become important. Ultimately, the rate is independentmore » of the DR physics and is in good agreement with kinetic results.« less
NASA Astrophysics Data System (ADS)
Lapenta, Giovanni; Goldman, Martin; Newman, David
2015-11-01
Magnetic reconnection is a mechanism to convert magnetic energy to particle energy in the form of heat and directed flows. We study here where reconnection and particle energization are found in full 3D models of a reconnecting plasma sheet. Three regions emerge as the primary loci of energy conversion: the separatrices, the dipolarization fronts and the electron diffusion region near x-points. We consider two scenarios: one where the exhaust from multiple x-lines forms a plasmoid (a flux rope in 3D) and one where the exhaust encounters pristine unreconnected plasma and forms a pile-up front. A key process intrinsically 3D, not present in 2D, is the development of an instability in the outflow leading to the formation of secondary reconnection sites that further enhance energy conversion. The MMS mission of NASA was launched on March 12 of this year with the stated goal of finding these regions. We will soon know if we are right in predicting these additional regions of dissipation in the reconnection outflow. Work supported by the NASA MMS Program and by the BOF funds of the KU Leuven. Simulations conducted on NASA Computing facilities, PRACE Tier-0 computing and at NERSC (DOE Office of Science User Contract No. DE-AC02-05CH11231).
3D instabilities connected with reconnection in full 3D PIC simulations
NASA Astrophysics Data System (ADS)
Lapenta, Giovanni
2013-10-01
Kinetic reconnection is characterized by a distinct behavior of electrons and ions with regions of strong relative speeds between the species. Electrons can flow at great speed relative to ions and can be characterized by a strong non-gyrotropy and anisotropy. When studied in full three dilensions, these electron peculiar properties can drive numerous instabilities that have been investigated by the suggested speaker and his collaborators in a number of recent published papers. Two regions have received most attention: 1) the separatrices where instabilities are caused by the electron flow and the electron phase space features, 2) the downstream fronts where an interchange instability leads to strong energy exchanges and secondary reconnection. In both situations the ions are demagnitezed but the electrons are not and their behaviour is rich in full kinetic processes. At the separatrices, two types of instabilities have been observed. The electron phase space is characterized by multiple populations at relative drifts (electron beams) and the whole electron species is drifting with respect to the ions. This condition is subject to different streaming instabilities. Additionally, the separatrices are regions of intense density and flow shear, with free energy available to drive Kelvin-Helmholtz-type instabilities. In the downstream fronts of reconnection, a density gradient develops in conditions where the acceleration is directed unfavourably for stability, leading to ballooning and interchange-type instabilities. Both cases are of great importance for the upcoming Magnetospheric Multiscale Mission that is bent on finding and analyzing the regions where the electron scale physics is dominant. The processes discussed above can provide key information for the operation of the mission and the interpretation of its results. Collaboration between the University of Colorado NASA-MMSIDS team (M. Goldman, D. Newman, L. Anderson, S. Erikson) and the KULeuven Swiff team
Experimental onset threshold and magnetic pressure pileup for 3D Sweet-Parker reconnection
Intrator, Thomas P; Sun, Xuan; Lapenta, Giovanni; Furno, Ivo
2008-01-01
In space, astrophysical and laboratory plasmas, magnetic reconnect ion converts magnetic into particle energy during unsteady, explosive events. The abrupt onset and cessation has been a long standing puzzle. We show the first three-dimensional (3D) laboratory example of onset and stagnation of Sweet-Parker type magnetic reconnection between magnetized and parallel current (flux) ropes driven by magnetohydrodynamic (MHD) attraction and 3D instability. Mutually attracting flux ropes advect and merge oppositely directed magnetic fields. Magnetic flux is annihilated, but reaches soon a threshold where magnetic flux and pressure pile up, and reconnection magnetic topology appears. This occurs when inflow speeds exceed the SweetParker speed v{sub SP} = v{sub A} / S{sup 1/2}, where v{sub A} is the Alfven speed and S is the Lundquist number for the reconnection layer, as magnetic flux arrives faster than flux annihilation can process it. Finally piled up fields generate MHD reaction forces that stall the inflow and the reconnection process.
FastScript3D - A Companion to Java 3D
NASA Technical Reports Server (NTRS)
Koenig, Patti
2005-01-01
FastScript3D is a computer program, written in the Java 3D(TM) programming language, that establishes an alternative language that helps users who lack expertise in Java 3D to use Java 3D for constructing three-dimensional (3D)-appearing graphics. The FastScript3D language provides a set of simple, intuitive, one-line text-string commands for creating, controlling, and animating 3D models. The first word in a string is the name of a command; the rest of the string contains the data arguments for the command. The commands can also be used as an aid to learning Java 3D. Developers can extend the language by adding custom text-string commands. The commands can define new 3D objects or load representations of 3D objects from files in formats compatible with such other software systems as X3D. The text strings can be easily integrated into other languages. FastScript3D facilitates communication between scripting languages [which enable programming of hyper-text markup language (HTML) documents to interact with users] and Java 3D. The FastScript3D language can be extended and customized on both the scripting side and the Java 3D side.
Magnetic reconnection in 3D magnetosphere models: magnetic separators and open flux production
NASA Astrophysics Data System (ADS)
Glocer, A.; Dorelli, J.; Toth, G.; Komar, C. M.; Cassak, P.
2014-12-01
There are multiple competing definitions of magnetic reconnection in 3D (e.g., Hesse and Schindler [1988], Lau and Finn [1990], and Boozer [2002]). In this work we focus on separator reconnection. A magnetic separator can be understood as the 3D analogue of a 2D x line with a guide field, and is defined by the line corresponding to the intersection of the separatrix surfaces associated with the magnetic nulls. A separator in the magnetosphere represents the intersection of four distinct magnetic topologies: solar wind, closed, open connected to the northern hemisphere, and open connected to the southern hemisphere. The integral of the parallel electric field along the separator defines the rate of open flux production, and is one measure of the reconnection rate. We present three methods for locating magnetic separators and apply them to 3D resistive MHD simulations of the Earth's magnetosphere using the BATS-R-US code. The techniques for finding separators and determining the reconnection rate are insensitive to IMF clock angle and can in principle be applied to any magnetospheric model. The present work examines cases of high and low resistivity, for two clock angles. We also examine the separator during Flux Transfer Events (FTEs) and Kelvin-Helmholtz instability.
NASA Astrophysics Data System (ADS)
Lapenta, Giovanni; Goldman, Martin; Newman, David; olshevskyi, Vyacheslav; Markidis, Stefano
2016-04-01
Dipolarization fronts (DF) are formed by reconnection outflows interacting with the pre-existing environment. These regions are host of important energy exchanges [1], particle acceleration [2] and a complex structure and evolution [3]. Our recent work has investigated these regions via fully kinetic 3D simulations [4]. As reported recently on Nature Physics [3], based on 3D fully kinetic simulations started with a well defined x-line, we observe that in the DF reconnection transitions towards a more chaotic regime. In the fronts an instability devel- ops caused by the local gradients of the density and by the unfavourable acceleration and field line curvature. The consequence is the break up of the fronts in a fashion similar to the classical fluid Rayleigh-Taylor instability with the formation of "fingers" of plasma and embedded magnetic fields. These fingers interact and produce secondary reconnection sites. We present several different diagnostics that prove the existence of these secondary reconnection sites. Each site is surrounded by its own electron diffusion region. At the fronts the ions are generally not magnetized and considerable ion slippage is present. The discovery we present is that electrons are also slipping, forming localized diffusion regions near secondary reconnection sites [1]. The consequence of this discovery is twofold. First, the instability in the fronts has strong energetic implications. We observe that the energy transfer locally is very strong, an order of magnitude stronger than in the "X" line. However, this energy transfer is of both signs as it is natural for a wavy rippling with regions of magnetic to kinetic and regions of kinetic to magnetic energy conversion. Second, and most important for this session, is that MMS should not limit the search for electron diffusion regions to the location marked with X in all reconnection cartoons. Our simulations predict more numerous and perhaps more easily measurable electron diffusion
Catastrophic onset of fast magnetic reconnection with a guide field
NASA Astrophysics Data System (ADS)
Cassak, P. A.; Drake, J. F.; Shay, M. A.
2007-05-01
It was recently shown that the slow (collisional) Sweet-Parker and the fast (collisionless) Hall magnetic reconnection solutions simultaneously exist for a wide range of resistivities; reconnection is bistable [Cassak, Shay, and Drake, Phys. Rev. Lett., 95, 235002 (2005)]. When the thickness of the dissipation region becomes smaller than a critical value, the Sweet-Parker solution disappears and fast reconnection ensues, potentially explaining how large amounts of magnetic free energy can accrue without significant release before the onset of fast reconnection. Two-fluid numerical simulations extending the previous results for anti-parallel reconnection (where the critical thickness is the ion skin depth) to component reconnection with a large guide field (where the critical thickness is the thermal ion Larmor radius) are presented. Applications to laboratory experiments of magnetic reconnection and the sawtooth crash are discussed.
3D palmprint data fast acquisition and recognition
NASA Astrophysics Data System (ADS)
Wang, Xiaoxu; Huang, Shujun; Gao, Nan; Zhang, Zonghua
2014-11-01
This paper presents a fast 3D (Three-Dimension) palmprint capturing system and develops an efficient 3D palmprint feature extraction and recognition method. In order to fast acquire accurate 3D shape and texture of palmprint, a DLP projector triggers a CCD camera to realize synchronization. By generating and projecting green fringe pattern images onto the measured palm surface, 3D palmprint data are calculated from the fringe pattern images. The periodic feature vector can be derived from the calculated 3D palmprint data, so undistorted 3D biometrics is obtained. Using the obtained 3D palmprint data, feature matching test have been carried out by Gabor filter, competition rules and the mean curvature. Experimental results on capturing 3D palmprint show that the proposed acquisition method can fast get 3D shape information of palmprint. Some initial experiments on recognition show the proposed method is efficient by using 3D palmprint data.
NASA Astrophysics Data System (ADS)
Wilder, F. D.; Ergun, R.; Goodrich, K.; Malaspina, D.; Eriksson, S.; Stawarz, J. E.; Sturner, A. P.; Holmes, J.; Burch, J. L.; Torbert, R. B.; Phan, T.; Le Contel, O.; Goldman, M. V.; Newman, D. L.; Lindqvist, P. A.; Khotyaintsev, Y. V.; Strangeway, R. J.; Russell, C. T.; Giles, B. L.; Pollock, C. J.
2015-12-01
The phenomenon of magnetic reconnection, especially at electron scales, is still poorly understood. One process that warrants further investigation is the role of wave phenomenon in mediating magnetic reconnection. Previous observations have shown the presence of electrostatic solitary waves (ESWs) as well as whistler mode waves near the dayside reconnection site. Additionally, recent simulations have suggested that whistler waves might be generated by electron phase space holes associated with ESWs as they propagate along the magnetic separatrix towards the diffusion region. Other observations have shown ESWs with distinct speeds and time scales, suggesting that different instabilities generate the ESWs. NASA's recently launched Magnetospheric Multiscale (MMS) mission presents a unique opportunity to investigate the roles of wave phenomena, such as ESWs and whistlers, in asymmetric reconnection at the dayside magnetopause. We will present 3-D electric and magnetic field data from magnetopause crossings by MMS during its first dayside science phase. Burst mode wave data and electron distributions from all four spacecraft will be analyzed to investigate the origin of these wave phenomena, as well as their impact on the reconnection electric field.
Existence of two MHD reconnection modes in a solar 3D magnetic null point topology
NASA Astrophysics Data System (ADS)
Pariat, Etienne; Antiochos, Spiro; DeVore, C. Richard; Dalmasse, Kévin
2012-07-01
Magnetic topologies with a 3D magnetic null point are common in the solar atmosphere and occur at different spatial scales: such structures can be associated with some solar eruptions, with the so-called pseudo-streamers, and with numerous coronal jets. We have recently developed a series of numerical experiments that model magnetic reconnection in such configurations in order to study and explain the properties of jet-like features. Our model uses our state-of-the-art adaptive-mesh MHD solver ARMS. Energy is injected in the system by line-tied motion of the magnetic field lines in a corona-like configuration. We observe that, in the MHD framework, two reconnection modes eventually appear in the course of the evolution of the system. A very impulsive one, associated with a highly dynamic and fully 3D current sheet, is associated with the energetic generation of a jet. Before and after the generation of the jet, a quasi-steady reconnection mode, more similar to the standard 2D Sweet-Parker model, presents a lower global reconnection rate. We show that the geometry of the magnetic configuration influences the trigger of one or the other mode. We argue that this result carries important implications for the observed link between observational features such as solar jets, solar plumes, and the emission of coronal bright points.
Kinetic turbulence in 3D collisionless magnetic reconnection with a guide magnetic field
NASA Astrophysics Data System (ADS)
Alejandro Munoz Sepulveda, Patricio; Kilian, Patrick; Jain, Neeraj; Büchner, Jörg
2016-04-01
The features of kinetic plasma turbulence developed during non-relativistic 3D collisionless magnetic reconnection are still not fully understood. This is specially true under the influence of a strong magnetic guide field, a scenario common in space plasmas such as in the solar corona and also in laboratory experiments such as MRX or VINETA II. Therefore, we study the mechanisms and micro-instabilities leading to the development of turbulence during 3D magnetic reconnection with a fully kinetic PIC code, emphasizing the role of the guide field with an initial setup suitable for the aforementioned environments. We also clarify the relations between these processes and the generation of non-thermal populations and particle acceleration.
3D meshes of carbon nanotubes guide functional reconnection of segregated spinal explants.
Usmani, Sadaf; Aurand, Emily Rose; Medelin, Manuela; Fabbro, Alessandra; Scaini, Denis; Laishram, Jummi; Rosselli, Federica B; Ansuini, Alessio; Zoccolan, Davide; Scarselli, Manuela; De Crescenzi, Maurizio; Bosi, Susanna; Prato, Maurizio; Ballerini, Laura
2016-07-01
In modern neuroscience, significant progress in developing structural scaffolds integrated with the brain is provided by the increasing use of nanomaterials. We show that a multiwalled carbon nanotube self-standing framework, consisting of a three-dimensional (3D) mesh of interconnected, conductive, pure carbon nanotubes, can guide the formation of neural webs in vitro where the spontaneous regrowth of neurite bundles is molded into a dense random net. This morphology of the fiber regrowth shaped by the 3D structure supports the successful reconnection of segregated spinal cord segments. We further observed in vivo the adaptability of these 3D devices in a healthy physiological environment. Our study shows that 3D artificial scaffolds may drive local rewiring in vitro and hold great potential for the development of future in vivo interfaces. PMID:27453939
3D meshes of carbon nanotubes guide functional reconnection of segregated spinal explants
Usmani, Sadaf; Aurand, Emily Rose; Medelin, Manuela; Fabbro, Alessandra; Scaini, Denis; Laishram, Jummi; Rosselli, Federica B.; Ansuini, Alessio; Zoccolan, Davide; Scarselli, Manuela; De Crescenzi, Maurizio; Bosi, Susanna; Prato, Maurizio; Ballerini, Laura
2016-01-01
In modern neuroscience, significant progress in developing structural scaffolds integrated with the brain is provided by the increasing use of nanomaterials. We show that a multiwalled carbon nanotube self-standing framework, consisting of a three-dimensional (3D) mesh of interconnected, conductive, pure carbon nanotubes, can guide the formation of neural webs in vitro where the spontaneous regrowth of neurite bundles is molded into a dense random net. This morphology of the fiber regrowth shaped by the 3D structure supports the successful reconnection of segregated spinal cord segments. We further observed in vivo the adaptability of these 3D devices in a healthy physiological environment. Our study shows that 3D artificial scaffolds may drive local rewiring in vitro and hold great potential for the development of future in vivo interfaces. PMID:27453939
FAST MAGNETIC RECONNECTION AND SPONTANEOUS STOCHASTICITY
Eyink, Gregory L.; Lazarian, A.; Vishniac, E. T.
2011-12-10
Magnetic field lines in astrophysical plasmas are expected to be frozen-in at scales larger than the ion gyroradius. The rapid reconnection of magnetic-flux structures with dimensions vastly larger than the gyroradius requires a breakdown in the standard Alfven flux-freezing law. We attribute this breakdown to ubiquitous MHD plasma turbulence with power-law scaling ranges of velocity and magnetic energy spectra. Lagrangian particle trajectories in such environments become 'spontaneously stochastic', so that infinitely many magnetic field lines are advected to each point and must be averaged to obtain the resultant magnetic field. The relative distance between initial magnetic field lines which arrive at the same final point depends upon the properties of two-particle turbulent dispersion. We develop predictions based on the phenomenological Goldreich and Sridhar theory of strong MHD turbulence and on weak MHD turbulence theory. We recover the predictions of the Lazarian and Vishniac theory for the reconnection rate of large-scale magnetic structures. Lazarian and Vishniac also invoked 'spontaneous stochasticity', but of the field lines rather than of the Lagrangian trajectories. More recent theories of fast magnetic reconnection appeal to microscopic plasma processes that lead to additional terms in the generalized Ohm's law, such as the collisionless Hall term. We estimate quantitatively the effect of such processes on the inertial-range turbulence dynamics and find them to be negligible in most astrophysical environments. For example, the predictions of the Lazarian and Vishniac theory are unchanged in Hall MHD turbulence with an extended inertial range, whenever the ion skin depth {delta}{sub i} is much smaller than the turbulent integral length or injection-scale L{sub i} .
Physical conditions for fast reconnection evolution in space plasmas
Ugai, M.
2012-07-15
The present paper studies physical conditions for fast reconnection mechanism involving slow shocks to evolve spontaneously in space (high-temperature) plasmas. This is fundamental for onset mechanisms of geomagnetic substorms and solar flares. It is demonstrated that reconnection evolution strongly depends on effective resistivity available in space plasmas as well as on dimensions of initial current sheet. If a current sheet is sufficiently thin, fast reconnection spontaneously evolves only when resistivity is locally enhanced around X reconnection point. This is because in space plasmas reconnection flows cause vital current concentration locally around X point. For current-driven anomalous resistivity, the resulting resistivity is automatically localized around X point, so fast reconnection mechanism can be realized. On the other hand, for uniform or Spitzer resistivity, any fast reconnection cannot grow; in particular, Spitzer resistivity is reduced around X point because of Joule heating. Regarding reconnection simulations (either fluid or particle), unless numerical resistivities are made negligibly small, they seriously mask the effects of physical resistivity, leading to a misleading conclusion that reconnection evolution is little influenced by plasma resistivity.
Measurements of Fast Magnetic Reconnection Driven by Relativistic Electrons
NASA Astrophysics Data System (ADS)
Raymond, Anthony; McKelvey, Andrew; Zulick, Calvin; Chuanfei, Dong; Maksimchuk, Anatoly; Thomas, Alexander; Yanovsky, Victor; Krushelnick, Karl; Willingale, Louise; Chykov, Vladimir; Nilson, Phil; Chen, Hui; Williams, Gerald; Bhattacharjee, Amitava; Fox, Will
2015-11-01
Magnetic reconnection is a process whereby opposing magnetic field lines are forced together and topologically rearrange, resulting in lower magnetic potential energy and in corresponding plasma heating. Such occurrences are ubiquitous in astrophysics as well as appearing in laboratory plasmas such as in ICF in the form of instabilities. We report measurements in the domain of ultra-fast, ultra-intense lasers, in which the mechanism responsible follows from radially expanding surface electrons with v ~ c . Results are compared from two laser facilities (HERCULES and Omega EP), both of which produced two relativistic intensity pulses focused within close proximity onto copper foils. A spherical X-ray crystal was used to image the Kα radiation induced by electron currents, revealing the midplane diffusion region wherein electrons are accelerated into the target by the electric field generated during reconnection. The characteristics of this signal are studied as a function of the focal spot separation, laser energy, and pulse duration. The results are then compared to 3D PIC simulations.
Slip versus Field-Line Mapping in Describing 3D Reconnection of Coronal Magnetic Fields
NASA Astrophysics Data System (ADS)
Titov, V. S.; Mikic, Z.; Torok, T.; Downs, C.; Lionello, R.; Linker, J.
2015-12-01
We demonstrate two techniques for describing the structure of the coronal magnetic field and its evolution due to reconnection in numerical 3D simulations of the solar corona and CMEs. These techniques employ two types of mapping of the boundary of the computational domain on itself. One of them is defined at a given time moment via connections of the magnetic field lines to their opposite endpoints. The other mapping, called slip mapping, relates field line endpoints at two different time moments and allows one to identify the slippage of plasma elements due to resistivity across field lines for a given time interval (Titov et al. 2009). The distortion of each of these mappings can be measured by using the so-called squashing factor Q (Titov 2007). The high-Q layers computed for the first and second mappings define, respectively, (quasi-)separatrix surfaces and reconnection fronts in evolving magnetic configurations. Analyzing these structural features, we are able to reveal topologically different domains and reconnected flux systems in the configurations, in particular, open, closed and disconnected magnetic flux tubes, as well as quantify the related magnetic flux transfer. Comparison with observations makes it possible also to relate these features to observed morphological elements such as flare loops and ribbons, and EUV dimmings. We illustrate these general techniques by applying them to particular data-driven MHD simulations. *Research supported by NASA's HSR and LWS Programs, and NSF/SHINE and NSF/FESD.
ONSET OF FAST MAGNETIC RECONNECTION IN PARTIALLY IONIZED GASES
Malyshkin, Leonid M.; Zweibel, Ellen G. E-mail: zweibel@astro.wisc.edu
2011-10-01
We consider quasi-stationary two-dimensional magnetic reconnection in a partially ionized incompressible plasma. We find that when the plasma is weakly ionized and the collisions between the ions and the neutral particles are significant, the transition to fast collisionless reconnection due to the Hall effect in the generalized Ohm's law is expected to occur at much lower values of the Lundquist number, as compared to a fully ionized plasma case. We estimate that these conditions for fast reconnection are satisfied in molecular clouds and in protostellar disks.
Fast Magnetic Reconnection: Bridging Laboratory and Space Plasma Physics
Bhattacharjee, Amitava
2012-02-16
Recent developments in experimental and theoretical studies of magnetic reconnection hold promise for providing solutions to outstanding problems in laboratory and space plasma physics. Examples include sawtooth crashes in tokamaks, substorms in the Earth’s Magnetosphere, eruptive solar flares, and more recently, fast reconnection in laser-produced high energy density plasmas. In each of these examples, a common and long-standing challenge has been to explain why fast reconnection proceeds rapidly from a relatively quiescent state. In this talk, we demonstrate the advantages of viewing these problems and their solutions from a common perspective. We focus on some recent, surprising discoveries regarding the role of secondary plasmoid instabilities of thin current sheets. Nonlinearly, these instabilities lead to fast reconnection rates that are very weakly dependent on the Lundquist number of the plasma.
NASA Astrophysics Data System (ADS)
Shuster, J. R.; Torbert, R. B.; Vaith, H.; Argall, M. R.; Li, G.; Chen, L. J.; Ergun, R. E.; Lindqvist, P. A.; Marklund, G. T.; Khotyaintsev, Y. V.; Russell, C. T.; Magnes, W.; Le Contel, O.; Pollock, C. J.; Giles, B. L.
2015-12-01
The electron drift instruments (EDIs) onboard each MMS spacecraft are designed with large geometric factors (~0.01cm2 str) to facilitate detection of weak (~100 nA) electron beams fired and received by the two gun-detector units (GDUs) when EDI is in its "electric field mode" to determine the local electric and magnetic fields. A consequence of the large geometric factor is that "ambient mode" electron flux measurements (500 eV electrons having 0°, 90°, or 180° pitch angle) can vary depending on the orientation of the EDI instrument with respect to the magnetic field, a nonphysical effect that requires a correction. Here, we present determinations of the θ- and ø-dependent correction factors for the eight EDI GDUs, where θ (ø) is the polar (azimuthal) angle between the GDU symmetry axis and the local magnetic field direction, and compare the corrected fluxes with those measured by the fast plasma instrument (FPI). Using these corrected, high time resolution (~1,000 samples per second) ambient electron fluxes, combined with the unprecedentedly high resolution 3D electric field measurements taken by the spin-plane and axial double probes (SDP and ADP), we are equipped to accurately detect electron-scale current layers and electric field waves associated with the non-Maxwellian (anisotropic and agyrotropic) particle distribution functions predicted to exist in the reconnection diffusion region. We compare initial observations of the diffusion region with distributions and wave analysis from PIC simulations of asymmetric reconnection applicable for modeling reconnection at the Earth's magnetopause, where MMS will begin Science Phase 1 as of September 1, 2015.
Fast magnetic reconnection in three-dimensional magnetohydrodynamics simulations
Pang Bijia; Pen, U.-L.; Vishniac, Ethan T.
2010-10-15
A constructive numerical example of fast magnetic reconnection in a three-dimensional periodic box is presented. Reconnection is initiated by a strong, localized perturbation to the field lines. The solution is intrinsically three-dimensional and its gross properties do not depend on the details of the simulations. {approx}30% of the magnetic energy is released in an event which lasts about one Alfven time, but only after a delay during which the field lines evolve into a critical configuration. The physical picture of the process is presented. The reconnection regions are dynamical and mutually interacting. In the comoving frame of these regions, reconnection occurs through a x-like point, analogous to Petschek reconnection. The dynamics appear to be driven by global flows, not local processes.
3D fast wavelet network model-assisted 3D face recognition
NASA Astrophysics Data System (ADS)
Said, Salwa; Jemai, Olfa; Zaied, Mourad; Ben Amar, Chokri
2015-12-01
In last years, the emergence of 3D shape in face recognition is due to its robustness to pose and illumination changes. These attractive benefits are not all the challenges to achieve satisfactory recognition rate. Other challenges such as facial expressions and computing time of matching algorithms remain to be explored. In this context, we propose our 3D face recognition approach using 3D wavelet networks. Our approach contains two stages: learning stage and recognition stage. For the training we propose a novel algorithm based on 3D fast wavelet transform. From 3D coordinates of the face (x,y,z), we proceed to voxelization to get a 3D volume which will be decomposed by 3D fast wavelet transform and modeled after that with a wavelet network, then their associated weights are considered as vector features to represent each training face . For the recognition stage, an unknown identity face is projected on all the training WN to obtain a new vector features after every projection. A similarity score is computed between the old and the obtained vector features. To show the efficiency of our approach, experimental results were performed on all the FRGC v.2 benchmark.
NASA Astrophysics Data System (ADS)
Krebs, Isabel; Jardin, Stephen C.; Igochine, Valentin; Guenter, Sibylle; Hoelzl, Matthias; ASDEX Upgrade Team
2014-10-01
We study sawtooth reconnection in ASDEX Upgrade tokamak plasmas by means of 3D non-linear two-fluid MHD simulations in toroidal geometry using the high-order finite element code M3D-C1. Parameters and equilibrium of the simulations are based on typical sawtoothing ASDEX Upgrade discharges. The simulation results are compared to features of the experimental observations such as the sawtooth crash time and frequency, the evolution of the safety factor profile and the 3D evolution of the temperature. 2D ECE imaging measurements during sawtooth crashes in ASDEX Upgrade indicate that the heat is transported out of the core through a narrow poloidally localized region. We investigate if incomplete sawtooth reconnection can be seen in the simulations which is suggested by soft X-ray tomography measurements in ASDEX Upgrade showing that an (m = 1, n = 1) perturbation is typically observed to survive the sawtooth crash and approximately maintain its radial position.
Fast 3D shape measurements using laser speckle projection
NASA Astrophysics Data System (ADS)
Schaffer, Martin; Grosse, Marcus; Harendt, Bastian; Kowarschik, Richard
2011-05-01
3D measurement setups based on structured light projection are widely used for many industrial applications. Due to intense research in the past the accuracy is comparably high in connection with relatively low cost of the equipment. But facing higher acquisition rates in industries especially for chain assembling lines there are still hurdles to take when accelerating 3D measurements and at the same time retaining accuracies. We developed a projection technique that uses laser speckles to enable fast 3D measurements with statistically structured light patterns. In combination with a temporal correlation technique dense and accurate 3D reconstructions at nearly video rate can be achieved.
Increased Speed: 3D Silicon Sensors. Fast Current Amplifiers
Parker, Sherwood; Kok, Angela; Kenney, Christopher; Jarron, Pierre; Hasi, Jasmine; Despeisse, Matthieu; Da Via, Cinzia; Anelli, Giovanni; /CERN
2012-05-07
The authors describe techniques to make fast, sub-nanosecond time resolution solid-state detector systems using sensors with 3D electrodes, current amplifiers, constant-fraction comparators or fast wave-form recorders, and some of the next steps to reach still faster results.
Fast sawtooth reconnection at realistic Lundquist numbers
NASA Astrophysics Data System (ADS)
Günter, S.; Yu, Q.; Lackner, K.; Bhattacharjee, A.; Huang, Y.-M.
2015-01-01
Magnetic reconnection, a ubiquitous phenomenon in astrophysics, space science and magnetic confinement research, frequently proceeds much faster than predicted by simple resistive MHD theory. Acceleration can result from the break-up of the thin Sweet-Parker current sheet into plasmoids, or from two-fluid effects decoupling mass and magnetic flux transport over the ion inertial length {{v}A}/{ωci} or the drift scale \\sqrt{{{T}e}/{{m}i}}/{ωci}, depending on the absence or presence of a strong magnetic guide field. We describe new results on the modelling of sawtooth reconnection in a simple tokamak geometry (circular cylindrical equilibrium) pushed to realistic Lundquist numbers for present day tokamaks. For the resistive MHD case, the onset criteria and the influence of plasmoids on the reconnection process agree well with earlier results found in the case of vanishing magnetic guide fields. While plasmoids are also observed in two-fluid calculations, they do not dominate the reconnection process for the range of plasma parameters considered in this study. In the two-fluid case they form as a transient phenomenon only. The reconnection times become weakly dependent on the S-value and for the most complete model—including two-fluid effects and equilibrium temperature and density gradients—agree well with those experimentally found on ASDEX Upgrade ≤ft(≤slant 100 μ s\\right).
Trigger of Fast Reconnection via Collapsing Current Sheets
NASA Astrophysics Data System (ADS)
Tenerani, A.; Velli, M.; Rappazzo, A. F.; Pucci, F.
2015-12-01
It has been widely believed that reconnection is the underlying mechanism of many explosive processes observed both in astrophysical and laboratory plasmas. However, both the questions of how magnetic reconnection is triggered in high Lundquist (S) and Reynolds (R) number plasmas, and how it can then occur on fast, ideal, time-scales remain open. Indeed, it has been argued that fast reconnection rates could be achieved once kinetic scales are reached, or, alternatively, by the onset of the so-called plasmoid instability within Sweet-Parker current sheets. However, it has been shown recently that a tearing mode instability (the "ideal tearing") can grow on an ideal, i.e., S-independent, timescale once the width a of a current sheet becomes thin enough with respect to its macroscopic length L, a/L ~ S-1/3. This suggests that current sheet thinning down to such a threshold aspect ratio —much larger, for S>>1, than the Sweet-Parker one that scales as a/L ~ S-1/2— might provide the trigger for fast reconnection even within the fluid plasma framework. Here we discuss the transition to fast reconnection by studying with visco-resistive MHD simulations the onset and evolution of the tearing instability within a single collapsing current sheet. We indeed show that the transition to a fast tearing mode instability takes place when an inverse aspect ratio of the order of the threshold a/L ~ S-1/3 is reached, and that the secondary current sheets forming nonlinearly become the source of a succession of recursive tearing instabilities. The latter is reminiscent of the fractal reconnection model of flares, which we modify in the light of the "ideal tearing" scenario.
NASA Astrophysics Data System (ADS)
Buechner, J.; Jain, N.; Sharma, A.
2013-12-01
The four s/c of the Magnetospheric Multiscale (MMS) mission, to be launched in 2014, will use the Earth's magnetosphere as a laboratory to study the microphysics of three fundamental plasma processes. One of them is magnetic reconnection, an essentially multi-scale process. While laboratory experiments and past theoretical investigations have shown that important processes necessary to understand magnetic reconnection take place at electron scales the MMS mission for the first time will be able to resolve these scales by in space observations. For the measurement strategy of MMS it is important to make specific predictions of the behavior of current sheets with a thickness of the order of the electron skin depth which play an important role in the evolution of collisionless magnetic reconnection. Since these processes are highly nonlinear and non-local numerical simulation is needed to specify the current sheet evolution. Here we present new results about the nonlinear evolution of electron-scale current sheets starting from the linear stage and using 3-D electron-magnetohydrodynamic (EMHD) simulations. The growth rates of the simulated instabilities compared well with the growth rates obtained from linear theory. Mechanisms and conditions of the formation of flux ropes and of current filamentation will be discussed in comparison with the results of fully kinetic simulations. In 3D the X- and O-point configurations of the magnetic field formed in reconnection planes alternate along the out-of-reconnection-plane direction with the wavelength of the unstable mode. In the presence of multiple reconnection sites, the out-of-plane magnetic field can develop nested structure of quadrupoles in reconnection planes, similar to the 2-D case, but now with variations in the out-of-plane direction. The structures of the electron flow and magnetic field in 3-D simulations will be compared with those in 2-D simulations to discriminate the essentially 3D features. We also discuss
Fast and precise 3D fluorophore localization by gradient fitting
NASA Astrophysics Data System (ADS)
Ma, Hongqiang; Xu, Jianquan; Jin, Jingyi; Gao, Ying; Lan, Li; Liu, Yang
2016-02-01
Astigmatism imaging is widely used to encode the 3D position of fluorophore in single-particle tracking and super-resolution localization microscopy. Here, we present a fast and precise localization algorithm based on gradient fitting to decode the 3D subpixel position of the fluorophore. This algorithm determines the center of the emitter by finding the position with the best-fit gradient direction distribution to the measured point spread function (PSF), and can retrieve the 3D subpixel position of the emitter in a single iteration. Through numerical simulation and experiments with mammalian cells, we demonstrate that our algorithm yields comparable localization precision to the traditional iterative Gaussian function fitting (GF) based method, while exhibits over two orders-of-magnitude faster execution speed. Our algorithm is a promising online reconstruction method for 3D super-resolution microscopy.
Instability, Turbulence, and 3D Magnetic Reconnection in a Line-Tied, Zero Net Current Screw Pinch
NASA Astrophysics Data System (ADS)
Brookhart, Matthew I.; Stemo, Aaron; Zuberbier, Amanda; Zweibel, Ellen; Forest, Cary B.
2015-04-01
This Letter reports the first experimental investigation into a line-tied plasma with a reversed current profile. Discrete current sources create a cylindrical plasma equilibrium with an axial field and zero net current. Detailed magnetic measurements show that an internal m =1 mode with no external character grows exponentially. The nonlinear evolution of the mode drives 3D reconnection events that reorganize the plasma equilibrium. The plasma is turbulent and exhibits reconnection events on a range of scales. These data are consistent with recent simulations of coronal loops and the nanoflare coronal heating mechanism.
Connecting Global Measures of 3D Magnetic Reconnection to Local Kinetic Physics
Daughton, William Scott
2015-07-16
After giving the motivation for the work, slides present the topic under the following headings: Description of LAPD experiment; Actual simulation setup; Simple kinetic theory of ined-tied tearing; Diagnostics to characterizing 3D reconnection; Example #1 - short-tied system; and Example #2 - long line-tied system. Colorful simulations are shown for quasipotential vs field line exponentiation, field line integrated Ohms Law, and correlation with agyrotopy & energy conversion for example #1; and evolution of current density for largest case, field exponentiation vs quasi-potential, and time evolution of magnetic field lines for example #2. To satisfy line-tied boundary conditions, there is need for superposition of oblique modes--the simple two-mode approximation works surprisingly well. For force-free layers with b_{g} >1, the fastest growing periodic modes are oblique with k_{x}λ ~0.5. This implies a minimum length of L_{y} > 2πλb_{g}. There are strong correlations between σ → Ξ → A_{0e} (observable with spacecraft). Electron pressure tensor is the dominant non-ideal term.
Fast Mode Decision for 3D-HEVC Depth Intracoding
Li, Nana; Wu, Qinggang
2014-01-01
The emerging international standard of high efficiency video coding based 3D video coding (3D-HEVC) is a successor to multiview video coding (MVC). In 3D-HEVC depth intracoding, depth modeling mode (DMM) and high efficiency video coding (HEVC) intraprediction mode are both employed to select the best coding mode for each coding unit (CU). This technique achieves the highest possible coding efficiency, but it results in extremely large encoding time which obstructs the 3D-HEVC from practical application. In this paper, a fast mode decision algorithm based on the correlation between texture video and depth map is proposed to reduce 3D-HEVC depth intracoding computational complexity. Since the texture video and its associated depth map represent the same scene, there is a high correlation among the prediction mode from texture video and depth map. Therefore, we can skip some specific depth intraprediction modes rarely used in related texture CU. Experimental results show that the proposed algorithm can significantly reduce computational complexity of 3D-HEVC depth intracoding while maintaining coding efficiency. PMID:24963512
Analysis of Fast Reconnection in MRX
NASA Astrophysics Data System (ADS)
Yamada, M.; Ji, H.; Kulsrud, R.; Kuritsyn, A.; Terry, S.
2002-11-01
In MRX (Magnetic Reconnection Experiment), the detailed magnetic field structure and the plasma profiles of the neutral sheet have been measured and evidence for non-MHD physics mechanisms has been obtained(M. Yamada et al., Phys. Plasmas, 7), 1781 (2000). As the electron-ion collision frequency is reduced, the reconnection rate has been observed to be significantly enhanced over the classical value. Both electrostatic and magnetic turbulence of lower hybrid frequency range has been observed in the low collisionality regime. The amplitude of magnetic fluctuation is the largest at the center of neutral sheet while electropotential fluctuations peak at the edge of the sheath. To describe the physics of this region, much thought has been given to the `` generalized" Ohm's law where electrons and ions are treated separately. In this `` E-MHD" region electrons but not ions are considered to be magnetized and the lower hybrid wave can be connected continuously to an obliquely propagating Whistler wave(R.L. Stenzel and J.M. Urrutia, Phys. Plasmas, 7), 4450 (2000). In this paper we examine the entire spectrum of waves from the electron cyclotron frequency down to the ion cyclotron frequency, focusing on comparisons with the recent experimental data from the MRX neutral sheet plasmas(H. Ji, et al., this meeting). The physical mechanisms of the observed enhanced resistivity due to these waves will be discussed. This work is supported by DoE, NASA and NSF.
NASA Astrophysics Data System (ADS)
Daldorff, L. K. S.; Toth, G.; Borovikov, D.; Gombosi, T. I.; Lapenta, G.
2014-12-01
With the new modeling capability in the Space Weather Modeling Framework (SWMF) of embedding an implicit Particle-in-Cell (PIC) model iPIC3D into the BATS-R-US magnetohydrodynamics model (Daldorff et al. 2014, JCP, 268, 236) we are ready to locally handle the full physics of the reconnection and its implications on the full system where globally, away from the reconnection region, a magnetohydrodynamic description is satisfactory. As magnetic reconnection is one of the main drivers in magnetospheric and heliospheric plasma dynamics, the self-consistent description of the electron dynamics in the coupled MHD-EPIC model is well suited for investigating the nature of these systems. We will compare the new embedded MHD-EPIC model with pure MHD and Hall MHD simulations of the Earth's magnetosphere.
Fast Wave Trains Associated with Solar Eruptions: Insights from 3D Thermodynamic MHD Simulations
NASA Astrophysics Data System (ADS)
Downs, C.; Liu, W.; Torok, T.; Linker, J.; Mikic, Z.; Ofman, L.
2015-12-01
EUV imaging observations during the SDO/AIA era have provided new insights into a variety of wave phenomena occurring in the low solar corona. One example is the observation of quasi-periodic, fast-propagating wave trains that are associated with solar eruptions, including flares and CMEs. While there has been considerable progress in understanding such waves from both an observational and theoretical perspective, it remains a challenge to pin down their physical origin. In this work, we detail our results from a case-study 3D thermodynamic MHD simulation of a coronal mass ejection where quasi-periodic wave trains are generated during the simulated eruption. We find a direct correlation between the onset of non-steady reconnection in the flare current sheet and the generation of quasi-periodic wave train signatures when patchy, collimated downflows interact with the flare arcade. Via forward modeling of SDO/AIA observables, we explore how the appearance of the wave trains is affected by line-of-sight integration and the multi-thermal nature of the coronal medium. We also examine how the wave trains themselves are channeled by natural waveguides formed in 3D by the non-uniform background magnetic field. While the physical association of the reconnection dynamics to the generation of quasi-periodic wave trains appears to be a compelling result, unanswered questions posed from recent observations as well as future prospects will be discussed.
a Fast Method for Measuring the Similarity Between 3d Model and 3d Point Cloud
NASA Astrophysics Data System (ADS)
Zhang, Zongliang; Li, Jonathan; Li, Xin; Lin, Yangbin; Zhang, Shanxin; Wang, Cheng
2016-06-01
This paper proposes a fast method for measuring the partial Similarity between 3D Model and 3D point Cloud (SimMC). It is crucial to measure SimMC for many point cloud-related applications such as 3D object retrieval and inverse procedural modelling. In our proposed method, the surface area of model and the Distance from Model to point Cloud (DistMC) are exploited as measurements to calculate SimMC. Here, DistMC is defined as the weighted distance of the distances between points sampled from model and point cloud. Similarly, Distance from point Cloud to Model (DistCM) is defined as the average distance of the distances between points in point cloud and model. In order to reduce huge computational burdens brought by calculation of DistCM in some traditional methods, we define SimMC as the ratio of weighted surface area of model to DistMC. Compared to those traditional SimMC measuring methods that are only able to measure global similarity, our method is capable of measuring partial similarity by employing distance-weighted strategy. Moreover, our method is able to be faster than other partial similarity assessment methods. We demonstrate the superiority of our method both on synthetic data and laser scanning data.
Simulations of 3D Magnetic Merging: Resistive Scalings for Null Point and QSL Reconnection
NASA Astrophysics Data System (ADS)
Effenberger, Frederic; Craig, I. J. D.
2016-01-01
Starting from an exact, steady-state, force-free solution of the magnetohydrodynamic (MHD) equations, we investigate how resistive current layers are induced by perturbing line-tied three-dimensional magnetic equilibria. This is achieved by the superposition of a weak perturbation field in the domain, in contrast to studies where the boundary is driven by slow motions, like those present in photospheric active regions. Our aim is to quantify how the current structures are altered by the contribution of so-called quasi-separatrix layers (QSLs) as the null point is shifted outside the computational domain. Previous studies based on magneto-frictional relaxation have indicated that despite the severe field line gradients of the QSL, the presence of a null is vital in maintaining fast reconnection. Here, we explore this notion using highly resolved simulations of the full MHD evolution. We show that for the null-point configuration, the resistive scaling of the peak current density is close to J˜η^{-1}, while the scaling is much weaker, i.e. J˜η^{-0.4}, when only the QSL connectivity gradients provide a site for the current accumulation.
A new fast reconnection model in a collisionless regime
Tsiklauri, David
2008-11-15
Based on the first principles [i.e., (i) by balancing the magnetic field advection with the term containing electron pressure tensor nongyrotropic components in the generalized Ohm's law; (ii) using the conservation of mass; and (iii) assuming that the weak magnetic field region width, where electron meandering motion supports electron pressure tensor off-diagonal (nongyrotropic) components, is of the order of electron Larmor radius] a simple model of magnetic reconnection in a collisionless regime is formulated. The model is general, resembling its collisional Sweet-Parker analog in that it is not specific to any initial configuration, e.g., Harris-type tearing unstable current sheet, X-point collapse or otherwise. In addition to its importance from the fundamental point of view, the collisionless reconnection model offers a much faster reconnection rate [M{sub c{sup '}}{sub less}=(c/{omega}{sub pe}){sup 2}/(r{sub L,e}L)] than Sweet-Parker's classical one (M{sub sp}=S{sup -1/2}). The width of the diffusion region (current sheet) in the collisionless regime is found to be {delta}{sub c{sup '}}{sub less}=(c/{omega}{sub pe}){sup 2}/r{sub L,e}, which is independent of the global reconnection scale L and is only prescribed by microphysics (electron inertial length, c/{omega}{sub pe}, and electron Larmor radius, r{sub L,e}). Amongst other issues, the fastness of the reconnection rate alleviates, e.g., the problem of interpretation of solar flares by means of reconnection, as for the typical solar coronal parameters the obtained collisionless reconnection time can be a few minutes, as opposed to Sweet-Parker's equivalent value of less than a day. The new theoretical reconnection rate is compared to the Magnetic Reconnection Experiment device experimental data by Yamada et al. [Phys. Plasmas 13, 052119 (2006)] and Ji et al. [Geophys. Res. Lett. 35, 13106 (2008)], and a good agreement is obtained.
Fast 3D stereo flood simulations in urban areas
NASA Astrophysics Data System (ADS)
Hoes, O.; de Haan, G.; Stelling, G.; van Leeuwen, E.; van Dam, A.; Pleumeekers, O.; Schuurmans, W.
2012-04-01
Flood propagation models are essential to study floods as it is problematic to collect data during actual floods. Moreover, models are needed to explore the consequences of additional scenarios above the actual flood itself. The results of these model studies are generally graphs with water levels over time for certain locations or maps with the flood extent in an area for different return periods. A main difficulty with these depictions of flood information is that they do not reflect the seriousness of flood impacts well in terms of life-like human experience. Typically, one needs a (near) flood before measures are implemented. Apparently, a graph or map is not the proper material to convince politicians and policy makers, even if they live in the threatened area. The recent introduction of commercially available 3D stereo projectors and high resolution elevation data make it possible to build life-like visualizations of simulations. In our research we explored using 3D stereo, the recently collected elevation data of the Netherlands (20 laser points per m2!) in combination with aerial photographs, and a new fast 2D flood propagation calculation scheme. This scheme (under construction) is able to simulate floods using such high amounts of data points. The model simulates flood propagation on an irregular grid; at locations with large elevation differences (e.g. in urban areas) and fast flowing water, smaller cells are used compared to flat surfaces where the water is not or hardly flowing. The result of our combination is a very detailed flood simulation model that can be used to simulate floods within a fraction of the current calculation time. The opportunities of models and their results increase enormously with fast calculations and visualizations combined. For example, the model allows on the spot exploration of measures during a flood, with the 3D visualization ensuring that flood impacts become clear for decision makers. We will show the preliminary
Fast magnetic reconnection due to anisotropic electron pressure
NASA Astrophysics Data System (ADS)
Cassak, Paul; Baylor, Robert; Fermo, Raymond; Beidler, Matthew; Shay, Michael; Swisdak, Marc; Drake, James; Karimabadi, Homa
2015-11-01
A new regime of fast magnetic reconnection with an out-of-plane (guide) magnetic field is reported in which the key role is played by an electron pressure anisotropy described by the Chew-Goldberger-Low gyrotropic equations of state in the generalized Ohm's law, which even dominates the Hall term. A description of the physical cause of this behavior is provided and two-dimensional fluid simulations are used to confirm the results. The electron pressure anisotropy causes the out-of-plane magnetic field to develop a quadrupole structure of opposite polarity to the Hall magnetic field and gives rise to dispersive waves. In addition to being important for understanding what causes reconnection to be fast, this mechanism should dominate in plasmas with low plasma beta and a high in-plane plasma beta with electron temperature comparable to or larger than ion temperature, so it could be relevant in the solar wind and some tokamaks.
Fast magnetic reconnection due to anisotropic electron pressure
NASA Astrophysics Data System (ADS)
Cassak, P. A.; Baylor, R. N.; Fermo, R. L.; Beidler, M. T.; Shay, M. A.; Swisdak, M.; Drake, J. F.; Karimabadi, H.
2015-02-01
A new regime of fast magnetic reconnection with an out-of-plane (guide) magnetic field is reported in which the key role is played by an electron pressure anisotropy described by the Chew-Goldberger-Low gyrotropic equations of state in the generalized Ohm's law, which even dominates the Hall term. A description of the physical cause of this behavior is provided and two-dimensional fluid simulations are used to confirm the results. The electron pressure anisotropy causes the out-of-plane magnetic field to develop a quadrupole structure of opposite polarity to the Hall magnetic field and gives rise to dispersive waves. In addition to being important for understanding what causes reconnection to be fast, this mechanism should dominate in plasmas with low plasma beta and a high in-plane plasma beta with electron temperature comparable to or larger than ion temperature, so it could be relevant in the solar wind and some tokamaks.
Fast magnetic reconnection due to anisotropic electron pressure
Cassak, P. A.; Baylor, R. N.; Fermo, R. L.; Beidler, M. T.; Shay, M. A.; Swisdak, M.; Drake, J. F.; Karimabadi, H.
2015-02-15
A new regime of fast magnetic reconnection with an out-of-plane (guide) magnetic field is reported in which the key role is played by an electron pressure anisotropy described by the Chew-Goldberger-Low gyrotropic equations of state in the generalized Ohm's law, which even dominates the Hall term. A description of the physical cause of this behavior is provided and two-dimensional fluid simulations are used to confirm the results. The electron pressure anisotropy causes the out-of-plane magnetic field to develop a quadrupole structure of opposite polarity to the Hall magnetic field and gives rise to dispersive waves. In addition to being important for understanding what causes reconnection to be fast, this mechanism should dominate in plasmas with low plasma beta and a high in-plane plasma beta with electron temperature comparable to or larger than ion temperature, so it could be relevant in the solar wind and some tokamaks.
3D magnetic field configuration of small-scale reconnection events in the solar plasma atmosphere
Shimizu, T.
2015-10-15
The outer solar atmosphere, i.e., the corona and the chromosphere, is replete with small energy-release events, which are accompanied by transient brightening and jet-like ejections. These events are considered to be magnetic reconnection events in the solar plasma, and their dynamics have been studied using recent advanced observations from the Hinode spacecraft and other observatories in space and on the ground. These events occur at different locations in the solar atmosphere and vary in their morphology and amount of the released energy. The magnetic field configurations of these reconnection events are inferred based on observations of magnetic fields at the photospheric level. Observations suggest that these magnetic configurations can be classified into two groups. In the first group, two anti-parallel magnetic fields reconnect to each other, yielding a 2D emerging flux configuration. In the second group, helical or twisted magnetic flux tubes are parallel or at a relative angle to each other. Reconnection can occur only between anti-parallel components of the magnetic flux tubes and may be referred to as component reconnection. The latter configuration type may be more important for the larger class of small-scale reconnection events. The two types of magnetic configurations can be compared to counter-helicity and co-helicity configurations, respectively, in laboratory plasma collision experiments.
3D Dynamics of Magnetopause Reconnection Using Hall-MHD Global Simulations
NASA Astrophysics Data System (ADS)
Maynard, K.; Germaschewski, K.; Raeder, J.; Bhattacharjee, A.
2011-12-01
Magnetic reconnection at Earth's magnetopause and in the magnetotail is of crucial importance for the dynamics of the global magnetosphere and space weather. Even though the plasma conditions in the magnetosphere are largely in the collisionless regime, most of the existing research using global computational models employ single-fluid magnetohydrodynamics (MHD) with artificial resistivity. Studies of reconnection in simplified, two-dimensional geometries have established that two-fluid and kinetic effects can dramatically alter dynamics and reconnection rates when compared with single-fluid models. These enhanced models also introduce particular signatures, for example a quadrupolar out-of-plane magnetic field component that has already been observed in space by satellite measurements. However, results from simplified geometries cannot be translated directly to the dynamics of three-dimensional magnetospheric reconnection. For instance, magnetic flux originating from the solar wind and arriving at the magnetopause can either reconnect or be advected around the magnetosphere. In this study, we use a new version of the OpenGGCM code that incorporates the Hall term in a Generalized Ohm's Law to study magnetopause reconnection under synthetic solar wind conditions and investigate how reconnection rates and dynamics of flux transfer events depend on the strength of the Hall term. The OpenGGCM, a global model of Earth's magnetosphere, has recently been ported to exploit modern computing architectures like the Cell processor and SIMD capabilities of conventional processors using an automatic code generator. These enhancements provide us with the performance needed to include the computationally expensive Hall physics.
Plasmoid dynamics in 3D resistive MHD simulations of magnetic reconnection
NASA Astrophysics Data System (ADS)
Samtaney, R.; Loureiro, N. F.; Uzdensky, D. A.; Schekochihin, A. A.
2012-04-01
Magnetic reconnection is a well known plasma process believed to lie at the heart of a variety of phenomena such as sub-storms in the Earth's magnetosphere, solar/stellar and accretion-disk flares, sawteeth activity in fusion devices, etc. During reconnection, the global magnetic field topology changes rapidly, leading to the violent release of magnetic energy. Over the past few years, the basic understanding of this fundamental process has undergone profound changes. The validity of the most basic, and widely accepted, reconnection paradigm - the famous Sweet-Parker (SP) model, which predicts that, in MHD, reconnection is extremely slow, its rate scaling as S-1/2, where S is the Lundquist number of the system - has been called into question as it was analytically demonstrated that, for S ≫ 1, SP-like current sheets are violently unstable to the formation of a large number of secondary islands, or plasmoids. Subsequent numerical simulations in 2D have confirmed the validity of the linear theory, and shown that plasmoids quickly grow to become wider than the thickness of the original SP current sheet, thus effectively changing the underlying reconnection geometry. Ensuing numerical work has revealed that the process of plasmoid formation, coalescence and ejection from the sheet drastically modifies the steady state picture assumed by Sweet and Parker, and leads to the unexpected result that MHD reconnection is independent of S. In this talk, we review these recent developments and present results from three-dimensional simulations of high-Lundquist number reconnection in the presence of a guide field. A parametric study varying the strength of the guide field is presented. Plasmoid flux and width distribution functions are quantified and compared with corresponding two dimensional simulations.
NASA Astrophysics Data System (ADS)
Giles, B. L.; Pollock, C. J.; Avanov, L. A.; Barrie, A. C.; Burch, J. L.; Chandler, M. O.; Clark, G. B.; Coffey, V. N.; Dickson, C.; Dorelli, J.; Ergun, R. E.; Fuselier, S. A.; Gershman, D. J.; Gliese, U.; Holland, M. P.; Jacques, A. D.; Kreisler, S.; Lavraud, B.; MacDonald, E.; Mauk, B.; Moore, T. E.; Mukai, T.; Nakamura, R.; Paterson, W. R.; Rager, A. C.; Saito, Y.; Salo, C.; Sauvaud, J. A.; Torbert, R. B.; Vinas, A. F.; Yokota, S.
2015-12-01
The primary focus of the Magnetospheric Multiscale (MMS) mission, launched in March 2015, is magnetic reconnection and associated processes. Understanding hinges critically on the kinetic physics that allows reconnection to take place. The Fast Plasma Investigation (FPI) provides electron and ion distribution functions at 4.5s cadence and, for select periods of time, at cadences of 30ms for electrons and 150ms for ions. These select time periods are chosen after in situ acquisition based on inspection of the low resolution data. Thus the FPI provides, independent of spacecraft spin rate, the time resolution needed to resolve the small, fast-moving reconnection diffusion regions. The first mission phase focuses on the dayside magnetopause and this presentation is intended to demonstrate the capabilities of FPI to resolve the important spatial scales relevant to the reconnection process. Magnetopause and other boundary crossings will be examined and the phase-space trajectories identified at the tetrahedral satellite locations through analysis of the 3D distribution functions.
Undersampling k-space using fast progressive 3D trajectories.
Spiniak, Juan; Guesalaga, Andres; Mir, Roberto; Guarini, Marcelo; Irarrazaval, Pablo
2005-10-01
In 3D MRI, sampling k-space with traditional trajectories can be excessively time-consuming. Fast imaging trajectories are used in an attempt to efficiently cover the k-space and reduce the scan time without significantly affecting the image quality. In many applications, further reductions in scan time can be achieved via undersampling of the k-space; however, no clearly optimal method exists. In most 3D trajectories the k-space is divided into regions that are sampled with shots that share a common geometry (e.g., spirals). A different approach is to design trajectories that gradually but uniformly cover the k-space. In the current work, successive shots progressively add sampled regions to the 3D frequency space. By cutting the sequence short, a natural undersampled method is obtained. This can be particularly efficient because in these types of trajectories the contribution of new information by later shots is less significant. In this work the performance of progressive trajectories for different degrees of undersampling is assessed with trajectories based on missile guidance (MG) ideas. The results show that the approach can be efficient in terms of reducing the scan time, and performs better than the stack of spirals (SOS) technique, particularly under nonideal conditions. PMID:16142719
Fast magnetic reconnection with plasmoid / current sheet ejection events in laboratory experiments
NASA Astrophysics Data System (ADS)
Inomoto, Michiaki; Ono, Yasushi; Hayashi, Yoshinori
2012-07-01
Non-steady and fast magnetic reconnections due to plasmoid or current sheet ejection events have been investigated in laboratory experiments using TS-3, TS-4 and UTST plasma merging devices in the University of Tokyo. In these devices, magnetic reconnection is induced by two different schemes, a) push reconnection driven by flux injection from the upstream region, b) pull reconnection driven by flux extraction to the downstream region. Current sheet or plasmoid ejection events are observed in these reconnection experiments particularly with strong guide magnetic field parallel to the reconnection electric field. In push reconnection experiments, anomalous resistivity is induced by the ion's kinetic effect (meandering motion) when the current sheet width is compressed shorter than the ion gyroradius by the strongly injected inflow flux. This fast reconnection regime does not involve plasmoid / current sheet ejection events. On the other hand, the guide field reduces the ion gyroradius and suppresses the onset of the anomalous resistivity, providing slow and steady magnetic reconnection. Impulsive fast reconnection with strong guide field develops, nevertheless, due to plasmoid / current sheet ejection events in pull and push reconnection experiments with extremely large external driving forces. In such a situation, the inflow flux is forcedly pushed into the reconnection region even faster than the maximal reconnection rate, resulting in flux pile up in front of the diffusion region. This piled flux induces large current density inside the current sheet in which plasmoid structure with closed flux surface is formed in pull reconnection case. The induced large current density or plasmoid is then ejected from the diffusion region with significant increase of reconnection electric field. As a result, magnetic reconnection condition with even larger reconnection rate than that obtained by anomalous resistivity was achieved under strong guide field and large external
Generation of field-aligned currents and Alfven waves by 3D magnetic reconnection
Ma, Z.W.; Lee, L.C.; Otto, A.
1995-07-01
The authors have carried out a three-dimensional compressible MHD simulation to study the generation of field-aligned currents (FAC`s) and Alfven waves by magnetic reconnection for locally antiparallel magnetic fields across the current sheet. Reconnection is triggered by a localized resistivity. The results indicate that both FAC`s and Alfven waves are generated by the three-dimensional reconnection process. Two pairs of FAC`s are generated on each side of current sheet. The polarities of the resulting FAC pair in the leading bulge region are opposite to those of a FAC pair in the trailing quasi-steady region. It is further found that a large portion of the FAC`s ({approximately}40%) is located in the closed field line region. They examine the Walen relation between FAC and parallel vorticity and find that Alfven waves are generated and propagate away from the reconnection site. They discuss the relevance of the results to the observed Region 1 FAC`s at noon. 15 refs., 4 figs.
Fast fully 3-D image reconstruction in PET using planograms.
Brasse, D; Kinahan, P E; Clackdoyle, R; Defrise, M; Comtat, C; Townsend, D W
2004-04-01
We present a method of performing fast and accurate three-dimensional (3-D) backprojection using only Fourier transform operations for line-integral data acquired by planar detector arrays in positron emission tomography. This approach is a 3-D extension of the two-dimensional (2-D) linogram technique of Edholm. By using a special choice of parameters to index a line of response (LOR) for a pair of planar detectors, rather than the conventional parameters used to index a LOR for a circular tomograph, all the LORs passing through a point in the field of view (FOV) lie on a 2-D plane in the four-dimensional (4-D) data space. Thus, backprojection of all the LORs passing through a point in the FOV corresponds to integration of a 2-D plane through the 4-D "planogram." The key step is that the integration along a set of parallel 2-D planes through the planogram, that is, backprojection of a plane of points, can be replaced by a 2-D section through the origin of the 4-D Fourier transform of the data. Backprojection can be performed as a sequence of Fourier transform operations, for faster implementation. In addition, we derive the central-section theorem for planogram format data, and also derive a reconstruction filter for both backprojection-filtering and filtered-backprojection reconstruction algorithms. With software-based Fourier transform calculations we provide preliminary comparisons of planogram backprojection to standard 3-D backprojection and demonstrate a reduction in computation time by a factor of approximately 15. PMID:15084067
3-D lookup: Fast protein structure database searches
Holm. L.; Sander, C.
1995-12-31
There are far fewer classes of three-dimensional protein folds than sequence families but the problem of detecting three-dimensional similarities is NP-complete. We present a novel heuristic for identifying 3-D similarities between a query structure and the database of known protein structures. Many methods for structure alignment use a bottom-up approach, identifying first local matches and then solving a combinatorial problem in building up larger clusters of matching substructures. Here the top-down approach is to start with the global comparison and select a rough superimposition using a fast 3-D lookup of secondary structure motifs. The superimposition is then extended to an alignment of C{sup {alpha}} atoms by an iterative dynamic programming step. An all-against-all comparison of 385-representative proteins (150,000 pair comparisons) took 1 day of computer time on a single R8000 processor. In other words, one query structure is scanned against the database in a matter of minutes. The method is rated at 90% reliability at capturing statistically significant similarities. It is useful as a rapid preprocessor to a comprehensive protein structure database search system.
Models of coronal heating, turbulence and fast reconnection.
Velli, M; Pucci, F; Rappazzo, F; Tenerani, A
2015-05-28
Coronal heating is at the origin of the EUV and X-ray emission and mass loss from the sun and many other stars. While different scenarios have been proposed to explain the heating of magnetically confined and open regions of the corona, they must all rely on the transfer, storage and dissipation of the abundant energy present in photospheric motions, which, coupled to magnetic fields, give rise to the complex phenomenology seen at the chromosphere and transition region (i.e. spicules, jets, 'tornadoes'). Here we discuss models and numerical simulations which rely on magnetic fields and electric currents both for energy transfer and for storage in the corona. We will revisit the sources and frequency spectrum of kinetic and electromagnetic energies, the role of boundary conditions, and the routes to small scales required for effective dissipation. Because reconnection in current sheets has been, and still is, one of the most important processes for coronal heating, we will also discuss recent aspects concerning the triggering of reconnection instabilities and the transition to fast reconnection. PMID:25897086
Models of coronal heating, turbulence and fast reconnection
NASA Astrophysics Data System (ADS)
Velli, M.; Pucci, F.; Rappazzo, F.; Tenerani, A.
2015-04-01
Coronal heating is at the origin of the EUV and X-ray emission and mass loss from the sun and many other stars. While different scenarios have been proposed to explain the heating of magnetically confined and open regions of the corona, they must all rely on the transfer, storage and dissipation of the abundant energy present in photospheric motions, which, coupled to magnetic fields, give rise to the complex phenomenology seen at the chromosphere and transition region (i.e. spicules, jets, 'tornadoes'). Here we discuss models and numerical simulations which rely on magnetic fields and electric currents both for energy transfer and for storage in the corona. We will revisit the sources and frequency spectrum of kinetic and electromagnetic energies, the role of boundary conditions, and the routes to small scales required for effective dissipation. Because reconnection in current sheets has been, and still is, one of the most important processes for coronal heating, we will also discuss recent aspects concerning the triggering of reconnection instabilities and the transition to fast reconnection.
Fast reconnection in relativistic plasmas: the magnetohydrodynamics tearing instability revisited
NASA Astrophysics Data System (ADS)
Del Zanna, L.; Papini, E.; Landi, S.; Bugli, M.; Bucciantini, N.
2016-08-01
Fast reconnection operating in magnetically dominated plasmas is often invoked in models for magnetar giant flares, for magnetic dissipation in pulsar winds, or to explain the gamma-ray flares observed in the Crab nebula, hence its investigation is of paramount importance in high-energy astrophysics. Here we study, by means of two dimensional numerical simulations, the linear phase and the subsequent nonlinear evolution of the tearing instability within the framework of relativistic resistive magnetohydrodynamics, as appropriate in situations where the Alfven velocity approaches the speed of light. It is found that the linear phase of the instability closely matches the analysis in classical MHD, where the growth rate scales with the Lundquist number S as S^-1/2, with the only exception of an enhanced inertial term due to the thermal and magnetic energy contributions. In addition, when thin current sheets of inverse aspect ratio scaling as S^-1/3 are considered, the so-called "ideal" tearing regime is retrieved, with modes growing independently on S and extremely fast, on only a few light crossing times of the sheet length. The overall growth of fluctuations is seen to solely depend on the value of the background Alfven velocity. In the fully nonlinear stage we observe an inverse cascade towards the fundamental mode, with Petschek-type supersonic jets propagating at the external Alfven speed from the X-point, and a fast reconnection rate at the predicted value R~(ln S)^-1.
Global Extended MHD Studies of Fast Magnetic Reconnection
Breslau J.A.; Jardin, S.C.
2002-09-18
Recent experimental and theoretical results have led to two lines of thought regarding the physical processes underlying fast magnetic reconnection. One is based on the traditional Sweet-Parker model but replaces the Spitzer resistivity with an enhanced resistivity caused by electron scattering by ion acoustic turbulence. The other includes the finite gyroradius effects that enter Ohm's law through the Hall and electron pressure gradient terms. A 2-D numerical study, conducted with a new implicit parallel two-fluid code, has helped to clarify the similarities and differences in predictions between these two models and provides some insight into their respective ranges of validity.
Research of Fast 3D Imaging Based on Multiple Mode
NASA Astrophysics Data System (ADS)
Chen, Shibing; Yan, Huimin; Ni, Xuxiang; Zhang, Xiuda; Wang, Yu
2016-02-01
Three-dimensional (3D) imaging has received increasingly extensive attention and has been widely used currently. Lots of efforts have been put on three-dimensional imaging method and system study, in order to meet fast and high accurate requirement. In this article, we realize a fast and high quality stereo matching algorithm on field programmable gate array (FPGA) using the combination of time-of-flight (TOF) camera and binocular camera. Images captured from the two cameras own a same spatial resolution, letting us use the depth maps taken by the TOF camera to figure initial disparity. Under the constraint of the depth map as the stereo pairs when comes to stereo matching, expected disparity of each pixel is limited within a narrow search range. In the meanwhile, using field programmable gate array (FPGA, altera cyclone IV series) concurrent computing we can configure multi core image matching system, thus doing stereo matching on embedded system. The simulation results demonstrate that it can speed up the process of stereo matching and increase matching reliability and stability, realize embedded calculation, expand application range.
Fast reconnection in relativistic plasmas: the magnetohydrodynamics tearing instability revisited
NASA Astrophysics Data System (ADS)
Del Zanna, L.; Papini, E.; Landi, S.; Bugli, M.; Bucciantini, N.
2016-08-01
Fast reconnection operating in magnetically dominated plasmas is often invoked in models for magnetar giant flares, for magnetic dissipation in pulsar winds, or to explain the gamma-ray flares observed in the Crab nebula; hence, its investigation is of paramount importance in high-energy astrophysics. Here we study, by means of two-dimensional numerical simulations, the linear phase and the subsequent non-linear evolution of the tearing instability within the framework of relativistic resistive magnetohydrodynamics (MHD), as appropriate in situations where the Alfvén velocity approaches the speed of light. It is found that the linear phase of the instability closely matches the analysis in classical MHD, where the growth rate scales with the Lundquist number S as S-1/2, with the only exception of an enhanced inertial term due to the thermal and magnetic energy contributions. In addition, when thin current sheets of inverse aspect ratio scaling as S-1/3 are considered, the so-called ideal tearing regime is retrieved, with modes growing independently of S and extremely fast, on only a few light crossing times of the sheet length. The overall growth of fluctuations is seen to solely depend on the value of the background Alfvén velocity. In the fully non-linear stage, we observe an inverse cascade towards the fundamental mode, with Petschek-type supersonic jets propagating at the external Alfvén speed from the X-point, and a fast reconnection rate at the predicted value {R}˜ (ln S)^{-1}.
Three-dimensional fast magnetic reconnection driven by relativistic ultraintense femtosecond lasers.
Ping, Y L; Zhong, J Y; Sheng, Z M; Wang, X G; Liu, B; Li, Y T; Yan, X Q; He, X T; Zhang, J; Zhao, G
2014-03-01
Three-dimensional fast magnetic reconnection driven by two ultraintense femtosecond laser pulses is investigated by relativistic particle-in-cell simulation, where the two paralleled incident laser beams are shot into a near-critical plasma layer to form a magnetic reconnection configuration in self-generated magnetic fields. A reconnection X point and out-of-plane quadrupole field structures associated with magnetic reconnection are formed. The reconnection rate is found to be faster than that found in previous two-dimensional Hall magnetohydrodynamic simulations and electrostatic turbulence contribution to the reconnection electric field plays an essential role. Both in-plane and out-of-plane electron and ion accelerations up to a few MeV due to the magnetic reconnection process are also obtained. PMID:24730781
Fast magnetic reconnection supported by sporadic small-scale Petschek-type shocks
Shibayama, Takuya Nakabou, Takashi; Kusano, Kanya; Miyoshi, Takahiro; Vekstein, Grigory
2015-10-15
Standard magnetohydrodynamic (MHD) theory predicts reconnection rate that is far too slow to account for a wide variety of reconnection events observed in space and laboratory plasmas. Therefore, it was commonly accepted that some non-MHD (kinetic) effects play a crucial role in fast reconnection. A recently renewed interest in simple MHD models is associated with the so-called plasmoid instability of reconnecting current sheets. Although it is now evident that this effect can significantly enhance the rate of reconnection, many details of the underlying multiple-plasmoid process still remain controversial. Here, we report results of a high-resolution computer simulation which demonstrate that fast albeit intermittent magnetic reconnection is sustained by numerous small-scale Petschek-type shocks spontaneously formed in the current sheet due to its plasmoid instability.
3D RISM theory with fast reciprocal-space electrostatics
Heil, Jochen; Kast, Stefan M.
2015-03-21
The calculation of electrostatic solute-solvent interactions in 3D RISM (“three-dimensional reference interaction site model”) integral equation theory is recast in a form that allows for a computational treatment analogous to the “particle-mesh Ewald” formalism as used for molecular simulations. In addition, relations that connect 3D RISM correlation functions and interaction potentials with thermodynamic quantities such as the chemical potential and average solute-solvent interaction energy are reformulated in a way that calculations of expensive real-space electrostatic terms on the 3D grid are completely avoided. These methodical enhancements allow for both, a significant speedup particularly for large solute systems and a smoother convergence of predicted thermodynamic quantities with respect to box size, as illustrated for several benchmark systems.
3D RISM theory with fast reciprocal-space electrostatics.
Heil, Jochen; Kast, Stefan M
2015-03-21
The calculation of electrostatic solute-solvent interactions in 3D RISM ("three-dimensional reference interaction site model") integral equation theory is recast in a form that allows for a computational treatment analogous to the "particle-mesh Ewald" formalism as used for molecular simulations. In addition, relations that connect 3D RISM correlation functions and interaction potentials with thermodynamic quantities such as the chemical potential and average solute-solvent interaction energy are reformulated in a way that calculations of expensive real-space electrostatic terms on the 3D grid are completely avoided. These methodical enhancements allow for both, a significant speedup particularly for large solute systems and a smoother convergence of predicted thermodynamic quantities with respect to box size, as illustrated for several benchmark systems. PMID:25796231
Fast vision-based catheter 3D reconstruction.
Moradi Dalvand, Mohsen; Nahavandi, Saeid; Howe, Robert D
2016-07-21
Continuum robots offer better maneuverability and inherent compliance and are well-suited for surgical applications as catheters, where gentle interaction with the environment is desired. However, sensing their shape and tip position is a challenge as traditional sensors can not be employed in the way they are in rigid robotic manipulators. In this paper, a high speed vision-based shape sensing algorithm for real-time 3D reconstruction of continuum robots based on the views of two arbitrary positioned cameras is presented. The algorithm is based on the closed-form analytical solution of the reconstruction of quadratic curves in 3D space from two arbitrary perspective projections. High-speed image processing algorithms are developed for the segmentation and feature extraction from the images. The proposed algorithms are experimentally validated for accuracy by measuring the tip position, length and bending and orientation angles for known circular and elliptical catheter shaped tubes. Sensitivity analysis is also carried out to evaluate the robustness of the algorithm. Experimental results demonstrate good accuracy (maximum errors of ±0.6 mm and ±0.5 deg), performance (200 Hz), and robustness (maximum absolute error of 1.74 mm, 3.64 deg for the added noises) of the proposed high speed algorithms. PMID:27352011
Fast 3D fluid registration of brain magnetic resonance images
NASA Astrophysics Data System (ADS)
Leporé, Natasha; Chou, Yi-Yu; Lopez, Oscar L.; Aizenstein, Howard J.; Becker, James T.; Toga, Arthur W.; Thompson, Paul M.
2008-03-01
Fluid registration is widely used in medical imaging to track anatomical changes, to correct image distortions, and to integrate multi-modality data. Fluid mappings guarantee that the template image deforms smoothly into the target, without tearing or folding, even when large deformations are required for accurate matching. Here we implemented an intensity-based fluid registration algorithm, accelerated by using a filter designed by Bro-Nielsen and Gramkow. We validated the algorithm on 2D and 3D geometric phantoms using the mean square difference between the final registered image and target as a measure of the accuracy of the registration. In tests on phantom images with different levels of overlap, varying amounts of Gaussian noise, and different intensity gradients, the fluid method outperformed a more commonly used elastic registration method, both in terms of accuracy and in avoiding topological errors during deformation. We also studied the effect of varying the viscosity coefficients in the viscous fluid equation, to optimize registration accuracy. Finally, we applied the fluid registration algorithm to a dataset of 2D binary corpus callosum images and 3D volumetric brain MRIs from 14 healthy individuals to assess its accuracy and robustness.
Fast vision-based catheter 3D reconstruction
NASA Astrophysics Data System (ADS)
Moradi Dalvand, Mohsen; Nahavandi, Saeid; Howe, Robert D.
2016-07-01
Continuum robots offer better maneuverability and inherent compliance and are well-suited for surgical applications as catheters, where gentle interaction with the environment is desired. However, sensing their shape and tip position is a challenge as traditional sensors can not be employed in the way they are in rigid robotic manipulators. In this paper, a high speed vision-based shape sensing algorithm for real-time 3D reconstruction of continuum robots based on the views of two arbitrary positioned cameras is presented. The algorithm is based on the closed-form analytical solution of the reconstruction of quadratic curves in 3D space from two arbitrary perspective projections. High-speed image processing algorithms are developed for the segmentation and feature extraction from the images. The proposed algorithms are experimentally validated for accuracy by measuring the tip position, length and bending and orientation angles for known circular and elliptical catheter shaped tubes. Sensitivity analysis is also carried out to evaluate the robustness of the algorithm. Experimental results demonstrate good accuracy (maximum errors of ±0.6 mm and ±0.5 deg), performance (200 Hz), and robustness (maximum absolute error of 1.74 mm, 3.64 deg for the added noises) of the proposed high speed algorithms.
Fast 3D Surface Extraction 2 pages (including abstract)
Sewell, Christopher Meyer; Patchett, John M.; Ahrens, James P.
2012-06-05
Ocean scientists searching for isosurfaces and/or thresholds of interest in high resolution 3D datasets required a tedious and time-consuming interactive exploration experience. PISTON research and development activities are enabling ocean scientists to rapidly and interactively explore isosurfaces and thresholds in their large data sets using a simple slider with real time calculation and visualization of these features. Ocean Scientists can now visualize more features in less time, helping them gain a better understanding of the high resolution data sets they work with on a daily basis. Isosurface timings (512{sup 3} grid): VTK 7.7 s, Parallel VTK (48-core) 1.3 s, PISTON OpenMP (48-core) 0.2 s, PISTON CUDA (Quadro 6000) 0.1 s.
Fast collisionless reconnection and electron heating in strongly magnetized plasmas.
Loureiro, N F; Schekochihin, A A; Zocco, A
2013-07-12
Magnetic reconnection in strongly magnetized (low-beta), weakly collisional plasmas is investigated by using a novel fluid-kinetic model [Zocco and Schekochihin, Phys. Plasmas 18, 102309 (2011)] which retains nonisothermal electron kinetics. It is shown that electron heating via Landau damping (linear phase mixing) is the dominant dissipation mechanism. In time, electron heating occurs after the peak of the reconnection rate; in space, it is concentrated along the separatrices of the magnetic island. For sufficiently large systems, the peak reconnection rate is cE(∥)(max) ≈ 0.2v(A)B(y,0), where v(A) is the Alfvén speed based on the reconnecting field B(y,0). The island saturation width is the same as in magnetohydrodynamics models except for small systems, when it becomes comparable to the kinetic scales. PMID:23889411
Scanning fringe projection for fast 3D inspection
NASA Astrophysics Data System (ADS)
Honegger, Marc; Kahl, Michael; Trunz, Sandra; Rinner, Stefan; Ettemeyer, Andreas; Lambelet, Patrick
2013-04-01
In an earlier paper we have described a concept for high speed 3D inspection using fringe projection techniques. We use a special CMOS camera with 300 x 300 px which can calculate the phase on board. The focus of the first step of development had been a fringe projector, which was able to modulate the projected fringes with up to 250 kHz. In the second step the image acquisition part of the system was developed. In case of 3D measurement with a matrix camera, the camera resp. the measuring object has to be moved stepwise in the lateral direction to cover multiple acquisition areas of the measurement object. Between each image the lateral movement has to correspond to the field of view of the camera. At the intended very high image acquisition rates the high acceleration of the system between each image will lead to inacceptable mechanical forces. In order to obtain a continuous scanning procedure and at the same time to use the performance of a matrix camera, a special lens system was developed. A measurement field 120 mm long and 3 mm wide is imaged onto the camera. The width of the measuring field is imaged onto the 300 rows of the camera, giving a lateral resolution of 10 μm. In the longitudinal direction the 120 mm object length is divided into 12'000 lines to reach the same resolution of 10 μm. Due to technical reasons that will be described later only 240 of the 300 pixel rows on the chip were used. Consequently, with each camera shot 240 separated lines are imaged onto the chip. Between each of these 240 lines there is a large empty space, which is not imaged by the camera. In principle, the camera is operating as 240 single line cameras. Therefore, if the camera is shifted in an inclined direction to the camera orientation over the object, the empty spaces can be recorded as well. In an optimum alignment, the complete measuring object can be scanned in a continuous movement, covering the total length of 120 mm. In this paper we will describe this image
Temporal-spatial modeling of fast-moving and deforming 3D objects
NASA Astrophysics Data System (ADS)
Wu, Xiaoliang; Wei, Youzhi
1998-09-01
This paper gives a brief description of the method and techniques developed for the modeling and reconstruction of fast moving and deforming 3D objects. A new approach using close-range digital terrestrial photogrammetry in conjunction with high speed photography and videography is proposed. A sequential image matching method (SIM) has been developed to automatically process pairs of images taken continuously of any fast moving and deforming 3D objects. Using the SIM technique a temporal-spatial model (TSM) of any fast moving and deforming 3D objects can be developed. The TSM would include a series of reconstructed surface models of the fast moving and deforming 3D object in the form of 3D images. The TSM allows the 3D objects to be visualized and analyzed in sequence. The SIM method, specifically the left-right matching and forward-back matching techniques are presented in the paper. An example is given which deals with the monitoring of a typical blast rock bench in a major open pit mine in Australia. With the SIM approach and the TSM model it is possible to automatically and efficiently reconstruct the 3D images of the blasting process. This reconstruction would otherwise be impossible to achieve using a labor intensive manual processing approach based on 2D images taken from conventional high speed cameras. The case study demonstrates the potential of the SIM approach and the TSM for the automatic identification, tracking and reconstruction of any fast moving and deforming 3D targets.
Fast 3D subsurface imaging with stepped-frequency GPR
NASA Astrophysics Data System (ADS)
Masarik, Matthew P.; Burns, Joseph; Thelen, Brian T.; Sutter, Lena
2015-05-01
This paper investigates an algorithm for forming 3D images of the subsurface using stepped-frequency GPR data. The algorithm is specifically designed for a handheld GPR and therefore accounts for the irregular sampling pattern in the data and the spatially-variant air-ground interface by estimating an effective "ground-plane" and then registering the data to the plane. The algorithm efficiently solves the 4th-order polynomial for the Snell reflection points using a fully vectorized iterative scheme. The forward operator is implemented efficiently using an accelerated nonuniform FFT (Greengard and Lee, 2004); the adjoint operator is implemented efficiently using an interpolation step coupled with an upsampled FFT. The imaging is done as a linearized version of the full inverse problem, which is regularized using a sparsity constraint to reduce sidelobes and therefore improve image localization. Applying an appropriate sparsity constraint, the algorithm is able to eliminate most the surrounding clutter and sidelobes, while still rendering valuable image properties such as shape and size. The algorithm is applied to simulated data, controlled experimental data (made available by Dr. Waymond Scott, Georgia Institute of Technology), and government-provided data with irregular sampling and air-ground interface.
Fast algorithm of 3D median filter for medical image despeckling
NASA Astrophysics Data System (ADS)
Xiong, Chengyi; Hou, Jianhua; Gao, Zhirong; He, Xiang; Chen, Shaoping
2007-12-01
Three-dimensional (3-D) median filtering is very useful to eliminate speckle noise from a medical imaging source, such as functional magnetic resonance imaging (fMRI) and ultrasonic imaging. 3-D median filtering is characterized by its higher computation complexity. N 3(N 3-1)/2 comparison operations would be required for 3-D median filtering with N×N×N window if the conventional bubble-sorting algorithm is adopted. In this paper, an efficient fast algorithm for 3-D median filtering was presented, which considerably reduced the computation complexity for extracting the median of a 3-D data array. Compared to the state-of-the-art, the proposed method could reduce the computation complexity of 3-D median filtering by 33%. It results in efficiently reducing the system delay of the 3-D median filter by software implementation, and the system cost and power consumption by hardware implementation.
Spontaneous three-dimensional magnetic reconnection in merging toroidal plasma experiment
Ii, Toru; Ono, Yasushi
2013-01-15
We investigated a new phenomenon of three-dimensional (3D) magnetic reconnection in TS-4 torus plasma merging experiments by directly measuring the 3D structures of the current sheet. Removal of all toroidal asymmetry of the device reveals that a strong external drive of reconnection inflow increases the toroidal asymmetry of the current sheet only during the reconnection. This spontaneous 3D deformation of the current sheet increases the reconnection outflow as well as the reconnection electric field, probably because local compression of the current sheet to a thickness less than the ion gyroradius triggers its strong dissipation of the current sheet, responsible for the onset of 3D reconnection. These mechanisms indicate that the 3D reconnection is a newly observed spontaneous process of fast reconnection.
Spontaneous three-dimensional magnetic reconnection in merging toroidal plasma experiment
NASA Astrophysics Data System (ADS)
Ii, Toru; Ono, Yasushi
2013-01-01
We investigated a new phenomenon of three-dimensional (3D) magnetic reconnection in TS-4 torus plasma merging experiments by directly measuring the 3D structures of the current sheet. Removal of all toroidal asymmetry of the device reveals that a strong external drive of reconnection inflow increases the toroidal asymmetry of the current sheet only during the reconnection. This spontaneous 3D deformation of the current sheet increases the reconnection outflow as well as the reconnection electric field, probably because local compression of the current sheet to a thickness less than the ion gyroradius triggers its strong dissipation of the current sheet, responsible for the onset of 3D reconnection. These mechanisms indicate that the 3D reconnection is a newly observed spontaneous process of fast reconnection.
Fast Magnetic Reconnection in the Plasmoid-Dominated Regime
Uzdensky, D. A.; Loureiro, N. F.; Schekochihin, A. A.
2010-12-03
A conceptual model of resistive magnetic reconnection via a stochastic plasmoid chain is proposed. The global reconnection rate is shown to be independent of the Lundquist number. The distribution of fluxes in the plasmoids is shown to be an inverse-square law. It is argued that there is a finite probability of emergence of abnormally large plasmoids, which can disrupt the chain (and may be responsible for observable large abrupt events in solar flares and sawtooth crashes). A criterion for the transition from the resistive magnetohydrodynamic to the collisionless regime is provided.
Fast magnetic reconnection in the plasmoid-dominated regime.
Uzdensky, D A; Loureiro, N F; Schekochihin, A A
2010-12-01
A conceptual model of resistive magnetic reconnection via a stochastic plasmoid chain is proposed. The global reconnection rate is shown to be independent of the Lundquist number. The distribution of fluxes in the plasmoids is shown to be an inverse-square law. It is argued that there is a finite probability of emergence of abnormally large plasmoids, which can disrupt the chain (and may be responsible for observable large abrupt events in solar flares and sawtooth crashes). A criterion for the transition from the resistive magnetohydrodynamic to the collisionless regime is provided. PMID:21231473
A fast 3D reconstruction system with a low-cost camera accessory
Zhang, Yiwei; Gibson, Graham M.; Hay, Rebecca; Bowman, Richard W.; Padgett, Miles J.; Edgar, Matthew P.
2015-01-01
Photometric stereo is a three dimensional (3D) imaging technique that uses multiple 2D images, obtained from a fixed camera perspective, with different illumination directions. Compared to other 3D imaging methods such as geometry modeling and 3D-scanning, it comes with a number of advantages, such as having a simple and efficient reconstruction routine. In this work, we describe a low-cost accessory to a commercial digital single-lens reflex (DSLR) camera system allowing fast reconstruction of 3D objects using photometric stereo. The accessory consists of four white LED lights fixed to the lens of a commercial DSLR camera and a USB programmable controller board to sequentially control the illumination. 3D images are derived for different objects with varying geometric complexity and results are presented, showing a typical height error of <3 mm for a 50 mm sized object. PMID:26057407
A fast 3D reconstruction system with a low-cost camera accessory
NASA Astrophysics Data System (ADS)
Zhang, Yiwei; Gibson, Graham M.; Hay, Rebecca; Bowman, Richard W.; Padgett, Miles J.; Edgar, Matthew P.
2015-06-01
Photometric stereo is a three dimensional (3D) imaging technique that uses multiple 2D images, obtained from a fixed camera perspective, with different illumination directions. Compared to other 3D imaging methods such as geometry modeling and 3D-scanning, it comes with a number of advantages, such as having a simple and efficient reconstruction routine. In this work, we describe a low-cost accessory to a commercial digital single-lens reflex (DSLR) camera system allowing fast reconstruction of 3D objects using photometric stereo. The accessory consists of four white LED lights fixed to the lens of a commercial DSLR camera and a USB programmable controller board to sequentially control the illumination. 3D images are derived for different objects with varying geometric complexity and results are presented, showing a typical height error of <3 mm for a 50 mm sized object.
ON THE ROLE OF FAST MAGNETIC RECONNECTION IN ACCRETING BLACK HOLE SOURCES
Singh, C. B.; De Gouveia Dal Pino, E. M.; Kadowaki, L. H. S. E-mail: dalpino@iag.usp.br
2015-01-30
We attempt to explain the observed radio and gamma-ray emission produced in the surroundings of black holes by employing a magnetically dominated accretion flow model and fast magnetic reconnection triggered by turbulence. In earlier work, a standard disk model was used and we refine the model by focusing on the sub-Eddington regime to address the fundamental plane of black hole activity. The results do not change substantially with regard to previous work, ensuring that the details of accretion physics are not relevant in the magnetic reconnection process occurring in the corona. Rather, our work puts fast magnetic reconnection events as a powerful mechanism operating in the core region near the jet base of black hole sources on more solid ground. For microquasars and low-luminosity active galactic nuclei, the observed correlation between radio emission and the mass of the sources can be explained by this process. The corresponding gamma-ray emission also seems to be produced in the same core region. On the other hand, emission from blazars and gamma-ray bursts cannot be correlated to core emission based on fast reconnection.
NASA Astrophysics Data System (ADS)
Del Sarto, Daniele; Pucci, Fulvia; Tenerani, Anna; Velli, Marco
2016-03-01
This paper discusses the transition to fast growth of the tearing instability in thin current sheets in the collisionless limit where electron inertia drives the reconnection process. It has been previously suggested that in resistive MHD there is a natural maximum aspect ratio (ratio of sheet length and breadth to thickness) which may be reached for current sheets with a macroscopic length L, the limit being provided by the fact that the tearing mode growth time becomes of the same order as the Alfvén time calculated on the macroscopic scale. For current sheets with a smaller aspect ratio than critical the normalized growth rate tends to zero with increasing Lundquist number S, while for current sheets with an aspect ratio greater than critical the growth rate diverges with S. Here we carry out a similar analysis but with electron inertia as the term violating magnetic flux conservation: previously found scalings of critical current sheet aspect ratios with the Lundquist number are generalized to include the dependence on the ratio de2/L2, where de is the electron skin depth, and it is shown that there are limiting scalings which, as in the resistive case, result in reconnecting modes growing on ideal time scales. Finite Larmor radius effects are then included, and the rescaling argument at the basis of "ideal" reconnection is proposed to explain secondary fast reconnection regimes naturally appearing in numerical simulations of current sheet evolution.
Fast 3D shape measurement using Fourier transform profilometry without phase unwrapping
NASA Astrophysics Data System (ADS)
Song, Kechen; Hu, Shaopeng; Wen, Xin; Yan, Yunhui
2016-09-01
This paper presents a novel, simple, yet fast 3D shape measurement method using Fourier transform profilometry. Different from the conventional Fourier transform profilometry, this proposed method introduces the binocular stereo vision and employs two image pairs (i.e., original image pairs and fringe image pairs) to restructure 3D shape. In this proposed method, instead of phase unwrapping algorithm, a coarse disparity map is adopted as a constraint condition to realize phase matching using wrapped phase. Since the local phase matching and sub-pixel disparity refinement are proposed to obtain high measuring accuracy, high-quality phase is not required. The validity of the proposed method is verified by experiments.
Fast Reconnection Rates Based on Group Velocity Cones: Whistler Regime and Pair Plasmas
NASA Astrophysics Data System (ADS)
Singh, N.
2009-05-01
Based on the group velocity vector of the whistler mode, we predict the range of whistler-regime reconnection rate depending on the half width (w) of the current sheet (CS. During the reconnection process electromagnetic perturbations (EMPs) are generated in the localized diffusion region (DR, which acts like an antenna and radiates whistler waves for certain range of CS widths. The reconnection structure (exhaust) is approximately the radiation pattern of the DR antenna and it is determined by the group velocity directions. Since the whistler waves originate from the electromagnetic perturbations (EMPs) localized in the DR, we calculate R over a range of the discrete values of the perpendicular wave number (k'') contained in the Fourier spectrum of the EMPs. We have used such calculations to determine the reconnection rates
Frank, Lawrence R.; Jung, Youngkyoo; Inati, Souheil; Tyszka, J. Michael; Wong, Eric C.
2009-01-01
We present an acquisition and reconstruction method designed to acquire high resolution 3D fast spin echo diffusion tensor images while mitigating the major sources of artifacts in DTI - field distortions, eddy currents and motion. The resulting images, being 3D, are of high SNR, and being fast spin echoes, exhibit greatly reduced field distortions. This sequence utilizes variable density spiral acquisition gradients, which allow for the implementation of a self-navigation scheme by which both eddy current and motion artifacts are removed. The result is that high resolution 3D DTI images are produced without the need for eddy current compensating gradients or B0 field correction. In addition, a novel method for fast and accurate reconstruction of the non-Cartesian data is employed. Results are demonstrated in the brains of normal human volunteers. PMID:19778618
A fast method to measure the 3D surface of the human heart
NASA Astrophysics Data System (ADS)
Cao, Yiping; Su, Xianyu; Xiang, Liqun; Chen, Wenjing; Zhang, Qican
2003-12-01
Three-dimensional (3-D) automatic measurement of an object is widely used in many fields. In Biology and Medicine society, it can be applicable for surgery, orthopedics, viscera disease analysis and diagnosis etc. Here a new fast method to measure the 3D surface of human heart is proposed which can provide doctors a lot of information, such as the size of heart profile, the sizes of the left or right heart ventricle, and the curvature center and radius of heart ventricle, to fully analyze and diagnose pathobiology of human heart. The new fast method is optically and noncontacted and based upon the Phase Measurement Profilometry (PMP), which has higher measuring precision. A human heart specimen experiment has verified our method.
Fast Probabilistic Fusion of 3d Point Clouds via Occupancy Grids for Scene Classification
NASA Astrophysics Data System (ADS)
Kuhn, Andreas; Huang, Hai; Drauschke, Martin; Mayer, Helmut
2016-06-01
High resolution consumer cameras on Unmanned Aerial Vehicles (UAVs) allow for cheap acquisition of highly detailed images, e.g., of urban regions. Via image registration by means of Structure from Motion (SfM) and Multi View Stereo (MVS) the automatic generation of huge amounts of 3D points with a relative accuracy in the centimeter range is possible. Applications such as semantic classification have a need for accurate 3D point clouds, but do not benefit from an extremely high resolution/density. In this paper, we, therefore, propose a fast fusion of high resolution 3D point clouds based on occupancy grids. The result is used for semantic classification. In contrast to state-of-the-art classification methods, we accept a certain percentage of outliers, arguing that they can be considered in the classification process when a per point belief is determined in the fusion process. To this end, we employ an octree-based fusion which allows for the derivation of outlier probabilities. The probabilities give a belief for every 3D point, which is essential for the semantic classification to consider measurement noise. For an example point cloud with half a billion 3D points (cf. Figure 1), we show that our method can reduce runtime as well as improve classification accuracy and offers high scalability for large datasets.
Ugai, M.
2009-11-15
It is well known that magnetic pulsations of long periods impulsively occur in accordance with the sudden onset of geomagnetic substorms and drastic enhancement of electrojets in the ionosphere. On the basis of the spontaneous fast reconnection model, the present paper examines the physical mechanism by which both magnetic pulsations and strong electrojets are impulsively driven by the fast (Alfvenic) reconnection jet. When a large-scale plasmoid [or traveling compression region (TCR)], directly caused by the fast reconnection jet, collides with the magnetic loop footpoint, strong electrojets are impulsively driven in a finite extent in the loop footpoint in accordance with the evolution of the current wedge and the generator current circuit. Simultaneously, magnetohydrodynamic (Alfven) waves, accompanied by the TCR, are reflected from the electrojet layer, leading to impulsive magnetic pulsations ahead of the loop footpoint because of the interaction (or resonance) between the reflected waves and the waves traveling toward the footpoint. The pulsations propagate outward in all directions from the source region of the wave reflection, and the pulsation periods are typically estimated to be of several tens of seconds.
3D WKB solution for fast magnetoacoustic wave behaviour around an X-line
NASA Astrophysics Data System (ADS)
McLaughlin, J. A.; Botha, G. J. J.; Régnier, S.; Spoors, D. L.
2016-06-01
Context. We study the propagation of a fast magnetoacoustic wave in a 3D magnetic field created from two magnetic dipoles. The magnetic topology contains an X-line. Aims: We aim to contribute to the overall understanding of MHD wave propagation within inhomogeneous media, specifically around X-lines. Methods: We investigate the linearised, 3D MHD equations under the assumptions of ideal and cold plasma. We utilise the WKB approximation and Charpit's method during our investigation. Results: It is found that the behaviour of the fast magnetoacoustic wave is entirely dictated by the local, inhomogeneous, equilibrium Alfvén speed profile. All parts of the wave experience refraction during propagation, where the magnitude of the refraction effect depends on the location of an individual wave element within the inhomogeneous magnetic field. The X-line, along which the Alfvén speed is identically zero, acts as a focus for the refraction effect. There are two main types of wave behaviour: part of the wave is either trapped by the X-line or escapes the system, and there exists a critical starting region around the X-line that divides these two types of behaviour. For the set-up investigated, it is found that 15.5% of the fast wave energy is trapped by the X-line. Conclusions: We conclude that linear, β = 0 fast magnetoacoustic waves can accumulate along X-lines and thus these will be specific locations of fast wave energy deposition and thus preferential heating. The work here highlights the importance of understanding the magnetic topology of a system. We also demonstrate how the 3D WKB technique described in this paper can be applied to other magnetic configurations.
Fast-shock formation in line-tied magnetic reconnection models of solar flares
NASA Technical Reports Server (NTRS)
Forbes, T. G.
1986-01-01
In a previous study by the author, an approximately stationary fast shock was tentatively identified in a numerical experiment designed to study line-tied magnetic reconnection. Here the evidence for the occurrence of a stationary fast shock is reexamined, and the previous identification is confirmed. In the numerical experiment, line-tied reconnection is modeled by a configuration which produces two supermagnetosonic outflow jets - one directed upward, away from the photosphere, and one directed downward, toward an arcade of closed magnetic loops tied to the photosphere. The fast shock occurs when the downward-directed jet encounters the obstacle formed by the closed loops. Although the existence of a stationary, or nearly stationary, fast shock is confirmed, the transition from the supermagnetosonic flow region upstream of the shock to the nearly static region downstream of the shock is more complicated than was previously thought. Immediately downstream of the shock, there exists a deflection sheath in which the submagnetosonic flow coming out of the shock is diverted around the region of static closed loops. The MHD jump conditions are used to investigate the characteristics of the fast shock and to show that a stationary shock cannot exist unless accompanied by a deflection sheath. Analysis of the shock's location and dimensions suggests that such fast shocks may contribute to particle acceleration and to thermal condensation in flares.
Fast extraction of minimal paths in 3D images and applications to virtual endoscopy.
Deschamps, T; Cohen, L D
2001-12-01
The aim of this article is to build trajectories for virtual endoscopy inside 3D medical images, using the most automatic way. Usually the construction of this trajectory is left to the clinician who must define some points on the path manually using three orthogonal views. But for a complex structure such as the colon, those views give little information on the shape of the object of interest. The path construction in 3D images becomes a very tedious task and precise a priori knowledge of the structure is needed to determine a suitable trajectory. We propose a more automatic path tracking method to overcome those drawbacks: we are able to build a path, given only one or two end points and the 3D image as inputs. This work is based on previous work by Cohen and Kimmel [Int. J. Comp. Vis. 24 (1) (1997) 57] for extracting paths in 2D images using Fast Marching algorithm. Our original contribution is twofold. On the first hand, we present a general technical contribution which extends minimal paths to 3D images and gives new improvements of the approach that are relevant in 2D as well as in 3D to extract linear structures in images. It includes techniques to make the path extraction scheme faster and easier, by reducing the user interaction. We also develop a new method to extract a centered path in tubular structures. Synthetic and real medical images are used to illustrate each contribution. On the other hand, we show that our method can be efficiently applied to the problem of finding a centered path in tubular anatomical structures with minimum interactivity, and that this path can be used for virtual endoscopy. Results are shown in various anatomical regions (colon, brain vessels, arteries) with different 3D imaging protocols (CT, MR). PMID:11731307
Fast 3D Spatial EPR Imaging Using Spiral Magnetic Field Gradient
Deng, Yuanmu; Petryakov, Sergy; He, Guanglong; Kesselring, Eric; Kuppusamy, Periannan; Zweier, Jay L.
2007-01-01
Electron paramagnetic resonance imaging (EPRI) provides direct detection and mapping of free radicals. The continuous wave (CW) EPRI technique, in particular, has been widely used in a variety of applications in the fields of biology and medicine due to its high sensitivity and applicability to a wide range of free radicals and paramagnetic species. However, the technique requires long image acquisition periods, and this limits its use for many in vivo applications where relatively rapid changes occur in the magnitude and distribution of spins. Therefore, there has been a great need to develop fast EPRI techniques. We report the development of a fast 3D CW EPRI technique using spiral magnetic field gradient. By spiraling the magnetic field gradient and stepping the main magnetic field, this approach acquires a 3D image in one sweep of the main magnetic field, enabling significant reduction of the imaging time. A direct one-stage 3D image reconstruction algorithm, modified for reconstruction of the EPR images from the projections acquired with the spiral magnetic field gradient, was used. We demonstrated using a home-built L-band EPR system that the spiral magnetic field gradient technique enabled a 4 to 7-fold accelerated acquisition of projections. This technique has great potential for in vivo studies of free radicals and their metabolism. PMID:17267252
Linear theory for fast collisionless magnetic reconnection in the lower-hybrid frequency range
NASA Astrophysics Data System (ADS)
Jovanović, D.; Shukla, P. K.
2005-05-01
A linear theory is presented for the interplay between the fast collisionless magnetic reconnection and the lower-hybrid waves that has been observed in recent computer simulations [J. F. Drake, M. Swisdak, C. Cattell et al., Science 299, 873 (2003)]. In plasma configurations with a strong guide field and anisotropic electron temperature, the electron dynamics is described within the framework of standard electron magnetohydrodynamic equations, accounting also for the effects of the electron polarization and ion motions in the presence of perpendicular electric fields. In the linear phase, we find two types of instabilities of a thin current sheet with steep edges, corresponding to its filamentation (or tearing) and bending. Using a surface-wave formalism for the perturbations whose wavelength is larger than the thickness of the current sheet, the corresponding growth rates are calculated as the contributions of singularities in the plasma dispersion function. These are governed by the electron inertia and the linear coupling of the reconnecting magnetic field with local plasma modes propagating in the perpendicular direction that are subject to the Buneman instability. The linear surface wave instability may be particularly important as a secondary instability, dissipating the thin current sheets that develop in the course of the fast reconnection in the shear-Alfvén and kinetic-Alfvén regimes, and providing the anomalous resistivity for the growth of magnetic islands beyond the shear-Alfvén and kinetic-Alfvén scales.
Comparison of Parallel MRI Reconstruction Methods for Accelerated 3D Fast Spin-Echo Imaging
Xiao, Zhikui; Hoge, W. Scott; Mulkern, R.V.; Zhao, Lei; Hu, Guangshu; Kyriakos, Walid E.
2014-01-01
Parallel MRI (pMRI) achieves imaging acceleration by partially substituting gradient-encoding steps with spatial information contained in the component coils of the acquisition array. Variable-density subsampling in pMRI was previously shown to yield improved two-dimensional (2D) imaging in comparison to uniform subsampling, but has yet to be used routinely in clinical practice. In an effort to reduce acquisition time for 3D fast spin-echo (3D-FSE) sequences, this work explores a specific nonuniform sampling scheme for 3D imaging, subsampling along two phase-encoding (PE) directions on a rectilinear grid. We use two reconstruction methods—2D-GRAPPA-Operator and 2D-SPACE RIP—and present a comparison between them. We show that high-quality images can be reconstructed using both techniques. To evaluate the proposed sampling method and reconstruction schemes, results via simulation, phantom study, and in vivo 3D human data are shown. We find that fewer artifacts can be seen in the 2D-SPACE RIP reconstructions than in 2D-GRAPPA-Operator reconstructions, with comparable reconstruction times. PMID:18727083
Algorithms for Accurate and Fast Plotting of Contour Surfaces in 3D Using Hexahedral Elements
NASA Astrophysics Data System (ADS)
Singh, Chandan; Saini, Jaswinder Singh
2016-07-01
In the present study, Fast and accurate algorithms for the generation of contour surfaces in 3D are described using hexahedral elements which are popular in finite element analysis. The contour surfaces are described in the form of groups of boundaries of contour segments and their interior points are derived using the contour equation. The locations of contour boundaries and the interior points on contour surfaces are as accurate as the interpolation results obtained by hexahedral elements and thus there are no discrepancies between the analysis and visualization results.
Algorithms for Accurate and Fast Plotting of Contour Surfaces in 3D Using Hexahedral Elements
NASA Astrophysics Data System (ADS)
Singh, Chandan; Saini, Jaswinder Singh
2016-05-01
In the present study, Fast and accurate algorithms for the generation of contour surfaces in 3D are described using hexahedral elements which are popular in finite element analysis. The contour surfaces are described in the form of groups of boundaries of contour segments and their interior points are derived using the contour equation. The locations of contour boundaries and the interior points on contour surfaces are as accurate as the interpolation results obtained by hexahedral elements and thus there are no discrepancies between the analysis and visualization results.
Thermal analysis of 3D composites by a new fast multipole hybrid boundary node method
NASA Astrophysics Data System (ADS)
Miao, Yu; Wang, Qiao; Zhu, Hongping; Li, Yinping
2014-01-01
This paper applies the hybrid boundary node method (Hybrid BNM) for the thermal analysis of 3D composites. A new formulation is derived for the inclusion-based composites. In the new formulation, the unknowns of the interfaces are assembled only once in the final system equation, which can reduce nearly one half of degrees of freedom (DOFs) compared with the conventional multi-domain solver when there are lots of inclusions. A new version of the fast multipole method (FMM) is also coupled with the new formulation and the technique is applied to thermal analysis of composites with many inclusions. In the new fast multipole hybrid boundary node method (FM-HBNM), a diagonal form for translation operators is used and the method presented can be applied to the computation of more than 1,000,000 DOFs on a personal computer. Numerical examples are presented to analyze the thermal behavior of composites with many inclusions.
3D Fast Automatic Segmentation of Kidney Based on Modified AAM and Random Forest.
Jin, Chao; Shi, Fei; Xiang, Dehui; Jiang, Xueqing; Zhang, Bin; Wang, Ximing; Zhu, Weifang; Gao, Enting; Chen, Xinjian
2016-06-01
In this paper, a fully automatic method is proposed to segment the kidney into multiple components: renal cortex, renal column, renal medulla and renal pelvis, in clinical 3D CT abdominal images. The proposed fast automatic segmentation method of kidney consists of two main parts: localization of renal cortex and segmentation of kidney components. In the localization of renal cortex phase, a method which fully combines 3D Generalized Hough Transform (GHT) and 3D Active Appearance Models (AAM) is applied to localize the renal cortex. In the segmentation of kidney components phase, a modified Random Forests (RF) method is proposed to segment the kidney into four components based on the result from localization phase. During the implementation, a multithreading technology is applied to speed up the segmentation process. The proposed method was evaluated on a clinical abdomen CT data set, including 37 contrast-enhanced volume data using leave-one-out strategy. The overall true-positive volume fraction and false-positive volume fraction were 93.15%, 0.37% for renal cortex segmentation; 83.09%, 0.97% for renal column segmentation; 81.92%, 0.55% for renal medulla segmentation; and 80.28%, 0.30% for renal pelvis segmentation, respectively. The average computational time of segmenting kidney into four components took 20 seconds. PMID:26742124
Fast 3D modeling in complex environments using a single Kinect sensor
NASA Astrophysics Data System (ADS)
Yue, Haosong; Chen, Weihai; Wu, Xingming; Liu, Jingmeng
2014-02-01
Three-dimensional (3D) modeling technology has been widely used in inverse engineering, urban planning, robot navigation, and many other applications. How to build a dense model of the environment with limited processing resources is still a challenging topic. A fast 3D modeling algorithm that only uses a single Kinect sensor is proposed in this paper. For every color image captured by Kinect, corner feature extraction is carried out first. Then a spiral search strategy is utilized to select the region of interest (ROI) that contains enough feature corners. Next, the iterative closest point (ICP) method is applied to the points in the ROI to align consecutive data frames. Finally, the analysis of which areas can be walked through by human beings is presented. Comparative experiments with the well-known KinectFusion algorithm have been done and the results demonstrate that the accuracy of the proposed algorithm is the same as KinectFusion but the computing speed is nearly twice of KinectFusion. 3D modeling of two scenes of a public garden and traversable areas analysis in these regions further verified the feasibility of our algorithm.
Fast and robust 3D ultrasound registration--block and game theoretic matching.
Banerjee, Jyotirmoy; Klink, Camiel; Peters, Edward D; Niessen, Wiro J; Moelker, Adriaan; van Walsum, Theo
2015-02-01
Real-time 3D US has potential for image guidance in minimally invasive liver interventions. However, motion caused by patient breathing makes it hard to visualize a localized area, and to maintain alignment with pre-operative information. In this work we develop a fast affine registration framework to compensate in real-time for liver motion/displacement due to breathing. The affine registration of two consecutive ultrasound volumes in time is performed using block-matching. For a set of evenly distributed points in one volume and their correspondences in the other volume, we propose a robust outlier rejection method to reject false matches. The inliers are then used to determine the affine transformation. The approach is evaluated on 13 4D ultrasound sequences acquired from 8 subjects. For 91 pairs of 3D ultrasound volumes selected from these sequences, a mean registration error of 1.8mm is achieved. A graphics processing unit (GPU) implementation runs the 3D US registration at 8 Hz. PMID:25484018
Plasmoid Instabilities Mediated Three-Dimensional Magnetohydrodynamic Turbulent Reconnection
Huang, Yi-min; Guo, Fan
2015-07-21
After some introductory remarks on fast reconnection in resistive MHD due to plasmoid instability, oblique tearing modes in 3D, and previous studies on 3D turbulent reconnection, the subject is presented under the following topics: 3D simulation setup, time evolution of the 3D simulation, comparison with Sweet-Parker and 2D plasmoid reconnection, and diagnostics of the turbulent state (decomposition of mean fields and fluctuations, power spectra of energy fluctuations, structure function and eddy anisotropy with respect to local magnetic field). Three primary conclusions were reached: (1) The results suggest that 3D plasmoid instabilities can lead to self-generated turbulent reconnection (evidence of energy cascade and development of inertial range, energy fluctuations preferentially align with the local magnetic field, which is one of the characteristics of MHD turbulence); (2) The turbulence is highly inhomogeneous, due to the presence of magnetic shear and outflow jets (conventional MHD turbulence theories or phenomenologies may not be applicable – e.g. scale-dependent anisotropy as predicted by Goldreich & Sridhar is not found); (3) 3D turbulent reconnection is different from 2D plasmoid-dominated reconnection in many aspects. However, in fully developed state, reconnection rates in 2D and 3D are comparable — this result needs to be further checked in higher S.
Capabilities of a Global 3D MHD Model for Monitoring Extremely Fast CMEs
NASA Astrophysics Data System (ADS)
Wu, C. C.; Plunkett, S. P.; Liou, K.; Socker, D. G.; Wu, S. T.; Wang, Y. M.
2015-12-01
Since the start of the space era, spacecraft have recorded many extremely fast coronal mass ejections (CMEs) which have resulted in severe geomagnetic storms. Accurate and timely forecasting of the space weather effects of these events is important for protecting expensive space assets and astronauts and avoiding communications interruptions. Here, we will introduce a newly developed global, three-dimensional (3D) magnetohydrodynamic (MHD) model (G3DMHD). The model takes the solar magnetic field maps at 2.5 solar radii (Rs) and intepolates the solar wind plasma and field out to 18 Rs using the algorithm of Wang and Sheeley (1990, JGR). The output is used as the inner boundary condition for a 3D MHD model. The G3DMHD model is capable of simulating (i) extremely fast CME events with propagation speeds faster than 2500 km/s; and (ii) multiple CME events in sequence or simultaneously. We will demonstrate the simulation results (and comparison with in-situ observation) for the fastest CME in record on 23 July 2012, the shortest transit time in March 1976, and the well-known historic Carrington 1859 event.
FAST MAGNETIC RECONNECTION IN THE SOLAR CHROMOSPHERE MEDIATED BY THE PLASMOID INSTABILITY
Ni, Lei; Kliem, Bernhard; Lin, Jun; Wu, Ning
2015-01-20
Magnetic reconnection in the partially ionized solar chromosphere is studied in 2.5 dimensional magnetohydrodynamic simulations including radiative cooling and ambipolar diffusion. A Harris current sheet with and without a guide field is considered. Characteristic values of the parameters in the middle chromosphere imply a high magnetic Reynolds number of ∼10{sup 6}-10{sup 7} in the present simulations. Fast magnetic reconnection then develops as a consequence of the plasmoid instability without the need to invoke anomalous resistivity enhancements. Multiple levels of the instability are followed as it cascades to smaller scales, which approach the ion inertial length. The reconnection rate, normalized to the asymptotic values of magnetic field and Alfvén velocity in the inflow region, reaches values in the range ∼0.01-0.03 throughout the cascading plasmoid formation and for zero as well as for strong guide field. The outflow velocity reaches ≈40 km s{sup –1}. Slow-mode shocks extend from the X-points, heating the plasmoids up to ∼8 × 10{sup 4} K. In the case of zero guide field, the inclusion of both ambipolar diffusion and radiative cooling causes a rapid thinning of the current sheet (down to ∼30 m) and early formation of secondary islands. Both of these processes have very little effect on the plasmoid instability for a strong guide field. The reconnection rates, temperature enhancements, and upward outflow velocities from the vertical current sheet correspond well to their characteristic values in chromospheric jets.
XuvTools: free, fast and reliable stitching of large 3D datasets.
Emmenlauer, M; Ronneberger, O; Ponti, A; Schwarb, P; Griffa, A; Filippi, A; Nitschke, R; Driever, W; Burkhardt, H
2009-01-01
Current biomedical research increasingly requires imaging large and thick 3D structures at high resolution. Prominent examples are the tracking of fine filaments over long distances in brain slices, or the localization of gene expression or cell migration in whole animals like Caenorhabditis elegans or zebrafish. To obtain both high resolution and a large field of view (FOV), a combination of multiple recordings ('tiles') is one of the options. Although hardware solutions exist for fast and reproducible acquisition of multiple 3D tiles, generic software solutions are missing to assemble ('stitch') these tiles quickly and accurately. In this paper, we present a framework that achieves fully automated recombination of tiles recorded at arbitrary positions in 3D space, as long as some small overlap between tiles is provided. A fully automated 3D correlation between all tiles is achieved such that no manual interaction or prior knowledge about tile positions is needed. We use (1) phase-only correlation in a multi-scale approach to estimate the coarse positions, (2) normalized cross-correlation of small patches extracted at salient points to obtain the precise matches, (3) find the globally optimal placement for all tiles by a singular value decomposition and (4) accomplish a nearly seamless stitching by a bleaching correction at the tile borders. If the dataset contains multiple channels, all channels are used to obtain the best matches between tiles. For speedup we employ a heuristic method to prune unneeded correlations, and compute all correlations via the fast Fourier transform (FFT), thereby achieving very good runtime performance. We demonstrate the successful application of the proposed framework to a wide range of different datasets from whole zebrafish embryos and C. elegans, mouse and rat brain slices and fine plant hairs (trichome). Further, we compare our stitching results to those of other commercially and freely available software solutions. The
Park, Sung-Hee; Ryu, Keun Ho; Gilbert, David
2005-08-01
Similarity search for protein 3D structures become complex and computationally expensive due to the fact that the size of protein structure databases continues to grow tremendously. Recently, fast structural similarity search systems have been required to put them into practical use in protein structure classification whilst existing comparison systems do not provide comparison results on time. Our approach uses multi-step processing that composes of a preprocessing step to represent geometry of protein structures with spatial objects, a filter step to generate a small candidate set using approximate topological string matching, and a refinement step to compute a structural alignment. This paper describes the preprocessing and filtering for fast similarity search using the discovery of topological patterns of secondary structure elements based on spatial relations. Our system is fully implemented by using Oracle 8i spatial. We have previously shown that our approach has the advantage of speed of performance compared with other approach such as DALI. This work shows that the discovery of topological relations of secondary structure elements in protein structures by using spatial relations of spatial databases is practical for fast structural similarity search for proteins. PMID:16187404
Miyazaki, Jun; Iida, Tadatsune; Tanaka, Shinji; Hayashi-Takagi, Akiko; Kasai, Haruo; Okabe, Shigeo; Kobayashi, Takayoshi
2016-05-01
A fast, high-sensitivity photothermal microscope was developed by implementing a spatially segmented balanced detection scheme into a laser scanning microscope. We confirmed a 4.9 times improvement in signal-to-noise ratio in the spatially segmented balanced detection compared with that of conventional detection. The system demonstrated simultaneous bi-modal photothermal and confocal fluorescence imaging of transgenic mouse brain tissue with a pixel dwell time of 20 μs. The fluorescence image visualized neurons expressing yellow fluorescence proteins, while the photothermal signal detected endogenous chromophores in the mouse brain, allowing 3D visualization of the distribution of various features such as blood cells and fine structures probably due to lipids. This imaging modality was constructed using compact and cost-effective laser diodes, and will thus be widely useful in the life and medical sciences. PMID:27231615
Fast-response switchable lens for 3D and wearable displays.
Lee, Yun-Han; Peng, Fenglin; Wu, Shin-Tson
2016-01-25
We report a switchable lens in which a twisted nematic (TN) liquid crystal cell is utilized to control the input polarization. Different polarization state leads to different path length in the proposed optical system, which in turn results in different focal length. This type of switchable lens has advantages in fast response time, low operation voltage, and inherently lower chromatic aberration. Using a pixelated TN panel, we can create depth information to the selected pixels and thus add depth information to a 2D image. By cascading three such device structures together, we can generate 8 different focuses for 3D displays, wearable virtual/augmented reality, and other head mounted display devices. PMID:26832545
Fast phase-added stereogram algorithm for generation of photorealistic 3D content.
Kang, Hoonjong; Stoykova, Elena; Yoshikawa, Hiroshi
2016-01-20
A new phase-added stereogram algorithm for accelerated computation of holograms from a point cloud model is proposed. The algorithm relies on the hologram segmentation, sampling of directional information, and usage of the fast Fourier transform with a finer grid in the spatial frequency domain than is provided by the segment size. The algorithm gives improved quality of reconstruction due to new phase compensation introduced in the segment fringe patterns. The result is finer beam steering leading to high peak intensity and a large peak signal-to-noise ratio in reconstruction. The feasibility of the algorithm is checked by the generation of 3D contents for a color wavefront printer. PMID:26835945
Miyazaki, Jun; Iida, Tadatsune; Tanaka, Shinji; Hayashi-Takagi, Akiko; Kasai, Haruo; Okabe, Shigeo; Kobayashi, Takayoshi
2016-01-01
A fast, high-sensitivity photothermal microscope was developed by implementing a spatially segmented balanced detection scheme into a laser scanning microscope. We confirmed a 4.9 times improvement in signal-to-noise ratio in the spatially segmented balanced detection compared with that of conventional detection. The system demonstrated simultaneous bi-modal photothermal and confocal fluorescence imaging of transgenic mouse brain tissue with a pixel dwell time of 20 μs. The fluorescence image visualized neurons expressing yellow fluorescence proteins, while the photothermal signal detected endogenous chromophores in the mouse brain, allowing 3D visualization of the distribution of various features such as blood cells and fine structures probably due to lipids. This imaging modality was constructed using compact and cost-effective laser diodes, and will thus be widely useful in the life and medical sciences. PMID:27231615
Fall, Mandiaye; Boutami, Salim; Glière, Alain; Stout, Brian; Hazart, Jerome
2013-06-01
A combination of the multilevel fast multipole method (MLFMM) and boundary element method (BEM) can solve large scale photonics problems of arbitrary geometry. Here, MLFMM-BEM algorithm based on a scalar and vector potential formulation, instead of the more conventional electric and magnetic field formulations, is described. The method can deal with multiple lossy or lossless dielectric objects of arbitrary geometry, be they nested, in contact, or dispersed. Several examples are used to demonstrate that this method is able to efficiently handle 3D photonic scatterers involving large numbers of unknowns. Absorption, scattering, and extinction efficiencies of gold nanoparticle spheres, calculated by the MLFMM, are compared with Mie's theory. MLFMM calculations of the bistatic radar cross section (RCS) of a gold sphere near the plasmon resonance and of a silica coated gold sphere are also compared with Mie theory predictions. Finally, the bistatic RCS of a nanoparticle gold-silver heterodimer calculated with MLFMM is compared with unmodified BEM calculations. PMID:24323115
ALE3D Simulation of Heating and Violence in a Fast Cookoff Experiment with LX-10
McClelland, M A; Maienschein, J L; Howard, W M; Nichols, A L; deHaven, M R; Strand, O T
2006-06-26
We performed a computational and experimental analysis of fast cookoff of LX-10 (94.7% HMX, 5.3% Viton A) confined in a 2 kbar steel tube with reinforced end caps. A Scaled-Thermal-Explosion-eXperiment (STEX) was completed in which three radiant heaters were used to heat the vessel until ignition, resulting in a moderately violent explosion after 20.4 minutes. Thermocouple measurements showed tube temperatures as high as 340 C at ignition and LX-10 surface temperatures as high as 279 C, which is near the melting point of HMX. Three micro-power radar systems were used to measure mean fragment velocities of 840 m/s. Photonics Doppler Velocimeters (PDVs) showed a rapid acceleration of fragments over 80 {micro}s. A one-dimensional ALE3D cookoff model at the vessel midplane was used to simulate the heating, thermal expansion, LX-10 decomposition composition, and closing of the gap between the HE (High Explosive) and vessel wall. Although the ALE3D simulation terminated before ignition, the model provided a good representation of heat transfer through the case and across the dynamic gap to the explosive.
ALE3D Simulation and Measurement of Violence in a Fast Cookoff Experiment with LX-10
McClelland, M A; Maienschein, J L; Howard, W M; deHaven, M R
2006-11-22
We performed a computational and experimental analysis of fast cookoff of LX-10 (94.7% HMX, 5.3% Viton A) confined in a 2 kbar steel tube with reinforced end caps. A Scaled-Thermal-Explosion-eXperiment (STEX) was completed in which three radiant heaters were used to heat the vessel until ignition, resulting in a moderately violent explosion after 20.4 minutes. Thermocouple measurements showed tube temperatures as high as 340 C at ignition and LX-10 surface temperatures as high as 279 C, which is near the melting point of HMX. Three micro-power radar systems were used to measure mean fragment velocities of 840 m/s. Photonics Doppler Velocimeters (PDVs) showed a rapid acceleration of fragments over 80 {micro}s. A one-dimensional ALE3D cookoff model at the vessel midplane was used to simulate the heating, thermal expansion, LX-10 decomposition composition, and closing of the gap between the HE (High Explosive) and vessel wall. Although the ALE3D simulation terminated before ignition, the model provided a good representation of heat transfer through the case and across the dynamic gap to the explosive.
NASA Astrophysics Data System (ADS)
Wang, Zhongren; Li, Bo; Zhou, Yuebin
2014-11-01
Fast 3D reconstruction of tool wear from 2D images has great importance to 3D measuring and objective evaluating tool wear condition, determining accurate tool change and insuring machined part's quality. Extracting 3D information of tool wear zone based on monocular multi-color structured light can realize fast recovery of surface topography of tool wear, which overcomes the problems of traditional methods such as solution diversity and slow convergence when using SFS method and stereo match when using 3D reconstruction from multiple images. In this paper, a kind of new multi-color structured light illuminator was put forward. An information mapping model was established among illuminator's structure parameters, surface morphology and color images. The mathematical model to reconstruct 3D morphology based on monocular multi-color structured light was presented. Experimental results show that this method is effective and efficient to reconstruct the surface morphology of tool wear zone.
ALE3D Simulation and Measurement of Violence in a Fast Cookoff Experiment for LX-10
McClelland, M A; Maienschein, J L; Howard, W M; deHaven, M R
2006-05-23
Fast cookoff is of interest in the areas of fire hazard reduction and the development of directed energy systems for defense. During a fast cookoff (thermal explosion), high heat fluxes cause rapid temperature increases and ignition in thin boundary layers. We are developing ALE3D models to describe the thermal, chemical, and mechanical behavior during the heating, ignition, and explosive phases. The candidate models and numerical strategies are being evaluated using benchmark cookoff experiments. Fast cookoff measurements were made in a Scaled-Thermal-Explosion-eXperiment (STEX) for LX-10 (94.7% HMX, 5.3% Viton A) confined in a 4130 steel tube with reinforced end caps. Gaps were present at the side and top of the explosive charge to allow for thermal expansion. The explosive was heated until explosion using radiant heaters. Temperatures were measured using thermocouples positioned on the tube wall and in the explosive. During the explosion, the tube expansion and fragment velocities were measured with strain gauges, Photonic-Doppler-Velocimeters (PDVs), and micropower radar units. A fragment size distribution was constructed from fragments captured in Lexan panels. ALE3D models for chemical, thermal, and mechanical behavior were developed for the heating and explosive processes. A multi-step chemical kinetics model is employed for the HMX while a one-step model is used for the Viton. A pressure-dependent deflagration model is employed during the expansion. A Steinberg-Guinan model represents the mechanical behavior of the solid constituents while polynomial and gamma-law expressions are used for the equation of state of the solid and gas species, respectively. Parameters for the kinetics model were specified using measurements of the One-Dimensional-Time-to-Explosion (ODTX), while measurements for burn rate were employed to determine parameters in the burn front model. The simulations include radiative and conductive transport across the dynamic gaps between the
Fast 3D Net Expeditions: Tools for Effective Scientific Collaboration on the World Wide Web
NASA Technical Reports Server (NTRS)
Watson, Val; Chancellor, Marisa K. (Technical Monitor)
1996-01-01
Two new technologies, the FASTexpedition and Remote FAST, have been developed that provide remote, 3D (three dimensional), high resolution, dynamic, interactive viewing of scientific data. The FASTexpedition permits one to access scientific data from the World Wide Web, take guided expeditions through the data, and continue with self controlled expeditions through the data. Remote FAST permits collaborators at remote sites to simultaneously view an analysis of scientific data being controlled by one of the collaborators. Control can be transferred between sites. These technologies are now being used for remote collaboration in joint university, industry, and NASA projects. Also, NASA Ames Research Center has initiated a project to make scientific data and guided expeditions through the data available as FASTexpeditions on the World Wide Web for educational purposes. Previously, remote visualization of dynamic data was done using video format (transmitting pixel information) such as video conferencing or MPEG (Motion Picture Expert Group) movies on the Internet. The concept for this new technology is to send the raw data (e.g., grids, vectors, and scalars) along with viewing scripts over the Internet and have the pixels generated by a visualization tool running on the viewers local workstation. The visualization tool that is currently used is FAST (Flow Analysis Software Toolkit). The advantages of this new technology over using video format are: (1) The visual is much higher in resolution (1280x1024 pixels with 24 bits of color) than typical video format transmitted over the network. (2) The form of the visualization can be controlled interactively (because the viewer is interactively controlling the visualization tool running on his workstation). (3) A rich variety of guided expeditions through the data can be included easily. (4) A capability is provided for other sites to see a visual analysis of one site as the analysis is interactively performed. Control of
3-D Ground Displacement Monitoring of very fast-moving Landslides in Emergency Scenario
NASA Astrophysics Data System (ADS)
Casu, Francesco; Manconi, Andrea; Bonano, Manuela; De Luca, Claudio; Elefante, Stefano
2014-05-01
On December 3rd, 2013, a large and fast-moving landslide phenomena, which occurred in South-West of Montescaglioso town (southern Italy) after some days of intense raining, caused ground displacements on the order of several meters. The mass wasting involved an important freeway connection disrupting more than 500 meters of the route and some isolated buildings. In this work we present a case study of application of SAR remote sensing techniques for retrieving ground displacement field in a landslide emergency scenario. To this aim, thanks to the availability of ascending and descending COSMO-SkyMed (CSK) satellite acquisitions, we first applied the DInSAR technique (Massonnet et al., 1993) to both datasets, for generating differential interferograms across the investigated event. In particular, two data pairs (one ascending and one descending) involving pre- and post-event epochs and approximately spanning the same time interval were identified. Unfortunately, the DInSAR analysis produced unsatisfactory results, because of the excessive phase noise within the area of interest, mainly related to the fast-moving deformation pattern (several meters) and also to the presence of vegetation. To overcome the above mentioned limitations, the amplitude-based Pixel Offset (PO) technique (Fialko and Simons, 2001) was applied to the previous identified CSK data pairs. In this case, the PO technique allowed us to retrieve the projection of the surface displacements across and along the satellite's track (range and azimuth, respectively) for both the ascending and descending orbits. Then, by properly combining these 2-D maps of the measured surface movements, we also retrieved the 3-D ground deformation pattern, i.e. the North, East and Vertical displacement components. The ground displacements have a main SSE component, with values exceeding 10 meters. Moreover, large subsidence values were identified in those areas experiencing the largest damages, as well as a clear uplift
NASA Astrophysics Data System (ADS)
Chai, Kil-Byoung; Zhai, Xiang; Bellan, Paul M.
2016-03-01
A spatially localized energetic extreme ultra-violet (EUV) burst is imaged at the presumed position of fast magnetic reconnection in a plasma jet produced by a coaxial helicity injection source; this EUV burst indicates strong localized electron heating. A circularly polarized high frequency magnetic field perturbation is simultaneously observed at some distance from the reconnection region indicating that the reconnection emits whistler waves and that Hall dynamics likely governs the reconnection. Spectroscopic measurement shows simultaneous fast ion heating. The electron heating is consistent with Ohmic dissipation, while the ion heating is consistent with ion trajectories becoming stochastic.
Does the resistive tearing instability nonlinearly evolve to a fast reconnection mode
NASA Technical Reports Server (NTRS)
Steinolfson, R. S.
1988-01-01
A fundamental problem in applying linear tearing instability theory to the rapid processes (particle acceleration, heating) in flares was the characteristically slow rate of reconnection. This problem can be at least partially overcome if the tearin mode nonlinearly evolves to a regime in which the reconnection rate is substantially enhanced, such as that for the Petschek configuration. This possibility was often suggested, and some numerical simulations appear to provide support for such a view. Numerical simulation are used to study the nonlinear evolution of the tearing stability and show that a fast Petschek-like regime may not be achieved. This conclusion follows when there are sufficient grid points within the diffusion region to completely resolve the nonlinear dynamic interactions in the diffusion layer. When the numerical resolution is not adequate, the solution does appear to approach a Petschek configuration. The resolved solution contains reverse flow vortices and current sheets, terminated with a current reversal, similar to those obtained by Syrovatsky (JEPT, 33, 933, 1971).
Fast 3D inversion of airborne gravity-gradiometry data using Lanczos bidiagonalization method
NASA Astrophysics Data System (ADS)
Meng, Zhaohai; Li, Fengting; Zhang, Dailei; Xu, Xuechun; Huang, Danian
2016-09-01
We developed a new fast inversion method for to process and interpret airborne gravity gradiometry data, which was based on Lanczos bidiagonalization algorithm. Here, we describe the application of this new 3D gravity gradiometry inversion method to recover a subsurface density distribution model from the airborne measured gravity gradiometry anomalies. For this purpose, the survey area is divided into a large number of rectangular cells with each cell possessing a constant unknown density. It is well known that the solution of large linear gravity gradiometry is an ill-posed problem since using the smoothest inversion method is considerably time consuming. We demonstrate that the Lanczos bidiagonalization method can be an appropriate algorithm to solve a Tikhonov solver time cost function for resolving the large equations within a short time. Lanczos bidiagonalization is designed to make the very large gravity gradiometry forward modeling matrices to become low-rank, which will considerably reduce the running time of the inversion method. We also use a weighted generalized cross validation method to choose the appropriate Tikhonov parameter to improve inversion results. The inversion incorporates a model norm that allows us to attain the smoothing and depth of the solution; in addition, the model norm counteracts the natural decay of the kernels, which concentrate at shallow depths. The method is applied on noise-contaminated synthetic gravity gradiometry data to demonstrate its suitability for large 3D gravity gradiometry data inversion. The airborne gravity gradiometry data from the Vinton Salt Dome, USE, were considered as a case study. The validity of the new method on real data is discussed with reference to the Vinton Dome inversion result. The intermediate density values in the constructed model coincide well with previous results and geological information. This demonstrates the validity of the gravity gradiometry inversion method.
Explosive turbulent magnetic reconnection.
Higashimori, K; Yokoi, N; Hoshino, M
2013-06-21
We report simulation results for turbulent magnetic reconnection obtained using a newly developed Reynolds-averaged magnetohydrodynamics model. We find that the initial Harris current sheet develops in three ways, depending on the strength of turbulence: laminar reconnection, turbulent reconnection, and turbulent diffusion. The turbulent reconnection explosively converts the magnetic field energy into both kinetic and thermal energy of plasmas, and generates open fast reconnection jets. This fast turbulent reconnection is achieved by the localization of turbulent diffusion. Additionally, localized structure forms through the interaction of the mean field and turbulence. PMID:23829741
Numerical 3D models support two distinct hydrothermal circulation systems at fast spreading ridges
NASA Astrophysics Data System (ADS)
Hasenclever, Jörg; Theissen-Krah, Sonja; Rüpke, Lars
2013-04-01
We present 3D numerical calculations of hydrothermal fluid flow at fast spreading ridges. The setup of the 3D models is based our previous 2D studies, in which we have coupled numerical models for crustal accretion and hydrothermal fluid flow. One result of these calculations is a crustal permeability field that leads to a thermal structure in the crust that matches seismic tomography data of the East Pacific Rise (EPR). The 1000°C isotherm obtained from the 2D results is now used as the lower boundary of the 3D model domain, while the upper boundary is a smoothed bathymetry of the EPR. The same permeability field as in the 2D models is used, with the highest permeability at the ridge axis and a decrease with both depth and distance to the ridge. Permeability is also reduced linearly between 600 and 1000°C. Using a newly developed parallel finite element code written in Matlab that solves for thermal evolution, fluid pressure and Darcy flow, we simulate the flow patterns of hydrothermal circulation in a segment of 5000m along-axis, 10000m across-axis and up to 5000m depth. We observe two distinct hydrothermal circulation systems: An on-axis system forming a series of vents with a spacing ranging from 100 to 500m that is recharged by nearby (100-200m) downflows on both sides of the ridge axis. Simultaneously a second system with much broader extensions both laterally and vertically exists off-axis. It is recharged by fluids intruding between 1500m to 5000m off-axis and sampling both upper and lower crust. These fluids are channeled in the deepest and hottest regions with high permeability and migrate up-slope following the 600°C isotherm until reaching the edge of the melt lens. Depending on the width of the melt lens these off-axis fluids either merge with the on-axis hydrothermal system or form separate vents. We observe separate off-axis vent fields if the magma lens half-width exceeds 1000m and confluence of both systems for half-widths smaller than 500m. For
Fast Semantic Segmentation of 3d Point Clouds with Strongly Varying Density
NASA Astrophysics Data System (ADS)
Hackel, Timo; Wegner, Jan D.; Schindler, Konrad
2016-06-01
We describe an effective and efficient method for point-wise semantic classification of 3D point clouds. The method can handle unstructured and inhomogeneous point clouds such as those derived from static terrestrial LiDAR or photogammetric reconstruction; and it is computationally efficient, making it possible to process point clouds with many millions of points in a matter of minutes. The key issue, both to cope with strong variations in point density and to bring down computation time, turns out to be careful handling of neighborhood relations. By choosing appropriate definitions of a point's (multi-scale) neighborhood, we obtain a feature set that is both expressive and fast to compute. We evaluate our classification method both on benchmark data from a mobile mapping platform and on a variety of large, terrestrial laser scans with greatly varying point density. The proposed feature set outperforms the state of the art with respect to per-point classification accuracy, while at the same time being much faster to compute.
NASA Astrophysics Data System (ADS)
Cheng, Chio Z.; Ono, Yasushi; Yang, Ya-Hui; Choe, Gwangson
2012-10-01
Impulsively fast magnetic reconnection has been shown to be the major mechanism responsible for explosive flare non-thermal energy release and acceleration of coronal mass ejection (CME) motion. It has been observed that for most large solar flares non-thermal emissions in hard X-rays (HXR) and millimeter/submillimeter waves impulsively rise and decade during the soft X-ray (SXR) emission rise phase. Moreover, the filament/CME upward motion is accelerated temporally in correlation with the impulsive enhancement of flare non-thermal emission and reconnection electric field in the current sheet in both simulations and observations. The peak reconnection electric field during flare impulsive phase is on the order of a few kV/m for X-class flares. Here, we demonstrated for the first time in laboratory plasma merging experiments the correlation of the magnetic reconnection rate with the acceleration of plasmoid ejected from the current sheet using the TS-4 device of the Tokyo University. Moreover, we have also found that the electron heating occurs in the current sheet and the ion heating occurs in the down-stream outflow region. Thus, we conclude that the plasmoid/CME acceleration is a key mechanism for the impulsive enhancement of magnetic reconnection rate (electric field).
Ultra-fast hybrid CPU-GPU multiple scatter simulation for 3-D PET.
Kim, Kyung Sang; Son, Young Don; Cho, Zang Hee; Ra, Jong Beom; Ye, Jong Chul
2014-01-01
Scatter correction is very important in 3-D PET reconstruction due to a large scatter contribution in measurements. Currently, one of the most popular methods is the so-called single scatter simulation (SSS), which considers single Compton scattering contributions from many randomly distributed scatter points. The SSS enables a fast calculation of scattering with a relatively high accuracy; however, the accuracy of SSS is dependent on the accuracy of tail fitting to find a correct scaling factor, which is often difficult in low photon count measurements. To overcome this drawback as well as to improve accuracy of scatter estimation by incorporating multiple scattering contribution, we propose a multiple scatter simulation (MSS) based on a simplified Monte Carlo (MC) simulation that considers photon migration and interactions due to photoelectric absorption and Compton scattering. Unlike the SSS, the MSS calculates a scaling factor by comparing simulated prompt data with the measured data in the whole volume, which enables a more robust estimation of a scaling factor. Even though the proposed MSS is based on MC, a significant acceleration of the computational time is possible by using a virtual detector array with a larger pitch by exploiting that the scatter distribution varies slowly in spatial domain. Furthermore, our MSS implementation is nicely fit to a parallel implementation using graphic processor unit (GPU). In particular, we exploit a hybrid CPU-GPU technique using the open multiprocessing and the compute unified device architecture, which results in 128.3 times faster than using a single CPU. Overall, the computational time of MSS is 9.4 s for a high-resolution research tomograph (HRRT) system. The performance of the proposed MSS is validated through actual experiments using an HRRT. PMID:24403412
NASA Astrophysics Data System (ADS)
Pucci, Fulvia; Del Sarto, Daniele; Tenerani, Anna; Velli, Marco
2015-04-01
By examining sheets with thicknesses scaling as different powers of the Lundquist number S, we previously showed (Pucci and Velli, 2014) that the growth rate of the tearing mode increases as current sheets thin and, once the inverse aspect ratio reaches a scaling a/L = S-1/3, the time-scale for the instability to develop becomes of the order of the Alfvén time. That means that a fast instability sets in well before Sweet-Parker type current sheets can form. In addition, such an instability produces many islands in the sheet, leading to fast nonlinear evolution and most probably a turbulent disruption of the sheet itself. This has fundamental implications for magnetically driven reconnection throughout the corona, and in particular for coronal heating and the triggering of coronal mass ejections. Here we extend the study of reconnection instabilities to magnetic fields of grater complexity, displaying different current structures such as, for example, multiple or asymmetric current layers. We also consider the possibility of a Δ' dependence on wave-number k-p for different values of p, studying analogies and variations of the trigger scaling relation a/L ~ S-1/3 with respect to the Harris current sheet equilibrium. At large Lundquist numbers in typical Heliospheric plasmas kinetic effects become more important in Ohm's law: we consider the effects of electron skin depth reconnection, showing that we can define a trigger relation similar to the resistive case. The results are important to the transition to fast reconnection in the solar corona, solar wind, magnetosphere as well as laboratory plasmas. F. Pucci and M. Velli, "Reconnection of quasi-singular current sheets: the 'ideal" tearing mode" ApJ 780:L19, 2014.
Fast Susceptibility-Weighted Imaging (SWI) with 3D Short-Axis Propeller (SAP)-EPI
Holdsworth, Samantha J.; Yeom, Kristen W.; Moseley, Michael E.; Skare, S.
2014-01-01
Purpose Susceptibility-Weighted Imaging (SWI) in neuroimaging can be challenging due to long scan times of 3D Gradient Recalled Echo (GRE), while faster techniques such as 3D interleaved EPI (iEPI) are prone to motion artifacts. Here we outline and implement a 3D Short-Axis Propeller Echo-Planar Imaging (SAP-EPI) trajectory as a faster, motion-correctable approach for SWI. Methods Experiments were conducted on a 3T MRI system. 3D SAP-EPI, 3D iEPI, and 3D GRE SWI scans were acquired on two volunteers. Controlled motion experiments were conducted to test the motion-correction capability of 3D SAP-EPI. 3D SAP-EPI SWI data were acquired on two pediatric patients as a potential alternative to 2D GRE used clinically. Results 3D GRE images had a better target resolution (0.47 × 0.94 × 2mm, scan time = 5min), iEPI and SAP-EPI images (resolution = 0.94 × 0.94 × 2mm) were acquired in a faster scan time (1:52min) with twice the brain coverage. SAP-EPI showed motion-correction capability and some immunity to undersampling from rejected data. Conclusion While 3D SAP-EPI suffers from some geometric distortion, its short scan time and motion-correction capability suggest that SAP-EPI may be a useful alternative to GRE and iEPI for use in SWI, particularly in uncooperative patients. PMID:24956237
Dynamic Response of Magnetic Reconnection Due to Current Sheet Variability
NASA Astrophysics Data System (ADS)
George, D. E.; Jahn, J. M.; Burch, J. L.; Hesse, M.; Pollock, C. J.
2014-12-01
Magnetic reconnection is a process which regulates the interaction between regions of magnetized plasma. While many factors have an impact on the evolution of this process, there still remains a lack of understanding of the key behaviors involved in the triggering of fast reconnection. Despite an abundance of in-situ measurements, indicating the high degree of variability in the thickness, density and composition along the current sheet, no simulation studies exist which account for such current sheet variations. 2D and 3D simulations have a periodic boundary in the dimension along the current sheet and so tend to neglect these variations in the current sheet originating external to the modeled reconnection region. Here we focus on the effects on reconnection due to the variability in the thickness and density of the current sheet. Using 2.5D kinetic simulations of 2-species plasma, we isolate and explore the dynamic effects on reconnection associated with variations in the current sheet originating externally to the reconnection region. While periodic boundary conditions are still used, in the direction along the current sheet, a step-change perturbation in thickness or density of the current sheet is introduced once a stable reconnection rate is reached. The dynamic response of the overall system, after introducing the perturbation, is then evaluated, with a focus on the reconnection rate. When the reconnection rate is slowed significantly over time, loading of the inflow region occurs (a build-up of plasma and magnetic energy/pressure. This state is indicated by an asymptotic behavior in the reconnection rate over time. If a sudden variation in the current sheet is introduced under these conditions, a resultant triggering of fast reconnection may occur, which could lead to an episode of fast reconnection, saw-tooth-crash condition or even act as a trigger for sub-storms.
Fast and Precise 3D Fluorophore Localization based on Gradient Fitting
NASA Astrophysics Data System (ADS)
Ma, Hongqiang; Xu, Jianquan; Jin, Jingyi; Gao, Ying; Lan, Li; Liu, Yang
2015-09-01
Astigmatism imaging approach has been widely used to encode the fluorophore’s 3D position in single-particle tracking and super-resolution localization microscopy. Here, we present a new high-speed localization algorithm based on gradient fitting to precisely decode the 3D subpixel position of the fluorophore. This algebraic algorithm determines the center of the fluorescent emitter by finding the position with the best-fit gradient direction distribution to the measured point spread function (PSF), and can retrieve the 3D subpixel position of the fluorophore in a single iteration. Through numerical simulation and experiments with mammalian cells, we demonstrate that our algorithm yields comparable localization precision to the traditional iterative Gaussian function fitting (GF) based method, while exhibits over two orders-of-magnitude faster execution speed. Our algorithm is a promising high-speed analyzing method for 3D particle tracking and super-resolution localization microscopy.
Fast and Precise 3D Fluorophore Localization based on Gradient Fitting
Ma, Hongqiang; Xu, Jianquan; Jin, Jingyi; Gao, Ying; Lan, Li; Liu, Yang
2015-01-01
Astigmatism imaging approach has been widely used to encode the fluorophore’s 3D position in single-particle tracking and super-resolution localization microscopy. Here, we present a new high-speed localization algorithm based on gradient fitting to precisely decode the 3D subpixel position of the fluorophore. This algebraic algorithm determines the center of the fluorescent emitter by finding the position with the best-fit gradient direction distribution to the measured point spread function (PSF), and can retrieve the 3D subpixel position of the fluorophore in a single iteration. Through numerical simulation and experiments with mammalian cells, we demonstrate that our algorithm yields comparable localization precision to the traditional iterative Gaussian function fitting (GF) based method, while exhibits over two orders-of-magnitude faster execution speed. Our algorithm is a promising high-speed analyzing method for 3D particle tracking and super-resolution localization microscopy. PMID:26390959
Self-generated Turbulence in Magnetic Reconnection
NASA Astrophysics Data System (ADS)
Oishi, Jeffrey S.; Mac Low, Mordecai-Mark; Collins, David C.; Tamura, Moeko
2015-06-01
Classical Sweet-Parker models of reconnection predict that reconnection rates depend inversely on the resistivity, usually parameterized using the dimensionless Lundquist number (S). We describe magnetohydrodynamic (MHD) simulations using a static, nested grid that show the development of a three-dimensional (3D) instability in the plane of a current sheet between reversing field lines without a guide field. The instability leads to rapid reconnection of magnetic field lines at a rate independent of S over at least the range 3.2× {{10}3}≲ S≲ 3.2× {{10}5} resolved by the simulations. We find that this instability occurs even for cases with S≲ {{10}4} that in our models appear stable to the recently described, two-dimensional, plasmoid instability. Our results suggest that 3D, MHD processes alone produce fast (resistivity independent) reconnection without recourse to kinetic effects or external turbulence. The unstable reconnection layers provide a self-consistent environment in which the extensively studied turbulent reconnection process can occur.
Three-Dimensional Localized Magnetic Reconnection in Torus Plasma Merging Device TS-4
NASA Astrophysics Data System (ADS)
Ii, Toru; Hayashi, Yoshinori; Inomoto, Michiaki; Ono, Yasushi
Three-dimensional (3-D) localized magnetic reconnection was studied experimentally using torus plasma merging device TS-4. The direct measurements of 3-D structures of current sheet revealed two unsteady and fast reconnection mechanisms: 3-D deformation of current sheet and mass ejection. When strong compression force IAcc∼60kA was applied to two plasma toroids with low guide field Bt/B¦¦∼1, toroidal modes n=1-3 of current sheet were observed to grow only during their reconnection and to disappear after the reconnection. This 3-D deformation promoted mass ejection from the current sheet, increasing the reconnection rate as well as reconnection (toroidal) electric field and outflow. On the other hand, the reconnection rate was maintained low under the high guide field Bt/B¦¦∼7 and weak compression IAcc∼0kA. These phenomena suggest that local compression of current sheet triggers its strong dissipation as well as plasma mass ejection, which are responsible for the onset of 3-D localized reconnection.
Song, Liang; Maslov, Konstantin; Bitton, Rachel; Shung, K. Kirk; Wang, Lihong V.
2009-01-01
We present an in vivo dark-field reflection-mode photoacoustic microscopy system that performs cross-sectional (B-scan) imaging at 50 Hz with realtime beamforming and 3D imaging consisting of 166 B-scan frames at 1 Hz with post-beamforming. To our knowledge, this speed is currently the fastest in photoacoustic imaging. A custom-designed light delivery system is integrated with a 30-MHz ultrasound linear array to realize dark-field reflection-mode imaging. Linear mechanical scanning of the array produces 3D images. The system has axial, lateral, and elevational resolutions of 25, 70, and 200 μm, respectively, and can image 3 mm deep in scattering biological tissues. Volumetric images of subcutaneous vasculature in rats are demonstrated in vivo. Fast 3D photoacoustic microscopy is anticipated to facilitate applications of photoacoustic imaging in biomedical studies that involve dynamics and clinical procedures that demand immediate diagnosis. PMID:19021408
Boulanger, Jérôme; Gueudry, Charles; Münch, Daniel; Cinquin, Bertrand; Paul-Gilloteaux, Perrine; Bardin, Sabine; Guérin, Christophe; Senger, Fabrice; Blanchoin, Laurent; Salamero, Jean
2014-01-01
Total internal reflection fluorescence microscopy (TIRFM) is the method of choice to visualize a variety of cellular processes in particular events localized near the plasma membrane of live adherent cells. This imaging technique not relying on particular fluorescent probes provides a high sectioning capability. It is, however, restricted to a single plane. We present here a method based on a versatile design enabling fast multiwavelength azimuthal averaging and incidence angles scanning to computationally reconstruct 3D images sequences. We achieve unprecedented 50-nm axial resolution over a range of 800 nm above the coverslip. We apply this imaging modality to obtain structural and dynamical information about 3D actin architectures. We also temporally decipher distinct Rab11a-dependent exocytosis events in 3D at a rate of seven stacks per second. PMID:25404337
Masaaki Yamada, Russell Kulsrud and Hantao Ji
2009-09-17
We review the fundamental physics of magnetic reconnection in laboratory and space plasmas, by discussing results from theory, numerical simulations, observations from space satellites, and the recent results from laboratory plasma experiments. After a brief review of the well-known early work, we discuss representative recent experimental and theoretical work and attempt to interpret the essence of significant modern findings. In the area of local reconnection physics, many significant findings have been made with regard to two- uid physics and are related to the cause of fast reconnection. Profiles of the neutral sheet, Hall currents, and the effects of guide field, collisions, and micro-turbulence are discussed to understand the fundamental processes in a local reconnection layer both in space and laboratory plasmas. While the understanding of the global reconnection dynamics is less developed, notable findings have been made on this issue through detailed documentation of magnetic self-organization phenomena in fusion plasmas. Application of magnetic reconnection physics to astrophysical plasmas is also brie y discussed.
Courau, T.; Moustafa, S.; Plagne, L.; Poncot, A.
2013-07-01
As part of its activity, EDF R and D is developing a new nuclear core simulation code named COCAGNE. This code relies on DIABOLO, a Simplified PN (SPN) method to compute the neutron flux inside the core for eigenvalue calculations. In order to assess the accuracy of SPN calculations, we have developed DOMINO, a new 3D Cartesian SN solver. The parallel implementation of DOMINO is very efficient and allows to complete an eigenvalue calculation involving around 300 x 10{sup 9} degrees of freedom within a few hours on a single shared-memory supercomputing node. This computation corresponds to a 26-group S{sub 8} 3D PWR core model used to assess the SPN accuracy. At the pin level, the maximal error for the SP{sub 5} DIABOLO fission production rate is lower than 0.2% compared to the S{sub 8} DOMINO reference for this 3D PWR core model. (authors)
Bettadapura, Radhakrishna; Rasheed, Muhibur; Vollrath, Antje; Bajaj, Chandrajit
2015-10-01
There continue to be increasing occurrences of both atomistic structure models in the PDB (possibly reconstructed from X-ray diffraction or NMR data), and 3D reconstructed cryo-electron microscopy (3D EM) maps (albeit at coarser resolution) of the same or homologous molecule or molecular assembly, deposited in the EMDB. To obtain the best possible structural model of the molecule at the best achievable resolution, and without any missing gaps, one typically aligns (match and fits) the atomistic structure model with the 3D EM map. We discuss a new algorithm and generalized framework, named PF(2) fit (Polar Fast Fourier Fitting) for the best possible structural alignment of atomistic structures with 3D EM. While PF(2) fit enables only a rigid, six dimensional (6D) alignment method, it augments prior work on 6D X-ray structure and 3D EM alignment in multiple ways: Scoring. PF(2) fit includes a new scoring scheme that, in addition to rewarding overlaps between the volumes occupied by the atomistic structure and 3D EM map, rewards overlaps between the volumes complementary to them. We quantitatively demonstrate how this new complementary scoring scheme improves upon existing approaches. PF(2) fit also includes two scoring functions, the non-uniform exterior penalty and the skeleton-secondary structure score, and implements the scattering potential score as an alternative to traditional Gaussian blurring. Search. PF(2) fit utilizes a fast polar Fourier search scheme, whose main advantage is the ability to search over uniformly and adaptively sampled subsets of the space of rigid-body motions. PF(2) fit also implements a new reranking search and scoring methodology that considerably improves alignment metrics in results obtained from the initial search. PMID:26469938
Bettadapura, Radhakrishna; Rasheed, Muhibur; Vollrath, Antje; Bajaj, Chandrajit
2015-01-01
There continue to be increasing occurrences of both atomistic structure models in the PDB (possibly reconstructed from X-ray diffraction or NMR data), and 3D reconstructed cryo-electron microscopy (3D EM) maps (albeit at coarser resolution) of the same or homologous molecule or molecular assembly, deposited in the EMDB. To obtain the best possible structural model of the molecule at the best achievable resolution, and without any missing gaps, one typically aligns (match and fits) the atomistic structure model with the 3D EM map. We discuss a new algorithm and generalized framework, named PF2 fit (Polar Fast Fourier Fitting) for the best possible structural alignment of atomistic structures with 3D EM. While PF2 fit enables only a rigid, six dimensional (6D) alignment method, it augments prior work on 6D X-ray structure and 3D EM alignment in multiple ways: Scoring. PF2 fit includes a new scoring scheme that, in addition to rewarding overlaps between the volumes occupied by the atomistic structure and 3D EM map, rewards overlaps between the volumes complementary to them. We quantitatively demonstrate how this new complementary scoring scheme improves upon existing approaches. PF2 fit also includes two scoring functions, the non-uniform exterior penalty and the skeleton-secondary structure score, and implements the scattering potential score as an alternative to traditional Gaussian blurring. Search. PF2 fit utilizes a fast polar Fourier search scheme, whose main advantage is the ability to search over uniformly and adaptively sampled subsets of the space of rigid-body motions. PF2 fit also implements a new reranking search and scoring methodology that considerably improves alignment metrics in results obtained from the initial search. PMID:26469938
A fast rebinning algorithm for 3D positron emission tomography using John's equation
NASA Astrophysics Data System (ADS)
Defrise, Michel; Liu, Xuan
1999-08-01
Volume imaging in positron emission tomography (PET) requires the inversion of the three-dimensional (3D) x-ray transform. The usual solution to this problem is based on 3D filtered-backprojection (FBP), but is slow. Alternative methods have been proposed which factor the 3D data into independent 2D data sets corresponding to the 2D Radon transforms of a stack of parallel slices. Each slice is then reconstructed using 2D FBP. These so-called rebinning methods are numerically efficient but are approximate. In this paper a new exact rebinning method is derived by exploiting the fact that the 3D x-ray transform of a function is the solution to the second-order partial differential equation first studied by John. The method is proposed for two sampling schemes, one corresponding to a pair of infinite plane detectors and another one corresponding to a cylindrical multi-ring PET scanner. The new FORE-J algorithm has been implemented for this latter geometry and was compared with the approximate Fourier rebinning algorithm FORE and with another exact rebinning algorithm, FOREX. Results with simulated data demonstrate a significant improvement in accuracy compared to FORE, while the reconstruction time is doubled. Compared to FOREX, the FORE-J algorithm is slightly less accurate but more than three times faster.
A fast and accurate method to predict 2D and 3D aerodynamic boundary layer flows
NASA Astrophysics Data System (ADS)
Bijleveld, H. A.; Veldman, A. E. P.
2014-12-01
A quasi-simultaneous interaction method is applied to predict 2D and 3D aerodynamic flows. This method is suitable for offshore wind turbine design software as it is a very accurate and computationally reasonably cheap method. This study shows the results for a NACA 0012 airfoil. The two applied solvers converge to the experimental values when the grid is refined. We also show that in separation the eigenvalues remain positive thus avoiding the Goldstein singularity at separation. In 3D we show a flow over a dent in which separation occurs. A rotating flat plat is used to show the applicability of the method for rotating flows. The shown capabilities of the method indicate that the quasi-simultaneous interaction method is suitable for design methods for offshore wind turbine blades.
Faganello, M; Califano, F; Pegoraro, F
2008-09-01
We give evidence for the first time of the onset of undriven fast, collisionless magnetic reconnection during the evolution of an initially homogeneous magnetic field advected in a sheared velocity field. We consider the interaction of the solar wind with the magnetospheric plasma at low latitude and show that reconnection takes place in the layer between adjacent vortices generated by the Kelvin-Helmholtz instability. This process generates coherent magnetic structures with a size comparable to the ion inertial scale, much smaller than the system dimensions but much larger than the electron inertial scale. These magnetic structures are further advected in the plasma in a complex pattern but remain stable over a time interval much longer than their formation time. These results can be crucial for the interpretation of satellite data showing coherent magnetic structures in the Earth's magnetosheath or the magnetotail. PMID:18851219
The fast and accurate 3D-face scanning technology based on laser triangle sensors
NASA Astrophysics Data System (ADS)
Wang, Jinjiang; Chang, Tianyu; Ge, Baozhen; Tian, Qingguo; Chen, Yang; Kong, Bin
2013-08-01
A laser triangle scanning method and the structure of 3D-face measurement system were introduced. In presented system, a liner laser source was selected as an optical indicated signal in order to scanning a line one times. The CCD image sensor was used to capture image of the laser line modulated by human face. The system parameters were obtained by system calibrated calculated. The lens parameters of image part of were calibrated with machine visual image method and the triangle structure parameters were calibrated with fine wire paralleled arranged. The CCD image part and line laser indicator were set with a linear motor carry which can achieve the line laser scanning form top of the head to neck. For the nose is ledge part and the eyes are sunk part, one CCD image sensor can not obtain the completed image of laser line. In this system, two CCD image sensors were set symmetric at two sides of the laser indicator. In fact, this structure includes two laser triangle measure units. Another novel design is there laser indicators were arranged in order to reduce the scanning time for it is difficult for human to keep static for longer time. The 3D data were calculated after scanning. And further data processing include 3D coordinate refine, mesh calculate and surface show. Experiments show that this system has simply structure, high scanning speed and accurate. The scanning range covers the whole head of adult, the typical resolution is 0.5mm.
Fast error simulation of optical 3D measurements at translucent objects
NASA Astrophysics Data System (ADS)
Lutzke, P.; Kühmstedt, P.; Notni, G.
2012-09-01
The scan results of optical 3D measurements at translucent objects deviate from the real objects surface. This error is caused by the fact that light is scattered in the objects volume and is not exclusively reflected at its surface. A few approaches were made to separate the surface reflected light from the volume scattered. For smooth objects the surface reflected light is dominantly concentrated in specular direction and could only be observed from a point in this direction. Thus the separation either leads to measurement results only creating data for near specular directions or provides data from not well separated areas. To ensure the flexibility and precision of optical 3D measurement systems for translucent materials it is necessary to enhance the understanding of the error forming process. For this purpose a technique for simulating the 3D measurement at translucent objects is presented. A simple error model is shortly outlined and extended to an efficient simulation environment based upon ordinary raytracing methods. In comparison the results of a Monte-Carlo simulation are presented. Only a few material and object parameters are needed for the raytracing simulation approach. The attempt of in-system collection of these material and object specific parameters is illustrated. The main concept of developing an error-compensation method based on the simulation environment and the collected parameters is described. The complete procedure is using both, the surface reflected and the volume scattered light for further processing.
WaveQ3D: Fast and accurate acoustic transmission loss (TL) eigenrays, in littoral environments
NASA Astrophysics Data System (ADS)
Reilly, Sean M.
This study defines a new 3D Gaussian ray bundling acoustic transmission loss model in geodetic coordinates: latitude, longitude, and altitude. This approach is designed to lower the computation burden of computing accurate environmental effects in sonar training application by eliminating the need to transform the ocean environment into a collection of Nx2D Cartesian radials. This approach also improves model accuracy by incorporating real world 3D effects, like horizontal refraction, into the model. This study starts with derivations for a 3D variant of Gaussian ray bundles in this coordinate system. To verify the accuracy of this approach, acoustic propagation predictions of transmission loss, time of arrival, and propagation direction are compared to analytic solutions and other models. To validate the model's ability to predict real world phenomena, predictions of transmission loss and propagation direction are compared to at-sea measurements, in an environment where strong horizontal refraction effect have been observed. This model has been integrated into U.S. Navy active sonar training system applications, where testing has demonstrated its ability to improve transmission loss calculation speed without sacrificing accuracy.
Progressive attenuation fields: Fast 2D-3D image registration without precomputation
Rohlfing, Torsten; Russakoff, Daniel B.; Denzler, Joachim; Mori, Kensaku; Maurer, Calvin R. Jr.
2005-09-15
Computation of digitally reconstructed radiograph (DRR) images is the rate-limiting step in most current intensity-based algorithms for the registration of three-dimensional (3D) images to two-dimensional (2D) projection images. This paper introduces and evaluates the progressive attenuation field (PAF), which is a new method to speed up DRR computation. A PAF is closely related to an attenuation field (AF). A major difference is that a PAF is constructed on the fly as the registration proceeds; it does not require any precomputation time, nor does it make any prior assumptions of the patient pose or limit the permissible range of patient motion. A PAF effectively acts as a cache memory for projection values once they are computed, rather than as a lookup table for precomputed projections like standard AFs. We use a cylindrical attenuation field parametrization, which is better suited for many medical applications of 2D-3D registration than the usual two-plane parametrization. The computed attenuation values are stored in a hash table for time-efficient storage and access. Using clinical gold-standard spine image data sets from five patients, we demonstrate consistent speedups of intensity-based 2D-3D image registration using PAF DRRs by a factor of 10 over conventional ray casting DRRs with no decrease of registration accuracy or robustness.
Fast imaging of laboratory core floods using 3D compressed sensing RARE MRI.
Ramskill, N P; Bush, I; Sederman, A J; Mantle, M D; Benning, M; Anger, B C; Appel, M; Gladden, L F
2016-09-01
Three-dimensional (3D) imaging of the fluid distributions within the rock is essential to enable the unambiguous interpretation of core flooding data. Magnetic resonance imaging (MRI) has been widely used to image fluid saturation in rock cores; however, conventional acquisition strategies are typically too slow to capture the dynamic nature of the displacement processes that are of interest. Using Compressed Sensing (CS), it is possible to reconstruct a near-perfect image from significantly fewer measurements than was previously thought necessary, and this can result in a significant reduction in the image acquisition times. In the present study, a method using the Rapid Acquisition with Relaxation Enhancement (RARE) pulse sequence with CS to provide 3D images of the fluid saturation in rock core samples during laboratory core floods is demonstrated. An objective method using image quality metrics for the determination of the most suitable regularisation functional to be used in the CS reconstructions is reported. It is shown that for the present application, Total Variation outperforms the Haar and Daubechies3 wavelet families in terms of the agreement of their respective CS reconstructions with a fully-sampled reference image. Using the CS-RARE approach, 3D images of the fluid saturation in the rock core have been acquired in 16min. The CS-RARE technique has been applied to image the residual water saturation in the rock during a water-water displacement core flood. With a flow rate corresponding to an interstitial velocity of vi=1.89±0.03ftday(-1), 0.1 pore volumes were injected over the course of each image acquisition, a four-fold reduction when compared to a fully-sampled RARE acquisition. Finally, the 3D CS-RARE technique has been used to image the drainage of dodecane into the water-saturated rock in which the dynamics of the coalescence of discrete clusters of the non-wetting phase are clearly observed. The enhancement in the temporal resolution that has
Fast imaging of laboratory core floods using 3D compressed sensing RARE MRI
NASA Astrophysics Data System (ADS)
Ramskill, N. P.; Bush, I.; Sederman, A. J.; Mantle, M. D.; Benning, M.; Anger, B. C.; Appel, M.; Gladden, L. F.
2016-09-01
Three-dimensional (3D) imaging of the fluid distributions within the rock is essential to enable the unambiguous interpretation of core flooding data. Magnetic resonance imaging (MRI) has been widely used to image fluid saturation in rock cores; however, conventional acquisition strategies are typically too slow to capture the dynamic nature of the displacement processes that are of interest. Using Compressed Sensing (CS), it is possible to reconstruct a near-perfect image from significantly fewer measurements than was previously thought necessary, and this can result in a significant reduction in the image acquisition times. In the present study, a method using the Rapid Acquisition with Relaxation Enhancement (RARE) pulse sequence with CS to provide 3D images of the fluid saturation in rock core samples during laboratory core floods is demonstrated. An objective method using image quality metrics for the determination of the most suitable regularisation functional to be used in the CS reconstructions is reported. It is shown that for the present application, Total Variation outperforms the Haar and Daubechies3 wavelet families in terms of the agreement of their respective CS reconstructions with a fully-sampled reference image. Using the CS-RARE approach, 3D images of the fluid saturation in the rock core have been acquired in 16 min. The CS-RARE technique has been applied to image the residual water saturation in the rock during a water-water displacement core flood. With a flow rate corresponding to an interstitial velocity of vi = 1.89 ± 0.03 ft day-1, 0.1 pore volumes were injected over the course of each image acquisition, a four-fold reduction when compared to a fully-sampled RARE acquisition. Finally, the 3D CS-RARE technique has been used to image the drainage of dodecane into the water-saturated rock in which the dynamics of the coalescence of discrete clusters of the non-wetting phase are clearly observed. The enhancement in the temporal resolution
Fast magnetic reconnection in low-density electron-positron plasmas
Bessho, Naoki; Bhattacharjee, A.
2010-10-15
Two-dimensional particle-in-cell simulations have been performed to study magnetic reconnection in low-density electron-positron plasmas without a guide magnetic field. Impulsive reconnection rates become of the order of unity when the background density is much smaller than 10% of the density in the initial current layer. It is demonstrated that the outflow speed is less than the upstream Alfven speed, and that the time derivative of the density must be taken into account in the definition of the reconnection rate. The reconnection electric fields in the low-density regime become much larger than the ones in the high-density regime, and it is possible to accelerate the particles to high energies more efficiently. The inertial term in the generalized Ohm's law is the most dominant term that supports a large reconnection electric field. An effective collisionless resistivity is produced and tracks the extension of the diffusion region in the late stage of the reconnection dynamics, and significant broadening of the diffusion region is observed. Because of the broadening of the diffusion region, no secondary islands, which have been considered to play a role to limit the diffusion region, are generated during the extension of the diffusion region in the outflow direction.
Fast quantitative susceptibility mapping using 3D EPI and total generalized variation.
Langkammer, Christian; Bredies, Kristian; Poser, Benedikt A; Barth, Markus; Reishofer, Gernot; Fan, Audrey Peiwen; Bilgic, Berkin; Fazekas, Franz; Mainero, Caterina; Ropele, Stefan
2015-05-01
Quantitative susceptibility mapping (QSM) allows new insights into tissue composition and organization by assessing its magnetic property. Previous QSM studies have already demonstrated that magnetic susceptibility is highly sensitive to myelin density and fiber orientation as well as to para- and diamagnetic trace elements. Image resolution in QSM with current approaches is limited by the long acquisition time of 3D scans and the need for high signal to noise ratio (SNR) to solve the dipole inversion problem. We here propose a new total-generalized-variation (TGV) based method for QSM reconstruction, which incorporates individual steps of phase unwrapping, background field removal and dipole inversion in a single iteration, thus yielding a robust solution to the reconstruction problem. This approach has beneficial characteristics for low SNR data, allowing for phase data to be rapidly acquired with a 3D echo planar imaging (EPI) sequence. The proposed method was evaluated with a numerical phantom and in vivo at 3 and 7 T. Compared to total variation (TV), TGV-QSM enforced higher order smoothness which yielded solutions closer to the ground truth and prevented stair-casing artifacts. The acquisition time for images with 1mm isotropic resolution and whole brain coverage was 10s on a clinical 3 Tesla scanner. In conclusion, 3D EPI acquisition combined with single-step TGV reconstruction yields reliable QSM images of the entire brain with 1mm isotropic resolution in seconds. The short acquisition time combined with the robust reconstruction may enable new QSM applications in less compliant populations, clinical susceptibility tensor imaging, and functional resting state examinations. PMID:25731991
NASA Astrophysics Data System (ADS)
Chakraborty, Bidisha; Heyde, Brecht; Alessandrini, Martino; D'hooge, Jan
2016-04-01
Image registration techniques using free-form deformation models have shown promising results for 3D myocardial strain estimation from ultrasound. However, the use of this technique has mostly been limited to research institutes due to the high computational demand, which is primarily due to the computational load of the regularization term ensuring spatially smooth cardiac strain estimates. Indeed, this term typically requires evaluating derivatives of the transformation field numerically in each voxel of the image during every iteration of the optimization process. In this paper, we replace this time-consuming step with a closed-form solution directly associated with the transformation field resulting in a speed up factor of ~10-60,000, for a typical 3D B-mode image of 2503 and 5003 voxels, depending upon the size and the parametrization of the transformation field. The performance of the numeric and the analytic solutions was contrasted by computing tracking and strain accuracy on two realistic synthetic 3D cardiac ultrasound sequences, mimicking two ischemic motion patterns. Mean and standard deviation of the displacement errors over the cardiac cycle for the numeric and analytic solutions were 0.68+/-0.40 mm and 0.75+/-0.43 mm respectively. Correlations for the radial, longitudinal and circumferential strain components at end-systole were 0.89, 0.83 and 0.95 versus 0.90, 0.88 and 0.92 for the numeric and analytic regularization respectively. The analytic solution matched the performance of the numeric solution as no statistically significant differences (p>0.05) were found when expressed in terms of bias or limits-of-agreement.
Robust 2D/3D registration for fast-flexion motion of the knee joint using hybrid optimization.
Ohnishi, Takashi; Suzuki, Masahiko; Kobayashi, Tatsuya; Naomoto, Shinji; Sukegawa, Tomoyuki; Nawata, Atsushi; Haneishi, Hideaki
2013-01-01
Previously, we proposed a 2D/3D registration method that uses Powell's algorithm to obtain 3D motion of a knee joint by 3D computed-tomography and bi-plane fluoroscopic images. The 2D/3D registration is performed consecutively and automatically for each frame of the fluoroscopic images. This method starts from the optimum parameters of the previous frame for each frame except for the first one, and it searches for the next set of optimum parameters using Powell's algorithm. However, if the flexion motion of the knee joint is fast, it is likely that Powell's algorithm will provide a mismatch because the initial parameters are far from the correct ones. In this study, we applied a hybrid optimization algorithm (HPS) combining Powell's algorithm with the Nelder-Mead simplex (NM-simplex) algorithm to overcome this problem. The performance of the HPS was compared with the separate performances of Powell's algorithm and the NM-simplex algorithm, the Quasi-Newton algorithm and hybrid optimization algorithm with the Quasi-Newton and NM-simplex algorithms with five patient data sets in terms of the root-mean-square error (RMSE), target registration error (TRE), success rate, and processing time. The RMSE, TRE, and the success rate of the HPS were better than those of the other optimization algorithms, and the processing time was similar to that of Powell's algorithm alone. PMID:23138929
NASA Astrophysics Data System (ADS)
Cheng, Cheng; Bodin, Thomas; Allen, Richard M.
2016-02-01
We present a novel 3D pre-stack Kirchhoff depth migration (PKDM) method for teleseismic receiver functions. The proposed algorithm considers the effects of diffraction, scattering, and travel time alteration caused by 3D volumetric heterogeneities. It is therefore particularly useful for imaging complex 3D structures such as dipping discontinuities, which is hard to accomplish with traditional methods. The scheme is based on the acoustic wave migration principle, where at each time step of the receiver function, the energy is migrated back to the ensemble of potential conversion points in the image, given a smooth 3D reference model. Travel times for P and S waves are computed with an efficient Eikonal solver, the Fast Marching Method. We also consider elastic scattering patterns, where the amplitude of converted S waves depends on the angle between the incident P wave, and the scattered S wave. Synthetic experiments demonstrate the validity of the method for a variety of dipping angle discontinuities. Comparison with the widely used Common Conversion Point (CCP) stacking method reveals that our migration shows considerable improvement. For example, the effect of multiple reflections that usually produce apparent discontinuities is avoided. The proposed approach is practical, computationally efficient, and is therefore a potentially powerful alternative to standard CCP methods for imaging large-scale continental structure under dense networks.
Fast and efficient particle reconstruction on a 3D grid using sparsity
NASA Astrophysics Data System (ADS)
Cornic, P.; Champagnat, F.; Cheminet, A.; Leclaire, B.; Le Besnerais, G.
2015-03-01
We propose an approach for efficient localization and intensity reconstruction of particles on a 3D grid based on sparsity principles. The computational complexity of the method is limited by using the particle volume reconstruction paradigm (Champagnat et al. in Meas Sci Technol 25, 2014) and a reduction in the problem dimension. Tests on synthetic and experimental data show that the proposed method leads to more efficient detections and to reconstructions of higher quality than classical tomoPIV approaches on a large range of seeding densities, up to ppp ≈ 0.12.
Atlas Toolkit: Fast registration of 3D morphological datasets in the absence of landmarks
Grocott, Timothy; Thomas, Paul; Münsterberg, Andrea E.
2016-01-01
Image registration is a gateway technology for Developmental Systems Biology, enabling computational analysis of related datasets within a shared coordinate system. Many registration tools rely on landmarks to ensure that datasets are correctly aligned; yet suitable landmarks are not present in many datasets. Atlas Toolkit is a Fiji/ImageJ plugin collection offering elastic group-wise registration of 3D morphological datasets, guided by segmentation of the interesting morphology. We demonstrate the method by combinatorial mapping of cell signalling events in the developing eyes of chick embryos, and use the integrated datasets to predictively enumerate Gene Regulatory Network states. PMID:26864723
3D Multislab, Multishot Acquisition for Fast, Whole-Brain MR Elastography with High SNR Efficiency
Johnson, Curtis L; Holtrop, Joseph L; McGarry, Matthew DJ; Weaver, John B; Paulsen, Keith D; Georgiadis, John G; Sutton, Bradley P
2014-01-01
Purpose To develop an acquisition scheme for generating magnetic resonance elastography (MRE) displacement data with whole-brain coverage, high spatial resolution, and adequate signal-to-noise ratio (SNR) in a short scan time. Theory and Methods A 3D multislab, multishot acquisition for whole-brain MRE with 2.0 mm isotropic spatial resolution is proposed. The multislab approach allowed for the use of short repetition time to achieve very high SNR efficiency. High SNR efficiency allowed for a reduced acquisition time of only six minutes while the minimum SNR needed for inversion was maintained. Results The mechanical property maps estimated from whole-brain displacement data with nonlinear inversion (NLI) demonstrated excellent agreement with neuroanatomical features, including the cerebellum and brainstem. A comparison with an equivalent 2D acquisition illustrated the improvement in SNR efficiency of the 3D multislab acquisition. The flexibility afforded by the high SNR efficiency allowed for higher resolution with a 1.6 mm isotropic voxel size, which generated higher estimates of brainstem stiffness compared with the 2.0 mm isotropic acquisition. Conclusions The acquisition presented allows for the capture of whole-brain MRE displacement data in a short scan time, and may be used to generate local mechanical property estimates of neuroanatomical features throughout the brain. PMID:24347237
Fast and Memory-Efficient Topological Denoising of 2D and 3D Scalar Fields.
Günther, David; Jacobson, Alec; Reininghaus, Jan; Seidel, Hans-Peter; Sorkine-Hornung, Olga; Weinkauf, Tino
2014-12-01
Data acquisition, numerical inaccuracies, and sampling often introduce noise in measurements and simulations. Removing this noise is often necessary for efficient analysis and visualization of this data, yet many denoising techniques change the minima and maxima of a scalar field. For example, the extrema can appear or disappear, spatially move, and change their value. This can lead to wrong interpretations of the data, e.g., when the maximum temperature over an area is falsely reported being a few degrees cooler because the denoising method is unaware of these features. Recently, a topological denoising technique based on a global energy optimization was proposed, which allows the topology-controlled denoising of 2D scalar fields. While this method preserves the minima and maxima, it is constrained by the size of the data. We extend this work to large 2D data and medium-sized 3D data by introducing a novel domain decomposition approach. It allows processing small patches of the domain independently while still avoiding the introduction of new critical points. Furthermore, we propose an iterative refinement of the solution, which decreases the optimization energy compared to the previous approach and therefore gives smoother results that are closer to the input. We illustrate our technique on synthetic and real-world 2D and 3D data sets that highlight potential applications. PMID:26356972
Contrast-Based 3D/2D Registration of the Left Atrium: Fast versus Consistent
Kowalewski, Christopher; Kurzidim, Klaus; Strobel, Norbert; Hornegger, Joachim
2016-01-01
For augmented fluoroscopy during cardiac ablation, a preoperatively acquired 3D model of a patient's left atrium (LA) can be registered to X-ray images recorded during a contrast agent (CA) injection. An automatic registration method that works also for small amounts of CA is desired. We propose two similarity measures: The first focuses on edges of the patient anatomy. The second computes a contrast agent distribution estimate (CADE) inside the 3D model and rates its consistency with the CA as seen in biplane fluoroscopic images. Moreover, temporal filtering on the obtained registration results of a sequence is applied using a Markov chain framework. Evaluation was performed on 11 well-contrasted clinical angiographic sequences and 10 additional sequences with less CA. For well-contrasted sequences, the error for all 73 frames was 7.9 ± 6.3 mm and it dropped to 4.6 ± 4.0 mm when registering to an automatically selected, well enhanced frame in each sequence. Temporal filtering reduced the error for all frames from 7.9 ± 6.3 mm to 5.7 ± 4.6 mm. The error was typically higher if less CA was used. A combination of both similarity measures outperforms a previously proposed similarity measure. The mean accuracy for well contrasted sequences is in the range of other proposed manual registration methods. PMID:27051412
Radial subsampling for fast cost function computation in intensity-based 3D image registration
NASA Astrophysics Data System (ADS)
Boettger, Thomas; Wolf, Ivo; Meinzer, Hans-Peter; Celi, Juan Carlos
2007-03-01
Image registration is always a trade-off between accuracy and speed. Looking towards clinical scenarios the time for bringing two or more images into registration should be around a few seconds only. We present a new scheme for subsampling 3D-image data to allow for efficient computation of cost functions in intensity-based image registration. Starting from an arbitrary center point voxels are sampled along scan lines which do radially extend from the center point. We analyzed the characteristics of different cost functions computed on the sub-sampled data and compared them to known cost functions with respect to local optima. Results show the cost functions are smooth and give high peaks at the expected optima. Furthermore we investigated capture range of cost functions computed under the new subsampling scheme. Capture range was remarkably better for the new scheme compared to metrics using all voxels or different subsampling schemes and high registration accuracy was achieved as well. The most important result is the improvement in terms of speed making this scheme very interesting for clinical scenarios. We conclude using the new subsampling scheme intensity-based 3D image registration can be performed much faster than using other approaches while maintaining high accuracy. A variety of different extensions of the new approach is conceivable, e.g. non-regular distribution of the scan lines or not to let the scan lines start from a center point only, but from the surface of an organ model for example.
Contrast-Based 3D/2D Registration of the Left Atrium: Fast versus Consistent.
Hoffmann, Matthias; Kowalewski, Christopher; Maier, Andreas; Kurzidim, Klaus; Strobel, Norbert; Hornegger, Joachim
2016-01-01
For augmented fluoroscopy during cardiac ablation, a preoperatively acquired 3D model of a patient's left atrium (LA) can be registered to X-ray images recorded during a contrast agent (CA) injection. An automatic registration method that works also for small amounts of CA is desired. We propose two similarity measures: The first focuses on edges of the patient anatomy. The second computes a contrast agent distribution estimate (CADE) inside the 3D model and rates its consistency with the CA as seen in biplane fluoroscopic images. Moreover, temporal filtering on the obtained registration results of a sequence is applied using a Markov chain framework. Evaluation was performed on 11 well-contrasted clinical angiographic sequences and 10 additional sequences with less CA. For well-contrasted sequences, the error for all 73 frames was 7.9 ± 6.3 mm and it dropped to 4.6 ± 4.0 mm when registering to an automatically selected, well enhanced frame in each sequence. Temporal filtering reduced the error for all frames from 7.9 ± 6.3 mm to 5.7 ± 4.6 mm. The error was typically higher if less CA was used. A combination of both similarity measures outperforms a previously proposed similarity measure. The mean accuracy for well contrasted sequences is in the range of other proposed manual registration methods. PMID:27051412
Ruh, Dominic; Tränkle, Benjamin; Rohrbach, Alexander
2011-10-24
Multi-dimensional, correlated particle tracking is a key technology to reveal dynamic processes in living and synthetic soft matter systems. In this paper we present a new method for tracking micron-sized beads in parallel and in all three dimensions - faster and more precise than existing techniques. Using an acousto-optic deflector and two quadrant-photo-diodes, we can track numerous optically trapped beads at up to tens of kHz with a precision of a few nanometers by back-focal plane interferometry. By time-multiplexing the laser focus, we can calibrate individually all traps and all tracking signals in a few seconds and in 3D. We show 3D histograms and calibration constants for nine beads in a quadratic arrangement, although trapping and tracking is easily possible for more beads also in arbitrary 2D arrangements. As an application, we investigate the hydrodynamic coupling and diffusion anomalies of spheres trapped in a 3 × 3 arrangement. PMID:22109012
NASA Astrophysics Data System (ADS)
Huang, Qun-xing; Yan, Fei Wang Jian-hua; Chi, Yong
2013-04-01
A new approach to inverse radiation analysis is presented for non-intrusive 3-D flame temperature reconstruction using flame emission images from four CCD camera detectors installed on the furnace wall. The scattering from participating medium in the flame was considered by combining the discrete radiative transfer method with the discrete ordinate method. A modified minimum residual algorithm was employed to calculate the least squares solution of the ill-conditioned inverse problem. A numerical test problem simulating real temperature measurements in an industrial furnace was used to assess the performance of the proposed method. These assessments indicate that this method is capable of reconstructing 3-D temperature distributions fast and accurately, even with noisy flame emission data. Such a capability has potential in real-time temperature measurement for combustion optimization and pollution emission control.
Jia, Jia; Liu, Juan; Jin, Guofan; Wang, Yongtian
2014-09-20
Occlusion culling is an important process that produces correct depth cues for observers in holographic displays, whereas current methods suffer from occlusion errors or high computational loads. We propose a fast and effective method for occlusion culling based on multiple light-point sampling planes and an inverse orthographic projection technique. Multiple light-point sampling planes are employed to remove the hidden surfaces for each direction of the view of the three-dimensional (3D) scene by forward orthographic projection, and the inverse orthographic projection technique is used to determine the effective sampling points of the 3D scene. A numerical simulation and an optical experiment are performed. The results show that this approach can realize accurate occlusion effects, smooth motion parallax, and continuous depth using low angular sampling without any extra computation costs. PMID:25322109
Very fast road database verification using textured 3D city models obtained from airborne imagery
NASA Astrophysics Data System (ADS)
Bulatov, Dimitri; Ziems, Marcel; Rottensteiner, Franz; Pohl, Melanie
2014-10-01
Road databases are known to be an important part of any geodata infrastructure, e.g. as the basis for urban planning or emergency services. Updating road databases for crisis events must be performed quickly and with the highest possible degree of automation. We present a semi-automatic algorithm for road verification using textured 3D city models, starting from aerial or even UAV-images. This algorithm contains two processes, which exchange input and output, but basically run independently from each other. These processes are textured urban terrain reconstruction and road verification. The first process contains a dense photogrammetric reconstruction of 3D geometry of the scene using depth maps. The second process is our core procedure, since it contains various methods for road verification. Each method represents a unique road model and a specific strategy, and thus is able to deal with a specific type of roads. Each method is designed to provide two probability distributions, where the first describes the state of a road object (correct, incorrect), and the second describes the state of its underlying road model (applicable, not applicable). Based on the Dempster-Shafer Theory, both distributions are mapped to a single distribution that refers to three states: correct, incorrect, and unknown. With respect to the interaction of both processes, the normalized elevation map and the digital orthophoto generated during 3D reconstruction are the necessary input - together with initial road database entries - for the road verification process. If the entries of the database are too obsolete or not available at all, sensor data evaluation enables classification of the road pixels of the elevation map followed by road map extraction by means of vectorization and filtering of the geometrically and topologically inconsistent objects. Depending on the time issue and availability of a geo-database for buildings, the urban terrain reconstruction procedure has semantic models
Mapping of the spontaneous deletion in the Ap3d1 gene of mocha mice: fast and reliable genotyping
Drasbek, Kim Ryun; Holm, Mai Marie; Delenclos, Marion; Jensen, Kimmo
2008-01-01
Background The mocha mouse carries a spontaneous deletion in the Ap3d1 gene, encoding the delta 1 subunit of the adaptor related protein complex 3, (Ap3d1), and subsequently lack the expression of functional AP-3. This leads to a deficiency in vesicle transport and storage, which affects neurotransmitter vesicle turnover and release in the central nervous system. Since the genomic sequence of the Ap3d1 gene of mocha mouse is not known, precise mapping of the deletion as well as reliable genotyping protocols are lacking. Findings We sequenced the Ap3d1 gene (HGNC GeneID: 8943) around the deletion site in the mocha mouse and revealed a 10639 bp deletion covering exon 2 to 6. Subsequently, new PCR primers were designed yielding a reliable genotyping protocol of both newborn and adult tissue. To examine the genotypes further, hippocampal neurons were cultured from mocha and control mice. Patch-clamp recordings showed that mocha neurons had a higher input resistance, and that autaptic EPSC in mocha cultures depressed faster and stronger as compared with control cultures. Conclusion Our study reports the sequence of the deleted part of the Ap3d1 gene in mocha mice, as well as a reliable PCR-based genotyping protocol. We cultured hippocampal neurons from control and mocha mice, and found a difference in input resistance of the neurons, and in the synaptic short-term plasticity of glutamatergic autapses showing a larger synaptic depression than controls. The described procedures may be useful for the future utilization of the mocha mouse as a model of defective vesicle biogenesis. Importantly, as genotyping by eye color is complicated in newborn mice, the designed protocol is so fast and reliable that newborn mice could rapidly be genotyped and hippocampal neurons dissociated and cultured, which is normally best done at P0-P2. PMID:19032734
Fast and Robust Sixth Order Multigrid Computation for 3D Convection Diffusion Equation.
Wang, Yin; Zhang, Jun
2010-10-15
We present a sixth order explicit compact finite difference scheme to solve the three dimensional (3D) convection diffusion equation. We first use multiscale multigrid method to solve the linear systems arising from a 19-point fourth order discretization scheme to compute the fourth order solutions on both the coarse grid and the fine grid. Then an operator based interpolation scheme combined with an extrapolation technique is used to approximate the sixth order accurate solution on the fine grid. Since the multigrid method using a standard point relaxation smoother may fail to achieve the optimal grid independent convergence rate for solving convection diffusion equation with a high Reynolds number, we implement the plane relaxation smoother in the multigrid solver to achieve better grid independency. Supporting numerical results are presented to demonstrate the efficiency and accuracy of the sixth order compact scheme (SOC), compared with the previously published fourth order compact scheme (FOC). PMID:21151737
Fast DRR generation for 2D to 3D registration on GPUs
Tornai, Gabor Janos; Cserey, Gyoergy
2012-08-15
Purpose: The generation of digitally reconstructed radiographs (DRRs) is the most time consuming step on the CPU in intensity based two-dimensional x-ray to three-dimensional (CT or 3D rotational x-ray) medical image registration, which has application in several image guided interventions. This work presents optimized DRR rendering on graphical processor units (GPUs) and compares performance achievable on four commercially available devices. Methods: A ray-cast based DRR rendering was implemented for a 512 Multiplication-Sign 512 Multiplication-Sign 72 CT volume. The block size parameter was optimized for four different GPUs for a region of interest (ROI) of 400 Multiplication-Sign 225 pixels with different sampling ratios (1.1%-9.1% and 100%). Performance was statistically evaluated and compared for the four GPUs. The method and the block size dependence were validated on the latest GPU for several parameter settings with a public gold standard dataset (512 Multiplication-Sign 512 Multiplication-Sign 825 CT) for registration purposes. Results: Depending on the GPU, the full ROI is rendered in 2.7-5.2 ms. If sampling ratio of 1.1%-9.1% is applied, execution time is in the range of 0.3-7.3 ms. On all GPUs, the mean of the execution time increased linearly with respect to the number of pixels if sampling was used. Conclusions: The presented results outperform other results from the literature. This indicates that automatic 2D to 3D registration, which typically requires a couple of hundred DRR renderings to converge, can be performed quasi on-line, in less than a second or depending on the application and hardware in less than a couple of seconds. Accordingly, a whole new field of applications is opened for image guided interventions, where the registration is continuously performed to match the real-time x-ray.
NASA Technical Reports Server (NTRS)
Hung, Stephen H. Y.
1989-01-01
A fast 3-D object recognition algorithm that can be used as a quick-look subsystem to the vision system for the Special-Purpose Dexterous Manipulator (SPDM) is described. Global features that can be easily computed from range data are used to characterize the images of a viewer-centered model of an object. This algorithm will speed up the processing by eliminating the low level processing whenever possible. It may identify the object, reject a set of bad data in the early stage, or create a better environment for a more powerful algorithm to carry the work further.
A fast technique applied to the analysis of Resistive Wall Modes with 3D conducting structures
Rubinacci, Guglielmo Liu, Yueqiang
2009-03-20
This paper illustrates the development of a 'fast' technique for the analysis of Resistive Wall Modes (RWMs) in fusion devices with three-dimensional conducting structures, by means of the recently developed CarMa code. Thanks to its peculiar features, the computational cost scales almost linearly with the number of discrete unknowns. Some large scale problems are solved in configurations of interest for the International Thermonuclear Experimental Reactor (ITER)
A fast nested dissection solver for Cartesian 3D elliptic problems using hierarchical matrices
NASA Astrophysics Data System (ADS)
Schmitz, Phillip G.; Ying, Lexing
2014-02-01
We present a fast algorithm for solutions to linear systems arising from three dimensional elliptic problems on a regular Cartesian mesh. We follow the approach of Schmitz and Ying (2012) on combining the nested dissection matrix factorization method with hierarchical matrices in two dimensions and extend it to the three dimensional case. A theoretical linear time complexity is derived and a more practical variant with slightly worse scaling is demonstrated.
Poulin, Eric; Racine, Emmanuel; Beaulieu, Luc; Binnekamp, Dirk
2015-03-15
Purpose: In high dose rate brachytherapy (HDR-B), current catheter reconstruction protocols are relatively slow and error prone. The purpose of this technical note is to evaluate the accuracy and the robustness of an electromagnetic (EM) tracking system for automated and real-time catheter reconstruction. Methods: For this preclinical study, a total of ten catheters were inserted in gelatin phantoms with different trajectories. Catheters were reconstructed using a 18G biopsy needle, used as an EM stylet and equipped with a miniaturized sensor, and the second generation Aurora{sup ®} Planar Field Generator from Northern Digital Inc. The Aurora EM system provides position and orientation value with precisions of 0.7 mm and 0.2°, respectively. Phantoms were also scanned using a μCT (GE Healthcare) and Philips Big Bore clinical computed tomography (CT) system with a spatial resolution of 89 μm and 2 mm, respectively. Reconstructions using the EM stylet were compared to μCT and CT. To assess the robustness of the EM reconstruction, five catheters were reconstructed twice and compared. Results: Reconstruction time for one catheter was 10 s, leading to a total reconstruction time inferior to 3 min for a typical 17-catheter implant. When compared to the μCT, the mean EM tip identification error was 0.69 ± 0.29 mm while the CT error was 1.08 ± 0.67 mm. The mean 3D distance error was found to be 0.66 ± 0.33 mm and 1.08 ± 0.72 mm for the EM and CT, respectively. EM 3D catheter trajectories were found to be more accurate. A maximum difference of less than 0.6 mm was found between successive EM reconstructions. Conclusions: The EM reconstruction was found to be more accurate and precise than the conventional methods used for catheter reconstruction in HDR-B. This approach can be applied to any type of catheters and applicators.
Observations and models of magnetic reconnection
NASA Astrophysics Data System (ADS)
Barta, Miroslav
2015-08-01
Magnetic reconnection is now almost unanimously considered to be a key plasma process for energy release in solar and stellar flares. Recent decade have seen rapid development in the theory, simulations and searching for observational evidences of magnetic reconnection being in action in the core of flares. Modern modeling approach involves many realistic aspects of magnetic reconnection such as intrinsically 3D nature of the process and, namely, its highly dynamic character connected with violent formation of plasmoids at many scales. The cascade of plasmoid formation represents natural process of fast, turbulent energy transfer to the kinetic dissipation scale. This concept, revealed by numerical simulations, has found its ground in the theory of (ideal) plasmoid instability in current layers with high aspect ratio. The plasmoid dominated reconnection regime is capable to account for many puzzling dilemmas in the flare physics ranging from the observation-demanded energy release rate vs. standard reconnection-regime timescales, observed organized large-scale structures vs. signatures of fragmented energy release etc. The talk aims at reviewing recent theoretical and simulation development in this direction and observational support for the concept of plasmoid-driven reconnection cascade namely in solar flares.
Development of a 3D numerical methodology for fast prediction of gun blast induced loading
NASA Astrophysics Data System (ADS)
Costa, E.; Lagasco, F.
2014-05-01
In this paper, the development of a methodology based on semi-empirical models from the literature to carry out 3D prediction of pressure loading on surfaces adjacent to a weapon system during firing is presented. This loading is consequent to the impact of the blast wave generated by the projectile exiting the muzzle bore. When exceeding a pressure threshold level, loading is potentially capable to induce unwanted damage to nearby hard structures as well as frangible panels or electronic equipment. The implemented model shows the ability to quickly predict the distribution of the blast wave parameters over three-dimensional complex geometry surfaces when the weapon design and emplacement data as well as propellant and projectile characteristics are available. Considering these capabilities, the use of the proposed methodology is envisaged as desirable in the preliminary design phase of the combat system to predict adverse effects and then enable to identify the most appropriate countermeasures. By providing a preliminary but sensitive estimate of the operative environmental loading, this numerical means represents a good alternative to more powerful, but time consuming advanced computational fluid dynamics tools, which use can, thus, be limited to the final phase of the design.
Fast simulation of solid tumors thermal ablation treatments with a 3D reaction diffusion model.
Bertaccini, Daniele; Calvetti, Daniela
2007-08-01
An efficient computational method for near real-time simulation of thermal ablation of tumors via radio frequencies is proposed. Model simulations of the temperature field in a 3D portion of tissue containing the tumoral mass for different patterns of source heating can be used to design the ablation procedure. The availability of a very efficient computational scheme makes it possible to update the predicted outcome of the procedure in real time. In the algorithms proposed here a discretization in space of the governing equations is followed by an adaptive time integration based on implicit multistep formulas. A modification of the ode15s MATLAB function which uses Krylov space iterative methods for the solution of the linear systems arising at each integration step makes it possible to perform the simulations on standard desktop for much finer grids than using the built-in ode15s. The proposed algorithm can be applied to a wide class of nonlinear parabolic differential equations. PMID:17173888
FAST SIMULATION OF SOLID TUMORS THERMAL ABLATION TREATMENTS WITH A 3D REACTION DIFFUSION MODEL *
BERTACCINI, DANIELE; CALVETTI, DANIELA
2007-01-01
An efficient computational method for near real-time simulation of thermal ablation of tumors via radio frequencies is proposed. Model simulations of the temperature field in a 3D portion of tissue containing the tumoral mass for different patterns of source heating can be used to design the ablation procedure. The availability of a very efficient computational scheme makes it possible update the predicted outcome of the procedure in real time. In the algorithms proposed here a discretization in space of the governing equations is followed by an adaptive time integration based on implicit multistep formulas. A modification of the ode15s MATLAB function which uses Krylov space iterative methods for the solution of for the linear systems arising at each integration step makes it possible to perform the simulations on standard desktop for much finer grids than using the built-in ode15s. The proposed algorithm can be applied to a wide class of nonlinear parabolic differential equations. PMID:17173888
Sakhalkar, H S; Oldham, M
2008-01-01
This study introduces a charge coupled device (CCD) area detector based optical-computed tomography (optical-CT) scanner for comprehensive verification of radiation dose distributions recorded in nonscattering radiochromic dosimeters. Defining characteristics include: (i) a very fast scanning time of approximately 5 min to acquire a complete three-dimensional (3D) dataset, (ii) improved image formation through the use of custom telecentric optics, which ensures accurate projection images and minimizes artifacts from scattered and stray-light sources, and (iii) high resolution (potentially 50 microm) isotropic 3D dose readout. The performance of the CCD scanner for 3D dose readout was evaluated by comparison with independent 3D readout from the single laser beam OCTOPUS-scanner for the same PRESAGE dosimeters. The OCTOPUS scanner was considered the "gold standard" technique in light of prior studies demonstrating its accuracy. Additional comparisons were made against calculated dose distributions from the ECLIPSE treatment-planning system. Dose readout for the following treatments were investigated: (i) a single rectangular beam irradiation to investigate small field and very steep dose gradient dosimetry away from edge effects, (ii) a 2-field open beam parallel-opposed irradiation to investigate dosimetry along steep dose gradients, and (iii) a 7-field intensity modulated radiation therapy (IMRT) irradiation to investigate dosimetry for complex treatment delivery involving modulation of fluence and for dosimetry along moderate dose gradients. Dose profiles, dose-difference plots, and gamma maps were employed to evaluate quantitative estimates of agreement between independently measured and calculated dose distributions. Results indicated that dose readout from the CCD scanner was in agreement with independent gold-standard readout from the OCTOPUS-scanner as well as the calculated ECLIPSE dose distribution for all treatments, except in regions within a few
NCC-RANSAC: A Fast Plane Extraction Method for 3-D Range Data Segmentation
Qian, Xiangfei; Ye, Cang
2015-01-01
This paper presents a new plane extraction (PE) method based on the random sample consensus (RANSAC) approach. The generic RANSAC-based PE algorithm may over-extract a plane, and it may fail in case of a multistep scene where the RANSAC procedure results in multiple inlier patches that form a slant plane straddling the steps. The CC-RANSAC PE algorithm successfully overcomes the latter limitation if the inlier patches are separate. However, it fails if the inlier patches are connected. A typical scenario is a stairway with a stair wall where the RANSAC plane-fitting procedure results in inliers patches in the tread, riser, and stair wall planes. They connect together and form a plane. The proposed method, called normal-coherence CC-RANSAC (NCC-RANSAC), performs a normal coherence check to all data points of the inlier patches and removes the data points whose normal directions are contradictory to that of the fitted plane. This process results in separate inlier patches, each of which is treated as a candidate plane. A recursive plane clustering process is then executed to grow each of the candidate planes until all planes are extracted in their entireties. The RANSAC plane-fitting and the recursive plane clustering processes are repeated until no more planes are found. A probabilistic model is introduced to predict the success probability of the NCC-RANSAC algorithm and validated with real data of a 3-D time-of-flight camera–SwissRanger SR4000. Experimental results demonstrate that the proposed method extracts more accurate planes with less computational time than the existing RANSAC-based methods. PMID:24771605
NASA Astrophysics Data System (ADS)
Zhang, Lisha
We present fast and robust numerical algorithms for 3-D scattering from perfectly electrical conducting (PEC) and dielectric random rough surfaces in microwave remote sensing. The Coifman wavelets or Coiflets are employed to implement Galerkin's procedure in the method of moments (MoM). Due to the high-precision one-point quadrature, the Coiflets yield fast evaluations of the most off-diagonal entries, reducing the matrix fill effort from O(N2) to O( N). The orthogonality and Riesz basis of the Coiflets generate well conditioned impedance matrix, with rapid convergence for the conjugate gradient solver. The resulting impedance matrix is further sparsified by the matrix-formed standard fast wavelet transform (SFWT). By properly selecting multiresolution levels of the total transformation matrix, the solution precision can be enhanced while matrix sparsity and memory consumption have not been noticeably sacrificed. The unified fast scattering algorithm for dielectric random rough surfaces can asymptotically reduce to the PEC case when the loss tangent grows extremely large. Numerical results demonstrate that the reduced PEC model does not suffer from ill-posed problems. Compared with previous publications and laboratory measurements, good agreement is observed.
Fast particles identification in programmable form at level-0 trigger by means of the 3D-Flow system
Crosetto, Dario B.
1998-10-30
The 3D-Flow Processor system is a new, technology-independent concept in very fast, real-time system architectures. Based on either an FPGA or an ASIC implementation, it can address, in a fully programmable manner, applications where commercially available processors would fail because of throughput requirements. Possible applications include filtering-algorithms (pattern recognition) from the input of multiple sensors, as well as moving any input validated by these filtering-algorithms to a single output channel. Both operations can easily be implemented on a 3D-Flow system to achieve a real-time processing system with a very short lag time. This system can be built either with off-the-shelf FPGAs or, for higher data rates, with CMOS chips containing 4 to 16 processors each. The basic building block of the system, a 3D-Flow processor, has been successfully designed in VHDL code written in ''Generic HDL'' (mostly made of reusable blocks that are synthesizable in different technologies, or FPGAs), to produce a netlist for a four-processor ASIC featuring 0.35 micron CBA (Ceil Base Array) technology at 3.3 Volts, 884 mW power dissipation at 60 MHz and 63.75 mm sq. die size. The same VHDL code has been targeted to three FPGA manufacturers (Altera EPF10K250A, ORCA-Lucent Technologies 0R3T165 and Xilinx XCV1000). A complete set of software tools, the 3D-Flow System Manager, equally applicable to ASIC or FPGA implementations, has been produced to provide full system simulation, application development, real-time monitoring, and run-time fault recovery. Today's technology can accommodate 16 processors per chip in a medium size die, at a cost per processor of less than $5 based on the current silicon die/size technology cost.
Suppression law of quantum states in a 3D photonic fast Fourier transform chip.
Crespi, Andrea; Osellame, Roberto; Ramponi, Roberta; Bentivegna, Marco; Flamini, Fulvio; Spagnolo, Nicolò; Viggianiello, Niko; Innocenti, Luca; Mataloni, Paolo; Sciarrino, Fabio
2016-01-01
The identification of phenomena able to pinpoint quantum interference is attracting large interest. Indeed, a generalization of the Hong-Ou-Mandel effect valid for any number of photons and optical modes would represent an important leap ahead both from a fundamental perspective and for practical applications, such as certification of photonic quantum devices, whose computational speedup is expected to depend critically on multi-particle interference. Quantum distinctive features have been predicted for many particles injected into multimode interferometers implementing the Fourier transform over the optical modes. Here we develop a scalable approach for the implementation of the fast Fourier transform algorithm using three-dimensional photonic integrated interferometers, fabricated via femtosecond laser writing technique. We observe the suppression law for a large number of output states with four- and eight-mode optical circuits: the experimental results demonstrate genuine quantum interference between the injected photons, thus offering a powerful tool for diagnostic of photonic platforms. PMID:26843135
Fast spatial beam shaping by acousto-optic diffraction for 3D non-linear microscopy.
Akemann, Walther; Léger, Jean-François; Ventalon, Cathie; Mathieu, Benjamin; Dieudonné, Stéphane; Bourdieu, Laurent
2015-11-01
Acousto-optic deflection (AOD) devices offer unprecedented fast control of the entire spatial structure of light beams, most notably their phase. AOD light modulation of ultra-short laser pulses, however, is not straightforward to implement because of intrinsic chromatic dispersion and non-stationarity of acousto-optic diffraction. While schemes exist to compensate chromatic dispersion, non-stationarity remains an obstacle. In this work we demonstrate an efficient AOD light modulator for stable phase modulation using time-locked generation of frequency-modulated acoustic waves at the full repetition rate of a high power laser pulse amplifier of 80 kHz. We establish the non-local relationship between the optical phase and the generating acoustic frequency function and verify the system for temporal stability, phase accuracy and generation of non-linear two-dimensional phase functions. PMID:26561090
Suppression law of quantum states in a 3D photonic fast Fourier transform chip
Crespi, Andrea; Osellame, Roberto; Ramponi, Roberta; Bentivegna, Marco; Flamini, Fulvio; Spagnolo, Nicolò; Viggianiello, Niko; Innocenti, Luca; Mataloni, Paolo; Sciarrino, Fabio
2016-01-01
The identification of phenomena able to pinpoint quantum interference is attracting large interest. Indeed, a generalization of the Hong–Ou–Mandel effect valid for any number of photons and optical modes would represent an important leap ahead both from a fundamental perspective and for practical applications, such as certification of photonic quantum devices, whose computational speedup is expected to depend critically on multi-particle interference. Quantum distinctive features have been predicted for many particles injected into multimode interferometers implementing the Fourier transform over the optical modes. Here we develop a scalable approach for the implementation of the fast Fourier transform algorithm using three-dimensional photonic integrated interferometers, fabricated via femtosecond laser writing technique. We observe the suppression law for a large number of output states with four- and eight-mode optical circuits: the experimental results demonstrate genuine quantum interference between the injected photons, thus offering a powerful tool for diagnostic of photonic platforms. PMID:26843135
LBM-EP: Lattice-Boltzmann method for fast cardiac electrophysiology simulation from 3D images.
Rapaka, S; Mansi, T; Georgescu, B; Pop, M; Wright, G A; Kamen, A; Comaniciu, Dorin
2012-01-01
Current treatments of heart rhythm troubles require careful planning and guidance for optimal outcomes. Computational models of cardiac electrophysiology are being proposed for therapy planning but current approaches are either too simplified or too computationally intensive for patient-specific simulations in clinical practice. This paper presents a novel approach, LBM-EP, to solve any type of mono-domain cardiac electrophysiology models at near real-time that is especially tailored for patient-specific simulations. The domain is discretized on a Cartesian grid with a level-set representation of patient's heart geometry, previously estimated from images automatically. The cell model is calculated node-wise, while the transmembrane potential is diffused using Lattice-Boltzmann method within the domain defined by the level-set. Experiments on synthetic cases, on a data set from CESC'10 and on one patient with myocardium scar showed that LBM-EP provides results comparable to an FEM implementation, while being 10 - 45 times faster. Fast, accurate, scalable and requiring no specific meshing, LBM-EP paves the way to efficient and detailed models of cardiac electrophysiology for therapy planning. PMID:23286029
Field calibration of binocular stereo vision based on fast reconstruction of 3D control field
NASA Astrophysics Data System (ADS)
Zhang, Haijun; Liu, Changjie; Fu, Luhua; Guo, Yin
2015-08-01
Construction of high-speed railway in China has entered a period of rapid growth. To accurately and quickly obtain the dynamic envelope curve of high-speed vehicle is an important guarantee for safe driving. The measuring system is based on binocular stereo vision. Considering the difficulties in field calibration such as environmental changes and time limits, carried out a field calibration method based on fast reconstruction of three-dimensional control field. With the rapid assembly of pre-calibrated three-dimensional control field, whose coordinate accuracy is guaranteed by manufacture accuracy and calibrated by V-STARS, two cameras take a quick shot of it at the same time. The field calibration parameters are then solved by the method combining linear solution with nonlinear optimization. Experimental results showed that the measurement accuracy can reach up to +/- 0.5mm, and more importantly, in the premise of guaranteeing accuracy, the speed of the calibration and the portability of the devices have been improved considerably.
Comparison of fast 3D simulation and actinic inspection for EUV masks with buries defects
Clifford, C. H.; Wiraatmadja, S.; Chan, T. T.; Neureuther, A. R.; Goldberg, K. A.; Mochi, I.; Liang, T.
2009-02-23
Aerial images for isolated defects and the interactions of defects with features are compared between the Actinic Inspection Tool (AIT) at Lawrence Berkeley National Laboratory (LBNL) and the fast EUV simulation program RADICAL. Comparisons between AIT images from August 2007 and RADICAL simulations are used to extract aberrations. At this time astigmatism was the dominant aberration with a value of 0.55 waves RMS. Significant improvements in the imaging performance of the AIT were made between August 2007 and December 2008. A good match will be shown between the most recent AIT images and RADICAL simulations without aberrations. These comparisons will demonstrate that a large defect, in this case 7nm tall on the surface, is still printable even if it is centered under the absorber line. These comparisons also suggest that the minimum defect size is between 1.5nm and 0.8nm surface height because a 1.5nm defect was printable but a 0.8nm was not. Finally, the image of a buried defect near an absorber line through focus will demonstrate an inversion in the effect of the defect from a protrusion of the dark line into the space to a protrusion of the space into the line.
Sanz-Requena, Roberto; Moratal, David; García-Sánchez, Diego Ramón; Bodí, Vicente; Rieta, José Joaquín; Sanchis, Juan Manuel
2007-03-01
Intravascular ultrasound (IVUS) imaging is used along with X-ray coronary angiography to detect vessel pathologies. Manual analysis of IVUS images is slow and time-consuming and it is not feasible for clinical purposes. A semi-automated method is proposed to generate 3D reconstructions from IVUS video sequences, so that a fast diagnose can be easily done, quantifying plaque length and severity as well as plaque volume of the vessels under study. The methodology described in this work has four steps: a pre-processing of IVUS images, a segmentation of media-adventitia contour, a detection of intima and plaque and a 3D reconstruction of the vessel. Preprocessing is intended to remove noise from the images without blurring the edges. Segmentation of media-adventitia contour is achieved using active contours (snakes). In particular, we use the gradient vector flow (GVF) as external force for the snakes. The detection of lumen border is obtained taking into account gray-level information of the inner part of the previously detected contours. A knowledge-based approach is used to determine which level of gray corresponds statistically to the different regions of interest: intima, plaque and lumen. The catheter region is automatically discarded. An estimate of plaque type is also given. Finally, 3D reconstruction of all detected regions is made. The suitability of this methodology has been verified for the analysis and visualization of plaque length, stenosis severity, automatic detection of the most problematic regions, calculus of plaque volumes and a preliminary estimation of plaque type obtaining for automatic measures of lumen and vessel area an average error smaller than 1mm(2) (equivalent aproximately to 10% of the average measure), for calculus of plaque and lumen volume errors smaller than 0.5mm(3) (equivalent approximately to 20% of the average measure) and for plaque type estimates a mismatch of less than 8% in the analysed frames. PMID:17215103
Fast, large field-of-view, telecentric optical-CT scanning system for 3D radiochromic dosimetry
Thomas, A; Oldham, M
2010-01-01
We describe initial experiences with an in-house, fast, large field-of-view optical-CT telecentric scanner (the Duke Large field of view Optical-CT Scanner (DLOS)). The DLOS system is designed to enable telecentric optical-CT imaging of dosimeters up to 24 cm in diameter with a spatial resolution of 1 mm3, in approximately 10 minutes. These capabilities render the DLOS system a unique device at present. The system is a scaled up version of early prototypes in our lab. This scaling introduces several challenges, including the accurate measurement of a greatly increased range of light attenuation within the dosimeter, and the need to reduce even minor reflections and scattered light within the imaging chain. We present several corrections and techniques that enable accurate, low noise, 3D dosimetery with the DLOS system. PMID:21218169
Fast, Accurate and Precise Mid-Sagittal Plane Location in 3D MR Images of the Brain
NASA Astrophysics Data System (ADS)
Bergo, Felipe P. G.; Falcão, Alexandre X.; Yasuda, Clarissa L.; Ruppert, Guilherme C. S.
Extraction of the mid-sagittal plane (MSP) is a key step for brain image registration and asymmetry analysis. We present a fast MSP extraction method for 3D MR images, based on automatic segmentation of the brain and on heuristic maximization of the cerebro-spinal fluid within the MSP. The method is robust to severe anatomical asymmetries between the hemispheres, caused by surgical procedures and lesions. The method is also accurate with respect to MSP delineations done by a specialist. The method was evaluated on 64 MR images (36 pathological, 20 healthy, 8 synthetic), and it found a precise and accurate approximation of the MSP in all of them with a mean time of 60.0 seconds per image, mean angular variation within a same image (precision) of 1.26o and mean angular difference from specialist delineations (accuracy) of 1.64o.
NASA Astrophysics Data System (ADS)
Lee, Kyungmoo; Zhang, Li; Abramoff, Michael D.; Sonka, Milan
2015-03-01
Image segmentation is important for quantitative analysis of medical image data. Recently, our research group has introduced a 3-D graph search method which can simultaneously segment optimal interacting surfaces with respect to the cost function in volumetric images. Although it provides excellent segmentation accuracy, it is computationally demanding (both CPU and memory) to simultaneously segment multiple surfaces from large volumetric images. Therefore, we propose a new, fast, and memory-efficient graph search method for intraretinal layer segmentation of 3-D macular optical coherence tomograpy (OCT) scans. The key idea is to reduce the size of a graph by combining the nodes with high costs based on the multiscale approach. The new approach requires significantly less memory and achieves significantly faster processing speeds (p < 0.01) with only small segmentation differences compared to the original graph search method. This paper discusses sub-optimality of this approach and assesses trade-off relationships between decreasing processing speed and increasing segmentation differences from that of the original method as a function of employed scale of the underlying graph construction.
HIFI-C: a robust and fast method for determining NMR couplings from adaptive 3D to 2D projections.
Cornilescu, Gabriel; Bahrami, Arash; Tonelli, Marco; Markley, John L; Eghbalnia, Hamid R
2007-08-01
We describe a novel method for the robust, rapid, and reliable determination of J couplings in multi-dimensional NMR coupling data, including small couplings from larger proteins. The method, "High-resolution Iterative Frequency Identification of Couplings" (HIFI-C) is an extension of the adaptive and intelligent data collection approach introduced earlier in HIFI-NMR. HIFI-C collects one or more optimally tilted two-dimensional (2D) planes of a 3D experiment, identifies peaks, and determines couplings with high resolution and precision. The HIFI-C approach, demonstrated here for the 3D quantitative J method, offers vital features that advance the goal of rapid and robust collection of NMR coupling data. (1) Tilted plane residual dipolar couplings (RDC) data are collected adaptively in order to offer an intelligent trade off between data collection time and accuracy. (2) Data from independent planes can provide a statistical measure of reliability for each measured coupling. (3) Fast data collection enables measurements in cases where sample stability is a limiting factor (for example in the presence of an orienting medium required for residual dipolar coupling measurements). (4) For samples that are stable, or in experiments involving relatively stronger couplings, robust data collection enables more reliable determinations of couplings in shorter time, particularly for larger biomolecules. As a proof of principle, we have applied the HIFI-C approach to the 3D quantitative J experiment to determine N-C' RDC values for three proteins ranging from 56 to 159 residues (including a homodimer with 111 residues in each subunit). A number of factors influence the robustness and speed of data collection. These factors include the size of the protein, the experimental set up, and the coupling being measured, among others. To exhibit a lower bound on robustness and the potential for time saving, the measurement of dipolar couplings for the N-C' vector represents a realistic
The ATLAS3D project - XIX. The hot gas content of early-type galaxies: fast versus slow rotators
NASA Astrophysics Data System (ADS)
Sarzi, Marc; Alatalo, Katherine; Blitz, Leo; Bois, Maxime; Bournaud, Frédéric; Bureau, Martin; Cappellari, Michele; Crocker, Alison; Davies, Roger L.; Davis, Timothy A.; de Zeeuw, P. T.; Duc, Pierre-Alain; Emsellem, Eric; Khochfar, Sadegh; Krajnović, Davor; Kuntschner, Harald; Lablanche, Pierre-Yves; McDermid, Richard M.; Morganti, Raffaella; Naab, Thorsten; Oosterloo, Tom; Scott, Nicholas; Serra, Paolo; Young, Lisa M.; Weijmans, Anne-Marie
2013-07-01
For early-type galaxies, the ability to sustain a corona of hot, X-ray-emitting gas could have played a key role in quenching their star formation history. A halo of hot gas may act as an effective shield against the acquisition of cold gas and can quickly absorb stellar mass loss material. Yet, since the discovery by the Einstein Observatory of such X-ray haloes around early-type galaxies, the precise amount of hot gas around these galaxies still remains a matter of debate. By combining homogeneously derived photometric and spectroscopic measurements for the early-type galaxies observed as part of the ATLAS3D integral field survey with measurements of their X-ray luminosity based on X-ray data of both low and high spatial resolution (for 47 and 19 objects, respectively) we conclude that the hot gas content of early-type galaxies can depend on their dynamical structure. Specifically, whereas slow rotators generally have X-ray haloes with luminosity LX, gas and temperature T values that are well in line with what is expected if the hot gas emission is sustained by the thermalization of the kinetic energy carried by the stellar mass loss material, fast rotators tend to display LX, gas values that fall consistently below the prediction of this model, with similar T values that do not scale with the stellar kinetic energy (traced by the stellar velocity dispersion) as observed in the case of slow rotators. Such a discrepancy between the hot gas content of slow and fast rotators would appear to reduce, or even disappear, for large values of the dynamical mass (above ˜3 × 1011 M⊙), with younger fast rotators displaying also somewhat larger LX, gas values possibly owing to the additional energy input from recent supernovae explosions. Considering that fast rotators are likely to be intrinsically flatter than slow rotators, and that the few LX, gas-deficient slow rotators also happen to be relatively flat, the observed LX, gas deficiency in these objects would support
Magnetic reconnection in space
Boozer, Allen H.
2012-09-15
Models of magnetic reconnection in space plasmas generally consider only a segment of the magnetic field lines. The consideration of only a segment of the lines is shown to lead to paradoxical results in which reconnection can be impossible even in a magnetic field constrained to be curl free or can be at an Alfven rate even when the plasma is a perfect conductor. A model of reconnecting magnetic fields is developed which shows the smallness of the interdiffusion distance {delta}{sub d} of magnetic field lines does not limit the speed of reconnection but does provide a reconnection trigger. When the reconnection region has a natural length L{sub r}, the spatial scale of the gradient of magnetic field across the magnetic field lines must reach L{sub g} Almost-Equal-To 0.3L{sub r}/ln(L{sub r}/{delta}{sub d}) for fast reconnection to be triggered, which implies a current density j Almost-Equal-To B/{mu}{sub 0}L{sub g} that is far lower than that usually thought required for fast reconnection. The relation between magnetic reconnection in space and in toroidal laboratory plasmas is also discussed.
NASA Astrophysics Data System (ADS)
Cooper, Wilfred A.; Pfefferle, David; Graves, Jonathan P.
2014-10-01
Designed to accurately solve the motion of energetic particles in the presence of 3D magnetic fields, the VENUS-LEVIS code leans on a non-canonical general coordinate Lagrangian formulation of the equations of motion. It switches between full-orbit particle following and guiding-centre tracing by verifying the perpendicular variation of magnetic vector field, not only including gradients and curvature terms but also the shearing of field-lines. The criteria is particularly relevant for the study of fast ion redistribution in the kinked core of hybrid plasmas, where the compression of flux-surfaces against the axisymmetric outer mantle creates strongly varying magnetic field-lines and large parallel currents. Slowing-down simulations of NBI fast ions show that co-passing particles helically align in the opposite side of the plasma deformation whereas counter-passing particles are barely affected by the kinked structure. Results are compared with experimental neutron camera traces and FIDA measurements during long-lived modes (LLM).
NASA Astrophysics Data System (ADS)
Czarske, J. W.; Kuschmierz, R.; Günther, P.
2013-06-01
Precise measurements of distance, eccentricity and 3D-shape of fast moving objects such as turning parts of lathes, gear shafts, magnetic bearings, camshafts, crankshafts and rotors of vacuum pumps are on the one hand important tasks. On the other hand they are big challenges, since contactless precise measurement techniques are required. Optical techniques are well suitable for distance measurements of non-moving surfaces. However, measurements of laterally fast moving surfaces are still challenging. For such tasks the laser Doppler distance sensor technique was invented by the TU Dresden some years ago. This technique has been realized by two mutually tilted interference fringe systems, where the distance is coded in the phase difference between the generated interference signals. However, due to the speckle effect different random envelopes and phase jumps of the interference signals occur. They disturb the phase difference estimation between the interference signals. In this paper, we will report on a scientific breakthrough on the measurement uncertainty budget which has been achieved recently. Via matching of the illumination and receiving optics the measurement uncertainty of the displacement and distance can be reduced by about one magnitude. For displacement measurements of a recurring rough surface a standard deviation of 110 nm were attained at lateral velocities of 5 m / s. Due to the additionally measured lateral velocity and the rotational speed, the two-dimensional shape of rotating objects is calculated. The three-dimensional shape can be conducted by employment of a line camera. Since the measurement uncertainty of the displacement, vibration, distance, eccentricity, and shape is nearly independent of the lateral surface velocity, this technique is predestined for fast-rotating objects. Especially it can be advantageously used for the quality control of workpieces inside of a lathe towards the reduction of process tolerances, installation times and
Mitsouras, Dimitris; Mulkern, Robert V.; Rybicki, Frank J.
2006-01-01
Fast Spin Echo (FSE) trains elicited by non-selective “hard” refocusing radio frequency (RF) pulses have been proposed as a means to enable application of FSE methods for high resolution 3D magnetic resonance imaging (MRI). Hard-pulse FSE (HPFSE) trains offer short (3–4 ms) echo spacings, but are unfortunately limited to imaging the entire sample within the coil sensitivity thus requiring lengthy imaging times, consequently limiting clinical application. In this work we formulate and analyze two general purpose combinations of 3D HPFSE with Inner Volume (IV) MR imaging to circumvent this limitation. The first method employs a 2D selective RF excitation followed by the HPFSE train, and focuses on required properties of the spatial excitation profile with respect to limiting RF pulse duration in the 5–6 ms range. The second method employs two orthogonally selective 1D RF excitations (a 90x°– 180y° pair) to generate an echo from magnetization within the volume defined by their intersection. Subsequent echoes are formed via the HPFSE train, placing the focus of the method on (a) avoiding spurious echoes that may arise from transverse magnetization located outside the slab intersection when it is unavoidably affected by the non-selective refocusing pulses, and (b) avoiding signal losses due to the necessarily different spacing (in time) of the RF pulse applications. The performance of each method is experimentally measured using Carr-Purcell-Meiboom-Gill (CPMG) multi-echo imaging, enabling examination of the magnetization evolution throughout the echo train. The methods as implemented achieve 95% to 97% outer volume signal suppression, and higher suppression appears to be well within reach, by further refinement of the selective RF excitations. Example images of the human brain and spine are presented with each technique. We conclude that the SNR effciency of volume imaging in conjunction with the short echo spacing afforded by hard pulse trains enable high
Optimized design of a LED-array-based TOF range imaging sensor for fast 3-D shape measurement
NASA Astrophysics Data System (ADS)
Wang, Huanqin; Wang, Ying; Xu, Jun; He, Deyong; Zhao, Tianpeng; Ming, Hai; Kong, Deyi
2011-06-01
A LED-array-based range imaging sensor using Time-of-Flight (TOF) distance measurement was developed to capture the depth information of three-dimensional (3-D) object. By time-division electronic scanning of the LED heterodyne phase-shift TOF range finders in array, the range images were fast obtained without any mechanical moving parts. The design of LED-array-based range imaging sensor was adequately described and a range imaging theoretical model based on photoelectric signal processing was built, which showed there was mutual restriction relationship among the measurement time of a depth pixel, the bandwidth of receiver and the sensor's signal-to-noise ratio (SNR). In order to improve the key parameters of sensor such as range resolution and measurement speed simultaneously, some optimized designs needed to be done for the proposed range imaging sensor, including choosing proper parameters for the filters in receiver, adopting special structure feedback automatic gain control (AGC) circuit with short response time, etc. The final experiment results showed the sensor after optimization could acquire the range images at a rate of 10 frames per second with a range resolution as high as +/-2mm in the range of 50-1200mm. The essential advantages of the proposed range imaging sensor were construction with simple structure, high range resolution, short measurement time and low cost, which was sufficient for many robotic and industrial automation applications.
Medical applications of fast 3D cameras in real-time image-guided radiotherapy (IGRT) of cancer
NASA Astrophysics Data System (ADS)
Li, Shidong; Li, Tuotuo; Geng, Jason
2013-03-01
Dynamic volumetric medical imaging (4DMI) has reduced motion artifacts, increased early diagnosis of small mobile tumors, and improved target definition for treatment planning. High speed cameras for video, X-ray, or other forms of sequential imaging allow a live tracking of external or internal movement useful for real-time image-guided radiation therapy (IGRT). However, none of 4DMI can track real-time organ motion and no camera has correlated with 4DMI to show volumetric changes. With a brief review of various IGRT techniques, we propose a fast 3D camera for live-video stereovision, an automatic surface-motion identifier to classify body or respiratory motion, a mechanical model for synchronizing the external surface movement with the internal target displacement by combination use of the real-time stereovision and pre-treatment 4DMI, and dynamic multi-leaf collimation for adaptive aiming the moving target. Our preliminary results demonstrate that the technique is feasible and efficient in IGRT of mobile targets. A clinical trial has been initiated for validation of its spatial and temporal accuracies and dosimetric impact for intensity-modulated RT (IMRT), volumetric-modulated arc therapy (VMAT), and stereotactic body radiotherapy (SBRT) of any mobile tumors. The technique can be extended for surface-guided stereotactic needle insertion in biopsy of small lung nodules.
Comparison of 3-D Modeling With Experimental Results on Fast Wave Antenna Loading in DIII-D
NASA Astrophysics Data System (ADS)
Pinsker, R. I.; Ryan, P. M.; Goulding, R. H.; Hanson, G. R.; Milanesio, D.; Maggiora, R.; Hosea, J. C.; Nagy, A.; Porkolab, M.; Zeng, L.
2011-10-01
In DIII-D and other tokamaks, with a fixed system voltage limit, the parameter that limits the ICRF power that can be coupled to H-mode plasmas is the antenna loading resistance RL. For a fixed antenna geometry and excitation (phasing), RL is determined by the electron density profile in the antenna near-field region. Quantitative understanding of the coupling physics is obtained by comparing the resistive (RL) and reactive components of the antenna loading, without and with plasma, to predictions of 3-D models of the antenna and the edge plasma (Microwave Studio and TOPICA). When measured density profiles from reflectometers are used, good agreement between predicted and measured values of RL is obtained without any adjustable parameters in the model. The improved understanding is applied to enhancement of RL in advanced scenarios in DIII-D to increase the coupled fast wave power. Supported in part by US DOE under DE-FC02-04ER54698, DE-AC05-00OR22725, DE-AC02-09CH11466, DE-FG02-08ER54984.
NASA Astrophysics Data System (ADS)
Windhari, Ayuty; Handayani, Gunawan
2015-04-01
The 3D inversion gravity anomaly to estimate topographical density using a matlab source code from gridded data provided by Parker Oldenburg algorithm based on fast Fourier transform was computed. We extend and improved the source code of 3DINVERT.M invented by Gomez Ortiz and Agarwal (2005) using the relationship between Fourier transform of the gravity anomaly and the sum of the Fourier transform from the topography density. We gave density contrast between the two media to apply the inversion. FFT routine was implemented to construct amplitude spectrum to the given mean depth. The results were presented as new graphics of inverted topography density, the gravity anomaly due to the inverted topography and the difference between the input gravity data and the computed ones. It terminates when the RMS error is lower than pre-assigned value used as convergence criterion or until maximum of iterations is reached. As an example, we used the matlab program on gravity data of Banten region, Indonesia.
2016-01-01
Objective. Our objective is to develop a computerized scheme for liver tumor segmentation in MR images. Materials and Methods. Our proposed scheme consists of four main stages. Firstly, the region of interest (ROI) image which contains the liver tumor region in the T1-weighted MR image series was extracted by using seed points. The noise in this ROI image was reduced and the boundaries were enhanced. A 3D fast marching algorithm was applied to generate the initial labeled regions which are considered as teacher regions. A single hidden layer feedforward neural network (SLFN), which was trained by a noniterative algorithm, was employed to classify the unlabeled voxels. Finally, the postprocessing stage was applied to extract and refine the liver tumor boundaries. The liver tumors determined by our scheme were compared with those manually traced by a radiologist, used as the “ground truth.” Results. The study was evaluated on two datasets of 25 tumors from 16 patients. The proposed scheme obtained the mean volumetric overlap error of 27.43% and the mean percentage volume error of 15.73%. The mean of the average surface distance, the root mean square surface distance, and the maximal surface distance were 0.58 mm, 1.20 mm, and 6.29 mm, respectively. PMID:27597960
Le, Trong-Ngoc; Bao, Pham The; Huynh, Hieu Trung
2016-01-01
Objective. Our objective is to develop a computerized scheme for liver tumor segmentation in MR images. Materials and Methods. Our proposed scheme consists of four main stages. Firstly, the region of interest (ROI) image which contains the liver tumor region in the T1-weighted MR image series was extracted by using seed points. The noise in this ROI image was reduced and the boundaries were enhanced. A 3D fast marching algorithm was applied to generate the initial labeled regions which are considered as teacher regions. A single hidden layer feedforward neural network (SLFN), which was trained by a noniterative algorithm, was employed to classify the unlabeled voxels. Finally, the postprocessing stage was applied to extract and refine the liver tumor boundaries. The liver tumors determined by our scheme were compared with those manually traced by a radiologist, used as the "ground truth." Results. The study was evaluated on two datasets of 25 tumors from 16 patients. The proposed scheme obtained the mean volumetric overlap error of 27.43% and the mean percentage volume error of 15.73%. The mean of the average surface distance, the root mean square surface distance, and the maximal surface distance were 0.58 mm, 1.20 mm, and 6.29 mm, respectively. PMID:27597960
NASA Astrophysics Data System (ADS)
Chang, Chenliang; Xia, Jun; Lei, Wei
2015-03-01
We proposed a new method to calculate the color computer generated hologram of three-dimensional object in holographic display. The three-dimensional object is composed of several tilted planes which are tilted from the hologram. The diffraction from each tilted plane to the hologram plane is calculated based on the coordinate rotation in Fourier spectrum domains. We used the nonuniform fast Fourier transformation (NUFFT) to calculate the nonuniform sampled Fourier spectrum on the tilted plane after coordinate rotation. By using the NUFFT, the diffraction calculation from tilted plane to the hologram plane with variable sampling rates can be achieved, which overcomes the sampling restriction of FFT in the conventional angular spectrum based method. The holograms of red, green and blue component of the polygon-based object are calculated separately by using our NUFFT based method. Then the color hologram is synthesized by placing the red, green and blue component hologram in sequence. The chromatic aberration caused by the wavelength difference can be solved effectively by restricting the sampling rate of the object in the calculation of each wavelength component. The computer simulation shows the feasibility of our method in calculating the color hologram of polygon-based object. The 3D object can be displayed in color with adjustable size and no chromatic aberration in holographic display system, which can be considered as an important application in the colorful holographic three-dimensional display.
Nugroho, Hendro; Widiyantoro, Sri; Nugraha, Andri Dian
2013-09-09
Determination of earthquake hypocenter in Indonesia conducted by the Meteorological, Climatological, and Geophysical Agency (MCGA) has still used a 1-D seismic velocity model. In this research, we have applied a Fast Grid Search (FGM) method and a 3-D velocity model resulting from tomographic imaging to relocate earthquakes in the Sumatran region. The data were taken from the MCGA data catalog from 2009 to 2011 comprising of subduction zone and on land fault earthquakes with magnitude greater than 4 Mw. Our preliminary results show some significant changes in the depths of the relocated earthquakes which are in general deeper than the depths of hypocenters from the MCGA data catalog. The residual times resulting from the relocation process are smaller than those prior to the relocation. Encouraged by these results, we will continue to conduct hypocenter relocation for all events from the MCGA data catalog periodically in order to produce a new data catalog with good quality. We hope that the new data catalog will be useful for further studies.
MO-G-BRF-07: Anomalously Fast Diffusion of Carbon Nanotubes Carriers in 3D Tissue Model
Wang, Y; Bahng, J; Kotov, N
2014-06-15
Purpose: We aim to investigate and understand diffusion process of carbon nanotubes (CNTs) and other nanoscale particles in tissue and organs. Methods: In this research, we utilized a 3D model tissue of hepatocellular carcinoma (HCC)cultured in inverted colloidal crystal (ICC) scaffolds to compare the diffusivity of CNTs with small molecules such as Rhodamine and FITC in vitro, and further investigated the transportation of CNTs with and without targeting ligand, TGFβ1. The real-time permeation profiles of CNTs in HCC tissue model with high temporal and spatial resolution was demonstrated by using standard confocal microscopy. Quantitative analysis of the diffusion process in 3D was carried out using luminescence intensity in a series of Z-stack images obtained for different time points of the diffusion process after initial addition of CNTs or small molecules to the cell culture and the image data was analyzed by software ImageJ and Mathematica. Results: CNTs display diffusion rate in model tissues substantially faster than small molecules of the similar charge such as FITC, and the diffusion rate of CNTs are significantly enhanced with targeting ligand, TGFβ1. Conclusion: In terms of the advantages of in-vitro model, we were able to have access to measuring the rate of CNT penetration at designed conditions with variable parameters. And the findings by using this model, changed our understanding about advantages of CNTs as nanoscale drug carriers and provides design principles for making new drug carriers for both treatment and diagnostics. Additionally the fast diffusion opens the discussion of the best possible drug carriers to reach deep parts of cancerous tissues, which is often a prerequisite for successful cancer treatment. This work was supported by the Center for Photonic and Multiscale Nanomaterials funded by National Science Foundation Materials Research Science and Engineering Center program DMR 1120923. The work was also partially supported by NSF
NASA Astrophysics Data System (ADS)
McLaughlin, J. A.; De Moortel, I.; Hood, A. W.; Brady, C. S.
2009-01-01
Context: This paper extends the models of Craig & McClymont (1991, ApJ, 371, L41) and McLaughlin & Hood (2004, A&A, 420, 1129) to include finite β and nonlinear effects. Aims: We investigate the nature of nonlinear fast magnetoacoustic waves about a 2D magnetic X-point. Methods: We solve the compressible and resistive MHD equations using a Lagrangian remap, shock capturing code (Arber et al. 2001, J. Comp. Phys., 171, 151) and consider an initial condition in {v}×{B} \\cdot {hat{z}} (a natural variable of the system). Results: We observe the formation of both fast and slow oblique magnetic shocks. The nonlinear wave deforms the X-point into a “cusp-like” point which in turn collapses to a current sheet. The system then evolves through a series of horizontal and vertical current sheets, with associated changes in connectivity, i.e. the system exhibits oscillatory reconnection. Our final state is non-potential (but in force balance) due to asymmetric heating from the shocks. Larger amplitudes in our initial condition correspond to larger values of the final current density left in the system. Conclusions: The inclusion of nonlinear terms introduces several new features to the system that were absent from the linear regime. A movie is available in electronic form at http://www.aanda.org
Radiative tearing - Magnetic reconnection on a fast thermal-instability time scale
NASA Technical Reports Server (NTRS)
Steinolfson, R. S.; Van Hoven, G.
1984-01-01
Two energy modification mechanisms which are known to occur in sheared magnetic fields are the tearing and thermal instabilities. These processes can be studied separately with formalisms incorporating just the effective driving mechanism of interest (finite resistivity for the tearing mode and unstable radiation for the thermal mode). A model which includes both effects, and a temperature-dependent resistivity, indicates that modified forms of these two instabilities may coexist for identical physical conditions. When they are isolated computationally, one can show that their limiting growth rates are approximately those of the uncoupled instabilities. The spatial structure and energy content of these two new hybrid processes are then individually examined and are found to differ considerably from those obtained from separate treatments of the driving mechanisms. The faster radiative instability, which has a hydromagnetically scaled growth rate like the condensation mode of the thermal instability, is shown to involve a substantial amount of magnetic field reconnection. This can be partially explained by a large temperature drop (or resistivity rise) at the X-point. The island width of the Coulomb-coupled radiative mode is 30 percent of that produced by a comparable level of the slower tearing instability. In addition, the perturbed magnetic energy in the radiative instability is 5 times that of the perturbed thermal energy, indicating an appreciable modification of the initial magnetic structure.
NASA Astrophysics Data System (ADS)
Canavesi, Cristina; Cogliati, Andrea; Hayes, Adam; Santhanam, Anand P.; Tankam, Patrice; Rolland, Jannick P.
2015-10-01
Fast, robust, nondestructive 3D imaging is needed for characterization of microscopic structures in industrial and clinical applications. A custom micro-electromechanical system (MEMS)-based 2D scanner system was developed to achieve 55 kHz A-scan acquisition in a Gabor-domain optical coherence microscopy (GD-OCM) instrument with a novel multilevel GPU architecture for high-speed imaging. GD-OCM yields high-definition volumetric imaging with dynamic depth of focusing through a bio-inspired liquid lens-based microscope design, which has no moving parts and is suitable for use in a manufacturing setting or in a medical environment. A dual-axis MEMS mirror was chosen to replace two single-axis galvanometer mirrors; as a result, the astigmatism caused by the mismatch between the optical pupil and the scanning location was eliminated and a 12x reduction in volume of the scanning system was achieved. Imaging at an invariant resolution of 2 μm was demonstrated throughout a volume of 1 × 1 × 0.6 mm3, acquired in less than 2 minutes. The MEMS-based scanner resulted in improved image quality, increased robustness and lighter weight of the system - all factors that are critical for on-field deployment. A custom integrated feedback system consisting of a laser diode and a position-sensing detector was developed to investigate the impact of the resonant frequency of the MEMS and the driving signal of the scanner on the movement of the mirror. Results on the metrology of manufactured materials and characterization of tissue samples with GD-OCM are presented.
None, None
2015-09-28
Coulomb interaction between charged particles inside a bunch is one of the most importance collective effects in beam dynamics, becoming even more significant as the energy of the particle beam is lowered to accommodate analytical and low-Z material imaging purposes such as in the time resolved Ultrafast Electron Microscope (UEM) development currently underway at Michigan State University. In addition, space charge effects are the key limiting factor in the development of ultrafast atomic resolution electron imaging and diffraction technologies and are also correlated with an irreversible growth in rms beam emittance due to fluctuating components of the nonlinear electron dynamics. In the short pulse regime used in the UEM, space charge effects also lead to virtual cathode formation in which the negative charge of the electrons emitted at earlier times, combined with the attractive surface field, hinders further emission of particles and causes a degradation of the pulse properties. Space charge and virtual cathode effects and their remediation are core issues for the development of the next generation of high-brightness UEMs. Since the analytical models are only applicable for special cases, numerical simulations, in addition to experiments, are usually necessary to accurately understand the space charge effect. In this paper we will introduce a grid-free differential algebra based multiple level fast multipole algorithm, which calculates the 3D space charge field for n charged particles in arbitrary distribution with an efficiency of O(n), and the implementation of the algorithm to a simulation code for space charge dominated photoemission processes.
None, None
2015-09-28
Coulomb interaction between charged particles inside a bunch is one of the most importance collective effects in beam dynamics, becoming even more significant as the energy of the particle beam is lowered to accommodate analytical and low-Z material imaging purposes such as in the time resolved Ultrafast Electron Microscope (UEM) development currently underway at Michigan State University. In addition, space charge effects are the key limiting factor in the development of ultrafast atomic resolution electron imaging and diffraction technologies and are also correlated with an irreversible growth in rms beam emittance due to fluctuating components of the nonlinear electron dynamics.more » In the short pulse regime used in the UEM, space charge effects also lead to virtual cathode formation in which the negative charge of the electrons emitted at earlier times, combined with the attractive surface field, hinders further emission of particles and causes a degradation of the pulse properties. Space charge and virtual cathode effects and their remediation are core issues for the development of the next generation of high-brightness UEMs. Since the analytical models are only applicable for special cases, numerical simulations, in addition to experiments, are usually necessary to accurately understand the space charge effect. In this paper we will introduce a grid-free differential algebra based multiple level fast multipole algorithm, which calculates the 3D space charge field for n charged particles in arbitrary distribution with an efficiency of O(n), and the implementation of the algorithm to a simulation code for space charge dominated photoemission processes.« less
Experimental Investigation of the Trigger Problem in Magnetic Reconnection
NASA Astrophysics Data System (ADS)
Egedal, Jan
2012-07-01
Magnetic reconnection releases magnetic energy not only in steady-state, but also in time-dependent and often explosive events requiring a transition from slow reconnection to fast. The question of what causes this transition is known as the ``trigger problem'' and is not well understood. We address the trigger problem using the Versatile Toroidal Facility (VTF) at MIT. We observe spontaneous reconnection events [1] with exponentially growing reconnection rates, and we characterize the 3D dynamics of these events using multiple internal probes. The observed reconnection is asymmetric: it begins at one toroidal location and propagates around in both directions. The spontaneous onset is facilitated by an interaction between the x-line current channel and a global mode in the electrostatic potential. It is this mode which breaks axisymmetry and enables a localized decrease in x-line current. We model the onset using an empirical Ohm's law and current continuity, which is maintained by ion polarization currents associated with the mode. The model reproduces the exponential growth of the reconnection electric field, and the model growth rate agrees well with the experimentally measured growth rate. The onset location is likely determined by small asymmetries in the in-vessel coils. To further investigate this conjecture new coils have been installed, which allow for controlled changes in toroidal asymmetry. The observations are suggestive of solar flare dynamics and are relevant to tokamak research [3]. [1] Egedal, J. et al. Laboratory Observations of Spontaneous Magnetic Reconnection. Phys. Rev. Lett. 98, 015003 (2007). [2] Katz, N. et al. Laboratory Observation of Localized Onset of Magnetic Reconnection. Phys. Rev. Lett. 104, 255004 (2010). [3] Park, H.K. et al. Self-organized Te redistribution during driven reconnection processes in high-temperature plasmas. Phys. Plasmas 13, 055907 (2006).
Ai, H; Pan, T; Hwang, K
2014-06-15
Purpose: To determine the feasibility of identifying cortical bone on MR images with a short-TE 3D fast-GRE sequence for attenuation correction of PET data in PET/MR. Methods: A water-fat-bone phantom was constructed with two pieces of beef shank. MR scans were performed on a 3T MR scanner (GE Discovery™ MR750). A 3D GRE sequence was first employed to measure the level of residual signal in cortical bone (TE{sub 1}/TE{sub 2}/TE{sub 3}=2.2/4.4/6.6ms, TR=20ms, flip angle=25°). For cortical bone segmentation, a 3D fast-GRE sequence (TE/TR=0.7/1.9ms, acquisition voxel size=2.5×2.5×3mm{sup 3}) was implemented along with a 3D Dixon sequence (TE{sub 1}/TE{sub 2}/TR=1.2/2.3/4.0ms, acquisition voxel size=1.25×1.25×3mm{sup 3}) for water/fat imaging. Flip angle (10°), acquisition bandwidth (250kHz), FOV (480×480×144mm{sup 3}) and reconstructed voxel size (0.94×0.94×1.5mm{sup 3}) were kept the same for both sequences. Soft tissue and fat tissue were first segmented on the reconstructed water/fat image. A tissue mask was created by combining the segmented water/fat masks, which was then applied on the fast-GRE image (MRFGRE). A second mask was created to remove the Gibbs artifacts present in regions in close vicinity to the phantom. MRFGRE data was smoothed with a 3D anisotropic diffusion filter for noise reduction, after which cortical bone and air was separated using a threshold determined from the histogram. Results: There is signal in the cortical bone region in the 3D GRE images, indicating the possibility of separating cortical bone and air based on signal intensity from short-TE MR image. The acquisition time for the 3D fast-GRE sequence was 17s, which can be reduced to less than 10s with parallel imaging. The attenuation image created from water-fat-bone segmentation is visually similar compared to reference CT. Conclusion: Cortical bone and air can be separated based on intensity in MR image with a short-TE 3D fast-GRE sequence. Further research is required
The polyGeVero® software for fast and easy computation of 3D radiotherapy dosimetry data
NASA Astrophysics Data System (ADS)
Kozicki, Marek; Maras, Piotr
2015-01-01
The polyGeVero® software package was elaborated for calculations of 3D dosimetry data such as the polymer gel dosimetry. It comprises four workspaces designed for: i) calculating calibrations, ii) storing calibrations in a database, iii) calculating dose distribution 3D cubes, iv) comparing two datasets e.g. a measured one with a 3D dosimetry with a calculated one with the aid of a treatment planning system. To accomplish calculations the software was equipped with a number of tools such as the brachytherapy isotopes database, brachytherapy dose versus distance calculation based on the line approximation approach, automatic spatial alignment of two 3D dose cubes for comparison purposes, 3D gamma index, 3D gamma angle, 3D dose difference, Pearson's coefficient, histograms calculations, isodoses superimposition for two datasets, and profiles calculations in any desired direction. This communication is to briefly present the main functions of the software and report on the speed of calculations performed by polyGeVero®.
3D multiple optical tweezers based on time-shared scanning with a fast focus tunable lens
NASA Astrophysics Data System (ADS)
Tanaka, Yoshio
2013-02-01
Three-dimensional controlled manipulation of individual micro-objects requires multiple optical tweezers that can be independently controlled in a 3D working space with high spatiotemporal resolution. Here, the author presents 3D multiple optical tweezers based on a time-shared scanning technique with an electrically focus tunable lens for axial steering and a two-axis steering mirror for lateral steering. Four typical examples of 3D controlled manipulation, including the rotation of a single bead on its axis, are demonstrated in real time. The optical system design and the control method are also described.
NASA Astrophysics Data System (ADS)
Jo, Hong Li; Song, Yo Han; Park, Jinho; Jo, Eun-Jung; Goh, Yeongchang; Shin, Kyujin; Kim, Min-Gon; Lee, Kang Taek
2015-11-01
We report on the development of a three-dimensional (3D) live-cell imaging technique with high spatiotemporal resolution using lanthanide-doped upconverting nanoparticles (UCNPs). It employs the sectioning capability of confocal microscopy except that the two-dimensional (2D) section images are acquired by wide-field epi-fluorescence microscopy. Although epi-fluorescence images are contaminated with the out-of-focus background in general, the near-infrared (NIR) excitation used for the excitation of UCNPs does not generate any autofluorescence, which helps to lower the background. Moreover, the image blurring due to defocusing was naturally eliminated in the image reconstruction process. The 3D images were used to investigate the cellular dynamics such as nuclear uptake and single-particle tracking that require 3D description.We report on the development of a three-dimensional (3D) live-cell imaging technique with high spatiotemporal resolution using lanthanide-doped upconverting nanoparticles (UCNPs). It employs the sectioning capability of confocal microscopy except that the two-dimensional (2D) section images are acquired by wide-field epi-fluorescence microscopy. Although epi-fluorescence images are contaminated with the out-of-focus background in general, the near-infrared (NIR) excitation used for the excitation of UCNPs does not generate any autofluorescence, which helps to lower the background. Moreover, the image blurring due to defocusing was naturally eliminated in the image reconstruction process. The 3D images were used to investigate the cellular dynamics such as nuclear uptake and single-particle tracking that require 3D description. Electronic supplementary information (ESI) available: Methods for synthetic procedure, characterization, live-cell imaging with UCNPs, and two video files for 3D imaging. See DOI: 10.1039/c5nr05875a
3D Ion Temperature Reconstruction
NASA Astrophysics Data System (ADS)
Tanabe, Hiroshi; You, Setthivoine; Balandin, Alexander; Inomoto, Michiaki; Ono, Yasushi
2009-11-01
The TS-4 experiment at the University of Tokyo collides two spheromaks to form a single high-beta compact toroid. Magnetic reconnection during the merging process heats and accelerates the plasma in toroidal and poloidal directions. The reconnection region has a complex 3D topology determined by the pitch of the spheromak magnetic fields at the merging plane. A pair of multichord passive spectroscopic diagnostics have been established to measure the ion temperature and velocity in the reconnection volume. One setup measures spectral lines across a poloidal plane, retrieving velocity and temperature from Abel inversion. The other, novel setup records spectral lines across another section of the plasma and reconstructs velocity and temperature from 3D vector and 2D scalar tomography techniques. The magnetic field linking both measurement planes is determined from in situ magnetic probe arrays. The ion temperature is then estimated within the volume between the two measurement planes and at the reconnection region. The measurement is followed over several repeatable discharges to follow the heating and acceleration process during the merging reconnection.
Rupprecht, Peter; Prendergast, Andrew; Wyart, Claire; Friedrich, Rainer W
2016-01-01
There is a high demand for 3D multiphoton imaging in neuroscience and other fields but scanning in axial direction presents technical challenges. We developed a focusing technique based on a remote movable mirror that is conjugate to the specimen plane and translated by a voice coil motor. We constructed cost-effective z-scanning modules from off-the-shelf components that can be mounted onto standard multiphoton laser scanning microscopes to extend scan patterns from 2D to 3D. Systems were designed for large objectives and provide high resolution, high speed and a large z-scan range (>300 μm). We used these systems for 3D multiphoton calcium imaging in the adult zebrafish brain and measured odor-evoked activity patterns across >1500 neurons with single-neuron resolution and high signal-to-noise ratio. PMID:27231612
Rupprecht, Peter; Prendergast, Andrew; Wyart, Claire; Friedrich, Rainer W
2016-05-01
There is a high demand for 3D multiphoton imaging in neuroscience and other fields but scanning in axial direction presents technical challenges. We developed a focusing technique based on a remote movable mirror that is conjugate to the specimen plane and translated by a voice coil motor. We constructed cost-effective z-scanning modules from off-the-shelf components that can be mounted onto standard multiphoton laser scanning microscopes to extend scan patterns from 2D to 3D. Systems were designed for large objectives and provide high resolution, high speed and a large z-scan range (>300 μm). We used these systems for 3D multiphoton calcium imaging in the adult zebrafish brain and measured odor-evoked activity patterns across >1500 neurons with single-neuron resolution and high signal-to-noise ratio. PMID:27231612
David L. Alumbaugh
2006-07-03
This report addresses the effects of electrical anisotropy on the 3D inversion of single-well induction logging data when anisotropy is not considered. Of concern are possible artifacts that may lead to an incorrect interpretation of the formation about the borehole. Comparison is made of 3D isotropic inversion on a suite of model data, with and without anisotropy, consisting of an infinite layer and layer terminated at the borehole. In both cases, the layer dip (or well deviation) is varied. Inversion of the anisotropic data result in an overestimate of the layer conductivity, and the lateral extent of the layer about the borehole.
NASA Astrophysics Data System (ADS)
Zhang, M. Q.
1989-09-01
A new Monte Carlo algorithm for 3D Kawasaki spin-exchange simulations and its implementation on a CDC CYBER 205 is presented. This approach is applicable to lattices with sizes between 4×4×4 and 256×L2×L3 ((L2+2)(L3+4)/4⩽65535) and periodic boundary conditions. It is adjustable to various kinetic models in which the total magnetization is conserved. Maximum speed on 10 million steps per second can be reached for 3-D Ising model with Metropolis rate.
NASA Astrophysics Data System (ADS)
Stover, E.; Berger, G.; Wendel, M.; Petter, J.
2015-10-01
A method for non-contact 3D form testing of aspheric surfaces including determination of decenter and wedge errors and lens thickness is presented. The principle is based on the absolute measurement capability of multi-wavelength interferometry (MWLI). The approach produces high density 3D shape information and geometric parameters at high accuracy in short measurement times. The system allows inspection of aspheres without restrictions in terms of spherical departures, of segmented and discontinuous optics. The optics can be polished or ground and made of opaque or transparent materials.
VIRO 3D: fast three-dimensional full-body scanning for humans and other living objects
NASA Astrophysics Data System (ADS)
Stein, Norbert; Minge, Bernhard
1998-03-01
The development of a family of partial and whole body scanners provides a complete technology for fully three-dimensional and contact-free scans on human bodies or other living objects within seconds. This paper gives insight into the design and the functional principles of the whole body scanner VIRO 3D operating on the basis of the laser split-beam method. The arrangement of up to 24 camera/laser combinations, thus dividing the area into different camera fields and an all- around sensor configuration travelling in vertical direction allow the complete 360-degree-scan of an object within 6 - 20 seconds. Due to a special calibration process the different sensors are matched and the measured data are combined. Up to 10 million 3D measuring points with a resolution of approximately 1 mm are processed in all coordinate axes to generate a 3D model. By means of high-performance processors in combination with real-time image processing chips the image data from almost any number of sensors can be recorded and evaluated synchronously in video real-time. VIRO 3D scanning systems have already been successfully implemented in various applications and will open up new perspectives in different other fields, ranging from industry, orthopaedic medicine, plastic surgery to art and photography.
3D in the Fast Lane: Render as You Go with the Latest OpenGL Boards.
ERIC Educational Resources Information Center
Sauer, Jeff; Murphy, Sam
1997-01-01
NT OpenGL hardware allows modelers and animators to work at relatively inexpensive NT workstations in their own offices or homes previous to shared space and workstation time in expensive studios. Rates seven OpenGL boards and two QuickDraw 3D accelerator boards for Mac users on overall value, wireframe and texture rendering, 2D acceleration, and…
NASA Astrophysics Data System (ADS)
Bates, Jason; Schmitt, Andrew; Klapisch, Marcel; Karasik, Max; Obenschain, Steve
2013-10-01
Modifications to the FAST3D code have been made to enhance its ability to simulate the dynamics of plastic ICF targets with high-Z overcoats. This class of problems is challenging computationally due in part to plasma conditions that are not in a state of local thermodynamic equilibrium and to the presence of mixed computational cells containing more than one material. Recently, new opacity tables for gold, palladium and plastic have been generated with an improved version of the STA code. These improved tables provide smoother, higher-fidelity opacity data over a wider range of temperature and density states than before, and contribute to a more accurate treatment of radiative transfer processes in FAST3D simulations. Furthermore, a new, more efficient subroutine known as ``MMEOS'' has been installed in the FAST3D code for determining pressure and temperature equilibrium conditions within cells containing multiple materials. We will discuss these topics, and present new simulation results for high-Z planar-target experiments performed recently on the NIKE Laser Facility. Work supported by DOE/NNSA.
3D pyCloudy modelling of bipolar planetary nebulae: Evidence for fast fading of the lobes
NASA Astrophysics Data System (ADS)
Gesicki, K.; Zijlstra, A. A.; Morisset, C.
2016-01-01
Aims: The origin and evolution of the shapes of bipolar planetary nebulae are poorly understood. We postulate that their history can be traced through their internal velocity fields in a procedure similar to the one well established for spherical objects. Such an analysis requires 3D photoionization and kinematical modelling that is computationally very time consuming. We apply an axially symmetric pseudo-3D photoionization model, pyCloudy, to derive the structures of six bipolar nebulae and two suggested post-bipolars in an attempt to constrain the bipolar planetary nebulae evolution. Methods: HST images and VLT/UVES spectroscopy are used for the modelling. The targets are located in the direction of the Galactic bulge. A 3D model structure is used as input to the photoionization code in order to fit the HST images. Line profiles of different ions constrain the velocity field. The model and associated velocity fields allow us to derive masses, velocities, and ages. Results: The 3D models find much lower ionized masses than required in 1D models: ionized masses are reduced by factors of 2-7. The selected bi-lobed planetary nebulae show a narrow range of ages: the averaged radii and velocities result in values between 1300 and 2000 yr. The lobes are fitted well with velocities that increase linearly with radius. These Hubble-type flows have been found before, and suggest that the lobes form at a defined point in time. The lobes appear to be slightly younger, by ~500 yr, than the main (host) nebulae; they seem to form at an early phase of PN evolution and fade after 1-2 kyr. We find that 30-35% of bulge PNe pass through a bipolar phase. Based on observations collected at the European Organisation for Astronomical Research in the Southern Hemisphere, Chile (proposal 075.D-0104) and HST (program 9356).
Simulating Magnetotail Reconnection in the Presence of a Significant Guide Field
NASA Astrophysics Data System (ADS)
Ashour-Abdalla, Maha; Lapenta, Giovanni; Walker, Raymond; El-Alaoui, Mostafa; Liang, Haoming; Zhou, Meng; Berchem, Jean; Goldstein, Melvyn
2016-04-01
For a large number of solar wind driven substorms a significant guide field is present in the magnetotail. We have developed a new simulation approach for magnetotail reconnection in which we combine results from global magnetohydrodynamic simulations of the solar wind, magnetosphere and ionosphere system with large-scale 3D particle-in-cell simulations of the magnetotail. In this approach we first run a global MHD simulation driven by upstream solar wind observations. Then we use the results from the MHD simulation to set the initial and boundary conditions of a large-scale PIC simulation. We have selected the iPIC3D kinetic code to go along with the UCLA global simulation model. Our approach allows us to investigate indicators of reconnection in a realistic magnetospheric configuration, which is essential to determine the dynamics of the substorm process, in particular the evolution of reconnection and dipolarization fronts. In preparation for the MMS phase in the magnetotail, we have developed an analysis system with multiprobes so we can directly compare the simulation results with MMS observations and have applied our multiscale simulation approach to the solar wind driven substorm observed on March 1, 2008. This substorm was characterized by a significant guide field. In the PIC simulation fast but episodic reconnection is found in the near-Earth tail. The reconnection has a complex 3D structure. When normalized by the Alfvén speed and magnetic field magnitude, the reconnection rate although highly variable oscillates about the expected value of 0.1. The location of reconnection also is highly variable. In the simulation it moves several RE across the tail in less than two minutes. We have evaluated several indicators of reconnection and find that agyrotropy (or non-gyrotropy) is the most reliable indicator.
Fast multi-core based multimodal registration of 2D cross-sections and 3D datasets
2010-01-01
Background Solving bioinformatics tasks often requires extensive computational power. Recent trends in processor architecture combine multiple cores into a single chip to improve overall performance. The Cell Broadband Engine (CBE), a heterogeneous multi-core processor, provides power-efficient and cost-effective high-performance computing. One application area is image analysis and visualisation, in particular registration of 2D cross-sections into 3D image datasets. Such techniques can be used to put different image modalities into spatial correspondence, for example, 2D images of histological cuts into morphological 3D frameworks. Results We evaluate the CBE-driven PlayStation 3 as a high performance, cost-effective computing platform by adapting a multimodal alignment procedure to several characteristic hardware properties. The optimisations are based on partitioning, vectorisation, branch reducing and loop unrolling techniques with special attention to 32-bit multiplies and limited local storage on the computing units. We show how a typical image analysis and visualisation problem, the multimodal registration of 2D cross-sections and 3D datasets, benefits from the multi-core based implementation of the alignment algorithm. We discuss several CBE-based optimisation methods and compare our results to standard solutions. More information and the source code are available from http://cbe.ipk-gatersleben.de. Conclusions The results demonstrate that the CBE processor in a PlayStation 3 accelerates computational intensive multimodal registration, which is of great importance in biological/medical image processing. The PlayStation 3 as a low cost CBE-based platform offers an efficient option to conventional hardware to solve computational problems in image processing and bioinformatics. PMID:20064262
Bao Yidong; Hu Sibo; Lang Zhikui; Hu Ping
2005-08-05
A fast simulation scheme for 3D curved binder flanging and blank shape prediction of sheet metal based on one-step inverse finite element method is proposed, in which the total plasticity theory and proportional loading assumption are used. The scheme can be actually used to simulate 3D flanging with complex curve binder shape, and suitable for simulating any type of flanging model by numerically determining the flanging height and flanging lines. Compared with other methods such as analytic algorithm and blank sheet-cut return method, the prominent advantage of the present scheme is that it can directly predict the location of the 3D flanging lines when simulating the flanging process. Therefore, the prediction time of flanging lines will be obviously decreased. Two typical 3D curve binder flanging including stretch and shrink characters are simulated in the same time by using the present scheme and incremental FE non-inverse algorithm based on incremental plasticity theory, which show the validity and high efficiency of the present scheme.
NASA Astrophysics Data System (ADS)
Huang, Guo-Jiao; Bai, Chao-Ying; Greenhalgh, Stewart
2013-09-01
The traditional grid/cell-based wavefront expansion algorithms, such as the shortest path algorithm, can only find the first arrivals or multiply reflected (or mode converted) waves transmitted from subsurface interfaces, but cannot calculate the other later reflections/conversions having a minimax time path. In order to overcome the above limitations, we introduce the concept of a stationary minimax time path of Fermat's Principle into the multistage irregular shortest path method. Here we extend it from Cartesian coordinates for a flat earth model to global ray tracing of multiple phases in a 3-D complex spherical earth model. The ray tracing results for 49 different kinds of crustal, mantle and core phases show that the maximum absolute traveltime error is less than 0.12 s and the average absolute traveltime error is within 0.09 s when compared with the AK135 theoretical traveltime tables for a 1-D reference model. Numerical tests in terms of computational accuracy and CPU time consumption indicate that the new scheme is an accurate, efficient and a practical way to perform 3-D multiphase arrival tracking in regional or global traveltime tomography.
Magnetic reconnection: from the Sweet–Parker model to stochastic plasmoid chains
NASA Astrophysics Data System (ADS)
Loureiro, N. F.; Uzdensky, D. A.
2016-01-01
Magnetic reconnection is the topological reconfiguration of the magnetic field in a plasma, accompanied by the violent release of energy and particle acceleration. Reconnection is as ubiquitous as plasmas themselves, with solar flares perhaps the most popular example. Other fascinating processes where reconnection plays a key role include the magnetic dynamo, geomagnetic storms and the sawtooth crash in tokamaks. Over the last few years, the theoretical understanding of magnetic reconnection in large-scale fluid systems has undergone a major paradigm shift. The steady-state model of reconnection described by the famous Sweet–Parker (SP) theory, which dominated the field for ∼50 years, has been replaced with an essentially time-dependent, bursty picture of the reconnection layer, dominated by the continuous formation and ejection of multiple secondary islands (plasmoids). Whereas in the SP model reconnection was predicted to be slow, a major implication of this new paradigm is that reconnection in fluid systems is fast (i.e. independent of the Lundquist number), provided that the system is large enough. This conceptual shift hinges on the realization that SP-like current layers are violently unstable to the plasmoid (tearing) instability—implying, therefore, that such current sheets are super-critically unstable and thus can never form in the first place. This suggests that the formation of a current sheet and the subsequent reconnection process cannot be decoupled, as is commonly assumed. This paper provides an introductory-level overview of the recent developments in reconnection theory and simulations that led to this essentially new framework. We briefly discuss the role played by the plasmoid instability in selected applications, and describe some of the outstanding challenges that remain at the frontier of this subject. Amongst these are the analytical and numerical extension of the plasmoid instability to (i) 3D and (ii) non
Fluxon Modeling of Eruptive Events With and Without Reconnection
NASA Astrophysics Data System (ADS)
DeForest, Craig; Rachmeler, L.; Davey, A.; Kankelborg, C.
2007-05-01
Fluxon MHD models represent the coronal magnetic field as a "skeleton" of discretized field lines. This quasi-Lagrangian approach eliminates numerical resistivity and allows 3-D time-dependent plasma simulation in a desktop workstation.Using our fluxon code, FLUX, we have demonstrated that ideal MHD instabilities can drive fast eruptive events even in the complete absence of magnetic reconnection. The mechanism ("herniation") is probably not the main driver of fast CMEs but may be applicable to microjets, macrospicules, or other small scale events where vortical flows are present in the solar atmosphere. In this presentation, we use time-dependent simulations to demonstrate energy release in several idealized plasma systems with and without magnetic reconnection.This work was funded by NASA's LWS and SHP-SR&T programs.
NASA Astrophysics Data System (ADS)
Khochfar, Sadegh; Emsellem, Eric; Serra, Paolo; Bois, Maxime; Alatalo, Katherine; Bacon, R.; Blitz, Leo; Bournaud, Frédéric; Bureau, M.; Cappellari, Michele; Davies, Roger L.; Davis, Timothy A.; de Zeeuw, P. T.; Duc, Pierre-Alain; Krajnović, Davor; Kuntschner, Harald; Lablanche, Pierre-Yves; McDermid, Richard M.; Morganti, Raffaella; Naab, Thorsten; Oosterloo, Tom; Sarzi, Marc; Scott, Nicholas; Weijmans, Anne-Marie; Young, Lisa M.
2011-10-01
We propose a simple model for the origin of fast and slow rotator early-type galaxies (ETG) within the hierarchical Λcold dark matter (ΛCDM) scenario, that is based on the assumption that the mass fraction of stellar discs in ETGs is a proxy for the specific angular momentum expressed via λR. Within our model we reproduce the fraction of fast and slow rotators as a function of magnitude in the ATLAS3D survey, assuming that fast-rotating ETGs have at least 10 per cent of their total stellar mass in a disc component. In agreement with ATLAS3D observations we find that slow rotators are predominantly galaxies with M* > 1010.5 M⊙ contributing ˜20 per cent to the overall ETG population. We show in detail that the growth histories of fast and slow rotators are different, supporting the classification of ETGs into these two categories. Slow rotators accrete between ˜50 and 90 per cent of their stellar mass from satellites and their most massive progenitors have on average up to three major mergers during their evolution. Fast rotators in contrast accrete less than 50 per cent and have on average less than one major merger in their past. We find that the underlying physical reason for the different growth histories is the slowing down and ultimately complete shut-down of gas cooling in massive galaxies. Once cooling and associated star formation in disc stop, galaxies grow via infall from satellites. Frequent minor mergers thereby destroy existing stellar discs via violent relaxation and also tend to lower the specific angular momentum of the main stellar body, lowering λR into the slow rotator regime. On average, the last gas-rich major merger interaction in slow rotators happens at z > 1.5, followed by a series of minor mergers. These results support the idea that kinematically decoupled cores (KDC) form during gas-rich major mergers at high z followed by minor mergers, which build-up the outer layers of the remnant, and make remnants that are initially too flat
Mori, Shinichiro; Kumagai, Motoki; Miki, Kentaro; Fukuhara, Riki; Haneishi, Hideaki
2015-09-01
To improve treatment workflow, we developed a graphic processing unit (GPU)-based patient positional verification software application and integrated it into carbon-ion scanning beam treatment. Here, we evaluated the basic performance of the software. The algorithm provides 2D/3D registration matching using CT and orthogonal X-ray flat panel detector (FPD) images. The participants were 53 patients with tumors of the head and neck, prostate or lung receiving carbon-ion beam treatment. 2D/3D-ITchi-Gime (ITG) calculation accuracy was evaluated in terms of computation time and registration accuracy. Registration calculation was determined using the similarity measurement metrics gradient difference (GD), normalized mutual information (NMI), zero-mean normalized cross-correlation (ZNCC), and their combination. Registration accuracy was dependent on the particular metric used. Representative examples were determined to have target registration error (TRE) = 0.45 ± 0.23 mm and angular error (AE) = 0.35 ± 0.18° with ZNCC + GD for a head and neck tumor; TRE = 0.12 ± 0.07 mm and AE = 0.16 ± 0.07° with ZNCC for a pelvic tumor; and TRE = 1.19 ± 0.78 mm and AE = 0.83 ± 0.61° with ZNCC for lung tumor. Calculation time was less than 7.26 s.The new registration software has been successfully installed and implemented in our treatment process. We expect that it will improve both treatment workflow and treatment accuracy. PMID:26081313
Mori, Shinichiro; Kumagai, Motoki; Miki, Kentaro; Fukuhara, Riki; Haneishi, Hideaki
2015-01-01
To improve treatment workflow, we developed a graphic processing unit (GPU)-based patient positional verification software application and integrated it into carbon-ion scanning beam treatment. Here, we evaluated the basic performance of the software. The algorithm provides 2D/3D registration matching using CT and orthogonal X-ray flat panel detector (FPD) images. The participants were 53 patients with tumors of the head and neck, prostate or lung receiving carbon-ion beam treatment. 2D/3D-ITchi-Gime (ITG) calculation accuracy was evaluated in terms of computation time and registration accuracy. Registration calculation was determined using the similarity measurement metrics gradient difference (GD), normalized mutual information (NMI), zero-mean normalized cross-correlation (ZNCC), and their combination. Registration accuracy was dependent on the particular metric used. Representative examples were determined to have target registration error (TRE) = 0.45 ± 0.23 mm and angular error (AE) = 0.35 ± 0.18° with ZNCC + GD for a head and neck tumor; TRE = 0.12 ± 0.07 mm and AE = 0.16 ± 0.07° with ZNCC for a pelvic tumor; and TRE = 1.19 ± 0.78 mm and AE = 0.83 ± 0.61° with ZNCC for lung tumor. Calculation time was less than 7.26 s.The new registration software has been successfully installed and implemented in our treatment process. We expect that it will improve both treatment workflow and treatment accuracy. PMID:26081313
Bessho, Naoki; Bhattacharjee, A.
2012-05-10
Magnetic reconnection and particle acceleration in relativistic Harris sheets in low-density electron-positron plasmas with no guide field have been studied by means of two-dimensional particle-in-cell simulations. Reconnection rates are of the order of one when the background density in a Harris sheet is of the order of 1% of the density in the current sheet, which is consistent with previous results in the non-relativistic regime. It has been demonstrated that the increase of the Lorentz factors of accelerated particles significantly enhances the collisionless resistivity needed to sustain a large reconnection electric field. It is shown analytically and numerically that the energy spectrum of accelerated particles near the X-line is the product of a power law and an exponential function of energy, {gamma}{sup -1/4}exp (- a{gamma}{sup 1/2}), where {gamma} is the Lorentz factor and a is a constant. However, in the low-density regime, while the most energetic particles are produced near X-lines, many more particles are energized within magnetic islands. Particles are energized in contracting islands by multiple reflection, but the mechanism is different from Fermi acceleration in magnetic islands for magnetized particles in the presence of a guide field. In magnetic islands, strong core fields are generated and plasma beta values are reduced. As a consequence, the fire-hose instability condition is not satisfied in most of the island region, and island contraction and particle acceleration can continue. In island coalescence, reconnection between two islands can accelerate some particles, however, many particles are decelerated and cooled, which is contrary to what has been discussed in the literature on particle acceleration due to reconnection in non-relativistic hydrogen plasmas.
NASA Astrophysics Data System (ADS)
Aghdasi, Nava; Li, Yangming; Berens, Angelique; Moe, Kris S.; Bly, Randall A.; Hannaford, Blake
2015-03-01
Minimally invasive neuroendoscopic surgery provides an alternative to open craniotomy for many skull base lesions. These techniques provides a great benefit to the patient through shorter ICU stays, decreased post-operative pain and quicker return to baseline function. However, density of critical neurovascular structures at the skull base makes planning for these procedures highly complex. Furthermore, additional surgical portals are often used to improve visualization and instrument access, which adds to the complexity of pre-operative planning. Surgical approach planning is currently limited and typically involves review of 2D axial, coronal, and sagittal CT and MRI images. In addition, skull base surgeons manually change the visualization effect to review all possible approaches to the target lesion and achieve an optimal surgical plan. This cumbersome process relies heavily on surgeon experience and it does not allow for 3D visualization. In this paper, we describe a rapid pre-operative planning system for skull base surgery using the following two novel concepts: importance-based highlight and mobile portal. With this innovation, critical areas in the 3D CT model are highlighted based on segmentation results. Mobile portals allow surgeons to review multiple potential entry portals in real-time with improved visualization of critical structures located inside the pathway. To achieve this we used the following methods: (1) novel bone-only atlases were manually generated, (2) orbits and the center of the skull serve as features to quickly pre-align the patient's scan with the atlas, (3) deformable registration technique was used for fine alignment, (4) surgical importance was assigned to each voxel according to a surgical dictionary, and (5) pre-defined transfer function was applied to the processed data to highlight important structures. The proposed idea was fully implemented as independent planning software and additional
SU-E-J-251: Fast MR-Based DRR Generation Using Highly Undersampled 3D Radial Trajectories
Pereira, G; Traughber, B; Traughber, M; Hu, L; Su, K; Muzic, R
2014-06-01
Purpose: The construction of a digitally reconstructed radiograph (DRR) from a magnetic resonance image (MRI) is possible if the cortical bone signal can be acquired and separated from air and soft tissue. This may be accomplished by subtracting a long echo-time, in-phase, gradient echo (GRE) image volume from an ultra-short echo time free induction decay (FID) image to produce a bone-enhanced (BE) image that reveals cortical bone. One limitation of this approach is the length of time required for data acquisition, which can limit the quality of the DRRs due to patient and organ motion. This study aimed to significantly reduce the acquisition time without compromising DRR quality. Methods: Brain data were acquired from two volunteers using a 3T MR scanner (Ingenia, Philips Healthcare). The FID and GRE images were acquired in a single acquisition using a 3D radial readout sequence with the following parameters: TE1=0.142ms (ultra-short), TE2=2.197ms (nearly in-phase), 2*2*2mm3 isotropic voxels, 250*250*250mm3 FOV. To reduce the acquisition time, k-space was sampled at 75, 50 and 25% of a full 3D sphere . The TE2 image was subtracted from the TE1 image to generate the BE images. The BE images were used to generate DRRs using the Pinnacle treatment planning system (Philips-version 9.2). The quality of the DRRs was evaluated qualitatively by 5 board certified medical physicists for clinical usefulness. Results: The acquisition time for 75, 50 and 25% sampling schemes were 219s, 146s, and 73s, respectively, the latter of which was a four-fold reduction in scan time compared to a 300s fully-sampled acquisition. All DRRs obtained were of acceptable quality and were shown to have sufficient information for clinical 2D image matching. Conclusion: Undersampling k-space while maintaining the same range of frequency information results in significantly reduced scan time and clinically acceptable DRR image quality. Drs. B Traughber and R Muzic have research support from Philips
Jiang, Chao; Jia, Shuhai; Dong, Jun; Bao, Qingchen; Yang, Jia; Lian, Qin; Li, Dichen
2015-09-21
We propose a novel multi-frequency color-marked fringe projection profilometry approach to measure the 3D shape of objects with depth discontinuities. A digital micromirror device projector is used to project a color map consisting of a series of different-frequency color-marked fringe patterns onto the target object. We use a chromaticity curve to calculate the color change caused by the height of the object. The related algorithm to measure the height is also described in this paper. To improve the measurement accuracy, a chromaticity curve correction method is presented. This correction method greatly reduces the influence of color fluctuations and measurement error on the chromaticity curve and the calculation of the object height. The simulation and experimental results validate the utility of our method. Our method avoids the conventional phase shifting and unwrapping process, as well as the independent calculation of the object height required by existing techniques. Thus, it can be used to measure complex and dynamic objects with depth discontinuities. These advantages are particularly promising for industrial applications. PMID:26406621
NASA Astrophysics Data System (ADS)
Chixiao, Chen; Jixuan, Xiang; Huabin, Chen; Jun, Xu; Fan, Ye; Ning, Li; Junyan, Ren
2015-05-01
Asynchronous successive approximation register (SAR) analog-to-digital converters (ADC) feature high energy efficiency but medium performance. From the point of view of speed, the key bottleneck is the unit capacitor size. In this paper, a small size three-dimensional (3-D) metal—oxide—metal (MOM) capacitor is proposed. The unit capacitor has a capacitance of 1-fF. It shapes as an umbrella, which is designed for fast settling consideration. A comparison among the proposed capacitor with other 3-D MOM capacitors is also given in the paper. To demonstrate the effectiveness of the MOM capacitor, a 6-b capacitive DAC is implemented in TSMC 1P9M 65 nm LP CMOS technology. The DAC consumes a power dissipation of 0.16 mW at the rate of 100 MS/s, excluding a source-follower based output buffer. Static measurement result shows that INL is less than ±1 LSB and DNL is less than ±0.5 LSB. In addition, a 100 MS/s 9-bit SAR ADC with the proposed 3-D capacitor is simulated.
Lukas, Carsten; Hahn, Horst K; Bellenberg, Barbara; Rexilius, Jan; Schmid, Gebhard; Schimrigk, Sebastian K; Przuntek, Horst; Köster, Odo; Peitgen, Heinz-Otto
2004-11-01
Fast, reliable and easy-to-use methods to quantify brain atrophy are of increasing importance in clinical studies on neuro-degenerative diseases. Here, ILAB 4, a new volumetry software that uses a fast semi-automated 3D segmentation of thin-slice T1-weighted 3D MR images based on a modified watershed transform and an automatic histogram analysis was evaluated. It provides the cerebral volumes: whole brain, white matter, gray matter and intracranial cavity. Inter- and intra-rater reliability and scan-rescan reproducibility were excellent in measuring whole brain volumes (coefficients of variation below 0.5%) of volunteers and patients. However, gray and white matter volumes were more susceptible to image quality. High accuracy of the absolute volume results (+/-5 ml) were shown by phantom and preparation measurements. Analysis times were 6 min for processing of 128 slices. The proposed technique is reliable and highly suitable for quantitative studies of brain atrophy, e.g., in multiple sclerosis. PMID:15536555
NASA Astrophysics Data System (ADS)
Bates, Jason; Schmitt, Andrew; Zalesak, Steve
2015-11-01
The ablative Rayleigh-Taylor (RT) instability is a key factor in the performance of directly-drive inertial-confinement-fusion (ICF) targets. Although this subject has been studied for quite some time, the accurate simulation of the ablative RT instability has proven to be a challenging task for many radiation hydrodynamics codes, particularly when it comes to capturing the ablatively-stabilized region of the linear dispersion spectrum and modeling ab initio perturbations. In this poster, we present results from recent two-dimensional numerical simulations of the ablative RT instability that were performed using the Eulerian code FastRad3D at the U.S. Naval Research Laboratory. We consider both planar and spherical geometries, low and moderate-Z target materials, different laser wavelengths and where possible, compare our findings with experiment data, linearized theory and/or results from other radiation hydrodynamics codes. Overall, we find that FastRad3D is capable of simulating the ablative RT instability quite accurately, although some uncertainties/discrepancies persist. We discuss these issues, as well as some of the numerical challenges associated with modeling this class of problems. Work supported by U.S. DOE/NNSA.
Magnetohydrodynamic simulations of turbulent magnetic reconnection
Fan Quanlin; Feng Xueshang; Xiang Changqing
2004-12-01
Turbulent reconnection process in a one-dimensional current sheet is investigated by means of a two-dimensional compressible one-fluid magnetohydrodynamic simulation with spatially uniform, fixed resistivity. Turbulence is set up by adding to the sheet pinch small but finite level of broadband random-phased magnetic field components. To clarify the nonlinear spatial-temporal nature of the turbulent reconnection process the reconnection system is treated as an unforced initial value problem without any anomalous resistivity model adopted. Numerical results demonstrate the duality of turbulent reconnection, i.e., a transition from Sweet-Parker-like slow reconnection to Petschek-like fast reconnection in its nonlinear evolutionary process. The initial slow reconnection phase is characterized by many independent microreconnection events confined within the sheet region and a global reconnection rate mainly dependent on the initially added turbulence and insensitive to variations of the plasma {beta} and resistivity. The formation and amplification of the major plasmoid leads the following reconnection process to a rapid reconnection stage with a fast reconnection rate of the order of 0.1 or even larger, drastically changing the topology of the global magnetic field. That is, the presence of magnetohydrodynamic turbulence in large-scale current sheets can raise the reconnection rate from small values on the order of the Sweet-Parker rate to high values on the order of the Petscheck rate through triggering an evolution toward fast magnetic reconnection. Meanwhile, the backward coupling between the small- and large-scale magnetic field dynamics has been properly represented through the present high resolution simulation. The undriven turbulent reconnection model established here expresses a solid numerical basis for the previous schematic two-step magnetic reconnection models and a possible explanation of two-stage energy release process of solar explosives.
NASA Astrophysics Data System (ADS)
Otake, Y.; Schafer, S.; Stayman, J. W.; Zbijewski, W.; Kleinszig, G.; Graumann, R.; Khanna, A. J.; Siewerdsen, J. H.
2012-02-01
Localization of target vertebrae is an essential step in minimally invasive spine surgery, with conventional methods relying on "level counting" - i.e., manual counting of vertebrae under fluoroscopy starting from readily identifiable anatomy (e.g., the sacrum). The approach requires an undesirable level of radiation, time, and is prone to counting errors due to the similar appearance of vertebrae in projection images; wrong-level surgery occurs in 1 of every ~3000 cases. This paper proposes a method to automatically localize target vertebrae in x-ray projections using 3D-2D registration between preoperative CT (in which vertebrae are preoperatively labeled) and intraoperative fluoroscopy. The registration uses an intensity-based approach with a gradient-based similarity metric and the CMA-ES algorithm for optimization. Digitally reconstructed radiographs (DRRs) and a robust similarity metric are computed on GPU to accelerate the process. Evaluation in clinical CT data included 5,000 PA and LAT projections randomly perturbed to simulate human variability in setup of mobile intraoperative C-arm. The method demonstrated 100% success for PA view (projection error: 0.42mm) and 99.8% success for LAT view (projection error: 0.37mm). Initial implementation on GPU provided automatic target localization within about 3 sec, with further improvement underway via multi-GPU. The ability to automatically label vertebrae in fluoroscopy promises to streamline surgical workflow, improve patient safety, and reduce wrong-site surgeries, especially in large patients for whom manual methods are time consuming and error prone.