Azimuthally Anisotropic 3D Velocity Continuation
Burnett, William; Fomel, Sergey
2011-01-01
We extend time-domain velocity continuation to the zero-offset 3D azimuthally anisotropic case. Velocity continuation describes how a seismic image changes given a change in migration velocity. This description turns out to be of a wave propagation process, in which images change along a velocity axis. In the anisotropic case, the velocity model is multiparameter. Therefore, anisotropic image propagation is multidimensional. We use a three-parameter slowness model, which is related to azimuthal variations in velocity, as well as their principal directions. This information is useful for fracture and reservoir characterization from seismic data. We provide synthetic diffraction imaging examples to illustratemore » the concept and potential applications of azimuthal velocity continuation and to analyze the impulse response of the 3D velocity continuation operator.« less
MOSSFRAC: An anisotropic 3D fracture model
Moss, W C; Levatin, J L
2006-08-14
Despite the intense effort for nearly half a century to construct detailed numerical models of plastic flow and plastic damage accumulation, models for describing fracture, an equally important damage mechanism still cannot describe basic fracture phenomena. Typical fracture models set the stress tensor to zero for tensile fracture and set the deviatoric stress tensor to zero for compressive fracture. One consequence is that the simple case of the tensile fracture of a cylinder under combined compressive radial and tensile axial loads is not modeled correctly. The experimental result is a cylinder that can support compressive radial loads, but no axial load, whereas, the typical numerical result is a cylinder with all stresses equal to zero. This incorrect modeling of fracture locally also has a global effect, because material that is fracturing produces stress release waves, which propagate from the fracture and influence the surrounding material. Consequently, it would be useful to have a model that can describe the stress relief and the resulting anisotropy due to fracture. MOSSFRAC is a material model that simulates three-dimensional tensile and shear fracture in initially isotropic elastic-plastic materials, although its framework is also amenable to initially anisotropic materials. It differs from other models by accounting for the effects of cracks on the constitutive response of the material, so that the previously described experiment, as well as complicated fracture scenarios are simulated more accurately. The model is implemented currently in the LLNL hydrocodes DYNA3D, PARADYN, and ALE3D. The purpose of this technical note is to present a complete qualitative description of the model and quantitative descriptions of salient features.
Vertical Flow Lithography for Fabrication of 3D Anisotropic Particles.
Habasaki, Shohei; Lee, Won Chul; Yoshida, Shotaro; Takeuchi, Shoji
2015-12-22
A microfluidics-based method for the 3D fabrication of anisotropic particles is reported. The method uses a vertical microchannel where tunable light patterns solidify photocurable resins for stacking multiple layers of the resins, thus enabling an application of stereolithography concepts to conventional flow lithography. Multilayered, tapered, and angular compartmental microparticles are demonstrated. PMID:26551590
3D, 9-C anisotropic seismic modeling and inversion
NASA Astrophysics Data System (ADS)
Rusmanugroho, Herurisa
The most complete representation of an elastic medium consists of an elastic tensor with 21 independent moduli. All 21 can be estimated from compressional and shear wave polarization and slowness vectors corresponding to wide apertures of polar and azimuth angles. In isotropic media, when seismic source and receiver components have the same orientation (such as XX and YY), the reflection amplitude contours align approximately perpendicular to the particle motions. The mixed components (such as XY and YX) have amplitude patterns that are in symmetrical pairs of either the same, or of opposite, polarity on either side of the diagonal of the 9-C response matrix. In anisotropic media, amplitude variations with azimuth show the same basic patterns and symmetries as for isotropic, but with a superimposed tendency for alignment parallel to the strike of the vertical cracks. Solutions for elastic tensor elements from synthetic slowness and polarization data calculated directly from the Christoffel equation are more sensitive to the polar angle aperture than to the azimuth aperture. Nine-component synthetic elastic vertical seismic profile data for a model with triclinic symmetry calculated by finite-differencing allows estimation of the elastic 21 tensor elements in the vicinity of a three-component borehole receiver. Wide polar angle and azimuth apertures are needed for accurately estimating the elastic tensor elements. The tensor elements become less independent as the data apertures decrease. Results obtained by extracting slowness and polarization data from the corresponding synthetic seismograms show similar results. The inversion algorithm has produced good results from field vertical seismic profile data set from the Weyburn Field in Southern Saskatchewan in Canada. Synthetic nine-component seismograms calculated from the extracted tensor are able to explain most of the significant features in the field data. The inverted stiffness elastic tensor shows orthorhombic
FPGA-based real-time anisotropic diffusion filtering of 3D ultrasound images
NASA Astrophysics Data System (ADS)
Castro-Pareja, Carlos R.; Dandekar, Omkar S.; Shekhar, Raj
2005-02-01
Three-dimensional ultrasonic imaging, especially the emerging real-time version of it, is particularly valuable in medical applications such as echocardiography, obstetrics and surgical navigation. A known problem with ultrasound images is their high level of speckle noise. Anisotropic diffusion filtering has been shown to be effective in enhancing the visual quality of 3D ultrasound images and as preprocessing prior to advanced image processing. However, due to its arithmetic complexity and the sheer size of 3D ultrasound images, it is not possible to perform online, real-time anisotropic diffusion filtering using standard software implementations. We present an FPGA-based architecture that allows performing anisotropic diffusion filtering of 3D images at acquisition rates, thus enabling the use of this filtering technique in real-time applications, such as visualization, registration and volume rendering.
Lagrange and average interpolation over 3D anisotropic elements
NASA Astrophysics Data System (ADS)
Acosta, Gabriel
2001-10-01
An average interpolation is introduced for 3-rectangles and tetrahedra, and optimal order error estimates in the H1 norm are proved. The constant in the estimate depends "weakly" (improving the results given in Durán (Math. Comp. 68 (1999) 187-199) on the uniformity of the mesh in each direction. For tetrahedra, the constant also depends on the maximum angle of the element. On the other hand, merging several known results (Acosta and Durán, SIAM J. Numer. Anal. 37 (1999) 18-36; Durán, Math. Comp. 68 (1999) 187-199; Krízek, SIAM J. Numer. Anal. 29 (1992) 513-520; Al Shenk, Math. Comp. 63 (1994) 105-119), we prove optimal order error for the -Lagrange interpolation in W1,p, p>2, with a constant depending on p as well as the maximum angle of the element. Again, under the maximum angle condition, optimal order error estimates are obtained in the H1 norm for higher degree interpolations.
Towards an Anisotropic Whole Mantle 3D Elastic Velocity Model
NASA Astrophysics Data System (ADS)
Panning, M. P.; Romanowicz, B.; Gung, Y.
2001-12-01
Many studies have documented the existence of anisotropy in the earth's upper mantle, concentrated in the top 200 km. This evidence comes from the study of surface waves as well as shear wave splitting. There is also evidence for shear wave splitting in D", at least in well sampled regions. There are some hints of anisotropy at the base of the transition zone. Tomographic models of the upper mantle have been developed with simplifying assumptions about the nature of the anisotropy, in order to minimize the number of free parameters in the inversions. Some assume transverse isotropy (e.g Ekström and Dziewonski, 1997), others include additional degrees of freedom with some realistic constraints on mineralogy (e.g. Montagner and Tanimoto, 1991). Our goal is to investigate anisotropy in the whole mantle, using the framework of waveform inversion, and the nonlinear asymptotic mode coupling theory (NACT), previously developed and applied to the construction of whole-mantle SH velocity models (Li and Romanowicz, 1996; Mégnin and Romanowicz, 2000). For this we require a 3 component dataset, and we have extended our automatic transverse (T) component wavepicking procedures to the vertical (Z) and longitudinal (L) component - a non-trivial task given the large number of phases present in the coupled P-SV system. A useful initial assumption, for which the theory has been readily adapted, is that of transverse isotropy. As a first step towards this, we have been investigating inversions using T component and Z,L component data separately. In particular, this allows us to explore the sampling that can be achieved with Z,L component data alone in the deepest part of the mantle. Indeed, D" is in general much better sampled in SH than in SV, owing to the availability of SHdiff at large distances, while SVdiff decays more rapidly due to mantle-core coupling. We present the results of our resolution experiments and discuss the differences between the 3D SV model obtained in well
Qian Chen
2008-08-18
The generation, motion, and interaction of dislocations play key roles during the plastic deformation process of crystalline solids. 3D Dislocation Dynamics has been employed as a mesoscale simulation algorithm to investigate the collective and cooperative behavior of dislocations. Most current research on 3D Dislocation Dynamics is based on the solutions available in the framework of classical isotropic elasticity. However, due to some degree of elastic anisotropy in almost all crystalline solids, it is very necessary to extend 3D Dislocation Dynamics into anisotropic elasticity. In this study, first, the details of efficient and accurate incorporation of the fully anisotropic elasticity into 3D discrete Dislocation Dynamics by numerically evaluating the derivatives of Green's functions are described. Then the intrinsic properties of perfect dislocations, including their stability, their core properties and disassociation characteristics, in newly discovered rare earth-based intermetallics and in conventional intermetallics are investigated, within the framework of fully anisotropic elasticity supplemented with the atomistic information obtained from the ab initio calculations. Moreover, the evolution and interaction of dislocations in these intermetallics as well as the role of solute segregation are presented by utilizing fully anisotropic 3D dislocation dynamics. The results from this work clearly indicate the role and the importance of elastic anisotropy on the evolution of dislocation microstructures, the overall ductility and the hardening behavior in these systems.
3D time-domain airborne EM modeling for an arbitrarily anisotropic earth
NASA Astrophysics Data System (ADS)
Yin, Changchun; Qi, Yanfu; Liu, Yunhe
2016-08-01
Time-domain airborne EM data is currently interpreted based on an isotropic model. Sometimes, it can be problematic when working in the region with distinct dipping stratifications. In this paper, we simulate the 3D time-domain airborne EM responses over an arbitrarily anisotropic earth with topography by edge-based finite-element method. Tetrahedral meshes are used to describe the abnormal bodies with complicated shapes. We further adopt the Backward Euler scheme to discretize the time-domain diffusion equation for electric field, obtaining an unconditionally stable linear equations system. We verify the accuracy of our 3D algorithm by comparing with 1D solutions for an anisotropic half-space. Then, we switch attentions to effects of anisotropic media on the strengths and the diffusion patterns of time-domain airborne EM responses. For numerical experiments, we adopt three typical anisotropic models: 1) an anisotropic anomalous body embedded in an isotropic half-space; 2) an isotropic anomalous body embedded in an anisotropic half-space; 3) an anisotropic half-space with topography. The modeling results show that the electric anisotropy of the subsurface media has big effects on both the strengths and the distribution patterns of time-domain airborne EM responses; this effect needs to be taken into account when interpreting ATEM data in areas with distinct anisotropy.
Anisotropic thermal expansion of a 3D metal–organic framework with hydrophilic and hydrophobic pores
Kondo, Atsushi Maeda, Kazuyuki
2015-01-15
A 3D flexible metal–organic framework (MOF) with 1D hydrophilic and hydrophobic pores shows anisotropic thermal expansion with relatively large thermal expansion coefficient (α{sub a}=−21×10{sup −6} K{sup −1} and α{sub c}=79×10{sup −6} K{sup −1}) between 133 K and 383 K. Temperature change gives deformation of both pores, which expand in diameter and elongate in length on cooling and vice versa. The thermally induced structural change should be derived from a unique framework topology like “lattice fence”. Silica accommodation changes not only the nature of the MOF but also thermal responsiveness of the MOF. Since the hydrophobic pores in the material are selectively blocked by the silica, the MOF with the silica is considered as a hydrophilic microporous material. Furthermore, inclusion of silica resulted in a drastic pore contraction in diameter and anisotropically changed the thermal responsiveness of the MOF. - Graphical abstract: A 3D metal–organic framework with hydrophilic and hydrophobic pores shows anisotropic thermal expansion behavior. The influence of silica filler in the hydrophobic pore was investigated. - Highlights: • Thermally induced structural change of a 3D MOF with a lattice fence topology was investigated. • The structural change was analyzed by synchrotron X-ray diffraction patterns. • Temperature change induces anisotropic thermal expansion/contraction of the MOF. • Silica inclusion anisotropically changes the thermal responsiveness of the MOF.
Howison, Mark
2010-05-06
We compare the performance of hand-tuned CUDA implementations of bilateral and anisotropic diffusion filters for denoising 3D MRI datasets. Our tests sweep comparable parameters for the two filters and measure total runtime, memory bandwidth, computational throughput, and mean squared errors relative to a noiseless reference dataset.
3D digital breast tomosynthesis image reconstruction using anisotropic total variation minimization.
Seyyedi, Saeed; Yildirim, Isa
2014-01-01
This paper presents a compressed sensing based reconstruction method for 3D digital breast tomosynthesis (DBT) imaging. Algebraic reconstruction technique (ART) has been in use in DBT imaging by minimizing the isotropic total variation (TV) of the reconstructed image. The resolution in DBT differs in sagittal and axial directions which should be encountered during the TV minimization. In this study we develop a 3D anisotropic TV (ATV) minimization by considering the different resolutions in different directions. A customized 3D Shepp-logan phantom was generated to mimic a real DBT image by considering the overlapping tissue and directional resolution issues. Results of the ART, ART+3D TV and ART+3D ATV are compared using structural similarity (SSIM) diagram. PMID:25571377
Stress distribution in a premolar 3D model with anisotropic and isotropic enamel.
Munari, Laís S; Cornacchia, Tulimar P M; Moreira, Allyson N; Gonçalves, Jason B; De Las Casas, Estevam B; Magalhães, Cláudia S
2015-08-01
The aim of this study was to compare the areas of stress concentration in a three-dimensional (3D) premolar tooth model with anisotropic or isotropic enamel using the finite element method. A computed tomography was imported to an image processing program to create the tooth model which was exported to a 3D modeling program. The mechanical properties and loading conditions were prescribed in Abaqus. In order to evaluate stresses, axial and oblique loads were applied simulating realistic conditions. Compression stress was observed on the side of load application, and tensile stress was observed on the opposite side. Tensile stress was concentrated mainly in the cervical region and in the alveolar insertion bone. Although stress concentration analyses of the isotropic 3D models produced similar stress distribution results when compared to the anisotropic models, tensile stress values shown by anisotropic models were smaller than the isotropic models. Oblique loads resulted in higher values of tensile stresses, which concentrate mainly in the cervical area of the tooth and in the alveolar bone insertion. Anisotropic properties must be utilized in enamel stress evaluation in non-carious cervical lesions. PMID:25850984
Optimal magnetic susceptibility matching in 3D.
Jia, Feng; Kumar, Rajesh; Korvink, Jan G
2013-04-01
When an object is inserted into the strong homogeneous magnetic field of a magnetic resonance magnet, its intrinsic relative susceptibility can cause unwanted local magnetic field inhomogeneities in the space surrounding the object. As is known, this effect can be partially countered by selectively adding material layers with opposing sign in susceptibility to the part. The determination of an optimal magnetic susceptibility distribution is an inverse problem, in which the susceptibility-induced inhomogeneity of the magnetic field inside a region of interest is reduced by redistributing the placement of materials in the design domain. This article proposes an efficient numerical topology optimization method for obtaining an optimal magnetic susceptibility distribution, in particular, for which the induced spatial magnetic field inhomogeneity is minimized. Using a material density function as a design variable, the value of the magnetic field inside a computational domain is determined using a finite element method. The first-order sensitivity of the objective function is calculated using an adjoint equation method. Numerical examples on a variety of design domain geometries illustrate the effectiveness of the optimization method. The method is of specific interest for the design of interventional magnetic resonance devices. It is a particularly useful method if passive shimming of magnetic resonance equipment is aimed for. PMID:22576319
Time domain topology optimization of 3D nanophotonic devices
NASA Astrophysics Data System (ADS)
Elesin, Y.; Lazarov, B. S.; Jensen, J. S.; Sigmund, O.
2014-02-01
We present an efficient parallel topology optimization framework for design of large scale 3D nanophotonic devices. The code shows excellent scalability and is demonstrated for optimization of broadband frequency splitter, waveguide intersection, photonic crystal-based waveguide and nanowire-based waveguide. The obtained results are compared to simplified 2D studies and we demonstrate that 3D topology optimization may lead to significant performance improvements.
NASA Astrophysics Data System (ADS)
Muench, Thomas; Koch, Manfred; Schlittenhardt, Jörg
2010-05-01
There is now ample evidence from both refraction seismic studies, done already a quarter century ago and from more recent local earthquake traveltime analysis of some of the authors above that large sections of the upper mantle underneath Europe and Germany, in particular, are anisotropic. Employing a modified version of the method of simultaneous inversion for structure and hypocenters (SSH) of the first author, including a priori known upper mantle anisotropy, the investigations of Song et al. [2001] and Song et al. [2004] by means of a 1D time-term analysis and a full 2D Pn anisotropic inversion of regional traveltime data are extended here to a full 3D SSH-inversion underneath Germany. Regional traveltimes from local events occurring between 1975 - 2003 are used which, after application of several selection criteria, results in ~1300 events with a total of ~30000 P- and S-phases for the SSH inversion. Because many of the recorded events appear to suffer from relatively poor hypocentral depth locations a full SSH analysis becomes an intricate undertaking. To alleviate the problem the SSH procedure is carried out in several incremental steps of increasing complexity. First of all improved vertically inhomogeneous velocity (1D) models are derived assuming an isotropic as well as an anisotropic upper mantle. In addition of a slightly better model fit for the anisotropic than for the isotropic model, the latter gives also a somewhat lower Pn-velocity of ~7.90 km/s, compared with ~8.0 km/s for the former. This indicates that inclusion of upper mantle anisotropy into the model is required to obtain physically reasonable Pn-velocities. The results for the P-velocity in the lower crustal layer of the model are less clear, as there appears to be some trade-off in the velocity of that layer and that of the upper mantle. During the course of this part of the study the 3D models have been increasingly refined, starting with a lateral discretization into 15 x 15 blocs
Optimization Techniques for 3D Graphics Deployment on Mobile Devices
NASA Astrophysics Data System (ADS)
Koskela, Timo; Vatjus-Anttila, Jarkko
2015-03-01
3D Internet technologies are becoming essential enablers in many application areas including games, education, collaboration, navigation and social networking. The use of 3D Internet applications with mobile devices provides location-independent access and richer use context, but also performance issues. Therefore, one of the important challenges facing 3D Internet applications is the deployment of 3D graphics on mobile devices. In this article, we present an extensive survey on optimization techniques for 3D graphics deployment on mobile devices and qualitatively analyze the applicability of each technique from the standpoints of visual quality, performance and energy consumption. The analysis focuses on optimization techniques related to data-driven 3D graphics deployment, because it supports off-line use, multi-user interaction, user-created 3D graphics and creation of arbitrary 3D graphics. The outcome of the analysis facilitates the development and deployment of 3D Internet applications on mobile devices and provides guidelines for future research.
3-D Isotropic and Anisotropic S-velocity Structure in the North American Upper Mantle
NASA Astrophysics Data System (ADS)
Yuan, H.; Marone, F.; Romanowicz, B.; Abt, D.; Fischer, K.
2008-12-01
The tectonic diversity of the North American continent has led to a number of geological, tectonic and geodynamical models, many of which can be better tested with high resolution 3-d tomographic models of the isotropic and anisotropic mantle structure of the continent. In the framework of non-linear asymptotic coupling theory (NACT), we recently developed tools to invert long period seismic waveforms combined with SKS splitting data, for both isotropic and radial and azimuthal anisotropic S-wave velocity structure in the upper mantle at the continental scale (Marone et al., 2007; Marone and Romanowicz, 2007). Striking differences in both isotropic and anisotropic velocity structure were observed: beneath the high velocity stable cratonic region a distinct two-layer anisotropic domain is present, with the bottom layer fast axis direction aligned with the absolute plate motion, and a shallower lithospheric layer with north pointing fast axis most likely showing records of past tectonic history; under the active western US the direction of tomographically inferred anisotropy is stable with depth and compatible with the absolute plate motion direction. Here we present an updated model which includes nearly five more years of data, including data from newly operative USArray stations, and a somewhat more extended frequency band. Our new model confirms our previous results, and reveals greater yet complex details of the anisotropic velocity structure beneath the western U.S.. We also show initial results of incorporating constraints on the depth to the lithosphere-asthenosphere boundary (LAB) using teleseismic receiver functions. We discuss the different anisotropic domains resolved both laterally and in depth, in the context of tectonic history of the north American continent.
Parallelization of Program to Optimize Simulated Trajectories (POST3D)
NASA Technical Reports Server (NTRS)
Hammond, Dana P.; Korte, John J. (Technical Monitor)
2001-01-01
This paper describes the parallelization of the Program to Optimize Simulated Trajectories (POST3D). POST3D uses a gradient-based optimization algorithm that reaches an optimum design point by moving from one design point to the next. The gradient calculations required to complete the optimization process, dominate the computational time and have been parallelized using a Single Program Multiple Data (SPMD) on a distributed memory NUMA (non-uniform memory access) architecture. The Origin2000 was used for the tests presented.
Graph-regularized 3D shape reconstruction from highly anisotropic and noisy images
Heinrich, Stephanie; Drewe, Philipp; Lou, Xinghua; Umrania, Shefali; Rätsch, Gunnar
2014-01-01
Analysis of microscopy images can provide insight into many biological processes. One particularly challenging problem is cellular nuclear segmentation in highly anisotropic and noisy 3D image data. Manually localizing and segmenting each and every cellular nucleus is very time-consuming, which remains a bottleneck in large-scale biological experiments. In this work, we present a tool for automated segmentation of cellular nuclei from 3D fluorescent microscopic data. Our tool is based on state-of-the-art image processing and machine learning techniques and provides a user-friendly graphical user interface. We show that our tool is as accurate as manual annotation and greatly reduces the time for the registration. PMID:25866587
Synthetic Seismograms for Realistic 3D Earth Model with Anisotropic Inner Core
NASA Astrophysics Data System (ADS)
Tsuboi, S.; Tono, Y.
2006-12-01
We have demonstrated that we can calculate global theoretical seismograms for realistic 3D Earth models based upon the combination of a precise numerical technique (the spectral-element method) and a sufficiently fast supercomputer (the Earth Simulator) [Tsuboi et al, 2003]. Here we have calculated synthetic seismograms by using model S20RTS of the mantle (Ritsema et al., 1999), model CRUST2.0 of the crust (Basin et al., 2000), topography and bathymetry model ETOPO5, and anisotropic inner core model (Ishii 2002). The calculations are performed on 4056 processors, which require 507 out of 640 nodes of the Earth Simulator. These synthetics are computed by using SPECFEM3D(Komatitsch and Tromp, 2002) and are accurate up to 3.5 seconds. We have calculated these synthetics with aisotropic inner core model for several earthquakes and compared with the synthetics which are calculated for isotropic inner core model. Preliminary comparison shows that the travel time differences between anisotropic inner core model and isotropic core model for PKPab phases are at most a few seconds. There seems to be no significant differences in waveforms of PKP phases. These differences in travel times may help us to improve inner core fine structure by comparing these synthetics with observation.
Measurements of the solid-body rotation of anisotropic particles in 3D turbulence
NASA Astrophysics Data System (ADS)
Marcus, Guy G.; Parsa, Shima; Kramel, Stefan; Ni, Rui; Voth, Greg A.
2014-10-01
We introduce a new method to measure Lagrangian vorticity and the rotational dynamics of anisotropic particles in a turbulent fluid flow. We use 3D printing technology to fabricate crosses (two perpendicular rods) and jacks (three mutually perpendicular rods). Time-resolved measurements of their orientation and solid-body rotation rate are obtained from four video images of their motion in a turbulent flow between oscillating grids with {{R}λ } = 91. The advected particles have a largest dimension of 6 times the Kolmogorov length, making them a good approximation to anisotropic tracer particles. Crosses rotate like disks and jacks rotate like spheres, so these measurements, combined with previous measurements of tracer rods, allow experimental study of axisymmetric ellipsoids across the full range of aspect ratios. The measured mean square tumbling rate, < {{\\dot{p}}i}{{\\dot{p}}i}> , confirms previous direct numerical simulations that indicate that disks tumble much more rapidly than rods. Measurements of the alignment of a unit vector defining the orientation of crosses with the direction of their solid-body rotation rate vector provide the first direct observation of the alignment of anisotropic particles by the velocity gradients in a turbulent flow.
An analysis of 3D anisotropic-viscoelastic forward modeling and dissipation
NASA Astrophysics Data System (ADS)
Yang, Chunying; Li, Xiang-Yang; Wang, Yun
2015-12-01
The anisotropic-viscoelastic wave equation is a generalized expression of a transversely isotropic-viscoelastic medium. The viscoelastic horizontal-tranverse isotropic (HTI) medium is a special case. Considering the horizontal symmetry axis of viscoelastic HTI media, we develop elasticities and wave equations for preserving its transverse isotropy. To model anisotropic-viscoelastic wave propagation, we apply the fourth-order Runge-Kutta and Fourier pseudospectral method to discretize the wave equation. The convolutional perfectly matched layer (CPML) absorbed boundary condition is applied to the 3D modeling algorithm and the result shows that it absorbs reflected energy efficiently. We present three models to investigate anisotropic-viscoelastic waves. Two half-space models demonstrate the azimuthal attenuation, which appears both on the PP and PSv wave. The quality factor is greater along the fracture direction for both the PP and PSv waves. Analysis shows that frequency-dependent amplitude attenuation behaves differently along the fracture azimuth. This demonstates that we can use this property for reservoir detection.
Intelligent speckle reducing anisotropic diffusion algorithm for automated 3-D ultrasound images.
Wu, Jun; Wang, Yuanyuan; Yu, Jinhua; Shi, Xinling; Zhang, Junhua; Chen, Yue; Pang, Yun
2015-02-01
A novel 3-D filtering method is presented for speckle reduction and detail preservation in automated 3-D ultrasound images. First, texture features of an image are analyzed by using the improved quadtree (QT) decomposition. Then, the optimal homogeneous and the obvious heterogeneous regions are selected from QT decomposition results. Finally, diffusion parameters and diffusion process are automatically decided based on the properties of these two selected regions. The computing time needed for 2-D speckle reduction is very short. However, the computing time required for 3-D speckle reduction is often hundreds of times longer than 2-D speckle reduction. This may limit its potential application in practice. Because this new filter can adaptively adjust the time step of iteration, the computation time is reduced effectively. Both synthetic and real 3-D ultrasound images are used to evaluate the proposed filter. It is shown that this filter is superior to other methods in both practicality and efficiency. PMID:26366596
Simulation of a true-triaxial deformation test on anisotropic gneiss using FLAC3D
NASA Astrophysics Data System (ADS)
Ye, Shenghua; Sehizadeh, Mehdi; Nasseri, Mohammad; Young, Paul
2016-04-01
A series of true-triaxial experiments have been carried out at the University of Toronto's Rock Fracture Dynamics Laboratory. Isotropic pegmatite and gneiss have been used to systematically study the effect of anisotropy on the strength, behaviour and seismic response. Samples were loaded under true-triaxial stress conditions and subjected to complex loading and unloading histories associated with rock deformation around underground openings. The results show expected patterns of weakness from preferentially oriented samples and highlight the importance of unloading history under true-triaxial conditions on the deformation and seismic response of the samples. These tests have been used to validate a synthetic simulation using the Itasca FLAC3D numerical code. The paper describes the FLAC3D simulations of the complex true-triaxial loading and unloading history of the different anisotropic samples. Various parameters were created to describe the physico-mechanical properties of the synthetic rock samples. Foliation planes of preferential orientations with respect to the primary loading direction were added to the synthetic rock samples to reflect the anisotropy of the gneiss. These synthetic rock samples were subjected to the same loading and unloading paths as the real rock samples, and failed in the same mechanism as what was observed from the experiments, and thus it proved the validity of this numerical simulation with FLAC3D.
NASA Astrophysics Data System (ADS)
Chong, J.; Yuan, H.; French, S. W.; Romanowicz, B. A.; Ni, S.
2011-12-01
Southeast Asia as a special region in the world which is seismically active and is surrounded by active tectonic belts, such as the Himalaya collision zone, western Pacific subduction zones and the Tianshan- Baikal tectonic belt. Seismic anisotropic tomography can shade light on the complex crust and upper mantle dynamics of this region, which is the subject of much debate. In this study, we applied full waveform time domain tomography to image 3D isotropic and anisotropic upper mantle shear velocity structure of Southeast Asia. Three component waveforms of teleseismic and far regional events (15 degree ≤ Δ≤ 165 degree) with magnitude ranges from Mw6.0 to Mw7.0 are collected from 91 permanent and 438 temporary broadband seismic stations in SE Asia. Wavepackets of both fundamental and overtone modes, filtered between 60 and 400 sec, are selected automatically according to the similarity between data and synthetic waveforms (Panning & Romanowicz, 2006). Wavepackets corresponding to event-station paths that sample the region considered are weighted according to path redundancy and signal to noise ratio. Higher modes and fundamental mode wavepackets are weighted separately in order to enhance the contribution of higher modes which are more sensitive to deeper structure compared to the fundamental mode. Synthetic waveforms and broadband sensitivity kernels are computed using normal mode asymptotic coupling theory (NACT, Li & Romanowicz, 1995). As a starting model, we consider a global anisotropic upper mantle shear velocity model based on waveform inversion using the Spectral Element Method (Lekic & Romanowicz, 2011), updated for more realistic crustal thickness (French et al., 2011) as our starting model, we correct waveforms for the effects of 3D structure outside of the region, and invert them for perturbations in the 3D structure of the target region only. We start with waveform inversion down to 60sec and after several iterations, we include shorter period
An optimal performance control scheme for a 3D crane
NASA Astrophysics Data System (ADS)
Maghsoudi, Mohammad Javad; Mohamed, Z.; Husain, A. R.; Tokhi, M. O.
2016-01-01
This paper presents an optimal performance control scheme for control of a three dimensional (3D) crane system including a Zero Vibration shaper which considers two control objectives concurrently. The control objectives are fast and accurate positioning of a trolley and minimum sway of a payload. A complete mathematical model of a lab-scaled 3D crane is simulated in Simulink. With a specific cost function the proposed controller is designed to cater both control objectives similar to a skilled operator. Simulation and experimental studies on a 3D crane show that the proposed controller has better performance as compared to a sequentially tuned PID-PID anti swing controller. The controller provides better position response with satisfactory payload sway in both rail and trolley responses. Experiments with different payloads and cable lengths show that the proposed controller is robust to changes in payload with satisfactory responses.
3D affine registration using teaching-learning based optimization
NASA Astrophysics Data System (ADS)
Jani, Ashish; Savsani, Vimal; Pandya, Abhijit
2013-09-01
3D image registration is an emerging research field in the study of computer vision. In this paper, two effective global optimization methods are considered for the 3D registration of point clouds. Experiments were conducted by applying each algorithm and their performance was evaluated with respect to rigidity, similarity and affine transformations. Comparison of algorithms and its effectiveness was tested for the average performance to find the global solution for minimizing the error in the terms of distance between the model cloud and the data cloud. The parameters for the transformation matrix were considered as the design variables. Further comparisons of the considered methods were done for the computational effort, computational time and the convergence of the algorithm. The results reveal that the use of TLBO was outstanding for image processing application involving 3D registration. [Figure not available: see fulltext.
Computation of optimized arrays for 3-D electrical imaging surveys
NASA Astrophysics Data System (ADS)
Loke, M. H.; Wilkinson, P. B.; Uhlemann, S. S.; Chambers, J. E.; Oxby, L. S.
2014-12-01
3-D electrical resistivity surveys and inversion models are required to accurately resolve structures in areas with very complex geology where 2-D models might suffer from artefacts. Many 3-D surveys use a grid where the number of electrodes along one direction (x) is much greater than in the perpendicular direction (y). Frequently, due to limitations in the number of independent electrodes in the multi-electrode system, the surveys use a roll-along system with a small number of parallel survey lines aligned along the x-direction. The `Compare R' array optimization method previously used for 2-D surveys is adapted for such 3-D surveys. Offset versions of the inline arrays used in 2-D surveys are included in the number of possible arrays (the comprehensive data set) to improve the sensitivity to structures in between the lines. The array geometric factor and its relative error are used to filter out potentially unstable arrays in the construction of the comprehensive data set. Comparisons of the conventional (consisting of dipole-dipole and Wenner-Schlumberger arrays) and optimized arrays are made using a synthetic model and experimental measurements in a tank. The tests show that structures located between the lines are better resolved with the optimized arrays. The optimized arrays also have significantly better depth resolution compared to the conventional arrays.
Parallel Optimization of 3D Cardiac Electrophysiological Model Using GPU
Xia, Yong; Wang, Kuanquan; Zhang, Henggui
2015-01-01
Large-scale 3D virtual heart model simulations are highly demanding in computational resources. This imposes a big challenge to the traditional computation resources based on CPU environment, which already cannot meet the requirement of the whole computation demands or are not easily available due to expensive costs. GPU as a parallel computing environment therefore provides an alternative to solve the large-scale computational problems of whole heart modeling. In this study, using a 3D sheep atrial model as a test bed, we developed a GPU-based simulation algorithm to simulate the conduction of electrical excitation waves in the 3D atria. In the GPU algorithm, a multicellular tissue model was split into two components: one is the single cell model (ordinary differential equation) and the other is the diffusion term of the monodomain model (partial differential equation). Such a decoupling enabled realization of the GPU parallel algorithm. Furthermore, several optimization strategies were proposed based on the features of the virtual heart model, which enabled a 200-fold speedup as compared to a CPU implementation. In conclusion, an optimized GPU algorithm has been developed that provides an economic and powerful platform for 3D whole heart simulations. PMID:26581957
Eigenfrequency optimized 3D continua, with possibility for cavities
NASA Astrophysics Data System (ADS)
Pedersen, Pauli; Pedersen, Niels L.
2015-04-01
Eigenfrequency optimization for 3D continua is formulated and exemplified by the geometry and boundary conditions of a thick plate. Numerical finite element models are based on four node tetrahedra and results from subspace iterations give directly the basis for the continuum redesign. The 3D modeling with a large number of elements has the possibility in optimal design to obtain (as found) not only holes but also cavities in the continuum. Sensitivity analysis is presented on the element level with simple physical interpretation of the involved terms. This general result has general value for control of eigenfrequencies. It is found that in the combination of partial differentiation with the chain rule of differentiation, a specific notation is needed and a suggestion is presented. The optimization method is based on a derived optimality criterion, and as such the maximization problem change to a problem of determining a design with uniform values of this criterion. Nonlinear stiffness interpolation may be a physical reality. A two parameter interpolation function is incorporated analytical, also in the sensitivity analysis and the optimality criterion, but without focusing on 1-0 optimal solutions. Two cases of boundary conditions, two cases of total amount of material, and cases of linear and nonlinear stiffness interpolation are studied.
3D Anisotropic structure of the south-central Mongolia from Rayleigh and Love wave tomography
NASA Astrophysics Data System (ADS)
Yu, D.; Wu, Q.; Montagner, J. P.
2014-12-01
A better understanding of the geodynamics of the crust and mantle below Baikal-Mongolia is required to identify the role of mantle processes versus that of far-field tectonic effects from India-Asia collision. Anisotropy tomography can provide new perspective to the continental growth mechanism. In order to study the 3D anisotropic structure of the upper mantle in the south-central Mongolia, we collected the vertical and transverse components of seismograms recorded at 69 broadband seismic stations. We have measured inter-station phase velocities of 7181 Rayleigh waves and 901 Love waves using the frequency-time analysis of wavelet transformation method for the fundamental mode at period range 10~80s. The lateral phase velocity variations are computed by using a regionalization method. These phase velocities have been inverted to obtain the first anisotropic model including Sv velocities, azimuthal and radial anisotropy. The Middle Gobi is associated with low velocity. Based on the distribution of the Cenozoic basalts in the Middle Gobi, it refers that the low velocity anomaly is related to the Cenozoic volcanism. In the northern domain, near to Baikal zone, the azimuthal anisotropy is normal to the Baikal rift and consistent with the fast direction of previous SKS splitting measurements. In the South Gobi, north to Main Mongolian Lineament, the azimuthal anisotropy is NEE-SWW in the crust and NW-SE in the mantle. It indicates that the crust and mantle are decoupled. We propose that the crustal deformation is related to the far-field effects of India-Asia collision and that the mantle flow is correlated with the Baikal rift activity. Further study in process will provide more evidence and insight to better understand the geodynamics in this region.
Optimized 3D Ultrashort Echo Time Pulmonary MRI
Johnson, Kevin M.; Fain, Sean B.; Schiebler, Mark L.; Nagle, Scott
2012-01-01
Purpose To optimize 3D radial ultrashort echo time MRI for high resolution whole-lung imaging. Methods 3D radial ultrashort echo time was implemented on a 3T scanner to investigate the effects of: (1) limited field-of-view excitation, (2) variable density readouts, and (3) radial oversampling. Improvements in noise performance and spatial resolution were assessed through simulation and phantom studies. Their effects on lung and airway visualization in five healthy male human subjects (mean age 32 years) were compared qualitatively through blinded ordinal scoring by two cardiothoracic radiologists using a nonparametric Friedman test (P < 0.05). Relative signal difference between endobronchial air and adjacent lung tissue, normalized to nearby vessel, was used as a surrogate for lung tissue signal. Quantitative measures were compared using the paired Student's t-test (P < 0.05). Finally, clinical feasibility was investigated in a patient with interstitial fibrosis. Results Simulation and phantom studies showed up to 67% improvement in SNR and reduced blurring for short T2* species using all three optimizations. In vivo images showed decreased artifacts and improved lung tissue and airway visualization both qualitatively and quantitatively. Conclusion The use of limited field-of-view excitation, variable readout gradients, and radial oversampling significantly improve the technical quality of 3D radial ultrashort echo time lung images. PMID:23213020
Transient 3D elastodynamic field in an embedded multilayered anisotropic plate.
Mora, Pierric; Ducasse, Eric; Deschamps, Marc
2016-07-01
The aim of this paper is to study the ultrasonic response to a transient source that radiates ultrasonic waves in a 3D embedded multilayered anisotropic and dissipative plate. The source can be inside the plate or outside, in a fluid loading the plate for example. In the context of Non-Destructive Testing applied to composite materials, our goal is to create a robust algorithm to calculate ultrasonic field, irrespective of the source and receiver positions. The principle of the method described in this paper is well-established. This method is based on time analysis using the Laplace transform. In the present work, it has been customized for computing ultrasonic source interactions with multilayered dissipative anisotropic plates. The fields are transformed in the 2D Fourier wave-vector domain for the space variables related to the plate surface, and they are expressed in the partial-wave basis. Surprisingly, this method has been very little used in the ultrasonic community, while it is a useful tool which complements the much used technique based on generalized Lamb wave decomposition. By avoiding mode analysis - which can be problematic in some cases - exact numerical calculations (i.e., approximations by truncating infinite series that may be poorly convergent are not needed) can be made in a relatively short time for immersed plates and viscoelastic layers. Even for 3D cases, numerical costs are relatively low. Special attention is given to separate up- and down-going waves, which is a simple matter when using the Laplace transform. Numerical results show the effectiveness of this method. Three examples are presented here to investigate the quality of the model and the robustness of the algorithm: first, a comparison of experiment and simulation for a monolayer carbon-epoxy plate, where the diffracted field is due to a source located on the first free surface of the sample, for both dissipative and non-dissipative cases; second, the basic configuration of an
PDE constrained optimization of electrical defibrillation in a 3D ventricular slice geometry.
Chamakuri, Nagaiah; Kunisch, Karl; Plank, Gernot
2016-04-01
A computational study of an optimal control approach for cardiac defibrillation in a 3D geometry is presented. The cardiac bioelectric activity at the tissue and bath volumes is modeled by the bidomain model equations. The model includes intramural fiber rotation, axially symmetric around the fiber direction, and anisotropic conductivity coefficients, which are extracted from a histological image. The dynamics of the ionic currents are based on the regularized Mitchell-Schaeffer model. The controls enter in the form of electrodes, which are placed at the boundary of the bath volume with the goal of dampening undesired arrhythmias. The numerical optimization is based on Newton techniques. We demonstrated the parallel architecture environment for the computation of potentials on multidomains and for the higher order optimization techniques. Copyright © 2015 John Wiley & Sons, Ltd. PMID:26249168
Anisotropic intrinsic anomalous Hall effect in ordered 3dPt alloys
NASA Astrophysics Data System (ADS)
Zhang, Hongbin; Blügel, Stefan; Mokrousov, Yuriy
2011-07-01
By performing first-principles calculations, we investigate the intrinsic anomalous Hall conductivity (AHC) and its anisotropy in ordered L10 FePt, CoPt, and NiPt ferromagnets and their intermediate alloys. We demonstrate that the AHC in this family of compounds depends strongly on the direction of the magnetization M in the crystal. We predict that such pronounced orientational dependence in combination with the general decreasing trend of the AHC when going from FePt to NiPt leads to a sign change of the AHC upon rotating the magnetization direction in the crystal of CoPt alloy. We also suggest that, for a range of concentration x in CoxNi1-xPt and FexCo1-xPt alloys, it is possible to achieve a complete quenching of the anomalous Hall current for a certain direction of the magnetization in the crystal. By analyzing the spin-resolved AHC in 3dPt alloys, we endeavor to relate the overall trend of the AHC in these compounds to the changes in their densities of d states around the Fermi energy upon varying the atomic number. Moreover, we show the generality of the phenomenon of anisotropic anomalous Hall effect by demonstrating its occurrence within the three-band tight-binding model.
Improved hybrid optimization algorithm for 3D protein structure prediction.
Zhou, Changjun; Hou, Caixia; Wei, Xiaopeng; Zhang, Qiang
2014-07-01
A new improved hybrid optimization algorithm - PGATS algorithm, which is based on toy off-lattice model, is presented for dealing with three-dimensional protein structure prediction problems. The algorithm combines the particle swarm optimization (PSO), genetic algorithm (GA), and tabu search (TS) algorithms. Otherwise, we also take some different improved strategies. The factor of stochastic disturbance is joined in the particle swarm optimization to improve the search ability; the operations of crossover and mutation that are in the genetic algorithm are changed to a kind of random liner method; at last tabu search algorithm is improved by appending a mutation operator. Through the combination of a variety of strategies and algorithms, the protein structure prediction (PSP) in a 3D off-lattice model is achieved. The PSP problem is an NP-hard problem, but the problem can be attributed to a global optimization problem of multi-extremum and multi-parameters. This is the theoretical principle of the hybrid optimization algorithm that is proposed in this paper. The algorithm combines local search and global search, which overcomes the shortcoming of a single algorithm, giving full play to the advantage of each algorithm. In the current universal standard sequences, Fibonacci sequences and real protein sequences are certified. Experiments show that the proposed new method outperforms single algorithms on the accuracy of calculating the protein sequence energy value, which is proved to be an effective way to predict the structure of proteins. PMID:25069136
3D Plasma Equilibrium and Stability with Hot Particle Anisotropic Pressure
Cooper, W. A.; Graves, J. P.; Hirshman, S. P.; Merkel, P.; Kisslinger, J.; Wobig, H. F. G.; Watanabe, K. Y.; Narushima, Y.
2008-11-01
The anisotropic pressure free-boundary three-dimsnsional (3D) equilibrium code ANI-MEC with nested magnetic flux surfaces has been developed as an extension of the VMEC2000 code. The preconditioning algorithm included is exploited to allow the computation of equilibrium states with radial force balance error improvements exceeding 4 orders of magnitude compared with the non-conditioned results. Large off-axis energetic particle deposition has been applied in a 2-field period quasiaxisymmetric stellarator reactor at <{beta}>{approx_equal}4.5% to test the limitations of the code. The hot particle pressures are roughly uniform around the flux surfaces when p{sub parallel}>p{sub perpendicular}. The fast particle perpendicular pressures localise in the region of deposition for p{sub perpendicular}>p{sub parallel}, while the energetic particle parallel pressures concentrate on the low-field side. Two anisotropic pressure models for global fluid stability implemented in the TERPSICHORE code have been applied to the LHD Heliotron for a sequence of equilibria with fixed <{beta}{sub dia}>{approx_equal}5%(<{beta}{sub th}>{approx_equal}3.5%) varying the fast particle temperature ratio T{sub parallel}/T{sub perpendicular}. Global magnetohydrodynamic modes are quasi-stable according to the model with rigid hot particle layers, while they become stabilised according to the fully interacting energetic particle model with increasing T{sub parallel}/T{sub perpendicular}. As T{sub parallel}/T{sub perpendicular} approaches 3, however, the n = 1 mode family becomes unstable. A transition from a nearly stable quasi-external ballooning-interchange structure to a weakly unstable internal kink mode takes place. The investigation of beam-driven fusion in a Heliotron system is broached. A background plasma with cold ions and warm electrons at <{beta}{sub ith}>{approx_equal}1% is examined with fixed T{sub parallel}/T{sub perpendicular} = 10 in which the hot particle contribution to <{beta
Anisotropic magnification distortion of the 3D galaxy correlation. I. Real space
Hui, Lam; LoVerde, Marilena; Gaztanaga, Enrique
2007-11-15
It has long been known that gravitational lensing, primarily via magnification bias, modifies the observed galaxy (or quasar) clustering. Such discussions have largely focused on the 2D angular correlation function. Here and in paper II [L. Hui, E. Gaztanaga, and M. LoVerde, arXiv:0710.4191] we explore how magnification bias distorts the 3D correlation function and power spectrum, as first considered by Matsubara [Astrophys. J. Lett. 537, L77 (2000).]. The interesting point is that the distortion is anisotropic. Magnification bias in general preferentially enhances the observed correlation in the line-of-sight (LOS) orientation, especially on large scales. For instance, at a LOS separation of {approx}100 Mpc/h, where the intrinsic galaxy-galaxy correlation is rather weak, the observed correlation can be enhanced by lensing by a factor of a few, even at a modest redshift of z{approx}0.35. This effect presents an interesting opportunity as well as a challenge. The opportunity: this lensing anisotropy is distinctive, making it possible to separately measure the galaxy-galaxy, galaxy-magnification, and magnification-magnification correlations, without measuring galaxy shapes. The anisotropy is distinguishable from the well-known distortion due to peculiar motions, as will be discussed in paper II. The challenge: the magnification distortion of the galaxy correlation must be accounted for in interpreting data as precision improves. For instance, the {approx}100 Mpc/h baryon acoustic oscillation scale in the correlation function is shifted by up to {approx}3% in the LOS orientation, and up to {approx}0.6% in the monopole, depending on the galaxy bias, redshift, and number count slope. The corresponding shifts in the inferred Hubble parameter and angular diameter distance, if ignored, could significantly bias measurements of the dark energy equation of state. Lastly, magnification distortion offers a plausible explanation for the well-known excess correlations seen in
NASA Astrophysics Data System (ADS)
Hofmann, I.; Boine-Frankenheim, O.
2014-12-01
The numerical noise inherent to particle-in-cell (PIC) simulation of 3d anisotropic high intensity bunched beams in periodic focusing is compared with the analytical model by Struckmeier [Part. Accel. 45, 229 (1994)]. The latter assumes that entropy growth can be related to Markov type stochastic processes due to temperature anisotropy and the artificial "collisions" caused by using macro-particles and calculating the space charge effect. The PIC simulations are carried out with the tracewin code widely used for high intensity beam simulation. The resulting noise can lead to growth of the six-dimensional rms emittance. The logarithm of the latter is shown to qualify as rms-based entropy. We confirm the dependence of this growth on the bunch temperature anisotropy as predicted by Struckmeier. However, we also find a grid and focusing dependent component of noise not predicted by Struckmeier. Although commonalities exist with well-established models for collision effects in PIC-simulation of extended plasmas, a distinctive feature is the presence of a periodic focusing potential, wherein the beam one-component plasma extends only over relatively few Debye lengths. Our findings are applied in particular to noise in high current linac beam simulation, where they help for optimization of the balance between the number of simulation particles and the grid resolution.
New X-Ray Tomography Method Based on the 3D Radon Transform Compatible with Anisotropic Sources
NASA Astrophysics Data System (ADS)
Vassholz, M.; Koberstein-Schwarz, B.; Ruhlandt, A.; Krenkel, M.; Salditt, T.
2016-02-01
In this work, we propose a novel computed tomography (CT) approach for three-dimensional (3D) object reconstruction, based on a generalized tomographic geometry with two-dimensional angular sampling (two angular degrees of freedom). The reconstruction is based on the 3D radon transform and is compatible with anisotropic beam conditions. This allows isotropic 3D imaging with a source, which can be extended along one direction for increased flux, while high resolution is achieved by a small source size only in the orthogonal direction. This novel scheme for analytical CT is demonstrated by numerical simulations and proof-of-concept experiments. In this way high resolution and coherence along a single direction determines the reconstruction quality of the entire 3D data set, opening up, for example, new opportunities to achieve nanoscale resolution and/or phase contrast with low brilliance sources such as laboratory x-ray or neutron sources.
New X-Ray Tomography Method Based on the 3D Radon Transform Compatible with Anisotropic Sources.
Vassholz, M; Koberstein-Schwarz, B; Ruhlandt, A; Krenkel, M; Salditt, T
2016-02-26
In this work, we propose a novel computed tomography (CT) approach for three-dimensional (3D) object reconstruction, based on a generalized tomographic geometry with two-dimensional angular sampling (two angular degrees of freedom). The reconstruction is based on the 3D Radon transform and is compatible with anisotropic beam conditions. This allows isotropic 3D imaging with a source, which can be extended along one direction for increased flux, while high resolution is achieved by a small source size only in the orthogonal direction. This novel scheme for analytical CT is demonstrated by numerical simulations and proof-of-concept experiments. In this way high resolution and coherence along a single direction determines the reconstruction quality of the entire 3D data set, opening up, for example, new opportunities to achieve nanoscale resolution and/or phase contrast with low brilliance sources such as laboratory x-ray or neutron sources. PMID:26967444
Optimizing Stellarators for Energetic Particle Confinement using BEAMS3D
NASA Astrophysics Data System (ADS)
Bolgert, Peter; Drevlak, Michael; Lazerson, Sam; Gates, David; White, Roscoe
2015-11-01
Energetic particle (EP) loss has been called the ``Achilles heel of stellarators,'' (Helander, Rep. Prog. Phys. 77 087001 (2014)) and there is a great need for magnetic configurations with improved EP confinement. In this study we utilize a newly developed capability of the stellarator optimization code STELLOPT: the ability to optimize EP confinement via an interface with guiding center code BEAMS3D (McMillan et al., Plasma Phys. Control. Fusion 56, 095019 (2014)). Using this new tool, optimizations of the W7-X experiment and ARIES-CS reactor are performed where the EP loss fraction is one of many target functions to be minimized. In W7-X, we simulate the experimental NBI system using realistic beam geometry and beam deposition physics. The goal is to find configurations with improved neutral beam deposition and energetic particle confinement. These calculations are compared to previous studies of W7-X NBI deposition. In ARIES-CS, we launch 3.5 MeV alpha particles from a near-axis flux surface using a uniform grid in toroidal and poloidal angle. As these particles are born from D-T reactions, we consider an isotropic distribution in velocity space. This research is supported by DoE Contract Number DE-AC02-09CH11466.
NASA Astrophysics Data System (ADS)
Greenhalgh, Stewart; Zhou, Bing; Maurer, Hansruedi
2010-05-01
We have developed a modified version of the spectral element method (SEM), called the Gaussian Quadrature Grid (GQG) approach, for frequency domain 3D seismic modelling in arbitrary heterogeneous, anisotropic media. The model may incorporate an arbitrary free-surface topography and irregular subsurface interfaces. Unlike the SEM ,it does not require a powerful mesh generator such as the Delauney Triangular or TetGen. Rather, the GQG approach replaces the element mesh with Gaussian quadrature abscissae to directly sample the physical properties of the model parameters and compute the weighted residual or variational integral. This renders the model discretisation simple and easily matched to the model topography, as well as direct control of the model paramterisation for subsequent inversion. In addition, it offers high accuracy in numerical modelling provided that an appropriate density of the Gaussian quadrature abscissae is employed. The second innovation of the GQG is the incorporation of a new implementation of perfectly matched layers to suppress artificial reflections from the domain margins. We employ PML model parameters (specified complex valued density and elastic moduli) rather than explicitly solving the governing wave equation with a complex co-ordinate system as in conventional approaches. Such an implementation is simple, general, effective and easily extendable to any class of anisotropy and other numerical modelling methods. The accuracy of the GQG approach is controlled by the number of Gaussian quadrature points per minimum wavelength, the so-called sampling density. The optimal sampling density should be the one which enables high definition of geological characteristics and high precision of the variational integral evaluation and spatial differentiation. Our experiments show that satisfactory results can be obtained using sampling densities of 5 points per minimum wavelength. Efficiency of the GQG approach mainly depends on the linear
NASA Astrophysics Data System (ADS)
Okaya, D. A.; Van Avendonk, H. J.
2013-12-01
Recent anisotropy studies at scales ranging from crust to full mantle have recognized the importance of 3D anisotropy geometry and heterogeneity as well as variability in anisotropic symmetry and orientation (tilt) of the Earth. The strong relationship between seismic anisotropy and geodynamic processes highlights the need to construct realistic Earth models that can explain observations of anisotropy in modern seismic data sets. For example, ray paths through a mantle slab window or a mountain belt may show that the crust or mantle exhibits low-order anisotropy due to a history of deformation and the development of tectonic fabrics. Observed traveltimes might not be fit with simple Transverse Isotropy (TI), so realistic calculations require an Earth model that accurately describes the wave speeds of compressional and shear waves. We have developed an anisotropic traveltime solver that allows for full 3D heterogeneity of anisotropy tensors, degrees of symmetry, and arbitrary orientation. This traveltime solver is based on the robust shortest path method (SPM) and a ray-bending algorithm that were previously applied to isotropic media (e.g., Van Avendonk et al., 2001). Instead of using an isotropic description of the seismic wave velocity, we define the full elastic tensor at each location in the model. The directional seismic velocity can subsequently be extracted using solutions of the Christoffel equations. For computational efficiency, we calculate all directional seismic velocities at each model node before the start of ray tracing. As we calculate a new ray segment, this information is quickly retrieved. We use these directional velocity maps to separately describe the propagation of compressional (P) and shear (S) body waves in anisotropic media and to subsequently calculate their traveltimes. Patterns within the velocity maps represent tensor symmetries and tilts, allowing for the construction of discretized large-scale 3D LPO flow fields or fabric
Debbaut, Charlotte; Vierendeels, Jan; Siggers, Jennifer H; Repetto, Rodolfo; Monbaliu, Diethard; Segers, Patrick
2014-01-01
The hepatic blood circulation is complex, particularly at the microcirculatory level. Previously, 2D liver lobule models using porous media and a 3D model using real sinusoidal geometries have been developed. We extended these models to investigate the role of vascular septa (VS) and anisotropic permeability. The lobule was modelled as a hexagonal prism (with or without VS) and the tissue was treated as a porous medium (isotropic or anisotropic permeability). Models were solved using computational fluid dynamics. VS inclusion resulted in more spatially homogeneous perfusion. Anisotropic permeability resulted in a larger axial velocity component than isotropic permeability. A parameter study revealed that results are most sensitive to the lobule size and radial pressure drop. Our model provides insight into hepatic microhaemodynamics, and suggests that inclusion of VS in the model leads to perfusion patterns that are likely to reflect physiological reality. The model has potential for applications to unphysiological and pathological conditions. PMID:23237543
NASA Astrophysics Data System (ADS)
Burkett, Michael W.; Clancy, Sean P.; Maudlin, Paul J.; Holian, Kathleen S.
2002-07-01
Previously developed constitutive models and solution algorithms for anisotropic elastoplastic material strength have been implemented in the three-dimensional Conejo hydrodynamics code. The anisotropic constitutive modeling is posed in an unrotated material frame of reference using the theorem of polar decomposition to obtain rigid body rotation. Continuous quadratic yield functions fitted from polycrystal simulations for a metallic hexagonal-close-packed structure were utilized. Simple rectangular shear problems, R-Value problems, and Taylor cylinder impact data were used to verify and validate the implementation of the anisotropic model. A stretching rod problem (involving large strain and high strain-rate deformation) was selected to investigate the effects of material anisotropy. Conejo simulations of rod topology were compared for two anisotropic cases.
Optimizing prostate needle biopsy through 3D simulation
NASA Astrophysics Data System (ADS)
Zeng, Jianchao; Kaplan, Charles; Xuan, Jian Hua; Sesterhenn, Isabell A.; Lynch, John H.; Freedman, Matthew T.; Mun, Seong K.
1998-06-01
Prostate needle biopsy is used for the detection of prostate cancer. The protocol of needle biopsy that is currently routinely used in the clinical environment is the systematic sextant technique, which defines six symmetric locations on the prostate surface for needle insertion. However, this protocol has been developed based on the long-term observation and experience of urologists. Little quantitative or scientific evidence supports the use of this biopsy technique. In this research, we aim at developing a statistically optimized new prostate needle biopsy protocol to improve the quality of diagnosis of prostate cancer. This new protocol will be developed by using a three-dimensional (3-D) computer- based probability map of prostate cancer. For this purpose, we have developed a computer-based 3-D visualization and simulation system with prostate models constructed from the digitized prostate specimens, in which the process of prostate needle biopsy can be simulated automatically by the computer. In this paper, we first develop an interactive biopsy simulation mode in the system, and evaluate the performance of the automatic biopsy simulation with the sextant biopsy protocol by comparing the results by the urologist using the interactive simulation mode with respect to 53 prostate models. This is required to confirm that the automatic simulation is accurate and reliable enough for the simulation with respect to a large number of prostate models. Then we compare the performance of the existing protocols using the automatic biopsy simulation system with respect to 107 prostate models, which will statistically identify if one protocol is better than another. Since the estimation of tumor volume is extremely important in determining the significance of a tumor and in deciding appropriate treatment methods, we further investigate correlation between the tumor volume and the positive core volume with 89 prostate models. This is done in order to develop a method to
A 3-D elasticity theory based model for acoustic radiation from multilayered anisotropic plates.
Shen, C; Xin, F X; Lu, T J
2014-05-01
A theoretical model built upon three-dimensional elasticity theory is developed to investigate the acoustic radiation from multilayered anisotropic plates subjected to a harmonic point force excitation. Fourier transform technique and stationary phase method are combined to predict the far-field radiated sound pressure of one-side water immersed plate. Compared to equivalent single-layer plate models, the present model based on elasticity theory can differentiate radiated sound pressure between dry-side and wet-side excited cases, as well as discrepancies induced by different layer sequences for multilayered anisotropic plates. These results highlight the superiority of the present theoretical model especially for handling multilayered anisotropic structures. PMID:24815294
NASA Astrophysics Data System (ADS)
Burkett, Michael; Clancy, Sean; Maudlin, Paul; Holian, Kathleen
2001-06-01
: Previously developed constitutive models and solution algorithms for anisotropic elastoplastic material strength has been implemented in the three-dimensional CONEJO hydrodynamics code. CONEJO is an explicit, Eulerian continuum mechanics code that is utilized to predict formation processes associated with material deformation at elevated strain-rates and is a code development project under the Accelerated Strategic Computing Initiative (ASCI) program. Some special features of CONEJO include a high-order advection algorithm, a material interface tracking scheme, and van Leer monotonic advection-limiting. The anisotropic constitutive modeling is posed in an unrotated material frame using the theorem of polar decomposition to describe rigid body rotation. An Euler-Rodrigues description is used to quantify the rigid body rotations. Continuous quadratic yield functions fitted from polycrystal simulations for a metallic hexagonal-close-packed structure were utilized. Associative flow formulations incorporating these yield functions were solved using a geometric normal return method. Simple rectangular shear problems, "R-value" problems, and Taylor cylinder impact test data were utilized to verify and validate the implementation of the anisotropic model. A "stretching rod" problem (involving large strain and strain-rate deformation) was selected to investigate the effects of material anisotropy for this deformation process. The rod necking rate and topology was compared for CONEJO simulations using several isotropic and anisotropic descriptions that utilized the Mechanical Threshold Stress (MTS) model.
NASA Astrophysics Data System (ADS)
Gansen, A.; Hachemi, M. El; Belouettar, S.; Hassan, O.; Morgan, K.
2016-09-01
The standard Yee algorithm is widely used in computational electromagnetics because of its simplicity and divergence free nature. A generalization of the classical Yee scheme to 3D unstructured meshes is adopted, based on the use of a Delaunay primal mesh and its high quality Voronoi dual. This allows the problem of accuracy losses, which are normally associated with the use of the standard Yee scheme and a staircased representation of curved material interfaces, to be circumvented. The 3D dual mesh leapfrog-scheme which is presented has the ability to model both electric and magnetic anisotropic lossy materials. This approach enables the modelling of problems, of current practical interest, involving structured composites and metamaterials.
3D head model classification using optimized EGI
NASA Astrophysics Data System (ADS)
Tong, Xin; Wong, Hau-san; Ma, Bo
2006-02-01
With the general availability of 3D digitizers and scanners, 3D graphical models have been used widely in a variety of applications. This has led to the development of search engines for 3D models. Especially, 3D head model classification and retrieval have received more and more attention in view of their many potential applications in criminal identifications, computer animation, movie industry and medical industry. This paper addresses the 3D head model classification problem using 2D subspace analysis methods such as 2D principal component analysis (2D PCA[3]) and 2D fisher discriminant analysis (2DLDA[5]). It takes advantage of the fact that the histogram is a 2D image, and we can extract the most useful information from these 2D images to get a good result accordingingly. As a result, there are two main advantages: First, we can perform less calculation to obtain the same rate of classification; second, we can reduce the dimensionality more than PCA to obtain a higher efficiency.
Optimizing 3D image quality and performance for stereoscopic gaming
NASA Astrophysics Data System (ADS)
Flack, Julien; Sanderson, Hugh; Pegg, Steven; Kwok, Simon; Paterson, Daniel
2009-02-01
The successful introduction of stereoscopic TV systems, such as Samsung's 3D Ready Plasma, requires high quality 3D content to be commercially available to the consumer. Console and PC games provide the most readily accessible source of high quality 3D content. This paper describes innovative developments in a generic, PC-based game driver architecture that addresses the two key issues affecting 3D gaming: quality and speed. At the heart of the quality issue are the same considerations that studios face producing stereoscopic renders from CG movies: how best to perform the mapping from a geometric CG environment into the stereoscopic display volume. The major difference being that for game drivers this mapping cannot be choreographed by hand but must be automatically calculated in real-time without significant impact on performance. Performance is a critical issue when dealing with gaming. Stereoscopic gaming has traditionally meant rendering the scene twice with the associated performance overhead. An alternative approach is to render the scene from one virtual camera position and use information from the z-buffer to generate a stereo pair using Depth-Image-Based Rendering (DIBR). We analyze this trade-off in more detail and provide some results relating to both 3D image quality and render performance.
The numerical simulation for a 3D two-phase anisotropic medium based on BISQ model
NASA Astrophysics Data System (ADS)
Wang, Zhejiang; He, Qiaodeng; Wang, Deli
2008-03-01
Biot-flow and squirt-flow are the two most important fluid flow mechanisms in porous media containing fluids. Based on the BISQ (Biot-Squirt) model where the two mechanisms are treated simultaneously, the elastic wave-field simulation in the porous medium is limited to two-dimensions and two-components (2D2C) or two-dimensions and three-components (2D3C). There is no previous report on wave simulation in three-dimensions and three-components. Only through three dimensional numerical simulations can we have an overall understanding of wave field coupling relations and the spatial distribution characteristics between the solid and fluid phases in the dual-phase anisotropic medium. In this paper, based on the BISQ equation, we present elastic wave propagation in a three dimensional dual-phase anisotropic medium simulated by the staggered-grid high-order finite-difference method. We analyze the resulting wave fields and show that the results are an improvement.
Interface traction stress of 3D dislocation loop in anisotropic bimaterial
NASA Astrophysics Data System (ADS)
Wu, Wenwang; Lv, Cunjing; Zhang, JinHuan
2016-02-01
By applying discrete Fast Fourier Transformation (FFT), semi-analytical solutions are developed to calculate the interface elastic fields of anisotropic bimaterial systems with perfect bonding, dislocation-like, force-like and linear spring-like interface models. Interface elastic fields are the linear superimposition of bulk stress, free surface relaxation image stress and interface traction stress (ITS) fields. Interface image energy of perfect bonding bimaterials can be solved through area integral over the interface plane, including the contribution of several componential stress fields. Calculation examples on dislocation loops within Cu-Nb bimaterial are performed to demonstrate the efficiency of such approaches. Effects of Ku = [kiju] for the dislocation-like, Kt = [kijt] for the force-like and Ks = diag [KT , KN ] for the linear spring-like imperfect interface models are investigated. Differences between perfect bonding and imperfect interface models, isotropic and anisotropic models are also studied. It is found that interface conditions and anisotropy have drastic effects on the interface elastic fields.
Enhanced Optical Transmission Mediated by Localized Plasmons in Anisotropic, 3D Nanohole Arrays
Yang, Jiun-Chan; Gao, Hanwei; Suh, Jae Yong; Zhou, Wei; Lee, Min Hyung; Odom, Teri W.
2010-01-01
This paper describes 3D nanohole arrays whose high optical transmission is mediated more by localized surface plasmon (LSP) excitations than by surface plasmon polaritons (SPPs). First, LSPs on 3D hole arrays lead to optical transmission an order of magnitude higher than 2D planar hole arrays. Second, LSP-mediated transmission is broadband and more tunable than SPP-enhanced transmission which is restricted by Bragg coupling. Third, for the first time, two types of surface plasmons can be selectively excited and manipulated on the same plasmonic substrate. This new plasmonic substrate fabricated by high-throughput nanolithography techniques paves the way for cutting-edge optoelectronic and biomedical applications. PMID:20698633
Sonka, Milan; Abramoff, Michael D.
2013-01-01
In this paper, MMSE estimator is employed for noise-free 3D OCT data recovery in 3D complex wavelet domain. Since the proposed distribution for noise-free data plays a key role in the performance of MMSE estimator, a priori distribution for the pdf of noise-free 3D complex wavelet coefficients is proposed which is able to model the main statistical properties of wavelets. We model the coefficients with a mixture of two bivariate Gaussian pdfs with local parameters which are able to capture the heavy-tailed property and inter- and intrascale dependencies of coefficients. In addition, based on the special structure of OCT images, we use an anisotropic windowing procedure for local parameters estimation that results in visual quality improvement. On this base, several OCT despeckling algorithms are obtained based on using Gaussian/two-sided Rayleigh noise distribution and homomorphic/nonhomomorphic model. In order to evaluate the performance of the proposed algorithm, we use 156 selected ROIs from 650 × 512 × 128 OCT dataset in the presence of wet AMD pathology. Our simulations show that the best MMSE estimator using local bivariate mixture prior is for the nonhomomorphic model in the presence of Gaussian noise which results in an improvement of 7.8 ± 1.7 in CNR. PMID:24222760
Interslice interpolation of anisotropic 3D images using multiresolution contour correlation
NASA Astrophysics Data System (ADS)
Lee, Jiann-Der; Wan, Shu-Yen; Ma, Cherng-Min
2002-05-01
To visualize, manipulate and analyze the geometrical structure of anatomical changes, it is often required to perform three-dimensional (3-D) interpolation of the interested organ shape from a series of cross-sectional images obtained from various imaging modalities, such as ultrasound, computed tomography (CT), magnetic resonance imaging (MRI), etc. In this paper, a novel wavelet-based interpolation scheme, which consists of four algorithms are proposed to 3-D image reconstruction. The multi-resolution characteristics of wavelet transform (WT) is completely used in this approach, which consists of two stages, boundary extraction and contour interpolation. More specifically, a wavelet-based radial search method is first designed to extract the boundary of the target object. Next, the global information of the extracted boundary is analyzed for interpolation using WT with various bases and scales. By using six performance measures to evaluate the effectiveness of the proposed scheme, experimental results show that the performance of all proposed algorithms is superior to traditional contour-based methods, linear interpolation and B-spline interpolation. The satisfactory outcome of the proposed scheme provides its capability for serving as an essential part of image processing system developed for medical applications.
Full 3D dispersion curve solutions for guided waves in generally anisotropic media
NASA Astrophysics Data System (ADS)
Hernando Quintanilla, F.; Lowe, M. J. S.; Craster, R. V.
2016-02-01
Dispersion curves of guided waves provide valuable information about the physical and elastic properties of waves propagating within a given waveguide structure. Algorithms to accurately compute these curves are an essential tool for engineers working in non-destructive evaluation and for scientists studying wave phenomena. Dispersion curves are typically computed for low or zero attenuation and presented in two or three dimensional plots. The former do not always provide a clear and complete picture of the dispersion loci and the latter are very difficult to obtain when high values of attenuation are involved and arbitrary anisotropy is considered in single or multi-layered systems. As a consequence, drawing correct and reliable conclusions is a challenging task in the modern applications that often utilize multi-layered anisotropic viscoelastic materials. These challenges are overcome here by using a spectral collocation method (SCM) to robustly find dispersion curves in the most complicated cases of high attenuation and arbitrary anisotropy. Solutions are then plotted in three-dimensional frequency-complex wavenumber space, thus gaining much deeper insight into the nature of these problems. The cases studied range from classical examples, which validate this approach, to new ones involving materials up to the most general triclinic class for both flat and cylindrical geometry in multi-layered systems. The apparent crossing of modes within the same symmetry family in viscoelastic media is also explained and clarified by the results. Finally, the consequences of the centre of symmetry, present in every crystal class, on the solutions are discussed.
NASA Astrophysics Data System (ADS)
Persch, Nico; Elhayek, Ahmed; Welk, Martin; Bruhn, Andrés; Grewenig, Sven; Böse, Katharina; Kraegeloh, Annette; Weickert, Joachim
2013-12-01
This paper proposes an advanced image enhancement method that is specifically tailored towards 3-D confocal and STED microscopy imagery. Our approach unifies image denoising, deblurring and interpolation in one joint method to handle the typical weaknesses of these advanced microscopy techniques: out-of-focus blur, Poisson noise and low axial resolution. In detail, we propose the combination of (i) Richardson-Lucy deconvolution, (ii) image restoration and (iii) anisotropic inpainting in one single scheme. To this end, we develop a novel PDE-based model that realizes these three ideas. First we consider a basic variational image restoration functional that is turned into a joint interpolation scheme by extending the regularization domain. Next, we integrate the variational representation of Richardson-Lucy deconvolution into our model, and illustrate its relation to Poisson distributed noise. In the following step, we supplement the components of our model with sub-quadratic penalization strategies that increase the robustness of the overall method. Finally, we consider the associated minimality conditions, where we exchange the occurring scalar-valued diffusivity function by a so-called diffusion tensor. This leads to an anisotropic regularization that is aligned with structures in the evolving image. As a further contribution of this paper, we propose a more efficient and faster semi-implicit iteration scheme that also increases the stability. Our experiments on real data sets demonstrate that this joint model achieves a superior reconstruction quality of the recorded cell.
Anderson localization of matter waves in 3D anisotropic disordered potentials
NASA Astrophysics Data System (ADS)
Piraud, Marie
2016-05-01
We study quantum transport and Anderson localization of matterwaves in 3 dimensional correlated disorder, focusing on the effects of the anisotropy. Indeed, understanding the anisotropy effects is fundamental for experiments with ultracold atoms as well as for several other systems, such as electrons in MOSFETs, light in biological medium, liquid crystals. A major challenge is to understand whether the anisotropy of the diffusion tensor is altered by the interference terms at the origin of Anderson localization. In particular, its anisotropy at the mobility edge remains to be investigated. So far, all theoretical analysis have assumed - more or less implicitly - that the anisotropy of the diffusion tensor is preserved by interference effects, and have focussed on the vanishing of diffusion as a whole. In this talk, I will start by presenting the usual description of matterwave transport in disordered medium. I will then present our method to go beyond the standard self-consistent theory, which includes in particular the full anisotropic structure of the spectral function. It thus avoids the infrared divergence of the usual self-consistent theory and, most importantly, does not make any assumption on the anisotropy of the renormalized diffusion tensor when including quantum interference terms. Using a generic model of disorder with elongated correlations, we find that the diffusion tensor is strongly affected by the quantum interference terms and that the anisotropy strongly diminishes in the vicinity of the mobility edge. Our work paves the way to further investigation with speckle potentials, which are directly relevant to ultracold-atom experiments. It will permit comparison with previous predictions for the mobility edge and shed new light on ongoing experiments in the field of ultracold atoms.
NASA Astrophysics Data System (ADS)
Berenschot, Erwin J. W.; Jansen, Henri V.; Tas, Niels R.
2013-05-01
When it comes to high-performance filtration, separation, sunlight collection, surface charge storage or catalysis, the effective surface area is what counts. Highly regular fractal structures seem to be the perfect candidates, but manufacturing can be quite cumbersome. Here it is shown--for the first time—that complex 3D fractals can be engineered using a recursive operation in conventional micromachining of single crystalline silicon. The procedure uses the built-in capability of the crystal lattice to form self-similar octahedral structures with minimal interference of the constructor. The silicon fractal can be used directly or as a mold to transfer the shape into another material. Moreover, they can be dense, porous, or like a wireframe. We demonstrate, after four levels of processing, that the initial number of octahedral structures is increased by a factor of 625. Meanwhile the size decreases 16 times down to 300 nm. At any level, pores of less than 100 nm can be fabricated at the octahedral vertices of the fractal. The presented technique supports the design of fractals with Hausdorff dimension D free of choice and up to D = 2.322.
Anisotropic magnification distortion of the 3D galaxy correlation. II. Fourier and redshift space
Hui Lam; Gaztanaga, Enrique; LoVerde, Marilena
2008-03-15
In paper I of this series we discuss how magnification bias distorts the 3D correlation function by enhancing the observed correlation in the line-of-sight (LOS) orientation, especially on large scales. This lensing anisotropy is distinctive, making it possible to separately measure the galaxy-galaxy, galaxy-magnification and magnification-magnification correlations. Here we extend the discussion to the power spectrum and also to redshift space. In real space, pairs oriented close to the LOS direction are not protected against nonlinearity even if the pair separation is large; this is because nonlinear fluctuations can enter through gravitational lensing at a small transverse separation (or i.e. impact parameter). The situation in Fourier space is different: by focusing on a small wave number k, as is usually done, linearity is guaranteed because both the LOS and transverse wave numbers must be small. This is why magnification distortion of the galaxy correlation appears less severe in Fourier space. Nonetheless, the effect is non-negligible, especially for the transverse Fourier modes, and should be taken into account in interpreting precision measurements of the galaxy power spectrum, for instance those that focus on the baryon oscillations. The lensing induced anisotropy of the power spectrum has a shape that is distinct from the more well-known redshift space anisotropies due to peculiar motions and the Alcock-Paczynski effect. The lensing anisotropy is highly localized in Fourier space while redshift space distortions are more spread out. This means that one could separate the magnification bias component in real observations, implying that potentially it is possible to perform a gravitational lensing measurement without measuring galaxy shapes.
Optimizing visual comfort for stereoscopic 3D display based on color-plus-depth signals.
Shao, Feng; Jiang, Qiuping; Fu, Randi; Yu, Mei; Jiang, Gangyi
2016-05-30
Visual comfort is a long-facing problem in stereoscopic 3D (S3D) display. In this paper, targeting to produce S3D content based on color-plus-depth signals, a general framework for depth mapping to optimize visual comfort for S3D display is proposed. The main motivation of this work is to remap the depth range of color-plus-depth signals to a new depth range that is suitable to comfortable S3D display. Towards this end, we first remap the depth range globally based on the adjusted zero disparity plane, and then present a two-stage global and local depth optimization solution to solve the visual comfort problem. The remapped depth map is used to generate the S3D output. We demonstrate the power of our approach on perceptually uncomfortable and comfortable stereoscopic images. PMID:27410090
Metos3D: a marine ecosystem toolkit for optimization and simulation in 3-D - Simulation Package v0.2
NASA Astrophysics Data System (ADS)
Piwonski, J.; Slawig, T.
2015-06-01
A general programming interface for parameter identification for marine ecosystem models is introduced. A comprehensive solver software for periodic steady-states is implemented that includes a fixed point iteration (spin-up) and a Newton solver. The software is based on the Portable, Extensible Toolkit for Scientific Computation (PETSc) library and uses transport matrices for efficient off-line simulation in 3-D. In addition to the usage of PETSc's parallel data structures and PETSc's Newton solver, an own load balancing algorithm is implemented. A simple verification is carried out using a well investigated biogeochemical model for phosphate (PO4) and dissolved organic phosphorous (DOP) with 7 parameters. The model is coupled via the interface to transport matrices that correspond to a longitudinal and latitudinal resolution of 2.8125° and 15 vertical layers. Initial tests show that both solvers and the load balancing algorithm work correctly. Further experiments demonstrate the robustness of the Newton solver with respect to parameter variations. Moreover, the numerical tests reveal that, with optimal control settings, the Newton solver converges at least 6 times faster towards a solution than the spin-up. However, additional twin experiments reveal differences between both solvers regarding a derivative-based black-box optimization. Whereas an optimization run with spin-up-based model evaluations is capable to identify model parameters of a reference solution, Newton-based model evaluations result in an inaccurate gradient approximation.
Optimizing fabrication of electrodeposited 3D surface features
Steffani, C., LLNL
1998-07-01
Selective electrodeposition onto specially masked surfaces (producing blind holes with non-conductive side walls) has been investigated. The thick masks are low cost polymer sheet into which holes and other patterns are machined. Plating with this type of masking system is one method of fabricating medium aspect ratio metallic structures (structures whose width to length ratio is greater than 10/1). The experiment was conducted with deposits of OFC copper from an acid sulfate solution, but similar results are expected with other metals. Structure diameter varied from 0.25 mm to 3.8 mm (0.01`` to 0.150``). The effect of current density, electrolyte concentration, solution temperature, and agitation on deposition rate were investigated. Deposit quality, deposition rate, and optimal plating parameters were evaluated.
3D modeling and optimization of the ITER ICRH antenna
NASA Astrophysics Data System (ADS)
Louche, F.; Dumortier, P.; Durodié, F.; Messiaen, A.; Maggiora, R.; Milanesio, D.
2011-12-01
The prediction of the coupling properties of the ITER ICRH antenna necessitates the accurate evaluation of the resistance and reactance matrices. The latter are mostly dependent on the geometry of the array and therefore a model as accurate as possible is needed to precisely compute these matrices. Furthermore simulations have so far neglected the poloidal and toroidal profile of the plasma, and it is expected that the loading by individual straps will vary significantly due to varying strap-plasma distance. To take this curvature into account, some modifications of the alignment of the straps with respect to the toroidal direction are proposed. It is shown with CST Microwave Studio® [1] that considering two segments in the toroidal direction, i.e. a "V-shaped" toroidal antenna, is sufficient. A new CATIA model including this segmentation has been drawn and imported into both MWS and TOPICA [2] codes. Simulations show a good agreement of the impedance matrices in vacuum. Various modifications of the geometry are proposed in order to further optimize the coupling. In particular we study the effect of the strap box parameters and the recess of the vertical septa.
NASA Astrophysics Data System (ADS)
Ghofrani Tabari, M.; Goodfellow, S. D.; Nasseri, M. B.; Young, R.
2013-12-01
A cubic specimen of water saturated Fontainebleau Sandstone is tested in the laboratory under true-triaxial loading where three different principal stresses are applied under drained conditions. Due to the loading arrangement, closure and opening of the pre-existing cracks in the rock, as well as creation and growth of the aligned cracks cause elliptical anisotropy and distributed heterogeneities. A Geophysical Imaging Cell equipped with an Acoustic Emission monitoring system is employed to image velocity structure of the sample during the experiment through repeated transducer to transducer non-destructive ultrasonic surveys. Apparent P-wave velocities along the rock body are calculated in different directions and shown in stereonet plots which demonstrate an overall anisotropy of the sample. The apparent velocities in the main three orthogonal cubic directions are used as raw data for building a mean spatial distribution model of anisotropy ratios. This approach is based on the concept of semi-principal axes in an elliptical anisotropic model and appointing two ratios between the three orthogonal velocities in each of the cubic grid cells. The spatial distribution model of anisotropy ratios are used to calculate the anisotropic ray-path segment matrix elements (Gij). These contain segment lengths of the ith ray in the jth cell in three dimensions where, length of each ray in each cell is computed for one principal direction based on the dip and strike of the ray and these lengths differ from the ones in an isotropic G Matrix. 3D strain of the squeezed rock and the consequent geometrical deformation is also included in the ray-path segment matrix. A Singular Value Decomposition (SVD) method is used for inversion from the data space of apparent velocities to the model space of P-wave propagation velocities in the three principal directions. Finally, spatial variation and temporal evolution of induced damages in the rock, representing uniformly distributed or
NASA Astrophysics Data System (ADS)
Weinmann, M.; Jutzi, B.; Mallet, C.
2014-08-01
3D scene analysis by automatically assigning 3D points a semantic label has become an issue of major interest in recent years. Whereas the tasks of feature extraction and classification have been in the focus of research, the idea of using only relevant and more distinctive features extracted from optimal 3D neighborhoods has only rarely been addressed in 3D lidar data processing. In this paper, we focus on the interleaved issue of extracting relevant, but not redundant features and increasing their distinctiveness by considering the respective optimal 3D neighborhood of each individual 3D point. We present a new, fully automatic and versatile framework consisting of four successive steps: (i) optimal neighborhood size selection, (ii) feature extraction, (iii) feature selection, and (iv) classification. In a detailed evaluation which involves 5 different neighborhood definitions, 21 features, 6 approaches for feature subset selection and 2 different classifiers, we demonstrate that optimal neighborhoods for individual 3D points significantly improve the results of scene interpretation and that the selection of adequate feature subsets may even further increase the quality of the derived results.
NASA Astrophysics Data System (ADS)
Marcus, Guy
2014-03-01
We have developed a general methodology to experimentally measure the time-resolved Lagrangian orientation and solid body rotation rate of anisotropic particles with arbitrary aspect ratio from standard stereoscopic video image data. We apply these techniques to particles advected in a Rλ ~ 110 fluid flow, where turbulence is generated by two grids oscillating in phase. We use 3D printing technology to design and fabricate neutrally buoyant rods, crosses (two perpendicular rods), and jacks (three mutually perpendicular rods) with a largest dimension of 7 times the Kolmogorov length scale, which makes them good approximations to tracer particles. We have measured the mean square rotation rate, p˙ip˙i , of particles spanning the full range of aspect ratios and obtained results that agree with direct numerical simulations. Our measurements of the full solid-body rotation of jacks, in particular, are of broad experimental relevance because they demonstrate a new and extensible way to directly probe the Lagrangian vorticity of a fluid. Lastly, we will present our direct measurements of the alignment of crosses with the direction of their solid body rotation rate vector, demonstrating how turbulence aligns particles along their longest dimension. This work was completed while at Wesleyan University.
3D Mesh optimization methods for unstructured polyhedra: A progress report
Miller, D.S.; Burton, D.E.
1994-11-22
A mesh optimization scheme allows a Lagrangian code to run problems with extreme mesh distortion by reconfiguring node and zone connectivity as the problem evolves. We have developed some 3D mesh optimization operations and criteria for applying them. These are demonstrated in a 3D Free Lagrange code being developed at LLNL. In the simplest case of a mesh or mesh subregion composed purely of tetrahedra we can maintain a Delaunay tetrahedralization. For more interesting meshes, made up of general polyhedra, a suite of optimization operations and their respective application criteria have been developed.
FPGA Implementation of Optimal 3D-Integer DCT Structure for Video Compression
Jacob, J. Augustin; Kumar, N. Senthil
2015-01-01
A novel optimal structure for implementing 3D-integer discrete cosine transform (DCT) is presented by analyzing various integer approximation methods. The integer set with reduced mean squared error (MSE) and high coding efficiency are considered for implementation in FPGA. The proposed method proves that the least resources are utilized for the integer set that has shorter bit values. Optimal 3D-integer DCT structure is determined by analyzing the MSE, power dissipation, coding efficiency, and hardware complexity of different integer sets. The experimental results reveal that direct method of computing the 3D-integer DCT using the integer set [10, 9, 6, 2, 3, 1, 1] performs better when compared to other integer sets in terms of resource utilization and power dissipation. PMID:26601120
FPGA Implementation of Optimal 3D-Integer DCT Structure for Video Compression.
Jacob, J Augustin; Kumar, N Senthil
2015-01-01
A novel optimal structure for implementing 3D-integer discrete cosine transform (DCT) is presented by analyzing various integer approximation methods. The integer set with reduced mean squared error (MSE) and high coding efficiency are considered for implementation in FPGA. The proposed method proves that the least resources are utilized for the integer set that has shorter bit values. Optimal 3D-integer DCT structure is determined by analyzing the MSE, power dissipation, coding efficiency, and hardware complexity of different integer sets. The experimental results reveal that direct method of computing the 3D-integer DCT using the integer set [10, 9, 6, 2, 3, 1, 1] performs better when compared to other integer sets in terms of resource utilization and power dissipation. PMID:26601120
Meta-Model Based Optimisation Algorithms for Robust Optimization of 3D Forging Sequences
Fourment, Lionel
2007-04-07
In order to handle costly and complex 3D metal forming optimization problems, we develop a new optimization algorithm that allows finding satisfactory solutions within less than 50 iterations (/function evaluation) in the presence of local extrema. It is based on the sequential approximation of the problem objective function by the Meshless Finite Difference Method (MFDM). This changing meta-model allows taking into account the gradient information, if available, or not. It can be easily extended to take into account uncertainties on the optimization parameters. This new algorithm is first evaluated on analytic functions, before being applied to a 3D forging benchmark, the preform tool shape optimization that allows minimizing the potential of fold formation during the two-stepped forging sequence.
Abràmoff, Michael D.; Kardon, Randy; Russell, Stephen R.; Wu, Xiaodong; Sonka, Milan
2008-01-01
Current techniques for segmenting macular optical coherence tomography (OCT) images have been 2-D in nature. Furthermore, commercially available OCT systems have only focused on segmenting a single layer of the retina, even though each intraretinal layer may be affected differently by disease. We report an automated approach for segmenting (anisotropic) 3-D macular OCT scans into five layers. Each macular OCT dataset consisted of six linear radial scans centered at the fovea. The six surfaces defining the five layers were identified on each 3-D composite image by transforming the segmentation task into that of finding a minimum-cost closed set in a geometric graph constructed from edge/regional information and a priori determined surface smoothness and interaction constraints. The method was applied to the macular OCT scans of 12 patients (24 3-D composite image datasets) with unilateral anterior ischemic optic neuropathy (AION). Using the average of three experts’ tracings as a reference standard resulted in an overall mean unsigned border positioning error of 6.1 ± 2.9 µm, a result comparable to the interobserver variability (6.9 ± 3.3 µm). Our quantitative analysis of the automated segmentation results from AION subject data revealed that the inner retinal layer thickness for the affected eye was 24.1 µm (21%) smaller on average than for the unaffected eye (P < 0.001), supporting the need for segmenting the layers separately. PMID:18815101
Topological evolutionary computing in the optimal design of 2D and 3D structures
NASA Astrophysics Data System (ADS)
Burczynski, T.; Poteralski, A.; Szczepanik, M.
2007-10-01
An application of evolutionary algorithms and the finite-element method to the topology optimization of 2D structures (plane stress, bending plates, and shells) and 3D structures is described. The basis of the topological evolutionary optimization is the direct control of the density material distribution (or thickness for 2D structures) by the evolutionary algorithm. The structures are optimized for stress, mass, and compliance criteria. The numerical examples demonstrate that this method is an effective technique for solving problems in computer-aided optimal design.
Liner Optimization Studies Using the Ducted Fan Noise Prediction Code TBIEM3D
NASA Technical Reports Server (NTRS)
Dunn, M. H.; Farassat, F.
1998-01-01
In this paper we demonstrate the usefulness of the ducted fan noise prediction code TBIEM3D as a liner optimization design tool. Boundary conditions on the interior duct wall allow for hard walls or a locally reacting liner with axially segmented, circumferentially uniform impedance. Two liner optimization studies are considered in which farfield noise attenuation due to the presence of a liner is maximized by adjusting the liner impedance. In the first example, the dependence of optimal liner impedance on frequency and liner length is examined. Results show that both the optimal impedance and attenuation levels are significantly influenced by liner length and frequency. In the second example, TBIEM3D is used to compare radiated sound pressure levels between optimal and non-optimal liner cases at conditions designed to simulate take-off. It is shown that significant noise reduction is achieved for most of the sound field by selecting the optimal or near optimal liner impedance. Our results also indicate that there is relatively large region of the impedance plane over which optimal or near optimal liner behavior is attainable. This is an important conclusion for the designer since there are variations in liner characteristics due to manufacturing imprecisions.
An optimal sensing strategy for recognition and localization of 3-D natural quadric objects
NASA Technical Reports Server (NTRS)
Lee, Sukhan; Hahn, Hernsoo
1991-01-01
An optimal sensing strategy for an optical proximity sensor system engaged in the recognition and localization of 3-D natural quadric objects is presented. The optimal sensing strategy consists of the selection of an optimal beam orientation and the determination of an optimal probing plane that compose an optimal data collection operation known as an optimal probing. The decision of an optimal probing is based on the measure of discrimination power of a cluster of surfaces on a multiple interpretation image (MII), where the measure of discrimination power is defined in terms of a utility function computing the expected number of interpretations that can be pruned out by a probing. An object representation suitable for active sensing based on a surface description vector (SDV) distribution graph and hierarchical tables is presented. Experimental results are shown.
Numerical Optimization Strategy for Determining 3D Flow Fields in Microfluidics
NASA Astrophysics Data System (ADS)
Eden, Alex; Sigurdson, Marin; Mezic, Igor; Meinhart, Carl
2015-11-01
We present a hybrid experimental-numerical method for generating 3D flow fields from 2D PIV experimental data. An optimization algorithm is applied to a theory-based simulation of an alternating current electrothermal (ACET) micromixer in conjunction with 2D PIV data to generate an improved representation of 3D steady state flow conditions. These results can be used to investigate mixing phenomena. Experimental conditions were simulated using COMSOL Multiphysics to solve the temperature and velocity fields, as well as the quasi-static electric fields. The governing equations were based on a theoretical model for ac electrothermal flows. A Nelder-Mead optimization algorithm was used to achieve a better fit by minimizing the error between 2D PIV experimental velocity data and numerical simulation results at the measurement plane. By applying this hybrid method, the normalized RMS velocity error between the simulation and experimental results was reduced by more than an order of magnitude. The optimization algorithm altered 3D fluid circulation patterns considerably, providing a more accurate representation of the 3D experimental flow field. This method can be generalized to a wide variety of flow problems. This research was supported by the Institute for Collaborative Biotechnologies through grant W911NF-09-0001 from the U.S. Army Research Office.
Toward 3D-guided prostate biopsy target optimization: an estimation of tumor sampling probabilities
NASA Astrophysics Data System (ADS)
Martin, Peter R.; Cool, Derek W.; Romagnoli, Cesare; Fenster, Aaron; Ward, Aaron D.
2014-03-01
Magnetic resonance imaging (MRI)-targeted, 3D transrectal ultrasound (TRUS)-guided "fusion" prostate biopsy aims to reduce the ~23% false negative rate of clinical 2D TRUS-guided sextant biopsy. Although it has been reported to double the positive yield, MRI-targeted biopsy still yields false negatives. Therefore, we propose optimization of biopsy targeting to meet the clinician's desired tumor sampling probability, optimizing needle targets within each tumor and accounting for uncertainties due to guidance system errors, image registration errors, and irregular tumor shapes. We obtained multiparametric MRI and 3D TRUS images from 49 patients. A radiologist and radiology resident contoured 81 suspicious regions, yielding 3D surfaces that were registered to 3D TRUS. We estimated the probability, P, of obtaining a tumor sample with a single biopsy. Given an RMS needle delivery error of 3.5 mm for a contemporary fusion biopsy system, P >= 95% for 21 out of 81 tumors when the point of optimal sampling probability was targeted. Therefore, more than one biopsy core must be taken from 74% of the tumors to achieve P >= 95% for a biopsy system with an error of 3.5 mm. Our experiments indicated that the effect of error along the needle axis on the percentage of core involvement (and thus the measured tumor burden) was mitigated by the 18 mm core length.
Optimization of a Hybrid Magnetic Bearing for a Magnetically Levitated Blood Pump via 3-D FEA.
Cheng, Shanbao; Olles, Mark W; Burger, Aaron F; Day, Steven W
2011-10-01
In order to improve the performance of a magnetically levitated (maglev) axial flow blood pump, three-dimensional (3-D) finite element analysis (FEA) was used to optimize the design of a hybrid magnetic bearing (HMB). Radial, axial, and current stiffness of multiple design variations of the HMB were calculated using a 3-D FEA package and verified by experimental results. As compared with the original design, the optimized HMB had twice the axial stiffness with the resulting increase of negative radial stiffness partially compensated for by increased current stiffness. Accordingly, the performance of the maglev axial flow blood pump with the optimized HMBs was improved: the maximum pump speed was increased from 6000 rpm to 9000 rpm (50%). The radial, axial and current stiffness of the HMB was found to be linear at nominal operational position from both 3-D FEA and empirical measurements. Stiffness values determined by FEA and empirical measurements agreed well with one another. The magnetic flux density distribution and flux loop of the HMB were also visualized via 3-D FEA which confirms the designers' initial assumption about the function of this HMB. PMID:22065892
Optimization of a Hybrid Magnetic Bearing for a Magnetically Levitated Blood Pump via 3-D FEA
Cheng, Shanbao; Olles, Mark W.; Burger, Aaron F.; Day, Steven W.
2011-01-01
In order to improve the performance of a magnetically levitated (maglev) axial flow blood pump, three-dimensional (3-D) finite element analysis (FEA) was used to optimize the design of a hybrid magnetic bearing (HMB). Radial, axial, and current stiffness of multiple design variations of the HMB were calculated using a 3-D FEA package and verified by experimental results. As compared with the original design, the optimized HMB had twice the axial stiffness with the resulting increase of negative radial stiffness partially compensated for by increased current stiffness. Accordingly, the performance of the maglev axial flow blood pump with the optimized HMBs was improved: the maximum pump speed was increased from 6000 rpm to 9000 rpm (50%). The radial, axial and current stiffness of the HMB was found to be linear at nominal operational position from both 3-D FEA and empirical measurements. Stiffness values determined by FEA and empirical measurements agreed well with one another. The magnetic flux density distribution and flux loop of the HMB were also visualized via 3-D FEA which confirms the designers’ initial assumption about the function of this HMB. PMID:22065892
3D nonrigid registration via optimal mass transport on the GPU.
Ur Rehman, Tauseef; Haber, Eldad; Pryor, Gallagher; Melonakos, John; Tannenbaum, Allen
2009-12-01
In this paper, we present a new computationally efficient numerical scheme for the minimizing flow approach for optimal mass transport (OMT) with applications to non-rigid 3D image registration. The approach utilizes all of the gray-scale data in both images, and the optimal mapping from image A to image B is the inverse of the optimal mapping from B to A. Further, no landmarks need to be specified, and the minimizer of the distance functional involved is unique. Our implementation also employs multigrid, and parallel methodologies on a consumer graphics processing unit (GPU) for fast computation. Although computing the optimal map has been shown to be computationally expensive in the past, we show that our approach is orders of magnitude faster then previous work and is capable of finding transport maps with optimality measures (mean curl) previously unattainable by other works (which directly influences the accuracy of registration). We give results where the algorithm was used to compute non-rigid registrations of 3D synthetic data as well as intra-patient pre-operative and post-operative 3D brain MRI datasets. PMID:19135403
3D nonrigid registration via optimal mass transport on the GPU
Rehman, Tauseef ur; Haber, Eldad; Pryor, Gallagher; Melonakos, John; Tannenbaum, Allen
2009-01-01
In this paper, we present a new computationally efficient numerical scheme for the minimizing flow approach for optimal mass transport (OMT) with applications to non-rigid 3D image registration. The approach utilizes all of the gray-scale data in both images, and the optimal mapping from image A to image B is the inverse of the optimal mapping from B to A. Further, no landmarks need to be specified, and the minimizer of the distance functional involved is unique. Our implementation also employs multigrid, and parallel methodologies on a consumer graphics processing unit (GPU) for fast computation. Although computing the optimal map has been shown to be computationally expensive in the past, we show that our approach is orders of magnitude faster then previous work and is capable of finding transport maps with optimality measures (mean curl) previously unattainable by other works (which directly influences the accuracy of registration). We give results where the algorithm was used to compute non-rigid registrations of 3D synthetic data as well as intra-patient pre-operative and post-operative 3D brain MRI datasets. PMID:19135403
Manifold Based Optimization for Single-Cell 3D Genome Reconstruction
Collas, Philippe
2015-01-01
The three-dimensional (3D) structure of the genome is important for orchestration of gene expression and cell differentiation. While mapping genomes in 3D has for a long time been elusive, recent adaptations of high-throughput sequencing to chromosome conformation capture (3C) techniques, allows for genome-wide structural characterization for the first time. However, reconstruction of "consensus" 3D genomes from 3C-based data is a challenging problem, since the data are aggregated over millions of cells. Recent single-cell adaptations to the 3C-technique, however, allow for non-aggregated structural assessment of genome structure, but data suffer from sparse and noisy interaction sampling. We present a manifold based optimization (MBO) approach for the reconstruction of 3D genome structure from chromosomal contact data. We show that MBO is able to reconstruct 3D structures based on the chromosomal contacts, imposing fewer structural violations than comparable methods. Additionally, MBO is suitable for efficient high-throughput reconstruction of large systems, such as entire genomes, allowing for comparative studies of genomic structure across cell-lines and different species. PMID:26262780
Optimization of 3D MP-RAGE for neonatal brain imaging at 3.0 T.
Williams, Lori-Anne; DeVito, Timothy J; Winter, Jeff D; Orr, Timothy N; Thompson, R Terry; Gelman, Neil
2007-10-01
Three-dimensional (3D) magnetic resonance imaging (MRI) has shown great potential for studying the impact of prematurity and pathology on brain development. We have investigated the potential of optimized T1-weighted 3D magnetization-prepared rapid gradient-echo imaging (MP-RAGE) for obtaining contrast between white matter (WM) and gray matter (GM) in neonates at 3 T. Using numerical simulations, we predicted that the inversion time (TI) for obtaining strongest contrast at 3 T is approximately 2 s for neonates, whereas for adults, this value is approximately 1.3 s. The optimal neonatal TI value was found to be insensitive to reasonable variations of the assumed T1 relaxation times. The maximum theoretical contrast for neonates was found to be approximately one third of that for adults. Using the optimized TI values, MP-RAGE images were obtained from seven neonates and seven adults at 3 T, and the contrast-to-noise ratio (CNR) was measured for WM versus five GM regions. Compared to adults, neonates exhibited lower CNR between cortical GM and WM and showed a different pattern of regional variation in CNR. These results emphasize the importance of sequence optimization specifically for neonates and demonstrate the challenge in obtaining strong contrast in neonatal brain with T1-weighted 3D imaging. PMID:17391887
Bauer, Carl A.; Werner, Gregory R.; Cary, John R.
2011-03-01
A new frequency-domain electromagnetics algorithm is developed for simulating curved interfaces between anisotropic dielectrics embedded in a Yee mesh with second-order error in resonant frequencies. The algorithm is systematically derived using the finite integration formulation of Maxwell's equations on the Yee mesh. Second-order convergence of the error in resonant frequencies is achieved by guaranteeing first-order error on dielectric boundaries and second-order error in bulk (possibly anisotropic) regions. Convergence studies, conducted for an analytically solvable problem and for a photonic crystal of ellipsoids with anisotropic dielectric constant, both show second-order convergence of frequency error; the convergence is sufficiently smooth that Richardson extrapolation yields roughly third-order convergence. The convergence of electric fields near the dielectric interface for the analytic problem is also presented.
Topology optimization of 3D structures with design-dependent loads
NASA Astrophysics Data System (ADS)
Zhang, Hui; Liu, Shu-Tian; Zhang, Xiong
2010-10-01
Topology optimization of continuum structures with design-dependent loads has long been a challenge. In this paper, the topology optimization of 3D structures subjected to design-dependent loads is investigated. A boundary search scheme is proposed for 3D problems, by means of which the load surface can be identified effectively and efficiently, and the difficulties arising in other approaches can be overcome. The load surfaces are made up of the boundaries of finite elements and the loads can be directly applied to corresponding element nodes, which leads to great convenience in the application of this method. Finally, the effectiveness and efficiency of the proposed method is validated by several numerical examples.
Simultaneous Aerodynamic and Structural Design Optimization (SASDO) for a 3-D Wing
NASA Technical Reports Server (NTRS)
Gumbert, Clyde R.; Hou, Gene J.-W.; Newman, Perry A.
2001-01-01
The formulation and implementation of an optimization method called Simultaneous Aerodynamic and Structural Design Optimization (SASDO) is shown as an extension of the Simultaneous Aerodynamic Analysis and Design Optimization (SAADO) method. It is extended by the inclusion of structure element sizing parameters as design variables and Finite Element Method (FEM) analysis responses as constraints. The method aims to reduce the computational expense. incurred in performing shape and sizing optimization using state-of-the-art Computational Fluid Dynamics (CFD) flow analysis, FEM structural analysis and sensitivity analysis tools. SASDO is applied to a simple. isolated, 3-D wing in inviscid flow. Results show that the method finds the saine local optimum as a conventional optimization method with some reduction in the computational cost and without significant modifications; to the analysis tools.
Analysis of the 3D acoustic cloaking problems using optimization method
NASA Astrophysics Data System (ADS)
Alekseev, G. V.; Spivak, Yu E.
2016-06-01
Control problems for the 3D model of acoustic scattering which describes scattering acoustic waves by a permeable obstacle with the form of a spherical layer are considered. These problems arise while developing the design technologies of acoustic cloaking devices using the wave flow method. The solvability of direct and control problems for the acoustic scattering model under study is proved. The sufficient conditions which provide local uniqueness and stability of optimal solutions are established.
NASA Astrophysics Data System (ADS)
Koch, M.; Muench, T.
2010-12-01
There is now ample evidence from both refraction seismic studies and from more recent local earthquake travel-time analysis of some of the authors that large sections of the upper mantle underneath Europe and Germany, in particular, are anisotropic. Employing a modified version of the method of simultaneous inversion for structure and hypocenters (SSH) of the first author, including a priori known upper mantle anisotropy, a full 3D SSH-inversion underneath Germany is carried out. Regional travel times from local events occurring between 1975 - 2003 are used which, after application of several selection criteria, results in ~1300 events with a total of ~30000 P- and S-phases for the SSH inversion. The SSH procedure is carried out in several incremental steps. First of all improved 1D seismic velocity models are derived assuming an isotropic as well as an anisotropic upper mantle. In addition of a slightly better model fit for the anisotropic than for the isotropic model, the latter gives also a somewhat lower Pn-velocity of ~7.90 km/s, compared with ~8.0 km/s for the former. This indicates that inclusion of upper mantle anisotropy into the SSH model is required to obtain physically reasonable Pn-velocities. The results for the P-velocity in the lower crust are less clear, because of some trade-off with the upper mantle layer. Increasingly refined 3D seismic models are then computed, starting with a lateral discretization into 15 x 15 blocs (=40 x 40 km per bloc) and finally going up to 35 × 35 blocs, (=16 x 16 km). For each of the models, inversion solutions for the isotropic, as well as the anisotropic case are examined. The quality of the solution is estimated by means of various tests for resolution, covariance and other trade-off characteristics of the data- and the model-space. Significant improvements for both the isotropic and anisotropic upper mantle cases are obtained for full 3D SSH inversion models. Similar to the 1D Pn-velocity models there are
Simultaneous Aerodynamic Analysis and Design Optimization (SAADO) for a 3-D Flexible Wing
NASA Technical Reports Server (NTRS)
Gumbert, Clyde R.; Hou, Gene J.-W.
2001-01-01
The formulation and implementation of an optimization method called Simultaneous Aerodynamic Analysis and Design Optimization (SAADO) are extended from single discipline analysis (aerodynamics only) to multidisciplinary analysis - in this case, static aero-structural analysis - and applied to a simple 3-D wing problem. The method aims to reduce the computational expense incurred in performing shape optimization using state-of-the-art Computational Fluid Dynamics (CFD) flow analysis, Finite Element Method (FEM) structural analysis and sensitivity analysis tools. Results for this small problem show that the method reaches the same local optimum as conventional optimization. However, unlike its application to the win,, (single discipline analysis), the method. as I implemented here, may not show significant reduction in the computational cost. Similar reductions were seen in the two-design-variable (DV) problem results but not in the 8-DV results given here.
Modelling 3D control of upright stance using an optimal control strategy.
Qu, Xingda; Nussbaum, Maury A
2012-01-01
A 3D balance control model of quiet upright stance is presented, based on an optimal control strategy, and evaluated in terms of its ability to simulate postural sway in both the anterior-posterior and medial-lateral directions. The human body was represented as a two-segment inverted pendulum. Several assumptions were made to linearise body dynamics, for example, that there was no transverse rotation during upright stance. The neural controller was presumed to be an optimal controller that generates ankle control torque and hip control torque according to certain performance criteria. An optimisation procedure was used to determine the values of unspecified model parameters including random disturbance gains and sensory delay times. This model was used to simulate postural sway behaviours characterised by centre-of-pressure (COP)-based measures. Confidence intervals for all normalised COP-based measures contained unity, indicating no significant differences between any of the simulated COP-based measures and corresponding experimental references. In addition, mean normalised errors for the traditional measures were < 8%, and those for most statistical mechanics measures were ∼3-66%. On the basis these results, the proposed 3D balance control model appears to have the ability to accurately simulate 3D postural sway behaviours. PMID:21598131
Optimal Image Stitching for Concrete Bridge Bottom Surfaces Aided by 3d Structure Lines
NASA Astrophysics Data System (ADS)
Liu, Yahui; Yao, Jian; Liu, Kang; Lu, Xiaohu; Xia, Menghan
2016-06-01
Crack detection for bridge bottom surfaces via remote sensing techniques is undergoing a revolution in the last few years. For such applications, a large amount of images, acquired with high-resolution industrial cameras close to the bottom surfaces with some mobile platform, are required to be stitched into a wide-view single composite image. The conventional idea of stitching a panorama with the affine model or the homographic model always suffers a series of serious problems due to poor texture and out-of-focus blurring introduced by depth of field. In this paper, we present a novel method to seamlessly stitch these images aided by 3D structure lines of bridge bottom surfaces, which are extracted from 3D camera data. First, we propose to initially align each image in geometry based on its rough position and orientation acquired with both a laser range finder (LRF) and a high-precision incremental encoder, and these images are divided into several groups with the rough position and orientation data. Secondly, the 3D structure lines of bridge bottom surfaces are extracted from the 3D cloud points acquired with 3D cameras, which impose additional strong constraints on geometrical alignment of structure lines in adjacent images to perform a position and orientation optimization in each group to increase the local consistency. Thirdly, a homographic refinement between groups is applied to increase the global consistency. Finally, we apply a multi-band blending algorithm to generate a large-view single composite image as seamlessly as possible, which greatly eliminates both the luminance differences and the color deviations between images and further conceals image parallax. Experimental results on a set of representative images acquired from real bridge bottom surfaces illustrate the superiority of our proposed approaches.
NASA Astrophysics Data System (ADS)
Chang, Yau-Zen; Hou, Jung-Fu; Tsao, Yi Hsiang; Lee, Shih-Tseng
2011-09-01
This paper proposes a scheme for finding the correspondence between uniformly spaced locations on the images of human face captured from different viewpoints at the same instant. The correspondence is dedicated for 3D reconstruction to be used in the registration procedure for neurosurgery where the exposure to projectors must be seriously restricted. The approach utilizes structured light to enhance patterns on the images and is initialized with the scale-invariant feature transform (SIFT). Successive locations are found according to spatial order using a parallel version of the particle swarm optimization algorithm. Furthermore, false locations are singled out for correction by searching for outliers from fitted curves. Case studies show that the scheme is able to correctly generate 456 evenly spaced 3D coordinate points in 23 seconds from a single shot of projected human face using a PC with 2.66 GHz Intel Q9400 CPU and 4GB RAM.
3D FEM Geometry and Material Flow Optimization of Porthole-Die Extrusion
Ceretti, Elisabetta; Mazzoni, Luca; Giardini, Claudio
2007-05-17
The aim of this work is to design and to improve the geometry of a porthole-die for the production of aluminum components by means of 3D FEM simulations. In fact, the use of finite element models will allow to investigate the effects of the die geometry (webs, extrusion cavity) on the material flow and on the stresses acting on the die so to reduce the die wear and to improve the tool life. The software used to perform the simulations was a commercial FEM code, Deform 3D. The technological data introduced in the FE model have been furnished by METRA S.p.A. Company, partner in this research. The results obtained have been considered valid and helpful by the Company for building a new optimized extrusion porthole-die.
OPTIMIZATION OF 3-D IMAGE-GUIDED NEAR INFRARED SPECTROSCOPY USING BOUNDARY ELEMENT METHOD
Srinivasan, Subhadra; Carpenter, Colin; Pogue, Brian W.; Paulsen, Keith D.
2010-01-01
Multimodality imaging systems combining optical techniques with MRI/CT provide high-resolution functional characterization of tissue by imaging molecular and vascular biomarkers. To optimize these hybrid systems for clinical use, faster and automatable algorithms are required for 3-D imaging. Towards this end, a boundary element model was used to incorporate tissue boundaries from MRI/CT into image formation process. This method uses surface rendering to describe light propagation in 3-D using diffusion equation. Parallel computing provided speedup of up to 54% in time of computation. Simulations showed that location of NIRS probe was crucial for quantitatively accurate estimation of tumor response. A change of up to 61% was seen between cycles 1 and 3 in monitoring tissue response to neoadjuvant chemotherapy. PMID:20523751
Sutradhar, Alok; Park, Jaejong; Carrau, Diana; Nguyen, Tam H; Miller, Michael J; Paulino, Glaucio H
2016-07-01
Large craniofacial defects require efficient bone replacements which should not only provide good aesthetics but also possess stable structural function. The proposed work uses a novel multiresolution topology optimization method to achieve the task. Using a compliance minimization objective, patient-specific bone replacement shapes can be designed for different clinical cases that ensure revival of efficient load transfer mechanisms in the mid-face. In this work, four clinical cases are introduced and their respective patient-specific designs are obtained using the proposed method. The optimized designs are then virtually inserted into the defect to visually inspect the viability of the design . Further, once the design is verified by the reconstructive surgeon, prototypes are fabricated using a 3D printer for validation. The robustness of the designs are mechanically tested by subjecting them to a physiological loading condition which mimics the masticatory activity. The full-field strain result through 3D image correlation and the finite element analysis implies that the solution can survive the maximum mastication of 120 lb. Also, the designs have the potential to restore the buttress system and provide the structural integrity. Using the topology optimization framework in designing the bone replacement shapes would deliver surgeons new alternatives for rather complicated mid-face reconstruction. PMID:26660897
Defining an optimal surface chemistry for pluripotent stem cell culture in 2D and 3D
NASA Astrophysics Data System (ADS)
Zonca, Michael R., Jr.
Surface chemistry is critical for growing pluripotent stem cells in an undifferentiated state. There is great potential to engineer the surface chemistry at the nanoscale level to regulate stem cell adhesion. However, the challenge is to identify the optimal surface chemistry of the substrata for ES cell attachment and maintenance. Using a high-throughput polymerization and screening platform, a chemically defined, synthetic polymer grafted coating that supports strong attachment and high expansion capacity of pluripotent stem cells has been discovered using mouse embryonic stem (ES) cells as a model system. This optimal substrate, N-[3-(Dimethylamino)propyl] methacrylamide (DMAPMA) that is grafted on 2D synthetic poly(ether sulfone) (PES) membrane, sustains the self-renewal of ES cells (up to 7 passages). DMAPMA supports cell attachment of ES cells through integrin beta1 in a RGD-independent manner and is similar to another recently reported polymer surface. Next, DMAPMA has been able to be transferred to 3D by grafting to synthetic, polymeric, PES fibrous matrices through both photo-induced and plasma-induced polymerization. These 3D modified fibers exhibited higher cell proliferation and greater expression of pluripotency markers of mouse ES cells than 2D PES membranes. Our results indicated that desirable surfaces in 2D can be scaled to 3D and that both surface chemistry and structural dimension strongly influence the growth and differentiation of pluripotent stem cells. Lastly, the feasibility of incorporating DMAPMA into a widely used natural polymer, alginate, has been tested. Novel adhesive alginate hydrogels have been successfully synthesized by either direct polymerization of DMAPMA and methacrylic acid blended with alginate, or photo-induced DMAPMA polymerization on alginate nanofibrous hydrogels. In particular, DMAPMA-coated alginate hydrogels support strong ES cell attachment, exhibiting a concentration dependency of DMAPMA. This research provides a
NASA Astrophysics Data System (ADS)
Nagarajan, Mahesh B.; De, Titas; Lochmüller, Eva-Maria; Eckstein, Felix; Wismüller, Axel
2014-04-01
The ability of Anisotropic Minkowski Functionals (AMFs) to capture local anisotropy while evaluating topological properties of the underlying gray-level structures has been previously demonstrated. We evaluate the ability of this approach to characterize local structure properties of trabecular bone micro-architecture in ex vivo proximal femur specimens, as visualized on multi-detector CT, for purposes of biomechanical bone strength prediction. To this end, volumetric AMFs were computed locally for each voxel of volumes of interest (VOI) extracted from the femoral head of 146 specimens. The local anisotropy captured by such AMFs was quantified using a fractional anisotropy measure; the magnitude and direction of anisotropy at every pixel was stored in histograms that served as a feature vectors that characterized the VOIs. A linear multi-regression analysis algorithm was used to predict the failure load (FL) from the feature sets; the predicted FL was compared to the true FL determined through biomechanical testing. The prediction performance was measured by the root mean square error (RMSE) for each feature set. The best prediction performance was obtained from the fractional anisotropy histogram of AMF Euler Characteristic (RMSE = 1.01 ± 0.13), which was significantly better than MDCT-derived mean BMD (RMSE = 1.12 ± 0.16, p<0.05). We conclude that such anisotropic Minkowski Functionals can capture valuable information regarding regional trabecular bone quality and contribute to improved bone strength prediction, which is important for improving the clinical assessment of osteoporotic fracture risk.
NASA Astrophysics Data System (ADS)
Yuan, X. H.; Zhao, G. P.; Yue, Ming; Ye, L. N.; Xia, J.; Zhang, X. C.; Chang, J.
2013-10-01
In this paper, the magnetic reversal process, hysteresis loops and energy products for exchange-coupled Nd2Fe14B/α-Fe bilayers are studied systematically by a three-dimensional (3D) model. The 3D calculations are numerically solved using the finite difference method, where the results are carefully compared with those calculated by one-dimensional (1D) model. It is found that the calculated hysteresis loops and energy products based on the two methods are consistent with each other. Both nucleation fields and coercivities decrease monotonically as the soft layer thickness Ls increases. In addition, the calculated spatial distributions of magnetization orientations in the thickness direction at various applied fields based on both methods signify a three-step magnetic reversal process, which are nucleation, growth and displacement of the domain wall. The calculated magnetic orientations within the film plane, however, are totally different according to the two methods. The 3D calculation exhibits a process of vortex formation and annihilation. On the other hand, the 1D calculation gives a quasi-coherent one, where magnetization orientation is coherent in the film plane and varies in the thickness direction. This new reversal mechanism displayed in the film plane has a systematic influence on the nucleation fields, coercivity and energy products.
NASA Astrophysics Data System (ADS)
Bai, Chao-Ying; Huang, Guo-Jiao; Li, Xiao-Ling; Zhou, Bing; Greenhalgh, Stewart
2013-11-01
To overcome the deficiency of some current grid-/cell-based ray tracing algorithms, which are only able to handle first arrivals or primary reflections (or conversions) in anisotropic media, we have extended the functionality of the multistage irregular shortest-path method to 2-D/3-D tilted transversely isotropic (TTI) media. The new approach is able to track multiple transmitted/reflected/converted arrivals composed of any kind of combinations of transmissions, reflections and mode conversions. The basic principle is that the seven parameters (five elastic parameters plus two polar angles defining the tilt of the symmetry axis) of the TTI media are sampled at primary nodes, and the group velocity values at secondary nodes are obtained by tri-linear interpolation of the primary nodes across each cell, from which the group velocities of the three wave modes (qP, qSV and qSH) are calculated. Finally, we conduct grid-/cell-based wave front expansion to trace multiple transmitted/reflected/converted arrivals from one region to the next. The results of calculations in uniform anisotropic media indicate that the numerical results agree with the analytical solutions except in directions of SV-wave triplications, at which only the lowest velocity value is selected at the singularity points by the multistage irregular shortest-path anisotropic ray tracing method. This verifies the accuracy of the methodology. Several simulation results show that the new method is able to efficiently and accurately approximate situations involving continuous velocity variations and undulating discontinuities, and that it is suitable for any combination of multiple transmitted/reflected/converted arrival tracking in TTI media of arbitrary strength and tilt. Crosshole synthetic traveltime tomographic tests have been performed, which highlight the importance of using such code when the medium is distinctly anisotropic.
Hedegaard, R.F.; Ho, J.; Eisert, J.
1996-12-31
Three-dimensional (3-D) geoscience volume modeling can be used to improve the efficiency of the environmental investigation and remediation process. At several unsaturated zone spill sites at two Superfund (CERCLA) sites (Military Installations) in California, all aspects of subsurface contamination have been characterized using an integrated computerized approach. With the aide of software such as LYNX GMS{trademark}, Wavefront`s Data Visualizer{trademark} and Gstools (public domain), the authors have created a central platform from which to map a contaminant plume, visualize the same plume three-dimensionally, and calculate volumes of contaminated soil or groundwater above important health risk thresholds. The developed methodology allows rapid data inspection for decisions such that the characterization process and remedial action design are optimized. By using the 3-D geoscience modeling and visualization techniques, the technical staff are able to evaluate the completeness and spatial variability of the data and conduct 3-D geostatistical predictions of contaminant and lithologic distributions. The geometry of each plume is estimated using 3-D variography on raw analyte values and indicator thresholds for the kriged model. Three-dimensional lithologic interpretation is based on either {open_quote}linked{close_quote} parallel cross sections or on kriged grid estimations derived from borehole data coded with permeability indicator thresholds. Investigative borings, as well as soil vapor extraction/injection wells, are sighted and excavation costs are estimated using these results. The principal advantages of the technique are the efficiency and rapidity with which meaningful results are obtained and the enhanced visualization capability which is a desirable medium to communicate with both the technical staff as well as nontechnical audiences.
Optimizing illumination in the greenhouse using a 3D model of tomato and a ray tracer.
de Visser, Pieter H B; Buck-Sorlin, Gerhard H; van der Heijden, Gerie W A M
2014-01-01
Reduction of energy use for assimilation lighting is one of the most urgent goals of current greenhouse horticulture in the Netherlands. In recent years numerous lighting systems have been tested in greenhouses, yet their efficiency has been very difficult to measure in practice. This simulation study evaluated a number of lighting strategies using a 3D light model for natural and artificial light in combination with a 3D model of tomato. The modeling platform GroIMP was used for the simulation study. The crop was represented by 3D virtual plants of tomato with fixed architecture. Detailed data on greenhouse architecture and lamp emission patterns of different light sources were incorporated in the model. A number of illumination strategies were modeled with the calibrated model. Results were compared to the standard configuration. Moreover, adaptation of leaf angles was incorporated for testing their effect on light use efficiency (LUE). A Farquhar photosynthesis model was used to translate the absorbed light for each leaf into a produced amount of carbohydrates. The carbohydrates produced by the crop per unit emitted light from sun or high pressure sodium lamps was the highest for horizontal leaf angles or slightly downward pointing leaves, and was less for more upward leaf orientations. The simulated leaf angles did not affect light absorption from inter-lighting LED modules, but the scenario with LEDs shining slightly upward (20(°)) increased light absorption and LUE relative to default horizontal beaming LEDs. Furthermore, the model showed that leaf orientation more perpendicular to the string of LEDs increased LED light interception. The combination of a ray tracer and a 3D crop model could compute optimal lighting of leaves by quantification of light fluxes and illustration by rendered lighting patterns. Results indicate that illumination efficiency increases when the lamp light is directed at most to leaves that have a high photosynthetic potential. PMID
Towards 3D lidar point cloud registration improvement using optimal neighborhood knowledge
NASA Astrophysics Data System (ADS)
Gressin, Adrien; Mallet, Clément; Demantké, Jérôme; David, Nicolas
2013-05-01
Automatic 3D point cloud registration is a main issue in computer vision and remote sensing. One of the most commonly adopted solution is the well-known Iterative Closest Point (ICP) algorithm. This standard approach performs a fine registration of two overlapping point clouds by iteratively estimating the transformation parameters, assuming good a priori alignment is provided. A large body of literature has proposed many variations in order to improve each step of the process (namely selecting, matching, rejecting, weighting and minimizing). The aim of this paper is to demonstrate how the knowledge of the shape that best fits the local geometry of each 3D point neighborhood can improve the speed and the accuracy of each of these steps. First we present the geometrical features that form the basis of this work. These low-level attributes indeed describe the neighborhood shape around each 3D point. They allow to retrieve the optimal size to analyze the neighborhoods at various scales as well as the privileged local dimension (linear, planar, or volumetric). Several variations of each step of the ICP process are then proposed and analyzed by introducing these features. Such variants are compared on real datasets with the original algorithm in order to retrieve the most efficient algorithm for the whole process. Therefore, the method is successfully applied to various 3D lidar point clouds from airborne, terrestrial, and mobile mapping systems. Improvement for two ICP steps has been noted, and we conclude that our features may not be relevant for very dissimilar object samplings.
Optimizing illumination in the greenhouse using a 3D model of tomato and a ray tracer
de Visser, Pieter H. B.; Buck-Sorlin, Gerhard H.; van der Heijden, Gerie W. A. M.
2014-01-01
Reduction of energy use for assimilation lighting is one of the most urgent goals of current greenhouse horticulture in the Netherlands. In recent years numerous lighting systems have been tested in greenhouses, yet their efficiency has been very difficult to measure in practice. This simulation study evaluated a number of lighting strategies using a 3D light model for natural and artificial light in combination with a 3D model of tomato. The modeling platform GroIMP was used for the simulation study. The crop was represented by 3D virtual plants of tomato with fixed architecture. Detailed data on greenhouse architecture and lamp emission patterns of different light sources were incorporated in the model. A number of illumination strategies were modeled with the calibrated model. Results were compared to the standard configuration. Moreover, adaptation of leaf angles was incorporated for testing their effect on light use efficiency (LUE). A Farquhar photosynthesis model was used to translate the absorbed light for each leaf into a produced amount of carbohydrates. The carbohydrates produced by the crop per unit emitted light from sun or high pressure sodium lamps was the highest for horizontal leaf angles or slightly downward pointing leaves, and was less for more upward leaf orientations. The simulated leaf angles did not affect light absorption from inter-lighting LED modules, but the scenario with LEDs shining slightly upward (20°) increased light absorption and LUE relative to default horizontal beaming LEDs. Furthermore, the model showed that leaf orientation more perpendicular to the string of LEDs increased LED light interception. The combination of a ray tracer and a 3D crop model could compute optimal lighting of leaves by quantification of light fluxes and illustration by rendered lighting patterns. Results indicate that illumination efficiency increases when the lamp light is directed at most to leaves that have a high photosynthetic potential. PMID
Bao, Qian; Jiang, Chenglong; Lin, Yun; Tan, Weixian; Wang, Zhirui; Hong, Wen
2016-01-01
With a short linear array configured in the cross-track direction, downward looking sparse linear array three-dimensional synthetic aperture radar (DLSLA 3-D SAR) can obtain the 3-D image of an imaging scene. To improve the cross-track resolution, sparse recovery methods have been investigated in recent years. In the compressive sensing (CS) framework, the reconstruction performance depends on the property of measurement matrix. This paper concerns the technique to optimize the measurement matrix and deal with the mismatch problem of measurement matrix caused by the off-grid scatterers. In the model of cross-track reconstruction, the measurement matrix is mainly affected by the configuration of antenna phase centers (APC), thus, two mutual coherence based criteria are proposed to optimize the configuration of APCs. On the other hand, to compensate the mismatch problem of the measurement matrix, the sparse Bayesian inference based method is introduced into the cross-track reconstruction by jointly estimate the scatterers and the off-grid error. Experiments demonstrate the performance of the proposed APCs' configuration schemes and the proposed cross-track reconstruction method. PMID:27556471
A multi-objective optimization framework to model 3D river and landscape evolution processes
NASA Astrophysics Data System (ADS)
Bizzi, Simone; Castelletti, Andrea; Cominola, Andrea; Mason, Emanuele; Paik, Kyungrock
2013-04-01
Water and sediment interactions shape hillslopes, regulate soil erosion and sedimentation, and organize river networks. Landscape evolution and river organization occur at various spatial and temporal scale and the understanding and modelling of them is highly complex. The idea of a least action principle governing river networks evolution has been proposed many times as a simpler approach among the ones existing in the literature. These theories assume that river networks, as observed in nature, self-organize and act on soil transportation in order to satisfy a particular "optimality" criterion. Accordingly, river and landscape weathering can be simulated by solving an optimization problem, where the choice of the criterion to be optimized becomes the initial assumption. The comparison between natural river networks and optimized ones verifies the correctness of this initial assumption. Yet, various criteria have been proposed in literature and there is no consensus on which is better able to explain river network features observed in nature like network branching and river bed profile: each one is able to reproduce some river features through simplified modelling of the natural processes, but it fails to characterize the whole complexity (3D and its dynamic) of the natural processes. Some of the criteria formulated in the literature partly conflict: the reason is that their formulation rely on mathematical and theoretical simplifications of the natural system that are suitable for specific spatial and temporal scale but fails to represent the whole processes characterizing landscape evolution. In an attempt to address some of these scientific questions, we tested the suitability of using a multi-objective optimization framework to describe river and landscape evolution in a 3D spatial domain. A synthetic landscape is used to this purpose. Multiple, alternative river network evolutions, corresponding to as many tradeoffs between the different and partly
Torque-consistent 3D force balance and optimization of non-resonant fields in tokamaks
NASA Astrophysics Data System (ADS)
Park, Jong-Kyu
2015-11-01
A non-axisymmetric magnetic perturbation in tokamaks breaks the toroidal symmetry and produces toroidal torque, which is well known as neoclassical toroidal viscosity (NTV) effects. Although NTV torque is second order, it is the first-order change in the pressure anisotropy that drives currents associated with local torques and thereby modifies the field penetration in force balance. The force operator becomes non-Hermitian, but can be directly solved using parallel, toroidal, and radial force balance, leading to a modified Euler-Lagrange equation. The general perturbed equilibrium code (GPEC), which has been successfully developed to solve the modified Euler-Lagrange equation, gives the torque-consistent 3D force balance as well as self-consistent NTV torque. The self-shielding of the torque becomes apparent in the solutions in high β, which was implied in recent MARS-K applications. Furthermore, the full response matrix including the torque in GPEC provides a new and systematic way of optimizing torque and non-resonant fields. Recently the optimization of 3D fields for torque has been actively studied using the stellarator optimizing tools, but the efficiency and accuracy can be greatly improved by directly incorporating the torque response matrix. There are salient features uncovered by response with the torque, as the response can become invisible in amplitudes but only significant in toroidal phase shift. A perturbation in backward helicity is an example, in which NTV can be induced substantially but quietly without measurable response in amplitudes. A number of other GPEC applications will also be discussed, including the multi-mode responses in high- β tokamak plasmas and the new non-axisymmetric control coil (NCC) design in NSTX-U. This work was supported by DOE Contract DE-AC02-09CH11466.
Optimizing penetration depth, contrast, and resolution in 3D dermatologic OCT
NASA Astrophysics Data System (ADS)
Aneesh, Alex; Považay, Boris; Hofer, Bernd; Zhang, Edward Z.; Kendall, Catherine; Laufer, Jan; Popov, Sergei; Glittenberg, Carl; Binder, Susanne; Stone, Nicholas; Beard, Paul C.; Drexler, Wolfgang
2010-02-01
High speed, three-dimensional optical coherence tomography (3D OCT) at 800nm, 1060nm and 1300nm with approximately 4μm, 7μm and 6μm axial and less than 15μm transverse resolution is demonstrated to investigate the optimum wavelength region for in vivo human skin imaging in terms of contrast, dynamic range and penetration depth. 3D OCT at 1300nm provides deeper penetration, while images obtained at 800nm were better in terms of contrast and speckle noise. 1060nm region was a compromise between 800nm and 1300nm in terms of penetration depth and image contrast. Optimizing sensitivity, penetration and contrast enabled unprecedented visualization of micro-structural morphology underneath the glabrous skin, hairy skin and in scar tissue. Higher contrast obtained at 800 nm appears to be critical in the in vitro tumor study. A multimodal approach combining OCT and PA helped to obtain morphological as well as vascular information from deeper regions of skin.
NASA Astrophysics Data System (ADS)
Wang, Shuai; Wang, Yu; Zi, Yanyang; He, Zhengjia
2015-12-01
A generalized and efficient model for rotating anisotropic rotor-bearing systems is presented in this paper with full considerations of the system's anisotropy in stiffness, inertia and damping. Based on the 3D finite element model and the model order reduction method, the effects of anisotropy in shaft and bearings on the forced response and whirling of anisotropic rotor-bearing systems are systematically investigated. First, the coefficients of journal bearings are transformed from the fixed frame to the rotating one. Due to the anisotropy in shaft and bearings, the motion is governed by differential equations with periodically time-variant coefficients. Then, a free-interface complex component mode synthesis (CMS) method is employed to generate efficient reduced-order models (ROM) for the periodically time-variant systems. In order to solve the obtained equations, a variant of Hill's method for systems with multiple harmonic excitations is developed. Four dimensionless parameters are defined to quantify the types and levels of anisotropy of bearings. Finally, the effects of the four types of anisotropy on the forced response and whirl orbits are studied. Numerical results show that the anisotropy of bearings in stiffness splits the sole resonant peak into two isolated ones, but the anisotropy of bearings in damping coefficients mainly affect the response amplitudes. Moreover, the whirl orbits become much more complex when the shaft and bearings are both anisotropic. In addition, the cross-coupling stiffness coefficients of bearings significantly affect the dynamic behaviors of the systems and cannot be neglected, though they are often much smaller than the principle stiffness terms.
External fixator configurations in tibia fractures: 1D optimization and 3D analysis comparison.
Roseiro, Luis M; Neto, M Augusta; Amaro, Ana; Leal, Rogerio P; Samarra, Miguel C
2014-01-01
The use of external fixation devices in orthopedic surgery is very common in open tibial fractures. A properly applied fixator may improve the healing process while one improperly applied might delay the healing process. The several external fixator systems used in clinical today, can be categorized into uniplanar-unilateral, uniplanar-bilateral, biplanar and multiplanar. The stability on the fracture focus and, therefore, the fracture healing process, is related with the type of external fixator configuration that is selected. The aim of this study is to discuss the principles for the successful application of unilateral-uniplanar external fixation, the assembly of its components, for the case of a transverse fractures using computational models. In this context, the fixation stiffness characteristics are evaluated using a simplified 1D finite element model for the tibia and external fixator. The beams are modeled with realistic cross-sectional geometry and material properties instead of a simplified model. The VABS (the Variational Asymptotic Beam Section analysis) methodology is used to compute the cross-sectional model for the generalized Timoshenko model, which was embedded in the finite element solver FEAP. The use of Timoshenko beam theory allows accounting for several kinds of loads, including torsion moments. Optimal design is performed with respect to the assembly of fixator components using a genetic algorithm. The optimization procedure is based on the evaluation of an objective function, which is dependent on the displacement at the fracture focus. The initial and optimal results are compared by performing a 3D analysis, for which different three-dimensional finite element models are created. The geometrical model of a tibia is created on the basis of data acquired by CAT scan, made for a healthy tibia of a 22 year old male. The 3D comparison of the 1D optimal results show a clear improvement on the objective function for the several load cases and
3-D Anisotropic Ambient Noise Tomography of Piton De La Fournaise Volcano (La Réunion Island)
NASA Astrophysics Data System (ADS)
Mordret, A.; Rivet, D. N.; Landes, M.; Shapiro, N.
2014-12-01
We cross-correlate four years of seismic noise continuously recorded by the seismic monitoring network of the Piton de la Fournaise volcano (La Réunion Island). The network is composed of 40 stations 27 of which have 3-component sensors. We use Vertical-to-Vertical (ZZ) cross-correlation components from all stations and Radial-to-Radial (RR) and Transverse-to-Transverse (TT) cross-correlations computed from 3-component records. The group velocity dispersion curves for Rayleigh and Love waves are measured using a Frequency-Time Analysis. We average measurements from ZZ and RR components to finally obtain 577 Rayleigh-wave dispersion curves. 395 Love-wave dispersion curves are obtained from the TT cross-correlations. We then regionalize the group velocities measurements to construct 2D dispersion maps at a set of periods between 0.4 and 8 s. Finally, we construct a 3D shear-velocity model down to 3 km below the sea level by jointly inverting the Rayleigh and Love wave group velocity maps with a Neighborhood Algorithm and with taking into account the radial anisotropy. The distribution of 3-D Voigt averaged S-wave velocities shows three distinct high-velocity anomalies surrounded by a low-velocity ring. The most western high-velocity anomaly is located below the actual "Plaine des Sables" and could be attributed to an old intrusive body at the location of the former volcanic center before it migrated toward its present location. The second high-velocity body is located below the summit of the volcano and likely corresponds to the actual preferential dyke intrusion zone as highlighted by the seismicity. The third high-velocity anomaly is located below the "Grandes Pentes" and the "Grand Brûlé" areas and is thought to be an imprint of the solidified magma chamber of the ancient dismantled "Les Alizé" volcano. The distribution of the radial anisotropy shows two main anomalies: a positive anisotropy (Vsh>Vsv) above sea level highlighting the recent edifice of Piton de
Structural optimization of 3D-printed synthetic spider webs for high strength
NASA Astrophysics Data System (ADS)
Qin, Zhao; Compton, Brett G.; Lewis, Jennifer A.; Buehler, Markus J.
2015-05-01
Spiders spin intricate webs that serve as sophisticated prey-trapping architectures that simultaneously exhibit high strength, elasticity and graceful failure. To determine how web mechanics are controlled by their topological design and material distribution, here we create spider-web mimics composed of elastomeric filaments. Specifically, computational modelling and microscale 3D printing are combined to investigate the mechanical response of elastomeric webs under multiple loading conditions. We find the existence of an asymptotic prey size that leads to a saturated web strength. We identify pathways to design elastomeric material structures with maximum strength, low density and adaptability. We show that the loading type dictates the optimal material distribution, that is, a homogeneous distribution is better for localized loading, while stronger radial threads with weaker spiral threads is better for distributed loading. Our observations reveal that the material distribution within spider webs is dictated by the loading condition, shedding light on their observed architectural variations.
Structural optimization of 3D-printed synthetic spider webs for high strength.
Qin, Zhao; Compton, Brett G; Lewis, Jennifer A; Buehler, Markus J
2015-01-01
Spiders spin intricate webs that serve as sophisticated prey-trapping architectures that simultaneously exhibit high strength, elasticity and graceful failure. To determine how web mechanics are controlled by their topological design and material distribution, here we create spider-web mimics composed of elastomeric filaments. Specifically, computational modelling and microscale 3D printing are combined to investigate the mechanical response of elastomeric webs under multiple loading conditions. We find the existence of an asymptotic prey size that leads to a saturated web strength. We identify pathways to design elastomeric material structures with maximum strength, low density and adaptability. We show that the loading type dictates the optimal material distribution, that is, a homogeneous distribution is better for localized loading, while stronger radial threads with weaker spiral threads is better for distributed loading. Our observations reveal that the material distribution within spider webs is dictated by the loading condition, shedding light on their observed architectural variations. PMID:25975372
3D transrectal ultrasound (TRUS) prostate segmentation based on optimal feature learning framework
NASA Astrophysics Data System (ADS)
Yang, Xiaofeng; Rossi, Peter J.; Jani, Ashesh B.; Mao, Hui; Curran, Walter J.; Liu, Tian
2016-03-01
We propose a 3D prostate segmentation method for transrectal ultrasound (TRUS) images, which is based on patch-based feature learning framework. Patient-specific anatomical features are extracted from aligned training images and adopted as signatures for each voxel. The most robust and informative features are identified by the feature selection process to train the kernel support vector machine (KSVM). The well-trained SVM was used to localize the prostate of the new patient. Our segmentation technique was validated with a clinical study of 10 patients. The accuracy of our approach was assessed using the manual segmentations (gold standard). The mean volume Dice overlap coefficient was 89.7%. In this study, we have developed a new prostate segmentation approach based on the optimal feature learning framework, demonstrated its clinical feasibility, and validated its accuracy with manual segmentations.
Structural optimization of 3D-printed synthetic spider webs for high strength
Qin, Zhao; Compton, Brett G.; Lewis, Jennifer A.; Buehler, Markus J.
2015-01-01
Spiders spin intricate webs that serve as sophisticated prey-trapping architectures that simultaneously exhibit high strength, elasticity and graceful failure. To determine how web mechanics are controlled by their topological design and material distribution, here we create spider-web mimics composed of elastomeric filaments. Specifically, computational modelling and microscale 3D printing are combined to investigate the mechanical response of elastomeric webs under multiple loading conditions. We find the existence of an asymptotic prey size that leads to a saturated web strength. We identify pathways to design elastomeric material structures with maximum strength, low density and adaptability. We show that the loading type dictates the optimal material distribution, that is, a homogeneous distribution is better for localized loading, while stronger radial threads with weaker spiral threads is better for distributed loading. Our observations reveal that the material distribution within spider webs is dictated by the loading condition, shedding light on their observed architectural variations. PMID:25975372
Bayesian Segmentation of Atrium Wall Using Globally-Optimal Graph Cuts on 3D Meshes
Veni, Gopalkrishna; Fu, Zhisong; Awate, Suyash P.; Whitaker, Ross T.
2014-01-01
Efficient segmentation of the left atrium (LA) wall from delayed enhancement MRI is challenging due to inconsistent contrast, combined with noise, and high variation in atrial shape and size. We present a surface-detection method that is capable of extracting the atrial wall by computing an optimal a-posteriori estimate. This estimation is done on a set of nested meshes, constructed from an ensemble of segmented training images, and graph cuts on an associated multi-column, proper-ordered graph. The graph/mesh is a part of a template/model that has an associated set of learned intensity features. When this mesh is overlaid onto a test image, it produces a set of costs which lead to an optimal segmentation. The 3D mesh has an associated weighted, directed multi-column graph with edges that encode smoothness and inter-surface penalties. Unlike previous graph-cut methods that impose hard constraints on the surface properties, the proposed method follows from a Bayesian formulation resulting in soft penalties on spatial variation of the cuts through the mesh. The novelty of this method also lies in the construction of proper-ordered graphs on complex shapes for choosing among distinct classes of base shapes for automatic LA segmentation. We evaluate the proposed segmentation framework on simulated and clinical cardiac MRI. PMID:24684007
GPU-Accelerated Denoising in 3D (GD3D)
Energy Science and Technology Software Center (ESTSC)
2013-10-01
The raw computational power GPU Accelerators enables fast denoising of 3D MR images using bilateral filtering, anisotropic diffusion, and non-local means. This software addresses two facets of this promising application: what tuning is necessary to achieve optimal performance on a modern GPU? And what parameters yield the best denoising results in practice? To answer the first question, the software performs an autotuning step to empirically determine optimal memory blocking on the GPU. To answer themore » second, it performs a sweep of algorithm parameters to determine the combination that best reduces the mean squared error relative to a noiseless reference image.« less
An Update on Design Tools for Optimization of CMC 3D Fiber Architectures
NASA Technical Reports Server (NTRS)
Lang, J.; DiCarlo, J.
2012-01-01
Objective: Describe and up-date progress for NASA's efforts to develop 3D architectural design tools for CMC in general and for SIC/SiC composites in particular. Describe past and current sequential work efforts aimed at: Understanding key fiber and tow physical characteristics in conventional 2D and 3D woven architectures as revealed by microstructures in the literature. Developing an Excel program for down-selecting and predicting key geometric properties and resulting key fiber-controlled properties for various conventional 3D architectures. Developing a software tool for accurately visualizing all the key geometric details of conventional 3D architectures. Validating tools by visualizing and predicting the Internal geometry and key mechanical properties of a NASA SIC/SIC panel with a 3D orthogonal architecture. Applying the predictive and visualization tools toward advanced 3D orthogonal SiC/SIC composites, and combining them into a user-friendly software program.
NASA Astrophysics Data System (ADS)
Ji, Songbai; Fan, Xiaoyao; Roberts, David W.; Hartov, Alex; Paulsen, Keith D.
2011-03-01
Compensating for brain shift as surgery progresses is important to ensure sufficient accuracy in patient-to-image registration in the operating room (OR) for reliable neuronavigation. Ultrasound has emerged as an important and practical imaging technique for brain shift compensation either by itself or through computational modeling that estimates whole-brain deformation. Using volumetric true 3D ultrasound (3DUS), it is possible to nonrigidly (e.g., based on B-splines) register two temporally different 3DUS images directly to generate feature displacement maps for data assimilation in the biomechanical model. Because of a large amount of data and number of degrees-of-freedom (DOFs) involved, however, a significant computational cost may be required that can adversely influence the clinical feasibility of the technique for efficiently generating model-updated MR (uMR) in the OR. This paper parametrically investigates three B-splines registration parameters and their influence on the computational cost and registration accuracy: number of grid nodes along each direction, floating image volume down-sampling rate, and number of iterations. A simulated rigid body displacement field was employed as a ground-truth against which the accuracy of displacements generated from the B-splines nonrigid registration was compared. A set of optimal parameters was then determined empirically that result in a registration computational cost of less than 1 min and a sub-millimetric accuracy in displacement measurement. These resulting parameters were further applied to a clinical surgery case to demonstrate their practical use. Our results indicate that the optimal set of parameters result in sufficient accuracy and computational efficiency in model computation, which is important for future application of the overall biomechanical modeling to generate uMR for image-guidance in the OR.
Multifractal anisotropic swimming: the optimal foraging behaviour of grouper larvae.
Mahjoub, M S; Dur, G; Souissi, S; Schmitt, F G; Hwang, J S
2016-05-01
It was hypothesized that the Malabar grouper Ephinephelus malabaricus larvae have developed search patterns adapted to the distribution of their prey to maximise their net energy intake per unit time. Analysis of the swimming behaviour of E. malabaricus larvae in both the presence and absence of Artemia sp. nauplii is presented to test this hypothesis. A method derived from turbulence studies (the moment function of the displacements) was used to characterize the behaviour. The results revealed that larval swimming pattern was multifractal (intermittent and long-range-correlated) and isotropic (i.e. uniform in all directions) in the presence of prey, but multifractal and anisotropic (i.e. more frequent long displacement on the vertical axis) in the absence of prey. It is suggested that the search behaviour observed in the absence of prey is an adaptive response to prey distribution pattern, which is often characterised by multifractality and anisotropy (i.e. larger patches on the horizontal axes). In the presence of prey, E. malabaricus shifted to intensive search behaviour. Other possible contributors to the observed patterns are discussed. It is concluded that multifractality and anisotropy of swimming patterns observed in the experiment are mainly explained in an optimal foraging theory framework. PMID:27021375
NASA Astrophysics Data System (ADS)
Marcus, Guy; Parsa, Shima; Kramel, Stefan; Ni, Rui; Voth, Greg
2013-11-01
We have developed a general methodology to experimentally measure the time-resolved Lagrangian orientation and solid body rotation rate of anisotropic particles with arbitrary aspect ratio from standard stereoscopic video image data. We apply these techniques to particles advected in a Rλ ~ 110 fluid flow, where turbulence is generated by two grids oscillating in phase. We use 3D printing technology to design and fabricate neutrally buoyant rods, crosses (two perpendicular rods), and jacks (three mutually perpendicular rods) with a largest dimension of 7 times the Kolmogorov length scale, which makes them good approximations to tracer particles. We have measured the mean square rotation rate, ṗiṗi , of particles spanning the full range of aspect ratios and obtained results that agree with direct numerical simulations. By measuring the full solid-body rotation of jacks, we provide a new, extensible way to directly probe the Lagrangian vorticity of a fluid. We also present direct measurements of the alignment of crosses with the direction of their solid body rotation rate vector--in agreement with direct numerical simulations. Supported by NSF grant DMR1208990.
Wang, Xi-fen; Zhou, Huai-chun
2005-01-01
The control of 3-D temperature distribution in a utility boiler furnace is essential for the safe, economic and clean operation of pc-fired furnace with multi-burner system. The development of the visualization of 3-D temperature distributions in pc-fired furnaces makes it possible for a new combustion control strategy directly with the furnace temperature as its goal to improve the control quality for the combustion processes. Studied in this paper is such a new strategy that the whole furnace is divided into several parts in the vertical direction, and the average temperature and its bias from the center in every cross section can be extracted from the visualization results of the 3-D temperature distributions. In the simulation stage, a computational fluid dynamics (CFD) code served to calculate the 3-D temperature distributions in a furnace, then a linear model was set up to relate the features of the temperature distributions with the input of the combustion processes, such as the flow rates of fuel and air fed into the furnaces through all the burners. The adaptive genetic algorithm was adopted to find the optimal combination of the whole input parameters which ensure to form an optimal 3-D temperature field in the furnace desired for the operation of boiler. Simulation results showed that the strategy could soon find the factors making the temperature distribution apart from the optimal state and give correct adjusting suggestions. PMID:16295911
Zherdetsky, Aleksej; Prakonina, Alena; Malony, Allen D.
2014-01-01
The Electrical Impedance Tomography (EIT) and electroencephalography (EEG) forward problems in anisotropic inhomogeneous media like the human head belongs to the class of the three-dimensional boundary value problems for elliptic equations with mixed derivatives. We introduce and explore the performance of several new promising numerical techniques, which seem to be more suitable for solving these problems. The proposed numerical schemes combine the fictitious domain approach together with the finite-difference method and the optimally preconditioned Conjugate Gradient- (CG-) type iterative method for treatment of the discrete model. The numerical scheme includes the standard operations of summation and multiplication of sparse matrices and vector, as well as FFT, making it easy to implement and eligible for the effective parallel implementation. Some typical use cases for the EIT/EEG problems are considered demonstrating high efficiency of the proposed numerical technique. PMID:24527060
Turovets, Sergei; Volkov, Vasily; Zherdetsky, Aleksej; Prakonina, Alena; Malony, Allen D
2014-01-01
The Electrical Impedance Tomography (EIT) and electroencephalography (EEG) forward problems in anisotropic inhomogeneous media like the human head belongs to the class of the three-dimensional boundary value problems for elliptic equations with mixed derivatives. We introduce and explore the performance of several new promising numerical techniques, which seem to be more suitable for solving these problems. The proposed numerical schemes combine the fictitious domain approach together with the finite-difference method and the optimally preconditioned Conjugate Gradient- (CG-) type iterative method for treatment of the discrete model. The numerical scheme includes the standard operations of summation and multiplication of sparse matrices and vector, as well as FFT, making it easy to implement and eligible for the effective parallel implementation. Some typical use cases for the EIT/EEG problems are considered demonstrating high efficiency of the proposed numerical technique. PMID:24527060
NASA Astrophysics Data System (ADS)
Guo, Y.; Ding, M. D.; Wiegelmann, T.; Li, H.
2008-06-01
The photospheric vector magnetic field of the active region NOAA 10930 was obtained with the Solar Optical Telescope (SOT) on board the Hinode satellite with a very high spatial resolution (about 0.3''). Observations of the two-ribbon flare on 2006 December 13 in this active region provide us a good sample to study the magnetic field configuration related to the occurrence of the flare. Using the optimization method for nonlinear force-free field (NLFFF) extrapolation proposed by Wheatland et al. and recently developed by Wiegelmann, we derive the three-dimensional (3D) vector magnetic field configuration associated with this flare. The general topology can be described as a highly sheared core field and a quasi-potential envelope arch field. The core field clearly shows some dips supposed to sustain a filament. Free energy release in the flare, calculated by subtracting the energy contained in the NLFFF and the corresponding potential field, is 2.4 × 1031 ergs, which is ~2% of the preflare potential field energy. We also calculate the shear angles, defined as the angles between the NLFFF and potential field, and find that they become larger at some particular sites in the lower atmosphere, while they become significantly smaller in most places, implying that the whole configuration gets closer to the potential field after the flare. The Ca II H line images obtained with the Broadband Filter Imager (BFI) of the SOT and the 1600 Å images with the Transition Region and Coronal Explorer (TRACE) show that the preflare heating occurs mainly in the core field. These results provide evidence in support of the tether-cutting model of solar flares.
Flynn, A A; Pedley, R B; Green, A J; Boxer, G M; Boden, R; Begent, R H
2001-10-01
The biological effect of radioimmunotherapy (RIT) is most commonly assessed in terms of the absorbed radiation dose. In tumor, conventional dosimetry methods assume a uniform radionuclide and calculate a mean dose throughout the tumor. However, the vasculature of solid tumors tends to be highly irregular and the systemic delivery of antibodies is therefore heterogeneous. Tumor-specific antibodies preferentially localize in the viable, radiosensitive parts of the tumor whereas non-specific antibodies can penetrate into the necrosis where the dose is wasted. As a result, the observed biological effect can be very different to the predicted effect from conventional dose estimates. The purpose of this study is to assess the potential for optimizing the biological effect of RIT by matching the dose-distribution with tumor structure through the selection of appropriate antibodies and radionuclides. Storage phosphor plate technology was used to acquire images of the antibody distribution in serial tumor sections. Images of the distributions of a trivalent (TFM), bivalent (A5B7-IgG), monovalent (MFE-23) and a non-specific antibody (MOPC) were obtained. These images were registered with corresponding images showing tumor morphology. Serial images were reconstructed to form 3D maps of the antibody distribution and tumor structure. Convolution of the image of antibody distribution with beta dose point kernals generated dose-rate distributions for 14C, 131I and 90Y. These were statistically compared with the tumor structure. The highest correlation was obtained for the multivalent antibodies combined with 131I, due to specific retention in viable areas of tumor coupled with the fact that much of the dose was deposted locally. With decreasing avidity the correlation also decreased and with the non-specific antibody this correlation was negative, indicating higher concentrations in the necrotic regions. In conclusion, the dose distribution can be optimized in tumor by selecting
NASA Astrophysics Data System (ADS)
Lin, F. C.; Schmandt, B.; Tsai, V. C.
2014-12-01
The EarthScope USArray Transportable Array (TA) has provided a great opportunity for imaging the detailed lithospheric structure beneath the continental US. In this presentation, we will report our recent progress on constructing detailed 3D isotropic and anisotropic crustal models of the contiguous US using Rayleigh wave phase velocity and ellipticity measurements across TA. In particular, we will discuss our recent methodology development of extracting short period Rayleigh wave ellipticity, or Rayleigh-wave H/V (horizontal to vertical) amplitude ratios, using multicomponent noise cross-correlations. To retain the amplitude ratio information between vertical and horizontal components, for each station, we perform daily noise pre-processing (temporal normalization and spectrum whitening) simultaneously for all three components. For each station pair, amplitude measurements between cross-correlations of different components (radial-radial, radial-vertical, vertical-radial and vertical-vertical) are then used to determine the Rayleigh-wave H/V ratios at the two station locations. Measurements from all available station pairs are used to determine isotropic and directionally dependent Rayleigh-wave H/V ratios at each location between 8- and 24-second period. The isotropic H/V ratio maps, combined with previous longer period Rayleigh-wave H/V ratio maps from earthquakes and Rayleigh-wave phase velocity maps from both ambient noise and earthquakes, are used to invert for a new 3-D isotropic crustal and upper-mantle model in the western United States. The new model has an outstanding vertical resolution in the upper crust and tradeoffs between different parameters are mitigated. A clear 180-degree periodicity is observed in the directionally dependent H/V ratio measurements for many locations where upper crustal anisotropy is likely strong. Across the US, good correlation is observed between the inferred fast directions in the upper crust and documented maximum
MSE optimal bit-rate allocation in JPEG2000 Part 2 compression applied to a 3D data set
NASA Astrophysics Data System (ADS)
Kosheleva, Olga M.; Cabrera, Sergio D.; Usevitch, Bryan E.; Aguirre, Alberto; Vidal, Edward, Jr.
2004-10-01
A bit rate allocation (BRA) strategy is needed to optimally compress three-dimensional (3-D) data on a per-slice basis, treating it as a collection of two-dimensional (2-D) slices/components. This approach is compatible with the framework of JPEG2000 Part 2 which includes the option of pre-processing the slices with a decorrelation transform in the cross-component direction so that slices of transform coefficients are compressed. In this paper, we illustrate the impact of a recently developed inter-slice rate-distortion optimal bit-rate allocation approach that is applicable to this compression system. The approach exploits the MSE optimality of all JPEG2000 bit streams for all slices when each is produced in the quality progressive mode. Each bit stream can be used to produce a rate-distortion curve (RDC) for each slice that is MSE optimal at each bit rate of interest. The inter-slice allocation approach uses all RDCs for all slices to optimally select an overall optimal set of bit rates for all the slices using a constrained optimization procedure. The optimization is conceptually similar to Post-Compression Rate-Distortion optimization that is used within JPEG2000 to optimize bit rates allocated to codeblocks. Results are presented for two types of data sets: a 3-D computed tomography (CT) medical image, and a 3-D metereological data set derived from a particular modeling program. For comparison purposes, compression results are also illustrated for the traditional log-variance approach and for a uniform allocation strategy. The approach is illustrated using two decorrelation tranforms (the Karhunen Loeve transform, and the discrete wavelet transform) for which the inter-slice allocation scheme has the most impact.
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
Development and Optimization of Viable Human Platforms through 3D Printing
Parker, Paul R.; Moya, Monica L.; Wheeler, Elizabeth K.
2015-08-21
3D printing technology offers a unique method for creating cell cultures in a manner far more conducive to accurate representation of human tissues and systems. Here we print cellular structures capable of forming vascular networks and exhibiting qualities of natural tissues and human systems. This allows for cheaper and readily available sources for further study of biological and pharmaceutical agents.
PRESAGETM - Development and optimization studies of a 3D radiochromic plastic dosimeter - Part 1
NASA Astrophysics Data System (ADS)
Adamovics, J.; Jordan, K.; Dietrich, J.
2006-12-01
This paper studies the polymerization of six different transparent plastics as potential 3D dosimeter matrices. In addition, six different leuco dyes and sixteen different free radical initiators were evaluated. Finally, the photoreactivity of the dosimeter was studied so that the effect of exposure to UV could be minimized.
Optimizing RetroICor and RetroKCor corrections for multi-shot 3D FMRI acquisitions
Tijssen, Rob H.N.; Jenkinson, Mark; Brooks, Jonathan C.W.; Jezzard, Peter; Miller, Karla L.
2015-01-01
Physiological noise, if unaccounted for, can drastically reduce the statistical significance of detected activation in FMRI. In this paper, we systematically optimize physiological noise regressions for multi-shot 3D FMRI data. First, we investigate whether 3D FMRI data are best corrected in image space (RetroICor) or k-space (RetroKCor), in which each k-space segment can be assigned its unique physiological phase. In addition, the optimal regressor set is determined using the Bayesian Information Criterion (BIC) for a variety of 3D acquisitions corresponding to different image contrasts and k-space readouts. Our simulations and experiments indicate that: (a) k-space corrections are more robust when performed on real/imaginary than magnitude/phase data; (b) k-space corrections do not outperform image-space corrections, despite the ability to synchronize physiological phase to acquisition time more accurately; and (c) the optimal model varied considerably between the various acquisition techniques. These results suggest the use of a tailored set of volume-wide regressors, determined by BIC or other selection criteria, that achieves optimal balance between variance reduction and potential over-fitting. PMID:24018307
NASA Astrophysics Data System (ADS)
van der Bom, I. M. J.; Klein, S.; Staring, M.; Homan, R.; Bartels, L. W.; Pluim, J. P. W.
2011-03-01
The advantage of 2D-3D image registration methods versus direct image-to-patient registration, is that these methods generally do not require user interaction (such as manual annotations), additional machinery or additional acquisition of 3D data. A variety of intensity-based similarity measures has been proposed and evaluated for different applications. These studies showed that the registration accuracy and capture range are influenced by the choice of similarity measure. However, the influence of the optimization method on intensity-based 2D-3D image registration has not been investigated. We have compared the registration performance of seven optimization methods in combination with three similarity measures: gradient difference, gradient correlation, and pattern intensity. Optimization methods included in this study were: regular step gradient descent, Nelder-Mead, Powell-Brent, Quasi-Newton, nonlinear conjugate gradient, simultaneous perturbation stochastic approximation, and evolution strategy. Registration experiments were performed on multiple patient data sets that were obtained during cerebral interventions. Various component combinations were evaluated on registration accuracy, capture range, and registration time. The results showed that for the same similarity measure, different registration accuracies and capture ranges were obtained when different optimization methods were used. For gradient difference, largest capture ranges were obtained with Powell-Brent and simultaneous perturbation stochastic approximation. Gradient correlation and pattern intensity had the largest capture ranges in combination with Powell-Brent, Nelder-Mead, nonlinear conjugate gradient, and Quasi-Newton. Average registration time, expressed in the number of DRRs required for convergence, was the lowest for Powell-Brent. Based on these results, we conclude that Powell-Brent is a reliable optimization method for intensity-based 2D-3D registration of x-ray images to CBCT
NASA Astrophysics Data System (ADS)
Pekşen, Ertan; Yas, Türker; Kıyak, Alper
2014-09-01
We examine the one-dimensional direct current method in anisotropic earth formation. We derive an analytic expression of a simple, two-layered anisotropic earth model. Further, we also consider a horizontally layered anisotropic earth response with respect to the digital filter method, which yields a quasi-analytic solution over anisotropic media. These analytic and quasi-analytic solutions are useful tests for numerical codes. A two-dimensional finite difference earth model in anisotropic media is presented in order to generate a synthetic data set for a simple one-dimensional earth. Further, we propose a particle swarm optimization method for estimating the model parameters of a layered anisotropic earth model such as horizontal and vertical resistivities, and thickness. The particle swarm optimization is a naturally inspired meta-heuristic algorithm. The proposed method finds model parameters quite successfully based on synthetic and field data. However, adding 5 % Gaussian noise to the synthetic data increases the ambiguity of the value of the model parameters. For this reason, the results should be controlled by a number of statistical tests. In this study, we use probability density function within 95 % confidence interval, parameter variation of each iteration and frequency distribution of the model parameters to reduce the ambiguity. The result is promising and the proposed method can be used for evaluating one-dimensional direct current data in anisotropic media.
Shape optimization of 3D continuum structures via force approximation techniques
NASA Technical Reports Server (NTRS)
Vanderplaats, Garret N.; Kodiyalam, Srinivas
1988-01-01
The existing need to develop methods whereby the shape design efficiency can be improved through the use of high quality approximation methods is addressed. An efficient approximation method for stress constraints in 3D shape design problems is proposed based on expanding the nodal forces in Taylor series with respect to shape variations. The significance of this new method is shown through elementary beam theory calculations and via numerical computations using 3D solid finite elements. Numerical examples including the classical cantilever beam structure and realistic automotive parts like the engine connecting rod are designed for optimum shape using the proposed method. The numerical results obtained from these methods are compared with other published results, to assess the efficiency and the convergence rate of the proposed method.
Optimization of spine surgery planning with 3D image templating tools
NASA Astrophysics Data System (ADS)
Augustine, Kurt E.; Huddleston, Paul M.; Holmes, David R., III; Shridharani, Shyam M.; Robb, Richard A.
2008-03-01
The current standard of care for patients with spinal disorders involves a thorough clinical history, physical exam, and imaging studies. Simple radiographs provide a valuable assessment but prove inadequate for surgery planning because of the complex 3-dimensional anatomy of the spinal column and the close proximity of the neural elements, large blood vessels, and viscera. Currently, clinicians still use primitive techniques such as paper cutouts, pencils, and markers in an attempt to analyze and plan surgical procedures. 3D imaging studies are routinely ordered prior to spine surgeries but are currently limited to generating simple, linear and angular measurements from 2D views orthogonal to the central axis of the patient. Complex spinal corrections require more accurate and precise calculation of 3D parameters such as oblique lengths, angles, levers, and pivot points within individual vertebra. We have developed a clinician friendly spine surgery planning tool which incorporates rapid oblique reformatting of each individual vertebra, followed by interactive templating for 3D placement of implants. The template placement is guided by the simultaneous representation of multiple 2D section views from reformatted orthogonal views and a 3D rendering of individual or multiple vertebrae enabling superimposition of virtual implants. These tools run efficiently on desktop PCs typically found in clinician offices or workrooms. A preliminary study conducted with Mayo Clinic spine surgeons using several actual cases suggests significantly improved accuracy of pre-operative measurements and implant localization, which is expected to increase spinal procedure efficiency and safety, and reduce time and cost of the operation.
Segmentation of bone structures in 3D CT images based on continuous max-flow optimization
NASA Astrophysics Data System (ADS)
Pérez-Carrasco, J. A.; Acha-Piñero, B.; Serrano, C.
2015-03-01
In this paper an algorithm to carry out the automatic segmentation of bone structures in 3D CT images has been implemented. Automatic segmentation of bone structures is of special interest for radiologists and surgeons to analyze bone diseases or to plan some surgical interventions. This task is very complicated as bones usually present intensities overlapping with those of surrounding tissues. This overlapping is mainly due to the composition of bones and to the presence of some diseases such as Osteoarthritis, Osteoporosis, etc. Moreover, segmentation of bone structures is a very time-consuming task due to the 3D essence of the bones. Usually, this segmentation is implemented manually or with algorithms using simple techniques such as thresholding and thus providing bad results. In this paper gray information and 3D statistical information have been combined to be used as input to a continuous max-flow algorithm. Twenty CT images have been tested and different coefficients have been computed to assess the performance of our implementation. Dice and Sensitivity values above 0.91 and 0.97 respectively were obtained. A comparison with Level Sets and thresholding techniques has been carried out and our results outperformed them in terms of accuracy.
Giżyńska, Marta K.; Kukołowicz, Paweł F.; Kordowski, Paweł
2014-01-01
Aim The aim of this work is to present a method of beam weight and wedge angle optimization for patients with prostate cancer. Background 3D-CRT is usually realized with forward planning based on a trial and error method. Several authors have published a few methods of beam weight optimization applicable to the 3D-CRT. Still, none on these methods is in common use. Materials and methods Optimization is based on the assumption that the best plan is achieved if dose gradient at ICRU point is equal to zero. Our optimization algorithm requires beam quality index, depth of maximum dose, profiles of wedged fields and maximum dose to femoral heads. The method was tested for 10 patients with prostate cancer, treated with the 3-field technique. Optimized plans were compared with plans prepared by 12 experienced planners. Dose standard deviation in target volume, and minimum and maximum doses were analyzed. Results The quality of plans obtained with the proposed optimization algorithms was comparable to that prepared by experienced planners. Mean difference in target dose standard deviation was 0.1% in favor of the plans prepared by planners for optimization of beam weights and wedge angles. Introducing a correction factor for patient body outline for dose gradient at ICRU point improved dose distribution homogeneity. On average, a 0.1% lower standard deviation was achieved with the optimization algorithm. No significant difference in mean dose–volume histogram for the rectum was observed. Conclusions Optimization shortens very much time planning. The average planning time was 5 min and less than a minute for forward and computer optimization, respectively. PMID:25337411
Li, Bai; Lin, Mu; Liu, Qiao; Li, Ya; Zhou, Changjun
2015-10-01
Protein folding is a fundamental topic in molecular biology. Conventional experimental techniques for protein structure identification or protein folding recognition require strict laboratory requirements and heavy operating burdens, which have largely limited their applications. Alternatively, computer-aided techniques have been developed to optimize protein structures or to predict the protein folding process. In this paper, we utilize a 3D off-lattice model to describe the original protein folding scheme as a simplified energy-optimal numerical problem, where all types of amino acid residues are binarized into hydrophobic and hydrophilic ones. We apply a balance-evolution artificial bee colony (BE-ABC) algorithm as the minimization solver, which is featured by the adaptive adjustment of search intensity to cater for the varying needs during the entire optimization process. In this work, we establish a benchmark case set with 13 real protein sequences from the Protein Data Bank database and evaluate the convergence performance of BE-ABC algorithm through strict comparisons with several state-of-the-art ABC variants in short-term numerical experiments. Besides that, our obtained best-so-far protein structures are compared to the ones in comprehensive previous literature. This study also provides preliminary insights into how artificial intelligence techniques can be applied to reveal the dynamics of protein folding. Graphical Abstract Protein folding optimization using 3D off-lattice model and advanced optimization techniques. PMID:26381910
NASA Astrophysics Data System (ADS)
Cipriani, L.; Fantini, F.; Bertacchi, S.
2014-06-01
Image-based modelling tools based on SfM algorithms gained great popularity since several software houses provided applications able to achieve 3D textured models easily and automatically. The aim of this paper is to point out the importance of controlling models parameterization process, considering that automatic solutions included in these modelling tools can produce poor results in terms of texture utilization. In order to achieve a better quality of textured models from image-based modelling applications, this research presents a series of practical strategies aimed at providing a better balance between geometric resolution of models from passive sensors and their corresponding (u,v) map reference systems. This aspect is essential for the achievement of a high-quality 3D representation, since "apparent colour" is a fundamental aspect in the field of Cultural Heritage documentation. Complex meshes without native parameterization have to be "flatten" or "unwrapped" in the (u,v) parameter space, with the main objective to be mapped with a single image. This result can be obtained by using two different strategies: the former automatic and faster, while the latter manual and time-consuming. Reverse modelling applications provide automatic solutions based on splitting the models by means of different algorithms, that produce a sort of "atlas" of the original model in the parameter space, in many instances not adequate and negatively affecting the overall quality of representation. Using in synergy different solutions, ranging from semantic aware modelling techniques to quad-dominant meshes achieved using retopology tools, it is possible to obtain a complete control of the parameterization process.
Density and level set-XFEM schemes for topology optimization of 3-D structures
NASA Astrophysics Data System (ADS)
Villanueva, Carlos H.; Maute, Kurt
2014-07-01
As the capabilities of additive manufacturing techniques increase, topology optimization provides a promising approach to design geometrically sophisticated structures. Traditional topology optimization methods aim at finding conceptual designs, but they often do not resolve sufficiently the geometry and the structural response such that the optimized designs can be directly used for manufacturing. To overcome these limitations, this paper studies the viability of the extended finite element method (XFEM) in combination with the level-set method (LSM) for topology optimization of three dimensional structures. The LSM describes the geometry by defining the nodal level set values via explicit functions of the optimization variables. The structural response is predicted by a generalized version of the XFEM. The LSM-XFEM approach is compared against results from a traditional Solid Isotropic Material with Penalization method for two-phase "solid-void" and "solid-solid" problems. The numerical results demonstrate that the LSM-XFEM approach describes crisply the geometry and predicts the structural response with acceptable accuracy even on coarse meshes.
Structure Design of the 3-D Braided Composite Based on a Hybrid Optimization Algorithm
NASA Astrophysics Data System (ADS)
Zhang, Ke
Three-dimensional braided composite has the better designable characteristic. Whereas wide application of hollow-rectangular-section three-dimensional braided composite in engineering, optimization design of the three-dimensional braided composite made by 4-step method were introduced. Firstly, the stiffness and damping characteristic analysis of the composite is presented. Then, the mathematical models for structure design of the three-dimensional braided composite were established. The objective functions are based on the specific damping capacity and stiffness of the composite. The design variables are the braiding parameters of the composites and sectional geometrical size of the composite. The optimization problem is solved by using ant colony optimization (ACO), contenting the determinate restriction. The results of numeral examples show that the better damping and stiffness characteristic could be obtained. The method proposed here is useful for the structure design of the kind of member and its engineering application.
A Case Study of Communication Optimizations on 3D Mesh Interconnects
NASA Astrophysics Data System (ADS)
Bhatelé, Abhinav; Bohm, Eric; Kalé, Laxmikant V.
Optimal network performance is critical to efficient parallel scaling for communication-bound applications on large machines. With wormhole routing, no-load latencies do not increase significantly with number of hops traveled. Yet, we, and others have recently shown that in presence of contention, message latencies can grow substantially large. Hence task mapping strategies should take the topology of the machine into account on large machines. In this paper, we present topology aware mapping as a technique to optimize communication on 3-dimensional mesh interconnects and hence improve performance.
Optimization of 3D laser scanning speed by use of combined variable step
NASA Astrophysics Data System (ADS)
Garcia-Cruz, X. M.; Sergiyenko, O. Yu.; Tyrsa, Vera; Rivas-Lopez, M.; Hernandez-Balbuena, D.; Rodriguez-Quiñonez, J. C.; Basaca-Preciado, L. C.; Mercorelli, P.
2014-03-01
The problem of 3D TVS slow functioning caused by constant small scanning step becomes its solution in the presented research. It can be achieved by combined scanning step application for the fast search of n obstacles in unknown surroundings. Such a problem is of keynote importance in automatic robot navigation. To maintain a reasonable speed robots must detect dangerous obstacles as soon as possible, but all known scanners able to measure distances with sufficient accuracy are unable to do it in real time. So, the related technical task of the scanning with variable speed and precise digital mapping only for selected spatial sectors is under consideration. A wide range of simulations in MATLAB 7.12.0 of several variants of hypothetic scenes with variable n obstacles in each scene (including variation of shapes and sizes) and scanning with incremented angle value (0.6° up to 15°) is provided. The aim of such simulation was to detect which angular values of interval still permit getting the maximal information about obstacles without undesired time losses. Three of such local maximums were obtained in simulations and then rectified by application of neuronal network formalism (Levenberg-Marquradt Algorithm). The obtained results in its turn were applied to MET (Micro-Electro-mechanical Transmission) design for practical realization of variable combined step scanning on an experimental prototype of our previously known laser scanner.
NASA Astrophysics Data System (ADS)
Bhowmik, Tanmoy; Liu, Hanli; Ye, Zhou; Oraintara, Soontorn
2016-03-01
Diffuse optical tomography (DOT) is a relatively low cost and portable imaging modality for reconstruction of optical properties in a highly scattering medium, such as human tissue. The inverse problem in DOT is highly ill-posed, making reconstruction of high-quality image a critical challenge. Because of the nature of sparsity in DOT, sparsity regularization has been utilized to achieve high-quality DOT reconstruction. However, conventional approaches using sparse optimization are computationally expensive and have no selection criteria to optimize the regularization parameter. In this paper, a novel algorithm, Dimensionality Reduction based Optimization for DOT (DRO-DOT), is proposed. It reduces the dimensionality of the inverse DOT problem by reducing the number of unknowns in two steps and thereby makes the overall process fast. First, it constructs a low resolution voxel basis based on the sensing-matrix properties to find an image support. Second, it reconstructs the sparse image inside this support. To compensate for the reduced sensitivity with increasing depth, depth compensation is incorporated in DRO-DOT. An efficient method to optimally select the regularization parameter is proposed for obtaining a high-quality DOT image. DRO-DOT is also able to reconstruct high-resolution images even with a limited number of optodes in a spatially limited imaging set-up.
Bhowmik, Tanmoy; Liu, Hanli; Ye, Zhou; Oraintara, Soontorn
2016-01-01
Diffuse optical tomography (DOT) is a relatively low cost and portable imaging modality for reconstruction of optical properties in a highly scattering medium, such as human tissue. The inverse problem in DOT is highly ill-posed, making reconstruction of high-quality image a critical challenge. Because of the nature of sparsity in DOT, sparsity regularization has been utilized to achieve high-quality DOT reconstruction. However, conventional approaches using sparse optimization are computationally expensive and have no selection criteria to optimize the regularization parameter. In this paper, a novel algorithm, Dimensionality Reduction based Optimization for DOT (DRO-DOT), is proposed. It reduces the dimensionality of the inverse DOT problem by reducing the number of unknowns in two steps and thereby makes the overall process fast. First, it constructs a low resolution voxel basis based on the sensing-matrix properties to find an image support. Second, it reconstructs the sparse image inside this support. To compensate for the reduced sensitivity with increasing depth, depth compensation is incorporated in DRO-DOT. An efficient method to optimally select the regularization parameter is proposed for obtaining a high-quality DOT image. DRO-DOT is also able to reconstruct high-resolution images even with a limited number of optodes in a spatially limited imaging set-up. PMID:26940661
Bhowmik, Tanmoy; Liu, Hanli; Ye, Zhou; Oraintara, Soontorn
2016-01-01
Diffuse optical tomography (DOT) is a relatively low cost and portable imaging modality for reconstruction of optical properties in a highly scattering medium, such as human tissue. The inverse problem in DOT is highly ill-posed, making reconstruction of high-quality image a critical challenge. Because of the nature of sparsity in DOT, sparsity regularization has been utilized to achieve high-quality DOT reconstruction. However, conventional approaches using sparse optimization are computationally expensive and have no selection criteria to optimize the regularization parameter. In this paper, a novel algorithm, Dimensionality Reduction based Optimization for DOT (DRO-DOT), is proposed. It reduces the dimensionality of the inverse DOT problem by reducing the number of unknowns in two steps and thereby makes the overall process fast. First, it constructs a low resolution voxel basis based on the sensing-matrix properties to find an image support. Second, it reconstructs the sparse image inside this support. To compensate for the reduced sensitivity with increasing depth, depth compensation is incorporated in DRO-DOT. An efficient method to optimally select the regularization parameter is proposed for obtaining a high-quality DOT image. DRO-DOT is also able to reconstruct high-resolution images even with a limited number of optodes in a spatially limited imaging set-up. PMID:26940661
NASA Astrophysics Data System (ADS)
Marquet, F.; Aubry, J. F.; Pernot, M.; Fink, M.; Tanter, M.
2011-11-01
Recent studies have demonstrated the feasibility of transcostal high intensity focused ultrasound (HIFU) treatment in liver. However, two factors limit thermal necrosis of the liver through the ribs: the energy deposition at focus is decreased by the respiratory movement of the liver and the energy deposition on the skin is increased by the presence of highly absorbing bone structures. Ex vivo ablations were conducted to validate the feasibility of a transcostal real-time 3D movement tracking and correction mode. Experiments were conducted through a chest phantom made of three human ribs immersed in water and were placed in front of a 300 element array working at 1 MHz. A binarized apodization law introduced recently in order to spare the rib cage during treatment has been extended here with real-time electronic steering of the beam. Thermal simulations have been conducted to determine the steering limits. In vivo 3D-movement detection was performed on pigs using an ultrasonic sequence. The maximum error on the transcostal motion detection was measured to be 0.09 ± 0.097 mm on the anterior-posterior axis. Finally, a complete sequence was developed combining real-time 3D transcostal movement correction and spiral trajectory of the HIFU beam, allowing the system to treat larger areas with optimized efficiency. Lesions as large as 1 cm in diameter have been produced at focus in excised liver, whereas no necroses could be obtained with the same emitted power without correcting the movement of the tissue sample.
NASA Astrophysics Data System (ADS)
Ivy, D. J.; Rigby, M. L.; Prinn, R. G.; Muhle, J.; Weiss, R. F.
2009-12-01
We present optimized annual global emissions from 1973-2008 of nitrogen trifluoride (NF3), a powerful greenhouse gas which is not currently regulated by the Kyoto Protocol. In the past few decades, NF3 production has dramatically increased due to its usage in the semiconductor industry. Emissions were estimated through the 'pulse-method' discrete Kalman filter using both a simple, flexible 2-D 12-box model used in the Advanced Global Atmospheric Gases Experiment (AGAGE) network and the Model for Ozone and Related Tracers (MOZART v4.5), a full 3-D atmospheric chemistry model. No official audited reports of industrial NF3 emissions are available, and with limited information on production, a priori emissions were estimated using both a bottom-up and top-down approach with two different spatial patterns based on semiconductor perfluorocarbon (PFC) emissions from the Emission Database for Global Atmospheric Research (EDGAR v3.2) and Semiconductor Industry Association sales information. Both spatial patterns used in the models gave consistent results, showing the robustness of the estimated global emissions. Differences between estimates using the 2-D and 3-D models can be attributed to transport rates and resolution differences. Additionally, new NF3 industry production and market information is presented. Emission estimates from both the 2-D and 3-D models suggest that either the assumed industry release rate of NF3 or industry production information is still underestimated.
APEnet+: a 3D Torus network optimized for GPU-based HPC Systems
NASA Astrophysics Data System (ADS)
Ammendola, R.; Biagioni, A.; Frezza, O.; Lo Cicero, F.; Lonardo, A.; Paolucci, P. S.; Rossetti, D.; Simula, F.; Tosoratto, L.; Vicini, P.
2012-12-01
In the supercomputing arena, the strong rise of GPU-accelerated clusters is a matter of fact. Within INFN, we proposed an initiative — the QUonG project — whose aim is to deploy a high performance computing system dedicated to scientific computations leveraging on commodity multi-core processors coupled with latest generation GPUs. The inter-node interconnection system is based on a point-to-point, high performance, low latency 3D torus network which is built in the framework of the APEnet+ project. It takes the form of an FPGA-based PCIe network card exposing six full bidirectional links running at 34 Gbps each that implements the RDMA protocol. In order to enable significant access latency reduction for inter-node data transfer, a direct network-to-GPU interface was built. The specialized hardware blocks, integrated in the APEnet+ board, provide support for GPU-initiated communications using the so called PCIe peer-to-peer (P2P) transactions. This development is made in close collaboration with the GPU vendor NVIDIA. The final shape of a complete QUonG deployment is an assembly of standard 42U racks, each one capable of 80 TFLOPS/rack of peak performance, at a cost of 5 k€/T F LOPS and for an estimated power consumption of 25 kW/rack. In this paper we report on the status of final rack deployment and on the R&D activities for 2012 that will focus on performance enhancement of the APEnet+ hardware through the adoption of new generation 28 nm FPGAs allowing the implementation of PCIe Gen3 host interface and the addition of new fault tolerance-oriented capabilities.
NASA Astrophysics Data System (ADS)
Li, Y.; Han, B.; Métivier, L.; Brossier, R.
2016-09-01
We investigate an optimal fourth-order staggered-grid finite-difference scheme for 3D frequency-domain viscoelastic wave modeling. An anti-lumped mass strategy is incorporated to minimize the numerical dispersion. The optimal finite-difference coefficients and the mass weighting coefficients are obtained by minimizing the misfit between the normalized phase velocities and the unity. An iterative damped least-squares method, the Levenberg-Marquardt algorithm, is utilized for the optimization. Dispersion analysis shows that the optimal fourth-order scheme presents less grid dispersion and anisotropy than the conventional fourth-order scheme with respect to different Poisson's ratios. Moreover, only 3.7 grid-points per minimum shear wavelength are required to keep the error of the group velocities below 1%. The memory cost is then greatly reduced due to a coarser sampling. A parallel iterative method named CARP-CG is used to solve the large ill-conditioned linear system for the frequency-domain modeling. Validations are conducted with respect to both the analytic viscoacoustic and viscoelastic solutions. Compared with the conventional fourth-order scheme, the optimal scheme generates wavefields having smaller error under the same discretization setups. Profiles of the wavefields are presented to confirm better agreement between the optimal results and the analytic solutions.
NASA Astrophysics Data System (ADS)
Ansola, R.; Veguería, E.; Alonso, C.; Querin, O. M.
2016-03-01
This work presents a sequential element rejection and admission (SERA) method for optimum topology design of three dimensional compliant actuators. The proposed procedure has been successfully applied to several topology optimization problems, but most investigations for compliant devices design have been focused on planar systems. This investigation aims to progress on this line, where a generalization of the method for three dimensional topology optimization is explored. The methodology described in this work is useful for the synthesis of high performance flexure based micro and nano manipulation applications demanding for both sensing and control of motion and force trajectories. In this case the goal of the topology optimization problem is to design an actuator that transfers work from the input point to the output port in a structurally efficient way. Here we will use the classical formulation where the displacement performed on a work piece modelled by a spring is maximized. The technique implemented works with two separate criteria for the rejection and admission of elements to efficiently achieve the optimum design and overcomes problems encountered by other evolutionary methods when dealing with compliant mechanisms design. The use of the algorithm is demonstrated through several numerical examples.
Tunesi, Marta; Fusco, Federica; Fiordaliso, Fabio; Corbelli, Alessandro; Biella, Gloria; Raimondi, Manuela T
2016-01-01
Frontotemporal lobar degeneration (FTLD) is a severe neurodegenerative disorder that is diagnosed with increasing frequency in clinical setting. Currently, no therapy is available and in addition the molecular basis of the disease are far from being elucidated. Consequently, it is of pivotal importance to develop reliable and cost-effective in vitro models for basic research purposes and drug screening. To this respect, recent results in the field of Alzheimer's disease have suggested that a tridimensional (3D) environment is an added value to better model key pathologic features of the disease. Here, we have tried to add complexity to the 3D cell culturing concept by using a microfluidic bioreactor, where cells are cultured under a continuous flow of medium, thus mimicking the interstitial fluid movement that actually perfuses the body tissues, including the brain. We have implemented this model using a neuronal-like cell line (SH-SY5Y), a widely exploited cell model for neurodegenerative disorders that shows some basic features relevant for FTLD modeling, such as the release of the FTLD-related protein progranulin (PRGN) in specific vesicles (exosomes). We have efficiently seeded the cells on 3D scaffolds, optimized a disease-relevant oxidative stress experiment (by targeting mitochondrial function that is one of the possible FTLD-involved pathological mechanisms) and evaluated cell metabolic activity in dynamic culture in comparison to static conditions, finding that SH-SY5Y cells cultured in 3D scaffold are susceptible to the oxidative damage triggered by a mitochondrial-targeting toxin (6-OHDA) and that the same cells cultured in dynamic conditions kept their basic capacity to secrete PRGN in exosomes once recovered from the bioreactor and plated in standard 2D conditions. We think that a further improvement of our microfluidic system may help in providing a full device where assessing basic FTLD-related features (including PRGN dynamic secretion) that may be
Optimal 3-D culture of primary articular chondrocytes for use in the Rotating Wall Vessel Bioreactor
Mellor, Liliana F.; Baker, Travis L.; Brown, Raquel J.; Catlin, Lindsey W.; Oxford, Julia Thom
2014-01-01
INTRODUCTION Reliable culturing methods for primary articular chondrocytes are essential to study the effects of loading and unloading on joint tissue at the cellular level. Due to the limited proliferation capacity of primary chondrocytes and their tendency to dedifferentiate in conventional culture conditions, long-term culturing conditions of primary chondrocytes can be challenging. The goal of this study was to develop a suspension culturing technique that not only would retain the cellular morphology but also maintain gene expression characteristics of primary articular chondrocytes. METHODS Three-dimensional culturing methods were compared and optimized for primary articular chondrocytes in the rotating wall vessel bioreactor, which changes the mechanical culture conditions to provide a form of suspension culture optimized for low shear and turbulence. We performed gene expression analysis and morphological characterization of cells cultured in alginate beads, Cytopore-2 microcarriers, primary monolayer culture, and passaged monolayer cultures using reverse transcription-PCR and laser scanning confocal microscopy. RESULTS Primary chondrocytes grown on Cytopore-2 microcarriers maintained the phenotypical morphology and gene expression pattern observed in primary bovine articular chondrocytes, and retained these characteristics for up to 9 days. DISCUSSION Our results provide a novel and alternative culturing technique for primary chondrocytes suitable for studies that require suspension such as those using the rotating wall vessel bioreactor. In addition, we provide an alternative culturing technique for primary chondrocytes that can impact future mechanistic studies of osteoarthritis progression, treatments for cartilage damage and repair, and cartilage tissue engineering. PMID:25199120
Jung, Jangwook P.; Hu, Dongjian; Domian, Ibrahim J.; Ogle, Brenda M.
2015-01-01
The extracellular matrix (ECM) impacts stem cell differentiation, but identifying formulations supportive of differentiation is challenging in 3D models. Prior efforts involving combinatorial ECM arrays seemed intuitively advantageous. We propose an alternative that suggests reducing sample size and technological burden can be beneficial and accessible when coupled to design of experiments approaches. We predict optimized ECM formulations could augment differentiation of cardiomyocytes derived in vitro. We employed native chemical ligation to polymerize 3D poly (ethylene glycol) hydrogels under mild conditions while entrapping various combinations of ECM and murine induced pluripotent stem cells. Systematic optimization for cardiomyocyte differentiation yielded a predicted solution of 61%, 24%, and 15% of collagen type I, laminin-111, and fibronectin, respectively. This solution was confirmed by increased numbers of cardiac troponin T, α-myosin heavy chain and α-sarcomeric actinin-expressing cells relative to suboptimum solutions. Cardiomyocytes of composites exhibited connexin43 expression, appropriate contractile kinetics and intracellular calcium handling. Further, adding a modulator of adhesion, thrombospondin-1, abrogated cardiomyocyte differentiation. Thus, the integrated biomaterial platform statistically identified an ECM formulation best supportive of cardiomyocyte differentiation. In future, this formulation could be coupled with biochemical stimulation to improve functional maturation of cardiomyocytes derived in vitro or transplanted in vivo. PMID:26687770
Jung, Jangwook P; Hu, Dongjian; Domian, Ibrahim J; Ogle, Brenda M
2015-01-01
The extracellular matrix (ECM) impacts stem cell differentiation, but identifying formulations supportive of differentiation is challenging in 3D models. Prior efforts involving combinatorial ECM arrays seemed intuitively advantageous. We propose an alternative that suggests reducing sample size and technological burden can be beneficial and accessible when coupled to design of experiments approaches. We predict optimized ECM formulations could augment differentiation of cardiomyocytes derived in vitro. We employed native chemical ligation to polymerize 3D poly (ethylene glycol) hydrogels under mild conditions while entrapping various combinations of ECM and murine induced pluripotent stem cells. Systematic optimization for cardiomyocyte differentiation yielded a predicted solution of 61%, 24%, and 15% of collagen type I, laminin-111, and fibronectin, respectively. This solution was confirmed by increased numbers of cardiac troponin T, α-myosin heavy chain and α-sarcomeric actinin-expressing cells relative to suboptimum solutions. Cardiomyocytes of composites exhibited connexin43 expression, appropriate contractile kinetics and intracellular calcium handling. Further, adding a modulator of adhesion, thrombospondin-1, abrogated cardiomyocyte differentiation. Thus, the integrated biomaterial platform statistically identified an ECM formulation best supportive of cardiomyocyte differentiation. In future, this formulation could be coupled with biochemical stimulation to improve functional maturation of cardiomyocytes derived in vitro or transplanted in vivo. PMID:26687770
3D prostate MR-TRUS non-rigid registration using dual optimization with volume-preserving constraint
NASA Astrophysics Data System (ADS)
Qiu, Wu; Yuan, Jing; Fenster, Aaron
2016-03-01
We introduce an efficient and novel convex optimization-based approach to the challenging non-rigid registration of 3D prostate magnetic resonance (MR) and transrectal ultrasound (TRUS) images, which incorporates a new volume preserving constraint to essentially improve the accuracy of targeting suspicious regions during the 3D TRUS guided prostate biopsy. Especially, we propose a fast sequential convex optimization scheme to efficiently minimize the employed highly nonlinear image fidelity function using the robust multi-channel modality independent neighborhood descriptor (MIND) across the two modalities of MR and TRUS. The registration accuracy was evaluated using 10 patient images by calculating the target registration error (TRE) using manually identified corresponding intrinsic fiducials in the whole prostate gland. We also compared the MR and TRUS manually segmented prostate surfaces in the registered images in terms of the Dice similarity coefficient (DSC), mean absolute surface distance (MAD), and maximum absolute surface distance (MAXD). Experimental results showed that the proposed method with the introduced volume-preserving prior significantly improves the registration accuracy comparing to the method without the volume-preserving constraint, by yielding an overall mean TRE of 2:0+/-0:7 mm, and an average DSC of 86:5+/-3:5%, MAD of 1:4+/-0:6 mm and MAXD of 6:5+/-3:5 mm.
Sharma, Pramod Kumar; Sharma, Praveen Kumar; Swamidas, Jamema V; Mahantshetty, Umesh; Deshpande, D. D.; Manjhi, Jayanand; Rai, D V
2014-01-01
The aim of this study was to evaluate the dose optimization in 3D image based gynecological interstitial brachytherapy using Martinez Universal Perineal Interstitial Template (MUPIT). Axial CT image data set of 20 patients of gynecological cancer who underwent external radiotherapy and high dose rate (HDR) interstitial brachytherapy using MUPIT was employed to delineate clinical target volume (CTV) and organs at risk (OARs). Geometrical and graphical optimization were done for optimum CTV coverage and sparing of OARs. Coverage Index (CI), dose homogeneity index (DHI), overdose index (OI), dose non-uniformity ratio (DNR), external volume index (EI), conformity index (COIN) and dose volume parameters recommended by GEC-ESTRO were evaluated. The mean CTV, bladder and rectum volume were 137 ± 47cc, 106 ± 41cc and 50 ± 25cc, respectively. Mean CI, DHI and DNR were 0.86 ± 0.03, 0.69 ± 0.11 and 0.31 ± 0.09, while the mean OI, EI, and COIN were 0.08 ± 0.03, 0.07 ± 0.05 and 0.79 ± 0.05, respectively. The estimated mean CTV D90 was 76 ± 11Gy and D100 was 63 ± 9Gy. The different dosimetric parameters of bladder D2cc, D1cc and D0.1cc were 76 ± 11Gy, 81 ± 14Gy, and 98 ± 21Gy and of rectum/recto-sigmoid were 80 ± 17Gy, 85 ± 13Gy, and 124 ± 37Gy, respectively. Dose optimization yields superior coverage with optimal values of indices. Emerging data on 3D image based brachytherapy with reporting and clinical correlation of DVH parameters outcome is enterprizing and provides definite assistance in improving the quality of brachytherapy implants. DVH parameter for urethra in gynecological implants needs to be defined further. PMID:25190999
Michaelsen, Kelly E; Krishnaswamy, Venkataramanan; Shi, Linxi; Vedantham, Srinivasan; Poplack, Steven P; Karellas, Andrew; Pogue, Brian W; Paulsen, Keith D
2015-12-01
Calibration of a three-dimensional multimodal digital breast tomosynthesis (DBT) x-ray and non-fiber based near infrared spectral tomography (NIRST) system is challenging but essential for clinical studies. Phantom imaging results yielded linear contrast recovery of total hemoglobin (HbT) concentration for cylindrical inclusions of 15 mm, 10 mm and 7 mm with a 3.5% decrease in the HbT estimate for each 1 cm increase in inclusion depth. A clinical exam of a patient's breast containing both benign and malignant lesions was successfully imaged, with greater HbT was found in the malignancy relative to the benign abnormality and fibroglandular regions (11 μM vs. 9.5 μM). Tools developed improved imaging system characterization and optimization of signal quality, which will ultimately improve patient selection and subsequent clinical trial results. PMID:26713210
NASA Astrophysics Data System (ADS)
White, Brian; Squires, Todd M.; Hain, Timothy C.; Stone, Howard A.
2003-11-01
Benign paroxysmal positional vertigo (BPPV) is a mechanical disorder of the vestibular system where micron-size crystals abnormally drift into the semicircular canals of the inner ear that sense angular motion of the head. Sedimentation of these crystals causes sensation of motion after true head motion has stopped: vertigo results. The usual clinical treatment is through a series of head maneuvers designed to move the particles into a less sensitive region of the canal system. We present a three-dimensional model to simulate treatment of BPPV by determining the complete hydrodynamic motion of the particles through the course of a therapeutic maneuver while using a realistic representation of the actual geometry. Analyses of clinical maneuvers show the parameter range for which they are effective, and indicate inefficiencies in current practice. In addition, an optimization process determines the most effective head maneuver, which significantly differs from those currently in practice.
3D electromagnetic optimization of the front face of the ITER ICRF antenna
NASA Astrophysics Data System (ADS)
Louche, F.; Dumortier, P.; Messiaen, A.; Durodié, F.
2011-10-01
In the framework of the ion cyclotron resonance heating (ICRH) antenna development for ITER, a design based on an external matching concept has been proposed [1]. We present in this work a series of electromagnetic simulations of this design performed with the commercial code CST Microwave Studio [2]. On the one hand, we explore how various geometrical modifications of some parts of the antenna (the straps and the four-port junction) can practically double the RF power coupled to the plasma. This optimization is supported by transmission line analysis. On the other hand, we treat the important question of the opportunity to tilt the straps in the toroidal direction to follow the plasma curvature as close as possible. We show that a configuration with two toroidal segments is sufficient and that further segmentation is not necessary. This work also underlines significant progress in the realism of ICRH antenna modelling and the importance of considering realistic load shaping in the models.
Research on transformation and optimization of large scale 3D modeling for real time rendering
NASA Astrophysics Data System (ADS)
Yan, Hu; Yang, Yongchao; Zhao, Gang; He, Bin; Shen, Guosheng
2011-12-01
During the simulation process of real-time three-dimensional scene, the popular modeling software and the real-time rendering platform are not compatible. The common solution is to create three-dimensional scene model by using modeling software and then transform the format supported by rendering platform. This paper takes digital campus scene simulation as an example, analyzes and solves the problems of surface loss; texture distortion and loss; model flicker and so on during the transformation from 3Ds Max to MultiGen Creator. Besides, it proposes the optimization strategy of model which is transformed. The operation results show that this strategy is a good solution to all kinds of problems existing in transformation and it can speed up the rendering speed of the model.
Michaelsen, Kelly E.; Krishnaswamy, Venkataramanan; Shi, Linxi; Vedantham, Srinivasan; Poplack, Steven P.; Karellas, Andrew; Pogue, Brian W.; Paulsen, Keith D.
2015-01-01
Calibration of a three-dimensional multimodal digital breast tomosynthesis (DBT) x-ray and non-fiber based near infrared spectral tomography (NIRST) system is challenging but essential for clinical studies. Phantom imaging results yielded linear contrast recovery of total hemoglobin (HbT) concentration for cylindrical inclusions of 15 mm, 10 mm and 7 mm with a 3.5% decrease in the HbT estimate for each 1 cm increase in inclusion depth. A clinical exam of a patient’s breast containing both benign and malignant lesions was successfully imaged, with greater HbT was found in the malignancy relative to the benign abnormality and fibroglandular regions (11 μM vs. 9.5 μM). Tools developed improved imaging system characterization and optimization of signal quality, which will ultimately improve patient selection and subsequent clinical trial results. PMID:26713210
Donner, René; Menze, Bjoern H; Bischof, Horst; Langs, Georg
2013-12-01
The accurate localization of anatomical landmarks is a challenging task, often solved by domain specific approaches. We propose a method for the automatic localization of landmarks in complex, repetitive anatomical structures. The key idea is to combine three steps: (1) a classifier for pre-filtering anatomical landmark positions that (2) are refined through a Hough regression model, together with (3) a parts-based model of the global landmark topology to select the final landmark positions. During training landmarks are annotated in a set of example volumes. A classifier learns local landmark appearance, and Hough regressors are trained to aggregate neighborhood information to a precise landmark coordinate position. A non-parametric geometric model encodes the spatial relationships between the landmarks and derives a topology which connects mutually predictive landmarks. During the global search we classify all voxels in the query volume, and perform regression-based agglomeration of landmark probabilities to highly accurate and specific candidate points at potential landmark locations. We encode the candidates' weights together with the conformity of the connecting edges to the learnt geometric model in a Markov Random Field (MRF). By solving the corresponding discrete optimization problem, the most probable location for each model landmark is found in the query volume. We show that this approach is able to consistently localize the model landmarks despite the complex and repetitive character of the anatomical structures on three challenging data sets (hand radiographs, hand CTs, and whole body CTs), with a median localization error of 0.80 mm, 1.19 mm and 2.71 mm, respectively. PMID:23664450
Cotter, Meghan M.; Whyms, Brian J.; Kelly, Michael P.; Doherty, Benjamin M.; Gentry, Lindell R.; Bersu, Edward T.; Vorperian, Houri K.
2015-01-01
The hyoid bone anchors and supports the vocal tract. Its complex shape is best studied in three dimensions, but it is difficult to capture on computed tomography (CT) images and three-dimensional volume renderings. The goal of this study was to determine the optimal CT scanning and rendering parameters to accurately measure the growth and developmental anatomy of the hyoid and to determine whether it is feasible and necessary to use these parameters in the measurement of hyoids from in vivo CT scans. Direct linear and volumetric measurements of skeletonized hyoid bone specimens were compared to corresponding CT images to determine the most accurate scanning parameters and three-dimensional rendering techniques. A pilot study was undertaken using in vivo scans from a retrospective CT database to determine feasibility of quantifying hyoid growth. Scanning parameters and rendering technique affected accuracy of measurements. Most linear CT measurements were within 10% of direct measurements; however, volume was overestimated when CT scans were acquired with a slice thickness greater than 1.25 mm. Slice-by-slice thresholding of hyoid images decreased volume overestimation. The pilot study revealed that the linear measurements tested correlate with age. A fine-tuned rendering approach applied to small slice thickness CT scans produces the most accurate measurements of hyoid bones. However, linear measurements can be accurately assessed from in vivo CT scans at a larger slice thickness. Such findings imply that investigation into the growth and development of the hyoid bone, and the vocal tract as a whole, can now be performed using these techniques. PMID:25810349
NASA Astrophysics Data System (ADS)
Cheng, Jieyu; Qiu, Wu; Yuan, Jing; Fenster, Aaron; Chiu, Bernard
2016-03-01
Registration of longitudinally acquired 3D ultrasound (US) images plays an important role in monitoring and quantifying progression/regression of carotid atherosclerosis. We introduce an image-based non-rigid registration algorithm to align the baseline 3D carotid US with longitudinal images acquired over several follow-up time points. This algorithm minimizes the sum of absolute intensity differences (SAD) under a variational optical-flow perspective within a multi-scale optimization framework to capture local and global deformations. Outer wall and lumen were segmented manually on each image, and the performance of the registration algorithm was quantified by Dice similarity coefficient (DSC) and mean absolute distance (MAD) of the outer wall and lumen surfaces after registration. In this study, images for 5 subjects were registered initially by rigid registration, followed by the proposed algorithm. Mean DSC generated by the proposed algorithm was 79:3+/-3:8% for lumen and 85:9+/-4:0% for outer wall, compared to 73:9+/-3:4% and 84:7+/-3:2% generated by rigid registration. Mean MAD of 0:46+/-0:08mm and 0:52+/-0:13mm were generated for lumen and outer wall respectively by the proposed algorithm, compared to 0:55+/-0:08mm and 0:54+/-0:11mm generated by rigid registration. The mean registration time of our method per image pair was 143+/-23s.
Tunesi, Marta; Fusco, Federica; Fiordaliso, Fabio; Corbelli, Alessandro; Biella, Gloria; Raimondi, Manuela T.
2016-01-01
Frontotemporal lobar degeneration (FTLD) is a severe neurodegenerative disorder that is diagnosed with increasing frequency in clinical setting. Currently, no therapy is available and in addition the molecular basis of the disease are far from being elucidated. Consequently, it is of pivotal importance to develop reliable and cost-effective in vitro models for basic research purposes and drug screening. To this respect, recent results in the field of Alzheimer’s disease have suggested that a tridimensional (3D) environment is an added value to better model key pathologic features of the disease. Here, we have tried to add complexity to the 3D cell culturing concept by using a microfluidic bioreactor, where cells are cultured under a continuous flow of medium, thus mimicking the interstitial fluid movement that actually perfuses the body tissues, including the brain. We have implemented this model using a neuronal-like cell line (SH-SY5Y), a widely exploited cell model for neurodegenerative disorders that shows some basic features relevant for FTLD modeling, such as the release of the FTLD-related protein progranulin (PRGN) in specific vesicles (exosomes). We have efficiently seeded the cells on 3D scaffolds, optimized a disease-relevant oxidative stress experiment (by targeting mitochondrial function that is one of the possible FTLD-involved pathological mechanisms) and evaluated cell metabolic activity in dynamic culture in comparison to static conditions, finding that SH-SY5Y cells cultured in 3D scaffold are susceptible to the oxidative damage triggered by a mitochondrial-targeting toxin (6-OHDA) and that the same cells cultured in dynamic conditions kept their basic capacity to secrete PRGN in exosomes once recovered from the bioreactor and plated in standard 2D conditions. We think that a further improvement of our microfluidic system may help in providing a full device where assessing basic FTLD-related features (including PRGN dynamic secretion) that may
Neal, Devin; Sakar, Mahmut Selman; Bashir, Rashid; Chan, Vincent; Asada, Haruhiko Harry
2015-06-01
In this study, we present a quantitative approach to construct effective 3D muscle tissues through shape optimization and load impedance matching with electrical and optical stimulation. We have constructed long, thin, fascicle-like skeletal muscle tissue and optimized its form factor through mechanical characterization. A new apparatus was designed and built, which allowed us to measure force-displacement characteristics with diverse load stiffnesses. We have found that (1) there is an optimal form factor that maximizes the muscle stress, (2) the energy transmitted to the load can be maximized with matched load stiffness, and (3) optical stimulation using channelrhodopsin2 in the muscle tissue can generate a twitch force as large as its electrical counterpart for well-developed muscle tissue. Using our tissue construct method, we found that an optimal initial diameter of 500 μm outperformed tissues using 250 μm by more than 60% and tissues using 760 μm by 105%. Using optimal load stiffness, our tissues have generated 12 pJ of energy per twitch at a peak generated stress of 1.28 kPa. Additionally, the difference in optically stimulated twitch performance versus electrically stimulated is a function of how well the overall tissue performs, with average or better performing strips having less than 10% difference. The unique mechanical characterization method used is generalizable to diverse load conditions and will be used to match load impedance to muscle tissue impedance for a wide variety of applications. PMID:25714129
Lao, Zhaoxin; Hu, Yanlei; Zhang, Chenchu; Yang, Liang; Li, Jiawen; Chu, Jiaru; Wu, Dong
2015-12-22
The hierarchical structures are the derivation of various functionalities in the natural world and have inspired broad practical applications in chemical systhesis and biological manipulation. However, traditional top-down fabrication approaches suffered from low complexity. We propose a laser printing capillary-assisted self-assembly (LPCS) strategy for fabricating regular periodic structures. Microscale pillars are first produced by the localized femtosecond laser polymerization and are subsequently self-assembled into periodic hierarchical architectures with the assistance of controlled capillary force. Moreover, based on anisotropic assemblies of micropillars, the LPCS method is further developed for the preparation of more complicated and advanced functional microstructures. Pillars cross section, height, and spatial arrangement can be tuned to guide capillary force, and diverse assemblies with different configurations are thus achieved. Finally, we developed a strategy for growing micro/nanoparticles in designed spatial locations through solution-evaporation self-assembly induced by morphology. Due to the high flexibility of LPCS method, the special arrangements, sizes, and distribution density of the micro/nanoparticles can be controlled readily. Our method will be employed not only to fabricate anisotropic hierarchical structures but also to design and manufacture organic/inorganic microparticles. PMID:26506428
Optimized design of a LED-array-based TOF range imaging sensor for fast 3-D shape measurement
NASA Astrophysics Data System (ADS)
Wang, Huanqin; Wang, Ying; Xu, Jun; He, Deyong; Zhao, Tianpeng; Ming, Hai; Kong, Deyi
2011-06-01
A LED-array-based range imaging sensor using Time-of-Flight (TOF) distance measurement was developed to capture the depth information of three-dimensional (3-D) object. By time-division electronic scanning of the LED heterodyne phase-shift TOF range finders in array, the range images were fast obtained without any mechanical moving parts. The design of LED-array-based range imaging sensor was adequately described and a range imaging theoretical model based on photoelectric signal processing was built, which showed there was mutual restriction relationship among the measurement time of a depth pixel, the bandwidth of receiver and the sensor's signal-to-noise ratio (SNR). In order to improve the key parameters of sensor such as range resolution and measurement speed simultaneously, some optimized designs needed to be done for the proposed range imaging sensor, including choosing proper parameters for the filters in receiver, adopting special structure feedback automatic gain control (AGC) circuit with short response time, etc. The final experiment results showed the sensor after optimization could acquire the range images at a rate of 10 frames per second with a range resolution as high as +/-2mm in the range of 50-1200mm. The essential advantages of the proposed range imaging sensor were construction with simple structure, high range resolution, short measurement time and low cost, which was sufficient for many robotic and industrial automation applications.
Braithwaite, D.; Bourgault, D.; Sulpice, A.; Barbut, J.M.; Tournier, R. l'Universite Joseph Fourier, Grenoble ); Monot, I.; Lepropre, M.; Provost, J.; Desgardin, G. )
1993-04-01
The dc transport critical current densities of melt texture grown and magnetically melt textured bulk YBaCuO have been measured at 77 K and in magnetic fields. A maximum value of over 31,000 A/cm[sup 2] is obtained with a field of 7 teslas applied parallel to the (a,b) planes. Over the rest of the angular range the critical current is shown to be determined mainly by the c-axis component of the applied field. Although this dependency is expected in the presence of two-dimensional vortices, in fact the data are shown to correspond better to the behavior expected of an anisotropic three-dimensional superconductor. These results are compared to magnetization measurements on the same samples. Results show that when the field is directed close to the c-axis, superconducting transport currents flow at fields well above the field at which the irreversible magnetization disappears.
Aganj, Iman; Lenglet, Christophe; Yacoub, Essa; Sapiro, Guillermo; Harel, Noam
2011-01-01
Hardware constraints, scanning time limitations, patient movement, and SNR considerations, restrict the slice-selection and the in-plane resolutions of MRI differently, generally resulting in anisotropic voxels. This non-uniform sampling can be problematic, especially in image segmentation and clinical examination. To alleviate this, the acquisition is divided into (two or) three separate scans, with higher in-plane resolutions and thick slices, yet orthogonal slice-selection directions. In this work, a non-iterative wavelet-based approach for combining the three orthogonal scans is adopted, and its advantages compared to other existing methods, such as Fourier techniques, are discussed, including the consideration of the actual pulse response of the MRI scanner, and its lower computational complexity. Experimental results are shown on simulated and real 7T MRI data. PMID:21761448
NASA Astrophysics Data System (ADS)
Park, Hye-Suk; Kim, Ye-Seul; Lee, Haeng-Hwa; Gang, Won-Suk; Kim, Hee-Joung; Choi, Young-Wook; Choi, JaeGu
2015-08-01
The purpose of this study is to determine the optimal non-uniform angular dose distribution to improve the quality of the 3D reconstructed images and to acquire extra 2D projection images. In this analysis, 7 acquisition sets were generated by using four different values for the number of projections (11, 15, 21, and 29) and total angular range (±14°, ±17.5°, ±21°, and ±24.5° ). For all acquisition sets, the zero-degree projection was used as the 2D image that was close to that of standard conventional mammography (CM). Exposures used were 50, 100, 150, and 200 mR for the zero-degree projection, and the remaining dose was distributed over the remaining projection angles. To quantitatively evaluate image quality, we computed the CNR (contrast-to-noise ratio) and the ASF (artifact spread function) for the same radiation dose. The results indicate that, for microcalcifications, acquisition sets with approximately 4 times higher exposure on the zero-degree projection than the average exposure for the remaining projection angles yielded higher CNR values and were 3% higher than the uniform distribution. However, very high dose concentrations toward the zero-degree projection may reduce the quality of the reconstructed images due to increasing noise in the peripheral views. The zero-degree projection of the non-uniform dose distribution offers a 2D image similar to that of standard CM, but with a significantly lower radiation dose. Therefore, we need to evaluate the diagnostic potential of extra 2D projection image when diagnose breast cancer by using 3D images with non-uniform angular dose distributions.
Lemmo, Stephanie; Atefi, Ehsan; Luker, Gary D.; Tavana, Hossein
2014-01-01
Tumor spheroids are three-dimensional clusters of cancer cells that exhibit characteristics of poorly perfused tumors and hence present a relevant model for testing the efficacy of anti-cancer compounds. The use of spheroids for drug screening is hindered by technological complexities for high throughput generation of consistent size spheroids individually addressable by drug compounds. Here we present and optimize a simple spheroid technology based on the use of an aqueous two-phase system. Cancer cells confined in a drop of the denser aqueous dextran phase are robotically dispensed into a microwell containing the immersion aqueous polyethylene glycol phase. Cells remain within the drop and form a viable spheroid, without a need for any external stimuli. The size of resulting spheroids is sensitive to volume variations of dispensed drops from the air displacement pipetting head of a commercial liquid handling robot. Therefore, we parametrically optimize the process of dispensing of dextran phase drops. For a given cell density, this optimization reproducibly generates consistent size spheroids in standard 96-well plates. In addition, we evaluate the use of a commercial biochemical assay to examine cellular viability of cancer cell spheroids. Spheroids show a dose-dependent response to cisplatin similar to a monolayer culture. However unlike their two-dimensional counterpart, spheroids exhibit resistance to paclitaxel treatment. This technology, which uses only commercially-available reagents and equipment, can potentially expedite anti-cancer drug discovery. Although the use of robotics makes the ATPS spheroid technology particularly useful for drug screening applications, this approach is compatible with simpler liquid handling techniques such as manual micropipetting and offers a straightforward method of 3D cell culture in research laboratories. PMID:25221631
NASA Astrophysics Data System (ADS)
Simpson, J. J.; Taflove, A.
2005-12-01
We report a finite-difference time-domain (FDTD) computational solution of Maxwell's equations [1] that models the possibility of detecting and characterizing ionospheric disturbances above seismic regions. Specifically, we study anomalies in Schumann resonance spectra in the extremely low frequency (ELF) range below 30 Hz as observed in Japan caused by a hypothetical cylindrical ionospheric disturbance above Taiwan. We consider excitation of the global Earth-ionosphere waveguide by lightning in three major thunderstorm regions of the world: Southeast Asia, South America (Amazon region), and Africa. Furthermore, we investigate varying geometries and characteristics of the ionospheric disturbance above Taiwan. The FDTD technique used in this study enables a direct, full-vector, three-dimensional (3-D) time-domain Maxwell's equations calculation of round-the-world ELF propagation accounting for arbitrary horizontal as well as vertical geometrical and electrical inhomogeneities and anisotropies of the excitation, ionosphere, lithosphere, and oceans. Our entire-Earth model grids the annular lithosphere-atmosphere volume within 100 km of sea level, and contains over 6,500,000 grid-points (63 km laterally between adjacent grid points, 5 km radial resolution). We use our recently developed spherical geodesic gridding technique having a spatial discretization best described as resembling the surface of a soccer ball [2]. The grid is comprised entirely of hexagonal cells except for a small fixed number of pentagonal cells needed for completion. Grid-cell areas and locations are optimized to yield a smoothly varying area difference between adjacent cells, thereby maximizing numerical convergence. We compare our calculated results with measured data prior to the Chi-Chi earthquake in Taiwan as reported by Hayakawa et. al. [3]. Acknowledgement This work was suggested by Dr. Masashi Hayakawa, University of Electro-Communications, Chofugaoka, Chofu Tokyo. References [1] A
NASA Astrophysics Data System (ADS)
Gosálvez, M. A.; Pal, Prem; Ferrando, N.; Hida, H.; Sato, K.
2011-12-01
This is part I of a series of two papers dedicated to the presentation of a novel, large throughput, experimental procedure to determine the three-dimensional distribution of the etch rate of silicon in a wide range of anisotropic etchants, including a total of 30 different etching conditions in KOH, KOH+IPA, TMAH and TMAH+Triton solutions at various concentrations and temperatures. The method is based on the use of previously reported, vertically micromachined wagon wheels (WWs) (Wind and Hines 2000 Surf. Sci. 460 21-38 Nguyen and Elwenspoek 2007 J. Electrochem. Soc. 154 D684-91), focusing on speeding up the etch rate extraction process for each WW by combining macrophotography and image processing procedures. The proposed procedure positions the WWs as a realistic alternative to the traditional hemispherical specimen. The obtained, extensive etch rate database is used to perform wet etching simulations of advanced systems, showing good agreement with the experimental counterparts. In part II of this series (Gosálvez et al J. Micromech. Microeng. 21 125008), we provide a theoretical analysis of the etched spoke shapes, a detailed comparison to the etch rates from previous studies and a self-consistency study of the measured etch rates against maximum theoretical values derived from the spoke shape analysis.
Baino, Francesco; Ferraris, Monica; Bretcanu, Oana; Verné, Enrica; Vitale-Brovarone, Chiara
2013-03-01
Fabrication of 3-D highly porous, bioactive, and mechanically competent scaffolds represents a significant challenge of bone tissue engineering. In this work, Bioglass®-derived glass-ceramic scaffolds actually fulfilling this complex set of requirements were successfully produced through the sponge replication method. Scaffold processing parameters and sintering treatment were carefully designed in order to obtain final porous bodies with pore content (porosity above 70 %vol), trabecular architecture and mechanical properties (compressive strength up to 3 MPa) analogous to those of the cancellous bone. Influence of the Bioglass® particles size on the structural and mechanical features of the sintered scaffolds was considered and discussed. Relationship between porosity and mechanical strength was investigated and modeled. Three-dimensional architecture, porosity, mechanical strength and in vitro bioactivity of the optimized Bioglass®-derived scaffolds were also compared to those of CEL2-based glass-ceramic scaffolds (CEL2 is an experimental bioactive glass originally developed by the authors at Politecnico di Torino) fabricated by the same processing technique, in an attempt at understanding the role of different bioactive glass composition on the major features of scaffolds prepared by the same method. PMID:22207602
Highly optimized simulations on single- and multi-GPU systems of the 3D Ising spin glass model
NASA Astrophysics Data System (ADS)
Lulli, M.; Bernaschi, M.; Parisi, G.
2015-11-01
We present a highly optimized implementation of a Monte Carlo (MC) simulator for the three-dimensional Ising spin-glass model with bimodal disorder, i.e., the 3D Edwards-Anderson model running on CUDA enabled GPUs. Multi-GPU systems exchange data by means of the Message Passing Interface (MPI). The chosen MC dynamics is the classic Metropolis one, which is purely dissipative, since the aim was the study of the critical off-equilibrium relaxation of the system. We focused on the following issues: (i) the implementation of efficient memory access patterns for nearest neighbours in a cubic stencil and for lagged-Fibonacci-like pseudo-Random Numbers Generators (PRNGs); (ii) a novel implementation of the asynchronous multispin-coding Metropolis MC step allowing to store one spin per bit and (iii) a multi-GPU version based on a combination of MPI and CUDA streams. Cubic stencils and PRNGs are two subjects of very general interest because of their widespread use in many simulation codes.
NASA Astrophysics Data System (ADS)
Gotchev, Atanas; Smolic, Aljoscha; Jumisko-Pyykkö, Satu; Strohmeier, Dominik; Bozdagi Akar, Gozde; Merkle, Philipp; Daskalov, Nikolai
2009-02-01
A European consortium of six partners has been developing core technological components of a mobile 3D television system over DVB-H channel. In this overview paper, we present our current results on developing optimal methods for stereo-video content creation, coding and transmission and emphasize their significance for the power-constrained mobile platform, equipped with auto-stereoscopic display. We address the user requirements by applying modern usercentered approaches taking into account different user groups and usage contexts in contrast to the laboratory assessment methods which, though standardized, offer limited applicability to real applications. To this end, we have been aiming at developing a methodological framework for the whole system development process. One of our goals has been to further develop the user-centered approach towards experienced quality of critical system components. In this paper, we classify different research methods and technological solutions analyzing their pros and constraints. Based on this analysis we present the user-centered methodological framework being used throughout the whole development process of the system and aimed at achieving the best performance and quality appealing to the end user.
Combinatorial 3D Mechanical Metamaterials
NASA Astrophysics Data System (ADS)
Coulais, Corentin; Teomy, Eial; de Reus, Koen; Shokef, Yair; van Hecke, Martin
2015-03-01
We present a class of elastic structures which exhibit 3D-folding motion. Our structures consist of cubic lattices of anisotropic unit cells that can be tiled in a complex combinatorial fashion. We design and 3d-print this complex ordered mechanism, in which we combine elastic hinges and defects to tailor the mechanics of the material. Finally, we use this large design space to encode smart functionalities such as surface patterning and multistability.
NASA Astrophysics Data System (ADS)
Castaldo, Raffaele; De Novellis, Vincenzo; Lollino, Piernicola; Manunta, Michele; Tizzani, Pietro
2015-04-01
The new challenge that the research in slopes instabilities phenomena is going to tackle is the effective integration and joint exploitation of remote sensing measurements with in situ data and observations to study and understand the sub-surface interactions, the triggering causes, and, in general, the long term behaviour of the investigated landslide phenomenon. In this context, a very promising approach is represented by Finite Element (FE) techniques, which allow us to consider the intrinsic complexity of the mass movement phenomena and to effectively benefit from multi source observations and data. In this context, we perform a three dimensional (3D) numerical model of the Ivancich (Assisi, Central Italy) instability phenomenon. In particular, we apply an inverse FE method based on a Genetic Algorithm optimization procedure, benefitting from advanced DInSAR measurements, retrieved through the full resolution Small Baseline Subset (SBAS) technique, and an inclinometric array distribution. To this purpose we consider the SAR images acquired from descending orbit by the COSMO-SkyMed (CSK) X-band radar constellation, from December 2009 to February 2012. Moreover the optimization input dataset is completed by an array of eleven inclinometer measurements, from 1999 to 2006, distributed along the unstable mass. The landslide body is formed of debris material sliding on a arenaceous marl substratum, with a thin shear band detected using borehole and inclinometric data, at depth ranging from 20 to 60 m. Specifically, we consider the active role of this shear band in the control of the landslide evolution process. A large field monitoring dataset of the landslide process, including at-depth piezometric and geological borehole observations, were available. The integration of these datasets allows us to develop a 3D structural geological model of the considered slope. To investigate the dynamic evolution of a landslide, various physical approaches can be considered
Zhao, Yi; Cao, Xiangyu; Gao, Jun; Sun, Yu; Yang, Huanhuan; Liu, Xiao; Zhou, Yulong; Han, Tong; Chen, Wei
2016-01-01
We propose a new strategy to design broadband and wide angle diffusion metasurfaces. An anisotropic structure which has opposite phases under x- and y-polarized incidence is employed as the "0" and "1" elements base on the concept of coding metamaterial. To obtain a uniform backward scattering under normal incidence, Simulated Annealing algorithm is utilized in this paper to calculate the optimal layout. The proposed method provides an efficient way to design diffusion metasurface with a simple structure, which has been proved by both simulations and measurements. PMID:27034110
NASA Astrophysics Data System (ADS)
Zhao, Yi; Cao, Xiangyu; Gao, Jun; Sun, Yu; Yang, Huanhuan; Liu, Xiao; Zhou, Yulong; Han, Tong; Chen, Wei
2016-04-01
We propose a new strategy to design broadband and wide angle diffusion metasurfaces. An anisotropic structure which has opposite phases under x- and y-polarized incidence is employed as the “0” and “1” elements base on the concept of coding metamaterial. To obtain a uniform backward scattering under normal incidence, Simulated Annealing algorithm is utilized in this paper to calculate the optimal layout. The proposed method provides an efficient way to design diffusion metasurface with a simple structure, which has been proved by both simulations and measurements.
Zhao, Yi; Cao, Xiangyu; Gao, Jun; Sun, Yu; Yang, Huanhuan; Liu, Xiao; Zhou, Yulong; Han, Tong; Chen, Wei
2016-01-01
We propose a new strategy to design broadband and wide angle diffusion metasurfaces. An anisotropic structure which has opposite phases under x- and y-polarized incidence is employed as the “0” and “1” elements base on the concept of coding metamaterial. To obtain a uniform backward scattering under normal incidence, Simulated Annealing algorithm is utilized in this paper to calculate the optimal layout. The proposed method provides an efficient way to design diffusion metasurface with a simple structure, which has been proved by both simulations and measurements. PMID:27034110
The Tonga-Vanuatu Subduction Complex -- a Self-Optimized 3D Slab-Slab-Mantle Heat Pump
NASA Astrophysics Data System (ADS)
McCreary, J. A.
2008-12-01
Recently published geophysical and geochemical data and increasingly actualistic free subduction models prompted a fresh look at 2 classics hinting, in combination, that a coupled 3D slab-slab-upper mantle interaction (Scholz and Campos, 1995; full citations at URL below) might power the prodigious surface heat dissipation (Lagabrielle et al., 1997) characterizing one of Earth's most remarkable tectonomagmatic systems, the Tonga-Vanuatu Subduction Complex (TVSC). The 3D TVSC includes (1) the kinematically, magmatically, and bathymetrically distinct North Tonga (NT, 14-26° S) and South Vanuatu (SV, 16-23° S) trenches and slabs, (2) the shared NT-SV backarc, and (3) entrained mobile upper mantle (MUM). That Earth's greatest convergence, rollback, and spreading rates; most disseminated spreading (the North Fiji Basin (NFB) ridge swarm); and greatest concentration of aggregate active ridge length coincide in a 1,500 km TVSC can't be accidental. To the north and south, the respective active NT and SV trenches swing abruptly 90° counterclockwise into continuity with the Vitiaz and Hunter fossil trenches, both active in the Late Miocene but now sinistral strike-slip loci standing over long exposed PA and AU slab edges. These 2 active-fossil trench pairs bracket a hot, shallow and geophysically and geochemically exceptional TVSC interior consisting of 2 rapidly spreading backarcs set back-to-back in free sublithospheric communication: The Lau-Havre NT backarc on the east and the ridge-infested SV backarc (NFB) on the west. The NFB and adjacent North Fiji Plateau make up the unplatelike New Hebrides-Fiji Orogen (Bird, 2003). As in the western Aleutians, the NT-Vitiaz and SV-Hunter subduction-to-strike-slip transitions (SSSTs) stand above toroidal fluxes of hot, dry PA and AU MUM driven along-trench and around the free NT and SV slab edges from subslab to supraslab regions by dynamic pressure gradients powered by slab free-fall and induced viscous couplings. These edge
3D Ultrasonic Wave Simulations for Structural Health Monitoring
NASA Technical Reports Server (NTRS)
Campbell, Leckey Cara A/; Miler, Corey A.; Hinders, Mark K.
2011-01-01
Structural health monitoring (SHM) for the detection of damage in aerospace materials is an important area of research at NASA. Ultrasonic guided Lamb waves are a promising SHM damage detection technique since the waves can propagate long distances. For complicated flaw geometries experimental signals can be difficult to interpret. High performance computing can now handle full 3-dimensional (3D) simulations of elastic wave propagation in materials. We have developed and implemented parallel 3D elastodynamic finite integration technique (3D EFIT) code to investigate ultrasound scattering from flaws in materials. EFIT results have been compared to experimental data and the simulations provide unique insight into details of the wave behavior. This type of insight is useful for developing optimized experimental SHM techniques. 3D EFIT can also be expanded to model wave propagation and scattering in anisotropic composite materials.
NASA Astrophysics Data System (ADS)
Ali, Khuram; Khan, Sohail Aziz; Jafri, Mohd Zubir Mat
2012-10-01
Reduced surface reflectance and enhanced light trapping is required by any high efficiency solar cell. Anisotropic etching was done on silicon (1 0 0) by using tetramethyl ammonium hydroxide TMAH, (CH3)4NOH, solution at 85 °C. Process variables considered were solution concentration and time proposed by response surface methodology (RSM). An effective surface texture was resulted with reflectance less than 8% without antireflection coating. The antireflection mechanism was also co-related with the etch rate of Si. Optimized values predicted by RSM for time and TMAH concentration were 5 min and 3.50% respectively. The technique and optimization of parameters by using response surface methodology (RSM) could be valuable in the texturization process for high-efficiency Si solar cells.
Mijnheer, B; Mans, A; Olaciregui-Ruiz, I; Rozendaal, R; Spreeuw, H; Herk, M van
2014-06-15
Purpose: To develop a 3D in vivo dosimetry method that is able to substitute pre-treatment verification in an efficient way, and to terminate treatment delivery if the online measured 3D dose distribution deviates too much from the predicted dose distribution. Methods: A back-projection algorithm has been further developed and implemented to enable automatic 3D in vivo dose verification of IMRT/VMAT treatments using a-Si EPIDs. New software tools were clinically introduced to allow automated image acquisition, to periodically inspect the record-and-verify database, and to automatically run the EPID dosimetry software. The comparison of the EPID-reconstructed and planned dose distribution is done offline to raise automatically alerts and to schedule actions when deviations are detected. Furthermore, a software package for online dose reconstruction was also developed. The RMS of the difference between the cumulative planned and reconstructed 3D dose distributions was used for triggering a halt of a linac. Results: The implementation of fully automated 3D EPID-based in vivo dosimetry was able to replace pre-treatment verification for more than 90% of the patient treatments. The process has been fully automated and integrated in our clinical workflow where over 3,500 IMRT/VMAT treatments are verified each year. By optimizing the dose reconstruction algorithm and the I/O performance, the delivered 3D dose distribution is verified in less than 200 ms per portal image, which includes the comparison between the reconstructed and planned dose distribution. In this way it was possible to generate a trigger that can stop the irradiation at less than 20 cGy after introducing large delivery errors. Conclusion: The automatic offline solution facilitated the large scale clinical implementation of 3D EPID-based in vivo dose verification of IMRT/VMAT treatments; the online approach has been successfully tested for various severe delivery errors.
NASA Astrophysics Data System (ADS)
Bremner, P. M.; Panning, M. P.; Russo, R.; Mocanu, V. I.; Stanciu, A. C.; Torpey, M. E.; Hongsresawat, S.; VanDecar, J. C.
2015-12-01
We present new 3-D radially anisotropic and isotropic crustal velocity models beneath central Idaho and eastern Oregon. We produced the velocity models from Love and horizontal component Rayleigh wave group and phase velocity measurements on the IDaho/ORegon (IDOR) Passive seismic network, 86 broadband seismic stations, dataset using ambient noise tomography and the methods of Gallego et. al (2010) and Lin et. al (2008). We calculated inter-station group/phase velocities in narrow frequency bands from travel-time measurements of the rotated stacked horizontal component cross-correlations (bandpass filtered between 2 and 30 seconds), which we used to invert for velocity structure beneath the network. We derived group and phase velocity maps for each frequency band using the damped least-squares inversion method of Tarantola (2005), and then jointly inverted for velocity with depth. Moho depths are prescribed in the joint inversions based on receiver functions, also from the IDOR seismic data, and provides a starting crustal velocity model. Goals of our work include refining models of crustal structure in the accreted Blue Mountain terranes in the western study area; determining the depth extent of the Salmon River Suture/West Idaho Shear Zone (WISZ), which crosses north-south through the middle of the network; determining the architecture of the Idaho batholith, an extensive largely crustal-derived pluton; and examining the nature of the autochthonous (?) North American crust and lithosphere beneath and east of the batholith.
Filippi, Andrea Riccardo; Ragona, Riccardo; Piva, Cristina; Scafa, Davide; Fiandra, Christian; Fusella, Marco; Giglioli, Francesca Romana; Lohr, Frank; Ricardi, Umberto
2015-05-01
Purpose: The purpose of this study was to evaluate the risks of second cancers and cardiovascular diseases associated with an optimized volumetric modulated arc therapy (VMAT) planning solution in a selected cohort of stage I/II Hodgkin lymphoma (HL) patients treated with either involved-node or involved-site radiation therapy in comparison with 3-dimensional conformal radiation therapy (3D-CRT). Methods and Materials: Thirty-eight patients (13 males and 25 females) were included. Disease extent was mediastinum alone (n=8, 21.1%); mediastinum plus unilateral neck (n=19, 50%); mediastinum plus bilateral neck (n=11, 29.9%). Prescription dose was 30 Gy in 2-Gy fractions. Only 5 patients had mediastinal bulky disease at diagnosis (13.1%). Anteroposterior 3D-CRT was compared with a multiarc optimized VMAT solution. Lung, breast, and thyroid cancer risks were estimated by calculating a lifetime attributable risk (LAR), with a LAR ratio (LAR{sub VMAT}-to-LAR{sub 3D-CRT}) as a comparative measure. Cardiac toxicity risks were estimated by calculating absolute excess risk (AER). Results: The LAR ratio favored 3D-CRT for lung cancer induction risk in mediastinal alone (P=.004) and mediastinal plus unilateral neck (P=.02) presentations. LAR ratio for breast cancer was lower for VMAT in mediastinal plus bilateral neck presentations (P=.02), without differences for other sites. For thyroid cancer, no significant differences were observed, regardless of anatomical presentation. A significantly lower AER of cardiac (P=.038) and valvular diseases (P<.0001) was observed for VMAT regardless of disease extent. Conclusions: In a cohort of patients with favorable characteristics in terms of disease extent at diagnosis (large prevalence of nonbulky presentations without axillary involvement), optimized VMAT reduced heart disease risk with comparable risks of thyroid and breast cancer, with an increase in lung cancer induction probability. The results are however strongly influenced by
Unaldi, Numan; Temel, Samil; Asari, Vijayan K
2012-01-01
One of the most critical issues of Wireless Sensor Networks (WSNs) is the deployment of a limited number of sensors in order to achieve maximum coverage on a terrain. The optimal sensor deployment which enables one to minimize the consumed energy, communication time and manpower for the maintenance of the network has attracted interest with the increased number of studies conducted on the subject in the last decade. Most of the studies in the literature today are proposed for two dimensional (2D) surfaces; however, real world sensor deployments often arise on three dimensional (3D) environments. In this paper, a guided wavelet transform (WT) based deployment strategy (WTDS) for 3D terrains, in which the sensor movements are carried out within the mutation phase of the genetic algorithms (GAs) is proposed. The proposed algorithm aims to maximize the Quality of Coverage (QoC) of a WSN via deploying a limited number of sensors on a 3D surface by utilizing a probabilistic sensing model and the Bresenham's line of sight (LOS) algorithm. In addition, the method followed in this paper is novel to the literature and the performance of the proposed algorithm is compared with the Delaunay Triangulation (DT) method as well as a standard genetic algorithm based method and the results reveal that the proposed method is a more powerful and more successful method for sensor deployment on 3D terrains. PMID:22666078
Unaldi, Numan; Temel, Samil; Asari, Vijayan K.
2012-01-01
One of the most critical issues of Wireless Sensor Networks (WSNs) is the deployment of a limited number of sensors in order to achieve maximum coverage on a terrain. The optimal sensor deployment which enables one to minimize the consumed energy, communication time and manpower for the maintenance of the network has attracted interest with the increased number of studies conducted on the subject in the last decade. Most of the studies in the literature today are proposed for two dimensional (2D) surfaces; however, real world sensor deployments often arise on three dimensional (3D) environments. In this paper, a guided wavelet transform (WT) based deployment strategy (WTDS) for 3D terrains, in which the sensor movements are carried out within the mutation phase of the genetic algorithms (GAs) is proposed. The proposed algorithm aims to maximize the Quality of Coverage (QoC) of a WSN via deploying a limited number of sensors on a 3D surface by utilizing a probabilistic sensing model and the Bresenham's line of sight (LOS) algorithm. In addition, the method followed in this paper is novel to the literature and the performance of the proposed algorithm is compared with the Delaunay Triangulation (DT) method as well as a standard genetic algorithm based method and the results reveal that the proposed method is a more powerful and more successful method for sensor deployment on 3D terrains. PMID:22666078
3D braid scaffolds for regeneration of articular cartilage.
Ahn, Hyunchul; Kim, Kyoung Ju; Park, Sook Young; Huh, Jeong Eun; Kim, Hyun Jeong; Yu, Woong-Ryeol
2014-06-01
Regenerating articular cartilage in vivo from cultured chondrocytes requires that the cells be cultured and implanted within a biocompatible, biodegradable scaffold. Such scaffolds must be mechanically stable; otherwise chondrocytes would not be supported and patients would experience severe pain. Here we report a new 3D braid scaffold that matches the anisotropic (gradient) mechanical properties of natural articular cartilage and is permissive to cell cultivation. To design an optimal structure, the scaffold unit cell was mathematically modeled and imported into finite element analysis. Based on this analysis, a 3D braid structure with gradient axial yarn distribution was designed and manufactured using a custom-built braiding machine. The mechanical properties of the 3D braid scaffold were evaluated and compared with simulated results, demonstrating that a multi-scale approach consisting of unit cell modeling and continuum analysis facilitates design of scaffolds that meet the requirements for mechanical compatibility with tissues. PMID:24556323
NASA Astrophysics Data System (ADS)
Wilkie, William Keats
1997-12-01
An aeroelastic model suitable for control law and preliminary structural design of composite helicopter rotor blades incorporating embedded anisotropic piezoelectric actuator laminae is developed. The aeroelasticity model consists of a linear, nonuniform beam representation of the blade structure, including linear piezoelectric actuation terms, coupled with a nonlinear, finite-state unsteady aerodynamics model. A Galerkin procedure and numerical integration in the time domain are used to obtain a soluti An aeroelastic model suitable for control law and preliminary structural design of composite helicopter rotor blades incorporating embedded anisotropic piezoelectric actuator laminae is developed. The aeroelasticity model consists of a linear, nonuniform beam representation of the blade structure, including linear piezoelectric actuation terms, coupled with a nonlinear, finite-state unsteady aerodynamics model. A Galerkin procedure and numerical integration in the time domain are used to obtain amited additional piezoelectric material mass, it is shown that blade twist actuation approaches which exploit in-plane piezoelectric free-stain anisotropies are capable of producing amplitudes of oscillatory blade twisting sufficient for rotor vibration reduction applications. The second study examines the effectiveness of using embedded piezoelectric actuator laminae to alleviate vibratory loads due to retreating blade stall. A 10 to 15 percent improvement in dynamic stall limited forward flight speed, and a 5 percent improvement in stall limited rotor thrust were numerically demonstrated for the active twist rotor blade relative to a conventional blade design. The active twist blades are also demonstrated to be more susceptible than the conventional blades to dynamic stall induced vibratory loads when not operating with twist actuation. This is the result of designing the active twist blades with low torsional stiffness in order to maximize piezoelectric twist authority
Ayad, G.; Barriere, T.; Gelin, J. C.; Liu, B.
2007-05-17
The paper is concerned with optimization and parametric identification of Powder Injection Molding process that consists first in injection of powder mixture with polymer binder and then to the sintering of the resulting powders parts by solid state diffusion. In the first part, one describes an original methodology to optimize the injection stage based on the combination of Design Of Experiments and an adaptive Response Surface Modeling. Then the second part of the paper describes the identification strategy that one proposes for the sintering stage, using the identification of sintering parameters from dilatometer curves followed by the optimization of the sintering process. The proposed approaches are applied to the optimization for manufacturing of a ceramic femoral implant. One demonstrates that the proposed approach give satisfactory results.
KNUPP,PATRICK; MITCHELL,SCOTT A.
1999-11-01
In an attempt to automatically produce high-quality all-hex meshes, we investigated a mesh improvement strategy: given an initial poor-quality all-hex mesh, we iteratively changed the element connectivity, adding and deleting elements and nodes, and optimized the node positions. We found a set of hex reconnection primitives. We improved the optimization algorithms so they can untangle a negative-Jacobian mesh, even considering Jacobians on the boundary, and subsequently optimize the condition number of elements in an untangled mesh. However, even after applying both the primitives and optimization we were unable to produce high-quality meshes in certain regions. Our experiences suggest that many boundary configurations of quadrilaterals admit no hexahedral mesh with positive Jacobians, although we have no proof of this.
2009-01-01
Background Discovery of new bioactive molecules that could enter drug discovery programs or that could serve as chemical probes is a very complex and costly endeavor. Structure-based and ligand-based in silico screening approaches are nowadays extensively used to complement experimental screening approaches in order to increase the effectiveness of the process and facilitating the screening of thousands or millions of small molecules against a biomolecular target. Both in silico screening methods require as input a suitable chemical compound collection and most often the 3D structure of the small molecules has to be generated since compounds are usually delivered in 1D SMILES, CANSMILES or in 2D SDF formats. Results Here, we describe the new open source program DG-AMMOS which allows the generation of the 3D conformation of small molecules using Distance Geometry and their energy minimization via Automated Molecular Mechanics Optimization. The program is validated on the Astex dataset, the ChemBridge Diversity database and on a number of small molecules with known crystal structures extracted from the Cambridge Structural Database. A comparison with the free program Balloon and the well-known commercial program Omega generating the 3D of small molecules is carried out. The results show that the new free program DG-AMMOS is a very efficient 3D structure generator engine. Conclusion DG-AMMOS provides fast, automated and reliable access to the generation of 3D conformation of small molecules and facilitates the preparation of a compound collection prior to high-throughput virtual screening computations. The validation of DG-AMMOS on several different datasets proves that generated structures are generally of equal quality or sometimes better than structures obtained by other tested methods. PMID:19912625
NASA Astrophysics Data System (ADS)
Wang, Jian; Cheng, Shuang; Li, Wanfei; Zhang, Su; Li, Hongfei; Zheng, Zhaozhao; Li, Fujin; Shi, Liyi; Lin, Hongzhen; Zhang, Yuegang
2016-07-01
Lithium/sulfur (Li/S) battery is a promising next-generation energy storage system owing to its high theoretical energy density. However, for practical use there remains some key problems to be solved, such as low active material utilization and rapid capacity fading, especially at high areal sulfur loadings. Here, we report a facile one-pot method to prepare porous three-dimensional nitrogen, sulfur-codoped graphene through hydrothermal reduction of graphene oxide with multi-ion mixture modulation. We show solid evidence that the results of multi-ion mixture modulation can not only improve the surface affinity of the nanocarbons to polysulfides, but also alter their assembling manner and render the resultant 3D network a more favorable pore morphology for accommodating and confining sulfur. It also had an excellent rate performance and cycling stability, showing an initial capacity of 1304 mA h g-1 at 0.05C, 613 mA h g-1 at 5C and maintaining a reversible capacity of 462 mA h g-1 after 1500 cycles at 2C with capacity fading as low as 0.028% per cycle. Moreover, a high areal capacity of 5.1 mA h cm-2 at 0.2C is achieved at an areal sulfur loading of 6.3 mg cm-2, which are the best values reported so far for dual-doped sulfur cathodes.
Imani, Rana; Hojjati Emami, Shahriar; Fakhrzadeh, Hossein; Baheiraei, Nafiseh; Sharifi, Ali M
2012-04-01
The ultimate goal of tissue engineering is to design and fabricate functional human tissues that are similar to natural cells and are capable of regeneration. Preparation of cell aggregates is one of the important steps in 3D tissue engineering technology, particularly in organ printing. Two simple methods, hanging drop (HD) and conical tube (CT) were utilized to prepare cell aggregates. The size and viability of the aggregates obtained at different initial cell densities and pre-culture duration were compared. The proliferative ability of the cell aggregates and their ability to spread in culture plates were also investigated. In both methods, the optimum average size of the aggregates was less than 500 microm. CT aggregates were smaller than HD aggregates. 5,000 cells per drop HD aggregates showed a marked ability to attach and spread on the culture surface. The proliferative ability reduced when the initial cell density was increased. Comparing these methods, we found that the HD method having better size controlling ability as well as enhanced ability to maintain higher rates of viability, spreading, and proliferation. In conclusion, smaller HD aggregates might be a suitable choice as building blocks for making bioink particles in bioprinting technique. PMID:23173303
Zarpalas, Dimitrios; Gkontra, Polyxeni; Daras, Petros; Maglaveras, Nicos
2014-01-01
Assessing the structural integrity of the hippocampus (HC) is an essential step toward prevention, diagnosis, and follow-up of various brain disorders due to the implication of the structural changes of the HC in those disorders. In this respect, the development of automatic segmentation methods that can accurately, reliably, and reproducibly segment the HC has attracted considerable attention over the past decades. This paper presents an innovative 3-D fully automatic method to be used on top of the multiatlas concept for the HC segmentation. The method is based on a subject-specific set of 3-D optimal local maps (OLMs) that locally control the influence of each energy term of a hybrid active contour model (ACM). The complete set of the OLMs for a set of training images is defined simultaneously via an optimization scheme. At the same time, the optimal ACM parameters are also calculated. Therefore, heuristic parameter fine-tuning is not required. Training OLMs are subsequently combined, by applying an extended multiatlas concept, to produce the OLMs that are anatomically more suitable to the test image. The proposed algorithm was tested on three different and publicly available data sets. Its accuracy was compared with that of state-of-the-art methods demonstrating the efficacy and robustness of the proposed method. PMID:27170866
NASA Astrophysics Data System (ADS)
Zhang, Honghai; Abiose, Ademola K.; Campbell, Dwayne N.; Sonka, Milan; Martins, James B.; Wahle, Andreas
2010-03-01
Quantitative analysis of the left ventricular shape and motion patterns associated with left ventricular mechanical dyssynchrony (LVMD) is essential for diagnosis and treatment planning in congestive heart failure. Real-time 3D echocardiography (RT3DE) used for LVMD analysis is frequently limited by heavy speckle noise or partially incomplete data, thus a segmentation method utilizing learned global shape knowledge is beneficial. In this study, the endocardial surface of the left ventricle (LV) is segmented using a hybrid approach combining active shape model (ASM) with optimal graph search. The latter is used to achieve landmark refinement in the ASM framework. Optimal graph search translates the 3D segmentation into the detection of a minimum-cost closed set in a graph and can produce a globally optimal result. Various information-gradient, intensity distributions, and regional-property terms-are used to define the costs for the graph search. The developed method was tested on 44 RT3DE datasets acquired from 26 LVMD patients. The segmentation accuracy was assessed by surface positioning error and volume overlap measured for the whole LV as well as 16 standard LV regions. The segmentation produced very good results that were not achievable using ASM or graph search alone.
Gkontra, Polyxeni; Daras, Petros; Maglaveras, Nicos
2014-01-01
Assessing the structural integrity of the hippocampus (HC) is an essential step toward prevention, diagnosis, and follow-up of various brain disorders due to the implication of the structural changes of the HC in those disorders. In this respect, the development of automatic segmentation methods that can accurately, reliably, and reproducibly segment the HC has attracted considerable attention over the past decades. This paper presents an innovative 3-D fully automatic method to be used on top of the multiatlas concept for the HC segmentation. The method is based on a subject-specific set of 3-D optimal local maps (OLMs) that locally control the influence of each energy term of a hybrid active contour model (ACM). The complete set of the OLMs for a set of training images is defined simultaneously via an optimization scheme. At the same time, the optimal ACM parameters are also calculated. Therefore, heuristic parameter fine-tuning is not required. Training OLMs are subsequently combined, by applying an extended multiatlas concept, to produce the OLMs that are anatomically more suitable to the test image. The proposed algorithm was tested on three different and publicly available data sets. Its accuracy was compared with that of state-of-the-art methods demonstrating the efficacy and robustness of the proposed method. PMID:27170866
NASA Astrophysics Data System (ADS)
Lee, Anice; Lin, Chung-Yi; Chen, F. Y.; Chang, Wen-Hao; Hsu, Sean; Li, Allen; Luo, Ying; Wen, Youxian
2009-03-01
To fully characterize the lithography process, it is critical to have accurate CD and profile of photo resist structure at ADI stage. Traditionally, CDSEM can only provide limited profile information, and is extremely challenged to be integrated for real-time in line wafer level process control because of throughput issue. Over the past few years, optical digital profilometry (ODP(R)) developed by Timbre Technologies, Inc., has been adopted for real-time process control in Litho for optical CD and shape monitoring. In this paper, the integrated ODP(R) reflectometer is applied to study process signatures of 3D complicated ADI and AEI structures of a 70nm DRAM process. The DT structures from the 70nm node process studied in this paper, are elliptical photo resist via developed over a thin film stack at ADI, and via etched deeply through the thick (over 5um) dielectric film at AEI. At ADI, the ODP(R) library is qualified by careful cross check with CDSEM data. CD results from iODP(R) show very good correlation to that from CDSEM. The iODPÂ® measurement for a FEM wafer shows smoother and cleaner Bossung curves than the CDSEM does. At AEI, the library is then qualified for top CD measurement in comparison to CDSEM, and also to results at ADI. With implementing iODP(R) measurements into both ADI and AEI structures, their signature patterns from ADI to AEI for 70nm DT process can be matched successfully. Such a signature pattern match indicates the strong correlation between ADI and AEI processes and can be fully made use of for APC. It is a significant step toward IM APC control considering ADI and AEI process steps together.
NASA Astrophysics Data System (ADS)
Taccogna, F.; Minelli, P.; Cavenago, M.; Veltri, P.; Ippolito, N.
2016-02-01
The geometry of a single aperture in the extraction grid plays a relevant role for the optimization of negative ion transport and extraction probability in a hybrid negative ion source. For this reason, a three-dimensional particle-in-cell/Monte Carlo collision model of the extraction region around the single aperture including part of the source and part of the acceleration (up to the extraction grid (EG) middle) regions has been developed for the new aperture design prepared for negative ion optimization 1 source. Results have shown that the dimension of the flat and chamfered parts and the slope of the latter in front of the source region maximize the product of production rate and extraction probability (allowing the best EG field penetration) of surface-produced negative ions. The negative ion density in the plane yz has been reported.
Taccogna, F; Minelli, P; Cavenago, M; Veltri, P; Ippolito, N
2016-02-01
The geometry of a single aperture in the extraction grid plays a relevant role for the optimization of negative ion transport and extraction probability in a hybrid negative ion source. For this reason, a three-dimensional particle-in-cell/Monte Carlo collision model of the extraction region around the single aperture including part of the source and part of the acceleration (up to the extraction grid (EG) middle) regions has been developed for the new aperture design prepared for negative ion optimization 1 source. Results have shown that the dimension of the flat and chamfered parts and the slope of the latter in front of the source region maximize the product of production rate and extraction probability (allowing the best EG field penetration) of surface-produced negative ions. The negative ion density in the plane yz has been reported. PMID:26932027
NASA Astrophysics Data System (ADS)
Åsell, Mats; Hyödynmaa, Simo; Gustafsson, Anders; Brahme, Anders
1997-11-01
The possibilities of using simultaneous fluence and energy modulation techniques in electron beam therapy to shape the dose distribution and almost eliminate the influences of tissue inhomogeneities have been investigated. By using a radiobiologically based optimization algorithm the radiobiological properties of the tissues can be taken into account when trying to find the best possible dose delivery. First water phantoms with differently shaped surfaces were used to study the effect of surface irregularities. We also studied water phantoms with internal inhomogeneities consisting of air or cortical bone. It was possible to improve substantially the dose distribution by fluence modulation in these cases. In addition to the fluence modulation the most suitable single electron energy in each case was also determined. Finally, the simultaneous use of several preselected electron beam energies was also tested, each with an individually optimized fluence profile. One to six electron energies were used, resulting in a slow improvement in complication-free cure with increasing number of beam energies. To apply these techniques to a more clinically relevant situation a post-operative breast cancer patient was studied. For simplicity this patient was treated with only one anterior beam portal to clearly illustrate the effect of inhomogeneities like bone and lung on the dose distribution. It is shown that by using fluence modulation the influence of dose inhomogeneities can be significantly reduced. When two or more electron beam energies with individually optimized fluence profiles are used the dose conformality to the internal target volume is further increased, particularly for targets with complex shapes.
NASA Astrophysics Data System (ADS)
Darazi, R.; Gouze, A.; Macq, B.
2009-01-01
Reproducing a natural and real scene as we see in the real world everyday is becoming more and more popular. Stereoscopic and multi-view techniques are used for this end. However due to the fact that more information are displayed requires supporting technologies such as digital compression to ensure the storage and transmission of the sequences. In this paper, a new scheme for stereo image coding is proposed. The original left and right images are jointly coded. The main idea is to optimally exploit the existing correlation between the two images. This is done by the design of an efficient transform that reduces the existing redundancy in the stereo image pair. This approach was inspired by Lifting Scheme (LS). The novelty in our work is that the prediction step is been replaced by an hybrid step that consists in disparity compensation followed by luminance correction and an optimized prediction step. The proposed scheme can be used for lossless and for lossy coding. Experimental results show improvement in terms of performance and complexity compared to recently proposed methods.
NASA Astrophysics Data System (ADS)
Kohout, B.; Pirinen, J.; Ruiter, N. V.
2012-03-01
The established standard screening method to detect breast cancer is X-ray mammography. However X-ray mammography often has low contrast for tumors located within glandular tissue. A new approach is 3D Ultrasound Computer Tomography (USCT), which is expected to detect small tumors at an early stage. This paper describes the development, improvement and the results of Finite Element Method (FEM) simulations of the Transducer Array System (TAS) used in our 3D USCT. The focus of this work is on researching the influence of meshing and material parameters on the electrical impedance curves. Thereafter, these findings are used to optimize the simulation model. The quality of the simulation was evaluated by comparing simulated impedance characteristics with measured data of the real TAS. The resulting FEM simulation model is a powerful tool to analyze and optimize transducer array systems applied for USCT. With this simulation model, the behavior of TAS for different geometry modifications was researched. It provides a means to understand the acoustical performances inside of any ultrasound transducer represented by its electrical impedance characteristic.
Wang, Yugang; Wu, Xinjun; Sun, Pengfei; Li, Jian
2015-01-01
Electromagnetic acoustic transducers (EMATs) can generate non-dispersive T(0,1) mode guided waves in a metallic pipe for nondestructive testing (NDT) by using a periodic permanent magnet (PPM) EMAT circular array. In order to enhance the excitation efficiency of the sensor, the effects of varying the number of elements of the array on the excitation efficiency is studied in this paper. The transduction process of the PPM EMAT array is studied based on 3-D finite element method (FEM). The passing signal amplitude of the torsional wave is obtained to represent the excitation efficiency of the sensor. Models with different numbers of elements are established and the results are compared to obtain an optimal element number. The simulation result is verified by experiments. It is shown that after optimization, the amplitudes of both the passing signal and defect signal with the optimal element number are increased by 29%, which verifies the feasibility of this optimal method. The essence of the optimization is to find the best match between the static magnetic field and the eddy current field in a limited circumferential space to obtain the maximum circumferential Lorentz force. PMID:25654722
Wang, Yugang; Wu, Xinjun; Sun, Pengfei; Li, Jian
2015-01-01
Electromagnetic acoustic transducers (EMATs) can generate non-dispersive T(0,1) mode guided waves in a metallic pipe for nondestructive testing (NDT) by using a periodic permanent magnet (PPM) EMAT circular array. In order to enhance the excitation efficiency of the sensor, the effects of varying the number of elements of the array on the excitation efficiency is studied in this paper. The transduction process of the PPM EMAT array is studied based on 3-D finite element method (FEM). The passing signal amplitude of the torsional wave is obtained to represent the excitation efficiency of the sensor. Models with different numbers of elements are established and the results are compared to obtain an optimal element number. The simulation result is verified by experiments. It is shown that after optimization, the amplitudes of both the passing signal and defect signal with the optimal element number are increased by 29%, which verifies the feasibility of this optimal method. The essence of the optimization is to find the best match between the static magnetic field and the eddy current field in a limited circumferential space to obtain the maximum circumferential Lorentz force. PMID:25654722
NASA Astrophysics Data System (ADS)
Melchor, J. L., Jr.; Cabrera, S. D.; Aguirre, A.; Kosheleva, O. M.; Vidal, E., Jr.
2005-08-01
This paper describes an efficient algorithm and its Java implementation for a recently developed mean-squared error (MSE) rate-distortion optimal (RDO) inter-slice bit-rate allocation (BRA) scheme applicable to the JPEG2000 Part 2 (J2KP2) framework. Its performance is illustrated on hyperspectral imagery data using the J2KP2 with the Karhunen- Loeve transform (KLT) for decorrelation. The results are contrasted with those obtained using the traditional logvariance based BRA method and with the original RDO algorithm. The implementation has been developed as a Java plug-in to be incorporated into our evolving multi-dimensional data compression software tool denoted CompressMD. The RDO approach to BRA uses discrete rate distortion curves (RDCs) for each slice of transform coefficients. The generation of each point on a RDC requires a full decompression of that slice, therefore, the efficient version minimizes the number of RDC points needed from each slice by using a localized coarse-to-fine approach denoted RDOEfficient. The scheme is illustrated in detail using a subset of 10 bands of hyperspectral imagery data and is contrasted to the original RDO implementation and the traditional (log-variance) method of BRA showing that better results are obtained with the RDO methods. The three schemes are also tested on two hyperspectral imagery data sets with all bands present: the Cuprite radiance data from AVIRIS and a set derived from the Hyperion satellite. The results from the RDO and RDOEfficient are very close to each other in the MSE sense indicating that the adaptive approach can find almost the same BRA solution. Surprisingly, the traditional method also performs very close to the RDO methods, indicating that it is very close to being optimal for these types of data sets.
Parapat, Riny Y; Wijaya, Muliany; Schwarze, Michael; Selve, Sören; Willinger, Marc; Schomäcker, Reinhard
2016-04-01
Correction for 'Particle shape optimization by changing from an isotropic to an anisotropic nanostructure: preparation of highly active and stable supported Pt catalysts in microemulsions' by Riny Y. Parapat et al., Nanoscale, 2013, 5, 796-805. PMID:26961853
NASA Astrophysics Data System (ADS)
Parapat, Riny Y.; Wijaya, Muliany; Schwarze, Michael; Selve, Sören; Willinger, Marc; Schomäcker, Reinhard
2016-03-01
Correction for `Particle shape optimization by changing from an isotropic to an anisotropic nanostructure: preparation of highly active and stable supported Pt catalysts in microemulsions' by Riny Y. Parapat et al., Nanoscale, 2013, 5, 796-805.
NASA Astrophysics Data System (ADS)
Shahbazmohamadi, Sina; Jordan, Eric H.
2012-12-01
Creation of three-dimensional representations of surfaces from images taken at two or more view angles is a well-established technique applied to optical images and is frequently used in combination with scanning electron microscopy (SEM). The present work describes specific steps taken to optimize and enhance the repeatability of three-dimensional surfaces reconstructed from SEM images. The presented steps result in an approximately tenfold improvement in the repeatability of the surface reconstruction compared to more standard techniques. The enhanced techniques presented can be used with any SEM friendly samples. In this work the modified technique was developed in order to accurately quantify surface geometry changes in metallic bond coats used with thermal barrier coatings (TBCs) to provide improved turbine hot part durability. Bond coat surfaces are quite rough, and accurate determination of surface geometry change (rumpling) requires excellent repeatability. Rumpling is an important contributor to TBC failure, and accurate quantification of rumpling is important to better understanding of the failure behavior of TBCs.
[3-D ultrasound in gastroenterology].
Zoller, W G; Liess, H
1994-06-01
Three-dimensional (3D) sonography represents a development of noninvasive diagnostic imaging by real-time two-dimensional (2D) sonography. The use of transparent rotating scans, comparable to a block of glass, generates a 3D effect. The objective of the present study was to optimate 3D presentation of abdominal findings. Additional investigations were made with a new volumetric program to determine the volume of selected findings of the liver. The results were compared with the estimated volumes of 2D sonography and 2D computer tomography (CT). For the processing of 3D images, typical parameter constellations were found for the different findings, which facilitated processing of 3D images. In more than 75% of the cases examined we found an optimal 3D presentation of sonographic findings with respect to the evaluation criteria developed by us for the 3D imaging of processed data. Great differences were found for the estimated volumes of the findings of the liver concerning the three different techniques applied. 3D ultrasound represents a valuable method to judge morphological appearance in abdominal findings. The possibility of volumetric measurements enlarges its potential diagnostic significance. Further clinical investigations are necessary to find out if definite differentiation between benign and malign findings is possible. PMID:7919882
Barry, Robert L; Strother, Stephen C; Gatenby, J Christopher; Gore, John C
2011-04-01
Blood oxygenation level dependent (BOLD) functional magnetic resonance imaging (fMRI) is commonly performed using 2D single-shot echo-planar imaging (EPI). However, single-shot EPI at 7 Tesla (T) often suffers from significant geometric distortions (due to low bandwidth (BW) in the phase-encode (PE) direction) and amplified physiological noise. Recent studies have suggested that 3D multi-shot sequences such as PRESTO may offer comparable BOLD contrast-to-noise ratio with increased volume coverage and decreased geometric distortions. Thus, a four-way group-level comparison was performed between 2D and 3D acquisition sequences at two in-plane resolutions. The quality of fMRI data was evaluated via metrics of prediction and reproducibility using NPAIRS (Non-parametric Prediction, Activation, Influence and Reproducibility re-Sampling). Group activation maps were optimized for each acquisition strategy by selecting the number of principal components that jointly maximized prediction and reproducibility, and showed good agreement in sensitivity and specificity for positive BOLD changes. High-resolution EPI exhibited the highest z-scores of the four acquisition sequences; however, it suffered from the lowest BW in the PE direction (resulting in the worst geometric distortions) and limited spatial coverage, and also caused some subject discomfort through peripheral nerve stimulation (PNS). In comparison, PRESTO also had high z-scores (higher than EPI for a matched in-plane resolution), the highest BW in the PE direction (producing images with superior geometric fidelity), the potential for whole-brain coverage, and no reported PNS. This study provides evidence to support the use of 3D multi-shot acquisition sequences in lieu of single-shot EPI for ultra high field BOLD fMRI at 7T. PMID:21232613
Mishra, Subrata H; Harden, Bradley J
2014-01-01
Structure determination of proteins by solution NMR has become an established method, but challenges increase steeply with the size of proteins. Notably spectral crowding and signal overlap impair the analysis of cross-peaks in NOESY spectra that provide distance restraints for structural models. An optimal spectral resolution can alleviate overlap but requires prohibitively long experimental time with existing methods. Here we present a time-shared 3D experiment optimized for large proteins that provides 15N and 13C dispersed NOESY spectra in a single measurement. NOESY correlations appear in the detected dimension and hence benefit from the highest resolution achievable of all dimensions without increase in experimental time. By design, this experiment is inherently optimal for non-uniform sampling acquisition when compared to current alternatives. Thus, 15N and 13C dispersed NOESY spectra with ultra-high resolution in all dimensions were acquired in parallel within about 4 days instead of 80 days for a 52 kDa monomeric protein at a concentration of 350 μM. PMID:25381567
NASA Astrophysics Data System (ADS)
Marks, R.; Clarke, A.; Featherston, C.; Kawashita, L.; Paget, C.; Pullin, R.
2015-07-01
With the increasing complexity of aircraft structures and materials there is an essential need to continually monitor the structure for damage. This also drives the requirement for optimizing the location of sensors for damage detection to ensure full damage detection coverage of the structure whilst minimizing the number of sensors required, hence reducing costs, weight and data processing. An experiment was carried out to investigate the optimal sensor locations of an active sensor network for detecting adhesive disbonds of a stiffened panel. A piezoelectric transducer was coupled to two different stiffened aluminium panels; one healthy and one with a 25.4mm long disbond. The transducer was positioned at five individual locations to assess the effectiveness of damage detection at different transmission locations. One excitation frequency of 100kHz was used for this study. The panels were scanned with a 3D scanning laser vibrometer which represented a network of ‘ideal’ receiving transducers. The responses measured on the disbonded panel were cross- correlated with those measured on the healthy panel at a large number of potential sensor locations. This generated a cost surface which a genetic algorithm could interrogate in order to find the optimal sensor locations for a given size of sensor network. Probabilistic techniques were used to consider multiple disbond location scenarios, in order to optimise the sensor network for maximum probability of detection across a range of disbond locations.
Parameterization of real-time 3D speckle tracking framework for cardiac strain assessment.
Lorsakul, Auranuch; Duan, Qi; Po, Ming Jack; Angelini, Elsa; Homma, Shunichi; Laine, Andrew F
2011-01-01
Cross-correlation based 3D speckle tracking algorithm can be used to automatically track myocardial motion on three dimensional real-time (RT3D) echocardiography. The goal of this study was to experimentally investigate the effects of different parameters associated with such algorithm to ensure accurate cardiac strain measurements. The investigation was performed on 10 chronic obstructive pulmonary disease RT3DE cardiac ultrasound images. The following two parameters were investigated: 1) the gradient threshold of the anisotropic diffusion pre-filtering and 2) the window size of the cross correlation template matching in the speckle tracking. Results suggest that the optimal gradient threshold of the anisotropic filter depends on the average gradient of the background speckle noise, and that an optimal pair of template size and search window size can be identified determines the cross-correlation level and computational cost. PMID:22254887
Klein, Oliver J; Bhayana, Brijesh; Park, Yong Jin; Evans, Conor L
2012-11-01
Hypoxia and acidosis are widely recognized as major contributors to the development of treatment resistant cancer. For patients with disseminated metastatic lesions, such as most women with ovarian cancer (OvCa), the progression to treatment resistant disease is almost always fatal. Numerous therapeutic approaches have been developed to eliminate treatment resistant carcinoma, including novel biologic, chemo, radiation, and photodynamic therapy (PDT) regimens. Recently, PDT using the cationic photosensitizer EtNBS was found to be highly effective against therapeutically unresponsive hypoxic and acidic OvCa cellular populations in vitro. To optimize this treatment regimen, we developed a tiered, high-content, image-based screening approach utilizing a biologically relevant OvCa 3D culture model to investigate a small library of side-chain modified EtNBS derivatives. The uptake, localization, and photocytotoxicity of these compounds on both the cellular and nodular levels were observed to be largely mediated by their respective ethyl side chain chemical alterations. In particular, EtNBS and its hydroxyl-terminated derivative (EtNBS-OH) were found to have similar pharmacological parameters, such as their nodular localization patterns and uptake kinetics. Interestingly, these two molecules were found to induce dramatically different therapeutic outcomes: EtNBS was found to be more effective in killing the hypoxic, nodule core cells with superior selectivity, while EtNBS-OH was observed to trigger widespread structural degradation of nodules. This breakdown of the tumor architecture can improve the therapeutic outcome and is known to synergistically enhance the antitumor effects of front-line chemotherapeutic regimens. These results, which would not have been predicted or observed using traditional monolayer or in vivo animal screening techniques, demonstrate the powerful capabilities of 3D in vitro screening approaches for the selection and optimization of therapeutic
NASA Astrophysics Data System (ADS)
Pletinckx, D.
2011-09-01
The current 3D hype creates a lot of interest in 3D. People go to 3D movies, but are we ready to use 3D in our homes, in our offices, in our communication? Are we ready to deliver real 3D to a general public and use interactive 3D in a meaningful way to enjoy, learn, communicate? The CARARE project is realising this for the moment in the domain of monuments and archaeology, so that real 3D of archaeological sites and European monuments will be available to the general public by 2012. There are several aspects to this endeavour. First of all is the technical aspect of flawlessly delivering 3D content over all platforms and operating systems, without installing software. We have currently a working solution in PDF, but HTML5 will probably be the future. Secondly, there is still little knowledge on how to create 3D learning objects, 3D tourist information or 3D scholarly communication. We are still in a prototype phase when it comes to integrate 3D objects in physical or virtual museums. Nevertheless, Europeana has a tremendous potential as a multi-facetted virtual museum. Finally, 3D has a large potential to act as a hub of information, linking to related 2D imagery, texts, video, sound. We describe how to create such rich, explorable 3D objects that can be used intuitively by the generic Europeana user and what metadata is needed to support the semantic linking.
NIF Ignition Target 3D Point Design
Jones, O; Marinak, M; Milovich, J; Callahan, D
2008-11-05
We have developed an input file for running 3D NIF hohlraums that is optimized such that it can be run in 1-2 days on parallel computers. We have incorporated increasing levels of automation into the 3D input file: (1) Configuration controlled input files; (2) Common file for 2D and 3D, different types of capsules (symcap, etc.); and (3) Can obtain target dimensions, laser pulse, and diagnostics settings automatically from NIF Campaign Management Tool. Using 3D Hydra calculations to investigate different problems: (1) Intrinsic 3D asymmetry; (2) Tolerance to nonideal 3D effects (e.g. laser power balance, pointing errors); and (3) Synthetic diagnostics.
3d-3d correspondence revisited
NASA Astrophysics Data System (ADS)
Chung, Hee-Joong; Dimofte, Tudor; Gukov, Sergei; Sułkowski, Piotr
2016-04-01
In fivebrane compactifications on 3-manifolds, we point out the importance of all flat connections in the proper definition of the effective 3d {N}=2 theory. The Lagrangians of some theories with the desired properties can be constructed with the help of homological knot invariants that categorify colored Jones polynomials. Higgsing the full 3d theories constructed this way recovers theories found previously by Dimofte-Gaiotto-Gukov. We also consider the cutting and gluing of 3-manifolds along smooth boundaries and the role played by all flat connections in this operation.
Generic remeshing of 3D triangular meshes with metric-dependent discrete voronoi diagrams.
Valette, Sebastien; Chassery, Jean Marc; Prost, Rémy
2008-01-01
In this paper, we propose a generic framework for 3D surface remeshing. Based on a metric-driven Discrete Voronoi Diagram construction, our output is an optimized 3D triangular mesh with a user defined vertex budget. Our approach can deal with a wide range of applications, from high quality mesh generation to shape approximation. By using appropriate metric constraints the method generates isotropic or anisotropic elements. Based on point-sampling, our algorithm combines the robustness and theoretical strength of Delaunay criteria with the efficiency of entirely discrete geometry processing . Besides the general described framework, we show experimental results using isotropic, quadric-enhanced isotropic and anisotropic metrics which prove the efficiency of our method on large meshes, for a low computational cost. PMID:18192716
NASA Astrophysics Data System (ADS)
Ko, Hsu-Feng
Sustained delivery of therapeutic genes from a three-dimensional (3-D) scaffold and subsequent gene expression capable of triggering the regeneration of damaged tissues is a tissue engineering strategy that has been gaining increased attention. Nanostructured calcium phosphates (NanoCaPs) are biocompatible and non-toxic biomaterials. Furthermore, their efficient transfection in vitro have rendered them attractive gene delivery carriers compared to other viral- or lipid-based agents that tend to be immunogenic or cytotoxic, leading to undesirable responses when utilized above a critical threshold. However, NanoCaPs are typically characterized by variable transfection and short shelf life due to particle aggregation. A viable solution to this problem is the incorporation of NanoCaPs into 3-D scaffolds. The main objectives of this research are therefore two-fold: (1) Examination of the potential of achieving optimized transfection of NanoCaPs via anionic substitution and (2) high throughput synthesis and screening of non-viral gene delivery systems (GDS) comprised of naturally-derived polymers as scaffolds containing NanoCaPs gene carriers. Results indicated that in addition to the excellent transfection levels exhibited by NanoCaPs in vitro, an additional 20-30% increase was observed for NanoCaPs with 10-25 mol% anion substitution. In contrast, high anion substitution (>60%) yielded a drastic decline in transfection. Structural characterizations verified successful anion substitution with a noticeable increase in lattice parameters indicative of an expanded unit cell due to ionic substitution. All of the anion substituted calcium phosphates exhibited the primary phase of hydroxyapatite. For the first time, GDS composed of various concentrations of alginate (AA), fibronectin (FN), and NanoCaPs-DNA complexes were demonstrated. The presence of AA and FN was effective in immobilizing NanoCaPs and reducing the aggregation. High throughput synthesis and screening
Freed, Melanie; Park, Subok; Badano, Aldo
2010-01-01
Purpose: Accurate models of detector blur are crucial for performing meaningful optimizations of three-dimensional (3D) x-ray breast imaging systems as well as for developing reconstruction algorithms that faithfully reproduce the imaged object anatomy. So far, x-ray detector blur has either been ignored or modeled as a shift-invariant symmetric function for these applications. The recent development of a Monte Carlo simulation package called MANTIS has allowed detailed modeling of these detector blur functions and demonstrated the magnitude of the anisotropy for both tomosynthesis and breast CT imaging systems. Despite the detailed results that MANTIS produces, the long simulation times required make inclusion of these results impractical in rigorous optimization and reconstruction algorithms. As a result, there is a need for detector blur models that can be rapidly generated. Methods: In this study, the authors have derived an analytical model for deterministic detector blur functions, referred to here as point response functions (PRFs), of columnar CsI phosphor screens. The analytical model is x-ray energy and incidence angle dependent and draws on results from MANTIS to indirectly include complicated interactions that are not explicitly included in the mathematical model. Once the mathematical expression is derived, values of the coefficients are determined by a two-dimensional (2D) fit to MANTIS-generated results based on a figure-of-merit (FOM) that measures the normalized differences between the MANTIS and analytical model results averaged over a region of interest. A smaller FOM indicates a better fit. This analysis was performed for a monochromatic x-ray energy of 25 keV, a CsI scintillator thickness of 150 μm, and four incidence angles (0°, 15°, 30°, and 45°). Results: The FOMs comparing the analytical model to MANTIS for these parameters were 0.1951±0.0011, 0.1915±0.0014, 0.2266±0.0021, and 0.2416±0.0074 for 0°, 15°, 30°, and 45
NASA Astrophysics Data System (ADS)
Kobayashi, M.; Xu, Y.; Ida, K.; Corre, Y.; Feng, Y.; Schmitz, O.; Frerichs, H.; Tabares, F. L.; Evans, T. E.; Coenen, J. W.; Liang, Y.; Bader, A.; Itoh, K.; Yamada, H.; Ghendrih, Ph.; Ciraolo, G.; Tafalla, D.; Lopez-Fraguas, A.; Guo, H. Y.; Cui, Z. Y.; Reiter, D.; Asakura, N.; Wenzel, U.; Morita, S.; Ohno, N.; Peterson, B. J.; Masuzaki, S.
2015-10-01
This paper assesses the three-dimensional (3D) effects of the edge magnetic field structure on divertor/scrape-off layer transport, based on an inter-machine comparison of experimental data and on the recent progress of 3D edge transport simulation. The 3D effects are elucidated as a consequence of competition between transports parallel (\\parallel ) and perpendicular (\\bot ) to the magnetic field, in open field lines cut by divertor plates, or in magnetic islands. The competition has strong impacts on divertor functions, such as determination of the divertor density regime, impurity screening and detachment control. The effects of magnetic perturbation on the edge electric field and turbulent transport are also discussed. Parameterization to measure the 3D effects on the edge transport is attempted for the individual divertor functions. Based on the suggested key parameters, an operation domain of the 3D divertor configuration is discussed for future devices.
Su, M; Sura, S
2014-06-01
Purpose: To evaluate dosimetric impact of two axillary nodes (AX) boost techniques: (1) posterior-oblique optimized field boost (POB), (2) traditional posterior-anterior boost (PAB) with field optimization (O-PAB), for a postmastectomy breast patient with positive axillary lymph nodes. Methods: Five patients, 3 left and 2 right chest walls, were included in this study. All patients were simulated in 5mm CT slice thickness. Supraclavicular (SC) and level I/II/III AX were contoured based on the RTOG atlas guideline. Five treatment plans, (1) tangential chest wall, (2) oblique SC including AX, (3) PAB, O-PAB and POB, were created for each patient. Three plan sums (PS) were generated by sum one of (3) plan with plan (1) and (2). The field optimization was done through PS dose distribution, which included a field adjustment, a fractional dose, a calculation location and a gantry angle selection for POB. A dosimetric impact was evaluated by comparing a SC and AX coverage, a PS maximum dose, an irradiated area percentage volume received dose over 105% prescription dose (V105), an ipsi-laterial mean lung dose (MLD), an ipsi-laterial mean humeral head dose (MHHD), a mean heart dose (MHD) (for left case only) and their DVH amount these three technique. Results: O-PAB, POB and PAB dosimetric results showed that there was no significant different on SC and AX coverage (p>0.43) and MHD (p>0.16). The benefit of sparing lung irradiation from PAB to O-PAB to POB was significant (p<0.004). PAB showed a highest PS maximum dose (p<0.005), V105 (p<0.023) and MLD (compared with OPAB, p=0.055). MHHD showed very sensitive to the patient arm positioning and anatomy. O-PAB convinced a lower MHHD than PAB (p=0.03). Conclusion: 3D CT contouring plays main role in accuracy radiotherapy. Dosimetric advantage of POB and O-PAB was observed for a better normal tissue irradiation sparing.
NASA Astrophysics Data System (ADS)
Meulien Ohlmann, Odile
2013-02-01
Today the industry offers a chain of 3D products. Learning to "read" and to "create in 3D" becomes an issue of education of primary importance. 25 years professional experience in France, the United States and Germany, Odile Meulien set up a personal method of initiation to 3D creation that entails the spatial/temporal experience of the holographic visual. She will present some different tools and techniques used for this learning, their advantages and disadvantages, programs and issues of educational policies, constraints and expectations related to the development of new techniques for 3D imaging. Although the creation of display holograms is very much reduced compared to the creation of the 90ies, the holographic concept is spreading in all scientific, social, and artistic activities of our present time. She will also raise many questions: What means 3D? Is it communication? Is it perception? How the seeing and none seeing is interferes? What else has to be taken in consideration to communicate in 3D? How to handle the non visible relations of moving objects with subjects? Does this transform our model of exchange with others? What kind of interaction this has with our everyday life? Then come more practical questions: How to learn creating 3D visualization, to learn 3D grammar, 3D language, 3D thinking? What for? At what level? In which matter? for whom?
NASA Technical Reports Server (NTRS)
Broussard, J. R.; Halyo, N.
1984-01-01
This report contains the development of a digital outer-loop three dimensional radio navigation (3-D RNAV) flight control system for a small commercial jet transport. The outer-loop control system is designed using optimal stochastic limited state feedback techniques. Options investigated using the optimal limited state feedback approach include integrated versus hierarchical control loop designs, 20 samples per second versus 5 samples per second outer-loop operation and alternative Type 1 integration command errors. Command generator tracking techniques used in the digital control design enable the jet transport to automatically track arbitrary curved flight paths generated by waypoints. The performance of the design is demonstrated using detailed nonlinear aircraft simulations in the terminal area, frequency domain multi-input sigma plots, frequency domain single-input Bode plots and closed-loop poles. The response of the system to a severe wind shear during a landing approach is also presented.
Rilling, M; Goulet, M; Thibault, S; Archambault, L
2015-06-15
specifications. This work leads the way to improving the 3D dosimeter’s achievable resolution, efficiency and build for providing a quality assurance tool fully meeting clinical needs. M.R. is financially supported by a Master’s Canada Graduate Scholarship from the NSERC. This research is also supported by the NSERC Industrial Research Chair in Optical Design.
ERIC Educational Resources Information Center
Hastings, S. K.
2002-01-01
Discusses 3 D imaging as it relates to digital representations in virtual library collections. Highlights include X-ray computed tomography (X-ray CT); the National Science Foundation (NSF) Digital Library Initiatives; output peripherals; image retrieval systems, including metadata; and applications of 3 D imaging for libraries and museums. (LRW)
Stoian, Sebastian A; Paraschiv, Carmen; Kiritsakas, Nathalie; Lloret, Francesc; Münck, Eckard; Bominaar, Emile L; Andruh, Marius
2010-04-01
The synthesis and crystallographic characterization of a new family of M(mu-CN)Ln complexes are reported. Two structural series have been prepared by reacting in water rare earth nitrates (Ln(III) = La, Pr, Nd, Sm, Eu, Gd, Dy, Ho) with K(3)[M(CN)(6)] (M(III) = Fe, Co) in the presence of hexamethylenetetramine (hmt). The first series consists of six isomorphous heterobinuclear complexes, [(CN)(5)M-CN-Ln(H(2)O)(8)].2hmt ([FeLa] 1, [FePr] 2, [FeNd] 3, [FeSm] 4, [FeEu] 5, [FeGd] 6), while the second series consists of four isostructural ionic complexes, [M(CN)(6)][Ln(H(2)O)(8)].hmt ([FeDy] 7, [FeHo] 8, [CoEu] 9, [CoGd] 10). The hexamethylenetetramine molecules contribute to the stabilization of the crystals by participating in an extended network of hydrogen bond interactions. In both series the aqua ligands are hydrogen bonded to the nitrogen atoms from both the terminal CN(-) groups and the hmt molecules. The [FeGd] complex has been analyzed with (57)Fe Mossbauer spectroscopy and magnetic susceptibility measurements. We have also analyzed the [FeLa] complex, in which the paramagnetic Gd(III) is replaced by diamagnetic La(III), with (57)Fe Mossbauer spectroscopy, electron paramagnetic resonance (EPR), and magnetic susceptibility measurements, to obtain information about the low-spin Fe(III) site that is not accessible in the presence of a paramagnetic ion at the complementary site. For the same reason, the [CoGd] complex, containing diamagnetic Co(III), was studied with EPR and magnetic susceptibility measurements, which confirmed the S = 7/2 spin of Gd(III). Prior knowledge about the paramagnetic sites in [FeGd] allows a detailed analysis of the exchange interactions between them. In particular, the question of whether the exchange interaction in [FeGd] is isotropic or anisotropic has been addressed. Standard variable-temperature magnetic susceptibility measurements provide only the value for a linear combination of J(x), J(y), and J(z) but contain no information
Parallel Anisotropic Tetrahedral Adaptation
NASA Technical Reports Server (NTRS)
Park, Michael A.; Darmofal, David L.
2008-01-01
An adaptive method that robustly produces high aspect ratio tetrahedra to a general 3D metric specification without introducing hybrid semi-structured regions is presented. The elemental operators and higher-level logic is described with their respective domain-decomposed parallelizations. An anisotropic tetrahedral grid adaptation scheme is demonstrated for 1000-1 stretching for a simple cube geometry. This form of adaptation is applicable to more complex domain boundaries via a cut-cell approach as demonstrated by a parallel 3D supersonic simulation of a complex fighter aircraft. To avoid the assumptions and approximations required to form a metric to specify adaptation, an approach is introduced that directly evaluates interpolation error. The grid is adapted to reduce and equidistribute this interpolation error calculation without the use of an intervening anisotropic metric. Direct interpolation error adaptation is illustrated for 1D and 3D domains.
Crandall, K.R.
1987-08-01
TRACE 3-D is an interactive beam-dynamics program that calculates the envelopes of a bunched beam, including linear space-charge forces, through a user-defined transport system. TRACE 3-D provides an immediate graphics display of the envelopes and the phase-space ellipses and allows nine types of beam-matching options. This report describes the beam-dynamics calculations and gives detailed instruction for using the code. Several examples are described in detail.
Shri Ramaswamy, University of Minnesota; B.V. Ramarao, State University of New York
2004-08-29
the top and bottom sides of the porous material, i.e. "two-sidedness" due to processing and raw material characteristics may lead to differences in end-use performance. The measurements of surface structure characteristics include thickness distribution, surface volume distribution, contact fraction distribution and surface pit distribution. This complements our earlier method to analyze the bulk structure and Z-D structure of porous materials. As one would expect, the surface structure characteristics will be critically dependent on the quality and resolution of the images. This presents a useful tool to characterize and engineer the surface structure of porous materials such as paper and board tailored to specific end-use applications. This will also help troubleshoot problems related to manufacturing and end-use applications. This study attempted to identify the optimal resolution through a comparison between 3D images obtained by monochromatic synchrotron radiation X-CT in phase contrast mode (resolution 1 m) and polychromatic radiation X-CT in absorption mode (res.5 m). It was found that both resolutions have the ability to show the expected trends when comparing different paper samples. The low resolution technique shows fewer details resulting in lower specific surface area, larger pore channels, characterized as hydraulic radii, and lower tortuosities, where differences between samples and principal directions are more difficult to detect. The disadvantages of the high resolution images are high cost and limited availability of hard x-ray beam time as well as the small size of the sample volumes imaged. The results show that the low resolution images can be used for comparative studies, whereas the high resolution images may be better suited for fundamental research on the paper structure and its influence on paper properties, as one gets more accurate physical measurements. In addition, pore space diffusion model has been developed to simulate simultaneous
NASA Astrophysics Data System (ADS)
Oldham, Mark
2015-01-01
Radiochromic materials exhibit a colour change when exposed to ionising radiation. Radiochromic film has been used for clinical dosimetry for many years and increasingly so recently, as films of higher sensitivities have become available. The two principle advantages of radiochromic dosimetry include greater tissue equivalence (radiologically) and the lack of requirement for development of the colour change. In a radiochromic material, the colour change arises direct from ionising interactions affecting dye molecules, without requiring any latent chemical, optical or thermal development, with important implications for increased accuracy and convenience. It is only relatively recently however, that 3D radiochromic dosimetry has become possible. In this article we review recent developments and the current state-of-the-art of 3D radiochromic dosimetry, and the potential for a more comprehensive solution for the verification of complex radiation therapy treatments, and 3D dose measurement in general.
NASA Astrophysics Data System (ADS)
Chixiao, Chen; Jixuan, Xiang; Huabin, Chen; Jun, Xu; Fan, Ye; Ning, Li; Junyan, Ren
2015-05-01
Asynchronous successive approximation register (SAR) analog-to-digital converters (ADC) feature high energy efficiency but medium performance. From the point of view of speed, the key bottleneck is the unit capacitor size. In this paper, a small size three-dimensional (3-D) metal—oxide—metal (MOM) capacitor is proposed. The unit capacitor has a capacitance of 1-fF. It shapes as an umbrella, which is designed for fast settling consideration. A comparison among the proposed capacitor with other 3-D MOM capacitors is also given in the paper. To demonstrate the effectiveness of the MOM capacitor, a 6-b capacitive DAC is implemented in TSMC 1P9M 65 nm LP CMOS technology. The DAC consumes a power dissipation of 0.16 mW at the rate of 100 MS/s, excluding a source-follower based output buffer. Static measurement result shows that INL is less than ±1 LSB and DNL is less than ±0.5 LSB. In addition, a 100 MS/s 9-bit SAR ADC with the proposed 3-D capacitor is simulated.
Wood, Michael R; Noetzel, Meredith J; Engers, Julie L; Bollinger, Katrina A; Melancon, Bruce J; Tarr, James C; Han, Changho; West, Mary; Gregro, Alison R; Lamsal, Atin; Chang, Sichen; Ajmera, Sonia; Smith, Emery; Chase, Peter; Hodder, Peter S; Bubser, Michael; Jones, Carrie K; Hopkins, Corey R; Emmitte, Kyle A; Niswender, Colleen M; Wood, Michael W; Duggan, Mark E; Conn, P Jeffrey; Bridges, Thomas M; Lindsley, Craig W
2016-07-01
This Letter describes the chemical optimization of a novel series of M4 positive allosteric modulators (PAMs) based on a 5,6-dimethyl-4-(piperidin-1-yl)thieno[2,3-d]pyrimidine core, identified from an MLPCN functional high-throughput screen. The HTS hit was potent and selective, but not CNS penetrant. Potency was maintained, while CNS penetration was improved (rat brain:plasma Kp=0.74), within the original core after several rounds of optimization; however, the thieno[2,3-d]pyrimidine core was subject to extensive oxidative metabolism. Ultimately, we identified a 6-fluoroquinazoline core replacement that afforded good M4 PAM potency, muscarinic receptor subtype selectivity and CNS penetration (rat brain:plasma Kp>10). Moreover, this campaign provided fundamentally distinct M4 PAM chemotypes, greatly expanding the available structural diversity for this exciting CNS target. PMID:27185330
NASA Astrophysics Data System (ADS)
Iliesiu, Luca; Kos, Filip; Poland, David; Pufu, Silviu S.; Simmons-Duffin, David; Yacoby, Ran
2016-03-01
We study the conformal bootstrap for a 4-point function of fermions < ψψψψ> in 3D. We first introduce an embedding formalism for 3D spinors and compute the conformal blocks appearing in fermion 4-point functions. Using these results, we find general bounds on the dimensions of operators appearing in the ψ × ψ OPE, and also on the central charge C T . We observe features in our bounds that coincide with scaling dimensions in the GrossNeveu models at large N . We also speculate that other features could coincide with a fermionic CFT containing no relevant scalar operators.
3-D Technology Approaches for Biological Ecologies
NASA Astrophysics Data System (ADS)
Liu, Liyu; Austin, Robert; U. S-China Physical-Oncology Sciences Alliance (PS-OA) Team
Constructing three dimensional (3-D) landscapes is an inevitable issue in deep study of biological ecologies, because in whatever scales in nature, all of the ecosystems are composed by complex 3-D environments and biological behaviors. Just imagine if a 3-D technology could help complex ecosystems be built easily and mimic in vivo microenvironment realistically with flexible environmental controls, it will be a fantastic and powerful thrust to assist researchers for explorations. For years, we have been utilizing and developing different technologies for constructing 3-D micro landscapes for biophysics studies in in vitro. Here, I will review our past efforts, including probing cancer cell invasiveness with 3-D silicon based Tepuis, constructing 3-D microenvironment for cell invasion and metastasis through polydimethylsiloxane (PDMS) soft lithography, as well as explorations of optimized stenting positions for coronary bifurcation disease with 3-D wax printing and the latest home designed 3-D bio-printer. Although 3-D technologies is currently considered not mature enough for arbitrary 3-D micro-ecological models with easy design and fabrication, I hope through my talk, the audiences will be able to sense its significance and predictable breakthroughs in the near future. This work was supported by the State Key Development Program for Basic Research of China (Grant No. 2013CB837200), the National Natural Science Foundation of China (Grant No. 11474345) and the Beijing Natural Science Foundation (Grant No. 7154221).
NASA Astrophysics Data System (ADS)
Ayala, Conxi; Izquierdo-Llavall, Esther; Pueyo, Emilio Luis; Rubio, Félix; Rodríguez-Pintó, Adriana; María Casas, Antonio; Oliva-Urcía, Belén; Rey-Moral, Carmen
2015-04-01
Obtaining an accurate 3D image of the geometry and physical properties of geological structures in depth is a challenge regardless the scale and the aim of the investigation. In this framework, assessing the origin of the uncertainties and reducing them is a key issue when building a 3D reconstruction of a target area. Usually, this process involves an interdisciplinary approach and also the use of different software whose inputs and outputs have to be interoperable. We have designed a new workflow for 2.5D and 3D geological and potential field modelling, especially useful in areas where no seismic data is available. The final aim is to obtain a 3D geological model, at a regional or local scale, with the smaller uncertainty as possible. Once the study area and the working scale are is decided, the first obvious step is to compile all preexisting data and to determine its uncertainties. If necessary, a survey will be carried out to acquire additional data (e.g., gravity, magnetic or petrophysical data) to have an appropriated coverage of information and rock samples. A thorough study of the petrophysical properties is made to determine the density, magnetic susceptibility and remanence that will be assigned to each lithology, together with its corresponding uncertainty. Finally, the modelling process is started, and it includes a feedback between geology and potential fields in order to progressively refine the model until it fits all the existing data. The procedure starts with the construction of balanced geological cross sections from field work, available geological maps as well as data from stratigraphic columns, boreholes, etc. These geological cross sections are exported and imported in GMSYS software to carry out the 2.5D potential field modelling. The model improves and its uncertainty is reduced through the feedback between the geologists and the geophysicists. Once the potential field anomalies are well adjusted, the cross sections are exported into 3
Measuring the orientation and rotation rate of 3D printed particles in turbulent flow
NASA Astrophysics Data System (ADS)
Voth, Greg; Marcus, Guy G.; Parsa, Shima; Kramel, Stefan; Ni, Rui; Cole, Brendan
2014-11-01
The orientation distribution and rotations of anisotropic particles plays a key role in many applications ranging from icy clouds to papermaking and drag reduction in pipe flow. Experimental access to time resolved orientations of anisotropic particles has not been easy to achieve. We have found that 3D printing technology can be used to fabricate a wide range of particle shapes with smallest dimension down to 300 μm. So far we have studied rods, crosses, jacks, tetrads, and helical shapes. We extract the particle orientations from stereoscopic video images using a method of least squares optimization in Euler angle space. We find that in turbulence the orientation and rotation rate of many particles can be understood using a simple picture of alignment of both the vorticity and a long axis of the particle with the Lagrangian stretching direction of the flow. This research is supported by NSF Grant DMR-1208990.
Measuring the orientation and rotation rate of 3D printed particles in turbulent flow
NASA Astrophysics Data System (ADS)
Voth, Greg; Kramel, Stefan; Cole, Brendan
2015-03-01
The orientation distribution and rotations of anisotropic particles plays a key role in many applications ranging from icy clouds to papermaking and drag reduction in pipe flow. Experimental access to time resolved orientations of anisotropic particles has not been easy to achieve. We have found that 3D printing technology can be used to fabricate a wide range of particle shapes with smallest dimension down to 300 ?m. So far we have studied rods, crosses, jacks, tetrads, and helical shapes. We extract the particle orientations from stereoscopic video images using a method of least squares optimization in Euler angle space. We find that in turbulence the orientation and rotation rate of many particles can be understood using a simple picture of alignment of both the vorticity and a long axis of the particle with the Lagrangian stretching direction of the flow.
NASA Astrophysics Data System (ADS)
Iizuka, Keigo
2008-02-01
In order to circumvent the fact that only one observer can view the image from a stereoscopic microscope, an attachment was devised for displaying the 3D microscopic image on a large LCD monitor for viewing by multiple observers in real time. The principle of operation, design, fabrication, and performance are presented, along with tolerance measurements relating to the properties of the cellophane half-wave plate used in the design.
Anisotropy effects on 3D waveform inversion
NASA Astrophysics Data System (ADS)
Stekl, I.; Warner, M.; Umpleby, A.
2010-12-01
In the recent years 3D waveform inversion has become achievable procedure for seismic data processing. A number of datasets has been inverted and presented (Warner el al 2008, Ben Hadj at all, Sirgue et all 2010) using isotropic 3D waveform inversion. However the question arises will the results be affected by isotropic assumption. Full-wavefield inversion techniques seek to match field data, wiggle-for-wiggle, to synthetic data generated by a high-resolution model of the sub-surface. In this endeavour, correctly matching the travel times of the principal arrivals is a necessary minimal requirement. In many, perhaps most, long-offset and wide-azimuth datasets, it is necessary to introduce some form of p-wave velocity anisotropy to match the travel times successfully. If this anisotropy is not also incorporated into the wavefield inversion, then results from the inversion will necessarily be compromised. We have incorporated anisotropy into our 3D wavefield tomography codes, characterised as spatially varying transverse isotropy with a tilted axis of symmetry - TTI anisotropy. This enhancement approximately doubles both the run time and the memory requirements of the code. We show that neglect of anisotropy can lead to significant artefacts in the recovered velocity models. We will present inversion results of inverting anisotropic 3D dataset by assuming isotropic earth and compare them with anisotropic inversion result. As a test case Marmousi model extended to 3D with no velocity variation in third direction and with added spatially varying anisotropy is used. Acquisition geometry is assumed as OBC with sources and receivers everywhere at the surface. We attempted inversion using both 2D and full 3D acquisition for this dataset. Results show that if no anisotropy is taken into account although image looks plausible most features are miss positioned in depth and space, even for relatively low anisotropy, which leads to incorrect result. This may lead to
Estimation of the thermal conductivity of hemp based insulation material from 3D tomographic images
NASA Astrophysics Data System (ADS)
El-Sawalhi, R.; Lux, J.; Salagnac, P.
2016-08-01
In this work, we are interested in the structural and thermal characterization of natural fiber insulation materials. The thermal performance of these materials depends on the arrangement of fibers, which is the consequence of the manufacturing process. In order to optimize these materials, thermal conductivity models can be used to correlate some relevant structural parameters with the effective thermal conductivity. However, only a few models are able to take into account the anisotropy of such material related to the fibers orientation, and these models still need realistic input data (fiber orientation distribution, porosity, etc.). The structural characteristics are here directly measured on a 3D tomographic image using advanced image analysis techniques. Critical structural parameters like porosity, pore and fiber size distribution as well as local fiber orientation distribution are measured. The results of the tested conductivity models are then compared with the conductivity tensor obtained by numerical simulation on the discretized 3D microstructure, as well as available experimental measurements. We show that 1D analytical models are generally not suitable for assessing the thermal conductivity of such anisotropic media. Yet, a few anisotropic models can still be of interest to relate some structural parameters, like the fiber orientation distribution, to the thermal properties. Finally, our results emphasize that numerical simulations on 3D realistic microstructure is a very interesting alternative to experimental measurements.
Lossy compression of hyperspectral images using shearlet transform and 3D SPECK
NASA Astrophysics Data System (ADS)
Karami, A.
2015-10-01
In this paper, a new lossy compression method for hyperspectral images (HSI) is introduced. HSI are considered as a 3D dataset with two dimensions in the spatial and one dimension in the spectral domain. In the proposed method, first 3D multidirectional anisotropic shearlet transform is applied to the HSI. Because, unlike traditional wavelets, shearlets are theoretically optimal in representing images with edges and other geometrical features. Second, soft thresholding method is applied to the shearlet transform coefficients and finally the modified coefficients are encoded using Three Dimensional- Set Partitioned Embedded bloCK (3D SPECK). Our simulation results show that the proposed method, in comparison with well-known approaches such as 3D SPECK (using 3D wavelet) and combined PCA and JPEG2000 algorithms, provides a higher SNR (signal to noise ratio) for any given compression ratio (CR). It is noteworthy to mention that the superiority of proposed method is distinguishable as the value of CR grows. In addition, the effect of proposed method on the spectral unmixing analysis is also evaluated.