2D-3D rigid registration to compensate for prostate motion during 3D TRUS-guided biopsy
NASA Astrophysics Data System (ADS)
De Silva, Tharindu; Fenster, Aaron; Bax, Jeffrey; Gardi, Lori; Romagnoli, Cesare; Samarabandu, Jagath; Ward, Aaron D.
2012-02-01
Prostate biopsy is the clinical standard for prostate cancer diagnosis. To improve the accuracy of targeting suspicious locations, systems have been developed that can plan and record biopsy locations in a 3D TRUS image acquired at the beginning of the procedure. Some systems are designed for maximum compatibility with existing ultrasound equipment and are thus designed around the use of a conventional 2D TRUS probe, using controlled axial rotation of this probe to acquire a 3D TRUS reference image at the start of the biopsy procedure. Prostate motion during the biopsy procedure causes misalignments between the prostate in the live 2D TRUS images and the pre-acquired 3D TRUS image. We present an image-based rigid registration technique that aligns live 2D TRUS images, acquired immediately prior to biopsy needle insertion, with the pre-acquired 3D TRUS image to compensate for this motion. Our method was validated using 33 manually identified intrinsic fiducials in eight subjects and the target registration error was found to be 1.89 mm. We analysed the suitability of two image similarity metrics (normalized cross correlation and mutual information) for this task by plotting these metrics as a function of varying parameters in the six degree-of-freedom transformation space, with the ground truth plane obtained from registration as the starting point for the parameter exploration. We observed a generally convex behaviour of the similarity metrics. This encourages their use for this registration problem, and could assist in the design of a tool for the detection of misalignment, which could trigger the execution of a non-real-time registration, when needed during the procedure.
An experimental method for eliminating effect of rigid out-of-plane motion on 2D-DIC
NASA Astrophysics Data System (ADS)
Zhiqiang, Wang; Fengzhou, Fang; Bing, Liu; Zhiyong, Wang
2015-10-01
The out-of-plane motion is one of the most important factors that affect the precision of two-dimensional digital image correlation (2D-DIC). In this paper, a novel solution is presented to improve conventional 2D-DIC by eliminating the effect of out-of-plane motion, including translation and rotation. Firstly, an experimental technique using two projected laser strips is proposed to measure the out-of-plane motion of a planar specimen. A theoretical model is then established to predict the pseudostrains caused by out-of-plane motion based on the pin-hole imaging model. Using the measured out-of-plane displacement, the captured deformed images used in 2D-DIC are amended to eliminate the effect of out-of-plane motion by the theoretical model. Finally, two experiments were conducted to validate the effectiveness of the proposed method. Results indicate that application of the proposed method can effectively eliminate the errors caused by out-of-plane motion.
Spoerk, Jakob; Gendrin, Christelle; Weber, Christoph; Figl, Michael; Pawiro, Supriyanto Ardjo; Furtado, Hugo; Fabri, Daniella; Bloch, Christoph; Bergmann, Helmar; Gröller, Eduard; Birkfellner, Wolfgang
2012-01-01
A common problem in image-guided radiation therapy (IGRT) of lung cancer as well as other malignant diseases is the compensation of periodic and aperiodic motion during dose delivery. Modern systems for image-guided radiation oncology allow for the acquisition of cone-beam computed tomography data in the treatment room as well as the acquisition of planar radiographs during the treatment. A mid-term research goal is the compensation of tumor target volume motion by 2D/3D registration. In 2D/3D registration, spatial information on organ location is derived by an iterative comparison of perspective volume renderings, so-called digitally rendered radiographs (DRR) from computed tomography volume data, and planar reference x-rays. Currently, this rendering process is very time consuming, and real-time registration, which should at least provide data on organ position in less than a second, has not come into existence. We present two GPU-based rendering algorithms which generate a DRR of 512 × 512 pixels size from a CT dataset of 53 MB size at a pace of almost 100 Hz. This rendering rate is feasible by applying a number of algorithmic simplifications which range from alternative volume-driven rendering approaches – namely so-called wobbled splatting – to sub-sampling of the DRR-image by means of specialized raycasting techniques. Furthermore, general purpose graphics processing unit (GPGPU) programming paradigms were consequently utilized. Rendering quality and performance as well as the influence on the quality and performance of the overall registration process were measured and analyzed in detail. The results show that both methods are competitive and pave the way for fast motion compensation by rigid and possibly even non-rigid 2D/3D registration and, beyond that, adaptive filtering of motion models in IGRT. PMID:21782399
Flow past 2-D Hemispherical Rigid Canopies
NASA Astrophysics Data System (ADS)
Carnasciali, Maria-Isabel
2013-11-01
The flow past a 2-dimensional rigid hemispherical shape is investigated using PIV. Flow field measurements and images were generated with the use of a Thermoflow® apparatus. Results of this study are compared to prior work (APS DFD 2012 Session E9.00003) which employed CFD to investigate the flow in the near wake of hemispherical parachutes. The various sized gaps/open areas were positioned at distinct locations. The work presented here is part of a larger research project to investigate flow fields in deceleration devices and parachutes. Understanding the pitch-stability of parachutes is essential for accurate design and implementation of these deceleration devices but they present a difficult system to analyze. The flexibility of the parachute fabric results in large variations in the parachute geometry leading to complex fluid-structure interactions. Such flow, combined with flow through gaps and open areas, has been postulated to shed alternating vortices causing pitching/oscillations of the canopy. The results presented provide some insight into which geometric features affect vortex shedding and may enable the redesign of the baseline parachute to minimize instabilities.
Non-rigid target tracking in 2D ultrasound images using hierarchical grid interpolation
NASA Astrophysics Data System (ADS)
Royer, Lucas; Babel, Marie; Krupa, Alexandre
2014-03-01
In this paper, we present a new non-rigid target tracking method within 2D ultrasound (US) image sequence. Due to the poor quality of US images, the motion tracking of a tumor or cyst during needle insertion is considered as an open research issue. Our approach is based on well-known compression algorithm in order to make our method work in real-time which is a necessary condition for many clinical applications. Toward that end, we employed a dedicated hierarchical grid interpolation algorithm (HGI) which can represent a large variety of deformations compared to other motion estimation algorithms such as Overlapped Block Motion Compensation (OBMC), or Block Motion Algorithm (BMA). The sum of squared difference of image intensity is selected as similarity criterion because it provides a good trade-off between computation time and motion estimation quality. Contrary to the others methods proposed in the literature, our approach has the ability to distinguish both rigid and non-rigid motions which are observed in ultrasound image modality. Furthermore, this technique does not take into account any prior knowledge about the target, and limits the user interaction which usually complicates the medical validation process. Finally, a technique aiming at identifying the main phases of a periodic motion (e.g. breathing motion) is introduced. The new approach has been validated from 2D ultrasound images of real human tissues which undergo rigid and non-rigid deformations.
Rotating rigid motion in general relativity
Mason, D.P.; Pooe, C.A.
1987-11-01
Kinematic and dynamic expressions are derived for the Lie derivative of vorticity along a particle world line in a rigid motion. It is found that the evolution of vorticity in a rigid motion is governed by the electric part of the Weyl tensor. Necessary and sufficient kinematic and dynamic conditions are established for a rotating rigid motion to be isometric.
Topology-Preserving Rigid Transformation of 2D Digital Images.
Ngo, Phuc; Passat, Nicolas; Kenmochi, Yukiko; Talbot, Hugues
2014-02-01
We provide conditions under which 2D digital images preserve their topological properties under rigid transformations. We consider the two most common digital topology models, namely dual adjacency and well-composedness. This paper leads to the proposal of optimal preprocessing strategies that ensure the topological invariance of images under arbitrary rigid transformations. These results and methods are proved to be valid for various kinds of images (binary, gray-level, label), thus providing generic and efficient tools, which can be used in particular in the context of image registration and warping. PMID:26270925
Interactive initialization of 2D/3D rigid registration
Gong, Ren Hui; Güler, Özgür; Kürklüoglu, Mustafa; Lovejoy, John; Yaniv, Ziv
2013-12-15
Purpose: Registration is one of the key technical components in an image-guided navigation system. A large number of 2D/3D registration algorithms have been previously proposed, but have not been able to transition into clinical practice. The authors identify the primary reason for the lack of adoption with the prerequisite for a sufficiently accurate initial transformation, mean target registration error of about 10 mm or less. In this paper, the authors present two interactive initialization approaches that provide the desired accuracy for x-ray/MR and x-ray/CT registration in the operating room setting. Methods: The authors have developed two interactive registration methods based on visual alignment of a preoperative image, MR, or CT to intraoperative x-rays. In the first approach, the operator uses a gesture based interface to align a volume rendering of the preoperative image to multiple x-rays. The second approach uses a tracked tool available as part of a navigation system. Preoperatively, a virtual replica of the tool is positioned next to the anatomical structures visible in the volumetric data. Intraoperatively, the physical tool is positioned in a similar manner and subsequently used to align a volume rendering to the x-ray images using an augmented reality (AR) approach. Both methods were assessed using three publicly available reference data sets for 2D/3D registration evaluation. Results: In the authors' experiments, the authors show that for x-ray/MR registration, the gesture based method resulted in a mean target registration error (mTRE) of 9.3 ± 5.0 mm with an average interaction time of 146.3 ± 73.0 s, and the AR-based method had mTREs of 7.2 ± 3.2 mm with interaction times of 44 ± 32 s. For x-ray/CT registration, the gesture based method resulted in a mTRE of 7.4 ± 5.0 mm with an average interaction time of 132.1 ± 66.4 s, and the AR-based method had mTREs of 8.3 ± 5.0 mm with interaction times of 58 ± 52 s. Conclusions: Based on the
Computational Fluid Dynamics Demonstration of Rigid Bodies in Motion
NASA Technical Reports Server (NTRS)
Camarena, Ernesto; Vu, Bruce T.
2011-01-01
The Design Analysis Branch (NE-Ml) at the Kennedy Space Center has not had the ability to accurately couple Rigid Body Dynamics (RBD) and Computational Fluid Dynamics (CFD). OVERFLOW-D is a flow solver that has been developed by NASA to have the capability to analyze and simulate dynamic motions with up to six Degrees of Freedom (6-DOF). Two simulations were prepared over the course of the internship to demonstrate 6DOF motion of rigid bodies under aerodynamic loading. The geometries in the simulations were based on a conceptual Space Launch System (SLS). The first simulation that was prepared and computed was the motion of a Solid Rocket Booster (SRB) as it separates from its core stage. To reduce computational time during the development of the simulation, only half of the physical domain with respect to the symmetry plane was simulated. Then a full solution was prepared and computed. The second simulation was a model of the SLS as it departs from a launch pad under a 20 knot crosswind. This simulation was reduced to Two Dimensions (2D) to reduce both preparation and computation time. By allowing 2-DOF for translations and 1-DOF for rotation, the simulation predicted unrealistic rotation. The simulation was then constrained to only allow translations.
Long waves induced motions to rigid spheroids
NASA Astrophysics Data System (ADS)
Zhou, Hongkun; Hong, Lianjin
2015-05-01
Responses of unconstrained and rigid spheroidal bodies subjected to long sound waves are analyzed by means of approaching hydrodynamic method. It is shown that in the low-frequency approximation the amplitude of translational velocity is completely determined by the density as well as the acoustic added mass which is equal to hydrodynamic one associated with the body. The inconformity of responses to sound waves in virtue of geometric asymmetry is also presented. In addition, rotational movement engendered by acoustic oblique incidence is discussed, and it represents as the modulated angular oscillation similar to the beat-frequency vibration. All these analyses on acoustically induced motions provide a theoretical evidence for developing spheroidal inertial vector receivers.
Rapid determination of RMSDs corresponding to macromolecular rigid body motions.
Popov, Petr; Grudinin, Sergei
2014-05-01
Finding the root mean sum of squared deviations (RMSDs) between two coordinate vectors that correspond to the rigid body motion of a macromolecule is an important problem in structural bioinformatics, computational chemistry, and molecular modeling. Standard algorithms compute the RMSD with time proportional to the number of atoms in the molecule. Here, we present RigidRMSD, a new algorithm that determines a set of RMSDs corresponding to a set of rigid body motions of a macromolecule in constant time with respect to the number of atoms in the molecule. Our algorithm is particularly useful for rigid body modeling applications, such as rigid body docking, and also for high-throughput analysis of rigid body modeling and simulation results. We also introduce a constant-time rotation RMSD as a similarity measure for rigid molecules. A C++ implementation of our algorithm is available at http://nano-d.inrialpes.fr/software/RigidRMSD. PMID:24615729
Flipping and scooping of curved 2D rigid fibers in simple shear: The Jeffery equations
NASA Astrophysics Data System (ADS)
Crowdy, Darren
2016-05-01
The dynamical system governing the motion of a curved rigid two-dimensional circular-arc fiber in simple shear is derived in analytical form. This is achieved by finding the solution for the associated low-Reynolds-number flow around such a fiber using the methods of complex analysis. Solutions of the dynamical system display the "flipping" and "scooping" recently observed in computational studies of three-dimensional fibers using linked rigid rod and bead-shell models [J. Wang et al., "Flipping, scooping, and spinning: Drift of rigid curved nonchiral fibers in simple shear flows," Phys. Fluids 24, 123304 (2012)]. To complete the Jeffery-type equations for a curved fiber in a linear flow field we also derive its evolution equations in an extensional flow. It is expected that the equations derived here also govern the motion of slender, curved, three-dimensional rigid fibers when they evolve purely in the plane of shear or strain.
Rigid Body Motion in Stereo 3D Simulation
ERIC Educational Resources Information Center
Zabunov, Svetoslav
2010-01-01
This paper addresses the difficulties experienced by first-grade students studying rigid body motion at Sofia University. Most quantities describing the rigid body are in relations that the students find hard to visualize and understand. They also lose the notion of cause-result relations between vector quantities, such as the relation between…
Nonrigid 2D registration of fluoroscopic coronary artery image sequence with layered motion
NASA Astrophysics Data System (ADS)
Park, Taewoo; Jung, Hoyup; Yun, Il Dong
2016-03-01
We present a new method for nonrigid registration of coronary artery models with layered motion information. 2D nonrigid registration method is proposed that brings layered motion information into correspondence with fluoroscopic angiograms. The registered model is overlaid on top of interventional angiograms to provide surgical assistance during image-guided chronic total occlusion procedures. The proposed methodology is divided into two parts: layered structures alignments and local nonrigid registration. In the first part, inpainting method is used to estimate a layered rigid transformation that aligns layered motion information. In the second part, a nonrigid registration method is implemented and used to compensate for any local shape discrepancy. Experimental evaluation conducted on a set of 7 fluoroscopic angiograms results in a reduced target registration error, which showed the effectiveness of the proposed method over single layered approach.
On the inertial motions of liquid-filled rigid bodies
NASA Astrophysics Data System (ADS)
Mazzone, Giusy; Galdi, Giovanni; Zunino, Paolo
2013-11-01
We consider a rigid body with a cavity completely filled by a viscous liquid and study the inertial motions of the system liquid-filled rigid body S . The equations governing the motion of this coupled system are given by the Navier-Stokes equations and the equations of the balance of the total angular momentum of S in absence of external forces and torques. Given any initial motion to the coupled system, characterized by an initial relative velocity of the fluid and an initial total angular momentum, we give a complete description of the behavior that the system liquid-filled rigid body will show at large times. From both analytical and numerical viewpoints, we are able to prove a longstanding conjecture stated by Zhukovskii, namely that S will eventually reach a steady state which is a rigid body permanent rotation. In other words, the liquid goes to rest with respect to the rigid body and the coupled system will rotate as a whole rigid body, with a constant angular velocity that is directed along one of the principal axes of inertia of the system.
2-D linear motion system. Innovative technology summary report
1998-11-01
The US Department of Energy's (DOE's) nuclear facility decontamination and decommissioning (D and D) program requires buildings to be decontaminated, decommissioned, and surveyed for radiological contamination in an expeditious and cost-effective manner. Simultaneously, the health and safety of personnel involved in the D and D activities is of primary concern. D and D workers must perform duties high off the ground, requiring the use of manlifts or scaffolding, often, in radiologically or chemically contaminated areas or in areas with limited access. Survey and decontamination instruments that are used are sometimes heavy or awkward to use, particularly when the worker is operating from a manlift or scaffolding. Finding alternative methods of performing such work on manlifts or scaffolding is important. The 2-D Linear Motion System (2-D LMS), also known as the Wall Walker{trademark}, is designed to remotely position tools and instruments on walls for use in such activities as radiation surveys, decontamination, and painting. Traditional (baseline) methods for operating equipment for these tasks require workers to perform duties on elevated platforms, sometimes several meters above the ground surface and near potential sources of contamination. The Wall Walker 2-D LMS significantly improves health and safety conditions by facilitating remote operation of equipment. The Wall Walker 2-D LMS performed well in a demonstration of its precision, accuracy, maneuverability, payload capacity, and ease of use. Thus, this innovative technology is demonstrated to be a viable alternative to standard methods of performing work on large, high walls, especially those that have potential contamination concerns. The Wall Walker was used to perform a final release radiological survey on over 167 m{sup 2} of walls. In this application, surveying using a traditional (baseline) method that employs an aerial lift for manual access was 64% of the total cost of the improved technology
Particle Filters and Occlusion Handling for Rigid 2D-3D Pose Tracking.
Lee, Jehoon; Sandhu, Romeil; Tannenbaum, Allen
2013-08-01
In this paper, we address the problem of 2D-3D pose estimation. Specifically, we propose an approach to jointly track a rigid object in a 2D image sequence and to estimate its pose (position and orientation) in 3D space. We revisit a joint 2D segmentation/3D pose estimation technique, and then extend the framework by incorporating a particle filter to robustly track the object in a challenging environment, and by developing an occlusion detection and handling scheme to continuously track the object in the presence of occlusions. In particular, we focus on partial occlusions that prevent the tracker from extracting an exact region properties of the object, which plays a pivotal role for region-based tracking methods in maintaining the track. To this end, a dynamical choice of how to invoke the objective functional is performed online based on the degree of dependencies between predictions and measurements of the system in accordance with the degree of occlusion and the variation of the object's pose. This scheme provides the robustness to deal with occlusions of an obstacle with different statistical properties from that of the object of interest. Experimental results demonstrate the practical applicability and robustness of the proposed method in several challenging scenarios. PMID:24058277
A rigid motion correction method for helical computed tomography (CT)
NASA Astrophysics Data System (ADS)
Kim, J.-H.; Nuyts, J.; Kyme, A.; Kuncic, Z.; Fulton, R.
2015-03-01
We propose a method to compensate for six degree-of-freedom rigid motion in helical CT of the head. The method is demonstrated in simulations and in helical scans performed on a 16-slice CT scanner. Scans of a Hoffman brain phantom were acquired while an optical motion tracking system recorded the motion of the bed and the phantom. Motion correction was performed by restoring projection consistency using data from the motion tracking system, and reconstructing with an iterative fully 3D algorithm. Motion correction accuracy was evaluated by comparing reconstructed images with a stationary reference scan. We also investigated the effects on accuracy of tracker sampling rate, measurement jitter, interpolation of tracker measurements, and the synchronization of motion data and CT projections. After optimization of these aspects, motion corrected images corresponded remarkably closely to images of the stationary phantom with correlation and similarity coefficients both above 0.9. We performed a simulation study using volunteer head motion and found similarly that our method is capable of compensating effectively for realistic human head movements. To the best of our knowledge, this is the first practical demonstration of generalized rigid motion correction in helical CT. Its clinical value, which we have yet to explore, may be significant. For example it could reduce the necessity for repeat scans and resource-intensive anesthetic and sedation procedures in patient groups prone to motion, such as young children. It is not only applicable to dedicated CT imaging, but also to hybrid PET/CT and SPECT/CT, where it could also ensure an accurate CT image for lesion localization and attenuation correction of the functional image data.
Hurricane Balls: A rigid-body-motion project for undergraduates
NASA Astrophysics Data System (ADS)
Jackson, David P.; Mertens, David; Pearson, Brett J.
2015-11-01
We discuss a project on rigid-body motion that is appropriate for students in an upper-division course in classical mechanics. We analyze the motion of Hurricane Balls, two spheres that are welded (or glued) together so they act as a single object that can be spun like a top. The steady-state motion consists of purely rotational motion about the center of mass, such that only one ball is in contact with the table as it rolls without slipping. We give a qualitative explanation for why one ball rises into the air, and we theoretically analyze the system using multiple approaches. We also perform a high-speed video analysis to obtain experimental data on how the orientation depends on the spin rate, and find agreement within a few percent of the theory.
Rigid facial motion influences featural, but not holistic, face processing.
Xiao, Naiqi G; Quinn, Paul C; Ge, Liezhong; Lee, Kang
2012-03-15
We report three experiments in which we investigated the effect of rigid facial motion on face processing. Specifically, we used the face composite effect to examine whether rigid facial motion influences primarily featural or holistic processing of faces. In Experiments 1-3, participants were first familiarized with dynamic displays in which a target face turned from one side to another; then at test, participants judged whether the top half of a composite face (the top half of the target/foil face aligned or misaligned with the bottom half of a foil face) belonged to the target face. We compared performance in the dynamic condition to various static control conditions in Experiments 1-3, which differed from each other in terms of the display order of the multiple static images or the inter-stimulus interval (ISI) between the images. We found that the size of the face composite effect in the dynamic condition was significantly smaller than that in the static conditions. In other words, the dynamic face display influenced participants to process the target faces in a part-based manner and consequently their recognition of the upper portion of the composite face at test became less interfered with by the aligned lower part of the foil face. The findings from the present experiments provide the strongest evidence to date to suggest that the rigid facial motion mainly influences facial featural, but not holistic, processing. PMID:22342561
The ‘twin paradox’ in relativistic rigid motion
NASA Astrophysics Data System (ADS)
Ben-Ya’acov, Uri
2016-09-01
Relativistic rigid motion suggests a new version for the so-called ‘twin paradox’, comparing the ages of two astronauts on a very long spaceship. Although there is always an instantaneous inertial frame in which the whole spaceship, being rigid, is simultaneously at rest, the twins’ ages, measured as the proper-times along their individual world lines, are different when they are located at remote parts of the spaceship. The age, or proper-time, difference depends on the distance at rest between the astronauts and the rapidity difference between start to end. The relation of the age difference with the relative Doppler shift of light signals transmitted between the astronauts and implications for the possibility to assign a common age (proper-time) to complex, spatially extended, relativistic systems are also discussed.
``Meta''-rigid motions and frames of reference
NASA Astrophysics Data System (ADS)
Bel, Lluís; Llosa, Josep
1995-10-01
We define the “meta”-rigid motions as particular classes of time-like congruences which are solutions of intrinsically defined partial differential equations that generalize Born's conditions. We consider in particular two hierarchies of such congruences. The first one is a geometrically motivated direct generalization of the symmetry concept inherent in Born congruences. The second one is an indirect generalization based on the conditions which guarantee the existence of a particular class of adapted coordinates of space, named quo-harmonic coordinates, whose definition is akin to the definition of harmonic coordinates but which differs from it in an essential point.
Maguinness, Corrina; Newell, Fiona N
2015-04-01
There is growing evidence to suggest that facial motion is an important cue for face recognition. However, it is poorly understood whether motion is integrated with facial form information or whether it provides an independent cue to identity. To provide further insight into this issue, we compared the effect of motion on face perception in two developmental prosopagnosics and age-matched controls. Participants first learned faces presented dynamically (video), or in a sequence of static images, in which rigid (viewpoint) or non-rigid (expression) changes occurred. Immediately following learning, participants were required to match a static face image to the learned face. Test face images varied by viewpoint (Experiment 1) or expression (Experiment 2) and were learned or novel face images. We found similar performance across prosopagnosics and controls in matching facial identity across changes in viewpoint when the learned face was shown moving in a rigid manner. However, non-rigid motion interfered with face matching across changes in expression in both individuals with prosopagnosia compared to the performance of control participants. In contrast, non-rigid motion did not differentially affect the matching of facial expressions across changes in identity for either prosopagnosics (Experiment 3). Our results suggest that whilst the processing of rigid motion information of a face may be preserved in developmental prosopagnosia, non-rigid motion can specifically interfere with the representation of structural face information. Taken together, these results suggest that both form and motion cues are important in face perception and that these cues are likely integrated in the representation of facial identity. PMID:25737056
Percepts of rigid motion within and across apertures
NASA Technical Reports Server (NTRS)
Shiffrar, Maggie; Pavel, M.
1991-01-01
Humans consistently err in their percepts of rotational motion viewed through an aperture. Such errors provide insight into the constraints observers use to interpret retinal images. In the first of two experiments, the subjects consistently perceived the fixed center of rotation for an unmarked line viewed through an aperture as located on the line, regardless of its actual location. Accuracy greatly improved with visible line endings. This finding was extended to explain why a square appears nonrigid when it rotates behind a partial occluder. This illusion is theorized to result from observers misperceiving the center of rotation of the unmarked square sides. In this situation, the subjects seemed unable to apply an object rigidity constraint across apertures. These findings support a conceptualization of the visual system in which consistent local information must be clearly present before prior knowledge can be used to interpret retinal stimulation.
Percepts of rigid motion within and across apertures
NASA Technical Reports Server (NTRS)
Shiffrar, Maggie; Pavel, M.
1991-01-01
Humans consistently err in their percepts of rotational motion viewed through an aperture. Such errors provide insight into the constraints observers use to interpret retinal images. In the first of two experiments, Ss consistently perceived the fixed center of rotation for an unmarked line viewed through an aperture as located on the line, regardless of its actual location. Accuracy greatly improved with visible line endings. This finding was extended to explain why a square appears nonrigid when it rotates behind a partial occluder. This illusion may result from observers misperceiving the center of rotation of the unmarked square sides. In this situation, Ss seemed unable to apply an object rigidity constraint across apertures. These findings support a conceptualization of the visual system in which consistent local information must be clearly present before prior knowledge can be used to interpret retinal stimulation.
NASA Astrophysics Data System (ADS)
Wang, Zenghui; Feng, Philip X.-L.
2016-07-01
Atomic layer crystals are emerging building blocks for enabling new two-dimensional (2D) nanomechanical systems, whose motions can be coupled to other attractive physical properties in such 2D systems. Optical interferometry has been very effective in reading out the infinitesimal motions of these 2D structures and spatially resolving different modes. To quantitatively understand the detection efficiency and its dependence on the device parameters and interferometric conditions, here we present a systematic study of the intrinsic motion responsivity in 2D nanomechanical systems using a Fresnel-law-based model. We find that in monolayer to 14-layer structures, MoS2 offers the highest responsivity among graphene, h-BN, and MoS2 devices and for the three commonly used visible laser wavelengths (633, 532, and 405 nm). We also find that the vacuum gap resulting from the widely used 300 nm-oxide substrate in making 2D devices, fortunately, leads to close-to-optimal responsivity for a wide range of 2D flakes. Our results elucidate and graphically visualize the dependence of motion transduction responsivity upon 2D material type and number of layers, vacuum gap, oxide thickness, and detecting wavelength, thus providing design guidelines for constructing 2D nanomechanical systems with optimal optical motion readout.
Wang, Zenghui; Feng, Philip X-L
2016-01-01
Atomic layer crystals are emerging building blocks for enabling new two-dimensional (2D) nanomechanical systems, whose motions can be coupled to other attractive physical properties in such 2D systems. Optical interferometry has been very effective in reading out the infinitesimal motions of these 2D structures and spatially resolving different modes. To quantitatively understand the detection efficiency and its dependence on the device parameters and interferometric conditions, here we present a systematic study of the intrinsic motion responsivity in 2D nanomechanical systems using a Fresnel-law-based model. We find that in monolayer to 14-layer structures, MoS2 offers the highest responsivity among graphene, h-BN, and MoS2 devices and for the three commonly used visible laser wavelengths (633, 532, and 405 nm). We also find that the vacuum gap resulting from the widely used 300 nm-oxide substrate in making 2D devices, fortunately, leads to close-to-optimal responsivity for a wide range of 2D flakes. Our results elucidate and graphically visualize the dependence of motion transduction responsivity upon 2D material type and number of layers, vacuum gap, oxide thickness, and detecting wavelength, thus providing design guidelines for constructing 2D nanomechanical systems with optimal optical motion readout. PMID:27464908
Wang, Zenghui; Feng, Philip X.-L.
2016-01-01
Atomic layer crystals are emerging building blocks for enabling new two-dimensional (2D) nanomechanical systems, whose motions can be coupled to other attractive physical properties in such 2D systems. Optical interferometry has been very effective in reading out the infinitesimal motions of these 2D structures and spatially resolving different modes. To quantitatively understand the detection efficiency and its dependence on the device parameters and interferometric conditions, here we present a systematic study of the intrinsic motion responsivity in 2D nanomechanical systems using a Fresnel-law-based model. We find that in monolayer to 14-layer structures, MoS2 offers the highest responsivity among graphene, h-BN, and MoS2 devices and for the three commonly used visible laser wavelengths (633, 532, and 405 nm). We also find that the vacuum gap resulting from the widely used 300 nm-oxide substrate in making 2D devices, fortunately, leads to close-to-optimal responsivity for a wide range of 2D flakes. Our results elucidate and graphically visualize the dependence of motion transduction responsivity upon 2D material type and number of layers, vacuum gap, oxide thickness, and detecting wavelength, thus providing design guidelines for constructing 2D nanomechanical systems with optimal optical motion readout. PMID:27464908
Non-Iterative Rigid 2D/3D Point-Set Registration Using Semidefinite Programming
NASA Astrophysics Data System (ADS)
Khoo, Yuehaw; Kapoor, Ankur
2016-07-01
We describe a convex programming framework for pose estimation in 2D/3D point-set registration with unknown point correspondences. We give two mixed-integer nonlinear program (MINP) formulations of the 2D/3D registration problem when there are multiple 2D images, and propose convex relaxations for both of the MINPs to semidefinite programs (SDP) that can be solved efficiently by interior point methods. Our approach to the 2D/3D registration problem is non-iterative in nature as we jointly solve for pose and correspondence. Furthermore, these convex programs can readily incorporate feature descriptors of points to enhance registration results. We prove that the convex programs exactly recover the solution to the original nonconvex 2D/3D registration problem under noiseless condition. We apply these formulations to the registration of 3D models of coronary vessels to their 2D projections obtained from multiple intra-operative fluoroscopic images. For this application, we experimentally corroborate the exact recovery property in the absence of noise and further demonstrate robustness of the convex programs in the presence of noise.
Interactions between rigid-body and flexible-body motions in maneuvering spacecraft
NASA Technical Reports Server (NTRS)
Silverberg, Larry M.; Park, Sungtae
1990-01-01
The present consideration of the significant interactions between rigid-body and flexible-body motions in maneuvering spacecraft proceeds by distinguishing between the two types of motion on the basis of a tracking coordinate system which coincides with the rigid-body component of the motion, as well as by maintaining the motion relative to the tracking coordinate as orthogonal to the rigid-body motion. The elastic motion is excited by the rigid-body motion via Coriolis terms, angular acceleration terms, and centrifugal terms. These interactions are illustrated for spacecraft undergoing bidirectional elastic motions via the dynamics of constantly rotating free-free beams subject to combined bending and longitudinal vibration.
NASA Astrophysics Data System (ADS)
Stahr, Donald W.; Law, Richard D.
2014-11-01
We model the development of shape preferred orientation (SPO) of a large population of two- and three-dimensional (2D and 3D) rigid clasts suspended in a linear viscous matrix deformed by superposed steady and continuously non-steady plane strain flows to investigate the sensitivity of clasts to changing boundary conditions during a single or superposed deformation events. Resultant clast SPOs are compared to one developed by an identical initial population that experienced a steady flow history of constant kinematic vorticity and reached an identical finite strain state, allowing examination of SPO sensitivity to deformation path. Rotation paths of individual triaxial inclusions are complex, even for steady plane strain flow histories. It has been suggested that the 3D nature of the system renders predictions based on 2D models inadequate for applied clast-based kinematic vorticity gauges. We demonstrate that for a large population of clasts, simplification to a 2D model does provide a good approximation to the SPO predicted by full 3D analysis for steady and non-steady plane strain deformation paths. Predictions of shape fabric development from 2D models are not only qualitatively similar to the more complex 3D analysis, but they display the same limitations of techniques based on clast SPO commonly used as a quantitative kinematic vorticity gauge. Our model results from steady, superposed, and non-steady flow histories with a significant pure shearing component at a wide range of finite strain resemble predictions for an identical initial population that experienced a single steady simple shearing deformation. We conclude that individual 2D and 3D clasts respond instantaneously to changes in boundary conditions, however, in aggregate, the SPO of a population of rigid inclusions does not reflect the late-stage kinematics of deformation, nor is it an indicator of the unique 'mean' kinematic vorticity experienced by a deformed rock volume.
Evaluation of low-dose limits in 3D-2D rigid registration for surgical guidance
NASA Astrophysics Data System (ADS)
Uneri, A.; Wang, A. S.; Otake, Y.; Kleinszig, G.; Vogt, S.; Khanna, A. J.; Gallia, G. L.; Gokaslan, Z. L.; Siewerdsen, J. H.
2014-09-01
An algorithm for intensity-based 3D-2D registration of CT and C-arm fluoroscopy is evaluated for use in surgical guidance, specifically considering the low-dose limits of the fluoroscopic x-ray projections. The registration method is based on a framework using the covariance matrix adaptation evolution strategy (CMA-ES) to identify the 3D patient pose that maximizes the gradient information similarity metric. Registration performance was evaluated in an anthropomorphic head phantom emulating intracranial neurosurgery, using target registration error (TRE) to characterize accuracy and robustness in terms of 95% confidence upper bound in comparison to that of an infrared surgical tracking system. Three clinical scenarios were considered: (1) single-view image + guidance, wherein a single x-ray projection is used for visualization and 3D-2D guidance; (2) dual-view image + guidance, wherein one projection is acquired for visualization, combined with a second (lower-dose) projection acquired at a different C-arm angle for 3D-2D guidance; and (3) dual-view guidance, wherein both projections are acquired at low dose for the purpose of 3D-2D guidance alone (not visualization). In each case, registration accuracy was evaluated as a function of the entrance surface dose associated with the projection view(s). Results indicate that images acquired at a dose as low as 4 μGy (approximately one-tenth the dose of a typical fluoroscopic frame) were sufficient to provide TRE comparable or superior to that of conventional surgical tracking, allowing 3D-2D guidance at a level of dose that is at most 10% greater than conventional fluoroscopy (scenario #2) and potentially reducing the dose to approximately 20% of the level in a conventional fluoroscopically guided procedure (scenario #3).
Elastic image registration via rigid object motion induced deformation
NASA Astrophysics Data System (ADS)
Zheng, Xiaofen; Udupa, Jayaram K.; Hirsch, Bruce E.
2011-03-01
In this paper, we estimate the deformations induced on soft tissues by the rigid independent movements of hard objects and create an admixture of rigid and elastic adaptive image registration transformations. By automatically segmenting and independently estimating the movement of rigid objects in 3D images, we can maintain rigidity in bones and hard tissues while appropriately deforming soft tissues. We tested our algorithms on 20 pairs of 3D MRI datasets pertaining to a kinematic study of the flexibility of the ankle complex of normal feet as well as ankles affected by abnormalities in foot architecture and ligament injuries. The results show that elastic image registration via rigid object-induced deformation outperforms purely rigid and purely nonrigid approaches.
Personalized x-ray reconstruction of the proximal femur via a non-rigid 2D-3D registration
NASA Astrophysics Data System (ADS)
Yu, Weimin; Zysset, Philippe; Zheng, Guoyan
2015-03-01
In this paper we present a new approach for a personalized X-ray reconstruction of the proximal femur via a non-rigid registration of a 3D volumetric template to 2D calibrated C-arm images. The 2D-3D registration is done with a hierarchical two-stage strategy: the global scaled rigid registration stage followed by a regularized deformable b-spline registration stage. In both stages, a set of control points with uniform spacing are placed over the domain of the 3D volumetric template and the registrations are driven by computing updated positions of these control points, which then allows to accurately register the 3D volumetric template to the reference space of the C-arm images. Comprehensive experiments on simulated images, on images of cadaveric femurs and on clinical datasets are designed and conducted to evaluate the performance of the proposed approach. Quantitative and qualitative evaluation results are given, which demonstrate the efficacy of the present approach.
An incompressible smoothed particle hydrodynamics method for the motion of rigid bodies in fluids
NASA Astrophysics Data System (ADS)
Tofighi, N.; Ozbulut, M.; Rahmat, A.; Feng, J. J.; Yildiz, M.
2015-09-01
A two-dimensional incompressible smoothed particle hydrodynamics scheme is presented for simulation of rigid bodies moving through Newtonian fluids. The scheme relies on combined usage of the rigidity constraints and the viscous penalty method to simulate rigid body motion. Different viscosity ratios and interpolation schemes are tested by simulating a rigid disc descending in quiescent medium. A viscosity ratio of 100 coupled with weighted harmonic averaging scheme has been found to provide satisfactory results. The performance of the resulting scheme is systematically tested for cases with linear motion, rotational motion and their combination. The test cases include sedimentation of a single and a pair of circular discs, sedimentation of an elliptic disc and migration and rotation of a circular disc in linear shear flow. Comparison with previous results at various Reynolds numbers indicates that the proposed method captures the motion of rigid bodies driven by flow or external body forces accurately.
Learning Spatially-Smooth Mappings in Non-Rigid Structure from Motion.
Hamsici, Onur C; Gotardo, Paulo F U; Martinez, Aleix M
2012-01-01
Non-rigid structure from motion (NRSFM) is a classical underconstrained problem in computer vision. A common approach to make NRSFM more tractable is to constrain 3D shape deformation to be smooth over time. This constraint has been used to compress the deformation model and reduce the number of unknowns that are estimated. However, temporal smoothness cannot be enforced when the data lacks temporal ordering and its benefits are less evident when objects undergo abrupt deformations. This paper proposes a new NRSFM method that addresses these problems by considering deformations as spatial variations in shape space and then enforcing spatial, rather than temporal, smoothness. This is done by modeling each 3D shape coefficient as a function of its input 2D shape. This mapping is learned in the feature space of a rotation invariant kernel, where spatial smoothness is intrinsically defined by the mapping function. As a result, our model represents shape variations compactly using custom-built coefficient bases learned from the input data, rather than a pre-specified set such as the Discrete Cosine Transform. The resulting kernel-based mapping is a by-product of the NRSFM solution and leads to another fundamental advantage of our approach: for a newly observed 2D shape, its 3D shape is recovered by simply evaluating the learned function. PMID:23946937
NASA Astrophysics Data System (ADS)
Olmez, O.; Ozbulut, M.; Yildiz, M.; Goren, O.
2016-06-01
The present study investigates the vortical and nonlinear effects in the roll motion of a 2-D body with square cross-sections by using Smoothed Particle Hydrodynamics (SPH). A 2-D rigid body with square cross-section is taken into account for the benchmark study and subjected to the oscillatory roll motion with a given angular frequency. The governing equations are continuity equation and Euler's equation with artificial viscosity term. Weakly Compressible SPH (WCSPH) scheme is employed for the discretization of the governing equations. Velocities of the fluid particles are updated by means of XSPH+Artificial Particle Displacement (VXSPH+APD) algorithm. In this method only the free surface fluid particles are subjected to VXSPH algorithm while the APD algorithm is employed for the fully populated flow regions. The hybrid usage of numerical treatment keeps free surface particles together by creating an artificial surface tension on the free surface. VXSPH+APD is a proven numerical treatment to provide the most accurate results for this type of free surface flows (Ozbulut et al. 2014). The results of the present study are compared with those of the experimental studies as well as with those of the numerical methods obtained from the current literature.
Lorenz, Kevin S.; Salama, Paul; Dunn, Kenneth W.; Delp, Edward J.
2013-01-01
Digital image analysis is a fundamental component of quantitative microscopy. However, intravital microscopy presents many challenges for digital image analysis. In general, microscopy volumes are inherently anisotropic, suffer from decreasing contrast with tissue depth, lack object edge detail, and characteristically have low signal levels. Intravital microscopy introduces the additional problem of motion artifacts, resulting from respiratory motion and heartbeat from specimens imaged in vivo. This paper describes an image registration technique for use with sequences of intravital microscopy images collected in time-series or in 3D volumes. Our registration method involves both rigid and non-rigid components. The rigid registration component corrects global image translations, while the non-rigid component manipulates a uniform grid of control points defined by B-splines. Each control point is optimized by minimizing a cost function consisting of two parts: a term to define image similarity, and a term to ensure deformation grid smoothness. Experimental results indicate that this approach is promising based on the analysis of several image volumes collected from the kidney, lung, and salivary gland of living rodents. PMID:22092443
Knowledge-In-Action: An Example with Rigid Body Motion
ERIC Educational Resources Information Center
Da Costa, Sayonara Salvador Cabral; Moreira, Marco Antonio
2005-01-01
This paper reports the analysis of the resolution of a paper-and-pencil problem, by eight undergraduate students majoring in engineering (six) and physics (two) at the Pontifcia Universidade Catlica do Rio Grande do Sul, in Porto Alegre, Brazil. The problem concerns kinetics of a rigid body, and the analysis was done in the light of Johnson-Lairds…
Lorentz Contraction, Bell's Spaceships and Rigid Body Motion in Special Relativity
ERIC Educational Resources Information Center
Franklin, Jerrold
2010-01-01
The meaning of Lorentz contraction in special relativity and its connection with Bell's spaceships parable is discussed. The motion of Bell's spaceships is then compared with the accelerated motion of a rigid body. We have tried to write this in a simple form that could be used to correct students' misconceptions due to conflicting earlier…
Edge preserving motion estimation with occlusions correction for assisted 2D to 3D conversion
NASA Astrophysics Data System (ADS)
Pohl, Petr; Sirotenko, Michael; Tolstaya, Ekaterina; Bucha, Victor
2014-02-01
In this article we propose high quality motion estimation based on variational optical flow formulation with non-local regularization term. To improve motion in occlusion areas we introduce occlusion motion inpainting based on 3-frame motion clustering. Variational formulation of optical flow proved itself to be very successful, however a global optimization of cost function can be time consuming. To achieve acceptable computation times we adapted the algorithm that optimizes convex function in coarse-to-fine pyramid strategy and is suitable for modern GPU hardware implementation. We also introduced two simplifications of cost function that significantly decrease computation time with acceptable decrease of quality. For motion clustering based motion inpaitning in occlusion areas we introduce effective method of occlusion aware joint 3-frame motion clustering using RANSAC algorithm. Occlusion areas are inpainted by motion model taken from cluster that shows consistency in opposite direction. We tested our algorithm on Middlebury optical flow benchmark, where we scored around 20th position, but being one of the fastest method near the top. We also successfully used this algorithm in semi-automatic 2D to 3D conversion tool for spatio-temporal background inpainting, automatic adaptive key frame detection and key points tracking.
Sequential Non-Rigid Structure from Motion Using Physical Priors.
Agudo, Antonio; Moreno-Noguer, Francesc; Calvo, Begona; Montiel, Jose M Martinez
2016-05-01
We propose a new approach to simultaneously recover camera pose and 3D shape of non-rigid and potentially extensible surfaces from a monocular image sequence. For this purpose, we make use of the Extended Kalman Filter based Simultaneous Localization And Mapping (EKF-SLAM) formulation, a Bayesian optimization framework traditionally used in mobile robotics for estimating camera pose and reconstructing rigid scenarios. In order to extend the problem to a deformable domain we represent the object's surface mechanics by means of Navier's equations, which are solved using a Finite Element Method (FEM). With these main ingredients, we can further model the material's stretching, allowing us to go a step further than most of current techniques, typically constrained to surfaces undergoing isometric deformations. We extensively validate our approach in both real and synthetic experiments, and demonstrate its advantages with respect to competing methods. More specifically, we show that besides simultaneously retrieving camera pose and non-rigid shape, our approach is adequate for both isometric and extensible surfaces, does not require neither batch processing all the frames nor tracking points over the whole sequence and runs at several frames per second. PMID:27046840
Könik, Arda; Connolly, Caitlin M; Johnson, Karen L; Dasari, Paul; Segars, Paul W; Pretorius, P H; Lindsay, Clifford; Dey, Joyoni; King, Michael A
2014-01-01
The development of methods for correcting patient motion in emission tomography has been receiving increased attention. Often performance of these methods is evaluated through simulations using digital anthropomorphic phantoms, such as the commonly used XCAT phantom, which models both respiratory and cardiac motion based on human studies. However, non-rigid body motion, which is frequently seen in clinical studies, is not present in the standard XCAT phantom. In addition, respiratory motion in the standard phantom is limited to a single generic trend. In this work, to obtain more realistic representation of motion, we developed a series of individual-specific XCAT phantoms modeling non-rigid respiratory and non-rigid body motions derived from the MRI acquisitions of volunteers. Acquisitions were performed in the sagittal orientation using the Navigator methodology. Baseline (no motion) acquisitions at end-expiration were obtained at the beginning of each imaging session for each volunteer. For the body motion studies, MRI was again acquired only at end-expiration for five body motion poses (shoulder stretch, shoulder twist, lateral bend, side roll, and axial slide). For the respiratory motion studies, MRI was acquired during free/regular breathing. The MR slices were then retrospectively sorted into 14 amplitude-binned respiratory states, end-expiration, end-inspiration, six intermediary states during inspiration, and six during expiration using the recorded Navigator signal. XCAT phantoms were then generated based on these MRI data by interactive alignment of the organ contours of the XCAT with the MRI slices using a GUI. Thus far we have created 5 body motion and 5 respiratory motion XCAT phantoms from MRI acquisitions of 6 healthy volunteers (3 males and 3 females). Non-rigid motion exhibited by the volunteers was reflected in both respiratory and body motion phantoms with a varying extent and character for each individual. In addition to these phantoms, we
NASA Astrophysics Data System (ADS)
Könik, Arda; Connolly, Caitlin M.; Johnson, Karen L.; Dasari, Paul; Segars, Paul W.; Pretorius, P. H.; Lindsay, Clifford; Dey, Joyoni; King, Michael A.
2014-07-01
The development of methods for correcting patient motion in emission tomography has been receiving increased attention. Often the performance of these methods is evaluated through simulations using digital anthropomorphic phantoms, such as the commonly used extended cardiac torso (XCAT) phantom, which models both respiratory and cardiac motion based on human studies. However, non-rigid body motion, which is frequently seen in clinical studies, is not present in the standard XCAT phantom. In addition, respiratory motion in the standard phantom is limited to a single generic trend. In this work, to obtain a more realistic representation of motion, we developed a series of individual-specific XCAT phantoms, modeling non-rigid respiratory and non-rigid body motions derived from the magnetic resonance imaging (MRI) acquisitions of volunteers. Acquisitions were performed in the sagittal orientation using the Navigator methodology. Baseline (no motion) acquisitions at end-expiration were obtained at the beginning of each imaging session for each volunteer. For the body motion studies, MRI was again acquired only at end-expiration for five body motion poses (shoulder stretch, shoulder twist, lateral bend, side roll, and axial slide). For the respiratory motion studies, an MRI was acquired during free/regular breathing. The magnetic resonance slices were then retrospectively sorted into 14 amplitude-binned respiratory states, end-expiration, end-inspiration, six intermediary states during inspiration, and six during expiration using the recorded Navigator signal. XCAT phantoms were then generated based on these MRI data by interactive alignment of the organ contours of the XCAT with the MRI slices using a graphical user interface. Thus far we have created five body motion and five respiratory motion XCAT phantoms from the MRI acquisitions of six healthy volunteers (three males and three females). Non-rigid motion exhibited by the volunteers was reflected in both respiratory
Rigid 2D/3D registration of intraoperative digital x-ray images and preoperative CT and MR images
NASA Astrophysics Data System (ADS)
Tomazevic, Dejan; Likar, Bostjan; Pernus, Franjo
2002-05-01
This paper describes a novel approach to register 3D computed tomography (CT) or magnetic resonance (MR) images to a set of 2D X-ray images. Such a registration may be a valuable tool for intraoperative determination of the precise position and orientation of some anatomy of interest, defined in preoperative images. The registration is based solely on the information present in 2D and 3D images. It does not require fiducial markers, X-ray image segmentation, or construction of digitally reconstructed radiographs. The originality of the approach is in using normals to bone surfaces, preoperatively defined in 3D MR or CT data, and gradients of intraoperative X-ray images, which are back-projected towards the X-ray source. The registration is then concerned with finding that rigid transformation of a CT or MR volume, which provides the best match between surface normals and back projected gradients, considering their amplitudes and orientations. The method is tested on a lumbar spine phantom. Gold standard registration is obtained by fidicual markers attached to the phantom. Volumes of interest, containing single vertebrae, are registered to different pairs of X-ray images from different starting positions, chosen randomly and uniformly around the gold standard position. Target registration errors and rotation errors are in order of 0.3 mm and 0.35 degrees for the CT to X-ray registration and 1.3 mm and 1.5 degrees for MR to X-ray registration. The registration is shown to be fast and accurate.
Rupture dynamics and ground motions from earthquakes in 2-D heterogeneous media
NASA Astrophysics Data System (ADS)
Bydlon, Samuel A.; Dunham, Eric M.
2015-03-01
We perform 2-D simulations of earthquakes on rough faults in media with random heterogeneities (with von Karman distribution) to study the effects of geometric and material heterogeneity on the rupture process and resulting high-frequency ground motions in the near-fault region (out to ˜20 km). Variations in slip and rupture velocity can arise from material heterogeneity alone but are dominantly controlled by fault roughness. Scattering effects become appreciable beyond ˜3 km from the fault. Near-fault scattering extends the duration of incoherent, high-frequency ground motions and, at least in our 2-D simulations, elevates root-mean-square accelerations (i.e., Arias intensity) with negligible reduction in peak velocities. We also demonstrate that near-fault scattering typically occurs in the power law tail of the power spectral density function, quantified by the Hurst exponent and another parameter combining standard deviation and correlation length.
Estimation of motion parameters for a rigid body from its orthogonal projection
NASA Technical Reports Server (NTRS)
Ganguly, S.; Ghosh, B.; Tarn, T. J.; Bejczy, A. K.
1989-01-01
An estimate is presented of the motion parameters, namely, linear and angular velocities of a rigid body rotating and translating in three-dimensional-space. It is assumed that the velocities are constant and that only the orthogonal projection of the motion is observable. In particular, if (x, y, z) is the Cartesian coordinate, it is assumed that the projection of the motion on the x-y plane is observed and the information along the z coordinate is lost.
NASA Technical Reports Server (NTRS)
Fitzpatrick, P. M.; Harmon, G. R.; Cochran, J. E.; Shaw, W. A.
1974-01-01
Some methods of approaching a solution to the Hamilton-Jacobi equation are outlined and examples are given to illustrate particular methods. These methods may be used for cases where the Hamilton-Jacobi equation is not separable and have been particularly useful in solving the rigid body motion of an earth satellite subjected to gravity torques. These general applications may also have usefulness in studying the motion of satellites with aerodynamic torque and in studying space vehicle motion where thrusting is involved.
Hurricane Balls: A rigid-body-motion student project
NASA Astrophysics Data System (ADS)
Jackson, David; Mertens, David; Pearson, Brett
Hurricane Balls is a spinning-top toy that consists of two metal spheres that are welded (or glued) together. The motion of Hurricane Balls provides a beautiful example of rotational motion in which the angular velocity and angular momentum point in different directions. Because the motion is both captivating to students and extremely reproducible, this system is an ideal example to include in a classical mechanics course. Moreover, the excellent agreement between theory and experiment makes a detailed analysis of Hurricane Balls a perfect topic for an independent student project. This talk will give an overview of the system and will provide some tips on how to make such a project a successful student experience.
Rigid Structure from Motion from a Blind Source Separation Perspective
Fortuna, Jeff
2013-01-01
We present an information theoretic approach to define the problem of structure from motion (SfM) as a blind source separation one. Given that for almost all practical joint densities of shape points, the marginal densities are non-Gaussian, we show how higher-order statistics can be used to provide improvements in shape estimates over the methods of factorization via Singular Value Decomposition (SVD), bundle adjustment and Bayesian approaches. Previous techniques have either explicitly or implicitly used only second-order statistics in models of shape or noise. A further advantage of viewing SfM as a blind source problem is that it easily allows for the inclusion of noise and shape models, resulting in Maximum Likelihood (ML) or Maximum a Posteriori (MAP) shape and motion estimates. A key result is that the blind source separation approach has the ability to recover the motion and shape matrices without the need to explicitly know the motion or shape pdf. We demonstrate that it suffices to know whether the pdf is sub-or super-Gaussian (i.e., semi-parametric estimation) and derive a simple formulation to determine this from the data. We provide extensive experimental results on synthetic and real tracked points in order to quantify the improvement obtained from this technique. PMID:23682206
NASA Astrophysics Data System (ADS)
Mechler, S.; Pershan, P. S.; Yahel, E.; Stoltz, S. E.; Shpyrko, O. G.; Lin, B.; Meron, M.; Sellner, S.
2010-10-01
The structural and mechanical properties of 2D crystalline surface phases that form at the surface of liquid eutectic Au82Si18 are studied using synchrotron x-ray scattering over a large temperature range. In the vicinity of the eutectic temperature the surface consists of a 2D atomic bilayer crystalline phase that transforms into a 2D monolayer crystalline phase during heating. The latter phase eventually melts into a liquidlike surface on further heating. We demonstrate that the short wavelength capillary wave fluctuations are suppressed due to the bending rigidity of 2D crystalline phases. The corresponding reduction in the Debye-Waller factor allows for measured reflectivity to be explained in terms of an electron density profile that is consistent with the 2D surface crystals.
The use of 2D ultrasound elastography for measuring tendon motion and strain.
Chernak Slane, Laura; Thelen, Darryl G
2014-02-01
The goal of the current study was to investigate the fidelity of a 2D ultrasound elastography method for the measurement of tendon motion and strain. Ultrasound phantoms and ex vivo porcine flexor tendons were cyclically stretched to 4% strain while cine ultrasound radiofrequency (RF) data and video data were simultaneously collected. 2D ultrasound elastography was used to estimate tissue motion and strain from RF data, and surface tissue motion and strain were separately estimated using digital image correlation (DIC). There were strong correlations (R(2)>0.97) between DIC and RF measurements of phantom displacement and strain, and good agreement in estimates of peak phantom strain (DIC: 3.5±0.2%; RF: 3.7±0.1%). For tendon, elastographic estimates of displacement profiles also correlated well with DIC measurements (R(2)>0.92), and exhibited similar estimated peak tendon strain (DIC: 2.6±1.4%; RF: 2.2±1.3%). Elastographic tracking with B-Mode images tended to under-predict peak strain for both the phantom and tendon. This study demonstrates the capacity to use quantitative elastographic techniques to measure tendon displacement and strain within an ultrasound image window. The approach may be extendible to in vivo use on humans, which would allow for the non-invasive analysis of tendon deformation in both normal and pathological states. PMID:24388164
The Use of 2D Ultrasound Elastography for Measuring Tendon Motion and Strain
Slane, Laura Chernak; Thelen, Darryl G.
2014-01-01
The goal of the current study was to investigate the fidelity of a 2D ultrasound elastography method for the measurement of tendon motion and strain. Ultrasound phantoms and ex vivo porcine flexor tendons were cyclically stretched to 4% strain while cine ultrasound radiofrequency (RF) data and video data were simultaneously collected. 2D ultrasound elastography was used to estimate tissue motion and strain from RF data, and surface tissue motion and strain were separately estimated using digital image correlation (DIC). There were strong correlations (R2 > 0.97) between DIC and RF measurements of phantom displacement and strain, and good agreement in estimates of peak phantom strain (DIC: 3.5 ± 0.2%; RF: 3.7 ± 0.1%). For tendon, elastographic estimates of displacement profiles also correlated well with DIC measurements (R2 > 0.92), and exhibited similar estimated peak tendon strain (DIC: 2.6 ± 1.4%; RF: 2.2 ± 1.3%). Elastographic tracking with B-Mode images tended to under-predict peak strain for both the phantom and tendon. This study demonstrates the capacity to use quantitative elastographic techniques to measure tendon displacement and strain within an ultrasound image window. The approach may be extendible to in vivo use on humans, which would allow for the non-invasive analysis of tendon deformation in both normal and pathological states. PMID:24388164
Coronary arteries motion modeling on 2D x-ray images
NASA Astrophysics Data System (ADS)
Gao, Yang; Sundar, Hari
2012-02-01
During interventional procedures, 3D imaging modalities like CT and MRI are not commonly used due to interference with the surgery and radiation exposure concerns. Therefore, real-time information is usually limited and building models of cardiac motion are difficult. In such case, vessel motion modeling based on 2-D angiography images become indispensable. Due to issues with existing vessel segmentation algorithms and the lack of contrast in occluded vessels, manual segmentation of certain branches is usually necessary. In addition, such occluded branches are the most important vessels during coronary interventions and obtaining motion models for these can greatly help in reducing the procedure time and radiation exposure. Segmenting different cardiac phases independently does not guarantee temporal consistency and is not efficient for occluded branches required manual segmentation. In this paper, we propose a coronary motion modeling system which extracts the coronary tree for every cardiac phase, maintaining the segmentation by tracking the coronary tree during the cardiac cycle. It is able to map every frame to the specific cardiac phase, thereby inferring the shape information of the coronary arteries using the model corresponding to its phase. Our experiments show that our motion modeling system can achieve promising results with real-time performance.
Emergence of electromotive force in precession-less rigid motion of deformed domain wall
NASA Astrophysics Data System (ADS)
Farajollahpour, Tohid; Darmiani, Narges; Phirouznia, Arash
2016-08-01
Recently it has been recognized that the electromotive force (emf) can be induced just by the spin precession where the generation of the electromotive force has been considered as a real-space topological pumping effect. It has been shown that the amount of the electromotive force is independent of the functionality of the localized moments. It was also demonstrated that the rigid domain wall (DW) motion cannot generate electromotive force in the system. Based on real-space topological pumping approach in the current study we show that the electromotive force can be induced by rigid motion of a deformed DW. We also demonstrate that the generated electromotive force strongly depends on the DW bulging. Meanwhile results show that the DW bulging leads to generation of the electromotive force both along the axis of the DW motion and normal to the direction of motion.
On the motion of a heavy rigid body in an ideal fluid with circulation.
Borisov, Alexey V; Mamaev, Ivan S
2006-03-01
We consider Chaplygin's equations [Izd. Akad. Nauk SSSR 3, 3 (1933)] describing the planar motion of a rigid body in an unbounded volume of an ideal fluid while circulation around the body is not zero. Hamiltonian structures and new integrable cases are revealed; certain remarkable partial solutions are found and their stability is examined. The nonintegrability of the system describing the motion of a body in the field of gravity is proved and the chaotic behavior of the system is illustrated. PMID:16599749
Geometrical connection between catacaustics and kinematics of planar motion of a rigid solid.
Bellver-Cebreros, Consuelo; Rodríguez-Danta, Marcelo
2016-09-01
Unnoticed and hidden optomechanical analogies between kinematics of planar motion of a rigid solid and catacaustics generated by mirror reflection on smooth profiles in geometrical optics are discussed. A concise and self-consistent theory is developed, which intends to explain and clarify many partial aspects covered by the literature. PMID:27607500
ERIC Educational Resources Information Center
Stanford Univ., CA. School Mathematics Study Group.
The teacher's guide for the eleventh unit in this SMSG series covers the chapter on rigid motions and vectors and the chapter on computers and programs. The overall purpose for each of the chapters is described, the prerequisite knowledge needed by students is specified, the mathematical development of each chapter is detailed, behavioral…
Pan, Xiaochang; Liu, Ke; Shao, Jinghua; Gao, Jing; Huang, Lingyun; Bai, Jing; Luo, Jianwen
2015-11-01
Tissue motion estimation is widely used in many ultrasound techniques. Rigid-model-based and nonrigid-modelbased methods are two main groups of space-domain methods of tissue motion estimation. The affine model is one of the commonly used nonrigid models. The performances of the rigid model and affine model have not been compared on ultrasound RF signals, which have been demonstrated to obtain higher accuracy, precision, and resolution in motion estimation compared with B-mode images. In this study, three methods, i.e., the normalized cross-correlation method with rigid model (NCC), the optical flow method with rigid model (OFRM), and the optical flow method with affine model (OFAM), are compared using ultrasound RF signals, rather than the B-mode images used in previous studies. Simulations, phantom, and in vivo experiments are conducted to make the comparison. In the simulations, the root-mean-square errors (RMSEs) of axial and lateral displacements and strains are used to assess the accuracy of motion estimation, and the elastographic signal-tonoise ratio (SNRe) and contrast-to-noise ratio (CNRe) are used to evaluate the quality of axial strain images. In the phantom experiments, the registration error between the pre- and postdeformation RF signals, as well as the SNRe and CNRe of axial strain images, are utilized as the evaluation criteria. In the in vivo experiments, the registration error is used to evaluate the estimation performance. The results show that the affinemodel- based method (i.e., OFAM) obtains the lowest RMSE or registration error and the highest SNRe and CNRe among all the methods. The affine model is demonstrated to be superior to the rigid model in motion estimation based on RF signals. PMID:26559623
An Interface for Specifying Rigid-Body Motions for CFD Applications
NASA Technical Reports Server (NTRS)
Murman, Scott M.; Chan, William; Aftosmis, Michael; Meakin, Robert L.; Kwak, Dochan (Technical Monitor)
2003-01-01
An interface for specifying rigid-body motions for CFD applications is presented. This interface provides a means of describing a component hierarchy in a geometric configuration, as well as the motion (prescribed or six-degree-of-freedom) associated with any component. The interface consists of a general set of datatypes, along with rules for their interaction, and is designed to be flexible in order to evolve as future needs dictate. The specification is currently implemented with an XML file format which is portable across platforms and applications. The motion specification is capable of describing general rigid body motions, and eliminates the need to write and compile new code within the application software for each dynamic configuration, allowing client software to automate dynamic simulations. The interface is integrated with a GUI tool which allows rigid body motions to be prescribed and verified interactively, promoting access to non-expert users. Illustrative examples, as well as the raw XML source of the file specifications, are included.
Nonlinear state-space modeling of human motion using 2-D marker observations.
Vartiainen, Paavo; Bragge, Timo; Arokoski, Jari P; Karjalainen, Pasi A
2014-07-01
A novel method for the estimation of human kinematics, based on state-space modeling, is proposed. The state consists of the positions, orientations, velocities, and accelerations of an articulated model. Estimation is performed using the unscented Kalman filter (UKF) algorithm with a fixed-interval smoother. Impulsive acceleration at floor contact of the foot is estimated by implementing a contact constraint in the UKF evolution model. The constraint inserts an acceleration impulse into the model state. The estimation method was applied to marker-based motion analysis in a motion laboratory. Validation measurements were performed with a rigid test device and with human gait. A triaxial accelerometer was used to evaluate acceleration estimates. Comparison between the proposed method and the extended Kalman smoother showed a clear difference in the quality of estimates during impulsive accelerations. The proposed approach enables estimation of human kinematics during both continuous and transient accelerations. The approach provides a novel way of estimating acceleration at foot initial contact, and thus enables more accurate evaluation of loading from the beginning of the floor contact. PMID:24760898
Coupled motion of rigid bodies about their center of mass. [Shuttle/payload system
NASA Technical Reports Server (NTRS)
Jezewski, D. J.; Donaldson, J. D.
1979-01-01
Nontrivial analytical solutions for the coupled motion of two rigid bodies about their center of mass are obtained on the assumptions that the rigid bodies are coupled by a massless rigid boom and that no external forces are acting on the system. Both relative rotational and translational motions of the two bodies are considered. General equations of motion are derived by regarding the two bodies as consisting of two distinct systems of particles and by applying the principle of conservation of angular momentum. It is shown that a basic nontrivial solution can be obtained for the translational problem if an assumption is made concerning the relative orientation of one principal axis of inertia of each body and that fundamental nontrivial solutions are readily obtained for the rotational problem if an additional assumption is made with respect to the symmetry of one body. Certain stability criteria are found for some of these motions by defining regions of constraint for the relative translational and rotational elements.
Quantifying Rigid and Nonrigid Motion of Liver Tumors During Stereotactic Body Radiation Therapy
Xu, Qianyi; Hanna, George; Grimm, Jimm; Kubicek, Gregory; Pahlajani, Niraj; Asbell, Sucha; Fan, Jiajin; Chen, Yan; LaCouture, Tamara
2014-09-01
Purpose: To quantify rigid and nonrigid motion of liver tumors using reconstructed 3-dimensional (3D) fiducials from stereo imaging during CyberKnife-based stereotactic body radiation therapy (SBRT). Methods and Materials: Twenty-three liver patients treated with 3 fractions of SBRT were used in this study. After 2 orthogonal kilovoltage images were taken during treatment, the 3D locations of the fiducials were generated by the CyberKnife system and validated using geometric derivations. A total of 4824 pairs of kilovoltage images from start to end of treatment were analyzed. For rigid motion, the rotational angles and translational shifts were reported by aligning 3D fiducial groups from different image pairs, using least-squares fitting. For nonrigid motion, we quantified interfractional tumor volume variations by using the proportional volume derived from the fiducials, which correlates to the sum of interfiducial distances. The individual fiducial displacements were also reported (1) after rigid corrections and (2) without angle corrections. Results: The proportional volume derived by the fiducials demonstrated a volume-increasing trend in the second (101.9% ± 3.6%) and third (101.0 ± 5.9%) fractions among most patients, possibly due to radiation-induced edema. For all patients, the translational shifts in left-right, anteroposterior, and superoinferior directions were 2.1 ± 2.3 mm, 2.9 ± 2.8 mm, and 6.4 ± 5.5 mm, respectively. The greatest translational shifts occurred in the superoinferior direction, likely due to respiratory motion from the diaphragm. The rotational angles in roll, pitch, and yaw were 1.2° ± 1.8°, 1.8° ± 2.4°, and 1.7° ± 2.1°, respectively. The 3D individual fiducial displacements with rigid corrections were 0.2 ± 0.2 mm and increased to 0.5 ± 0.4 mm without rotational corrections. Conclusions: Accurate 3D locations of internal fiducials can be reconstructed from stereo imaging during treatment. As an
Chen, Chia-Hsiung; Azari, David; Hu, Yu Hen; Lindstrom, Mary J.; Thelen, Darryl; Yen, Thomas Y.; Radwin, Robert G.
2015-01-01
Objective Marker-less 2D video tracking was studied as a practical means to measure upper limb kinematics for ergonomics evaluations. Background Hand activity level (HAL) can be estimated from speed and duty cycle. Accuracy was measured using a cross correlation template-matching algorithm for tracking a region of interest on the upper extremities. Methods Ten participants performed a paced load transfer task while varying HAL (2, 4, and 5) and load (2.2 N, 8.9 N and 17.8 N). Speed and acceleration measured from 2D video were compared against ground truth measurements using 3D infrared motion capture. Results The median absolute difference between 2D video and 3D motion capture was 86.5 mm/s for speed, and 591 mm/s2 for acceleration, and less than 93 mm/s for speed and 656 mm/s2 for acceleration when camera pan and tilt were within ±30 degrees. Conclusion Single-camera 2D video had sufficient accuracy (< 100 mm/s) for evaluating HAL. Practitioner Summary This study demonstrated that 2D video tracking had sufficient accuracy to measure HAL for ascertaining the American Conference of Government Industrial Hygienists Threshold Limit Value® for repetitive motion when the camera is located within ±30 degrees off the plane of motion when compared against 3D motion capture for a simulated repetitive motion task. PMID:25978764
Chen, Chia-Hsiung; Azari, David P; Hu, Yu Hen; Lindstrom, Mary J; Thelen, Darryl; Yen, Thomas Y; Radwin, Robert G
2015-01-01
Marker-less 2D video tracking was studied as a practical means to measure upper limb kinematics for ergonomics evaluations. Hand activity level (HAL) can be estimated from speed and duty cycle. Accuracy was measured using a cross-correlation template-matching algorithm for tracking a region of interest on the upper extremities. Ten participants performed a paced load transfer task while varying HAL (2, 4, and 5) and load (2.2 N, 8.9 N and 17.8 N). Speed and acceleration measured from 2D video were compared against ground truth measurements using 3D infrared motion capture. The median absolute difference between 2D video and 3D motion capture was 86.5 mm/s for speed, and 591 mm/s(2) for acceleration, and less than 93 mm/s for speed and 656 mm/s(2) for acceleration when camera pan and tilt were within ± 30 degrees. Single-camera 2D video had sufficient accuracy (< 100 mm/s) for evaluating HAL. Practitioner Summary: This study demonstrated that 2D video tracking had sufficient accuracy to measure HAL for ascertaining the American Conference of Government Industrial Hygienists Threshold Limit Value(®) for repetitive motion when the camera is located within ± 30 degrees off the plane of motion when compared against 3D motion capture for a simulated repetitive motion task. PMID:25978764
Flexibility and rigidity of cross-linked Straight Fibrils under axial motion constraints.
Nagy Kem, Gyula
2016-09-01
The Straight Fibrils are stiff rod-like filaments and play a significant role in cellular processes as structural stability and intracellular transport. Introducing a 3D mechanical model for the motion of braced cylindrical fibrils under axial motion constraint; we provide some mechanism and a graph theoretical model for fibril structures and give the characterization of the flexibility and the rigidity of this bar-and-joint spatial framework. The connectedness and the circuit of the bracing graph characterize the flexibility of these structures. In this paper, we focus on the kinematical properties of hierarchical levels of fibrils and evaluate the number of the bracing elements for the rigidity and its computational complexity. The presented model is a good characterization of the frameworks of bio-fibrils such as microtubules, cellulose, which inspired this work. PMID:27289214
New conditional integrable cases of motion of a rigid body with Kovalevskaya's configuration
NASA Astrophysics Data System (ADS)
Yehia, H. M.; Elmandouh, A. A.
2011-01-01
We consider the general problem of motion of a rigid body about a fixed point under the action of an axisymmetric combination of potential and gyroscopic forces. We introduce two new cases of this problem which are integrable on the zero level of the cyclic integral. The new cases are combined generalizations of several previously known ones, namely those of Kovalevskaya, Yehia, Sokolov, Yehia and Bedweihi and Goriatchev, by the introduction of additional parameters to the structure of each.
Stabilization of the rotational motion of a rigid body on a vibrating base
NASA Astrophysics Data System (ADS)
Krementulo, V. V.
1984-12-01
The problem of the optimum (in a certain sense) stabilization of the permanent rotation of a heavy rigid body on a vibrating base is solved in the context of analytical control theory. Stabilization is achieved by means of a gimbal-suspended balanced gyroscope controlled by three moments. The control moments, obtained in explicit form, ensure the asymptotic stability of the rotational motion of the body along all of its phase coordinates and a minimum of a certain integral functional.
Flap-lag equations of motion of rigid, articulated rotor blades with three hinge sequences
NASA Technical Reports Server (NTRS)
Chen, Robert T. N.
1987-01-01
A derivation of coupled flap-lag equations of motion for a rigid articulated rotor with hinge springs and viscous dampers is reported. Three different flapping-lag-pitch hinge sequences are considered and the Lagrange method is used to derive the equations. The effects of the complete six degrees-of-freedom aircraft motions are included and all the inertia dynamic terms are retained; no small-angle assumptions are used in the development. Comparisons of the results with those available in the literature are made. Sources of terms missing in previous analyses, especially those of the inertia dynamics, are identified.
Computing 3-D structure of rigid objects using stereo and motion
NASA Technical Reports Server (NTRS)
Nguyen, Thinh V.
1987-01-01
Work performed as a step toward an intelligent automatic machine vision system for 3-D imaging is discussed. The problem considered is the quantitative 3-D reconstruction of rigid objects. Motion and stereo are the two clues considered in this system. The system basically consists of three processes: the low level process to extract image features, the middle level process to establish the correspondence in the stereo (spatial) and motion (temporal) modalities, and the high level process to compute the 3-D coordinates of the corner points by integrating the spatial and temporal correspondences.
Identification of motion parameters of a rigid body from its orthogonal and perspective projections
NASA Technical Reports Server (NTRS)
Ganguly, S.; Ghosh, B.; Tarn, T. J.; Bejczy, A. K.
1989-01-01
An estimate is made of the motion parameters, namely, linear and angular velocities, of a rigid body rotating and translating in three-space. The authors assume that the velocities are constant and that the motion is not completely observable. They consider two separate cases of partial observations corresponding to the orthogonal and the perspective projections, respectively. If (x, y, z) is the Cartesian coordinate of the three-space, the authors assume in the first case that the projection of the motion on the x-y plane is observed. If (r, theta, phi) is the polar coordinates of the three-space, they assume in the second case that the parameter vector (theta, phi) is observed. The use of both of these cases to estimate the motion parameters is discussed.
Time-correlation analysis of simulated water motion in flexible and rigid gramicidin channels.
Chiu, S W; Jakobsson, E; Subramaniam, S; McCammon, J A
1991-01-01
Molecular dynamics simulations have been done on a system consisting of the polypeptide membrane channel former gramicidin, plus water molecules in the channel and caps of waters at the two ends of the channel. In the absence of explicit simulation of the surrounding membrane, the helical form of the channel was maintained by artificial restraints on the peptide motion. The characteristic time constant of the artificial restraint was varied to assess the effect of the restraints on the channel structure and water motions. Time-correlation analysis was done on the motions of individual channel waters and on the motions of the center of mass of the channel waters. It is found that individual water molecules confined in the channel execute higher frequency motions than bulk water, for all degrees of channel peptide restraint. The center-of-mass motion of the chain of channel waters (which is the motion that is critical for transmembrane transport, due to the mandatory single filing of water in the channel) does not exhibit these higher frequency motions. The mobility of the water chain is dramatically reduced by holding the channel rigid. Thus permeation through the channel is not like flow through a rigid pipe; rather permeation is facilitated by peptide motion. For the looser restraints we used, the mobility of the water chain was not very much affected by the degree of restraint. Depending on which set of experiments is considered, the computed mobility of our water chain in the flexible channel is four to twenty times too high to account for the experimentally measured resistance of the gramicidin channel to water flow. From this result it appears likely that the peptide motions of an actual gramicidin channel embedded in a lipid membrane may be more restrained than in our flexible channel model, and that these restraints may be a significant modulator of channel permeability. For the completely rigid channel model the "trapping" of the water molecules in
NASA Technical Reports Server (NTRS)
Likins, P. W.
1974-01-01
Equations of motion are derived for use in simulating a spacecraft or other complex electromechanical system amenable to idealization as a set of hinge-connected rigid bodies of tree topology, with rigid axisymmetric rotors and nonrigid appendages attached to each rigid body in the set. In conjunction with a previously published report on finite-element appendage vibration equations, this report provides a complete minimum-dimension formulation suitable for generic programming for digital computer numerical integration.
NASA Technical Reports Server (NTRS)
1976-01-01
The two-particle, steady-state Schroedinger equation is transformed to center of mass and internuclear distance vector coordinates, leading to the free particle wave equation for the kinetic energy motion of the molecule and a decoupled wave equation for a single particle of reduced mass moving in a spherical potential field. The latter describes the vibrational and rotational energy modes of the diatomic molecule. For fixed internuclear distance, this becomes the equation of rigid rotator motion. The classical partition function for the rotator is derived and compared with the quantum expression. Molecular symmetry effects are developed from the generalized Pauli principle that the steady-state wave function of any system of fundamental particles must be antisymmetric. Nuclear spin and spin quantum functions are introduced and ortho- and para-states of rotators, along with their degeneracies, are defined. Effects of nuclear spin on entropy are deduced. Next, rigid polyatomic rotators are considered and the partition function for this case is derived. The patterns of rotational energy levels for nonlinear molecules are discussed for the spherical symmetric top, for the prolate symmetric top, for the oblate symmetric top, and for the asymmetric top. Finally, the equilibrium energy and specific heat of rigid rotators are derived.
Myocardial motion estimation in tagged MR sequences by using alphaMI-based non rigid registration.
Oubel, E; Tobon-Gomez, C; Hero, A O; Frangi, A F
2005-01-01
Tagged Magnetic Resonance Imaging (MRI) is currently the reference MR modality for myocardial motion and strain analysis. NMI-based non rigid registration has proven to be an accurate method to retrieve cardiac deformation fields. The use of alphaMI permits higher dimensional features to be implemented in myocardial deformation estimation through image registration. This paper demonstrates that this is feasible with a set of Haar wavelet features of high dimension. While we do not demonstrate performance improvement for this set of features, there is no significant degradation as compared to implementing the registration method with the traditional NMI metric. We use Entropic Spanning Graphs (ESGs) to estimate the alphaMI of the wavelet feature vectors WFVs since this is not possible with histograms. To the best of our knowledge, this is the first time that ESGs are used for non rigid registration. PMID:16685969
Gradient-Driven Vortex Motion in Nonneutral Plasmas and Ideal 2D Fluids
NASA Astrophysics Data System (ADS)
Schecter, David A.
2000-10-01
Two-dimensional (2D) turbulent flows can relax to metastable patterns without dissipation of kinetic energy. This ``rapid'' relaxation has been observed in computer simulations of ideal 2D fluids, and more recently in experiments with pure electron plasmas, which can obey similar dynamics. The late stage of relaxation often involves small vortices moving in a larger ``background'' shear-flow.(X.P. Huang et al., Phys. Rev. Lett. 74), 4424 (1995). In time, positive vortices (rotating counter-clockwise) move to peaks in background vorticity, whereas negative vortices (rotating clockwise) move to minima.(C.G. Rossby, J. Mar. Res. 7), 175 (1948); C.H. Liu and L. Ting, Comp. & Fluids 15, 77 (1987). In general, the rate of this migration increases with the magnitude of the background vorticity gradient, whereas it decreases as the background shear intensifies.\\vspace12pt Positive and negative vortices can also be classified as either prograde or retrograde, depending on whether they rotate with or against the local background shear. Surprisingly, a retrograde vortex moves up or down a background vorticity gradient orders of magnitude faster than a prograde vortex of equal strength.(D.A. Schecter and D.H.E. Dubin, Phys. Rev. Lett. 83), 2191 (1999). An accurate expression for the velocity of a weak retrograde vortex is obtained from an analytic calculation, in which the response of the background flow to the vortex is linearized. However, this linear theory fails for prograde vortices of any strength. Interestingly, the velocity of a prograde vortex can be obtained from a simple estimate, which accounts for the nonlinear ``trapping'' of background fluid around the vortex. The analytic expressions for the velocities of both prograde and retrograde vortices are in good quantitative agreement with vortex-in-cell simulations, and with electron plasma experiments, when the background shear is below a critical level. When the ratio of background shear to background vorticity
Unsteady hydrodynamic effect of rotation on steady rigid-body motion
NASA Astrophysics Data System (ADS)
Bhattacharya, S.
2005-08-01
Owing to the inertial effect of the flow, an unsteady hydrodynamic force will act on a particle of arbitrary shape undergoing a steady rigid-body motion with small but finite Reynolds number if the axis of rotation of the particle is not its axis of rotational symmetry. Unsteady flow field is generated owing to such rotation of the body and as a result the particle experiences a time-dependent translational resistance. In this paper, we analyse this time-dependent hydrodynamic force and obtain its higher-order correction by systematically expanding the Navier Stokes equation in small Reynolds number.
NASA Astrophysics Data System (ADS)
Zhang, N.; Zong, Z.
2011-11-01
The dynamic elastic response of a floating ship hull girder to an underwater explosion bubble is normally composed of two parts: rigid-body motion and elastic deformation. However, the effects of rigid-body motion have consistently been neglected in the current literature based on the assumption that they are small. In this paper, our focus is on the study of rigid-body motion effects on the hull girder's elastic deformation, also known as the ‘whipping response’. A theory of interaction between a gas bubble and a hull girder is presented. The bubble dynamic equations combined with the bubble migration, free surface effect and drag force considerations are solved numerically using the Runge-Kutta method. The rigid-body and elastic responses of the hull that are induced by the impulsive pressure of a bubble are calculated using the methods presented herein. Two different examples of real ships are given to demonstrate the effect of rigid-body motion on whipping responses. The numerical results show that rigid-body motions reduce the amplitudes and vibration natural periods of the bending moments of the hull girder. These effects can be ignored for slender hulls, but must be taken into consideration for shorter/wider hulls so as not to underestimate the longitudinal strength.
Chen, Mingqing; Bai, Junjie; Siochi, R Alfredo C
2013-02-01
To present a new method of estimating 3D positions of the ipsi-lateral hemi-diaphragm apex (IHDA) from 2D projection images of mega-voltage cone beam CT (MVCBCT). The detection framework reconstructs a 3D volume from all the 2D projection images. An initial estimated 3D IHDA position is determined in this volume based on an imaging processing pipeline, including Otsu thresholding, connected component labeling and template matching. This initial position is then projected onto each 2D projection image to create a region of interest (ROI). To accurately detect the IHDA position in 2D projection space, two methods, dynamic Hough transform (DHT) and a tracking approach based on a joint probability density function (PDF) are developed. Both methods utilize a double-parabola model to fit the 2D diaphragm boundary. The 3D IHDA motion in the superior-inferior (SI) direction is estimated from the initial static 3D position and the detected 2D positions in projection space. The two Hough-based detection methods are tested on 35 MVCBCT scans from 15 patients. The detection is compared to manually identified IHDA positions in 2D projection space by three clinicians. An average and standard deviation of 4.252 ± 3.354 and 2.485 ± 1.750 mm was achieved for DHT and tracking-based approaches respectively, compared with the inter-expert variance among three experts of 1.822 ± 1.106 mm. Based on the results of the scans, the PDF tracking-based approach appears more robust than the DHT. The combination of the automatic ROI localization and the tracking-based approach is a quicker and more accurate method of extracting 3D IHDA motion from 2D projection images. PMID:23321998
Analytic Theory and Control of the Motion of Spinning Rigid Bodies
NASA Technical Reports Server (NTRS)
Tsiotras, Panagiotis
1993-01-01
Numerical simulations are often resorted to, in order to understand the attitude response and control characteristics of a rigid body. However, this approach in performing sensitivity and/or error analyses may be prohibitively expensive and time consuming, especially when a large number of problem parameters are involved. Thus, there is an important role for analytical models in obtaining an understanding of the complex dynamical behavior. In this dissertation, new analytic solutions are derived for the complete attitude motion of spinning rigid bodies, under minimal assumptions. Hence, we obtain the most general solutions reported in the literature so far. Specifically, large external torques and large asymmetries are included in the problem statement. Moreover, problems involving large angular excursions are treated in detail. A new tractable formulation of the kinematics is introduced which proves to be extremely helpful in the search for analytic solutions of the attitude history of such kinds of problems. The main utility of the new formulation becomes apparent however, when searching for feedback control laws for stabilization and/or reorientation of spinning spacecraft. This is an inherently nonlinear problem, where standard linear control techniques fail. We derive a class of control laws for spin axis stabilization of symmetric spacecraft using only two pairs of gas jet actuators. Practically, this could correspond to a spacecraft operating in failure mode, for example. Theoretically, it is also an important control problem which, because of its difficulty, has received little, if any, attention in the literature. The proposed control laws are especially simple and elegant. A feedback control law that achieves arbitrary reorientation of the spacecraft is also derived, using ideas from invariant manifold theory. The significance of this research is twofold. First, it provides a deeper understanding of the fundamental behavior of rigid bodies subject to body
Dynamics and cortical distribution of neural responses to 2D and 3D motion in human
McKee, Suzanne P.; Norcia, Anthony M.
2013-01-01
The perception of motion-in-depth is important for avoiding collisions and for the control of vergence eye-movements and other motor actions. Previous psychophysical studies have suggested that sensitivity to motion-in-depth has a lower temporal processing limit than the perception of lateral motion. The present study used functional MRI-informed EEG source-imaging to study the spatiotemporal properties of the responses to lateral motion and motion-in-depth in human visual cortex. Lateral motion and motion-in-depth displays comprised stimuli whose only difference was interocular phase: monocular oscillatory motion was either in-phase in the two eyes (lateral motion) or in antiphase (motion-in-depth). Spectral analysis was used to break the steady-state visually evoked potentials responses down into even and odd harmonic components within five functionally defined regions of interest: V1, V4, lateral occipital complex, V3A, and hMT+. We also characterized the responses within two anatomically defined regions: the inferior and superior parietal cortex. Even harmonic components dominated the evoked responses and were a factor of approximately two larger for lateral motion than motion-in-depth. These responses were slower for motion-in-depth and were largely independent of absolute disparity. In each of our regions of interest, responses at odd-harmonics were relatively small, but were larger for motion-in-depth than lateral motion, especially in parietal cortex, and depended on absolute disparity. Taken together, our results suggest a plausible neural basis for reduced psychophysical sensitivity to rapid motion-in-depth. PMID:24198326
Taylor, L.M.; Preece, D.S.
1989-07-01
The computer program DMC (Distinct Motion Code) determines the two-dimensional planar rigid body motion of an arbitrary number of spherical shaped particles. The code uses an explicit central difference time integration algorithm to calculate the motion of the particles. Contact constraints between the particles are enforced using the penalty method. Coulomb friction and viscous damping are included in the collisions. The explicit time integration is conditionally stable with a time increment size which is dependent on the mass of the smallest particle in the mesh and the penalty stiffness used for the contact forces. The code chooses the spring stiffness based on the Young's modulus and Poisson's ratio of the material. The ability to tie spheres in pairs with a constraint condition is included in the code. The code has been written in an extremely efficient manner with particular emphasis placed on vector processing. While this does not impose any restrictions on non-vector processing computers, it does provide extremely fast results on vector processing computers. A bucket sorting or boxing algorithm is used to reduce the number of comparisons which must be made between spheres to determine the contact pairs. The sorting algorithm is completely algebraic and contains no logical branching. 13 refs., 14 figs., 4 tabs.
NASA Astrophysics Data System (ADS)
Taylor, L. M.; Preece, D. S.
1989-07-01
The computer program Distinct Motion Code (DMC) determines the two-dimensional planar rigid body motion of an arbitrary number of spherical shaped particles. The code uses an explicit central difference time integration algorithm to calculate the motion of the particles. Contact constraints between the particles are enforced using the penalty method. Coulomb friction and viscous damping are included in the collisions. The explicit time integration is conditionally stable with a time increment size which is dependent on the mass of the smallest particle in the mesh and the penalty stiffness used for the contact forces. The code chooses the spring stiffness based on the Young's modulus and Poisson's ratio of the material. The ability to tie spheres in pairs with a constraint condition is included in the code. The code has been written in an extremely efficient manner with particular emphasis placed on vector processing. While this does not impose any restrictions on non-vector processing computers, it does provide extremely fast results on vector processing computers. A bucket sorting or boxing algorithm is used to reduce the number of comparisons which must be made between spheres to determine the contact pairs. The sorting algorithm is completely algebraic and contains no logical branching.
Hydrodynamic lubrication of rigid nonconformal contacts in combined rolling and normal motion
NASA Technical Reports Server (NTRS)
Ghosh, M. K.; Hamrock, B. J.; Brewe, D. E.
1984-01-01
A numerical solution to the problem of hydrodynamic lubrication of rigid point contacts with an isoviscous, incompressible lubricant was obtained. The hydrodynamic load-carrying capacity under unsteady (or dynamic) conditions arising from the combined effects of squeeze motion superposed upon the entraining motion was determined for both normal approach and separation. Superposed normal motion considerably increases net load-carrying capacity during normal approach and substantially reduces net load-carrying capacity during separation. Geometry was also found to have a significant influence on the dynamic load-carrying capacity. The ratio of dynamic to steady state load-carrying capacity increases with increasing geometry parameter for normal approach and decreases during separation. The cavitation (film rupture) boundary is also influenced significantly by the normal motion, moving downstream during approach and upstream during separation. For sufficiently high normal separation velocity the rupture boundary may even move upstream of the minimum-film-thickness position. Sixty-three cases were used to derive a functional relationship for the ratio of the dynamic to steady state load-carrying capacity in terms of the dimensionless normal velocity parameter (incorporating normal velocity, entraining velocity, and film thickness) and the geometry parameter.
Dynamic Linkages Between the Transition Zone & Surface Plate Motions in 2D Models of Subduction
NASA Astrophysics Data System (ADS)
Arredondo, K.; Billen, M. I.
2012-12-01
Descending subducted slabs affect both plate tectonics at the surface and overall mantle flow (e.g. Conrad and Lithgow-Bertelloni, 2002). For time-dependent numerical models, the potential evolution of these slabs, ranging from immediate penetration into the lower mantle to prior buckling and stagnation, are affected by parameters such as the plate age, the viscosity jump into the lower mantle, the presence of phase transitions, trench motion and the chosen governing equation approximation (e.g. Billen and Hirth, 2007). Similarly, the overall deviatoric stress within the slab, especially where modified by the phase transitions, may explain the uneven distribution of deep earthquakes with depth (e.g. Bina, 1997). Better understanding of these processes may arise from a more realistic 2-D model that is fully-dynamic, with an overriding plate, freely-moving trench, compositionally-layered slab and seven major phase transitions, in addition to using the compressible (TALA) form of the governing equations. Though the thermodynamic parameters of certain phase transitions may be uncertain, this study aims to test the latest data and encourage further mineralogical research. We will present fully-dynamic models, which explore the importance of the phase transitions, especially those that have been previously excluded such as the wadsleyite to ringwoodite and the pyroxene and garnet phase transitions. These phase transitions, coupled with the modeled compositionally distinct crust, harzburgite, and pyrolite lithosphere layers, may produce new large-scale dynamic behavior not seen in past numerical models, as well as stress variations within the slab related to deep slab seismicity. Feedback from the compositionally complex slab to the dynamic trench may provide further insight on the mechanics of slab stagnation and behavior in the upper and lower mantle. Billen, M. I., and G. Hirth, Rheologic controls on slab dynamics, Geochemistry, Geophysics and Geosystems, 8 (Q08012
Wu Jian; Kim, Minho; Peters, Jorg; Chung, Heeteak; Samant, Sanjiv S.
2009-12-15
Purpose: Rigid 2D-3D registration is an alternative to 3D-3D registration for cases where largely bony anatomy can be used for patient positioning in external beam radiation therapy. In this article, the authors evaluated seven similarity measures for use in the intensity-based rigid 2D-3D registration using a variation in Skerl's similarity measure evaluation protocol. Methods: The seven similarity measures are partitioned intensity uniformity, normalized mutual information (NMI), normalized cross correlation (NCC), entropy of the difference image, pattern intensity (PI), gradient correlation (GC), and gradient difference (GD). In contrast to traditional evaluation methods that rely on visual inspection or registration outcomes, the similarity measure evaluation protocol probes the transform parameter space and computes a number of similarity measure properties, which is objective and optimization method independent. The variation in protocol offers an improved property in the quantification of the capture range. The authors used this protocol to investigate the effects of the downsampling ratio, the region of interest, and the method of the digitally reconstructed radiograph (DRR) calculation [i.e., the incremental ray-tracing method implemented on a central processing unit (CPU) or the 3D texture rendering method implemented on a graphics processing unit (GPU)] on the performance of the similarity measures. The studies were carried out using both the kilovoltage (kV) and the megavoltage (MV) images of an anthropomorphic cranial phantom and the MV images of a head-and-neck cancer patient. Results: Both the phantom and the patient studies showed the 2D-3D registration using the GPU-based DRR calculation yielded better robustness, while providing similar accuracy compared to the CPU-based calculation. The phantom study using kV imaging suggested that NCC has the best accuracy and robustness, but its slow function value change near the global maximum requires a
Kim, Young-Keun; Kim, Kyung-Soo
2014-10-15
Maritime transportation demands an accurate measurement system to track the motion of oscillating container boxes in real time. However, it is a challenge to design a sensor system that can provide both reliable and non-contact methods of 6-DOF motion measurements of a remote object for outdoor applications. In the paper, a sensor system based on two 2D laser scanners is proposed for detecting the relative 6-DOF motion of a crane load in real time. Even without implementing a camera, the proposed system can detect the motion of a remote object using four laser beam points. Because it is a laser-based sensor, the system is expected to be highly robust to sea weather conditions.
NASA Astrophysics Data System (ADS)
Kim, Young-Keun; Kim, Kyung-Soo
2014-10-01
Maritime transportation demands an accurate measurement system to track the motion of oscillating container boxes in real time. However, it is a challenge to design a sensor system that can provide both reliable and non-contact methods of 6-DOF motion measurements of a remote object for outdoor applications. In the paper, a sensor system based on two 2D laser scanners is proposed for detecting the relative 6-DOF motion of a crane load in real time. Even without implementing a camera, the proposed system can detect the motion of a remote object using four laser beam points. Because it is a laser-based sensor, the system is expected to be highly robust to sea weather conditions.
Xu, Q; Hanna, G; Kubicek, G; Asbell, S; Chen, Y; LaCouture, T; Grimm, J; Pahlajani, N; Fan, J
2014-06-01
Purpose: To quantitatively evaluate rigid and nonrigid motion of liver tumors based on fiducial tracking in 3D by stereo imaging during CyberKnife SBRT. Methods: Twenty-five liver patients previously treated with three-fractions of SBRT were retrospectively recruited in this study. During treatment, the 3D locations of fiducials were reported by the CyberKnife system after two orthogonal kV X-ray images were taken and further validated by geometry derivations. A total of 5004 pairs of X-ray images acquired during the course of treatment for all the patients, were analyzed. For rigid motion, the rotational angles and translational shifts by aligning 3D fiducial groups in different image pairs after least-square fitting were reported. For nonrigid motion, the relative interfractional tumor shape variations were reported and correlated to the sum of inter-fiducial distances. The individual fiducial displacements were also reported after rigid corrections and without angle corrections. Results: The relative tumor volume variation indicated by the inter-fiducial distances demonstrated an increasing trend in the second (101.6±3.4%) and third fraction (101.2±5.6%) among most patients. The cause could be possibly due to radiation-induced edema. For all the patients, the translational shift was 8.1±5.7 mm, with shifts in LR, AP and SI were 2.1±2.4 mm, 2.8±2.9 mm and 6.7±5.1 mm, respectively. The greatest translation shift occurred in SI, mainly due the breathing motion of diaphragm The rotational angles were 1.1±1.7°, 1.9±2.6° and 1.6±2.2°, in roll, pitch, and yaw, respectively. The 3D fiducial displacement with rigid corrections were 0.2±0.2 mm and increased to 0.6±0.3 mm without rotational corrections. Conclusion: The fiducial locations in 3D can be precisely reconstructed from CyberKnife stereo imaging system during treatment. The fiducials provide close estimation of both rigid and nonrigid motion of .liver tumors. The reported data could be further
Evaluating angular deflections from the digital gradient sensing method with rigid-motion deleted
NASA Astrophysics Data System (ADS)
Zhang, Rui
2016-06-01
The digital gradient sensing method is used for measuring small angular deflections of light rays due to local stresses in transparent planar solids. The method is based on two-dimensional (2D) digital image correlation (DIC) to measure the angular deflection of light rays; however, when a specimen is subjected to loading, deformation measurement from DIC is not perfect because of the existence of small in-plane and out-of-plane motions of the test sample surface that occurred after loading. These disadvantages will lead to errors in the measured angular deflections. The influence of unavoidable in-plane and out-of-plane motions was discussed, and a method to eliminate the influence to show the pure stress gradient of polymethy methacrylate is demonstrated.
Petasecca, M. Newall, M. K.; Aldosari, A. H.; Fuduli, I.; Espinoza, A. A.; Porumb, C. S.; Guatelli, S.; Metcalfe, P.; Lerch, M. L. F.; Rosenfeld, A. B.; Booth, J. T.; Colvill, E.; Duncan, M.; Cammarano, D.; Carolan, M.; Oborn, B.; Perevertaylo, V.; Keall, P. J.
2015-06-15
Purpose: Spatial and temporal resolutions are two of the most important features for quality assurance instrumentation of motion adaptive radiotherapy modalities. The goal of this work is to characterize the performance of the 2D high spatial resolution monolithic silicon diode array named “MagicPlate-512” for quality assurance of stereotactic body radiation therapy (SBRT) and stereotactic radiosurgery (SRS) combined with a dynamic multileaf collimator (MLC) tracking technique for motion compensation. Methods: MagicPlate-512 is used in combination with the movable platform HexaMotion and a research version of radiofrequency tracking system Calypso driving MLC tracking software. The authors reconstruct 2D dose distributions of small field square beams in three modalities: in static conditions, mimicking the temporal movement pattern of a lung tumor and tracking the moving target while the MLC compensates almost instantaneously for the tumor displacement. Use of Calypso in combination with MagicPlate-512 requires a proper radiofrequency interference shielding. Impact of the shielding on dosimetry has been simulated by GEANT4 and verified experimentally. Temporal and spatial resolutions of the dosimetry system allow also for accurate verification of segments of complex stereotactic radiotherapy plans with identification of the instant and location where a certain dose is delivered. This feature allows for retrospective temporal reconstruction of the delivery process and easy identification of error in the tracking or the multileaf collimator driving systems. A sliding MLC wedge combined with the lung motion pattern has been measured. The ability of the MagicPlate-512 (MP512) in 2D dose mapping in all three modes of operation was benchmarked by EBT3 film. Results: Full width at half maximum and penumbra of the moving and stationary dose profiles measured by EBT3 film and MagicPlate-512 confirm that motion has a significant impact on the dose distribution. Motion
NASA Astrophysics Data System (ADS)
Merckel, Loic; Nishida, Toyoaki
In this paper, we introduce a method for recognizing a subject complex object in real world environment. We use a three dimensional model described by line segments of the object and the data provided by a three-axis orientation sensor attached to the video camera. We assume that existing methods for finding line features in the image allow at least one model line segment to be detected as a single continuous segment. The method consists of two main steps: generation of pose hypotheses and then evaluation of each pose in order to select the most appropriate one. The first stage is three-fold: model visibility, line matching and pose estimation; the second stage aims to rank the poses by evaluating the similarity between the projected model lines and the image lines. Furthermore, we propose an additional step that consists of refining the best candidate pose by using the Lie group formalism of spatial rigid motions. Such a formalism provides an efficient local parameterization of the set of rigid rotation via the exponential map. A set of experiments demonstrating the robustness of this approach is presented.
NASA Technical Reports Server (NTRS)
Damodaran, Murali
1988-01-01
Unsteady inviscid transonic flow over airfoils in arbitrary rigid body motion is analyzed numerically by solving the two-dimensional unsteady Euler equations in integral form using a finite volume scheme. The solution procedure is based on an explicit Runge-Kutta time-stepping scheme wherein the spatial terms are central-differenced and a combination of second- and fourth-differences in the flow variables are used to form the numerical dissipation terms to stabilize the scheme. Unsteady calculations are started from converged steady-state solutions as initial conditions. Nonreflective boundary conditions are imposed on the far-field boundaries. Results are presented and, where possible, validated against available numerical and experimental data for airfoils subjected to a step change in angle of attack, airfoils oscillating and plunging in transonic flow, and airfoils immersed in a time-varying free stream.
The acoustic far-field of rigid bodies in arbitrary motion
NASA Technical Reports Server (NTRS)
Farassat, F.
1974-01-01
The far-field sound produced by a rigid body in arbitrary motion, with shock discontinuities close to the body, is studied. The analysis is based on the work of Ffowcs Williams and Hawkings (1969). An expression for the far-field sound pressure is obtained in the form of surface and line integrals carried out over a contracting sphere and its intersection with the body and shock surfaces. It is also found that in addition to the quadrupole distribution, the discontinuities in Lighthill stress at the shock, the fluid stresses at the body surface, and the curvatures (principal and mean) of the body and shock surfaces contribute to the sound field. Two examples are worked out.
Unsteady transonic flow past airfoils in rigid-body motion. [UFLO5
Chang, I C
1981-03-01
With the aim of developing a fast and accurate computer code for predicting the aerodynamic forces needed for a flutter analysis, some basic concepts in computational transonics are reviewed. The unsteady transonic flow past airfoils in rigid body motion is adequately described by the potential flow equation as long as the boundary layer remains attached. The two dimensional unsteady transonic potential flow equation in quasilinear form with first order radiation boundary conditions is solved by an alternating direction implicit scheme in an airfoil attached sheared parabolic coordinate system. Numerical experiments show that the scheme is very stable and is able to resolve the higher nonlinear transonic effects for filter analysis within the context of an inviscid theory.
The Accuracy of Webcams in 2D Motion Analysis: Sources of Error and Their Control
ERIC Educational Resources Information Center
Page, A.; Moreno, R.; Candelas, P.; Belmar, F.
2008-01-01
In this paper, we show the potential of webcams as precision measuring instruments in a physics laboratory. Various sources of error appearing in 2D coordinate measurements using low-cost commercial webcams are discussed, quantifying their impact on accuracy and precision, and simple procedures to control these sources of error are presented.…
NASA Astrophysics Data System (ADS)
Lohr, M. B.
2008-10-01
The rotational motion of a torque-free axisymmetric rigid body is precession. This motion has been expressed analytically in the literature given the body's initial orientation and rotational dynamics parameters, i.e. inertia ratio and initial angular velocities or precession parameters. The inverse problem of deriving these dynamics parameters given orientation in time has been implemented numerically but has not yet been solved analytically. If a rigid body is precessing, and its orientation with respect to an arbitrary inertial frame is provided at three equally spaced points in time such that the rotational motion is not undersampled, an analytical inverse solution is presented for the precession rate, relative spin rate, coning angle and angular velocities; if the precessional motion is due to inertial axisymmetry and torque-free motion, the inertia ratio is also derived. Additionally, an analytical methodology is presented to test for non-precessional motion. These techniques are applicable to various problems in space science and astronomy, where non-precessional motion or the rotational dynamics parameters of this type of rigid body must be accurately derived from its orientation or relative orientation in time.
Vehicular motion in 2D city traffic network with signals controlled by phase shift
NASA Astrophysics Data System (ADS)
Komada, Kazuhito; Kojima, Kengo; Nagatani, Takashi
2011-03-01
We study the dynamic behavior of vehicular traffic through the series of traffic lights controlled by phase shift in two-dimensional (2D) city traffic network. The nonlinear-map model is presented for the vehicular traffic. The city traffic network is made of one-way perpendicular streets arranged in a square lattice with traffic signals where vertical streets are oriented upwards and horizontal streets are oriented rightwards. There are two traffic lights for the movement to north or that to east at each crossing. The traffic lights are controlled by the cycle time, split, and phase shift. The vehicle moves through the series of signals on a path selected by the driver. The city traffic with a heterogeneous density distribution is also studied. The dependence of the arrival time on cycle time, split, phase shift, selected path, and density is clarified for 2D city traffic. It is shown that the vehicular traffic is efficiently controlled by the phase shift.
Parker, Jason G; Mair, Bernard A; Gilland, David R
2009-10-01
In this article, a new method is introduced for estimating the motion of the heart due to respiration in gated cardiac SPECT using a rigid-body model with rotation parametrized by a unit quaternion. The method is based on minimizing the sum of squared errors between the reference and the deformed frames resulting from the usual optical flow constraint by using an optimized conjugate gradient routine. This method does not require any user-defined parameters or penalty terms, which simplifies its use in a clinical setting. Using a mathematical phantom, the method was quantitatively compared to the principal axis method, as well as an iterative method in which the rotation matrix was represented by Euler angles. The quaternion-based method was shown to be substantially more accurate and robust across a wide range of extramyocardial activity levels than the principal axis method. Compared with the Euler angle representation, the quaternion-based method resulted in similar accuracy but a significant reduction in computation times. Finally, the quaternion-based method was investigated using a respiratory-gated cardiac SPECT acquisition of a human subject. The motion-corrected image has increased sharpness and myocardial uniformity compared to the uncorrected image. PMID:19928105
NASA Technical Reports Server (NTRS)
Ghosh, M. K.; Brewe, D. E.; Hamrock, B. J.
1987-01-01
The effect of inlet starvation on the hydrodynamic lubrication of lightly loaded rigid nonconformal contacts in combined rolling and normal motion is determined through a numerical solution of the Reynolds' equation for an isoviscous, incompressible lubricant. Starvation is effected by systematically reducing the fluid inlet level. The pressures are taken to be ambient at the inlet meniscus boundary and Reynolds' boundary condition is applied for film rupture in the exit region. Results are presented for the dynamic performance of the starved contacts in combined rolling and normal motion for both normal approach and separation. During normal approach the dynamic load ratio (i.e. ratio of dynamic to steady state load capacity) increases considerably with increase in the inlet starvation. The effect of starvation on the dynamic peak pressure ratio is relatively small. Further, it has been observed that with increasing starvation, film thickness effects become significant in the dynamic behavior of the nonconformal contacts. For significantly starved contacts the dynamic load ratio increases with increase in film thickness during normal approach and a similar reduction is observed during separation. A similar effect is noted for the dynamic peak pressure ratio.
NASA Technical Reports Server (NTRS)
Ghosh, M. K.; Hamrock, B. J.; Brewe, D. E.
1986-01-01
The effect of inlet starvation on the hydrodynamic lubrication of lightly loaded rigid nonconformal contacts in combined rolling and normal motion is determined through a numerical solution of the Reynolds' equation for an isoviscous, incompressible lubricant. Starvation is effected by systematically reducing the fluid inlet level. The pressures are taken to be ambient at the inlet meniscus boundary and Reynolds' boundary condition is applied for film rupture in the exit region. Results are presented for the dynamic performance of the starved contacts in combined rolling and normal motion for both normal approach and separation. During normal approach the dynamic load ratio (i.e. ratio of dynamic to steady state load capacity) increases considerably with increase in the inlet starvation. The effect of starvation on the dynamic peak pressure ratio is relatively small. Further, it has been observed that with increasing starvation, film thickness effects become significant in the dynamic behavior of the nonconformal contacts. For significantly starved contacts the dynamic load ratio increases with increase in film thickness during normal approach and a similar reduction is observed during separation. A similar effect is noted for the dynamic peak pressure ratio.
A vector-dyadic development of the equations of motion for N-coupled rigid bodies and point masses
NASA Technical Reports Server (NTRS)
Frisch, H. P.
1974-01-01
The equations of motion are derived, in vector-dyadic format, for a topological tree of coupled rigid bodies, point masses, and symmetrical momentum wheels. These equations were programmed, and form the basis for the general-purpose digital computer program N-BOD. A complete derivation of the equations of motion is included along with a description of the methods used for kinematics, constraint elimination, and for the inclusion of nongyroscope forces and torques acting external or internal to the system.
GPS Constraints on Lesser Antilles Forearc Motion and Rigid Caribbean Plate
NASA Astrophysics Data System (ADS)
López, A. M.; Stein, S.; Sella, G.; Dixon, T. H.; Calais, E.; Jansma, P. E.
2005-05-01
We are using a decade of Global Positioning System data to address two tectonic problems of the Caribbean (CA) plate; 1) Whether a forearc sliver exists along the Lesser Antilles forearc and if so what is its dynamics and location, and 2) Whether the Caribbean plate is deforming internally. We approach this problem by developing GPS-derived velocity vectors at sites within the CA plate and its boundaries and comparing them to four decades of earthquake data. In a number of subduction zones, misfits between slip vectors and predicted convergence azimuths from Euler vectors suggest the presence of a forearc sliver, where trench-parallel motion is accommodated along a strike-slip fault system. Such a situation may be occurring at the eastern boundary of the CA plate along the Lesser Antilles (LA) forearc, where the North America (NA) plate subducts obliquely. Comparing slip vectors of shallow (0-60 km) thrust events to the predicted motions of GPS-based Euler vectors show a systematic northerly misfit, suggesting a trench-parallel component of motion taken up by the forearc sliver. This possibility can be tested with GPS data from the forearc. In addition, we use new GPS data to constrain the internal rigidity of the plate. Previous GPS work yielded a possible upper bound on internal deformation of 4-6 mm/yr. With an expansion in the data set on critically located stations in the CA plate (SANA, ROJO, CRO1 and AVES), we have computed new sets of Euler vector pairs for the CA-NA and CA-South America plate pairs.
Multi-level model for 2D human motion analysis and description
NASA Astrophysics Data System (ADS)
Foures, Thomas; Joly, Philippe
2003-01-01
This paper deals with the proposition of a model for human motion analysis in a video. Its main caracteristic is to adapt itself automatically to the current resolution, the actual quality of the picture, or the level of precision required by a given application, due to its possible decomposition into several hierarchical levels. The model is region-based to address some analysis processing needs. The top level of the model is only defined with 5 ribbons, which can be cut into sub-ribbons regarding to a given (or an expected) level of details. Matching process between model and current picture consists in the comparison of extracted subject shape with a graphical rendering of the model built on the base of some computed parameters. The comparison is processed by using a chamfer matching algorithm. In our developments, we intend to realize a platform of interaction between a dancer and tools synthetizing abstract motion pictures and music in the conditions of a real-time dialogue between a human and a computer. In consequence, we use this model in a perspective of motion description instead of motion recognition: no a priori gestures are supposed to be recognized as far as no a priori application is specially targeted. The resulting description will be made following a Description Scheme compliant with the movement notation called "Labanotation".
NASA Technical Reports Server (NTRS)
Fitzpatrick, P. M.; Harmon, G. R.; Liu, J. J. F.; Cochran, J. E.
1974-01-01
The formalism for studying perturbations of a triaxial rigid body within the Hamilton-Jacobi framework is developed. The motion of a triaxial artificial earth satellite about its center of mass is studied. Variables are found which permit separation, and the Euler angles and associated conjugate momenta are obtained as functions of canonical constants and time.
1 kHz 2D Visual Motion Sensor Using 20 × 20 Silicon Retina Optical Sensor and DSP Microcontroller.
Liu, Shih-Chii; Yang, MinHao; Steiner, Andreas; Moeckel, Rico; Delbruck, Tobi
2015-04-01
Optical flow sensors have been a long running theme in neuromorphic vision sensors which include circuits that implement the local background intensity adaptation mechanism seen in biological retinas. This paper reports a bio-inspired optical motion sensor aimed towards miniature robotic and aerial platforms. It combines a 20 × 20 continuous-time CMOS silicon retina vision sensor with a DSP microcontroller. The retina sensor has pixels that have local gain control and adapt to background lighting. The system allows the user to validate various motion algorithms without building dedicated custom solutions. Measurements are presented to show that the system can compute global 2D translational motion from complex natural scenes using one particular algorithm: the image interpolation algorithm (I2A). With this algorithm, the system can compute global translational motion vectors at a sample rate of 1 kHz, for speeds up to ±1000 pixels/s, using less than 5 k instruction cycles (12 instructions per pixel) per frame. At 1 kHz sample rate the DSP is 12% occupied with motion computation. The sensor is implemented as a 6 g PCB consuming 170 mW of power. PMID:25879969
The oscillatory motion of a surfactant-laden liquid plug in a 2D-channel
NASA Astrophysics Data System (ADS)
Fujioka, Hideki; Grotberg, James B.
2004-11-01
Liquid plugs can form in the lung's small airways near the end of expiration. This happens more frequently when the amount of pulmonary surfactant is reduced. In medical treatments such as surfactant replacement therapy, partial liquid ventilation, and drug delivery, the formation of plugs in an airway is important to deliver the instilled liquid uniformly throughout the lung. In this study, we investigate numerically the oscillatory motion of a surfactant-laden liquid plug within a two-dimensional channel lined by a thin liquid film. The viscosity of both the left and right air phases is assumed to be negligible, so that the only fluid dynamics of the liquid phase is considered. The plug motion is regulated by the flow rate in the left air phase, which is prescribed as a sinusoidal function of time. The pressure drop between the left and right air phases varies for time with a different phase of the flow rate. The plug length and the film thickness oscillate with an average value during a cycle. These behaviors changes by system parameters, Reynolds number, Womersley number, Capillary number, and surfactant properties. The significance of this study on mechanical stresses acting on airway epithelial cells caused by the motion of a liquid plug during normal breath, conventional or high-frequency ventilation is discussed. Supported by NIH grant HL41126, NASA grant NAG3-2740.
Arterial Mechanical Motion Estimation Based on a Semi-Rigid Body Deformation Approach
Guzman, Pablo; Hamarneh, Ghassan; Ros, Rafael; Ros, Eduardo
2014-01-01
Arterial motion estimation in ultrasound (US) sequences is a hard task due to noise and discontinuities in the signal derived from US artifacts. Characterizing the mechanical properties of the artery is a promising novel imaging technique to diagnose various cardiovascular pathologies and a new way of obtaining relevant clinical information, such as determining the absence of dicrotic peak, estimating the Augmentation Index (AIx), the arterial pressure or the arterial stiffness. One of the advantages of using US imaging is the non-invasive nature of the technique unlike Intra Vascular Ultra Sound (IVUS) or angiography invasive techniques, plus the relative low cost of the US units. In this paper, we propose a semi rigid deformable method based on Soft Bodies dynamics realized by a hybrid motion approach based on cross-correlation and optical flow methods to quantify the elasticity of the artery. We evaluate and compare different techniques (for instance optical flow methods) on which our approach is based. The goal of this comparative study is to identify the best model to be used and the impact of the accuracy of these different stages in the proposed method. To this end, an exhaustive assessment has been conducted in order to decide which model is the most appropriate for registering the variation of the arterial diameter over time. Our experiments involved a total of 1620 evaluations within nine simulated sequences of 84 frames each and the estimation of four error metrics. We conclude that our proposed approach obtains approximately 2.5 times higher accuracy than conventional state-of-the-art techniques. PMID:24871987
Arterial mechanical motion estimation based on a semi-rigid body deformation approach.
Guzman, Pablo; Hamarneh, Ghassan; Ros, Rafael; Ros, Eduardo
2014-01-01
Arterial motion estimation in ultrasound (US) sequences is a hard task due to noise and discontinuities in the signal derived from US artifacts. Characterizing the mechanical properties of the artery is a promising novel imaging technique to diagnose various cardiovascular pathologies and a new way of obtaining relevant clinical information, such as determining the absence of dicrotic peak, estimating the Augmentation Index (AIx), the arterial pressure or the arterial stiffness. One of the advantages of using US imaging is the non-invasive nature of the technique unlike Intra Vascular Ultra Sound (IVUS) or angiography invasive techniques, plus the relative low cost of the US units. In this paper, we propose a semi rigid deformable method based on Soft Bodies dynamics realized by a hybrid motion approach based on cross-correlation and optical flow methods to quantify the elasticity of the artery. We evaluate and compare different techniques (for instance optical flow methods) on which our approach is based. The goal of this comparative study is to identify the best model to be used and the impact of the accuracy of these different stages in the proposed method. To this end, an exhaustive assessment has been conducted in order to decide which model is the most appropriate for registering the variation of the arterial diameter over time. Our experiments involved a total of 1620 evaluations within nine simulated sequences of 84 frames each and the estimation of four error metrics. We conclude that our proposed approach obtains approximately 2.5 times higher accuracy than conventional state-of-the-art techniques. PMID:24871987
Binocular Perception of 2D Lateral Motion and Guidance of Coordinated Motor Behavior.
Fath, Aaron J; Snapp-Childs, Winona; Kountouriotis, Georgios K; Bingham, Geoffrey P
2016-04-01
Zannoli, Cass, Alais, and Mamassian (2012) found greater audiovisual lag between a tone and disparity-defined stimuli moving laterally (90-170 ms) than for disparity-defined stimuli moving in depth or luminance-defined stimuli moving laterally or in depth (50-60 ms). We tested if this increased lag presents an impediment to visually guided coordination with laterally moving objects. Participants used a joystick to move a virtual object in several constant relative phases with a laterally oscillating stimulus. Both the participant-controlled object and the target object were presented using a disparity-defined display that yielded information through changes in disparity over time (CDOT) or using a luminance-defined display that additionally provided information through monocular motion and interocular velocity differences (IOVD). Performance was comparable for both disparity-defined and luminance-defined displays in all relative phases. This suggests that, despite lag, perception of lateral motion through CDOT is generally sufficient to guide coordinated motor behavior. PMID:26614099
Energy Exchange during Plunge/Surge Motions of a 2D Wing
NASA Astrophysics Data System (ADS)
Kerstens, Wesley; Choi, Jeesoon; Colonius, Tim; Williams, David
2011-11-01
The rate of energy transfer between an NACA-0006 wing and an unsteady flow is examined at pre-stall and post-stall conditions using numerical simulations and wind tunnel experiments. The plunge and surge motions simulate the fluctuating vertical (wz) and longitudinal (wx) velocity components of a wind gust. In a steady flow the wing loses energy to the flow through the drag power term, but in an unsteady flow the wing may gain energy from the fluctuating lift power and fluctuating drag power terms. The net energy transfer averaged over the period of oscillation depends on the phase angle between the plunge and surge motions. The largest increase of energy occurs when wx and wz are in-phase. When the fluctuations are large enough, then it is possible for the net energy gain to be positive. The numerical simulations conducted at Reynolds numbers near the critical value for vortex shedding show qualitative agreement with the experiments. The simulations highlight the role of vortex shedding in determining the optimal frequency and phase for energy extraction from the gust. Support of the AFOSR through grant FA9550-09-1-0189 managed by Dr. Douglas Smith is gratefully acknnowledged.
On the relative rotational motion between rigid fibers and fluid in turbulent channel flow
NASA Astrophysics Data System (ADS)
Marchioli, C.; Zhao, L.; Andersson, H. I.
2016-01-01
In this study, the rotation of small rigid fibers relative to the surrounding fluid in wall-bounded turbulence is examined by means of direct numerical simulations coupled with Lagrangian tracking. Statistics of the relative (fiber-to-fluid) angular velocity, referred to as slip spin in the present study, are evaluated by modelling fibers as prolate spheroidal particles with Stokes number, St, ranging from 1 to 100 and aspect ratio, λ, ranging from 3 to 50. Results are compared one-to-one with those obtained for spherical particles (λ = 1) to highlight effects due to fiber length. The statistical moments of the slip spin show that differences in the rotation rate of fibers and fluid are influenced by inertia, but depend strongly also on fiber length: Departures from the spherical shape, even when small, are associated with an increase of rotational inertia and prevent fibers from passively following the surrounding fluid. An increase of fiber length, in addition, decouples the rotational dynamics of a fiber from its translational dynamics suggesting that the two motions can be modelled independently only for long enough fibers (e.g., for aspect ratios of order ten or higher in the present simulations).
Inertial Motions of a Rigid Body with a Cavity Filled with a Viscous Liquid
NASA Astrophysics Data System (ADS)
Disser, Karoline; Galdi, Giovanni P.; Mazzone, Giusy; Zunino, Paolo
2016-07-01
We study inertial motions of the coupled system, {S}, constituted by a rigid body containing a cavity entirely filled with a viscous liquid. We show that for arbitrary initial data having only finite kinetic energy, every corresponding weak solution (à la Leray-Hopf) converges, as time goes to infinity, to a uniform rotation, unless two central moments of inertia of {S} coincide and are strictly greater than the third one. This corroborates a famous "conjecture" of N.Ye. Zhukovskii in several physically relevant cases. Moreover, we show that, in a known range of initial data, this rotation may only occur along the central axis of inertia of {S} with the larger moment of inertia. We also provide necessary and sufficient conditions for the rigorous nonlinear stability of permanent rotations, which improve and/or generalize results previously given by other authors under different types of approximation. Finally, we present results obtained by a targeted numerical simulation that, on the one hand, complement the analytical findings, whereas, on the other hand, point out new features that the analysis is yet not able to catch, and, as such, lay the foundation for interesting and challenging future investigation.
NASA Astrophysics Data System (ADS)
Khoshelham, Kourosh
2016-04-01
Registration is often a prerequisite step in processing point clouds. While planar surfaces are suitable features for registration, most of the existing plane-based registration methods rely on iterative solutions for the estimation of transformation parameters from plane correspondences. This paper presents a new closed-form solution for the estimation of a rigid motion from a set of point-plane correspondences. The role of normalization is investigated and its importance for accurate plane fitting and plane-based registration is shown. The paper also presents a thorough evaluation of the closed-form solutions and compares their performance with the iterative solution in terms of accuracy, robustness, stability and efficiency. The results suggest that the closed-form solution based on point-plane correspondences should be the method of choice in point cloud registration as it is significantly faster than the iterative solution, and performs as well as or better than the iterative solution in most situations. The normalization of the point coordinates is also recommended as an essential preprocessing step for point cloud registration. An implementation of the closed-form solutions in MATLAB is available at: http://people.eng.unimelb.edu.au/kkhoshelham/research.html#directmotion
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
Towards real-time 2D/3D registration for organ motion monitoring in image-guided radiation therapy
NASA Astrophysics Data System (ADS)
Gendrin, C.; Spoerk, J.; Bloch, C.; Pawiro, S. A.; Weber, C.; Figl, M.; Markelj, P.; Pernus, F.; Georg, D.; Bergmann, H.; Birkfellner, W.
2010-02-01
Nowadays, radiation therapy systems incorporate kV imaging units which allow for the real-time acquisition of intra-fractional X-ray images of the patient with high details and contrast. An application of this technology is tumor motion monitoring during irradiation. For tumor tracking, implanted markers or position sensors are used which requires an intervention. 2D/3D intensity based registration is an alternative, non-invasive method but the procedure must be accelerate to the update rate of the device, which lies in the range of 5 Hz. In this paper we investigate fast CT to a single kV X-ray 2D/3D image registration using a new porcine reference phantom with seven implanted fiducial markers. Several parameters influencing the speed and accuracy of the registrations are investigated. First, four intensity based merit functions, namely Cross-Correlation, Rank Correlation, Mutual Information and Correlation Ratio, are compared. Secondly, wobbled splatting and ray casting rendering techniques are implemented on the GPU and the influence of each algorithm on the performance of 2D/3D registration is evaluated. Rendering times for a single DRR of 20 ms were achieved. Different thresholds of the CT volume were also examined for rendering to find the setting that achieves the best possible correspondence with the X-ray images. Fast registrations below 4 s became possible with an inplane accuracy down to 0.8 mm.
Modeling Selective Local Interactions with Memory: Motion on a 2D Lattice.
Weinberg, Daniel; Levy, Doron
2014-06-15
We consider a system of particles that simultaneously move on a two-dimensional periodic lattice at discrete times steps. Particles remember their last direction of movement and may either choose to continue moving in this direction, remain stationary, or move toward one of their neighbors. The form of motion is chosen based on predetermined stationary probabilities. Simulations of this model reveal a connection between these probabilities and the emerging patterns and size of aggregates. In addition, we develop a reaction diffusion master equation from which we derive a system of ODEs describing the dynamics of the particles on the lattice. Simulations demonstrate that solutions of the ODEs may replicate the aggregation patterns produced by the stochastic particle model. We investigate conditions on the parameters that influence the locations at which particles prefer to aggregate. This work is a two-dimensional generalization of [Galante & Levy, Physica D, http://dx.doi.org/10.1016/j.physd.2012.10.010], in which the corresponding one-dimensional problem was studied. PMID:25045193
Dynamic Linkages Between the Transition Zone & Surface Plate Motions in 2D Models of Subduction
NASA Astrophysics Data System (ADS)
Arredondo, K.; Billen, M. I.
2013-12-01
While slab pull is considered the dominant force controlling plate motion and speed, its magnitude is controlled by slab behavior in the mantle, where tomographic studies show a wide range of possibilities from direct penetration to folding, or stagnation directly above the lower mantle (e.g. Fukao et al., 2009). Geodynamic studies have investigated various parameters, such as plate age and two phase transitions, to recreate observed behavior (e.g. Běhounková and Cízková, 2008). However, past geodynamic models have left out known slab characteristics that may have a large impact on slab behavior and our understanding of subduction processes. Mineral experiments and seismic observations have indicated the existence of additional phase transitions in the mantle transition zone that may produce buoyancy forces large enough to affect the descent of a subducting slab (e.g. Ricard et al., 2005). The current study systematically tests different common assumptions used in geodynamic models: kinematic versus free-slip boundary conditions, the effects of adiabatic heating, viscous dissipation and latent heat, compositional layering and a more complete suite of phase transitions. Final models have a complete energy equation, with eclogite, harzburgite and pyrolite lithosphere compositional layers, and seven composition-dependent phase transitions within the olivine, pyroxene and garnet polymorph minerals. Results show important feedback loops between different assumptions and new behavior from the most complete models. Kinematic models show slab weakening or breaking above the 660 km boundary and between compositional layers. The behavior in dynamic models with a free-moving trench and overriding plate is compared to the more commonly found kinematic models. The new behavior may have important implications for the depth distribution of deep earthquakes within the slab. Though the thermodynamic parameters of certain phase transitions may be uncertain, their presence and
Rigidity and definition of Caribbean plate motion from COCONet and campaign GPS observations
NASA Astrophysics Data System (ADS)
Mattioli, Glen; Miller, Jamie; DeMets, Charles; Jansma, Pamela
2014-05-01
observations, which implies that the Caribbean is undergoing modest (1-3 mm/yr) deformation within its interior. Some sites in the western Caribbean included in our analysis may be biased by small, but significant coseismic deformation, which has not been removed from the site velocities used in our inversion to define Caribbean motion and rigidity. Scenarios for possible east-west deformation accommodated across the Lower Nicaraguan Rise and Beata Ridge will be presented.
Rigidity and definition of Caribbean plate motion from COCONet and campaign GPS observations
NASA Astrophysics Data System (ADS)
Mattioli, G. S.; Miller, J. A.; DeMets, C.; Jansma, P. E.
2015-12-01
sites in the western Caribbean included in our analysis may be biased by small, but significant coseismic deformation, which has not been removed from the site velocities used in our inversion to define Caribbean motion and rigidity. Scenarios for possible east-west deformation accommodated across the Lower Nicaraguan Rise and Beata Ridge will be presented.
NASA Astrophysics Data System (ADS)
Wu, Shunguang; Hong, Lang
2008-04-01
A framework of simultaneously estimating the motion and structure parameters of a 3D object by using high range resolution (HRR) and ground moving target indicator (GMTI) measurements with template information is given. By decoupling the motion and structure information and employing rigid-body constraints, we have developed the kinematic and measurement equations of the problem. Since the kinematic system is unobservable by using only one scan HRR and GMTI measurements, we designed an architecture to run the motion and structure filters in parallel by using multi-scan measurements. Moreover, to improve the estimation accuracy in large noise and/or false alarm environments, an interacting multi-template joint tracking (IMTJT) algorithm is proposed. Simulation results have shown that the averaged root mean square errors for both motion and structure state vectors have been significantly reduced by using the template information.
NASA Astrophysics Data System (ADS)
Njoroge, M. W.; Malservisi, R.; Hugentobler, U.; Mokhtari, M.; Voytenko, D.
2014-12-01
Plate rigidity is one of the main paradigms of plate tectonics and a fundamental assumption in the definition of a global reference frame as ITRF. Although still far for optimal, the increased GPS instrumentation of the African region can allow us to understand how rigid one of the major plate can be. The presence of diffused band of seismicity, the Cameroon volcanic line, Pan African Kalahari orogenic belt and East Africa Rift suggest the possibility of relative motion among the different regions within the Nubia. The study focuses on the rigidity of Nubia plate. We divide the plate into three regions: Western (West Africa craton plus Nigeria), Central (approximately the region of the Congo craton) and Southern (Kalahari craton plus South Africa) and we utilize Euler Vector formulation to study internal rigidity and eventual relative motion. Developing five different reference frames with different combinations of the 3 regions, we try to understand the presence of the relative motion between the 3 cratons thus the stability of the Nubia plate as a whole. All available GPS stations from the regions are used separately or combined in creation of the reference frames. We utilize continuous stations with at least 2.5 years of data between 1994 and 2014. Given the small relative velocity, it is important to eliminate eventual biases in the analysis and to have a good estimation in the uncertainties of the observed velocities. For this reason we perform our analysis using both Bernese and Gipsy-oasis codes to generate time series for each station. Velocities and relative uncertainties are analyzed using the Allan variance of rate technique, taking in account for colored noise. An analysis of the color of the noise as function of latitude and climatic region is also performed to each time series. Preliminary results indicate a slight counter clockwise motion of West Africa craton with respect to South Africa Kalahari, and South Africa Kalahari-Congo Cratons. In addition
Technology Transfer Automated Retrieval System (TEKTRAN)
This paper presents a depth-averaged two-dimensional shallow water model for simulating long waves in vegetated water bodies under breaking and non-breaking conditions. The effects of rigid vegetation are modelled in the form of drag and inertia forces as sink terms in the momentum equations. The dr...
Efficient framework for deformable 2D-3D registration
NASA Astrophysics Data System (ADS)
Fluck, Oliver; Aharon, Shmuel; Khamene, Ali
2008-03-01
Using 2D-3D registration it is possible to extract the body transformation between the coordinate systems of X-ray and volumetric CT images. Our initial motivation is the improvement of accuracy of external beam radiation therapy, an effective method for treating cancer, where CT data play a central role in radiation treatment planning. Rigid body transformation is used to compute the correct patient setup. The drawback of such approaches is that the rigidity assumption on the imaged object is not valid for most of the patient cases, mainly due to respiratory motion. In the present work, we address this limitation by proposing a flexible framework for deformable 2D-3D registration consisting of a learning phase incorporating 4D CT data sets and hardware accelerated free form DRR generation, 2D motion computation, and 2D-3D back projection.
ERIC Educational Resources Information Center
Bellver-Cebreros, Consuelo; Rodriguez-Danta, Marcelo
2009-01-01
An apparently unnoticed analogy between the torque-free motion of a rotating rigid body about a fixed point and the propagation of light in anisotropic media is stated. First, a new plane construction for visualizing this torque-free motion is proposed. This method uses an intrinsic representation alternative to angular momentum and independent of…
Lorenz, K S; Salama, P; Dunn, K W; Delp, E J
2012-02-01
Digital image analysis is a fundamental component of quantitative microscopy. However, intravital microscopy presents many challenges for digital image analysis. In general, microscopy volumes are inherently anisotropic, suffer from decreasing contrast with tissue depth, lack object edge detail and characteristically have low signal levels. Intravital microscopy introduces the additional problem of motion artefacts, resulting from respiratory motion and heartbeat from specimens imaged in vivo. This paper describes an image registration technique for use with sequences of intravital microscopy images collected in time-series or in 3D volumes. Our registration method involves both rigid and nonrigid components. The rigid registration component corrects global image translations, whereas the nonrigid component manipulates a uniform grid of control points defined by B-splines. Each control point is optimized by minimizing a cost function consisting of two parts: a term to define image similarity, and a term to ensure deformation grid smoothness. Experimental results indicate that this approach is promising based on the analysis of several image volumes collected from the kidney, lung and salivary gland of living rodents. PMID:22092443
NASA Astrophysics Data System (ADS)
Jokinen, Olli
2013-10-01
The paper deals with measurement of human facial deformations from synchronized image sequences taken with multiple calibrated cameras from different viewpoints. SIFT (Scale Invariant Feature Transform) keypoints are utilized as image feature points in the first place to determine spatial and temporal correspondences between images. If no temporal match is found for an image point by keypoint matching, then the tracking of the point is switched to least squares matching provided the point has one or more spatial corresponding points in the other views of the previous frame. For this purpose, a new method based on affine multi-image least squares matching is proposed where multiple spatial and temporal template images are simultaneously matched against each search image and part of the spatial template images also change during adjustment. A new method based on analyzing temporal changes in the image coordinates of the tracked points in multiple views is then presented for detecting the 3-D points which move only rigidly between consecutive frames. These points are used to eliminate the effect of rigid motion of the head and to obtain the changes in the 3-D points and in the corresponding image points due to pure deformation of the face. The methods are thoroughly tested with three multi-image sequences of four cameras including also quite large changes of facial deformations. The test results prove that the proposed affine multi-image least squares matching yields better results than another method using only fixed templates of the previous frame. The elimination of the effect of rigid motion works well and the points where the face is deforming can be correctly detected and the true deformation estimated. A method based on a novel adaptive threshold is also proposed for automated extraction and tracking of circular targets on a moving calibration object.
Bi-planar 2D-to-3D registration in Fourier domain for stereoscopic x-ray motion tracking
NASA Astrophysics Data System (ADS)
Zosso, Dominique; Le Callennec, Benoît; Bach Cuadra, Meritxell; Aminian, Kamiar; Jolles, Brigitte M.; Thiran, Jean-Philippe
2008-03-01
In this paper we present a new method to track bone movements in stereoscopic X-ray image series of the knee joint. The method is based on two different X-ray image sets: a rotational series of acquisitions of the still subject knee that allows the tomographic reconstruction of the three-dimensional volume (model), and a stereoscopic image series of orthogonal projections as the subject performs movements. Tracking the movements of bones throughout the stereoscopic image series means to determine, for each frame, the best pose of every moving element (bone) previously identified in the 3D reconstructed model. The quality of a pose is reflected in the similarity between its theoretical projections and the actual radiographs. We use direct Fourier reconstruction to approximate the three-dimensional volume of the knee joint. Then, to avoid the expensive computation of digitally rendered radiographs (DRR) for pose recovery, we develop a corollary to the 3-dimensional central-slice theorem and reformulate the tracking problem in the Fourier domain. Under the hypothesis of parallel X-ray beams, the heavy 2D-to-3D registration of projections in the signal domain is replaced by efficient slice-to-volume registration in the Fourier domain. Focusing on rotational movements, the translation-relevant phase information can be discarded and we only consider scalar Fourier amplitudes. The core of our motion tracking algorithm can be implemented as a classical frame-wise slice-to-volume registration task. Results on both synthetic and real images confirm the validity of our approach.
NASA Astrophysics Data System (ADS)
Almesallmy, Mohammed
Methodologies are developed for dynamic analysis of mechanical systems with emphasis on inertial propulsion systems. This work adopted the Lagrangian methodology. Lagrangian methodology is the most efficient classical computational technique, which we call Equations of Motion Code (EOMC). The EOMC is applied to several simple dynamic mechanical systems for easier understanding of the method and to aid other investigators in developing equations of motion of any dynamic system. In addition, it is applied to a rigid multibody system, such as Thomson IPS [Thomson 1986]. Furthermore, a simple symbolic algorithm is developed using Maple software, which can be used to convert any nonlinear n-order ordinary differential equation (ODE) systems into 1st-order ODE system in ready format to be used in Matlab software. A side issue, but equally important, we have started corresponding with the U.S. Patent office to persuade them that patent applications, claiming gross linear motion based on inertial propulsion systems should be automatically rejected. The precedent is rejection of patent applications involving perpetual motion machines.
NASA Astrophysics Data System (ADS)
Varkhalev, Iu. P.
1988-07-01
The first Liapunov method is used to investigate sufficient conditions for the existence of asymptotically pendulum motions of a rigid body with a single fixed point. A region of permissible values of the parameters characterizing the body mass corresponding to the asymptotically pendulum motions is defined. The existence of such motions in the case of the Kovalevskii solution and their absence in the Goriachev-Chapolygin solution are demonstrated.
Regular and chaotic motions in applied dynamics of a rigid body.
Beletskii, V. V.; Pivovarov, M. L.; Starostin, E. L.
1996-06-01
Periodic and regular motions, having a predictable functioning mode, play an important role in many problems of dynamics. The achievements of mathematics and mechanics (beginning with Poincare) have made it possible to establish that such motion modes, generally speaking, are local and form "islands" of regularity in a "chaotic sea" of essentially unpredictable trajectories. The development of computer techniques together with theoretical investigations makes it possible to study the global structure of the phase space of many problems having applied significance. A review of a number of such problems, considered by the authors in the past four or five years, is given in this paper. These include orientation and rotation problems of artificial and natural celestial bodies and the problem of controlling the motion of a locomotion robot. The structure of phase space is investigated for these problems. The phase trajectories of the motion are constructed by a numerical implementation of the Poincare point map method. Distinctions are made between regular (or resonance), quasiregular (or conditionally periodic), and chaotic trajectories. The evolution of the phase picture as the parameters are varied is investigated. A large number of "phase portraits" gives a notion of the arrangement and size of the stability islands in the "sea" of chaotic motions, about the appearance and disappearance of these islands as the parameters are varied, etc. (c) 1996 American Institute of Physics. PMID:12780243
NASA Astrophysics Data System (ADS)
Otake, Y.; Murphy, R. J.; Kutzer, M. D.; Taylor, R. H.; Armand, M.
2014-03-01
Background: Snake-like dexterous manipulators may offer significant advantages in minimally-invasive surgery in areas not reachable with conventional tools. Precise control of a wire-driven manipulator is challenging due to factors such as cable deformation, unknown internal (cable friction) and external forces, thus requiring correcting the calibration intraoperatively by determining the actual pose of the manipulator. Method: A method for simultaneously estimating pose and kinematic configuration of a piecewise-rigid object such as a snake-like manipulator from a single x-ray projection is presented. The method parameterizes kinematics using a small number of variables (e.g., 5), and optimizes them simultaneously with the 6 degree-of-freedom pose parameter of the base link using an image similarity between digitally reconstructed radiographs (DRRs) of the manipulator's attenuation model and the real x-ray projection. Result: Simulation studies assumed various geometric magnifications (1.2-2.6) and out-of-plane angulations (0°-90°) in a scenario of hip osteolysis treatment, which demonstrated the median joint angle error was 0.04° (for 2.0 magnification, +/-10° out-of-plane rotation). Average computation time was 57.6 sec with 82,953 function evaluations on a mid-range GPU. The joint angle error remained lower than 0.07° while out-of-plane rotation was 0°-60°. An experiment using video images of a real manipulator demonstrated a similar trend as the simulation study except for slightly larger error around the tip attributed to accumulation of errors induced by deformation around each joint not modeled with a simple pin joint. Conclusions: The proposed approach enables high precision tracking of a piecewise-rigid object (i.e., a series of connected rigid structures) using a single projection image by incorporating prior knowledge about the shape and kinematic behavior of the object (e.g., each rigid structure connected by a pin joint parameterized by a
NASA Astrophysics Data System (ADS)
Monnier, F.; Vallet, B.; Paparoditis, N.; Papelard, J.-P.; David, N.
2013-10-01
This article presents a generic and efficient method to register terrestrial mobile data with imperfect location on a geographic database with better overall accuracy but less details. The registration method proposed in this paper is based on a semi-rigid point to plane ICP ("Iterative Closest Point"). The main applications of such registration is to improve existing geographic databases, particularly in terms of accuracy, level of detail and diversity of represented objects. Other applications include fine geometric modelling and fine façade texturing, object extraction such as trees, poles, road signs marks, facilities, vehicles, etc. The geopositionning system of mobile mapping systems is affected by GPS masks that are only partially corrected by an Inertial Navigation System (INS) which can cause an important drift. As this drift varies non-linearly, but slowly in time, it will be modelled by a translation defined as a piecewise linear function of time which variation over time will be minimized (rigidity term). For each iteration of the ICP, the drift is estimated in order to minimise the distance between laser points and planar model primitives (data attachment term). The method has been tested on real data (a scan of the city of Paris of 3.6 million laser points registered on a 3D model of approximately 71,400 triangles).
ERIC Educational Resources Information Center
Stanford Univ., CA. School Mathematics Study Group.
Transformation geometry topics are covered in one chapter of Unit 11 of this SMSG series. Work with translations, reflections, rotations, and composition of motions is included; vectors are briefly discussed. The chapter on computers and programming deals with recent history and uses of of the computer, organization of a digital computer, an…
Ladstein, Jarle; Evensmoen, Hallvard R.; Håberg, Asta K.; Kristoffersen, Anders; Goa, Pål E.
2016-01-01
Purpose: To compare 2D and 3D echo-planar imaging (EPI) in a higher cognitive level fMRI paradigm. In particular, to study the link between the presence of task-correlated physiological fluctuations and motion and the fMRI contrast estimates from either 2D EPI or 3D EPI datasets, with and without adding nuisance regressors to the model. A signal model in the presence of partly task-correlated fluctuations is derived, and predictions for contrast estimates with and without nuisance regressors are made. Materials and Methods: Thirty-one healthy volunteers were scanned using 2D EPI and 3D EPI during a virtual environmental learning paradigm. In a subgroup of 7 subjects, heart rate and respiration were logged, and the correlation with the paradigm was evaluated. FMRI analysis was performed using models with and without nuisance regressors. Differences in the mean contrast estimates were investigated by analysis-of-variance using Subject, Sequence, Day, and Run as factors. The distributions of group level contrast estimates were compared. Results: Partially task-correlated fluctuations in respiration, heart rate and motion were observed. Statistically significant differences were found in the mean contrast estimates between the 2D EPI and 3D EPI when using a model without nuisance regressors. The inclusion of nuisance regressors for cardiorespiratory effects and motion reduced the difference to a statistically non-significant level. Furthermore, the contrast estimate values shifted more when including nuisance regressors for 3D EPI compared to 2D EPI. Conclusion: The results are consistent with 3D EPI having a higher sensitivity to fluctuations compared to 2D EPI. In the presence partially task-correlated physiological fluctuations or motion, proper correction is necessary to get expectation correct contrast estimates when using 3D EPI. As such task-correlated physiological fluctuations or motion is difficult to avoid in paradigms exploring higher cognitive functions, 2
NASA Technical Reports Server (NTRS)
Longuski, J. M.
1980-01-01
Analytic expressions are found for Euler's Equations of Motion and for the Eulerian Angles for both symmetric and near symmetric rigid bodies under the influence of arbitrary constant body-fixed torques. These solutions provide the body-fixed angular velocities and the attitude of the body, respectively, as functions of time. They are of special interest in applications to spinning spacecraft (such as the Galileo Spacecraft to be launched in 1984) because they include the effect of time-varying spin rate. Thus they can be applied to spin-up and spin-down maneuvers as well as to error analysis for thruster misalignments. The solutions are given for arbitrary initial conditions in terms of Fresnel, Sine and Cosine Integrals. Numerical integration of the governing differential equations has verified that the approximate analytic solutions are very accurate in many physical situations of interest.
Goksel, Orcun; Zahiri-Azar, Reza; Salcudean, Septimiu E
2007-01-01
Motion estimation in sequences of ultrasound echo signals is essential for a wide range of applications. In time domain cross correlation, which is a common motion estimation technique, the displacements are typically not integral multiples of the sampling period. Therefore, to estimate the motion with sub-sample accuracy, 1D and 2D interpolation methods such as parabolic, cosine, and ellipsoid fitting have been introduced in the literature. In this paper, a simulation framework is presented in order to compare the performance of currently available techniques. First, the tissue deformation is modeled using the finite element method (FEM) and then the corresponding pre-/post-deformation radio-frequency (RF) signals are generated using Field II ultrasound simulation software. Using these simulated RF data of deformation, both axial and lateral tissue motion are estimated with sub-sample accuracy. The estimated displacements are then evaluated by comparing them to the known displacements computed by the FEM. This simulation approach was used to evaluate three different lateral motion estimation techniques employing (i) two separate 1D sub-sampling, (ii) two consecutive 1D sub-sampling, and (iii) 2D joint sub-sampling estimators. The estimation errors during two different tissue compression tests are presented with and without spatial filtering. Results show that RF signal processing methods involving tissue deformation can be evaluated using the proposed simulation technique, which employs accurate models. PMID:18002416
Response of a polymer network to the motion of a rigid sphere.
Diamant, Haim
2015-05-01
In view of recent microrheology experiments we re-examine the problem of a rigid sphere oscillating inside a dilute polymer network. The network and its solvent are treated using the two-fluid model. We show that the dynamics of the medium can be decomposed into two independent incompressible flows. The first, dominant at large distances and obeying the Stokes equation, corresponds to the collective flow of the two components as a whole. The other, governing the dynamics over an intermediate range of distances and following the Brinkman equation, describes the flow of the network and solvent relative to one another. The crossover between these two regions occurs at a dynamic length scale which is much larger than the network's mesh size. The analysis focuses on the spatial structure of the medium's response and the role played by the dynamic crossover length. We examine different boundary conditions at the sphere surface. The large-distance collective flow is shown to be independent of boundary conditions and network compressibility, establishing the robustness of two-point microrheology at large separations. The boundary conditions that fit the experimental results for inert spheres in entangled F-actin networks are those of a free network, which does not interact directly with the sphere. Closed-form expressions and scaling relations are derived, allowing for the extraction of material parameters from a combination of one- and two-point microrheology. We discuss a basic deficiency of the two-fluid model and a way to bypass it when analyzing microrheological data. PMID:25957176
Chen, Jun; Zhang, Qing; Liu, Zhi-Fa; Wang, Shuai-Hua; Xiao, Yu; Li, Rong; Xu, Jian-Gang; Zhao, Ya-Ping; Zheng, Fa-Kun; Guo, Guo-Cong
2015-06-01
Two new lead(II) coordination polymers, [Pb(NO3)(tzib)]n (1) and [Pb(tzib)2]n (2), were successfully synthesized from the reaction of a rigid ligand 1-tetrazole-4-imidazole-benzene (Htzib) and lead(II) nitrate in different solvents. The obtained polymers have been characterized by single-crystal X-ray diffraction analyses, which show that both polymers feature 2D layer structures. The inorganic anion nitrate in 1 shows a μ2-κO3:κO3 bridging mode to connect adjacent lead ions into a zigzag chain, and then the organic ligands tzib(-) join the neighboring chains into a 2D layer by a μ3-κN1:κN2:κN6 connection mode. In 2, there are two different bridging modes of the tzib(-) ligand: μ3-κN1:κN2:κN6 and μ3-κN1:κN6 to coordinate the lead ions into a 2D layer structure. Interestingly, both polymers displayed broadband emissions covering the entire visible spectra, which could be tunable to near white-light emission by varying excitation wavelengths. PMID:25952460
Bonanno, Gabriele; Puy, Gilles; Wiaux, Yves; van Heeswijk, Ruud B.; Piccini, Davide; Stuber, Matthias
2014-01-01
Purpose Respiratory motion correction remains a challenge in coronary magnetic resonance imaging (MRI) and current techniques, such as navigator gating, suffer from sub-optimal scan efficiency and ease-of-use. To overcome these limitations, an image-based self-navigation technique is proposed that uses “sub-images” and compressed sensing (CS) to obtain translational motion correction in 2D. The method was preliminarily implemented as a 2D technique and tested for feasibility for targeted coronary imaging. Methods During a 2D segmented radial k-space data acquisition, heavily undersampled sub-images were reconstructed from the readouts collected during each cardiac cycle. These sub-images may then be used for respiratory self-navigation. Alternatively, a CS reconstruction may be used to create these sub-images, so as to partially compensate for the heavy undersampling. Both approaches were quantitatively assessed using simulations and in vivo studies, and the resulting self-navigation strategies were then compared to conventional navigator gating. Results Sub-images reconstructed using CS showed a lower artifact level than sub-images reconstructed without CS. As a result, the final image quality was significantly better when using CS-assisted self-navigation as opposed to the non-CS approach. Moreover, while both self-navigation techniques led to a 69% scan time reduction (as compared to navigator gating), there was no significant difference in image quality between the CS-assisted self-navigation technique and conventional navigator gating, despite the significant decrease in scan time. Conclusions CS-assisted self-navigation using 2D translational motion correction demonstrated feasibility of producing coronary MRA data with image quality comparable to that obtained with conventional navigator gating, and does so without the use of additional acquisitions or motion modeling, while still allowing for 100% scan efficiency and an improved ease-of-use. In
Li, Wei
2016-05-01
Consider cooperative manipulation and transportation of a rigid body by multiple two-wheeled nonholonomic robotic agents that attached to it, the agents are then physically constrained to maintain rigid-formation-motion (RFM); thus the system has two physical motion-constraints at two levels: 1) the nonholonomic constraint at the individual level and 2) the RFM constraint at the system level. First, we provide a novel notion: the encapsulation of a category of control with certain constraints for one motion-mode as a control-law module (CLM), any concrete control law with such constraints is called an instance of the CLM; here two CLMs are provided as the examples. Then we provide an RFM control framework by decomposing a feasible RFM configuration-path as a concatenation of partitions, with one type of CLMs for each partition; thus any instance for each partition can be designed separately and incorporated easily with the interchangeable property, which makes the framework modular, flexible, and adaptive, to satisfy different kinematics requirements. As a result, the transportation is achieved by RFM control of agents. Also, the RFM framework implies a valuable rigid-closure-method for accurate rigid body manipulation even when agents are not attached to the body. PMID:27093718
NASA Astrophysics Data System (ADS)
Rupnik, Ewelina; Jansa, Josef
2013-04-01
Central to our investigation is determination of dynamic behaviour of a highly reflective platform floating on water, as well as derivation of parameters defining instantaneous water state. The employed imaging setup consists of three off-the-shelf dSLR cameras capable of video recording at a 30Hz frame rate. In order to observe a change, the non-rigid and non-diffuse bodies impose the adoption of artificial targetting and custom measurement algorithms. Attention will be given to an in-house software tool implemented to carry out point measurement, correspondence search, tracking and outlier detection methods in the presence of specular reflections and a multimedia scene. A methodology for retrieval of wave parameters in regular wave conditions is also automatically handled by the software and will be discussed. In the context of performed measurements and achieved results, we will point out the extent to which consumer grade camera can fulfil automation and accuracy demands of industrial applications and the pitfalls entailed. Lastly, we will elaborate on visual representation of computed motion and deformations.
NASA Astrophysics Data System (ADS)
Haris, L.; Khotimah, S. N.; Haryanto, F.; Viridi, S.
2014-02-01
Molecular dynamics has been widely used to numerically solve equation of motion of classical many-particle system. It can be used to simulate many systems including biophysics, whose complexity level is determined by the involved elements. Based on this method, a numerical model had been constructed to mimic the behaviour of malaria-infected red blood cells within capillary vessel. The model was governed by three forces namely Coulomb force, normal force, and Stokes force. By utilizing two dimensional four-cells scheme, theoretical observation was carried out to test its capability. Although the parameters were chosen deliberately, all of the quantities were given arbitrary value. Despite this fact, the results were quite satisfactory. Combined with the previous results, it can be said that the proposed model were sufficient enough to mimic the malaria-infected red blood cells motion within obstructed capillary vessel.
Real-time ultrasound-tagging to track the 2D motion of the common carotid artery wall in vivo
Zahnd, Guillaume; Salles, Sébastien; Liebgott, Hervé; Vray, Didier; Sérusclat, André; Moulin, Philippe
2015-02-15
Purpose: Tracking the motion of biological tissues represents an important issue in the field of medical ultrasound imaging. However, the longitudinal component of the motion (i.e., perpendicular to the beam axis) remains more challenging to extract due to the rather coarse resolution cell of ultrasound scanners along this direction. The aim of this study is to introduce a real-time beamforming strategy dedicated to acquire tagged images featuring a distinct pattern in the objective to ease the tracking. Methods: Under the conditions of the Fraunhofer approximation, a specific apodization function was applied to the received raw channel data, in real-time during image acquisition, in order to introduce a periodic oscillations pattern along the longitudinal direction of the radio frequency signal. Analytic signals were then extracted from the tagged images, and subpixel motion tracking of the intima–media complex was subsequently performed offline, by means of a previously introduced bidimensional analytic phase-based estimator. Results: The authors’ framework was applied in vivo on the common carotid artery from 20 young healthy volunteers and 6 elderly patients with high atherosclerosis risk. Cine-loops of tagged images were acquired during three cardiac cycles. Evaluated against reference trajectories manually generated by three experienced analysts, the mean absolute tracking error was 98 ± 84 μm and 55 ± 44 μm in the longitudinal and axial directions, respectively. These errors corresponded to 28% ± 23% and 13% ± 9% of the longitudinal and axial amplitude of the assessed motion, respectively. Conclusions: The proposed framework enables tagged ultrasound images of in vivo tissues to be acquired in real-time. Such unconventional beamforming strategy contributes to improve tracking accuracy and could potentially benefit to the interpretation and diagnosis of biomedical images.
NASA Astrophysics Data System (ADS)
Velioǧlu, Deniz; Cevdet Yalçıner, Ahmet; Zaytsev, Andrey
2016-04-01
Tsunamis are huge waves with long wave periods and wave lengths that can cause great devastation and loss of life when they strike a coast. The interest in experimental and numerical modeling of tsunami propagation and inundation increased considerably after the 2011 Great East Japan earthquake. In this study, two numerical codes, FLOW 3D and NAMI DANCE, that analyze tsunami propagation and inundation patterns are considered. Flow 3D simulates linear and nonlinear propagating surface waves as well as long waves by solving three-dimensional Navier-Stokes (3D-NS) equations. NAMI DANCE uses finite difference computational method to solve 2D depth-averaged linear and nonlinear forms of shallow water equations (NSWE) in long wave problems, specifically tsunamis. In order to validate these two codes and analyze the differences between 3D-NS and 2D depth-averaged NSWE equations, two benchmark problems are applied. One benchmark problem investigates the runup of long waves over a complex 3D beach. The experimental setup is a 1:400 scale model of Monai Valley located on the west coast of Okushiri Island, Japan. Other benchmark problem is discussed in 2015 National Tsunami Hazard Mitigation Program (NTHMP) Annual meeting in Portland, USA. It is a field dataset, recording the Japan 2011 tsunami in Hilo Harbor, Hawaii. The computed water surface elevation and velocity data are compared with the measured data. The comparisons showed that both codes are in fairly good agreement with each other and benchmark data. The differences between 3D-NS and 2D depth-averaged NSWE equations are highlighted. All results are presented with discussions and comparisons. Acknowledgements: Partial support by Japan-Turkey Joint Research Project by JICA on earthquakes and tsunamis in Marmara Region (JICA SATREPS - MarDiM Project), 603839 ASTARTE Project of EU, UDAP-C-12-14 project of AFAD Turkey, 108Y227, 113M556 and 213M534 projects of TUBITAK Turkey, RAPSODI (CONCERT_Dis-021) of CONCERT
NASA Astrophysics Data System (ADS)
Rank, Christopher M.; Heußer, Thorsten; Flach, Barbara; Brehm, Marcus; Kachelrieß, Marc
2015-03-01
We propose a new method for PET/MR respiratory motion compensation, which is based on a 3D-2D registration of strongly undersampled MR data and a) runs in parallel with the PET acquisition, b) can be interlaced with clinical MR sequences, and c) requires less than one minute of the total MR acquisition time per bed position. In our simulation study, we applied a 3D encoded radial stack-of-stars sampling scheme with 160 radial spokes per slice and an acquisition time of 38 s. Gated 4D MR images were reconstructed using a 4D iterative reconstruction algorithm. Based on these images, motion vector fields were estimated using our newly-developed 3D-2D registration framework. A 4D PET volume of a patient with eight hot lesions in the lungs and upper abdomen was simulated and MoCo 4D PET images were reconstructed based on the motion vector fields derived from MR. For evaluation, average SUVmean values of the artificial lesions were determined for a 3D, a gated 4D, a MoCo 4D and a reference (with ten-fold measurement time) gated 4D reconstruction. Compared to the reference, 3D reconstructions yielded an underestimation of SUVmean values due to motion blurring. In contrast, gated 4D reconstructions showed the highest variation of SUVmean due to low statistics. MoCo 4D reconstructions were only slightly affected by these two sources of uncertainty resulting in a significant visual and quantitative improvement in terms of SUVmean values. Whereas temporal resolution was comparable to the gated 4D images, signal-to-noise ratio and contrast-to-noise ratio were close to the 3D reconstructions.
NASA Astrophysics Data System (ADS)
Furtado, H.; Steiner, E.; Stock, M.; Georg, D.; Birkfellner, W.
2014-03-01
Intra-fractional respiratorymotion during radiotherapy is one of themain sources of uncertainty in dose application creating the need to extend themargins of the planning target volume (PTV). Real-time tumormotion tracking by 2D/3D registration using on-board kilo-voltage (kV) imaging can lead to a reduction of the PTV. One limitation of this technique when using one projection image, is the inability to resolve motion along the imaging beam axis. We present a retrospective patient study to investigate the impact of paired portal mega-voltage (MV) and kV images, on registration accuracy. We used data from eighteen patients suffering from non small cell lung cancer undergoing regular treatment at our center. For each patient we acquired a planning CT and sequences of kV and MV images during treatment. Our evaluation consisted of comparing the accuracy of motion tracking in 6 degrees-of-freedom(DOF) using the anterior-posterior (AP) kV sequence or the sequence of kV-MV image pairs. We use graphics processing unit rendering for real-time performance. Motion along cranial-caudal direction could accurately be extracted when using only the kV sequence but in AP direction we obtained large errors. When using kV-MV pairs, the average error was reduced from 3.3 mm to 1.8 mm and the motion along AP was successfully extracted. The mean registration time was of 190+/-35ms. Our evaluation shows that using kVMV image pairs leads to improved motion extraction in 6 DOF. Therefore, this approach is suitable for accurate, real-time tumor motion tracking with a conventional LINAC.
Feng, Yongqiang; Max, Ludo
2014-01-01
Purpose Studying normal or disordered motor control requires accurate motion tracking of the effectors (e.g., orofacial structures). The cost of electromagnetic, optoelectronic, and ultrasound systems is prohibitive for many laboratories, and limits clinical applications. For external movements (lips, jaw), video-based systems may be a viable alternative, provided that they offer high temporal resolution and sub-millimeter accuracy. Method We examined the accuracy and precision of 2D and 3D data recorded with a system that combines consumer-grade digital cameras capturing 60, 120, or 240 frames per second (fps), retro-reflective markers, commercially-available computer software (APAS, Ariel Dynamics), and a custom calibration device. Results Overall mean error (RMSE) across tests was 0.15 mm for static tracking and 0.26 mm for dynamic tracking, with corresponding precision (SD) values of 0.11 and 0.19 mm, respectively. The effect of frame rate varied across conditions, but, generally, accuracy was reduced at 240 fps. The effect of marker size (3 vs. 6 mm diameter) was negligible at all frame rates for both 2D and 3D data. Conclusion Motion tracking with consumer-grade digital cameras and the APAS software can achieve sub-millimeter accuracy at frame rates that are appropriate for kinematic analyses of lip/jaw movements for both research and clinical purposes. PMID:24686484
Feng, Yongqiang; Max, Ludo
2014-04-01
PURPOSE Studying normal or disordered motor control requires accurate motion tracking of the effectors (e.g., orofacial structures). The cost of electromagnetic, optoelectronic, and ultrasound systems is prohibitive for many laboratories and limits clinical applications. For external movements (lips, jaw), video-based systems may be a viable alternative, provided that they offer high temporal resolution and submillimeter accuracy. METHOD The authors examined the accuracy and precision of 2-D and 3-D data recorded with a system that combines consumer-grade digital cameras capturing 60, 120, or 240 frames per second (fps), retro-reflective markers, commercially available computer software (APAS, Ariel Dynamics), and a custom calibration device. RESULTS Overall root-mean-square error (RMSE) across tests was 0.15 mm for static tracking and 0.26 mm for dynamic tracking, with corresponding precision (SD) values of 0.11 and 0.19 mm, respectively. The effect of frame rate varied across conditions, but, generally, accuracy was reduced at 240 fps. The effect of marker size (3- vs. 6-mm diameter) was negligible at all frame rates for both 2-D and 3-D data. CONCLUSION Motion tracking with consumer-grade digital cameras and the APAS software can achieve submillimeter accuracy at frame rates that are appropriate for kinematic analyses of lip/jaw movements for both research and clinical purposes. PMID:24686484
Feng, Bing; Gifford, Howard C.; Beach, Richard D.; Boening, Guido; Gennert, Michael A.; King, Michael A.
2008-01-01
Due to the extended imaging times employed in SPECT and PET, patient motion during imaging is a common clinical occurrence. The fast and accurate correction of the three-dimensional (3D) translational and rotational patient motion in iterative reconstruction is thus necessary to address this important cause of artifacts. We propose a method of incorporating 3D Gaussian interpolation in the projector/backprojector pair to facilitate compensation for rigid-body motion in addition to attenuation and distance-dependent blurring. The method works as the interpolation step for moving the current emission voxel estimates and attenuation maps in the global coordinate system to the new patient location in the rotating coordinate system when calculating the expected projection. It also is employed for moving back the backprojection of the ratio of the measured projection to the expected projection and backprojection of the unit value (sensitivity factor) to the original location. MCAT simulations with known six-degree-of-freedom (6DOF) motion were employed to evaluate the accuracy of our method of motion compensation. We also tested the method with acquisitions of the Data Spectrum Anthropomorphic phantom where motion during SPECT acquisition was measured using the Polaris IR motion tracking system. No motion artifacts were seen on the reconstructions with the motion compensation. PMID:16827485
NASA Astrophysics Data System (ADS)
Connolly, C. M.; Konik, A.; Dasari, P. K. R.; Segars, P.; Zheng, S.; Johnson, K. L.; Dey, J.; King, M. A.
2011-03-01
Patient motion can cause artifacts, which can lead to difficulty in interpretation. The purpose of this study is to create 3D digital anthropomorphic phantoms which model the location of the structures of the chest and upper abdomen of human volunteers undergoing a series of clinically relevant motions. The 3D anatomy is modeled using the XCAT phantom and based on MRI studies. The NURBS surfaces of the XCAT are interactively adapted to fit the MRI studies. A detailed XCAT phantom is first developed from an EKG triggered Navigator acquisition composed of sagittal slices with a 3 x 3 x 3 mm voxel dimension. Rigid body motion states are then acquired at breath-hold as sagittal slices partially covering the thorax, centered on the heart, with 9 mm gaps between them. For non-rigid body motion requiring greater sampling, modified Navigator sequences covering the entire thorax with 3 mm gaps between slices are obtained. The structures of the initial XCAT are then adapted to fit these different motion states. Simultaneous to MRI imaging the positions of multiple reflective markers on stretchy bands about the volunteer's chest and abdomen are optically tracked in 3D via stereo imaging. These phantoms with combined position tracking will be used to investigate both imaging-data-driven and motion-tracking strategies to estimate and correct for patient motion. Our initial application will be to cardiacperfusion SPECT imaging where the XCAT phantoms will be used to create patient activity and attenuation distributions for each volunteer with corresponding motion tracking data from the markers on the body-surface. Monte Carlo methods will then be used to simulate SPECT acquisitions, which will be used to evaluate various motion estimation and correction strategies.
NASA Technical Reports Server (NTRS)
Liu, J. J. F.; Fitzpatrick, P. M.
1973-01-01
Variational equations were applied to the case of a rapidly spinning triaxial body moving in an elliptic orbit, in which the orbital plane is regressing at a constant rate. The explicit differential equations obtained in this application were integrated by the method of averaging to develop secular analytical expressions, which, to first-order in a small parameter, describe the complete space motions of the rigid body under the influence of nonresonant gravity-gradient perturbations. The effects of aerodynamic torque on the rotational motion of an orbiting satellite are studied, as another example of the application of the variational equations derived and the method of averaging.
NASA Technical Reports Server (NTRS)
Macala, G. A.
1983-01-01
A computer program is described that can automatically generate symbolic equations of motion for systems of hinge-connected rigid bodies with tree topologies. The dynamical formulation underlying the program is outlined, and examples are given to show how a symbolic language is used to code the formulation. The program is applied to generate the equations of motion for a four-body model of the Galileo spacecraft. The resulting equations are shown to be a factor of three faster in execution time than conventional numerical subroutines.
Wang, Zhirui; Xu, Jia; Huang, Zuzhen; Zhang, Xudong; Xia, Xiang-Gen; Long, Teng; Bao, Qian
2016-01-01
To detect and estimate ground slowly moving targets in airborne single-channel synthetic aperture radar (SAR), a road-aided ground moving target indication (GMTI) algorithm is proposed in this paper. First, the road area is extracted from a focused SAR image based on radar vision. Second, after stationary clutter suppression in the range-Doppler domain, a moving target is detected and located in the image domain via the watershed method. The target’s position on the road as well as its radial velocity can be determined according to the target’s offset distance and traffic rules. Furthermore, the target’s azimuth velocity is estimated based on the road slope obtained via polynomial fitting. Compared with the traditional algorithms, the proposed method can effectively cope with slowly moving targets partly submerged in a stationary clutter spectrum. In addition, the proposed method can be easily extended to a multi-channel system to further improve the performance of clutter suppression and motion estimation. Finally, the results of numerical experiments are provided to demonstrate the effectiveness of the proposed algorithm. PMID:26999140
Wang, Zhirui; Xu, Jia; Huang, Zuzhen; Zhang, Xudong; Xia, Xiang-Gen; Long, Teng; Bao, Qian
2016-01-01
To detect and estimate ground slowly moving targets in airborne single-channel synthetic aperture radar (SAR), a road-aided ground moving target indication (GMTI) algorithm is proposed in this paper. First, the road area is extracted from a focused SAR image based on radar vision. Second, after stationary clutter suppression in the range-Doppler domain, a moving target is detected and located in the image domain via the watershed method. The target's position on the road as well as its radial velocity can be determined according to the target's offset distance and traffic rules. Furthermore, the target's azimuth velocity is estimated based on the road slope obtained via polynomial fitting. Compared with the traditional algorithms, the proposed method can effectively cope with slowly moving targets partly submerged in a stationary clutter spectrum. In addition, the proposed method can be easily extended to a multi-channel system to further improve the performance of clutter suppression and motion estimation. Finally, the results of numerical experiments are provided to demonstrate the effectiveness of the proposed algorithm. PMID:26999140
Hou, Gary Y.; Provost, Jean; Grondin, Julien; Wang, Shutao; Marquet, Fabrice; Bunting, Ethan; Konofagou, Elisa E.
2015-01-01
Harmonic Motion Imaging for Focused Ultrasound (HMIFU) is a recently developed High-Intensity Focused Ultrasound (HIFU) treatment monitoring method. HMIFU utilizes an Amplitude-Modulated (fAM = 25 Hz) HIFU beam to induce a localized focal oscillatory motion, which is simultaneously estimated and imaged by confocally-aligned imaging transducer. HMIFU feasibilities have been previously shown in silico, in vitro, and in vivo in 1-D or 2-D monitoring of HIFU treatment. The objective of this study is to develop and show the feasibility of a novel fast beamforming algorithm for image reconstruction using GPU-based sparse-matrix operation with real-time feedback. In this study, the algorithm was implemented onto a fully integrated, clinically relevant HMIFU system composed of a 93-element HIFU transducer (fcenter = 4.5MHz) and coaxially-aligned 64-element phased array (fcenter = 2.5MHz) for displacement excitation and motion estimation, respectively. A single transmit beam with divergent beam transmit was used while fast beamforming was implemented using a GPU-based delay-and-sum method and a sparse-matrix operation. Axial HMI displacements were then estimated from the RF signals using a 1-D normalized cross-correlation method and streamed to a graphic user interface. The present work developed and implemented a sparse matrix beamforming onto a fully-integrated, clinically relevant system, which can stream displacement images up to 15 Hz using a GPU-based processing, an increase of 100 fold in rate of streaming displacement images compared to conventional CPU-based conventional beamforming and reconstruction processing. The achieved feedback rate is also currently the fastest and only approach that does not require interrupting the HIFU treatment amongst the acoustic radiation force based HIFU imaging techniques. Results in phantom experiments showed reproducible displacement imaging, and monitoring of twenty two in vitro HIFU treatments using the new 2D system showed a
NASA Astrophysics Data System (ADS)
Frew, D. A.; Scheffer, C.
2008-01-01
Accurate rotordynamic analysis is critical in the achievement of efficient rotary machine design, however the majority of models concern flexible shafts with concentrated supports. The modified Euler equations of motion are used in a numerical model to calculate the natural frequencies and whirl amplitudes of a rigid rotor supported by a single-aerostatic bearing. The bearing is modelled with a non-constant stiffness distribution along its length and a non-symmetric centre of gravity. The results are compared with experimental modal analysis (EMA).
NASA Astrophysics Data System (ADS)
Feng, Tao; Wang, Jizhe; Fung, George; Tsui, Benjamin
2016-01-01
Respiratory motion (RM) and cardiac motion (CM) degrade the quality and resolution in cardiac PET scans. We have developed non-rigid motion estimation methods to estimate both RM and CM based on 4D cardiac gated PET data alone, and compensate the dual respiratory and cardiac (R&C) motions after (MCAR), during (MCDR), and before (MCBR) image reconstruction. In all three R&C motion correction methods, attenuation-activity mismatch effect was modeled by using transformed attenuation maps using the estimated RM. The difference of using activity preserving and non-activity preserving models in R&C correction was also studied. Realistic Monte Carlo simulated 4D cardiac PET data using the 4D XCAT phantom and accurate models of the scanner design parameters and performance characteristics at different noise levels were employed as the known truth and for method development and evaluation. Results from the simulation study suggested that all three dual R&C motion correction methods provide substantial improvement in the quality of 4D cardiac gated PET images as compared with no motion correction. Specifically, the MCDR method yields the best performance for all different noise levels compared with the MCAR and MCBR methods. While MCBR reduces computational time dramatically but the resultant 4D cardiac gated PET images has overall inferior image quality when compared to that from the MCAR and MCDR approaches in the ‘almost’ noise free case. Also, the MCBR method has better noise handling properties when compared with MCAR and provides better quantitative results in high noise cases. When the goal is to reduce scan time or patient radiation dose, MCDR and MCBR provide a good compromise between image quality and computational times.
2D/3D Image Registration using Regression Learning
Chou, Chen-Rui; Frederick, Brandon; Mageras, Gig; Chang, Sha; Pizer, Stephen
2013-01-01
In computer vision and image analysis, image registration between 2D projections and a 3D image that achieves high accuracy and near real-time computation is challenging. In this paper, we propose a novel method that can rapidly detect an object’s 3D rigid motion or deformation from a 2D projection image or a small set thereof. The method is called CLARET (Correction via Limited-Angle Residues in External Beam Therapy) and consists of two stages: registration preceded by shape space and regression learning. In the registration stage, linear operators are used to iteratively estimate the motion/deformation parameters based on the current intensity residue between the target projec-tion(s) and the digitally reconstructed radiograph(s) (DRRs) of the estimated 3D image. The method determines the linear operators via a two-step learning process. First, it builds a low-order parametric model of the image region’s motion/deformation shape space from its prior 3D images. Second, using learning-time samples produced from the 3D images, it formulates the relationships between the model parameters and the co-varying 2D projection intensity residues by multi-scale linear regressions. The calculated multi-scale regression matrices yield the coarse-to-fine linear operators used in estimating the model parameters from the 2D projection intensity residues in the registration. The method’s application to Image-guided Radiation Therapy (IGRT) requires only a few seconds and yields good results in localizing a tumor under rigid motion in the head and neck and under respiratory deformation in the lung, using one treatment-time imaging 2D projection or a small set thereof. PMID:24058278
NASA Astrophysics Data System (ADS)
Arredondo, K.; Billen, M. I.
2015-12-01
Observations of seismicity and seismic tomography provide constraints on the geometry of slabs within mantle, while compression/tension axis derived from moment tensor solutions provide constraints on the internal deformation of slabs. However, since these observations provide only a somewhat blurred or incomplete snapshot of the slab in time, it is difficult to directly relate these observations to the evolution of the slab geometry and the forces acting on and within the slab. In contrast, plate tectonic reconstructions provide time-dependent constraints on the surface motion of plates and the trench at subduction zones, which are related to the dynamical evolution of the slab. We use 2D geodynamical simulations of subduction to explore the relationship between dynamical process within the deforming slab and the observations of surface plate motion and the state-of-stress in slabs. Specifically we utilize models that include the extended Boussinesq approximation (shear heating and latent heat terms in the energy equation), a layered lithosphere with pyrolite, harzburgite and basalt/eclogite, compositionally-dependent phase transitions, and a composite rheology with yielding. The models employ a weak crustal layer that decouples the overriding and subducting plates and allows for dynamically determined trench motion. Here we show that, 1) multiple phase transitions increase slab folding, 2) ridge push significantly increases trench retreat, and 3) strength of the weak crustal layer influences slab detachment. Compared to past studies a more realistic treatment of the phase transitions makes trench retreat more difficult to generate: a weaker plate may encourage slab retreat but detaches once the slab tip crosses into the transition zone due to the rapid increase in slab density. As suggested by previous studies, slab folding within the transition zone changes the direction of forces on the slab and causes periodic changes from trench retreat to trench advance. We
O’Halloran, R; Aksoy, M; Aboussouan, E; Peterson, E; Van, A; Bammer, R
2014-01-01
Purpose Diffusion contrast in diffusion-weighted steady state free precession MRI is generated through the constructive addition of signal from many coherence pathways. Motion-induced phase causes destructive interference which results in loss of signal magnitude and diffusion contrast. In this work, a 3D navigator-based real-time correction of the rigid-body-motion-induced phase errors is developed for diffusion-weighted steady state free precession MRI. Methods The efficacy of the real-time prospective correction method in preserving phase coherence of the steady-state is tested in 3D phantom experiments and 3D scans of healthy human subjects. Results In nearly all experiments, the signal magnitude in images obtained with proposed prospective correction was higher than the signal magnitude in images obtained with no correction. In the human subjects the mean magnitude signal in the data was up to 30 percent higher with prospective motion correction than without. Prospective correction never resulted in a decrease in mean signal magnitude in either the data or in the images. Conclusions The proposed prospective motion correction method is shown to preserve the phase coherence of the steady state in diffusion-weighted steady state free precession MRI, thus mitigating signal magnitude losses that would confound the desired diffusion contrast. PMID:24715414
Feng, Bing; King, Michael A
2006-11-01
We developed a unique method for estimating and compensating rigid-body translations and rotations from scatter and-attenuation-compensated projection data in iterative reconstruction when multiple projection angles are acquired at the same time. During reconstruction, both the non-attenuated and attenuated line-integrals are calculated. Their ratios are then multiplied to the scatter-corrected projection data to estimate scatter-and-attenuation- compensated projection data. At the end of each iteration, the sets of compensated projection data for the angles acquired at the same time are employed to calculate the center-of mass and the inertia tensor, which are used to estimate the location and orientation of the imaging object by the principle-axes method. The estimated motion is applied in the next iteration to reposition the estimated slices and attenuation map in the projector and back-projector to match the pose of the patient at time the projections were acquired. To evaluate our method, we simulated an acquisition of the MCAT phantom with a 3-head SPECT system and imaged the Data Spectrum anthropomorphic phantom on a 3-head IRIX SPECT system. In simulations the phantom translated and rotated by the same amount 9 times. A numerical projector modeling the motion, attenuation, and distance-dependent blurring was used to generate the projection data. Poisson noise was added and 30 noise-realizations were generated. In the experiment with the anthropomorphic phantom, four 360-degree acquisitions were performed with the phantom translated or rotated beforehand. A motion-present dataset was made by mixing the 4 acquisitions. For both the MCAT phantom simulations and anthropomorphic phantom experiment, the motion-present data were reconstructed with 10 iterations of the OSEM which estimates and corrects the motion as described above. Our method obtained visually artifact-free reconstructions, while the reconstruction with no motion correction showed severe artifacts
Vreeburg, J. P. B.
1999-01-22
The paper reports on a current line of research in accelerometry. Two subjects are addressed: the reconstruction of the location and attitude of a linear, or uni-axial, accelerometer from its output under a known motion, and the reconstruction of the acceleration field constituent vectors from the combined output of a known arrangement of linear accelerometers. The arrangement can be arbitrary and, consequently, does not require precision mounting. The component of the acceleration along the sensitive direction gives the ideal output of the accelerometer. When the motion that induces the acceleration is known, a set of five ideal measurement data may suffice to recover the location and attitude of the accelerometer. The formulas for this calculation are given. Their use is illustrated by simulation of an accelerometer and its output. The effects of errors are shown; noisy data are much less detrimental to the reconstruction calculations than systematic errors in the known motion. If the geometry of a set of accelerometers is known, their output can be combined for the reconstruction of the linear and angular motion components that induce the acceleration. Conventionally this is achieved by elimination of the contribution of the angular rate of the geometry to the acceleration field. Only special arrangements of accelerometers, discussed in the literature, allow elimination by elementary operations. A method, thought to be new, is presented for the elimination of the linear and angular acceleration contributions to the field sensed by an arbitrary arrangement of accelerometers, and the consequent recovery of the angular rate vector from the reduced data set. Particular difficulties are encountered in this process but it has been shown that successful reconstruction is possible when a redundant set of data is available. Various options are suggested for further analysis, with the goal to determine the minimum arrangement, identify system errors or improve data
Developing rigid constraint for the estimation of pose and structure from a single image.
Wei, Bao-Gang; Liu, Yong-Huai
2004-07-01
Pose and structure estimation from a single image is a fundamental problem in machine vision and multiple sensor fusion and integration. In this paper we propose using rigid constraints described in different coordinate frames to iteratively estimate structural and camera pose parameters. Using geometric properties of reflected correspondences we put forward a new concept, the reflected pole of a rigid transformation. The reflected pole represents a general analysis of transformations that can be applied to both 2D and 3D transformations. We demonstrate how the concept is applied to calibration by proposing an iterative method to estimate the structural parameters of objects. The method is based on a coarse-to-fine strategy in which initial estimation is obtained through a classical linear algorithm which is then refined by iteration. For a comparative study of performance, we also implemented an extended motion estimation algorithm (from 2D-2D to 3D-2D case) based on epipolar geometry. PMID:15495305
ERIC Educational Resources Information Center
Feng, Yongqiang; Max, Ludo
2014-01-01
Purpose: Studying normal or disordered motor control requires accurate motion tracking of the effectors (e.g., orofacial structures). The cost of electromagnetic, optoelectronic, and ultrasound systems is prohibitive for many laboratories and limits clinical applications. For external movements (lips, jaw), video-based systems may be a viable…
Rigidity of the abdomen ... is a sore area inside the belly or abdomen, the pain will get worse when a hand ... Causes can include: Abscess inside the abdomen Appendicitis ... small intestine, large bowel, or gallbladder ( gastrointestinal ...
Tracking of deformable target in 2D ultrasound images
NASA Astrophysics Data System (ADS)
Royer, Lucas; Marchal, Maud; Le Bras, Anthony; Dardenne, Guillaume; Krupa, Alexandre
2015-03-01
In this paper, we propose a novel approach for automatically tracking deformable target within 2D ultrasound images. Our approach uses only dense information combined with a physically-based model and has therefore the advantage of not using any fiducial marker nor a priori knowledge on the anatomical environment. The physical model is represented by a mass-spring damper system driven by different types of forces where the external forces are obtained by maximizing image similarity metric between a reference target and a deformed target across the time. This deformation is represented by a parametric warping model where the optimal parameters are estimated from the intensity variation. This warping function is well-suited to represent localized deformations in the ultrasound images because it directly links the forces applied on each mass with the motion of all the pixels in its vicinity. The internal forces constrain the deformation to physically plausible motions, and reduce the sensitivity to the speckle noise. The approach was validated on simulated and real data, both for rigid and free-form motions of soft tissues. The results are very promising since the deformable target could be tracked with a good accuracy for both types of motion. Our approach opens novel possibilities for computer-assisted interventions where deformable organs are involved and could be used as a new tool for interactive tracking of soft tissues in ultrasound images.
NASA Astrophysics Data System (ADS)
Gordon, R. G.; Kreemer, C.
2015-12-01
Plate rigidity is the central tenet of plate tectonics. Mounting evidence suggests, however, that significant intraplate deformation occurs in oceanic lithosphere due to horizontal thermal contraction, the rate of which decreases as ≈ 1/age [Kumar & Gordon 2009]. Support for this hypothesis comes from the azimuths of submarine transform faults, which are fit significantly better assuming shrinking plates than by assuming rigid plates [Mishra & Gordon 2015]. Previously we estimated the intraplate velocity field of the Pacific plate accounting for horizontal thermal contraction. The ≈2 mm/yr southeastward motion predicted for the northeastern part of the plate relative to the Pacific-Antarctic Rise may contribute to the non-closure of the Pacific-North America plate motion circuit. In a reference frame in which fix the oldest portion of the Pacific plate, some sites on the plate move up to ≈2 mm/yr [Kreemer & Gordon 2014]. Here we present intraplate velocity fields of the Cocos and Nazca plates and discuss their implications for the non-rigidity of plates and the non-closure of the Pacific-Cocos-Nazca plate circuit, which fails closure by a stunning 14 ±5 mm/yr [DeMets et al. 2010]. If we fix the oldest part of the Cocos plate, intraplate velocities of up to ≈2 mm/yr are estimated, with the fastest motion occurring at the northern end of the plate. If we fix the oldest part of the Nazca plate, displacement rates up to 2 mm/yr are estimated, with the fastest motion occurring in the northeasternmost portion of the plate. In the velocity fields for both plates, the lithosphere adjacent to transform faults along the East Pacific Rise tends to move to the south, which would skew the azimuths of the transform faults clockwise of the values expected for rigid plates, which is the same as the sense of misfit between observed azimuths of transform faults and the azimuths calculated from the MORVEL global set of relative angular velocities [DeMets et al. 2010]. Direct
ERIC Educational Resources Information Center
Gerhart, James B.; Nussbaum, Rudi H.
This monograph was written for the Conference on the New Instructional Materials in Physics held at the University of Washington in summer, 1965. It is intended for use in an introductory course in college physics. It consists of an extensive qualitative discussion of motion followed by a detailed development of the quantitative methods needed to…
ERIC Educational Resources Information Center
Brand, Judith, Ed.
2002-01-01
This issue of Exploratorium Magazine focuses on the topic of motion. Contents include: (1) "First Word" (Zach Tobias); (2) "Cosmic Collisions" (Robert Irion); (3) "The Mobile Cell" (Karen E. Kalumuck); (4) "The Paths of Paths" (Steven Vogel); (5) "Fragments" (Pearl Tesler); (6) "Moving Pictures" (Amy Snyder); (7) "Plants on the Go" (Katharine…
2-D Animation's Not Just for Mickey Mouse.
ERIC Educational Resources Information Center
Weinman, Lynda
1995-01-01
Discusses characteristics of two-dimensional (2-D) animation; highlights include character animation, painting issues, and motion graphics. Sidebars present Silicon Graphics animations tools and 2-D animation programs for the desktop computer. (DGM)
Inertial solvation in femtosecond 2D spectra
NASA Astrophysics Data System (ADS)
Hybl, John; Albrecht Ferro, Allison; Farrow, Darcie; Jonas, David
2001-03-01
We have used 2D Fourier transform spectroscopy to investigate polar solvation. 2D spectroscopy can reveal molecular lineshapes beneath ensemble averaged spectra and freeze molecular motions to give an undistorted picture of the microscopic dynamics of polar solvation. The transition from "inhomogeneous" to "homogeneous" 2D spectra is governed by both vibrational relaxation and solvent motion. Therefore, the time dependence of the 2D spectrum directly reflects the total response of the solvent-solute system. IR144, a cyanine dye with a dipole moment change upon electronic excitation, was used to probe inertial solvation in methanol and propylene carbonate. Since the static Stokes' shift of IR144 in each of these solvents is similar, differences in the 2D spectra result from solvation dynamics. Initial results indicate that the larger propylene carbonate responds more slowly than methanol, but appear to be inconsistent with rotational estimates of the inertial response. To disentangle intra-molecular vibrations from solvent motion, the 2D spectra of IR144 will be compared to the time-dependent 2D spectra of the structurally related nonpolar cyanine dye HDITCP.
Alraiyes, Abdul Hamid; Machuzak, Michael S
2014-12-01
The purpose of this article is to provide an introduction to rigid bronchoscopy (RB). We will briefly discuss its history, evolution, and resurgence while we highlight its versatility and usefulness for today's interventional pulmonologist and thoracic surgeon. Despite being one of the earliest pulmonary procedures described, RB is still an important technique. Advances in thoracic medicine have made this skill critical for a fully functional interventional pulmonary program. If the interventional pulmonologist of this century is to be successful, he or she should be facile in this technique. Despite the availability of RB for decades, the invention of flexible bronchoscopy in 1966 led to a significant downturn in its usage. The growth of the interventional pulmonology field brought RB back into the spot light. Apart from the historic role of RB in treatment of central airway lesions and mechanical debulking of endobronchial lesions, RB is the key instrument that can adapt modern therapeutic tools such as laser, argon plasma coagulation, electrocautery, cryotherapy, and stent deployment. Performing RB requires proper preprocedure preparation, exceptional understanding of upper airway anatomy, specific hand-eye coordination, and open communication between the bronchoscopist and the anesthesiologist. These skills can be primarily learned and maintained with repetition. This article will review information relevant to this technique and lay a foundation to be built upon for years to come. PMID:25463158
Benavides, G.L.; Burt, J.D.
1994-07-12
The invention relates to a clamp mechanism that can be used to attach or temporarily support objects inside of tubular goods. The clamp mechanism can also be modified so that it grips objects. The clamp has a self-centering feature to accommodate out-of-roundness or other internal defections in tubular objects such as pipe. A plurality of clamping shoes are expanded by a linkage which is preferably powered by a motor to contact the inside of a pipe. The motion can be reversed and jaw elements can be connected to the linkage so as to bring the jaws together to grab an object. 12 figs.
Benavides, Gilbert L.; Burt, Jack D.
1994-01-01
The invention relates to a clamp mechanism that can be used to attach or temporarily support objects inside of tubular goods. The clamp mechanism can also be modified so that it grips objects. The clamp has a self-centering feature to accommodate out-of-roundness or other internal defections in tubular objects such as pipe. A plurality of clamping shoes are expanded by a linkage which is preferably powered by a motor to contact the inside of a pipe. The motion can be reversed and jaw elements can be connected to the linkage so as to bring the jaws together to grab an object.
2D-3D Registration of CT Vertebra Volume to Fluoroscopy Projection: A Calibration Model Assessment
NASA Astrophysics Data System (ADS)
Bifulco, P.; Cesarelli, M.; Allen, R.; Romano, M.; Fratini, A.; Pasquariello, G.
2009-12-01
This study extends a previous research concerning intervertebral motion registration by means of 2D dynamic fluoroscopy to obtain a more comprehensive 3D description of vertebral kinematics. The problem of estimating the 3D rigid pose of a CT volume of a vertebra from its 2D X-ray fluoroscopy projection is addressed. 2D-3D registration is obtained maximising a measure of similarity between Digitally Reconstructed Radiographs (obtained from the CT volume) and real fluoroscopic projection. X-ray energy correction was performed. To assess the method a calibration model was realised a sheep dry vertebra was rigidly fixed to a frame of reference including metallic markers. Accurate measurement of 3D orientation was obtained via single-camera calibration of the markers and held as true 3D vertebra position; then, vertebra 3D pose was estimated and results compared. Error analysis revealed accuracy of the order of 0.1 degree for the rotation angles of about 1 mm for displacements parallel to the fluoroscopic plane, and of order of 10 mm for the orthogonal displacement.
Welsh, Mark F; Willing, Ryan T; Giles, Joshua W; Athwal, George S; Johnson, James A
2016-02-29
The purpose of this study was to employ subject-specific computer models to evaluate the interaction of glenohumeral range-of-motion and Hill-Sachs humeral head bone defect size on engagement and shoulder dislocation. We hypothesized that the rate of engagement would increase as defect size increased, and that greater shoulder ROM would engage smaller defects. Three dimensional computer models of 12 shoulders were created. For each shoulder, additional models were created with simulated Hill-Sachs defects of varying severities (XS=15%, S=22.5%, M=30%, L=37.5%, XL=45% and XXL=52.5% of the humeral head diameter, respectively). Rotational motion simulations without translation were conducted. The simulations ended if the defect engaged the anterior glenoid rim with resultant dislocation. The results showed that the rate of engagement was significantly different between defect sizes (0.001
motions are considered. Since engagement of XS and S size Hill-Sachs defects is believed to occur clinically, we suspect that some amount of joint translation may be occurring, causing these defects to engage. Therefore, further studies on clinical pre-operative joint laxity and ROM may enable the prediction of engagement. PMID:26862040
2D/3D Visual Tracker for Rover Mast
NASA Technical Reports Server (NTRS)
Bajracharya, Max; Madison, Richard W.; Nesnas, Issa A.; Bandari, Esfandiar; Kunz, Clayton; Deans, Matt; Bualat, Maria
2006-01-01
A visual-tracker computer program controls an articulated mast on a Mars rover to keep a designated feature (a target) in view while the rover drives toward the target, avoiding obstacles. Several prior visual-tracker programs have been tested on rover platforms; most require very small and well-estimated motion between consecutive image frames a requirement that is not realistic for a rover on rough terrain. The present visual-tracker program is designed to handle large image motions that lead to significant changes in feature geometry and photometry between frames. When a point is selected in one of the images acquired from stereoscopic cameras on the mast, a stereo triangulation algorithm computes a three-dimensional (3D) location for the target. As the rover moves, its body-mounted cameras feed images to a visual-odometry algorithm, which tracks two-dimensional (2D) corner features and computes their old and new 3D locations. The algorithm rejects points, the 3D motions of which are inconsistent with a rigid-world constraint, and then computes the apparent change in the rover pose (i.e., translation and rotation). The mast pan and tilt angles needed to keep the target centered in the field-of-view of the cameras (thereby minimizing the area over which the 2D-tracking algorithm must operate) are computed from the estimated change in the rover pose, the 3D position of the target feature, and a model of kinematics of the mast. If the motion between the consecutive frames is still large (i.e., 3D tracking was unsuccessful), an adaptive view-based matching technique is applied to the new image. This technique uses correlation-based template matching, in which a feature template is scaled by the ratio between the depth in the original template and the depth of pixels in the new image. This is repeated over the entire search window and the best correlation results indicate the appropriate match. The program could be a core for building application programs for systems
Cincu, Rafael; Lorente, Francisco de Asis; Gomez, Joaquin; Eiras, Jose; Agrawal, Amit
2014-01-01
Background: With the advancement of technologies there is more interest in the maintenance of the spine's biomechanical properties focusing on the preservation of the functional motion segment. In present article we describe our experience with 25 cases managed with artificial cervical discs with 28 Solis cage following cervical discectomy with a mean follow-up period of 7.5 year. Materials and Methods: All surgeries were performed by single surgeon from March 2004 to June 2005 with a follow-up till date. Patients with symptomatic single or multiple level diseases that had no prior cervical surgery were candidates for the study. Cohort demographics were comparable. Standardized clinical outcome measures and radiographic examinations were used at prescribed post-operative intervals to compare the treatment groups. Relief in radicular pain, cervical spine motion, and degenerative changes at follow-up were noted. Results: In a total 53 cases, the mean age in prosthesis group was 47 years (age range: 30-63 years) and mean age in cage group was 44 years (32-62 years). Mean hospital stay was 2.7 days in both the groups. At 4 weeks complete cervical movements could be achieved in 19 cases in artificial disc group. Maintenance of movement after 7.5 years was in 76% of these patients. Lordosis was maintained in all cases till date. There was no mortality or wound infection in our series. Conclusions: We conclude that artificial cervical disc could be an alternative to fixed spinal fusion as it represents the most physiological substitute of disc. However, there is need for further studies to support the use of artificial cervical disc prosthesis. PMID:25685218
Energy Science and Technology Software Center (ESTSC)
2005-07-01
Aniso2d is a two-dimensional seismic forward modeling code. The earth is parameterized by an X-Z plane in which the seismic properties Can have monoclinic with x-z plane symmetry. The program uses a user define time-domain wavelet to produce synthetic seismograms anrwhere within the two-dimensional media.
Nonrigid Autofocus Motion Correction for Coronary MR Angiography with a 3D Cones Trajectory
Ingle, R. Reeve; Wu, Holden H.; Addy, Nii Okai; Cheng, Joseph Y.; Yang, Phillip C.; Hu, Bob S.; Nishimura, Dwight G.
2014-01-01
Purpose: To implement a nonrigid autofocus motion correction technique to improve respiratory motion correction of free-breathing whole-heart coronary magnetic resonance angiography (CMRA) acquisitions using an image-navigated 3D cones sequence. Methods: 2D image navigators acquired every heartbeat are used to measure superior-inferior, anterior-posterior, and right-left translation of the heart during a free-breathing CMRA scan using a 3D cones readout trajectory. Various tidal respiratory motion patterns are modeled by independently scaling the three measured displacement trajectories. These scaled motion trajectories are used for 3D translational compensation of the acquired data, and a bank of motion-compensated images is reconstructed. From this bank, a gradient entropy focusing metric is used to generate a nonrigid motion-corrected image on a pixel-by-pixel basis. The performance of the autofocus motion correction technique is compared with rigid-body translational correction and no correction in phantom, volunteer, and patient studies. Results: Nonrigid autofocus motion correction yields improved image quality compared to rigid-body-corrected images and uncorrected images. Quantitative vessel sharpness measurements indicate superiority of the proposed technique in 14 out of 15 coronary segments from three patient and two volunteer studies. Conclusion: The proposed technique corrects nonrigid motion artifacts in free-breathing 3D cones acquisitions, improving image quality compared to rigid-body motion correction. PMID:24006292
Greg Flach, Frank Smith
2011-12-31
Mesh2d is a Fortran90 program designed to generate two-dimensional structured grids of the form [x(i),y(i,j)] where [x,y] are grid coordinates identified by indices (i,j). The x(i) coordinates alone can be used to specify a one-dimensional grid. Because the x-coordinates vary only with the i index, a two-dimensional grid is composed in part of straight vertical lines. However, the nominally horizontal y(i,j0) coordinates along index i are permitted to undulate or otherwise vary. Mesh2d also assigns an integer material type to each grid cell, mtyp(i,j), in a user-specified manner. The complete grid is specified through three separate input files defining the x(i), y(i,j), and mtyp(i,j) variations.
Energy Science and Technology Software Center (ESTSC)
2011-12-31
Mesh2d is a Fortran90 program designed to generate two-dimensional structured grids of the form [x(i),y(i,j)] where [x,y] are grid coordinates identified by indices (i,j). The x(i) coordinates alone can be used to specify a one-dimensional grid. Because the x-coordinates vary only with the i index, a two-dimensional grid is composed in part of straight vertical lines. However, the nominally horizontal y(i,j0) coordinates along index i are permitted to undulate or otherwise vary. Mesh2d also assignsmore » an integer material type to each grid cell, mtyp(i,j), in a user-specified manner. The complete grid is specified through three separate input files defining the x(i), y(i,j), and mtyp(i,j) variations.« less
NASA Astrophysics Data System (ADS)
Lotsch, Bettina V.
2015-07-01
Graphene's legacy has become an integral part of today's condensed matter science and has equipped a whole generation of scientists with an armory of concepts and techniques that open up new perspectives for the postgraphene area. In particular, the judicious combination of 2D building blocks into vertical heterostructures has recently been identified as a promising route to rationally engineer complex multilayer systems and artificial solids with intriguing properties. The present review highlights recent developments in the rapidly emerging field of 2D nanoarchitectonics from a materials chemistry perspective, with a focus on the types of heterostructures available, their assembly strategies, and their emerging properties. This overview is intended to bridge the gap between two major—yet largely disjunct—developments in 2D heterostructures, which are firmly rooted in solid-state chemistry or physics. Although the underlying types of heterostructures differ with respect to their dimensions, layer alignment, and interfacial quality, there is common ground, and future synergies between the various assembly strategies are to be expected.
Rigid particulate matter sensor
Hall, Matthew
2011-02-22
A sensor to detect particulate matter. The sensor includes a first rigid tube, a second rigid tube, a detection surface electrode, and a bias surface electrode. The second rigid tube is mounted substantially parallel to the first rigid tube. The detection surface electrode is disposed on an outer surface of the first rigid tube. The detection surface electrode is disposed to face the second rigid tube. The bias surface electrode is disposed on an outer surface of the second rigid tube. The bias surface electrode is disposed to face the detection surface electrode on the first rigid tube. An air gap exists between the detection surface electrode and the bias surface electrode to allow particulate matter within an exhaust stream to flow between the detection and bias surface electrodes.
Pneumatically erected rigid habitat
NASA Technical Reports Server (NTRS)
Salles, Bradley
1992-01-01
The pneumatically erected rigid habitat concept consists of a structure based on an overexpanded metal bellows. The basic concept incorporates the advantages of both inflatable and rigid structures. The design and erection detail are presented with viewgraphs.
Electron spin-echo techniques for the study of protein motion
NASA Astrophysics Data System (ADS)
Kar, Leela; Johnson, Michael E.; Bowman, Michael K.
Electron spin-echo (ESE) spectroscopy has been used to make the first direct measurements of spin-spin relaxation times of a spin-labeled protein at physiological temperatures. Results from experiments using maleimide-labeled deoxygenated hemoglobin (dHb) from individuals homozygous for sickle cell anemia (dHbS) have been compared with those from control experiments using dHb from normal adults (dHbA). Hb "immobilized" by ammonium sulfate precipitation and by siloxane polymer entrapment have been studied for a suitable "rigid" reference. Two-dimensional ESE (2D-ESE) experiments have been performed using all of these systems. The 2D contour plots show that 2D-ESE is sensitive to the slow motion of dHbS polymers and can differentiate it from both that of immobilized Hb and of HbA molecules in solution at the same temperature and concentration. More importantly, the 2D-ESE technique enables one to select for slower motion and thereby extract the dHbS polymer signal from the total signal generated by the heterogeneous system containing dHbS molecules in solution as well as in the polymer. Computer simulations using current slow motional theories show that detailed motional and structural information may be obtained by such studies. The considerable potential of 2D-ESE spectroscopy in the study of macromolecular motion is illustrated by comparing 2D-ESE with the nonlinear technique of saturation transfer electron paramagnetic resonance.
Lie-Poisson integrators for rigid body dynamics in the solar system
NASA Astrophysics Data System (ADS)
Touma, J.; Wisdom, J.
1994-03-01
The n-body mapping method of Wisdom & Holman (1991) is generalized to encompass rotational dynamics. The Lie-Poisson structure of rigid body dynamics is discussed. Integrators which preserve that structure are derived for the motion of a free rigid body and for the motion of rigid bodies interacting gravitationally with mass points.
Modeling and Control of 2-D Grasping of an Object with Arbitrary Shape under Rolling Contact
NASA Astrophysics Data System (ADS)
Arimoto, Suguru; Yoshida, Morio; Sekimoto, Masahiro; Tahara, Kenji
Modeling, control, and stabilization of dynamics of two-dimensional object grasping by using a pair of multi-joint robot fingers are investigated under rolling contact constraints and an arbitrary geometry of the object and fingertips. First, a fundamental testbed problem of modeling and control of rolling motion between 2-D rigid bodies with an arbitrary shape is treated under the assumption that the two contour curves coincide at the contact point and share the same tangent. The rolling constraint induces the Euler equation of motion that is parameterized by a common arclength parameter and constrained onto the kernel space orthogonally complemented to the image space spanned from the constraint gradient. By extending the analysis to the problem of stable grasp of a 2-D object with an arbitrary shape by a pair of robot fingers, the Euler-Lagrange equation of motion of the overall fingers/object system parametrized by arclength parameters is derived, together with a couple of first-order differential equations that express evolutions of contact points in terms of the second fundamental form. It is shown that 2-D rolling constraints are integrable in the sense of Frobonius even if their Pfaffian forms are characterized by arclength parameters. A control signal called “blind grasping” is introduced and shown to be effective in stabilization of grasping without using the details of the object shape and parameters or external sensing. An extension of the Dirichlet-Lagrange stability theorem to a class of systems with DOF-redundancy under constraints is suggested by using a Morse-Bott-Lyapunov function.
NASA Astrophysics Data System (ADS)
Mettier, Ralph; Pfiffner, O. Adrian
2010-05-01
Surface motion is, apart from the obvious topography, the most easily accessible and best quantifiable characteristic of a typical alpine-style orogen. While it is understood that several different processes, such as i.e. isostatic unloading and thermodynamic effects contribute to the overall motion, it is mostly unclear how large the individual contributions are, and how much of the observed motion is a consequence of ongoing tectonic shortening. A number of methods, such as enhanced GPS measurements, Interferometric Synthetic Aperture Radar (InSAR) and fission track (FT) dating, as well as precise leveling can now provide us with a good description of the vertical motion at present as well as in the fairly recent history of the orogen. This in turn, provides us with reliable, and often much needed, criteria for calibrating conceptual and numerical models of orogenesis and the involved processes. We present a series of finite element models, that attempt to reproduce the observed vertical surface motion on a roughly north-south cross section of the Swiss Alps in the 'ABAQUS' commercial FEM package. Unlike most comparable modeling approaches, we apply a fairly simple formulation of rheology, and focus on a highly complex geometrical representation of the cross section, constructed of individual tectonomorphic units such as the Aar- and Gotthard massifs, the Helvetic and Penninic nappe structures as well as the underlying subduction of the European crust. The models simulate a short timespan, with a fixed rate of shortening prescribed by the boundary conditions and the various interactions between the tectonomorphic units being the dominant adjustable parameters. The resulting motion at the surface of the model, as well as the internal deformation of the individual tectonomorphic units is then examined, interpreted and compared to their real-world counterparts. The models incorporate variations in the chosen physical descriptions of the materials, deforming in
Enhanced rigid-bond restraints
Thorn, Andrea; Dittrich, Birger; Sheldrick, George M.
2012-07-01
An extension is proposed to the rigid-bond description of atomic thermal motion in crystals. The rigid-bond model [Hirshfeld (1976 ▶). Acta Cryst. A32, 239–244] states that the mean-square displacements of two atoms are equal in the direction of the bond joining them. This criterion is widely used for verification (as intended by Hirshfeld) and also as a restraint in structure refinement as suggested by Rollett [Crystallographic Computing (1970 ▶), edited by F. R. Ahmed et al., pp. 167–181. Copenhagen: Munksgaard]. By reformulating this condition, so that the relative motion of the two atoms is required to be perpendicular to the bond, the number of restraints that can be applied per anisotropic atom is increased from about one to about three. Application of this condition to 1,3-distances in addition to the 1,2-distances means that on average just over six restraints can be applied to the six anisotropic displacement parameters of each atom. This concept is tested against very high resolution data of a small peptide and employed as a restraint for protein refinement at more modest resolution (e.g. 1.7 Å)
NASA Astrophysics Data System (ADS)
Najafi, Amin
2014-05-01
Using the Monte Carlo simulations, we have calculated mean-square fluctuations in statistical mechanics, such as those for colloids energy configuration are set on square 2D periodic substrates interacting via a long range screened Coulomb potential on any specific and fixed substrate. Random fluctuations with small deviations from the state of thermodynamic equilibrium arise from the granular structure of them and appear as thermal diffusion with Gaussian distribution structure as well. The variations are showing linear form of the Fluctuation-Dissipation Theorem on the energy of particles constitutive a canonical ensemble with continuous diffusion process of colloidal particle systems. The noise-like variation of the energy per particle and the order parameter versus the Brownian displacement of sum of large number of random steps of particles at low temperatures phase are presenting a markovian process on colloidal particles configuration, too.
NASA Astrophysics Data System (ADS)
Wang, Jin; Ma, Jianyong; Zhou, Changhe
2014-11-01
A 3×3 high divergent 2D-grating with period of 3.842μm at wavelength of 850nm under normal incidence is designed and fabricated in this paper. This high divergent 2D-grating is designed by the vector theory. The Rigorous Coupled Wave Analysis (RCWA) in association with the simulated annealing (SA) is adopted to calculate and optimize this 2D-grating.The properties of this grating are also investigated by the RCWA. The diffraction angles are more than 10 degrees in the whole wavelength band, which are bigger than the traditional 2D-grating. In addition, the small period of grating increases the difficulties of fabrication. So we fabricate the 2D-gratings by direct laser writing (DLW) instead of traditional manufacturing method. Then the method of ICP etching is used to obtain the high divergent 2D-grating.
Curve-based 2D-3D registration of coronary vessels for image guided procedure
NASA Astrophysics Data System (ADS)
Duong, Luc; Liao, Rui; Sundar, Hari; Tailhades, Benoit; Meyer, Andreas; Xu, Chenyang
2009-02-01
3D roadmap provided by pre-operative volumetric data that is aligned with fluoroscopy helps visualization and navigation in Interventional Cardiology (IC), especially when contrast agent-injection used to highlight coronary vessels cannot be systematically used during the whole procedure, or when there is low visibility in fluoroscopy for partially or totally occluded vessels. The main contribution of this work is to register pre-operative volumetric data with intraoperative fluoroscopy for specific vessel(s) occurring during the procedure, even without contrast agent injection, to provide a useful 3D roadmap. In addition, this study incorporates automatic ECG gating for cardiac motion. Respiratory motion is identified by rigid body registration of the vessels. The coronary vessels are first segmented from a multislice computed tomography (MSCT) volume and correspondent vessel segments are identified on a single gated 2D fluoroscopic frame. Registration can be explicitly constrained using one or multiple branches of a contrast-enhanced vessel tree or the outline of guide wire used to navigate during the procedure. Finally, the alignment problem is solved by Iterative Closest Point (ICP) algorithm. To be computationally efficient, a distance transform is computed from the 2D identification of each vessel such that distance is zero on the centerline of the vessel and increases away from the centerline. Quantitative results were obtained by comparing the registration of random poses and a ground truth alignment for 5 datasets. We conclude that the proposed method is promising for accurate 2D-3D registration, even for difficult cases of occluded vessel without injection of contrast agent.
Ambiguous fluidity and rigidity and diamonds that ooze!
Meyer, G E; Dougherty, T J
1990-01-01
If white hemicircles rotate over the edges of a black diamond, there occurs an ambiguity of rigidity and motion. As the hemicircles obscure the vertices of the diamond, the figure transforms from a diamond to a rotating, nonrigid cross made of a tar-like fluid. When the corners reappear, the stimulus again becomes a rigid, solid diamond. Visibility of the vertices implies rigidity. If white squares are rotated, fluidity is not perceived. If the diamond has sawtooth edges and the hemicircles are rotated, no fluidity is perceived. Similarly, if illusory contours suggest the amodal completion of the vertices, rigidity is maintained. PMID:2096367
Local Metric Learning in 2D/3D Deformable Registration With Application in the Abdomen
Chou, Chen-Rui; Mageras, Gig; Pizer, Stephen
2015-01-01
In image-guided radiotherapy (IGRT) of disease sites subject to respiratory motion, soft tissue deformations can affect localization accuracy. We describe the application of a method of 2D/3D deformable registration to soft tissue localization in abdomen. The method, called registration efficiency and accuracy through learning a metric on shape (REALMS), is designed to support real-time IGRT. In a previously developed version of REALMS, the method interpolated 3D deformation parameters for any credible deformation in a deformation space using a single globally-trained Riemannian metric for each parameter. We propose a refinement of the method in which the metric is trained over a particular region of the deformation space, such that interpolation accuracy within that region is improved. We report on the application of the proposed algorithm to IGRT in abdominal disease sites, which is more challenging than in lung because of low intensity contrast and nonrespiratory deformation. We introduce a rigid translation vector to compensate for nonrespiratory deformation, and design a special region-of-interest around fiducial markers implanted near the tumor to produce a more reliable registration. Both synthetic data and actual data tests on abdominal datasets show that the localized approach achieves more accurate 2D/3D deformable registration than the global approach. PMID:24771575
Efficient implementation of the rank correlation merit function for 2D/3D registration.
Figl, M; Bloch, C; Gendrin, C; Weber, C; Pawiro, S A; Hummel, J; Markelj, P; Pernus, F; Bergmann, H; Birkfellner, W
2010-10-01
A growing number of clinical applications using 2D/3D registration have been presented recently. Usually, a digitally reconstructed radiograph is compared iteratively to an x-ray image of the known projection geometry until a match is achieved, thus providing six degrees of freedom of rigid motion which can be used for patient setup in image-guided radiation therapy or computer-assisted interventions. Recently, stochastic rank correlation, a merit function based on Spearman's rank correlation coefficient, was presented as a merit function especially suitable for 2D/3D registration. The advantage of this measure is its robustness against variations in image histogram content and its wide convergence range. The considerable computational expense of computing an ordered rank list is avoided here by comparing randomly chosen subsets of the DRR and reference x-ray. In this work, we show that it is possible to omit the sorting step and to compute the rank correlation coefficient of the full image content as fast as conventional merit functions. Our evaluation of a well-calibrated cadaver phantom also confirms that rank correlation-type merit functions give the most accurate results if large differences in the histogram content for the DRR and the x-ray image are present. PMID:20844334
Quantum mechanics of a generalised rigid body
NASA Astrophysics Data System (ADS)
Gripaios, Ben; Sutherland, Dave
2016-05-01
We consider the quantum version of Arnold’s generalisation of a rigid body in classical mechanics. Thus, we quantise the motion on an arbitrary Lie group manifold of a particle whose classical trajectories correspond to the geodesics of any one-sided-invariant metric. We show how the derivation of the spectrum of energy eigenstates can be simplified by making use of automorphisms of the Lie algebra and (for groups of type I) by methods of harmonic analysis. We show how the method can be extended to cosets, generalising the linear rigid rotor. As examples, we consider all connected and simply connected Lie groups up to dimension 3. This includes the universal cover of the archetypical rigid body, along with a number of new exactly solvable models. We also discuss a possible application to the topical problem of quantising a perfect fluid.
Energy Science and Technology Software Center (ESTSC)
2004-08-01
AnisWave2D is a 2D finite-difference code for a simulating seismic wave propagation in fully anisotropic materials. The code is implemented to run in parallel over multiple processors and is fully portable. A mesh refinement algorithm has been utilized to allow the grid-spacing to be tailored to the velocity model, avoiding the over-sampling of high-velocity materials that usually occurs in fixed-grid schemes.
Bao, Zhaosheng; Hong, Jeong-Mo; Teran, Joseph; Fedkiw, Ronald
2007-01-01
We propose a novel approach to fracturing (and denting) brittle materials. To avoid the computational burden imposed by the stringent time step restrictions of explicit methods or with solving nonlinear systems of equations for implicit methods, we treat the material as a fully rigid body in the limit of infinite stiffness. In addition to a triangulated surface mesh and level set volume for collisions, each rigid body is outfitted with a tetrahedral mesh upon which finite element analysis can be carried out to provide a stress map for fracture criteria. We demonstrate that the commonly used stress criteria can lead to arbitrary fracture (especially for stiff materials) and instead propose the notion of a time averaged stress directly into the FEM analysis. When objects fracture, the virtual node algorithm provides new triangle and tetrahedral meshes in a straightforward and robust fashion. Although each new rigid body can be rasterized to obtain a new level set, small shards can be difficult to accurately resolve. Therefore, we propose a novel collision handling technique for treating both rigid bodies and rigid body thin shells represented by only a triangle mesh. PMID:17218752
Method to estimate center of rigidity using vibration recordings
Safak, Erdal; Celebi, Mehmet
1990-01-01
A method to estimate the center of rigidity of buildings by using vibration recordings is presented. The method is based on the criterion that the coherence of translational motions with the rotational motion is minimum at the center of rigidity. Since the coherence is a function of frequency, a gross but frequency-independent measure of the coherency is defined as the integral of the coherence function over the frequency. The center of rigidity is determined by minimizing this integral. The formulation is given for two-dimensional motions. Two examples are presented for the method; a rectangular building with ambient-vibration recordings, and a triangular building with earthquake-vibration recordings. Although the examples given are for buildings, the method can be applied to any structure with two-dimensional motions.
Hallquist, J.O.
1982-02-01
This revised report provides an updated user's manual for DYNA2D, an explicit two-dimensional axisymmetric and plane strain finite element code for analyzing the large deformation dynamic and hydrodynamic response of inelastic solids. A contact-impact algorithm permits gaps and sliding along material interfaces. By a specialization of this algorithm, such interfaces can be rigidly tied to admit variable zoning without the need of transition regions. Spatial discretization is achieved by the use of 4-node solid elements, and the equations-of motion are integrated by the central difference method. An interactive rezoner eliminates the need to terminate the calculation when the mesh becomes too distorted. Rather, the mesh can be rezoned and the calculation continued. The command structure for the rezoner is described and illustrated by an example.
Pal, Tanmoy; Bhattacharjee, Somendra M
2016-05-01
The temperature dependence of DNA flexibility is studied in the presence of stretching and unzipping forces. Two classes of models are considered. In one case the origin of elasticity is entropic due to the polymeric correlations, and in the other the double-stranded DNA is taken to have an intrinsic rigidity for bending. In both cases single strands are completely flexible. The change in the elastic constant for the flexible case due to thermally generated bubbles is obtained exactly. For the case of intrinsic rigidity, the elastic constant is found to be proportional to the square root of the bubble number fluctuation. PMID:27300825
NASA Astrophysics Data System (ADS)
Pal, Tanmoy; Bhattacharjee, Somendra M.
2016-05-01
The temperature dependence of DNA flexibility is studied in the presence of stretching and unzipping forces. Two classes of models are considered. In one case the origin of elasticity is entropic due to the polymeric correlations, and in the other the double-stranded DNA is taken to have an intrinsic rigidity for bending. In both cases single strands are completely flexible. The change in the elastic constant for the flexible case due to thermally generated bubbles is obtained exactly. For the case of intrinsic rigidity, the elastic constant is found to be proportional to the square root of the bubble number fluctuation.
Self-propulsion of a body with rigid surface and variable coefficient of lift in a perfect fluid
NASA Astrophysics Data System (ADS)
Ramodanov, Sergey M.; Tenenev, Valentin A.; Treschev, Dmitry V.
2012-11-01
We study the system of a 2D rigid body moving in an unbounded volume of incompressible, vortex-free perfect fluid which is at rest at infinity. The body is equipped with a gyrostat and a so-called Flettner rotor. Due to the latter the body is subject to a lifting force (Magnus effect). The rotational velocities of the gyrostat and the rotor are assumed to be known functions of time (control inputs). The equations of motion are presented in the form of the Kirchhoff equations. The integrals of motion are given in the case of piecewise continuous control. Using these integrals we obtain a (reduced) system of first-order differential equations on the configuration space. Then an optimal control problem for several types of the inputs is solved using genetic algorithms.
NASA Astrophysics Data System (ADS)
Baka, N.; Lelieveldt, B. P. F.; Schultz, C.; Niessen, W.; van Walsum, T.
2015-05-01
During percutaneous coronary interventions (PCI) catheters and arteries are visualized by x-ray angiography (XA) sequences, using brief contrast injections to show the coronary arteries. If we could continue visualizing the coronary arteries after the contrast agent passed (thus in non-contrast XA frames), we could potentially lower contrast use, which is advantageous due to the toxicity of the contrast agent. This paper explores the possibility of such visualization in mono-plane XA acquisitions with a special focus on respiratory based coronary artery motion estimation. We use the patient specific coronary artery centerlines from pre-interventional 3D CTA images to project on the XA sequence for artery visualization. To achieve this, a framework for registering the 3D centerlines with the mono-plane 2D + time XA sequences is presented. During the registration the patient specific cardiac and respiratory motion is learned. We investigate several respiratory motion estimation strategies with respect to accuracy, plausibility and ease of use for motion prediction in XA frames with and without contrast. The investigated strategies include diaphragm motion based prediction, and respiratory motion extraction from the guiding catheter tip motion. We furthermore compare translational and rigid respiratory based heart motion. We validated the accuracy of the 2D/3D registration and the respiratory and cardiac motion estimations on XA sequences of 12 interventions. The diaphragm based motion model and the catheter tip derived motion achieved 1.58 mm and 1.83 mm median 2D accuracy, respectively. On a subset of four interventions we evaluated the artery visualization accuracy for non-contrast cases. Both diaphragm, and catheter tip based prediction performed similarly, with about half of the cases providing satisfactory accuracy (median error < 2 mm).
Effective rigidity of membranes
NASA Astrophysics Data System (ADS)
Peliti, L.
1986-12-01
The role of thermal fluctuations of shape (undulations) in reducing the effective rigidity of membranes is reviewed. The consequences of this effect on vesicle size distribution and on the structure of microemulsions, as well as on other physical phenomena, are sketched.
Steckle, W.P. Jr.; Mitchell, M.A.; Aspen, P.G.
1998-12-31
This is the final report of a three-year, Laboratory Directed Research and Development (LDRD) project at the Los Alamos National Laboratory (LANL). Organic analogues to inorganic zeolites would be a significant step forward in engineered porous materials and would provide advantages in range, selectivity, tailorability, and processing. Rigid molecular foams or {open_quotes}organic zeolites{close_quotes} would not be crystalline materials and could be tailored over a broader range of pore sizes and volumes. A novel process for preparing hypercrosslinked polymeric foams has been developed via a Friedel-Crafts polycondensation reaction. A series of rigid hypercrosslinked foams have been prepared using simple rigid polyaromatic hydrocarbons including benzene, biphenyl, m-terphenyl, diphenylmethane, and polystyrene, with dichloroxylene (DCX) as the pore size. After drying the foams are robust and rigid. Densities of the resulting foams can range from 0.15 g/cc to 0.75 g/cc. Nitrogen adsorption studies have shown that by judiciously selecting monomers and the crosslinking agent along with the level of crosslinking and the cure time of the resulting gel, the pore size, pore size distribution, and the total surface area of the foam can be tailored. Surface areas range from 160 to 1,200 m{sup 2}/g with pore sizes ranging from 6 {angstrom} to 2,000 {angstrom}.
Electrostatics of Rigid Polyelectrolytes
Wong, G.C.L.
2009-06-04
The organization of rigid biological polyelectrolytes by multivalent ions and macroions are important for many fundamental problems in biology and biomedicine, such as cytoskeletal regulation and antimicrobial sequestration in cystic fibrosis. These polyelectrolytes have been used as model systems for understanding electrostatics in complex fluids. Here, we review some recent results in theory, simulations, and experiments.
NASA Astrophysics Data System (ADS)
Mayor, Louise
2016-05-01
Graphene might be the most famous example, but there are other 2D materials and compounds too. Louise Mayor explains how these atomically thin sheets can be layered together to create flexible “van der Waals heterostructures”, which could lead to a range of novel applications.
NASA Astrophysics Data System (ADS)
Hartschuh, R. D.; Wargacki, S. P.; Xiong, H.; Neiswinger, J.; Kisliuk, A.; Sihn, S.; Ward, V.; Vaia, R. A.; Sokolov, A. P.
2008-08-01
Viruses have traditionally been studied as pathogens, but in recent years they have been adapted for applications ranging from drug delivery and gene therapy to nanotechnology, photonics, and electronics. Although the structures of many viruses are known, most of their biophysical properties remain largely unexplored. Using Brillouin light scattering, we analyzed the mechanical rigidity, intervirion coupling, and vibrational eigenmodes of Wiseana iridovirus (WIV). We identified phonon modes propagating through the viral assemblies as well as the localized vibrational eigenmode of individual viruses. The measurements indicate a Young’s modulus of ˜7GPa for single virus particles and their assemblies, surprisingly high for “soft” materials. Mechanical modeling confirms that the DNA core dominates the WIV rigidity. The results also indicate a peculiar mechanical coupling during self-assembly of WIV particles.
Motion and deformation compensation for freehand prostate biopsies
NASA Astrophysics Data System (ADS)
Khallaghi, Siavash; Nouranian, Saman; Sojoudi, Samira; Ashab, Hussam A.; Machan, Lindsay; Chang, Silvia; Black, Peter; Gleave, Martin; Goldenberg, Larry; Abolmaesumi, Purang
2014-03-01
In this paper, we present a registration pipeline to compensate for prostate motion and deformation during targeted freehand prostate biopsies. We perform 2D-3D registration by reconstructing a thin-volume around the real-time 2D ultrasound imaging plane. Constrained Sum of Squared Differences (SSD) and gradient descent optimization are used to rigidly align the moving volume to the fixed thin-volume. Subsequently, B-spline de- formable registration is performed to compensate for remaining non-linear deformations. SSD and zero-bounded Limited memory Broyden Fletcher Goldfarb Shannon (LBFGS) optimizer are used to find the optimum B-spline parameters. Registration results are validated on five prostate biopsy patients. Initial experiments suggest thin- volume-to-volume registration to be more effective than slice-to-volume registration. Also, a minimum consistent 2 mm improvement of Target Registration Error (TRE) is achieved following the deformable registration.
Silt motion simulation using finite volume particle method
NASA Astrophysics Data System (ADS)
Jahanbakhsh, E.; Vessaz, C.; Avellan, F.
2014-03-01
In this paper, we present a 3-D FVPM which features rectangular top-hat kernels. With this method, interaction vectors are computed exactly and efficiently. We introduce a new method to enforce the no-slip boundary condition. With this boundary enforcement, the interaction forces between fluid and wall are computed accurately. We employ the boundary force to predict the motion of rigid spherical silt particles inside the fluid. To validate the model, we simulate the 2-D sedimentation of a single particle in viscous fluid tank and compare results with benchmark data. The particle resolution is verified by convergence study. We also simulate the sedimentation of two particles exhibiting drafting, kissing and tumbling phenomena in 2-D and 3-D. We compare the results with other numerical solutions.
Kitazaki, M; Shimojo, S
1996-01-01
The generic-view principle (GVP) states that given a 2-D image the visual system interprets it as a generic view of a 3-D scene when possible. The GVP was applied to 3-D-motion perception to show how the visual system decomposes retinal image motion into three components of 3-D motion: stretch/shrinkage, rotation, and translation. First, the optical process of retinal image motion was analyzed, and predictions were made based on the GVP in the inverse-optical process. Then experiments were conducted in which the subject judged perception of stretch/shrinkage, rotation in depth, and translation in depth for a moving bar stimulus. Retinal-image parameters-2-D stretch/shrinkage, 2-D rotation, and 2-D translation-were manipulated categorically and exhaustively. The results were highly consistent with the predictions. The GVP seems to offer a broad and general framework for understanding the ambiguity-solving process in motion perception. Its relationship to other constraints such as that of rigidity is discussed. PMID:8923550
NASA Technical Reports Server (NTRS)
Gasch, Matthew J.
2011-01-01
NASA Exploration Systems Mission Directorate s (ESMD) Entry, Descent, and Landing (EDL) Technology Development Project (TDP) and the NASA Aeronautics Research Mission Directorate s (ARMD) Hypersonics Project are developing new advanced rigid ablators in an effort to substantially increase reliability, decrease mass, and reduce life cycle cost of rigid aeroshell-based entry systems for multiple missions. Advanced Rigid Ablators combine ablation resistant top layers capable of high heat flux entry and enable high-speed EDL with insulating mass-efficient bottom that, insulate the structure and lower the areal weight. These materials may benefit Commercial Orbital Transportation Services (COTS) vendors and may potentially enable new NASA missions for higher velocity returns (e.g. asteroid, Mars). The materials have been thermally tested to 400-450 W/sq cm at the Laser Hardened Materials Evaluation Lab (LHMEL), Hypersonics Materials Evaluation Test System (HyMETS) and in arcjet facilities. Tested materials exhibit much lower backface temperatures and reduced recession over the baseline materials (PICA). Although the EDL project is ending in FY11, NASA in-house development of advanced ablators will continue with a focus on varying resin systems and fiber/resin interactions.
Sim, Jaehyun; Sim, Jun; Park, Eunsung; Lee, Julian
2015-06-01
Many proteins undergo large-scale motions where relatively rigid domains move against each other. The identification of rigid domains, as well as the hinge residues important for their relative movements, is important for various applications including flexible docking simulations. In this work, we develop a method for protein rigid domain identification based on an exhaustive enumeration of maximal rigid domains, the rigid domains not fully contained within other domains. The computation is performed by mapping the problem to that of finding maximal cliques in a graph. A minimal set of rigid domains are then selected, which cover most of the protein with minimal overlap. In contrast to the results of existing methods that partition a protein into non-overlapping domains using approximate algorithms, the rigid domains obtained from exact enumeration naturally contain overlapping regions, which correspond to the hinges of the inter-domain bending motion. The performance of the algorithm is demonstrated on several proteins. PMID:25820699
Energy Science and Technology Software Center (ESTSC)
2001-01-31
This software reduces the data from two-dimensional kSA MOS program, k-Space Associates, Ann Arbor, MI. Initial MOS data is recorded without headers in 38 columns, with one row of data per acquisition per lase beam tracked. The final MOSS 2d data file is reduced, graphed, and saved in a tab-delimited column format with headers that can be plotted in any graphing software.
NASA Astrophysics Data System (ADS)
D'Avino, G.; Maffettone, P. L.; Hulsen, M. A.; Peters, G. W. M.
2007-09-01
In this work a new numerical method for concentrated inertialess rigid particle suspensions in a planar elongational flow using a fixed mesh is presented. The main concept is to randomly relocate a particle on an inflow section of the domain when it crosses the outflow boundaries. A three-layer domain is considered in order to: (i) develop a small computational domain as the representative sample of the whole suspension, (ii) impose the elongational flow boundary conditions far from the particles, (iii) achieve a steady state (in a statistical meaning). Our scheme uses a time-independent fixed grid avoiding the difficulties involved in deforming meshes and remeshing of the domain. In this way, computations can proceed indefinitely and micro-structural fluctuations around a steady state can be studied. A fictitious domain is implemented in order to easily manage the rigid-body motion. The particles are described by their boundaries only (rigid-ring description) and the rigid-body motion is imposed through Lagrange multipliers. The bulk properties are recovered by using an averaging procedure where the traction forces on the particle surface are recovered by the Lagrange multipliers. The scheme has been combined with a standard velocity-pressure finite element formulation and 2D simulations of a large number (150 and 225) of particles in a Newtonian medium are performed. Local as well as bulk properties are evaluated and discussed. The results show very good agreement with dilute theory as well as with other numerical simulations in the literature for higher concentrations. Our formulation is well suited for viscoelastic suspensions and can be easily extended to 3D simulations.
NASA Astrophysics Data System (ADS)
Pan, Bing; Yu, Liping; Wu, Dafang
2014-02-01
The ideal pinhole imaging model commonly assumed for an ordinary two-dimensional digital image correlation (2D-DIC) system is neither perfect nor stable because of the existence of small out-of-plane motion of the test sample surface that occurred after loading, small out-of-plane motion of the sensor target due to temperature variation of a camera and unavoidable geometric distortion of an imaging lens. In certain cases, these disadvantages can lead to significant errors in the measured displacements and strains. Although a high-quality bilateral telecentric lens has been strongly recommended to be used in the 2D-DIC system as an essential optical component to achieve high-accuracy measurement, it is not generally applicable due to its fixed field of view, limited depth of focus and high cost. To minimize the errors associated with the imperfectness and instability of a common 2D-DIC system using a low-cost imaging lens, a generalized compensation method using a non-deformable reference sample is proposed in this work. With the proposed method, the displacement of the reference sample rigidly attached behind the test sample is first measured using 2D-DIC, and then it is fitted using a parametric model. The fitted parametric model is then used to correct the displacements of the deformed sample to remove the influences of these unfavorable factors. The validity of the proposed compensation method is first verified using out-of-plane translation, out-of-plane rotation, in-plane translation tests and their combinations. Uniaxial tensile tests of an aluminum specimen were also performed to quantitatively examine the strain accuracy of the proposed compensation method. Experiments show that the proposed compensation method is an easy-to-implement yet effective technique for achieving high-accuracy deformation measurement using an ordinary 2D-DIC system.
Nanoimprint lithography: 2D or not 2D? A review
NASA Astrophysics Data System (ADS)
Schift, Helmut
2015-11-01
Nanoimprint lithography (NIL) is more than a planar high-end technology for the patterning of wafer-like substrates. It is essentially a 3D process, because it replicates various stamp topographies by 3D displacement of material and takes advantage of the bending of stamps while the mold cavities are filled. But at the same time, it keeps all assets of a 2D technique being able to pattern thin masking layers like in photon- and electron-based traditional lithography. This review reports about 20 years of development of replication techniques at Paul Scherrer Institut, with a focus on 3D aspects of molding, which enable NIL to stay 2D, but at the same time enable 3D applications which are "more than Moore." As an example, the manufacturing of a demonstrator for backlighting applications based on thermally activated selective topography equilibration will be presented. This technique allows generating almost arbitrary sloped, convex and concave profiles in the same polymer film with dimensions in micro- and nanometer scale.
Finite Element Analysis of 2-D Elastic Contacts Involving FGMs
NASA Astrophysics Data System (ADS)
Abhilash, M. N.; Murthy, H.
2014-05-01
The response of elastic indenters in contact with Functionally Graded Material (FGM) coated homogeneous elastic half space has been presented in the current paper. Finite element analysis has been used due to its ability to handle complex geometry, material, and boundary conditions. Indenters of different typical surface profiles have been considered and the problem has been idealized as a two-dimensional (2D) plane strain problem considering only normal loads. Initially, indenters were considered to be rigid and the results were validated with the solutions presented in the literature. The analysis has then been extended to the case of elastic indenters on FGM-coated half spaces and the results are discussed.
Chiang, Ta-Kuan; Straub, Douglas L.; Dennis, Richard A.
2000-01-01
The present invention involves a porous rigid filter including a plurality of concentric filtration elements having internal flow passages and forming external flow passages there between. The present invention also involves a pressure vessel containing the filter for the removal of particulates from high pressure particulate containing gases, and further involves a method for using the filter to remove such particulates. The present filter has the advantage of requiring fewer filter elements due to the high surface area-to-volume ratio provided by the filter, requires a reduced pressure vessel size, and exhibits enhanced mechanical design properties, improved cleaning properties, configuration options, modularity and ease of fabrication.
FPCAS2D user's guide, version 1.0
NASA Astrophysics Data System (ADS)
Bakhle, Milind A.
1994-12-01
The FPCAS2D computer code has been developed for aeroelastic stability analysis of bladed disks such as those in fans, compressors, turbines, propellers, or propfans. The aerodynamic analysis used in this code is based on the unsteady two-dimensional full potential equation which is solved for a cascade of blades. The structural analysis is based on a two degree-of-freedom rigid typical section model for each blade. Detailed explanations of the aerodynamic analysis, the numerical algorithms, and the aeroelastic analysis are not given in this report. This guide can be used to assist in the preparation of the input data required by the FPCAS2D code. A complete description of the input data is provided in this report. In addition, four test cases, including inputs and outputs, are provided.
3D Reconstruction of Human Motion from Monocular Image Sequences.
Wandt, Bastian; Ackermann, Hanno; Rosenhahn, Bodo
2016-08-01
This article tackles the problem of estimating non-rigid human 3D shape and motion from image sequences taken by uncalibrated cameras. Similar to other state-of-the-art solutions we factorize 2D observations in camera parameters, base poses and mixing coefficients. Existing methods require sufficient camera motion during the sequence to achieve a correct 3D reconstruction. To obtain convincing 3D reconstructions from arbitrary camera motion, our method is based on a-priorly trained base poses. We show that strong periodic assumptions on the coefficients can be used to define an efficient and accurate algorithm for estimating periodic motion such as walking patterns. For the extension to non-periodic motion we propose a novel regularization term based on temporal bone length constancy. In contrast to other works, the proposed method does not use a predefined skeleton or anthropometric constraints and can handle arbitrary camera motion. We achieve convincing 3D reconstructions, even under the influence of noise and occlusions. Multiple experiments based on a 3D error metric demonstrate the stability of the proposed method. Compared to other state-of-the-art methods our algorithm shows a significant improvement. PMID:27093439
Prospective Motion Correction using Inductively-Coupled Wireless RF Coils
Ooi, Melvyn B.; Aksoy, Murat; Maclaren, Julian; Watkins, Ronald D.; Bammer, Roland
2013-01-01
Purpose A novel prospective motion correction technique for brain MRI is presented that uses miniature wireless radio-frequency (RF) coils, or “wireless markers”, for position tracking. Methods Each marker is free of traditional cable connections to the scanner. Instead, its signal is wirelessly linked to the MR receiver via inductive coupling with the head coil. Real-time tracking of rigid head motion is performed using a pair of glasses integrated with three wireless markers. A tracking pulse-sequence, combined with knowledge of the markers’ unique geometrical arrangement, is used to measure their positions. Tracking data from the glasses is then used to prospectively update the orientation and position of the image-volume so that it follows the motion of the head. Results Wireless-marker position measurements were comparable to measurements using traditional wired RF tracking coils, with the standard deviation of the difference < 0.01 mm over the range of positions measured inside the head coil. RF safety was verified with B1 maps and temperature measurements. Prospective motion correction was demonstrated in a 2D spin-echo scan while the subject performed a series of deliberate head rotations. Conclusion Prospective motion correction using wireless markers enables high quality images to be acquired even during bulk motions. Wireless markers are small, avoid RF safety risks from electrical cables, are not hampered by mechanical connections to the scanner, and require minimal setup times. These advantages may help to facilitate adoption in the clinic. PMID:23813444
Non-rigid alignment in electron tomography in materials science.
Printemps, Tony; Bernier, Nicolas; Bleuet, Pierre; Mula, Guido; Hervé, Lionel
2016-09-01
Electron tomography is a key technique that enables the visualization of an object in three dimensions with a resolution of about a nanometre. High-quality 3D reconstruction is possible thanks to the latest compressed sensing algorithms and/or better alignment and preprocessing of the 2D projections. Rigid alignment of 2D projections is routine in electron tomography. However, it cannot correct misalignments induced by (i) deformations of the sample due to radiation damage or (ii) drifting of the sample during the acquisition of an image in scanning transmission electron microscope mode. In both cases, those misalignments can give rise to artefacts in the reconstruction. We propose a simple-to-implement non-rigid alignment technique to correct those artefacts. This technique is particularly suited for needle-shaped samples in materials science. It is initiated by a rigid alignment of the projections and it is then followed by several rigid alignments of different parts of the projections. Piecewise linear deformations are applied to each projection to force them to simultaneously satisfy the rigid alignments of the different parts. The efficiency of this technique is demonstrated on three samples, an intermetallic sample with deformation misalignments due to a high electron dose typical to spectroscopic electron tomography, a porous silicon sample with an extremely thin end particularly sensitive to electron beam and another porous silicon sample that was drifting during image acquisitions. PMID:27018779
Ultrasound 2D Strain Estimator Based on Image Registration for Ultrasound Elastography
Yang, Xiaofeng; Torres, Mylin; Kirkpatrick, Stephanie; Curran, Walter J.; Liu, Tian
2015-01-01
In this paper, we present a new approach to calculate 2D strain through the registration of the pre- and post-compression (deformation) B-mode image sequences based on an intensity-based non-rigid registration algorithm (INRA). Compared with the most commonly used cross-correlation (CC) method, our approach is not constrained to any particular set of directions, and can overcome displacement estimation errors introduced by incoherent motion and variations in the signal under high compression. This INRA method was tested using phantom and in vivo data. The robustness of our approach was demonstrated in the axial direction as well as the lateral direction where the standard CC method frequently fails. In addition, our approach copes well under large compression (over 6%). In the phantom study, we computed the strain image under various compressions and calculated the signal-to-noise (SNR) and contrast-to-noise (CNS) ratios. The SNR and CNS values of the INRA method were much higher than those calculated from the CC-based method. Furthermore, the clinical feasibility of our approach was demonstrated with the in vivo data from patients with arm lymphedema. PMID:25914492
Ultrasound 2D strain estimator based on image registration for ultrasound elastography
NASA Astrophysics Data System (ADS)
Yang, Xiaofeng; Torres, Mylin; Kirkpatrick, Stephanie; Curran, Walter J.; Liu, Tian
2014-03-01
In this paper, we present a new approach to calculate 2D strain through the registration of the pre- and post-compression (deformation) B-mode image sequences based on an intensity-based non-rigid registration algorithm (INRA). Compared with the most commonly used cross-correlation (CC) method, our approach is not constrained to any particular set of directions, and can overcome displacement estimation errors introduced by incoherent motion and variations in the signal under high compression. This INRA method was tested using phantom and in vivo data. The robustness of our approach was demonstrated in the axial direction as well as the lateral direction where the standard CC method frequently fails. In addition, our approach copes well under large compression (over 6%). In the phantom study, we computed the strain image under various compressions and calculated the signal-to-noise (SNR) and contrast-to-noise (CNS) ratios. The SNR and CNS values of the INRA method were much higher than those calculated from the CC-based method. Furthermore, the clinical feasibility of our approach was demonstrated with the in vivo data from patients with arm lymphedema.
Interactive initialization for 2D/3D intra-operative registration using the Microsoft Kinect
NASA Astrophysics Data System (ADS)
Gong, Ren Hui; Güler, Özgur; Yaniv, Ziv
2013-03-01
All 2D/3D anatomy based rigid registration algorithms are iterative, requiring an initial estimate of the 3D data pose. Current initialization methods have limited applicability in the operating room setting, due to the constraints imposed by this environment or due to insufficient accuracy. In this work we use the Microsoft Kinect device to allow the surgeon to interactively initialize the registration process. A Kinect sensor is used to simulate the mouse-based operations in a conventional manual initialization approach, obviating the need for physical contact with an input device. Different gestures from both arms are detected from the sensor in order to set or switch the required working contexts. 3D hand motion provides the six degree-of-freedom controls for manipulating the pre-operative data in the 3D space. We evaluated our method for both X-ray/CT and X-ray/MR initialization using three publicly available reference data sets. Results show that, with initial target registration errors of 117:7 +/- 28:9 mm a user is able to achieve final errors of 5:9 +/- 2:6 mm within 158 +/- 65 sec using the Kinect-based approach, compared to 4:8+/-2:0 mm and 88+/-60 sec when using the mouse for interaction. Based on these results we conclude that this method is sufficiently accurate for initialization of X-ray/CT and X-ray/MR registration in the OR.
The coupling VIV analysis of SCRs with rigid swing
NASA Astrophysics Data System (ADS)
Liu, Juan; Huang, Weiping
2015-08-01
With the development of deepwater oil and gas exploration, Steel Catenary Risers (SCRs) become preferred risers for resource production, import and export. Vortex induced vibration (VIV) is the key problem encountered in the design of SCRs. In this study, a new model, the rigid swing model, is proposed based on the consideration of large curvature of SCRs. The sag bend of SCRs is assumed as a rigid swing system around the axis from the hanging point to the touch down point (TDP) in the model. The torque, produced by the lift force and the swing vector, provides the driving torque for the swing system, and the weight of SCRs provides the restoring torque. The simulated response of rigid swing is coupled with bending vibration, and then the coupling VIV model of SCRs is studied in consideration of bending vibration and rigid motion. The calculated results indicate that the rigid swing has a magnitude equal to that of bending vibration, and the rigid motion affects the dynamic response of SCRs and can not be neglected in the VIV analysis.
Generalized flexibility-rigidity index
NASA Astrophysics Data System (ADS)
Nguyen, Duc Duy; Xia, Kelin; Wei, Guo-Wei
2016-06-01
Flexibility-rigidity index (FRI) has been developed as a robust, accurate, and efficient method for macromolecular thermal fluctuation analysis and B-factor prediction. The performance of FRI depends on its formulations of rigidity index and flexibility index. In this work, we introduce alternative rigidity and flexibility formulations. The structure of the classic Gaussian surface is utilized to construct a new type of rigidity index, which leads to a new class of rigidity densities with the classic Gaussian surface as a special case. Additionally, we introduce a new type of flexibility index based on the domain indicator property of normalized rigidity density. These generalized FRI (gFRI) methods have been extensively validated by the B-factor predictions of 364 proteins. Significantly outperforming the classic Gaussian network model, gFRI is a new generation of methodologies for accurate, robust, and efficient analysis of protein flexibility and fluctuation. Finally, gFRI based molecular surface generation and flexibility visualization are demonstrated.
Projectile transverse motion and stability in electromagnetic induction launchers
Shokair, I.R.
1993-12-31
The transverse motion of a projectile in an electromagnetic induction launcher is considered. The equations of motion for translation and rotation are derived assuming a rigid projectile and a flyway restoring force per unit length that is proportional to the local displacement. Linearized transverse forces and torques due to energized coils are derived for displaced or tilted armature elements based on a first order perturbation method. The resulting equations of motion for a rigid projectile composed of multiple elements in a multi-coil launcher are analyzed as a coupled oscillator system of equations and a simple linear stability condition is derived. The equations of motion are incorporated into the 2-D Slingshot circuit code and numerical solutions for the transverse motion are obtained. For a launcher with a 10 cm bore radius with a 40 cm long solid armature, we find that stability is achieved with a restoring force (per unit length) constant of k {approx} 1 {times} 10{sup 8} N/m{sup 2}. For k = 1.5 {times} 10{sup 8} N/m{sup 2} and sample coil misalignment modeled as a sine wave of 1 mm amplitude at wavelengths of one or two meters, the projectile displacement grows to a maximum of 4 mm. This growth is due to resonance between the natural frequency of the projectile transverse motion and the coil displacement wavelength. This resonance does not persist because of the changing axial velocity. Random coil displacement is also found to cause roughly the same projectile displacement. For the maximum displacement a rough estimate of the transverse pressure is 50 bars. Results for a wound armature with uniform current density throughout show very similar displacements.
Shirvanyants, David; Alexandrova, Anastassia N.; Dokholyan, Nikolay V.
2011-01-01
Motivation: Identifying the location of binding sites on proteins is of fundamental importance for a wide range of applications, including molecular docking, de novo drug design, structure identification and comparison of functional sites. Here we present Erebus, a web server that searches the entire Protein Data Bank for a given substructure defined by a set of atoms of interest, such as the binding scaffolds for small molecules. The identified substructure contains atoms having the same names, belonging to same amino acids and separated by the same distances (within a given tolerance) as the atoms of the query structure. The accuracy of a match is measured by the root-mean-square deviation or by the normal weight with a given variance. Tests show that our approach can reliably locate rigid binding scaffolds of drugs and metal ions. Availability and Implementation: We provide this service through a web server at http://erebus.dokhlab.org. Contact: dokh@unc.edu PMID:21460026
NASA Astrophysics Data System (ADS)
Latka, Miroslaw; Glaubic-Latka, Marta; Latka, Dariusz; West, Bruce J.
2004-04-01
We study the middle cerebral artery blood flow velocity (MCAfv) in humans using transcranial Doppler ultrasonography (TCD). Scaling properties of time series of the axial flow velocity averaged over a cardiac beat interval may be characterized by two exponents. The short time scaling exponent (STSE) determines the statistical properties of fluctuations of blood flow velocities in short-time intervals while the Hurst exponent describes the long-term fractal properties. In many migraineurs the value of the STSE is significantly reduced and may approach that of the Hurst exponent. This change in dynamical properties reflects the significant loss of short-term adaptability and the overall hyperexcitability of the underlying cerebral blood flow control system. We call this effect fractal rigidity.
Rigid collapsible dish structure
NASA Technical Reports Server (NTRS)
Palmer, William B. (Inventor); Giebler, Martin M. (Inventor)
1982-01-01
A collapsible dish structure composed of a plurality of rows of rigid radial petal assemblies concentric with the axis of the dish. The petal assemblies consist of a center petal and two side petals, the center petal hinged on an axis tangent to a circle concentric with the axis of the dish and the side petals hinged to the center petal at their mating edge. The center petal is foldable inwardly and the side petals rotate about their hinges such that the collapsed dish structure occupies a much smaller volume than the deployed dish. Means of controlling the shape of the dish to compensate for differential expansion of the deployed dish are also provided.
NASA Astrophysics Data System (ADS)
Smania, Daniel
2007-07-01
We describe a new and robust method to prove rigidity results in complex dynamics. The new ingredient is the geometry of the critical puzzle pieces: under control of geometry and ``complex bounds'', two generalized polynomial-like maps which admit a topological conjugacy, quasiconformal outside the filled-in Julia set, are indeed quasiconformally conjugate. The proof uses a new abstract removability-type result for quasiconformal maps, following ideas of Heinonen and Koskela and of Kallunki and Koskela, optimized for applications in complex dynamics. We prove, as the first application of this new method, that, for even criticalities distinct from two, the period two cycle of the Fibonacci renormalization operator is hyperbolic with 1 -dimensional unstable manifold.
The physics of 2D microfluidic droplet ensembles
NASA Astrophysics Data System (ADS)
Beatus, Tsevi; Bar-Ziv, Roy H.; Tlusty, Tsvi
2012-07-01
We review non-equilibrium many-body phenomena in ensembles of 2D microfluidic droplets. The system comprises of continuous two-phase flow with disc-shaped droplets driven in a channel, at low Reynolds number of 10-4-10-3. The basic physics is that of an effective potential flow, governed by the 2D Laplace equation, with multiple, static and dynamic, boundaries of the droplets and the walls. The motion of the droplets induces dipolar flow fields, which mediate 1/r2 hydrodynamic interaction between the droplets. Summation of these long-range 2D forces over droplet ensembles converges, in contrast to the divergence of the hydrodynamic forces in 3D. In analogy to electrostatics, the strong effect of boundaries on the equations of motion is calculated by means of image dipoles. We first consider the dynamics of droplets flowing in a 1D crystal, which exhibits unique phonon-like excitations, and a variety of nonlinear instabilities-all stemming from the hydrodynamic interactions. Narrowing the channel results in hydrodynamic screening of the dipolar interactions, which changes salient features of the phonon spectra. Shifting from a 1D ordered crystal to 2D disordered ensemble, the hydrodynamic interactions induce collective density waves and shocks, which are superposed on single-droplet randomized motion and dynamic clustering. These collective modes originate from density-velocity coupling, whose outcome is a 1D Burgers equation. The rich observational phenomenology and the tractable theory render 2D droplet ensembles a suitable table-top system for studying non-equilibrium many-body physics with long-range interactions.
NASA Astrophysics Data System (ADS)
Searle, G. F. C.
2014-05-01
1. Elementary theory of harmonic motion; 2. Experimental work in harmonic motion; Experiment 1. Determination of g by a simple pendulum; Experiment 2. Harmonic motion of a body suspended by a spring; Experiment 3. Harmonic motion of a rigid body suspended by a torsion wire; Experiment 4. Study of a system with variable moment of inertia; Experiment 5. Dynamical determination of ratio of couple to twist for a torsion wire; Experiment 6. Comparison of the moments of inertia of two bodies; Experiment 7. Experiment with a pair of inertia bars; Experiment 8. Determination of the moment of inertia of a rigid pendulum; Experiment 9. Experiment on a pendulum with variable moment of inertia; Experiment 10. Determination of g by a rigid pendulum; Experiment 11. Pendulum on a yielding support; Experiment 12. Determination of the radius of curvature of a concave mirror by the oscillations of a sphere rolling in it; Experiment 13. Determination of g by the oscillations of a rod rolling on a cylinder; Experiment 14. Study of a vibrating system with two degrees of freedom; Note 1. On the vibration of a body suspended from a light spring; Note 2. Periodic time of a pendulum vibrating through a finite arc; Note 3. Periodic time for finite motion; Note 4. Periodic times of a pendulum with two degrees of freedom.
Direct numerical simulation of rigid bodies in multiphase flow within an Eulerian framework
NASA Astrophysics Data System (ADS)
Rauschenberger, P.; Weigand, B.
2015-06-01
A new method is presented to simulate rigid body motion in the Volume-of-Fluid based multiphase code Free Surface 3D. The specific feature of the new method is that it works within an Eulerian framework without the need for a Lagrangian representation of rigid bodies. Several test cases are shown to prove the validity of the numerical scheme. The technique is able to conserve the shape of arbitrarily shaped rigid bodies and predict terminal velocities of rigid spheres. The instability of a falling ellipsoid is captured. Multiple rigid bodies including collisions may be considered using only one Volume-of-Fluid variable which allows to simulate the drafting, kissing and tumbling phenomena of two rigid spheres. The method can easily be extended to rigid bodies undergoing phase change processes.
Electrically conductive rigid polyurethane foam
Neet, Thomas E.; Spieker, David A.
1985-03-19
A rigid, polyurethane foam comprises about 2-10 weight percent, based on the total foam weight, of a carbon black which is CONDUCTEX CC-40-220 or CONDUCTEX SC, whereby the rigid polyurethane foam is electrically conductive and has essentially the same mechanical properties as the same foam without carbon black added.
Electrically conductive rigid polyurethane foam
Neet, T.E.; Spieker, D.A.
1983-12-08
A rigid, moldable polyurethane foam comprises about 2 to 10 weight percent, based on the total foam weight, of a carbon black which is CONDUCTEX CC-40-220 or CONDUCTEX SC, whereby the rigid polyurethane foam is electrically conductive and has essentially the same mechanical properties as the same foam without carbon black added.
Intraventricular flow alterations due to dyssynchronous wall motion
NASA Astrophysics Data System (ADS)
Pope, Audrey M.; Lai, Hong Kuan; Samaee, Milad; Santhanakrishnan, Arvind
2015-11-01
Roughly 30% of patients with systolic heart failure suffer from left ventricular dyssynchrony (LVD), in which mechanical discoordination of the ventricle walls leads to poor hemodynamics and suboptimal cardiac function. There is currently no clear mechanistic understanding of how abnormalities in septal-lateral (SL) wall motion affects left ventricle (LV) function, which is needed to improve the treatment of LVD using cardiac resynchronization therapy. We use an experimental flow phantom with an LV physical model to study mechanistic effects of SL wall motion delay on LV function. To simulate mechanical LVD, two rigid shafts were coupled to two segments (apical and mid sections) along the septal wall of the LV model. Flow through the LV model was driven using a piston pump, and stepper motors coupled to the above shafts were used to locally perturb the septal wall segments relative to the pump motion. 2D PIV was used to examine the intraventricular flow through the LV physical model. Alterations to SL delay results in a reduction in the kinetic energy (KE) of the flow field compared to synchronous SL motion. The effect of varying SL motion delay from 0% (synchronous) to 100% (out-of-phase) on KE and viscous dissipation will be presented. This research was supported by the Oklahoma Center for Advancement of Science and Technology (HR14-022).
NASA Technical Reports Server (NTRS)
Reddy, T. S. R.; Srivastava, R.
1996-01-01
This guide describes the input data required for using MSAP2D (Multi Stage Aeroelastic analysis Program - Two Dimensional) computer code. MSAP2D can be used for steady, unsteady aerodynamic, and aeroelastic (flutter and forced response) analysis of bladed disks arranged in multiple blade rows such as those found in compressors, turbines, counter rotating propellers or propfans. The code can also be run for single blade row. MSAP2D code is an extension of the original NPHASE code for multiblade row aerodynamic and aeroelastic analysis. Euler equations are used to obtain aerodynamic forces. The structural dynamic equations are written for a rigid typical section undergoing pitching (torsion) and plunging (bending) motion. The aeroelastic equations are solved in time domain. For single blade row analysis, frequency domain analysis is also provided to obtain unsteady aerodynamic coefficients required in an eigen analysis for flutter. In this manual, sample input and output are provided for a single blade row example, two blade row example with equal and unequal number of blades in the blade rows.
Brownian dynamics of confined rigid bodies.
Delong, Steven; Balboa Usabiaga, Florencio; Donev, Aleksandar
2015-10-14
We introduce numerical methods for simulating the diffusive motion of rigid bodies of arbitrary shape immersed in a viscous fluid. We parameterize the orientation of the bodies using normalized quaternions, which are numerically robust, space efficient, and easy to accumulate. We construct a system of overdamped Langevin equations in the quaternion representation that accounts for hydrodynamic effects, preserves the unit-norm constraint on the quaternion, and is time reversible with respect to the Gibbs-Boltzmann distribution at equilibrium. We introduce two schemes for temporal integration of the overdamped Langevin equations of motion, one based on the Fixman midpoint method and the other based on a random finite difference approach, both of which ensure that the correct stochastic drift term is captured in a computationally efficient way. We study several examples of rigid colloidal particles diffusing near a no-slip boundary and demonstrate the importance of the choice of tracking point on the measured translational mean square displacement (MSD). We examine the average short-time as well as the long-time quasi-two-dimensional diffusion coefficient of a rigid particle sedimented near a bottom wall due to gravity. For several particle shapes, we find a choice of tracking point that makes the MSD essentially linear with time, allowing us to estimate the long-time diffusion coefficient efficiently using a Monte Carlo method. However, in general, such a special choice of tracking point does not exist, and numerical techniques for simulating long trajectories, such as the ones we introduce here, are necessary to study diffusion on long time scales. PMID:26472363
Brownian dynamics of confined rigid bodies
NASA Astrophysics Data System (ADS)
Delong, Steven; Balboa Usabiaga, Florencio; Donev, Aleksandar
2015-10-01
We introduce numerical methods for simulating the diffusive motion of rigid bodies of arbitrary shape immersed in a viscous fluid. We parameterize the orientation of the bodies using normalized quaternions, which are numerically robust, space efficient, and easy to accumulate. We construct a system of overdamped Langevin equations in the quaternion representation that accounts for hydrodynamic effects, preserves the unit-norm constraint on the quaternion, and is time reversible with respect to the Gibbs-Boltzmann distribution at equilibrium. We introduce two schemes for temporal integration of the overdamped Langevin equations of motion, one based on the Fixman midpoint method and the other based on a random finite difference approach, both of which ensure that the correct stochastic drift term is captured in a computationally efficient way. We study several examples of rigid colloidal particles diffusing near a no-slip boundary and demonstrate the importance of the choice of tracking point on the measured translational mean square displacement (MSD). We examine the average short-time as well as the long-time quasi-two-dimensional diffusion coefficient of a rigid particle sedimented near a bottom wall due to gravity. For several particle shapes, we find a choice of tracking point that makes the MSD essentially linear with time, allowing us to estimate the long-time diffusion coefficient efficiently using a Monte Carlo method. However, in general, such a special choice of tracking point does not exist, and numerical techniques for simulating long trajectories, such as the ones we introduce here, are necessary to study diffusion on long time scales.
Brownian dynamics of confined rigid bodies
Delong, Steven; Balboa Usabiaga, Florencio; Donev, Aleksandar
2015-10-14
We introduce numerical methods for simulating the diffusive motion of rigid bodies of arbitrary shape immersed in a viscous fluid. We parameterize the orientation of the bodies using normalized quaternions, which are numerically robust, space efficient, and easy to accumulate. We construct a system of overdamped Langevin equations in the quaternion representation that accounts for hydrodynamic effects, preserves the unit-norm constraint on the quaternion, and is time reversible with respect to the Gibbs-Boltzmann distribution at equilibrium. We introduce two schemes for temporal integration of the overdamped Langevin equations of motion, one based on the Fixman midpoint method and the other based on a random finite difference approach, both of which ensure that the correct stochastic drift term is captured in a computationally efficient way. We study several examples of rigid colloidal particles diffusing near a no-slip boundary and demonstrate the importance of the choice of tracking point on the measured translational mean square displacement (MSD). We examine the average short-time as well as the long-time quasi-two-dimensional diffusion coefficient of a rigid particle sedimented near a bottom wall due to gravity. For several particle shapes, we find a choice of tracking point that makes the MSD essentially linear with time, allowing us to estimate the long-time diffusion coefficient efficiently using a Monte Carlo method. However, in general, such a special choice of tracking point does not exist, and numerical techniques for simulating long trajectories, such as the ones we introduce here, are necessary to study diffusion on long time scales.
Registration of 2D cardiac images to real-time 3D ultrasound volumes for 3D stress echocardiography
NASA Astrophysics Data System (ADS)
Leung, K. Y. Esther; van Stralen, Marijn; Voormolen, Marco M.; van Burken, Gerard; Nemes, Attila; ten Cate, Folkert J.; Geleijnse, Marcel L.; de Jong, Nico; van der Steen, Antonius F. W.; Reiber, Johan H. C.; Bosch, Johan G.
2006-03-01
Three-dimensional (3D) stress echocardiography is a novel technique for diagnosing cardiac dysfunction, by comparing wall motion of the left ventricle under different stages of stress. For quantitative comparison of this motion, it is essential to register the ultrasound data. We propose an intensity based rigid registration method to retrieve two-dimensional (2D) four-chamber (4C), two-chamber, and short-axis planes from the 3D data set acquired in the stress stage, using manually selected 2D planes in the rest stage as reference. The algorithm uses the Nelder-Mead simplex optimization to find the optimal transformation of one uniform scaling, three rotation, and three translation parameters. We compared registration using the SAD, SSD, and NCC metrics, performed on four resolution levels of a Gaussian pyramid. The registration's effectiveness was assessed by comparing the 3D positions of the registered apex and mitral valve midpoints and 4C direction with the manually selected results. The registration was tested on data from 20 patients. Best results were found using the NCC metric on data downsampled with factor two: mean registration errors were 8.1mm, 5.4mm, and 8.0° in the apex position, mitral valve position, and 4C direction respectively. The errors were close to the interobserver (7.1mm, 3.8mm, 7.4°) and intraobserver variability (5.2mm, 3.3mm, 7.0°), and better than the error before registration (9.4mm, 9.0mm, 9.9°). We demonstrated that the registration algorithm visually and quantitatively improves the alignment of rest and stress data sets, performing similar to manual alignment. This will improve automated analysis in 3D stress echocardiography.
Supramolecular Synthons: Will Giant Rigid Superspheres Do?
2016-01-01
For the first time, the concept of supramolecular synthons was applied to giant rigid superspheres based on pentaphosphaferrocene [CpRFe(η5-P5)] (R = Me, Et) and Cu(I) halides, which reach 2.1–3.0 nm in diameter. Two supramolecular synthons, σ–π and π–π, are discovered based on halogen···CpR and Cp*···Cp* specific interactions, respectively. The geometry of the synthons is reproducible in a series of crystal structures of various supramolecules. The σ–π synthon alone is realized more frequently for Br-containing superspheres. A combination of the σ–π and π–π synthons is more typical for Cl-containing supramolecules. Each supramolecule can bear up to nine synthons to give mostly 2D and 3D architectures. PMID:27081373
NKG2D ligands as therapeutic targets
Spear, Paul; Wu, Ming-Ru; Sentman, Marie-Louise; Sentman, Charles L.
2013-01-01
The Natural Killer Group 2D (NKG2D) receptor plays an important role in protecting the host from infections and cancer. By recognizing ligands induced on infected or tumor cells, NKG2D modulates lymphocyte activation and promotes immunity to eliminate ligand-expressing cells. Because these ligands are not widely expressed on healthy adult tissue, NKG2D ligands may present a useful target for immunotherapeutic approaches in cancer. Novel therapies targeting NKG2D ligands for the treatment of cancer have shown preclinical success and are poised to enter into clinical trials. In this review, the NKG2D receptor and its ligands are discussed in the context of cancer, infection, and autoimmunity. In addition, therapies targeting NKG2D ligands in cancer are also reviewed. PMID:23833565
New generalizations of the integrable problems in rigid body dynamics
NASA Astrophysics Data System (ADS)
Yehia, H. M.
1997-10-01
We consider the general problem of motion of a rigid body about a fixed point under the action of an axisymmetric combination of potential and gyroscopic forces. We introduce six cases of this problem which are completely integrable for arbitrary initial conditions. The new cases generalize by several parameters all, but one, of the known results in the subject of rigid body dynamics. Namely, we generalize all the results due to Euler, Lagrange, Clebsch, Kovalevskaya, Brun and Lyapunov and also their subsequent generalizations by Rubanovsky and the present author.