3-D UNSTRUCTURED HEXAHEDRAL-MESH Sn TRANSPORT METHODS
J. MOREL; J. MCGHEE; ET AL
2000-11-01
This is the final report of a three-year, Laboratory-Directed Research and Development (LDRD) project at the Los Alamos National Laboratory (LANL). We have developed a method for solving the neutral-particle transport equation on 3-D unstructured hexahedral meshes using a S{sub n} discretization in angle in conjunction with a discontinuous finite-element discretization in space and a multigroup discretization in energy. Previous methods for solving this equation in 3-D have been limited to rectangular meshes. The unstructured-mesh method that we have developed is far more efficient for solving problems with complex 3-D geometric features than rectangular-mesh methods. In spite of having to make several compromises in our spatial discretization technique and our iterative solution technique, our method has been found to be both accurate and efficient for a broad class of problems.
GPU accelerated simulations of 3D deterministic particle transport using discrete ordinates method
Gong Chunye; Liu Jie; Chi Lihua; Huang Haowei; Fang Jingyue; Gong Zhenghu
2011-07-01
Graphics Processing Unit (GPU), originally developed for real-time, high-definition 3D graphics in computer games, now provides great faculty in solving scientific applications. The basis of particle transport simulation is the time-dependent, multi-group, inhomogeneous Boltzmann transport equation. The numerical solution to the Boltzmann equation involves the discrete ordinates (S{sub n}) method and the procedure of source iteration. In this paper, we present a GPU accelerated simulation of one energy group time-independent deterministic discrete ordinates particle transport in 3D Cartesian geometry (Sweep3D). The performance of the GPU simulations are reported with the simulations of vacuum boundary condition. The discussion of the relative advantages and disadvantages of the GPU implementation, the simulation on multi GPUs, the programming effort and code portability are also reported. The results show that the overall performance speedup of one NVIDIA Tesla M2050 GPU ranges from 2.56 compared with one Intel Xeon X5670 chip to 8.14 compared with one Intel Core Q6600 chip for no flux fixup. The simulation with flux fixup on one M2050 is 1.23 times faster than on one X5670.
GPU accelerated simulations of 3D deterministic particle transport using discrete ordinates method
NASA Astrophysics Data System (ADS)
Gong, Chunye; Liu, Jie; Chi, Lihua; Huang, Haowei; Fang, Jingyue; Gong, Zhenghu
2011-07-01
Graphics Processing Unit (GPU), originally developed for real-time, high-definition 3D graphics in computer games, now provides great faculty in solving scientific applications. The basis of particle transport simulation is the time-dependent, multi-group, inhomogeneous Boltzmann transport equation. The numerical solution to the Boltzmann equation involves the discrete ordinates ( Sn) method and the procedure of source iteration. In this paper, we present a GPU accelerated simulation of one energy group time-independent deterministic discrete ordinates particle transport in 3D Cartesian geometry (Sweep3D). The performance of the GPU simulations are reported with the simulations of vacuum boundary condition. The discussion of the relative advantages and disadvantages of the GPU implementation, the simulation on multi GPUs, the programming effort and code portability are also reported. The results show that the overall performance speedup of one NVIDIA Tesla M2050 GPU ranges from 2.56 compared with one Intel Xeon X5670 chip to 8.14 compared with one Intel Core Q6600 chip for no flux fixup. The simulation with flux fixup on one M2050 is 1.23 times faster than on one X5670.
An unstaggered constrained transport method for the 3D ideal magnetohydrodynamic equations
NASA Astrophysics Data System (ADS)
Helzel, Christiane; Rossmanith, James A.; Taetz, Bertram
2011-05-01
Numerical methods for solving the ideal magnetohydrodynamic (MHD) equations in more than one space dimension must either confront the challenge of controlling errors in the discrete divergence of the magnetic field, or else be faced with nonlinear numerical instabilities. One approach for controlling the discrete divergence is through a so-called constrained transport method, which is based on first predicting a magnetic field through a standard finite volume solver, and then correcting this field through the appropriate use of a magnetic vector potential. In this work we develop a constrained transport method for the 3D ideal MHD equations that is based on a high-resolution wave propagation scheme. Our proposed scheme is the 3D extension of the 2D scheme developed by Rossmanith [J.A. Rossmanith, An unstaggered, high-resolution constrained transport method for magnetohydrodynamic flows, SIAM J. Sci. Comput. 28 (2006) 1766], and is based on the high-resolution wave propagation method of Langseth and LeVeque [J.O. Langseth, R.J. LeVeque, A wave propagation method for threedimensional hyperbolic conservation laws, J. Comput. Phys. 165 (2000) 126]. In particular, in our extension we take great care to maintain the three most important properties of the 2D scheme: (1) all quantities, including all components of the magnetic field and magnetic potential, are treated as cell-centered; (2) we develop a high-resolution wave propagation scheme for evolving the magnetic potential; and (3) we develop a wave limiting approach that is applied during the vector potential evolution, which controls unphysical oscillations in the magnetic field. One of the key numerical difficulties that is novel to 3D is that the transport equation that must be solved for the magnetic vector potential is only weakly hyperbolic. In presenting our numerical algorithm we describe how to numerically handle this problem of weak hyperbolicity, as well as how to choose an appropriate gauge condition. The
NASA Astrophysics Data System (ADS)
Thomas, Justin W.
2006-12-01
The Numerical Nuclear Reactor (NNR) is a code suite that is being developed to provide high-fidelity multi-physics capability for the analysis of light water nuclear reactors. The focus of the work here is to extend the capability of the NNR by incorporation of the neutronics module, DeCART, for Boiling Water Reactor (BWR) applications. The DeCART code has been coupled to the NNR fluid mechanics and heat transfer module STAR-CD for light water reactor applications. The coupling has been accomplished via an interface program, which is responsible for mapping the STAR-CD and DeCART meshes, managing communication, and monitoring convergence. DeCART obtains the solution of the 3-D Boltzmann transport equation by performing a series of 2-D modular ray tracing-based method of characteristics problems that are coupled within the framework of 3-D coarse-mesh finite difference. The relatively complex geometry and increased axial heterogeneity found in BWRs are beyond the modeling capability of the original version of DeCART. In this work, DeCART is extended in three primary areas. First, the geometric capability is generalized by extending the modular ray tracing scheme and permitting an unstructured mesh in the global finite difference kernel. Second, numerical instabilities, which arose as a result of the severe axial heterogeneity found in BWR cores, have been resolved. Third, an advanced nodal method has been implemented to improve the accuracy of the axial flux distribution. In this semi-analytic nodal method, the analytic solution to the transverse-integrated neutron diffusion equation is obtained, where the nonhomogeneous neutron source was first approximated by a quartic polynomial. The successful completion of these three tasks has allowed the application of the coupled DeCART/STAR-CD code to practical BWR problems.
An iterative KP1 method for solving the transport equation in 3D domains on unstructured grids
NASA Astrophysics Data System (ADS)
Kokonkov, N. I.; Nikolaeva, O. V.
2015-10-01
A two-step iterative KP1 method for solving systems of grid equations that approximate the integro-differential transport equation in 3D domains on unstructured grids using nodal SN methods is described. Results of testing the efficiency of the proposed method in solving benchmark problems of reactor protection on tetrahedral grids are presented.
Reconstruction of 3D structure using stochastic methods: morphology and transport properties
NASA Astrophysics Data System (ADS)
Karsanina, Marina; Gerke, Kirill; Čapek, Pavel; Vasilyev, Roman; Korost, Dmitry; Skvortsova, Elena
2013-04-01
One of the main factors defining numerous flow phenomena in rocks, soils and other porous media, including fluid and solute movements, is pore structure, e.g., pore sizes and their connectivity. Numerous numerical methods were developed to quantify single and multi-phase flow in such media on microscale. Among most popular ones are: 1) a wide range of finite difference/element/volume solutions of Navier-Stokes equations and its simplifications; 2) lattice-Boltzmann method; 3) pore-network models, among others. Each method has some advantages and shortcomings, so that different research teams usually utilize more than one, depending on the study case. Recent progress in 3D imaging of internal structure, e.g., X-ray tomography, FIB-SEM and confocal microscopy, made it possible to obtain digitized input pore parameters for such models, however, a trade-off between resolution and sample size is usually unavoidable. There are situations then only standard two-dimensional information of porous structure is known due to tomography high cost or resolution limitations. However, physical modeling on microscale requires 3D information. There are three main approaches to reconstruct (using 2D cut(s) or some other limited information/properties) porous media: 1) statistical methods (correlation functions and simulated annealing, multi-point statistics, entropy methods), 2) sequential methods (sphere or other granular packs) and 3) morphological methods. Stochastic reconstructions using correlation functions possess some important advantage - they provide a statistical description of the structure, which is known to have relationships with all physical properties. In addition, this method is more flexible for other applications to characterize porous media. Taking different 3D scans of natural and artificial porous materials (sandstones, soils, shales, ceramics) we choose some 2D cut/s as sources of input correlation functions. Based on different types of correlation functions
Intercomparison of 3D pore-scale flow and solute transport simulation methods
Mehmani, Yashar; Schoenherr, Martin; Pasquali, Andrea; Perkins, William A.; Kim, Kyungjoo; Perego, Mauro; Parks, Michael L.; Balhoff, Matthew T.; Richmond, Marshall C.; Geier, Martin; et al
2015-09-28
Multiple numerical approaches have been developed to simulate porous media fluid flow and solute transport at the pore scale. These include 1) methods that explicitly model the three-dimensional geometry of pore spaces and 2) methods that conceptualize the pore space as a topologically consistent set of stylized pore bodies and pore throats. In previous work we validated a model of the first type, using computational fluid dynamics (CFD) codes employing a standard finite volume method (FVM), against magnetic resonance velocimetry (MRV) measurements of pore-scale velocities. Here we expand that validation to include additional models of the first type based onmore » the lattice Boltzmann method (LBM) and smoothed particle hydrodynamics (SPH), as well as a model of the second type, a pore-network model (PNM). The PNM approach used in the current study was recently improved and demonstrated to accurately simulate solute transport in a two-dimensional experiment. While the PNM approach is computationally much less demanding than direct numerical simulation methods, the effect of conceptualizing complex three-dimensional pore geometries on solute transport in the manner of PNMs has not been fully determined. We apply all four approaches (FVM-based CFD, LBM, SPH and PNM) to simulate pore-scale velocity distributions and (for capable codes) nonreactive solute transport, and intercompare the model results. Comparisons are drawn both in terms of macroscopic variables (e.g., permeability, solute breakthrough curves) and microscopic variables (e.g., local velocities and concentrations). Generally good agreement was achieved among the various approaches, but some differences were observed depending on the model context. The intercomparison work was challenging because of variable capabilities of the codes, and inspired some code enhancements to allow consistent comparison of flow and transport simulations across the full suite of methods. This paper provides support for
Intercomparison of 3D pore-scale flow and solute transport simulation methods
Mehmani, Yashar; Schoenherr, Martin; Pasquali, Andrea; Perkins, William A.; Kim, Kyungjoo; Perego, Mauro; Parks, Michael L.; Balhoff, Matthew T.; Richmond, Marshall C.; Geier, Martin; Krafczyk, Manfred; Luo, Li -Shi; Tartakovsky, Alexandre M.; Yang, Xiaofan; Scheibe, Timothy D.; Trask, Nathaniel
2015-09-28
Multiple numerical approaches have been developed to simulate porous media fluid flow and solute transport at the pore scale. These include 1) methods that explicitly model the three-dimensional geometry of pore spaces and 2) methods that conceptualize the pore space as a topologically consistent set of stylized pore bodies and pore throats. In previous work we validated a model of the first type, using computational fluid dynamics (CFD) codes employing a standard finite volume method (FVM), against magnetic resonance velocimetry (MRV) measurements of pore-scale velocities. Here we expand that validation to include additional models of the first type based on the lattice Boltzmann method (LBM) and smoothed particle hydrodynamics (SPH), as well as a model of the second type, a pore-network model (PNM). The PNM approach used in the current study was recently improved and demonstrated to accurately simulate solute transport in a two-dimensional experiment. While the PNM approach is computationally much less demanding than direct numerical simulation methods, the effect of conceptualizing complex three-dimensional pore geometries on solute transport in the manner of PNMs has not been fully determined. We apply all four approaches (FVM-based CFD, LBM, SPH and PNM) to simulate pore-scale velocity distributions and (for capable codes) nonreactive solute transport, and intercompare the model results. Comparisons are drawn both in terms of macroscopic variables (e.g., permeability, solute breakthrough curves) and microscopic variables (e.g., local velocities and concentrations). Generally good agreement was achieved among the various approaches, but some differences were observed depending on the model context. The intercomparison work was challenging because of variable capabilities of the codes, and inspired some code enhancements to allow consistent comparison of flow and transport simulations across the full suite of methods. This paper provides support for confidence
Benchmark Study of 3D Pore-scale Flow and Solute Transport Simulation Methods
NASA Astrophysics Data System (ADS)
Scheibe, T. D.; Yang, X.; Mehmani, Y.; Perkins, W. A.; Pasquali, A.; Schoenherr, M.; Kim, K.; Perego, M.; Parks, M. L.; Trask, N.; Balhoff, M.; Richmond, M. C.; Geier, M.; Krafczyk, M.; Luo, L. S.; Tartakovsky, A. M.
2015-12-01
Multiple numerical approaches have been developed to simulate porous media fluid flow and solute transport at the pore scale. These include 1) methods that explicitly model the three-dimensional geometry of pore spaces and 2) methods that conceptualize the pore space as a topologically consistent set of stylized pore bodies and pore throats. In previous work we validated a model of the first type, using computational fluid dynamics (CFD) codes employing standard finite volume method (FVM), against magnetic resonance velocimetry (MRV) measurements of pore-scale velocities. Here we expand that benchmark study to include additional models of the first type based on the immersed-boundary method (IMB), lattice Boltzmann method (LBM), and smoothed particle hydrodynamics (SPH), as well as a model of the second type, a pore-network model (PNM). While the PNM approach is computationally much less demanding than direct numerical simulation methods, the effect of conceptualizing complex three-dimensional pore geometries in the manner of PNMs has not been fully determined. We apply all five approaches (FVM-based CFD, IMB, LBM, SPH and PNM) to simulate pore-scale velocity distributions and nonreactive solute transport, and intercompare the model results. Comparisons are drawn both in terms of macroscopic variables (e.g., permeability, solute breakthrough curves) and microscopic variables (e.g., local velocities and concentrations). Generally good agreement was achieved among the various approaches, but some differences were observed depending on the model context. The benchmark study was challenging because of variable capabilities of the codes, and inspired some code enhancements to allow consistent comparison of flow and transport simulations across the full suite of methods. This study provides support for confidence in a variety of pore-scale modeling methods, and motivates further development and application of pore-scale simulation methods.
MT3D was first developed by Chunmiao Zheng in 1990 at S.S. Papadopulos & Associates, Inc. with partial support from the U.S. Environmental Protection Agency (USEPA). Starting in 1990, MT3D was released as a pubic domain code from the USEPA. Commercial versions with enhanced capab...
NASA Astrophysics Data System (ADS)
Wehrer, Markus; Slater, Lee
2015-04-01
flow fraction was observed to be independent of precipitation rate. This suggests the presence of a fingering process driven by textural heterogeneities. As a consequence, preferential transport of the conservative and the reactive tracer also occurred. We found that 3D ERT can serve to quantitatively characterize shape measures of both tracer breakthroughs and water content dynamics. In particular, shape measures influenced by the advective propagation of the tracer peak, like mean velocity and normalized first central moment, are highly correlated between ERT data and validation data (consisting of tracer measurements in seepage water samples). Using shape measures proved to be advantageous over interpretation of ERT data with spatially uncertain petrophysical functions for the characterization of heterogeneous flow and transport. Consequently, for future applications of ERT in soil hydrological modeling, the use of temporal moments is recommended.
3D Multigroup Sn Neutron Transport Code
Energy Science and Technology Software Center (ESTSC)
2001-02-14
ATTILA is a 3D multigroup transport code with arbitrary order ansotropic scatter. The transport equation is solved in first order form using a tri-linear discontinuous spatial differencing on an arbitrary tetrahedral mesh. The overall solution technique is source iteration with DSA acceleration of the scattering source. Anisotropic boundary and internal sources may be entered in the form of spherical harmonics moments. Alpha and k eigenvalue problems are allowed, as well as fixed source problems. Forwardmore » and adjoint solutions are available. Reflective, vacumn, and source boundary conditions are available. ATTILA can perform charged particle transport calculations using slowing down (CSD) terms. ATTILA can also be used to peform infra-red steady-state calculations for radiative transfer purposes.« less
3D Multigroup Sn Neutron Transport Code
McGee, John; Wareing, Todd; Pautz, Shawn
2001-02-14
ATTILA is a 3D multigroup transport code with arbitrary order ansotropic scatter. The transport equation is solved in first order form using a tri-linear discontinuous spatial differencing on an arbitrary tetrahedral mesh. The overall solution technique is source iteration with DSA acceleration of the scattering source. Anisotropic boundary and internal sources may be entered in the form of spherical harmonics moments. Alpha and k eigenvalue problems are allowed, as well as fixed source problems. Forward and adjoint solutions are available. Reflective, vacumn, and source boundary conditions are available. ATTILA can perform charged particle transport calculations using slowing down (CSD) terms. ATTILA can also be used to peform infra-red steady-state calculations for radiative transfer purposes.
New method of 3-D object recognition
NASA Astrophysics Data System (ADS)
He, An-Zhi; Li, Qun Z.; Miao, Peng C.
1991-12-01
In this paper, a new method of 3-D object recognition using optical techniques and a computer is presented. We perform 3-D object recognition using moire contour to obtain the object's 3- D coordinates, projecting drawings of the object in three coordinate planes to describe it and using a method of inquiring library of judgement to match objects. The recognition of a simple geometrical entity is simulated by computer and studied experimentally. The recognition of an object which is composed of a few simple geometrical entities is discussed.
Methods for comparing 3D surface attributes
NASA Astrophysics Data System (ADS)
Pang, Alex; Freeman, Adam
1996-03-01
A common task in data analysis is to compare two or more sets of data, statistics, presentations, etc. A predominant method in use is side-by-side visual comparison of images. While straightforward, it burdens the user with the task of discerning the differences between the two images. The user if further taxed when the images are of 3D scenes. This paper presents several methods for analyzing the extent, magnitude, and manner in which surfaces in 3D differ in their attributes. The surface geometry are assumed to be identical and only the surface attributes (color, texture, etc.) are variable. As a case in point, we examine the differences obtained when a 3D scene is rendered progressively using radiosity with different form factor calculation methods. The comparison methods include extensions of simple methods such as mapping difference information to color or transparency, and more recent methods including the use of surface texture, perturbation, and adaptive placements of error glyphs.
A non-conforming 3D spherical harmonic transport solver
Van Criekingen, S.
2006-07-01
A new 3D transport solver for the time-independent Boltzmann transport equation has been developed. This solver is based on the second-order even-parity form of the transport equation. The angular discretization is performed through the expansion of the angular neutron flux in spherical harmonics (PN method). The novelty of this solver is the use of non-conforming finite elements for the spatial discretization. Such elements lead to a discontinuous flux approximation. This interface continuity requirement relaxation property is shared with mixed-dual formulations such as the ones based on Raviart-Thomas finite elements. Encouraging numerical results are presented. (authors)
NASA Technical Reports Server (NTRS)
Bidwell, Colin S.; Pinella, David; Garrison, Peter
1999-01-01
Collection efficiency and ice accretion calculations were made for a commercial transport using the NASA Lewis LEWICE3D ice accretion code, the ICEGRID3D grid code and the CMARC panel code. All of the calculations were made on a Windows 95 based personal computer. The ice accretion calculations were made for the nose, wing, horizontal tail and vertical tail surfaces. Ice shapes typifying those of a 30 minute hold were generated. Collection efficiencies were also generated for the entire aircraft using the newly developed unstructured collection efficiency method. The calculations highlight the flexibility and cost effectiveness of the LEWICE3D, ICEGRID3D, CMARC combination.
Recognition methods for 3D textured surfaces
NASA Astrophysics Data System (ADS)
Cula, Oana G.; Dana, Kristin J.
2001-06-01
Texture as a surface representation is the subject of a wide body of computer vision and computer graphics literature. While texture is always associated with a form of repetition in the image, the repeating quantity may vary. The texture may be a color or albedo variation as in a checkerboard, a paisley print or zebra stripes. Very often in real-world scenes, texture is instead due to a surface height variation, e.g. pebbles, gravel, foliage and any rough surface. Such surfaces are referred to here as 3D textured surfaces. Standard texture recognition algorithms are not appropriate for 3D textured surfaces because the appearance of these surfaces changes in a complex manner with viewing direction and illumination direction. Recent methods have been developed for recognition of 3D textured surfaces using a database of surfaces observed under varied imaging parameters. One of these methods is based on 3D textons obtained using K-means clustering of multiscale feature vectors. Another method uses eigen-analysis originally developed for appearance-based object recognition. In this work we develop a hybrid approach that employs both feature grouping and dimensionality reduction. The method is tested using the Columbia-Utrecht texture database and provides excellent recognition rates. The method is compared with existing recognition methods for 3D textured surfaces. A direct comparison is facilitated by empirical recognition rates from the same texture data set. The current method has key advantages over existing methods including requiring less prior information on both the training and novel images.
Cheng, Candong; Lee, Joon-Ho; Lim, Kim Hwa; Massoud, Hisham Z.; Liu, Qing Huo
2007-01-01
A 3-D quantum transport solver based on the spectral element method (SEM) and perfectly matched layer (PML) is introduced to solve the 3-D Schrödinger equation with a tensor effective mass. In this solver, the influence of the environment is replaced with the artificial PML open boundary extended beyond the contact regions of the device. These contact regions are treated as waveguides with known incident waves from waveguide mode solutions. As the transmitted wave function is treated as a total wave, there is no need to decompose it into waveguide modes, thus significantly simplifying the problem in comparison with conventional open boundary conditions. The spectral element method leads to an exponentially improving accuracy with the increase in the polynomial order and sampling points. The PML region can be designed such that less than −100 dB outgoing waves are reflected by this artificial material. The computational efficiency of the SEM solver is demonstrated by comparing the numerical and analytical results from waveguide and plane-wave examples, and its utility is illustrated by multiple-terminal devices and semiconductor nanotube devices. PMID:18037971
3D Shape and Indirect Appearance by Structured Light Transport.
OToole, Matthew; Mather, John; Kutulakos, Kiriakos N
2016-07-01
We consider the problem of deliberately manipulating the direct and indirect light flowing through a time-varying, general scene in order to simplify its visual analysis. Our approach rests on a crucial link between stereo geometry and light transport: while direct light always obeys the epipolar geometry of a projector-camera pair, indirect light overwhelmingly does not. We show that it is possible to turn this observation into an imaging method that analyzes light transport in real time in the optical domain, prior to acquisition. This yields three key abilities that we demonstrate in an experimental camera prototype: (1) producing a live indirect-only video stream for any scene, regardless of geometric or photometric complexity; (2) capturing images that make existing structured-light shape recovery algorithms robust to indirect transport; and (3) turning them into one-shot methods for dynamic 3D shape capture. PMID:27295455
Challenges in Lagrangian transport and predictability in 3D flows
NASA Astrophysics Data System (ADS)
Branicki, M.; Wiggins, S.; Kirwan, A. D.; Malek-Madani, R.
2011-12-01
The interplay between the geometrical theory of dynamical systems and the trajectory-based description of aperiodically time-dependent fluid flows has led to significant advances in understanding the role of chaotic transport in geophysical flows at scales dominated by advection. Lagrangian transport analysis utilizing either the time-dependent geometry of intersecting stable and unstable manifolds of the so-called Distinguished Hyperbolic Trajectories (DHT), or ridges of finite-time Lyapunov exponent fields (LCS), provide a much needed and complementary insight into ephemeral mechanisms responsible for the existence of `leaky' transport barriers and 'leaky' mesoscale eddies. However, to date most oceanic applications have been confined to 2D analysis of near surface regions in 'perfect' flows not accounting for model or measurement error, and with the tacit assumption of negligible vertical velocities. I will systematically address issues concerning the regimes of applicability of two-dimensional analysis in 3D aperiodically time-dependent flows, as well as outstanding challenges in fully 3D Lagrangian transport analysis. Even for perfect horizontal velocities, little is known about the vertical extent of stable/unstable manifolds associated with DHTs and/or other special structures relevant to stratified 3D flows. In particular, their sensitivity to errors in the vertical velocities and data assimilation methods has been little studied. Rigorous results regarding the above issues will be illustrated by revealing and mathematically tractable toy models, as well as examples from a detailed study in an eddy-rich region from the Gulf of Mexico and the Mediterranean. New ways of quantifying the uncertainty in Lagrangian predictions will also be presented.
Quantum transport through 3D Dirac materials
Salehi, M.; Jafari, S.A.
2015-08-15
Bismuth and its alloys provide a paradigm to realize three dimensional materials whose low-energy effective theory is given by Dirac equation in 3+1 dimensions. We study the quantum transport properties of three dimensional Dirac materials within the framework of Landauer–Büttiker formalism. Charge carriers in normal metal satisfying the Schrödinger equation, can be split into four-component with appropriate matching conditions at the boundary with the three dimensional Dirac material (3DDM). We calculate the conductance and the Fano factor of an interface separating 3DDM from a normal metal, as well as the conductance through a slab of 3DDM. Under certain circumstances the 3DDM appears transparent to electrons hitting the 3DDM. We find that electrons hitting the metal-3DDM interface from metallic side can enter 3DDM in a reversed spin state as soon as their angle of incidence deviates from the direction perpendicular to interface. However the presence of a second interface completely cancels this effect.
Quantum transport through 3D Dirac materials
NASA Astrophysics Data System (ADS)
Salehi, M.; Jafari, S. A.
2015-08-01
Bismuth and its alloys provide a paradigm to realize three dimensional materials whose low-energy effective theory is given by Dirac equation in 3+1 dimensions. We study the quantum transport properties of three dimensional Dirac materials within the framework of Landauer-Büttiker formalism. Charge carriers in normal metal satisfying the Schrödinger equation, can be split into four-component with appropriate matching conditions at the boundary with the three dimensional Dirac material (3DDM). We calculate the conductance and the Fano factor of an interface separating 3DDM from a normal metal, as well as the conductance through a slab of 3DDM. Under certain circumstances the 3DDM appears transparent to electrons hitting the 3DDM. We find that electrons hitting the metal-3DDM interface from metallic side can enter 3DDM in a reversed spin state as soon as their angle of incidence deviates from the direction perpendicular to interface. However the presence of a second interface completely cancels this effect.
Puso, M A; Laursen, T A
2002-05-02
Smoothing of contact surfaces can be used to eliminate the chatter typically seen with node on facet contact and give a better representation of the actual contact surface. The latter affect is well demonstrated for problems with interference fits. In this work we present two methods for the smoothing of contact surfaces for 3D finite element contact. In the first method, we employ Gregory patches to smooth the faceted surface in a node on facet implementation. In the second method, we employ a Bezier interpolation of the faceted surface in a mortar method implementation of contact. As is well known, node on facet approaches can exhibit locking due to the failure of the Babuska-Brezzi condition and in some instances fail the patch test. The mortar method implementation is stable and provides optimal convergence in the energy of error. In the this work we demonstrate the superiority of the smoothed versus the non-smoothed node on facet implementations. We also show where the node on facet method fails and some results from the smoothed mortar method implementation.
Transport of 3D space charge dominated beams
NASA Astrophysics Data System (ADS)
Lü, Jian-Qin
2013-10-01
In this paper we present the theoretical analysis and the computer code design for the intense pulsed beam transport. Intense beam dynamics is a very important issue in low-energy high-current accelerators and beam transport systems. This problem affects beam transmission and beam qualities. Therefore, it attracts the attention of the accelerator physicists worldwide. The analysis and calculation for the intense beam dynamics are very complicated, because the state of particle motion is dominated not only by the applied electromagnetic fields, but also by the beam-induced electromagnetic fields (self-fields). Moreover, the self fields are related to the beam dimensions and particle distributions. So, it is very difficult to get the self-consistent solutions of particle motion analytically. For this reason, we combine the Lie algebraic method and the particle in cell (PIC) scheme together to simulate intense 3D beam transport. With the Lie algebraic method we analyze the particle nonlinear trajectories in the applied electromagnetic fields up to third order approximation, and with the PIC algorithm we calculate the space charge effects to the particle motion. Based on the theoretical analysis, we have developed a computer code, which calculates beam transport systems consisting of electrostatic lenses, electrostatic accelerating columns, solenoid lenses, magnetic and electric quadruples, magnetic sextupoles, octopuses and different kinds of electromagnetic analyzers. The optimization calculations and the graphic display for the calculated results are provided by the code.
Bekar, Kursat B; Azmy, Yousry
2009-01-01
Improved TORT solutions to the 3D transport codes' suite of benchmarks exercise are presented in this study. Preliminary TORT solutions to this benchmark indicate that the majority of benchmark quantities for most benchmark cases are computed with good accuracy, and that accuracy improves with model refinement. However, TORT fails to compute accurate results for some benchmark cases with aspect ratios drastically different from 1, possibly due to ray effects. In this work, we employ the standard approach of splitting the solution to the transport equation into an uncollided flux and a fully collided flux via the code sequence GRTUNCL3D and TORT to mitigate ray effects. The results of this code sequence presented in this paper show that the accuracy of most benchmark cases improved substantially. Furthermore, the iterative convergence problems reported for the preliminary TORT solutions have been resolved by bringing the computational cells' aspect ratio closer to unity and, more importantly, by using 64-bit arithmetic precision in the calculation sequence. Results of this study are also reported.
3D Face Modeling Using the Multi-Deformable Method
Hwang, Jinkyu; Yu, Sunjin; Kim, Joongrock; Lee, Sangyoun
2012-01-01
In this paper, we focus on the problem of the accuracy performance of 3D face modeling techniques using corresponding features in multiple views, which is quite sensitive to feature extraction errors. To solve the problem, we adopt a statistical model-based 3D face modeling approach in a mirror system consisting of two mirrors and a camera. The overall procedure of our 3D facial modeling method has two primary steps: 3D facial shape estimation using a multiple 3D face deformable model and texture mapping using seamless cloning that is a type of gradient-domain blending. To evaluate our method's performance, we generate 3D faces of 30 individuals and then carry out two tests: accuracy test and robustness test. Our method shows not only highly accurate 3D face shape results when compared with the ground truth, but also robustness to feature extraction errors. Moreover, 3D face rendering results intuitively show that our method is more robust to feature extraction errors than other 3D face modeling methods. An additional contribution of our method is that a wide range of face textures can be acquired by the mirror system. By using this texture map, we generate realistic 3D face for individuals at the end of the paper. PMID:23201976
Turbulence and transport in a 3D magnetic boundary
NASA Astrophysics Data System (ADS)
Agostini, Matteo; Carraro, Lorella; Ciaccio, Giovanni; de Masi, Gianluca; Rea, Cristina; Scarin, Paolo; Spizzo, Gianluca; Spolaore, Monica; Vianello, Nicola
2014-10-01
In present fusion devices the interaction between 3D magnetic field, edge kinetic properties and turbulence is a crucial issue; not only in intrinsically 3D configurations such as the stellarators, but also in tokamaks, where magnetic perturbations are applied to control ELMs and plasma wall interaction. In the RFX-mod reversed field pinch the spontaneous development at high plasma current of a helical magnetic state displays strong analogies with the aforementioned configurations. At the edge the presence of a stochastic layer and magnetic islands with a well-defined helical symmetry leads to a helical pattern of flow, pressure gradients and turbulent fluctuations: larger fluctuations and shorter correlation lengths are observed near the X-point of the magnetic island, where also a flow slowing-down occurs. Aim of this work is to study the effect of edge turbulence on particle transport in a 3D magnetic boundary, characterizing the properties of the edge blobs along the helical deformation. The magnetic topology also modifies kinetic properties, with higher pressure gradients observed close to the O-point of the island. The measurement of the time evolution of pressure gradient and blob characteristics, can clarify the mutual relation between these two quantities.
New Algorithms for Large-scale 3D Radiation Transport
NASA Astrophysics Data System (ADS)
Lentz, Eric J.
2009-05-01
Radiation transport is critical not only for analysis of astrophysical objects but also for the dynamical transport of energy within. Increased fidelity and dimensionality of the other components of such models requires a similar improvement in the radiation transport. Modern astrophysical simulations can be large enough that the values for a single variable for the entire computational domain cannot be stored on a single compute node. The natural solution is to decompose the physical domain into pieces with each node responsible for a single sub-domain. Using localized plus "ghost" zone data works well for problems like explicit hydrodynamics or nuclear reaction networks with modest impact from inter-process communication. Unfortunately, radiation transport is an inherently non-local process that couples the entire model domain together and efficient algorithms are needed to conquer this problem. In this poster, I present the early development of a new parallel, 3-D transport code using ray tracing to formally solve the transport equation across numerically decomposed domains. The algorithm model takes advantage of one-sided communication to develop a scalable, parallel formal solver. Other aspects and future direction of the parallel code development such as scalability and the inclusion of scattering will also be discussed.
a Fast Method for Measuring the Similarity Between 3d Model and 3d Point Cloud
NASA Astrophysics Data System (ADS)
Zhang, Zongliang; Li, Jonathan; Li, Xin; Lin, Yangbin; Zhang, Shanxin; Wang, Cheng
2016-06-01
This paper proposes a fast method for measuring the partial Similarity between 3D Model and 3D point Cloud (SimMC). It is crucial to measure SimMC for many point cloud-related applications such as 3D object retrieval and inverse procedural modelling. In our proposed method, the surface area of model and the Distance from Model to point Cloud (DistMC) are exploited as measurements to calculate SimMC. Here, DistMC is defined as the weighted distance of the distances between points sampled from model and point cloud. Similarly, Distance from point Cloud to Model (DistCM) is defined as the average distance of the distances between points in point cloud and model. In order to reduce huge computational burdens brought by calculation of DistCM in some traditional methods, we define SimMC as the ratio of weighted surface area of model to DistMC. Compared to those traditional SimMC measuring methods that are only able to measure global similarity, our method is capable of measuring partial similarity by employing distance-weighted strategy. Moreover, our method is able to be faster than other partial similarity assessment methods. We demonstrate the superiority of our method both on synthetic data and laser scanning data.
The COMET method in 3-D hexagonal geometry
Connolly, K. J.; Rahnema, F.
2012-07-01
The hybrid stochastic-deterministic coarse mesh radiation transport (COMET) method developed at Georgia Tech now solves reactor core problems in 3-D hexagonal geometry. In this paper, the method is used to solve three preliminary test problems designed to challenge the method with steep flux gradients, high leakage, and strong asymmetry and heterogeneity in the core. The test problems are composed of blocks taken from a high temperature test reactor benchmark problem. As the method is still in development, these problems and their results are strictly preliminary. Results are compared to whole core Monte Carlo reference solutions in order to verify the method. Relative errors are on the order of 50 pcm in core eigenvalue, and mean relative error in pin fission density calculations is less than 1% in these difficult test cores. The method requires the one-time pre-computation of a response expansion coefficient library, which may be compiled in a comparable amount of time to a single whole core Monte Carlo calculation. After the library has been computed, COMET may solve any number of core configurations on the order of an hour, representing a significant gain in efficiency over other methods for whole core transport calculations. (authors)
An implicit dispersive transport algorithm for the US Geological Survey MOC3D solute-transport model
Kipp, K.L., Jr.; Konikow, L.F.; Hornberger, G.Z.
1998-01-01
This report documents an extension to the U.S. Geological Survey MOC3D transport model that incorporates an implicit-in-time difference approximation for the dispersive transport equation, including source/sink terms. The original MOC3D transport model (Version 1) uses the method of characteristics to solve the transport equation on the basis of the velocity field. The original MOC3D solution algorithm incorporates particle tracking to represent advective processes and an explicit finite-difference formulation to calculate dispersive fluxes. The new implicit procedure eliminates several stability criteria required for the previous explicit formulation. This allows much larger transport time increments to be used in dispersion-dominated problems. The decoupling of advective and dispersive transport in MOC3D, however, is unchanged. With the implicit extension, the MOC3D model is upgraded to Version 2. A description of the numerical method of the implicit dispersion calculation, the data-input requirements and output options, and the results of simulator testing and evaluation are presented. Version 2 of MOC3D was evaluated for the same set of problems used for verification of Version 1. These test results indicate that the implicit calculation of Version 2 matches the accuracy of Version 1, yet is more efficient than the explicit calculation for transport problems that are characterized by a grid Peclet number less than about 1.0.
3D scanning modeling method application in ancient city reconstruction
NASA Astrophysics Data System (ADS)
Ren, Pu; Zhou, Mingquan; Du, Guoguang; Shui, Wuyang; Zhou, Pengbo
2015-07-01
With the development of optical engineering technology, the precision of 3D scanning equipment becomes higher, and its role in 3D modeling is getting more distinctive. This paper proposed a 3D scanning modeling method that has been successfully applied in Chinese ancient city reconstruction. On one hand, for the existing architectures, an improved algorithm based on multiple scanning is adopted. Firstly, two pieces of scanning data were rough rigid registered using spherical displacers and vertex clustering method. Secondly, a global weighted ICP (iterative closest points) method is used to achieve a fine rigid registration. On the other hand, for the buildings which have already disappeared, an exemplar-driven algorithm for rapid modeling was proposed. Based on the 3D scanning technology and the historical data, a system approach was proposed for 3D modeling and virtual display of ancient city.
How do hydrodynamic instabilities affect 3D transport in geophysical vortices?
NASA Astrophysics Data System (ADS)
Wang, Peng; Özgökmen, Tamay M.
2015-03-01
Three-dimensional (3D) transport within geophysical vortices (e.g. ocean eddies) is important in understanding processes at a variety of scales, ranging from plankton production to climate variability. 3D transport can be affected by hydrodynamic instabilities of geophysical vortices; however, how the instabilities affecting 3D transport is not clear. Focusing on barotropic, inertial and 3D instabilities, we investigate the joint impacts of instabilities on 3D transport by using analytical methods and direct numerical simulations. We discover for the first time that material can be exchanged through 3D pathways which link a family of vortices generated by the instabilities in a single, initially unstable vortex. We also show that instabilities can increase the magnitude of vertical velocity, mixing rate and vertical material exchange. Besides, we find that instabilities can cause the kinetic energy wavenumber spectrum to have a power-law regime different than the classic regimes of k - 5 / 3 and k-3, and propose a new energy spectrum to interpret the non-classic regime.
Moving from Batch to Field Using the RT3D Reactive Transport Modeling System
NASA Astrophysics Data System (ADS)
Clement, T. P.; Gautam, T. R.
2002-12-01
The public domain reactive transport code RT3D (Clement, 1997) is a general-purpose numerical code for solving coupled, multi-species reactive transport in saturated groundwater systems. The code uses MODFLOW to simulate flow and several modules of MT3DMS to simulate the advection and dispersion processes. RT3D employs the operator-split strategy which allows the code solve the coupled reactive transport problem in a modular fashion. The coupling between reaction and transport is defined through a separate module where the reaction equations are specified. The code supports a versatile user-defined reaction option that allows users to define their own reaction system through a Fortran-90 subroutine, known as the RT3D-reaction package. Further a utility code, known as BATCHRXN, allows the users to independently test and debug their reaction package. To analyze a new reaction system at a batch scale, users should first run BATCHRXN to test the ability of their reaction package to model the batch data. After testing, the reaction package can simply be ported to the RT3D environment to study the model response under 1-, 2-, or 3-dimensional transport conditions. This paper presents example problems that demonstrate the methods for moving from batch to field-scale simulations using BATCHRXN and RT3D codes. The first example describes a simple first-order reaction system for simulating the sequential degradation of Tetrachloroethene (PCE) and its daughter products. The second example uses a relatively complex reaction system for describing the multiple degradation pathways of Tetrachloroethane (PCA) and its daughter products. References 1) Clement, T.P, RT3D - A modular computer code for simulating reactive multi-species transport in 3-Dimensional groundwater aquifers, Battelle Pacific Northwest National Laboratory Research Report, PNNL-SA-28967, September, 1997. Available at: http://bioprocess.pnl.gov/rt3d.htm.
Light Attenuation Method for 3D data acquisition (LAM3D) of bottom particle deposits
NASA Astrophysics Data System (ADS)
Er, Jenn Wei; Law, Adrian W. K.; Adams, E. Eric; Yang, Yang
2015-11-01
We have developed a novel experimental technique, Light Attenuation Method for 3D data acquisition (LAM3D), to acquire three-dimensional spatial characteristics and temporal development of bottom particle deposits. The new technique performs data acquisition with higher spatial and temporal resolution than existing approaches with laser and ultrasonic 3D profilers, and is therefore ideal for laboratory investigations with fast varying changes in the sediment bed, such as the developing deposition profile from sediment clouds commonly formed during dredging or land reclamation projects and the dynamic evolution in movable bed processes in rivers. The principle of the technique is based on the analysis of the light attenuation due to multiple light scattering through the particle deposits layer compared to the clear water column. With appropriate calibration, the particles size and distribution thickness can be quantified by the transmitted light spectrum. In the presentation, we will first show our measurement setup with a light panel for calibrated illumination and a system of DSLR cameras for the light capturing. Subsequently, we shall present the experimental results of fast evolving deposition profile of a barge-disposed sediment cloud upon its bottom impact on the sea bed.
A 3D Level Set Method for Microwave Breast Imaging
Colgan, Timothy J.; Hagness, Susan C.; Van Veen, Barry D.
2015-01-01
Objective Conventional inverse-scattering algorithms for microwave breast imaging result in moderate resolution images with blurred boundaries between tissues. Recent 2D numerical microwave imaging studies demonstrate that the use of a level set method preserves dielectric boundaries, resulting in a more accurate, higher resolution reconstruction of the dielectric properties distribution. Previously proposed level set algorithms are computationally expensive and thus impractical in 3D. In this paper we present a computationally tractable 3D microwave imaging algorithm based on level sets. Methods We reduce the computational cost of the level set method using a Jacobian matrix, rather than an adjoint method, to calculate Frechet derivatives. We demonstrate the feasibility of 3D imaging using simulated array measurements from 3D numerical breast phantoms. We evaluate performance by comparing full 3D reconstructions to those from a conventional microwave imaging technique. We also quantitatively assess the efficacy of our algorithm in evaluating breast density. Results Our reconstructions of 3D numerical breast phantoms improve upon those of a conventional microwave imaging technique. The density estimates from our level set algorithm are more accurate than those of conventional microwave imaging, and the accuracy is greater than that reported for mammographic density estimation. Conclusion Our level set method leads to a feasible level of computational complexity for full 3D imaging, and reconstructs the heterogeneous dielectric properties distribution of the breast more accurately than conventional microwave imaging methods. Significance 3D microwave breast imaging using a level set method is a promising low-cost, non-ionizing alternative to current breast imaging techniques. PMID:26011863
3D face recognition based on a modified ICP method
NASA Astrophysics Data System (ADS)
Zhao, Kankan; Xi, Jiangtao; Yu, Yanguang; Chicharo, Joe F.
2011-11-01
3D face recognition technique has gained much more attention recently, and it is widely used in security system, identification system, and access control system, etc. The core technique in 3D face recognition is to find out the corresponding points in different 3D face images. The classic partial Iterative Closest Point (ICP) method is iteratively align the two point sets based on repetitively calculate the closest points as the corresponding points in each iteration. After several iterations, the corresponding points can be obtained accurately. However, if two 3D face images with different scale are from the same person, the classic partial ICP does not work. In this paper we propose a modified partial Iterative Closest Point (ICP) method in which the scaling effect is considered to achieve 3D face recognition. We design a 3x3 diagonal matrix as the scale matrix in each iteration of the classic partial ICP. The probing face image which is multiplied by the scale matrix will keep the similar scale with the reference face image. Therefore, we can accurately determine the corresponding points even the scales of probing image and reference image are different. 3D face images in our experiments are acquired by a 3D data acquisition system based on Digital Fringe Projection Profilometry (DFPP). A 3D database consists of 30 group images, three images with the same scale, which are from the same person with different views, are included in each group. And in different groups, the scale of the 3 images may be different from other groups. The experiment results show that our proposed method can achieve 3D face recognition, especially in the case that the scales of probing image and referent image are different.
3-D numerical simulations of volcanic ash transport and deposition
NASA Astrophysics Data System (ADS)
Suzuki, Y. J.; Koyaguchi, T.
2012-12-01
During an explosive volcanic eruption, volcanic gas and pyroclasts are ejected from the volcanic vent. The pyroclasts are carried up within a convective plume, advected by the surrounding wind field, and sediment on the ground depending on their terminal velocity. The fine ash are expected to have atmospheric residence, whereas the coarser particles form fall deposits. Accurate modeling of particle transport and deposition is of critical importance from the viewpoint of disaster prevention. Previously, some particle-tracking models (e.g., PUFF) and advection-diffusion models (e.g., TEPHRA2 and FALL3D) tried to forecast particle concentration in the atmosphere and particle loading at ground level. However, these models assumed source conditions (the grain-size distribution, plume height, and mass release location) based on the simple 1-D model of convective plume. In this study, we aim to develop a new 3-D model which reproduces both of the dynamics of convective plume and the ash transport. The model is designed to describe the injection of eruption cloud and marker particles from a circular vent above a flat surface into the stratified atmosphere. Because the advection is the predominant mechanism of particle transport near the volcano, the diffusive process is not taken into account in this model. The distribution of wind velocity is given as an initial condition. The model of the eruption cloud dynamics is based on the 3-D time-dependent model of Suzuki et al. (2005). We apply a pseudo-gas model to calculate the eruption cloud dynamics: the effect of particle separation on the cloud dynamics is not considered. In order to reproduce the drastic change of eruption cloud density, we change the effective gas constant and heat capacity of the mixture in the equation of state for ideal gases with the mixing ratio between the ejected material and entrained air. In order to calculate the location and movement of ash particles, the present model employs Lagrangian marker
The Transient 3-D Transport Coupled Code TORT-TD/ATTICA3D for High-Fidelity Pebble-Bed HTGR Analyses
NASA Astrophysics Data System (ADS)
Seubert, Armin; Sureda, Antonio; Lapins, Janis; Bader, Johannes; Laurien, Eckart
2012-01-01
This article describes the 3D discrete ordinates-based coupled code system TORT-TD/ATTICA3D that aims at steady state and transient analyses of pebble-bed high-temperature gas cooled reactors. In view of increasing computing power, the application of time-dependent neutron transport methods becomes feasible for best estimate evaluations of safety margins. The calculation capabilities of TORT-TD/ATTICA3D are presented along with the coupling approach, with focus on the time-dependent neutron transport features of TORT-TD. Results obtained for the OECD/NEA/NSC PBMR-400 benchmark demonstrate the transient capabilities of TORT-TD/ATTICA3D.
Reactive Transport from Path3D: A Stream Tube Approach for Heterogeneous Aquifers
NASA Astrophysics Data System (ADS)
LI, L.
2001-05-01
Path3D (Zheng, 1991) is a popular computer program run in series with MODFLOW. Remediation engineers and hydrogeologist use it to track contaminant paths and to estimate solute travel time at heterogeneous sites. In order to predict fate and transport of multiple species at heterogeneous sites, numerical modeling packages, such as MT3D (Zheng,1990) or RT3D (Clement 1997) are often employed. These packages also are based on pre-processing with MODFLOW. However, for complex cases with aquifer heterogeneity, MT3D and RT3D often require very long computer run times. This paper addresses a new, stream-tube, approach that is both highly efficient and accurate to predict multi-species reactive transport at heterogeneous sites with steady flow. Our application of the stream tube approach is different from other stream tube approaches that apply the advection-dispersion-reaction equation in each stream tube (such as Ginn, 2000, Yabusaki, 1998, Charbeneau, 2000). In this work, the authors make use of properties of a linear system,working with decoupled reaction and sorption processes and mixing processes described by residence time distributions (RTDs). RTDs are abstracted from Path3D particle-tracking results and additional temporal and spatial dispersion (not caused by aquifer heterogeneity) is ignored. Reactions, including first order reactions and linear, reversible sorption, are applied through analytical transfer functions (called kinetic response functions). Convolution can then be applied to determine contaminant concentrations at monitoring points, using the RTDs determined from Path3D, kinetic transfer functions (expressed analytically), and expected trends of the source concentration. We are currently testing the approach and noting significant computational advantages for problems in three-dimensions, with first order reaction pathways and different retardation factors. We will demonstrate the method with several examples and compare the performance with MT3D and
2D/1D approximations to the 3D neutron transport equation. I: Theory
Kelley, B. W.; Larsen, E. W.
2013-07-01
A new class of '2D/1D' approximations is proposed for the 3D linear Boltzmann equation. These approximate equations preserve the exact transport physics in the radial directions x and y and diffusion physics in the axial direction z. Thus, the 2D/1D equations are more accurate approximations of the 3D Boltzmann equation than the conventional 3D diffusion equation. The 2D/1D equations can be systematically discretized, to yield accurate simulation methods for 3D reactor core problems. The resulting solutions will be more accurate than 3D diffusion solutions, and less expensive to generate than standard 3D transport solutions. In this paper, we (i) show that the simplest 2D/1D equation has certain desirable properties, (ii) systematically discretize this equation, and (iii) derive a stable iteration scheme for solving the discrete system of equations. In a companion paper [1], we give numerical results that confirm the theoretical predictions of accuracy and iterative stability. (authors)
[An integrated segmentation method for 3D ultrasound carotid artery].
Yang, Xin; Wu, Huihui; Liu, Yang; Xu, Hongwei; Liang, Huageng; Cai, Wenjuan; Fang, Mengjie; Wang, Yujie
2013-07-01
An integrated segmentation method for 3D ultrasound carotid artery was proposed. 3D ultrasound image was sliced into transverse, coronal and sagittal 2D images on the carotid bifurcation point. Then, the three images were processed respectively, and the carotid artery contours and thickness were obtained finally. This paper tries to overcome the disadvantages of current computer aided diagnosis method, such as high computational complexity, easily introduced subjective errors et al. The proposed method could get the carotid artery overall information rapidly, accurately and completely. It could be transplanted into clinical usage for atherosclerosis diagnosis and prevention. PMID:24195385
How Do Hydrodynamic Instabilities Affect 3D Transport in Geophysical Vortices?
NASA Astrophysics Data System (ADS)
Wang, P.; Ozgokmen, T. M.
2014-12-01
Understanding three-dimensional (3D) transport in ocean eddies is important for processes at a variety of scales, ranging from plankton production to climate variability. It is well known that geophysical vortices are subject to various hydrodynamic instabilities. Yet the influence of these instabilities on 3D material transport in vortex systems is not well investigated. Focusing on barotropic, inertial and 3D instabilities, we analyze these instabilities with normal-mode method, and reproduce their characteristics via highly-resolved numerical simulations using a spectral element Navier-Stokes solver. By comparing the simulation results of stable and unstable vortices, we investigate the joint impacts of instabilities on 3D transport through three major aspects: (i) energy transfer, (ii) overturning transport of the secondary circulation, and (iii) rates of vertical exchange and mixing. It is found that instabilities can enhance local nonlinear interactions and cause the kinetic energy wavenumber spectrum to have slopes between the conventional -5/3 and -3 at inertial ranges. The cascade of a new quantity is proposed to explain these non-conventional slopes. One of our main results is the discovery of material exchange between the central vortex and satellite vortices through 3D pathways, called funnels. These funnels modify the concept of elliptic regions that can trap material when confined to 2D dynamics. Thus, we show that a family of vortices, created by the hydrodynamic instabilities of the initially unstable vortex, can still continue to operate in unity in order to complete the 3D transport in these systems. We also show that flow instabilities can double the magnitude of vertical velocity, increase the rate of vertical exchange by an order of magnitude and enhance mixing rate more than 100%.
Improving automated 3D reconstruction methods via vision metrology
NASA Astrophysics Data System (ADS)
Toschi, Isabella; Nocerino, Erica; Hess, Mona; Menna, Fabio; Sargeant, Ben; MacDonald, Lindsay; Remondino, Fabio; Robson, Stuart
2015-05-01
This paper aims to provide a procedure for improving automated 3D reconstruction methods via vision metrology. The 3D reconstruction problem is generally addressed using two different approaches. On the one hand, vision metrology (VM) systems try to accurately derive 3D coordinates of few sparse object points for industrial measurement and inspection applications; on the other, recent dense image matching (DIM) algorithms are designed to produce dense point clouds for surface representations and analyses. This paper strives to demonstrate a step towards narrowing the gap between traditional VM and DIM approaches. Efforts are therefore intended to (i) test the metric performance of the automated photogrammetric 3D reconstruction procedure, (ii) enhance the accuracy of the final results and (iii) obtain statistical indicators of the quality achieved in the orientation step. VM tools are exploited to integrate their main functionalities (centroid measurement, photogrammetric network adjustment, precision assessment, etc.) into the pipeline of 3D dense reconstruction. Finally, geometric analyses and accuracy evaluations are performed on the raw output of the matching (i.e. the point clouds) by adopting a metrological approach. The latter is based on the use of known geometric shapes and quality parameters derived from VDI/VDE guidelines. Tests are carried out by imaging the calibrated Portable Metric Test Object, designed and built at University College London (UCL), UK. It allows assessment of the performance of the image orientation and matching procedures within a typical industrial scenario, characterised by poor texture and known 3D/2D shapes.
Hong, X; Gao, H
2014-06-15
Purpose: The Linear Boltzmann Transport Equation (LBTE) solved through statistical Monte Carlo (MC) method provides the accurate dose calculation in radiotherapy. This work is to investigate the alternative way for accurately solving LBTE using deterministic numerical method due to its possible advantage in computational speed from MC. Methods: Instead of using traditional spherical harmonics to approximate angular scattering kernel, our deterministic numerical method directly computes angular scattering weights, based on a new angular discretization method that utilizes linear finite element method on the local triangulation of unit angular sphere. As a Result, our angular discretization method has the unique advantage in positivity, i.e., to maintain all scattering weights nonnegative all the time, which is physically correct. Moreover, our method is local in angular space, and therefore handles the anisotropic scattering well, such as the forward-peaking scattering. To be compatible with image-guided radiotherapy, the spatial variables are discretized on the structured grid with the standard diamond scheme. After discretization, the improved sourceiteration method is utilized for solving the linear system without saving the linear system to memory. The accuracy of our 3D solver is validated using analytic solutions and benchmarked with Geant4, a popular MC solver. Results: The differences between Geant4 solutions and our solutions were less than 1.5% for various testing cases that mimic the practical cases. More details are available in the supporting document. Conclusion: We have developed a 3D LBTE solver based on a new angular discretization method that guarantees the positivity of scattering weights for physical correctness, and it has been benchmarked with Geant4 for photon dose calculation.
PHT3D-UZF: A Reactive Transport Model for Variably-Saturated Porous Media.
Wu, Ming Zhi; Post, Vincent E A; Salmon, S Ursula; Morway, Eric D; Prommer, Henning
2016-01-01
A modified version of the MODFLOW/MT3DMS-based reactive transport model PHT3D was developed to extend current reactive transport capabilities to the variably-saturated component of the subsurface system and incorporate diffusive reactive transport of gaseous species. Referred to as PHT3D-UZF, this code incorporates flux terms calculated by MODFLOW's unsaturated-zone flow (UZF1) package. A volume-averaged approach similar to the method used in UZF-MT3DMS was adopted. The PHREEQC-based computation of chemical processes within PHT3D-UZF in combination with the analytical solution method of UZF1 allows for comprehensive reactive transport investigations (i.e., biogeochemical transformations) that jointly involve saturated and unsaturated zone processes. Intended for regional-scale applications, UZF1 simulates downward-only flux within the unsaturated zone. The model was tested by comparing simulation results with those of existing numerical models. The comparison was performed for several benchmark problems that cover a range of important hydrological and reactive transport processes. A 2D simulation scenario was defined to illustrate the geochemical evolution following dewatering in a sandy acid sulfate soil environment. Other potential applications include the simulation of biogeochemical processes in variably-saturated systems that track the transport and fate of agricultural pollutants, nutrients, natural and xenobiotic organic compounds and micropollutants such as pharmaceuticals, as well as the evolution of isotope patterns. PMID:25628017
PHT3D-UZF: A reactive transport model for variably-saturated porous media
Wu, Ming Zhi; Post, Vincent E. A.; Salmon, S. Ursula; Morway, Eric; Prommer, H.
2016-01-01
A modified version of the MODFLOW/MT3DMS-based reactive transport model PHT3D was developed to extend current reactive transport capabilities to the variably-saturated component of the subsurface system and incorporate diffusive reactive transport of gaseous species. Referred to as PHT3D-UZF, this code incorporates flux terms calculated by MODFLOW's unsaturated-zone flow (UZF1) package. A volume-averaged approach similar to the method used in UZF-MT3DMS was adopted. The PHREEQC-based computation of chemical processes within PHT3D-UZF in combination with the analytical solution method of UZF1 allows for comprehensive reactive transport investigations (i.e., biogeochemical transformations) that jointly involve saturated and unsaturated zone processes. Intended for regional-scale applications, UZF1 simulates downward-only flux within the unsaturated zone. The model was tested by comparing simulation results with those of existing numerical models. The comparison was performed for several benchmark problems that cover a range of important hydrological and reactive transport processes. A 2D simulation scenario was defined to illustrate the geochemical evolution following dewatering in a sandy acid sulfate soil environment. Other potential applications include the simulation of biogeochemical processes in variably-saturated systems that track the transport and fate of agricultural pollutants, nutrients, natural and xenobiotic organic compounds and micropollutants such as pharmaceuticals, as well as the evolution of isotope patterns.
MR image denoising method for brain surface 3D modeling
NASA Astrophysics Data System (ADS)
Zhao, De-xin; Liu, Peng-jie; Zhang, De-gan
2014-11-01
Three-dimensional (3D) modeling of medical images is a critical part of surgical simulation. In this paper, we focus on the magnetic resonance (MR) images denoising for brain modeling reconstruction, and exploit a practical solution. We attempt to remove the noise existing in the MR imaging signal and preserve the image characteristics. A wavelet-based adaptive curve shrinkage function is presented in spherical coordinates system. The comparative experiments show that the denoising method can preserve better image details and enhance the coefficients of contours. Using these denoised images, the brain 3D visualization is given through surface triangle mesh model, which demonstrates the effectiveness of the proposed method.
A method to fabricate disconnected silver nanostructures in 3D.
Vora, Kevin; Kang, SeungYeon; Mazur, Eric
2012-01-01
The standard nanofabrication toolkit includes techniques primarily aimed at creating 2D patterns in dielectric media. Creating metal patterns on a submicron scale requires a combination of nanofabrication tools and several material processing steps. For example, steps to create planar metal structures using ultraviolet photolithography and electron-beam lithography can include sample exposure, sample development, metal deposition, and metal liftoff. To create 3D metal structures, the sequence is repeated multiple times. The complexity and difficulty of stacking and aligning multiple layers limits practical implementations of 3D metal structuring using standard nanofabrication tools. Femtosecond-laser direct-writing has emerged as a pre-eminent technique for 3D nanofabrication.(1,2) Femtosecond lasers are frequently used to create 3D patterns in polymers and glasses.(3-7) However, 3D metal direct-writing remains a challenge. Here, we describe a method to fabricate silver nanostructures embedded inside a polymer matrix using a femtosecond laser centered at 800 nm. The method enables the fabrication of patterns not feasible using other techniques, such as 3D arrays of disconnected silver voxels.(8) Disconnected 3D metal patterns are useful for metamaterials where unit cells are not in contact with each other,(9) such as coupled metal dot(10,11)or coupled metal rod(12,13) resonators. Potential applications include negative index metamaterials, invisibility cloaks, and perfect lenses. In femtosecond-laser direct-writing, the laser wavelength is chosen such that photons are not linearly absorbed in the target medium. When the laser pulse duration is compressed to the femtosecond time scale and the radiation is tightly focused inside the target, the extremely high intensity induces nonlinear absorption. Multiple photons are absorbed simultaneously to cause electronic transitions that lead to material modification within the focused region. Using this approach, one can
Novel 3D Compression Methods for Geometry, Connectivity and Texture
NASA Astrophysics Data System (ADS)
Siddeq, M. M.; Rodrigues, M. A.
2016-06-01
A large number of applications in medical visualization, games, engineering design, entertainment, heritage, e-commerce and so on require the transmission of 3D models over the Internet or over local networks. 3D data compression is an important requirement for fast data storage, access and transmission within bandwidth limitations. The Wavefront OBJ (object) file format is commonly used to share models due to its clear simple design. Normally each OBJ file contains a large amount of data (e.g. vertices and triangulated faces, normals, texture coordinates and other parameters) describing the mesh surface. In this paper we introduce a new method to compress geometry, connectivity and texture coordinates by a novel Geometry Minimization Algorithm (GM-Algorithm) in connection with arithmetic coding. First, each vertex ( x, y, z) coordinates are encoded to a single value by the GM-Algorithm. Second, triangle faces are encoded by computing the differences between two adjacent vertex locations, which are compressed by arithmetic coding together with texture coordinates. We demonstrate the method on large data sets achieving compression ratios between 87 and 99 % without reduction in the number of reconstructed vertices and triangle faces. The decompression step is based on a Parallel Fast Matching Search Algorithm (Parallel-FMS) to recover the structure of the 3D mesh. A comparative analysis of compression ratios is provided with a number of commonly used 3D file formats such as VRML, OpenCTM and STL highlighting the performance and effectiveness of the proposed method.
Color dithering methods for LEGO-like 3D printing
NASA Astrophysics Data System (ADS)
Sun, Pei-Li; Sie, Yuping
2015-01-01
Color dithering methods for LEGO-like 3D printing are proposed in this study. The first method is work for opaque color brick building. It is a modification of classic error diffusion. Many color primaries can be chosen. However, RGBYKW is recommended as its image quality is good and the number of color primary is limited. For translucent color bricks, multi-layer color building can enhance the image quality significantly. A LUT-based method is proposed to speed the dithering proceeding and make the color distribution even smoother. Simulation results show the proposed multi-layer dithering method can really improve the image quality of LEGO-like 3D printing.
3-D PARTICLE TRANSPORT WITHIN THE HUMAN UPPER RESPIRATORY TRACT
In this study trajectories of inhaled particulate matter (PM) were simulated within a three-dimensional (3-D) computer model of the human upper respiratory tract (URT). The airways were described by computer-reconstructed images of a silicone rubber cast of the human head, throat...
Parallel 3-D method of characteristics in MPACT
Kochunas, B.; Dovvnar, T. J.; Liu, Z.
2013-07-01
A new parallel 3-D MOC kernel has been developed and implemented in MPACT which makes use of the modular ray tracing technique to reduce computational requirements and to facilitate parallel decomposition. The parallel model makes use of both distributed and shared memory parallelism which are implemented with the MPI and OpenMP standards, respectively. The kernel is capable of parallel decomposition of problems in space, angle, and by characteristic rays up to 0(104) processors. Initial verification of the parallel 3-D MOC kernel was performed using the Takeda 3-D transport benchmark problems. The eigenvalues computed by MPACT are within the statistical uncertainty of the benchmark reference and agree well with the averages of other participants. The MPACT k{sub eff} differs from the benchmark results for rodded and un-rodded cases by 11 and -40 pcm, respectively. The calculations were performed for various numbers of processors and parallel decompositions up to 15625 processors; all producing the same result at convergence. The parallel efficiency of the worst case was 60%, while very good efficiency (>95%) was observed for cases using 500 processors. The overall run time for the 500 processor case was 231 seconds and 19 seconds for the case with 15625 processors. Ongoing work is focused on developing theoretical performance models and the implementation of acceleration techniques to minimize the number of iterations to converge. (authors)
NASA Technical Reports Server (NTRS)
Davis, Anthony B.; Marshak, Alexander
2010-01-01
The interplay of sunlight with clouds is a ubiquitous and often pleasant visual experience, but it conjures up major challenges for weather, climate, environmental science and beyond. Those engaged in the characterization of clouds (and the clear air nearby) by remote sensing methods are even more confronted. The problem comes, on the one hand, from the spatial complexity of real clouds and, on the other hand, from the dominance of multiple scattering in the radiation transport. The former ingredient contrasts sharply with the still popular representation of clouds as homogeneous plane-parallel slabs for the purposes of radiative transfer computations. In typical cloud scenes the opposite asymptotic transport regimes of diffusion and ballistic propagation coexist. We survey the three-dimensional (3D) atmospheric radiative transfer literature over the past 50 years and identify three concurrent and intertwining thrusts: first, how to assess the damage (bias) caused by 3D effects in the operational 1D radiative transfer models? Second, how to mitigate this damage? Finally, can we exploit 3D radiative transfer phenomena to innovate observation methods and technologies? We quickly realize that the smallest scale resolved computationally or observationally may be artificial but is nonetheless a key quantity that separates the 3D radiative transfer solutions into two broad and complementary classes: stochastic and deterministic. Both approaches draw on classic and contemporary statistical, mathematical and computational physics.
NASA Astrophysics Data System (ADS)
Davis, Anthony B.; Marshak, Alexander
2010-02-01
The interplay of sunlight with clouds is a ubiquitous and often pleasant visual experience, but it conjures up major challenges for weather, climate, environmental science and beyond. Those engaged in the characterization of clouds (and the clear air nearby) by remote sensing methods are even more confronted. The problem comes, on the one hand, from the spatial complexity of real clouds and, on the other hand, from the dominance of multiple scattering in the radiation transport. The former ingredient contrasts sharply with the still popular representation of clouds as homogeneous plane-parallel slabs for the purposes of radiative transfer computations. In typical cloud scenes the opposite asymptotic transport regimes of diffusion and ballistic propagation coexist. We survey the three-dimensional (3D) atmospheric radiative transfer literature over the past 50 years and identify three concurrent and intertwining thrusts: first, how to assess the damage (bias) caused by 3D effects in the operational 1D radiative transfer models? Second, how to mitigate this damage? Finally, can we exploit 3D radiative transfer phenomena to innovate observation methods and technologies? We quickly realize that the smallest scale resolved computationally or observationally may be artificial but is nonetheless a key quantity that separates the 3D radiative transfer solutions into two broad and complementary classes: stochastic and deterministic. Both approaches draw on classic and contemporary statistical, mathematical and computational physics.
SAMA: A Method for 3D Morphological Analysis
Cerruti, Florent; Sonnenschein, Carlos; Soto, Ana M.
2016-01-01
Three-dimensional (3D) culture models are critical tools for understanding tissue morphogenesis. A key requirement for their analysis is the ability to reconstruct the tissue into computational models that allow quantitative evaluation of the formed structures. Here, we present Software for Automated Morphological Analysis (SAMA), a method by which epithelial structures grown in 3D cultures can be imaged, reconstructed and analyzed with minimum human intervention. SAMA allows quantitative analysis of key features of epithelial morphogenesis such as ductal elongation, branching and lumen formation that distinguish different hormonal treatments. SAMA is a user-friendly set of customized macros operated via FIJI (http://fiji.sc/Fiji), an open-source image analysis platform in combination with a set of functions in R (http://www.r-project.org/), an open-source program for statistical analysis. SAMA enables a rapid, exhaustive and quantitative 3D analysis of the shape of a population of structures in a 3D image. SAMA is cross-platform, licensed under the GPLv3 and available at http://montevil.theobio.org/content/sama. PMID:27035711
A method for building 3D models of barchan dunes
NASA Astrophysics Data System (ADS)
Nai, Yang; Li-lan, Su; Lin, Wan; Jie, Yang; Shi-yi, Chen; Wei-lu, Hu
2016-01-01
The distributions of barchan dunes are usually represented by digital terrain models (DTMs) overlaid with digital orthophoto maps. Given that most regions with barchan dues have low relief, a 3D map obtained from a DTM may ineffectively show the stereoscopic shape of each dune. The method of building 3D models of barchan dunes using existing modeling software seldom considers the geographical environment. As a result, barchan dune models are often inconsistent with actual DTMs and incompletely express the morphological characteristics of dunes. Manual construction of barchan dune models is also costly and time consuming. Considering these problems, the morphological characteristics of barchan dunes and the mathematical relationships between the morphological parameters of the dunes, such as length, height, and width, are analyzed in this study. The methods of extracting the morphological feature points of barchan dunes, calculating their morphological parameters and building dune outlines and skeleton lines based on the medial axes, are also presented. The dune outlines, skeleton lines, and part of the medial axes of dunes are used to construct a constrained triangulated irregular network. C# and ArcEngine are employed to build 3D models of barchan dunes automatically. Experimental results of a study conducted in Tengger Desert show that the method can be used to approximate the morphological characteristics of barchan dunes and is less time consuming than manual methods.
Breast tumour visualization using 3D quantitative ultrasound methods
NASA Astrophysics Data System (ADS)
Gangeh, Mehrdad J.; Raheem, Abdul; Tadayyon, Hadi; Liu, Simon; Hadizad, Farnoosh; Czarnota, Gregory J.
2016-04-01
Breast cancer is one of the most common cancer types accounting for 29% of all cancer cases. Early detection and treatment has a crucial impact on improving the survival of affected patients. Ultrasound (US) is non-ionizing, portable, inexpensive, and real-time imaging modality for screening and quantifying breast cancer. Due to these attractive attributes, the last decade has witnessed many studies on using quantitative ultrasound (QUS) methods in tissue characterization. However, these studies have mainly been limited to 2-D QUS methods using hand-held US (HHUS) scanners. With the availability of automated breast ultrasound (ABUS) technology, this study is the first to develop 3-D QUS methods for the ABUS visualization of breast tumours. Using an ABUS system, unlike the manual 2-D HHUS device, the whole patient's breast was scanned in an automated manner. The acquired frames were subsequently examined and a region of interest (ROI) was selected in each frame where tumour was identified. Standard 2-D QUS methods were used to compute spectral and backscatter coefficient (BSC) parametric maps on the selected ROIs. Next, the computed 2-D parameters were mapped to a Cartesian 3-D space, interpolated, and rendered to provide a transparent color-coded visualization of the entire breast tumour. Such 3-D visualization can potentially be used for further analysis of the breast tumours in terms of their size and extension. Moreover, the 3-D volumetric scans can be used for tissue characterization and the categorization of breast tumours as benign or malignant by quantifying the computed parametric maps over the whole tumour volume.
Optical Sensors and Methods for Underwater 3D Reconstruction
Massot-Campos, Miquel; Oliver-Codina, Gabriel
2015-01-01
This paper presents a survey on optical sensors and methods for 3D reconstruction in underwater environments. The techniques to obtain range data have been listed and explained, together with the different sensor hardware that makes them possible. The literature has been reviewed, and a classification has been proposed for the existing solutions. New developments, commercial solutions and previous reviews in this topic have also been gathered and considered. PMID:26694389
Optical Sensors and Methods for Underwater 3D Reconstruction.
Massot-Campos, Miquel; Oliver-Codina, Gabriel
2015-01-01
This paper presents a survey on optical sensors and methods for 3D reconstruction in underwater environments. The techniques to obtain range data have been listed and explained, together with the different sensor hardware that makes them possible. The literature has been reviewed, and a classification has been proposed for the existing solutions. New developments, commercial solutions and previous reviews in this topic have also been gathered and considered. PMID:26694389
Discrete Method of Images for 3D Radio Propagation Modeling
NASA Astrophysics Data System (ADS)
Novak, Roman
2016-09-01
Discretization by rasterization is introduced into the method of images (MI) in the context of 3D deterministic radio propagation modeling as a way to exploit spatial coherence of electromagnetic propagation for fine-grained parallelism. Traditional algebraic treatment of bounding regions and surfaces is replaced by computer graphics rendering of 3D reflections and double refractions while building the image tree. The visibility of reception points and surfaces is also resolved by shader programs. The proposed rasterization is shown to be of comparable run time to that of the fundamentally parallel shooting and bouncing rays. The rasterization does not affect the signal evaluation backtracking step, thus preserving its advantage over the brute force ray-tracing methods in terms of accuracy. Moreover, the rendering resolution may be scaled back for a given level of scenario detail with only marginal impact on the image tree size. This allows selection of scene optimized execution parameters for faster execution, giving the method a competitive edge. The proposed variant of MI can be run on any GPU that supports real-time 3D graphics.
Parallel 3D Mortar Element Method for Adaptive Nonconforming Meshes
NASA Technical Reports Server (NTRS)
Feng, Huiyu; Mavriplis, Catherine; VanderWijngaart, Rob; Biswas, Rupak
2004-01-01
High order methods are frequently used in computational simulation for their high accuracy. An efficient way to avoid unnecessary computation in smooth regions of the solution is to use adaptive meshes which employ fine grids only in areas where they are needed. Nonconforming spectral elements allow the grid to be flexibly adjusted to satisfy the computational accuracy requirements. The method is suitable for computational simulations of unsteady problems with very disparate length scales or unsteady moving features, such as heat transfer, fluid dynamics or flame combustion. In this work, we select the Mark Element Method (MEM) to handle the non-conforming interfaces between elements. A new technique is introduced to efficiently implement MEM in 3-D nonconforming meshes. By introducing an "intermediate mortar", the proposed method decomposes the projection between 3-D elements and mortars into two steps. In each step, projection matrices derived in 2-D are used. The two-step method avoids explicitly forming/deriving large projection matrices for 3-D meshes, and also helps to simplify the implementation. This new technique can be used for both h- and p-type adaptation. This method is applied to an unsteady 3-D moving heat source problem. With our new MEM implementation, mesh adaptation is able to efficiently refine the grid near the heat source and coarsen the grid once the heat source passes. The savings in computational work resulting from the dynamic mesh adaptation is demonstrated by the reduction of the the number of elements used and CPU time spent. MEM and mesh adaptation, respectively, bring irregularity and dynamics to the computer memory access pattern. Hence, they provide a good way to gauge the performance of computer systems when running scientific applications whose memory access patterns are irregular and unpredictable. We select a 3-D moving heat source problem as the Unstructured Adaptive (UA) grid benchmark, a new component of the NAS Parallel
SOLIDFELIX: a transportable 3D static volume display
NASA Astrophysics Data System (ADS)
Langhans, Knut; Kreft, Alexander; Wörden, Henrik Tom
2009-02-01
Flat 2D screens cannot display complex 3D structures without the usage of different slices of the 3D model. Volumetric displays like the "FELIX 3D-Displays" can solve the problem. They provide space-filling images and are characterized by "multi-viewer" and "all-round view" capabilities without requiring cumbersome goggles. In the past many scientists tried to develop similar 3D displays. Our paper includes an overview from 1912 up to today. During several years of investigations on swept volume displays within the "FELIX 3D-Projekt" we learned about some significant disadvantages of rotating screens, for example hidden zones. For this reason the FELIX-Team started investigations also in the area of static volume displays. Within three years of research on our 3D static volume display at a normal high school in Germany we were able to achieve considerable results despite minor funding resources within this non-commercial group. Core element of our setup is the display volume which consists of a cubic transparent material (crystal, glass, or polymers doped with special ions, mainly from the rare earth group or other fluorescent materials). We focused our investigations on one frequency, two step upconversion (OFTS-UC) and two frequency, two step upconversion (TFTSUC) with IR-Lasers as excitation source. Our main interest was both to find an appropriate material and an appropriate doping for the display volume. Early experiments were carried out with CaF2 and YLiF4 crystals doped with 0.5 mol% Er3+-ions which were excited in order to create a volumetric pixel (voxel). In addition to that the crystals are limited to a very small size which is the reason why we later investigated on heavy metal fluoride glasses which are easier to produce in large sizes. Currently we are using a ZBLAN glass belonging to the mentioned group and making it possible to increase both the display volume and the brightness of the images significantly. Although, our display is currently
NASA Astrophysics Data System (ADS)
Sharkawi, K.-H.; Abdul-Rahman, A.
2013-09-01
to LoD4. The accuracy and structural complexity of the 3D objects increases with the LoD level where LoD0 is the simplest LoD (2.5D; Digital Terrain Model (DTM) + building or roof print) while LoD4 is the most complex LoD (architectural details with interior structures). Semantic information is one of the main components in CityGML and 3D City Models, and provides important information for any analyses. However, more often than not, the semantic information is not available for the 3D city model due to the unstandardized modelling process. One of the examples is where a building is normally generated as one object (without specific feature layers such as Roof, Ground floor, Level 1, Level 2, Block A, Block B, etc). This research attempts to develop a method to improve the semantic data updating process by segmenting the 3D building into simpler parts which will make it easier for the users to select and update the semantic information. The methodology is implemented for 3D buildings in LoD2 where the buildings are generated without architectural details but with distinct roof structures. This paper also introduces hybrid semantic-geometric 3D segmentation method that deals with hierarchical segmentation of a 3D building based on its semantic value and surface characteristics, fitted by one of the predefined primitives. For future work, the segmentation method will be implemented as part of the change detection module that can detect any changes on the 3D buildings, store and retrieve semantic information of the changed structure, automatically updates the 3D models and visualize the results in a userfriendly graphical user interface (GUI).
Towards a 3D Space Radiation Transport Code
NASA Technical Reports Server (NTRS)
Wilson, J. W.; Tripathl, R. K.; Cicomptta, F. A.; Heinbockel, J. H.; Tweed, J.
2002-01-01
High-speed computational procedures for space radiation shielding have relied on asymptotic expansions in terms of the off-axis scatter and replacement of the general geometry problem by a collection of flat plates. This type of solution was derived for application to human rated systems in which the radius of the shielded volume is large compared to the off-axis diffusion limiting leakage at lateral boundaries. Over the decades these computational codes are relatively complete and lateral diffusion effects are now being added. The analysis for developing a practical full 3D space shielding code is presented.
System and method for 3D printing of aerogels
Worsley, Marcus A.; Duoss, Eric; Kuntz, Joshua; Spadaccini, Christopher; Zhu, Cheng
2016-03-08
A method of forming an aerogel. The method may involve providing a graphene oxide powder and mixing the graphene oxide powder with a solution to form an ink. A 3D printing technique may be used to write the ink into a catalytic solution that is contained in a fluid containment member to form a wet part. The wet part may then be cured in a sealed container for a predetermined period of time at a predetermined temperature. The cured wet part may then be dried to form a finished aerogel part.
Method and simulation to study 3D crosstalk perception
NASA Astrophysics Data System (ADS)
Khaustova, Dar'ya; Blondé, Laurent; Huynh-Thu, Quan; Vienne, Cyril; Doyen, Didier
2012-03-01
To various degrees, all modern 3DTV displays suffer from crosstalk, which can lead to a decrease of both visual quality and visual comfort, and also affect perception of depth. In the absence of a perfect 3D display technology, crosstalk has to be taken into account when studying perception of 3D stereoscopic content. In order to improve 3D presentation systems and understand how to efficiently eliminate crosstalk, it is necessary to understand its impact on human perception. In this paper, we present a practical method to study the perception of crosstalk. The approach consists of four steps: (1) physical measurements of a 3DTV, (2) building of a crosstalk surface based on those measurements and representing specifically the behavior of that 3TV, (3) manipulation of the crosstalk function and application on reference images to produce test images degraded by crosstalk in various ways, and (4) psychophysical tests. Our approach allows both a realistic representation of the behavior of a 3DTV and the easy manipulation of its resulting crosstalk in order to conduct psycho-visual experiments. Our approach can be used in all studies requiring the understanding of how crosstalk affects perception of stereoscopic content and how it can be corrected efficiently.
Single-camera fixed perspective 360-deg 3D method
NASA Astrophysics Data System (ADS)
Harding, Kevin G.; Fergan, Robert K.
1997-01-01
The use of 3D methods for such applications as feature locations within a wide field-of-view, such as for automated guided vehicles or large assembly work, offers some distinct challenges. The use of stereo viewing has often been the method of choice due to the wide area coverage and hardware simplicity. However, stereo based methods suffer from a loss of spatial position resolution for more distant object as compared to close objects due to the high demagnification needed to cover large fields-of-view. A long depth-of-field in such systems may also degrade the general ability to perform correlations due to poor focus. In addition, stereo looses distance resolution for features nearing the line of the two cameras, typically requiring movement of the cameras. The paper presents a novel method of obtaining 3D scene information as seen from the center of a cylindrical field. The method described uses a single camera with a view that is rotated through 360 degrees by means of a continuously rotating mirror. The viewing systems uses a constant field of view optical system that provides a constant X-Y resolution of features in the scene over depths of several meters. Comparing successive images with the readout from an encoder on the rotating mirror generates all locations of objects within a limited height cylinder. This paper will discuss the sources of errors and typical capabilities of this approach in light of a real-time part location tracking application useful in assembly systems.
3D reconstruction methods of coronal structures by radio observations
NASA Astrophysics Data System (ADS)
Aschwanden, Markus J.; Bastian, T. S.; White, Stephen M.
1992-11-01
The ability to carry out the three dimensional (3D) reconstruction of structures in the solar corona would represent a major advance in the study of the physical properties in active regions and in flares. Methods which allow a geometric reconstruction of quasistationary coronal structures (for example active region loops) or dynamic structures (for example flaring loops) are described: stereoscopy of multi-day imaging observations by the VLA (Very Large Array); tomography of optically thin emission (in radio or soft x-rays); multifrequency band imaging by the VLA; and tracing of magnetic field lines by propagating electron beams.
3D reconstruction methods of coronal structures by radio observations
NASA Technical Reports Server (NTRS)
Aschwanden, Markus J.; Bastian, T. S.; White, Stephen M.
1992-01-01
The ability to carry out the three dimensional (3D) reconstruction of structures in the solar corona would represent a major advance in the study of the physical properties in active regions and in flares. Methods which allow a geometric reconstruction of quasistationary coronal structures (for example active region loops) or dynamic structures (for example flaring loops) are described: stereoscopy of multi-day imaging observations by the VLA (Very Large Array); tomography of optically thin emission (in radio or soft x-rays); multifrequency band imaging by the VLA; and tracing of magnetic field lines by propagating electron beams.
Transport of iron oxide nanoparticles in saturated porous media: a large-scale 3D study
NASA Astrophysics Data System (ADS)
Velimirovic, Milica; Schmid, Doris; Micić, Vesna; Miyajima, Kumiko; Klaas, Norbert; Braun, Jürgen; Bosch, Julian; Meckenstock, Rainer; von der Kammer, Frank; Hofmann, Thilo
2016-04-01
Iron oxide nanoparticles (FeOxNp) have a high potential as electron acceptor for in situ microbial oxidation of a wide range of recalcitrant groundwater contaminants (Bosch et al., 2010). Tosco et al. (2012) reported on high colloidal stability of FeOxNp dispersed in water, their low deposition behavior, and consequently improved transport in column experiments compared to extensively studied zerovalent iron nanoparticles. However, determination of FeOxNp transport behavior at the field-relevant conditions has not been done before. The present work is aimed to evaluate different complementary methods for detection, quantification and transport characterization of FeOxNp in a large-scale three-dimensional (3D) model aquifer. Prior to that, batch-scale experiments were performed in order to elucidate the potential of the selected methods for direct and indirect characterization and detection of FeOxNp. Direct methods included measurements of particle size distribution, particle concentration, Fetot content and turbidity of the FeOxNp suspension. Indirect methods included measurements of particle zeta potential, as well as TOC content and pH of the FeOxNp suspension. The results of the batch experiments indicated that the most suitable approach for detecting and quantifying FeOxNp was measuring Fetot content and suspension turbidity, as well as particle size determined using dynamic light scattering principle. These complementary methods were further applied in a large-scale 3D study containing medium and coarse sand in order to 1) assess the transport of FeOxNp in saturated porous medium during injection (VFeOx = 6 m3, cparticle = 20 g/L, Qinj = 0.7 m3/h), and 2) illustrate their spatial distribution after injection. The outcomes of the large-scale 3D study confirmed that FeOxNp transport can be successfully investigated applying complementary methods. Monitoring data including Fetot content, turbidity and particle size showed the transport of particles towards the
A perceptual preprocess method for 3D-HEVC
NASA Astrophysics Data System (ADS)
Shi, Yawen; Wang, Yongfang; Wang, Yubing
2015-08-01
A perceptual preprocessing method for 3D-HEVC coding is proposed in the paper. Firstly we proposed a new JND model, which accounts for luminance contrast masking effect, spatial masking effect, and temporal masking effect, saliency characteristic as well as depth information. We utilize spectral residual approach to obtain the saliency map and built a visual saliency factor based on saliency map. In order to distinguish the sensitivity of objects in different depth. We segment each texture frame into foreground and background by a automatic threshold selection algorithm using corresponding depth information, and then built a depth weighting factor. A JND modulation factor is built with a linear combined with visual saliency factor and depth weighting factor to adjust the JND threshold. Then, we applied the proposed JND model to 3D-HEVC for residual filtering and distortion coefficient processing. The filtering process is that the residual value will be set to zero if the JND threshold is greater than residual value, or directly subtract the JND threshold from residual value if JND threshold is less than residual value. Experiment results demonstrate that the proposed method can achieve average bit rate reduction of 15.11%, compared to the original coding scheme with HTM12.1, while maintains the same subjective quality.
3D unstructured-mesh radiation transport codes
Morel, J.
1997-12-31
Three unstructured-mesh radiation transport codes are currently being developed at Los Alamos National Laboratory. The first code is ATTILA, which uses an unstructured tetrahedral mesh in conjunction with standard Sn (discrete-ordinates) angular discretization, standard multigroup energy discretization, and linear-discontinuous spatial differencing. ATTILA solves the standard first-order form of the transport equation using source iteration in conjunction with diffusion-synthetic acceleration of the within-group source iterations. DANTE is designed to run primarily on workstations. The second code is DANTE, which uses a hybrid finite-element mesh consisting of arbitrary combinations of hexahedra, wedges, pyramids, and tetrahedra. DANTE solves several second-order self-adjoint forms of the transport equation including the even-parity equation, the odd-parity equation, and a new equation called the self-adjoint angular flux equation. DANTE also offers three angular discretization options: $S{_}n$ (discrete-ordinates), $P{_}n$ (spherical harmonics), and $SP{_}n$ (simplified spherical harmonics). DANTE is designed to run primarily on massively parallel message-passing machines, such as the ASCI-Blue machines at LANL and LLNL. The third code is PERICLES, which uses the same hybrid finite-element mesh as DANTE, but solves the standard first-order form of the transport equation rather than a second-order self-adjoint form. DANTE uses a standard $S{_}n$ discretization in angle in conjunction with trilinear-discontinuous spatial differencing, and diffusion-synthetic acceleration of the within-group source iterations. PERICLES was initially designed to run on workstations, but a version for massively parallel message-passing machines will be built. The three codes will be described in detail and computational results will be presented.
A killer micro attack on 3D neutron transport
Dorr, M.R.; Ferguson, J.M.
1990-11-01
We describe the deterministic solution of the neutron transport equation and the computation of the effective criticality of three-dimensional assemblies using the BBN TC2000 killer micros. We observe that the performance of our research code PTRAN running on 48 processors of the TC2000 is competitive with the partially vectorizable version running on a single Cray Y/MP processor. This performance scales well with the number of processors on real problems, including those that are not load balanced a priori. To obtain this performance, we explicitly specify and exploit data locality and data dependence using domain decomposition and dynamic job scheduling. 3 refs., 4 figs., 2 tabs.
A method of PSF generation for 3D brightfield deconvolution.
Tadrous, P J
2010-02-01
This paper addresses the problem of 3D deconvolution of through focus widefield microscope datasets (Z-stacks). One of the most difficult stages in brightfield deconvolution is finding the point spread function. A theoretically calculated point spread function (called a 'synthetic PSF' in this paper) requires foreknowledge of many system parameters and still gives only approximate results. A point spread function measured from a sub-resolution bead suffers from low signal-to-noise ratio, compounded in the brightfield setting (by contrast to fluorescence) by absorptive, refractive and dispersal effects. This paper describes a method of point spread function estimation based on measurements of a Z-stack through a thin sample. This Z-stack is deconvolved by an idealized point spread function derived from the same Z-stack to yield a point spread function of high signal-to-noise ratio that is also inherently tailored to the imaging system. The theory is validated by a practical experiment comparing the non-blind 3D deconvolution of the yeast Saccharomyces cerevisiae with the point spread function generated using the method presented in this paper (called the 'extracted PSF') to a synthetic point spread function. Restoration of both high- and low-contrast brightfield structures is achieved with fewer artefacts using the extracted point spread function obtained with this method. Furthermore the deconvolution progresses further (more iterations are allowed before the error function reaches its nadir) with the extracted point spread function compared to the synthetic point spread function indicating that the extracted point spread function is a better fit to the brightfield deconvolution model than the synthetic point spread function. PMID:20096049
A killer micro attack on 3D neutron transport
Dorr, M.R.; Ferguson, J.M.
1990-11-16
In this paper, we describe the deterministic solution of the neutron transport equation and the computation of the effective criticality of three-dimensional assemblies using the BBN TC2000 killer micros. We observe that the performance of our research code PTRAN running on 48 processors of the TC2000 is competitive with the partially vectorizable version running on a single Cray Y/MP processor. This performance scales well with the number of processors on real problems, including those that are not load balanced a priori. To obtain this performance, we explicitly specify and exploit data locality and data dependence using domain decomposition and dynamic job scheduling. From the results obtained here, it appears that, at least for this application, a production machine based on the TC2000 architecture with more powerful processors and a commensurate increase in switch speed could yield a significant gain in our design capability. 2 refs., 5 figs., 2 tabs.
The 3D inelastic analysis methods for hot section components
NASA Technical Reports Server (NTRS)
Dame, L. T.; Mcknight, R. L.
1983-01-01
The objective of this research is to develop an analytical tool capable of economically evaluating the cyclic time dependent plasticity which occurs in hot section engine components in areas of strain concentration resulting from the combination of both mechanical and thermal stresses. The techniques developed must be capable of accommodating large excursions in temperatures with the associated variations in material properties including plasticity and creep. The overall objective of this proposed program is to develop advanced 3-D inelastic structural/stress analysis methods and solution strategies for more accurate and yet more cost effective analysis of combustors, turbine blades, and vanes. The approach will be to develop four different theories, one linear and three higher order with increasing complexities including embedded singularities.
The 3D inelastic analysis methods for hot section components
NASA Technical Reports Server (NTRS)
Mcknight, R. L.; Maffeo, R. J.; Tipton, M. T.; Weber, G.
1992-01-01
A two-year program to develop advanced 3D inelastic structural stress analysis methods and solution strategies for more accurate and cost effective analysis of combustors, turbine blades, and vanes is described. The approach was to develop a matrix of formulation elements and constitutive models. Three constitutive models were developed in conjunction with optimized iterating techniques, accelerators, and convergence criteria within a framework of dynamic time incrementing. Three formulation models were developed: an eight-noded midsurface shell element; a nine-noded midsurface shell element; and a twenty-noded isoparametric solid element. A separate computer program has been developed for each combination of constitutive model-formulation model. Each program provides a functional stand alone capability for performing cyclic nonlinear structural analysis. In addition, the analysis capabilities incorporated into each program can be abstracted in subroutine form for incorporation into other codes or to form new combinations.
On 3D inelastic analysis methods for hot section components
NASA Technical Reports Server (NTRS)
Mcknight, R. L.; Chen, P. C.; Dame, L. T.; Holt, R. V.; Huang, H.; Hartle, M.; Gellin, S.; Allen, D. H.; Haisler, W. E.
1986-01-01
Accomplishments are described for the 2-year program, to develop advanced 3-D inelastic structural stress analysis methods and solution strategies for more accurate and cost effective analysis of combustors, turbine blades and vanes. The approach was to develop a matrix of formulation elements and constitutive models. Three constitutive models were developed in conjunction with optimized iterating techniques, accelerators, and convergence criteria within a framework of dynamic time incrementing. Three formulations models were developed; an eight-noded mid-surface shell element, a nine-noded mid-surface shell element and a twenty-noded isoparametric solid element. A separate computer program was developed for each combination of constitutive model-formulation model. Each program provides a functional stand alone capability for performing cyclic nonlinear structural analysis. In addition, the analysis capabilities incorporated into each program can be abstracted in subroutine form for incorporation into other codes or to form new combinations.
Domain Decomposition PN Solutions to the 3D Transport Benchmark over a Range in Parameter Space
NASA Astrophysics Data System (ADS)
Van Criekingen, S.
2014-06-01
The objectives of this contribution are twofold. First, the Domain Decomposition (DD) method used in the parafish parallel transport solver is re-interpreted as a Generalized Schwarz Splitting as defined by Tang [SIAM J Sci Stat Comput, vol.13 (2), pp. 573-595, 1992]. Second, parafish provides spherical harmonic (i.e., PN) solutions to the NEA benchmark suite for 3D transport methods and codes over a range in parameter space. To the best of the author's knowledge, these are the first spherical harmonic solutions provided for this demanding benchmark suite. They have been obtained using 512 CPU cores of the JuRoPa machine installed at the Jülich Computing Center (Germany).
Lattice Boltzmann Method for 3-D Flows with Curved Boundary
NASA Technical Reports Server (NTRS)
Mei, Renwei; Shyy, Wei; Yu, Dazhi; Luo, Li-Shi
2002-01-01
In this work, we investigate two issues that are important to computational efficiency and reliability in fluid dynamics applications of the lattice, Boltzmann equation (LBE): (1) Computational stability and accuracy of different lattice Boltzmann models and (2) the treatment of the boundary conditions on curved solid boundaries and their 3-D implementations. Three athermal 3-D LBE models (D3QI5, D3Ql9, and D3Q27) are studied and compared in terms of efficiency, accuracy, and robustness. The boundary treatment recently developed by Filippova and Hanel and Met et al. in 2-D is extended to and implemented for 3-D. The convergence, stability, and computational efficiency of the 3-D LBE models with the boundary treatment for curved boundaries were tested in simulations of four 3-D flows: (1) Fully developed flows in a square duct, (2) flow in a 3-D lid-driven cavity, (3) fully developed flows in a circular pipe, and (4) a uniform flow over a sphere. We found that while the fifteen-velocity 3-D (D3Ql5) model is more prone to numerical instability and the D3Q27 is more computationally intensive, the 63Q19 model provides a balance between computational reliability and efficiency. Through numerical simulations, we demonstrated that the boundary treatment for 3-D arbitrary curved geometry has second-order accuracy and possesses satisfactory stability characteristics.
3D Wavelet-Based Filter and Method
Moss, William C.; Haase, Sebastian; Sedat, John W.
2008-08-12
A 3D wavelet-based filter for visualizing and locating structural features of a user-specified linear size in 2D or 3D image data. The only input parameter is a characteristic linear size of the feature of interest, and the filter output contains only those regions that are correlated with the characteristic size, thus denoising the image.
3D Neutron Transport PWR Full-core Calculation with RMC code
NASA Astrophysics Data System (ADS)
Qiu, Yishu; She, Ding; Fan, Xiao; Wang, Kan; Li, Zeguang; Liang, Jingang; Leroyer, Hadrien
2014-06-01
Nowadays, there are more and more interests in the use of Monte Carlo codes to calculate the detailed power density distributions in full-core reactors. With the Inspur TS1000 HPC Server of Tsinghua University, several calculations have been done based on the EDF 3D Neutron Transport PWR Full-core benchmark through large-scale parallelism. To investigate and compare the results of the deterministic method and Monte Carlo method, EDF R&D and Department of Engineering Physics of Tsinghua University are having a collaboration to make code to code verification. So in this paper, two codes are used. One is the code COCAGNE developed by the EDF R&D, a deterministic core code, and the other is the Monte Carlo code RMC developed by Department of Engineering Physics in Tsinghua University. First, the full-core model is described and a 26-group calculation was performed by these two codes using the same 26-group cross-section library provided by EDF R&D. Then the parallel and tally performance of RMC is discussed. RMC employs a novel algorithm which can cut down most of the communications. It can be seen clearly that the speedup ratio almost linearly increases with the nodes. Furthermore the cell-mapping method applied by RMC consumes little time to tally even millions of cells. The results of the codes COCAGNE and RMC are compared in three ways. The results of these two codes agree well with each other. It can be concluded that both COCAGNE and RMC are able to provide 3D-transport solutions associated with detailed power density distributions calculation in PWR full-core reactors. Finally, to investigate how many histories are needed to obtain a given standard deviation for a full 3D solution, the non-symmetrized condensed 2-group fluxes of RMC are discussed.
3D Compressible Melt Transport with Adaptive Mesh Refinement
NASA Astrophysics Data System (ADS)
Dannberg, Juliane; Heister, Timo
2015-04-01
Melt generation and migration have been the subject of numerous investigations, but their typical time and length-scales are vastly different from mantle convection, which makes it difficult to study these processes in a unified framework. The equations that describe coupled Stokes-Darcy flow have been derived a long time ago and they have been successfully implemented and applied in numerical models (Keller et al., 2013). However, modelling magma dynamics poses the challenge of highly non-linear and spatially variable material properties, in particular the viscosity. Applying adaptive mesh refinement to this type of problems is particularly advantageous, as the resolution can be increased in mesh cells where melt is present and viscosity gradients are high, whereas a lower resolution is sufficient in regions without melt. In addition, previous models neglect the compressibility of both the solid and the fluid phase. However, experiments have shown that the melt density change from the depth of melt generation to the surface leads to a volume increase of up to 20%. Considering these volume changes in both phases also ensures self-consistency of models that strive to link melt generation to processes in the deeper mantle, where the compressibility of the solid phase becomes more important. We describe our extension of the finite-element mantle convection code ASPECT (Kronbichler et al., 2012) that allows for solving additional equations describing the behaviour of silicate melt percolating through and interacting with a viscously deforming host rock. We use the original compressible formulation of the McKenzie equations, augmented by an equation for the conservation of energy. This approach includes both melt migration and melt generation with the accompanying latent heat effects. We evaluate the functionality and potential of this method using a series of simple model setups and benchmarks, comparing results of the compressible and incompressible formulation and
Methods for Geometric Data Validation of 3d City Models
NASA Astrophysics Data System (ADS)
Wagner, D.; Alam, N.; Wewetzer, M.; Pries, M.; Coors, V.
2015-12-01
Geometric quality of 3D city models is crucial for data analysis and simulation tasks, which are part of modern applications of the data (e.g. potential heating energy consumption of city quarters, solar potential, etc.). Geometric quality in these contexts is however a different concept as it is for 2D maps. In the latter case, aspects such as positional or temporal accuracy and correctness represent typical quality metrics of the data. They are defined in ISO 19157 and should be mentioned as part of the metadata. 3D data has a far wider range of aspects which influence their quality, plus the idea of quality itself is application dependent. Thus, concepts for definition of quality are needed, including methods to validate these definitions. Quality on this sense means internal validation and detection of inconsistent or wrong geometry according to a predefined set of rules. A useful starting point would be to have correct geometry in accordance with ISO 19107. A valid solid should consist of planar faces which touch their neighbours exclusively in defined corner points and edges. No gaps between them are allowed, and the whole feature must be 2-manifold. In this paper, we present methods to validate common geometric requirements for building geometry. Different checks based on several algorithms have been implemented to validate a set of rules derived from the solid definition mentioned above (e.g. water tightness of the solid or planarity of its polygons), as they were developed for the software tool CityDoctor. The method of each check is specified, with a special focus on the discussion of tolerance values where they are necessary. The checks include polygon level checks to validate the correctness of each polygon, i.e. closeness of the bounding linear ring and planarity. On the solid level, which is only validated if the polygons have passed validation, correct polygon orientation is checked, after self-intersections outside of defined corner points and edges
New method for 3D reconstruction in digital tomosynthesis
NASA Astrophysics Data System (ADS)
Claus, Bernhard E. H.; Eberhard, Jeffrey W.
2002-05-01
Digital tomosynthesis mammography is an advanced x-ray application that can provide detailed 3D information about the imaged breast. We introduce a novel reconstruction method based on simple backprojection, which yields high contrast reconstructions with reduced artifacts at a relatively low computational complexity. The first step in the proposed reconstruction method is a simple backprojection with an order statistics-based operator (e.g., minimum) used for combining the backprojected images into a reconstructed slice. Accordingly, a given pixel value does generally not contribute to all slices. The percentage of slices where a given pixel value does not contribute, as well as the associated reconstructed values, are collected. Using a form of re-projection consistency constraint, one now updates the projection images, and repeats the order statistics backprojection reconstruction step, but now using the enhanced projection images calculated in the first step. In our digital mammography application, this new approach enhances the contrast of structures in the reconstruction, and allows in particular to recover the loss in signal level due to reduced tissue thickness near the skinline, while keeping artifacts to a minimum. We present results obtained with the algorithm for phantom images.
The 3D Shape of the Dendrite by WDT Method
NASA Astrophysics Data System (ADS)
Tang, Chao; Mitobe, Kazutaka; Yoshimura, Noboru
The purpose of this study is use of a three dimension (3D) measuring system that can automatically measure surface condition. We applied the WDT method that is one of the migration acceleration testing methods, to calculate the spatial variation of the electrodes of ion immigration on a glass epoxy printed wiring board. We also investigated the spatial shape and its variation of dendrite after short circuit for the cases of uniform and nonuniform field strength. As a result the phenomenon of immigration peak of separated matter from cathode to anode due to nonuniform was reported.The moving of the peak of the separated matter is supposed to be due to Cu(OH)2's change in accumulation status. Under the nonuniform and uniform situation, the behavior of separated matter will change after occurring short circuit between the electrodes. Therefore in order to avoid the progress of ion immigration, it is necessary to pay attention to the field strength in hardwiring and the curvature so that the field strength of the wiring pattern cannot be very high.
Two new methods for simulating photolithography development in 3D
Helmsen, J.; Colella, P.; Dorr, M.; Puckett, E.G.
1997-01-30
Two methods are presented for simulating the development of photolithographic profiles during the resist dissolution phase. These algorithms are the volume-of-fluid algorithm, and the steady level-set algorithm. They are compared with the ray-trace, cell, and level-set techniques employed in SAMPLE-3D. The volume-of-fluid algorithm employs an Euclidean Grid with volume fractions. At each time step, the surface is reconstructed by computing an approximation of the tangent plane of the surface in each cell that contains a value between 0 and 1. The geometry constructed in this manner is used to determine flow velocity vectors and the flux across each edge. The material is then advanced by a split advection scheme. The steady Level Set algorithm is an extension of the Iterative Level Set algorithm. The steady Level Set algorithm combines Fast Level Set concepts and a technique for finding zero residual solutions to the ( ) function. The etch time for each cell is calculated in a time ordered manner. Use of heap sorting data structures allows the algorithm to execute extremely quickly. Comparisons of the methods have been performed and results shown.
OS3D/GIMRT software for modeling multicomponent-multidimensional reactive transport
CI Steefel; SB Yabusaki
2000-05-17
OS3D/GIMRT is a numerical software package for simulating multicomponent reactive transport in porous media. The package consists of two principal components: (1) the code OS3D (Operator Splitting 3-Dimensional Reactive Transport) which simulates reactive transport by either splitting the reaction and transport steps in time, i.e., the classic time or operator splitting approach, or by iterating sequentially between reactions and transport, and (2) the code GIMRT (Global Implicit Multicomponent Reactive Transport) which treats up to two dimensional reactive transport with a one step or global implicit approach. Although the two codes do not yet have totally identical capabilities, they can be run from the same input file, allowing comparisons to be made between the two approaches in many cases. The advantages and disadvantages of the two approaches are discussed more fully below, but in general OS3D is designed for simulation of transient concentration fronts, particularly under high Peclet number transport conditions, because of its use of a total variation diminishing or TVD transport algorithm. GIMRT is suited for simulating water-rock alteration over long periods of time where the aqueous concentration field is at or close to a quasi-stationary state and the numerical transport errors are less important. Where water-rock interaction occurs over geological periods of time, GIMRT may be preferable to OS3D because of its ability to take larger time steps.
A support-operator method for 3-D rupture dynamics
NASA Astrophysics Data System (ADS)
Ely, Geoffrey P.; Day, Steven M.; Minster, Jean-Bernard
2009-06-01
We present a numerical method to simulate spontaneous shear crack propagation within a heterogeneous, 3-D, viscoelastic medium. Wave motions are computed on a logically rectangular hexahedral mesh, using the generalized finite-difference method of Support Operators (SOM). This approach enables modelling of non-planar surfaces and non-planar fault ruptures. Our implementation, the Support Operator Rupture Dynamics (SORD) code, is highly scalable, enabling large-scale, multiprocessors calculations. The fault surface is modelled by coupled double nodes, where rupture occurs as dictated by the local stress conditions and a frictional failure law. The method successfully performs test problems developed for the Southern California Earthquake Center (SCEC)/U.S. Geological Survey (USGS) dynamic earthquake rupture code validation exercise, showing good agreement with semi-analytical boundary integral method results. We undertake further dynamic rupture tests to quantify numerical errors introduced by shear deformations to the hexahedral mesh. We generate a family of meshes distorted by simple shearing, in the along-strike direction, up to a maximum of 73°. For SCEC/USGS validation problem number 3, grid-induced errors increase with mesh shear angle, with the logarithm of error approximately proportional to angle over the range tested. At 73°, rms misfits are about 10 per cent for peak slip rate, and 0.5 per cent for both rupture time and total slip, indicating that the method (which, up to now, we have applied mainly to near-vertical strike-slip faulting) is also capable of handling geometries appropriate to low-angle surface-rupturing thrust earthquakes. Additionally, we demonstrate non-planar rupture effects, by modifying the test geometry to include, respectively, cylindrical curvature and sharp kinks.
Ash3d: A finite-volume, conservative numerical model for ash transport and tephra deposition
Schwaiger, Hans F.; Denlinger, Roger P.; Mastin, Larry G.
2012-01-01
We develop a transient, 3-D Eulerian model (Ash3d) to predict airborne volcanic ash concentration and tephra deposition during volcanic eruptions. This model simulates downwind advection, turbulent diffusion, and settling of ash injected into the atmosphere by a volcanic eruption column. Ash advection is calculated using time-varying pre-existing wind data and a robust, high-order, finite-volume method. Our routine is mass-conservative and uses the coordinate system of the wind data, either a Cartesian system local to the volcano or a global spherical system for the Earth. Volcanic ash is specified with an arbitrary number of grain sizes, which affects the fall velocity, distribution and duration of transport. Above the source volcano, the vertical mass distribution with elevation is calculated using a Suzuki distribution for a given plume height, eruptive volume, and eruption duration. Multiple eruptions separated in time may be included in a single simulation. We test the model using analytical solutions for transport. Comparisons of the predicted and observed ash distributions for the 18 August 1992 eruption of Mt. Spurr in Alaska demonstrate to the efficacy and efficiency of the routine.
NASA Astrophysics Data System (ADS)
Lague, D.; Brodu, N.; Leroux, J.
2012-12-01
Ground based lidar and photogrammetric techniques are increasingly used to track the evolution of natural surfaces in 3D at an unprecedented resolution and precision. The range of applications encompass many type of natural surfaces with different geometries and roughness characteristics (landslides, cliff erosion, river beds, bank erosion,....). Unravelling surface change in these contexts requires to compare large point clouds in 2D or 3D. The most commonly used method in geomorphology is based on a 2D difference of the gridded point clouds. Yet this is hardly adapted to many 3D natural environments such as rivers (with horizontal beds and vertical banks), while gridding complex rough surfaces is a complex task. On the other hand, tools allowing to perform 3D comparison are scarce and may require to mesh the point clouds which is difficult on rough natural surfaces. Moreover, existing 3D comparison tools do not provide an explicit calculation of confidence intervals that would factor in registration errors, roughness effects and instrument related position uncertainties. To unlock this problem, we developed the first algorithm combining a 3D measurement of surface change directly on point clouds with an estimate of spatially variable confidence intervals (called M3C2). The method has two steps : (1) surface normal estimation and orientation in 3D at a scale consistent with the local roughness ; (2) measurement of mean surface change along the normal direction with explicit calculation of a local confidence interval. Comparison with existing 3D methods based on a closest-point calculation demonstrates the higher precision of the M3C2 method when mm changes needs to be detected. The M3C2 method is also simple to use as it does not require surface meshing or gridding, and is not sensitive to missing data or change in point density. We also present a 3D classification tool (CANUPO) for vegetation removal based on a new geometrical measure: the multi
3-D dynamic rupture simulations by a finite volume method
NASA Astrophysics Data System (ADS)
Benjemaa, M.; Glinsky-Olivier, N.; Cruz-Atienza, V. M.; Virieux, J.
2009-07-01
Dynamic rupture of a 3-D spontaneous crack of arbitrary shape is investigated using a finite volume (FV) approach. The full domain is decomposed in tetrahedra whereas the surface, on which the rupture takes place, is discretized with triangles that are faces of tetrahedra. First of all, the elastodynamic equations are described into a pseudo-conservative form for an easy application of the FV discretization. Explicit boundary conditions are given using criteria based on the conservation of discrete energy through the crack surface. Using a stress-threshold criterion, these conditions specify fluxes through those triangles that have suffered rupture. On these broken surfaces, stress follows a linear slip-weakening law, although other friction laws can be implemented. For The Problem Version 3 of the dynamic-rupture code verification exercise conducted by the SCEC/USGS, numerical solutions on a planar fault exhibit a very high convergence rate and are in good agreement with the reference one provided by a finite difference (FD) technique. For a non-planar fault of parabolic shape, numerical solutions agree satisfactorily well with those obtained with a semi-analytical boundary integral method in terms of shear stress amplitudes, stopping phases arrival times and stress overshoots. Differences between solutions are attributed to the low-order interpolation of the FV approach, whose results are particularly sensitive to the mesh regularity (structured/unstructured). We expect this method, which is well adapted for multiprocessor parallel computing, to be competitive with others for solving large scale dynamic ruptures scenarios of seismic sources in the near future.
Predicting longshore gradients in longshore transport: the CERC formula compared to Delft3D
List, Jeffrey H.; Hanes, Daniel M.; Ruggiero, Peter
2007-01-01
The prediction of longshore transport gradients is critical for forecasting shoreline change. We employ simple test cases consisting of shoreface pits at varying distances from the shoreline to compare the longshore transport gradients predicted by the CERC formula against results derived from the process-based model Delft3D. Results show that while in some cases the two approaches give very similar results, in many cases the results diverge greatly. Although neither approach is validated with field data here, the Delft3D-based transport gradients provide much more consistent predictions of erosional and accretionary zones as the pit location varies across the shoreface.
Pore-scale intermittent velocity structure underpinning anomalous transport through 3-D porous media
NASA Astrophysics Data System (ADS)
Kang, Peter K.; Anna, Pietro; Nunes, Joao P.; Bijeljic, Branko; Blunt, Martin J.; Juanes, Ruben
2014-09-01
We study the nature of non-Fickian particle transport in 3-D porous media by simulating fluid flow in the intricate pore space of real rock. We solve the full Navier-Stokes equations at the same resolution as the 3-D micro-CT (computed tomography) image of the rock sample and simulate particle transport along the streamlines of the velocity field. We find that transport at the pore scale is markedly anomalous: longitudinal spreading is superdiffusive, while transverse spreading is subdiffusive. We demonstrate that this anomalous behavior originates from the intermittent structure of the velocity field at the pore scale, which in turn emanates from the interplay between velocity heterogeneity and velocity correlation. Finally, we propose a continuous time random walk model that honors this intermittent structure at the pore scale and captures the anomalous 3-D transport behavior at the macroscale.
Methods For Electronic 3-D Moving Pictures Without Glasses
NASA Astrophysics Data System (ADS)
Collender, Robert B.
1987-06-01
This paper describes implementation approaches in image acquisition and playback for 3-D computer graphics, 3-D television and 3-D theatre movies without special glasses. Projection lamps, spatial light modulators, CRT's and dynamic scanning are all eliminated by the application of an active image array, all static components and a semi-specular screen. The resulting picture shows horizontal parallax with a wide horizontal view field (up to 360 de-grees) giving a holographic appearance in full color with smooth continuous viewing without speckle. Static component systems are compared with dynamic component systems using both linear and circular arrays. Implementation of computer graphic systems are shown that allow complex shaded color images to extend from the viewer's eyes to infinity. Large screen systems visible by hundreds of people are feasible by the use of low f-stops and high gain screens in projection. Screen geometries and special screen properties are shown. Viewing characteristics offer no restrictions in view-position over the entire view-field and have a "look-around" feature for all the categories of computer graphics, television and movies. Standard video cassettes and optical discs can also interface the system to generate a 3-D window viewable without glasses. A prognosis is given for technology application to 3-D pictures without glasses that replicate the daily viewing experience. Super-position of computer graphics on real-world pictures is shown feasible.
Bailey, Ryan T; Morway, Eric D; Niswonger, Richard G; Gates, Timothy K
2013-01-01
A numerical model was developed that is capable of simulating multispecies reactive solute transport in variably saturated porous media. This model consists of a modified version of the reactive transport model RT3D (Reactive Transport in 3 Dimensions) that is linked to the Unsaturated-Zone Flow (UZF1) package and MODFLOW. Referred to as UZF-RT3D, the model is tested against published analytical benchmarks as well as other published contaminant transport models, including HYDRUS-1D, VS2DT, and SUTRA, and the coupled flow and transport modeling system of CATHY and TRAN3D. Comparisons in one-dimensional, two-dimensional, and three-dimensional variably saturated systems are explored. While several test cases are included to verify the correct implementation of variably saturated transport in UZF-RT3D, other cases are included to demonstrate the usefulness of the code in terms of model run-time and handling the reaction kinetics of multiple interacting species in variably saturated subsurface systems. As UZF1 relies on a kinematic-wave approximation for unsaturated flow that neglects the diffusive terms in Richards equation, UZF-RT3D can be used for large-scale aquifer systems for which the UZF1 formulation is reasonable, that is, capillary-pressure gradients can be neglected and soil parameters can be treated as homogeneous. Decreased model run-time and the ability to include site-specific chemical species and chemical reactions make UZF-RT3D an attractive model for efficient simulation of multispecies reactive transport in variably saturated large-scale subsurface systems. PMID:23131109
Bailey, Ryan T.; Morway, Eric D.; Niswonger, Richard G.; Gates, Timothy K.
2013-01-01
A numerical model was developed that is capable of simulating multispecies reactive solute transport in variably saturated porous media. This model consists of a modified version of the reactive transport model RT3D (Reactive Transport in 3 Dimensions) that is linked to the Unsaturated-Zone Flow (UZF1) package and MODFLOW. Referred to as UZF-RT3D, the model is tested against published analytical benchmarks as well as other published contaminant transport models, including HYDRUS-1D, VS2DT, and SUTRA, and the coupled flow and transport modeling system of CATHY and TRAN3D. Comparisons in one-dimensional, two-dimensional, and three-dimensional variably saturated systems are explored. While several test cases are included to verify the correct implementation of variably saturated transport in UZF-RT3D, other cases are included to demonstrate the usefulness of the code in terms of model run-time and handling the reaction kinetics of multiple interacting species in variably saturated subsurface systems. As UZF1 relies on a kinematic-wave approximation for unsaturated flow that neglects the diffusive terms in Richards equation, UZF-RT3D can be used for large-scale aquifer systems for which the UZF1 formulation is reasonable, that is, capillary-pressure gradients can be neglected and soil parameters can be treated as homogeneous. Decreased model run-time and the ability to include site-specific chemical species and chemical reactions make UZF-RT3D an attractive model for efficient simulation of multispecies reactive transport in variably saturated large-scale subsurface systems.
X3D moving grid methods for semiconductor applications
Kuprat, A.; Cartwright, D.; Gammel, J.T.; George, D.; Kendrick, B.; Kilcrease, D.; Trease, H.; Walker, R.
1997-11-01
The Los Alamos 3D grid toolbox handles grid maintenance chores and provides access to a sophisticated set of optimization algorithms for unstructured grids. The application of these tools to semiconductor problems is illustrated in three examples: grain growth, topographic deposition and electrostatics. These examples demonstrate adaptive smoothing, front tracking, and automatic, adaptive refinement/derefinement.
Object-oriented urban 3D spatial data model organization method
NASA Astrophysics Data System (ADS)
Li, Jing-wen; Li, Wen-qing; Lv, Nan; Su, Tao
2015-12-01
This paper combined the 3d data model with object-oriented organization method, put forward the model of 3d data based on object-oriented method, implemented the city 3d model to quickly build logical semantic expression and model, solved the city 3d spatial information representation problem of the same location with multiple property and the same property with multiple locations, designed the space object structure of point, line, polygon, body for city of 3d spatial database, and provided a new thought and method for the city 3d GIS model and organization management.
Coupling 2-D cylindrical and 3-D x-y-z transport computations
Abu-Shumays, I.K.; Yehnert, C.E.; Pitcairn, T.N.
1998-06-30
This paper describes a new two-dimensional (2-D) cylindrical geometry to three-dimensional (3-D) rectangular x-y-z splice option for multi-dimensional discrete ordinates solutions to the neutron (photon) transport equation. Of particular interest are the simple transformations developed and applied in order to carry out the required spatial and angular interpolations. The spatial interpolations are linear and equivalent to those applied elsewhere. The angular interpolations are based on a high order spherical harmonics representation of the angular flux. Advantages of the current angular interpolations over previous work are discussed. An application to an intricate streaming problem is provided to demonstrate the advantages of the new method for efficient and accurate prediction of particle behavior in complex geometries.
3D Structure of Saharan Dust Transport Towards Europe as Seen by CALIPSO
NASA Astrophysics Data System (ADS)
Marinou, Eleni; Amiridis, Vassilis; Tsekeri, Alexandra; Solomos, Stavros; Kokkalis, Panos; Proestakis, Emmanouil; Kottas, Michael; Binietoglou, Ioannis; Zanis, Prodromos; Kazadzis, Stelios; Wandinger, Ulla; Ansmann, Albert
2016-06-01
We present a 3D multi-year monthly mean climatology of Saharan dust advection over Europe using an area-optimized pure dust CALIPSO product. The product has been developed by applying EARLINET-measured dust lidar ratios and depolarization-based dust discrimination methods and it is shown to have a very good agreement in terms of AOD when compared to AERONET over Europe/North Africa and MODIS over Mediterranean. The processing of such purely observational data reveals the certain seasonal patterns of dust transportation towards Europe and the Atlantic Ocean. The physical and optical properties of the dust layer are identified for several areas near the Saharan sources, over the Mediterranean and over continental Europe.
3D range scan enhancement using image-based methods
NASA Astrophysics Data System (ADS)
Herbort, Steffen; Gerken, Britta; Schugk, Daniel; Wöhler, Christian
2013-10-01
This paper addresses the problem of 3D surface scan refinement, which is desirable due to noise, outliers, and missing measurements being present in the 3D surfaces obtained with a laser scanner. We present a novel algorithm for the fusion of absolute laser scanner depth profiles and photometrically estimated surface normal data, which yields a noise-reduced and highly detailed depth profile with large scale shape robustness. In contrast to other approaches published in the literature, the presented algorithm (1) regards non-Lambertian surfaces, (2) simultaneously computes surface reflectance (i.e. BRDF) parameters required for 3D reconstruction, (3) models pixelwise incident light and viewing directions, and (4) accounts for interreflections. The algorithm as such relies on the minimization of a three-component error term, which penalizes intensity deviations, integrability deviations, and deviations from the known large-scale surface shape. The solution of the error minimization is obtained iteratively based on a calculus of variations. BRDF parameters are estimated by initially reducing and then iteratively refining the optical resolution, which provides the required robust data basis. The 3D reconstruction of concave surface regions affected by interreflections is improved by compensating global illumination in the image data. The algorithm is evaluated based on eight objects with varying albedos and reflectance behaviors (diffuse, specular, metallic). The qualitative evaluation shows a removal of outliers and a strong reduction of noise, while the large scale shape is preserved. Fine surface details Which are previously not contained in the surface scans, are incorporated through using image data. The algorithm is evaluated with respect to its absolute accuracy using two caliper objects of known shape, and based on synthetically generated data. The beneficial effect of interreflection compensation on the reconstruction accuracy is evaluated quantitatively in a
Calibration Methods for a 3D Triangulation Based Camera
NASA Astrophysics Data System (ADS)
Schulz, Ulrike; Böhnke, Kay
A sensor in a camera takes a gray level image (1536 x 512 pixels), which is reflected by a reference body. The reference body is illuminated by a linear laser line. This gray level image can be used for a 3D calibration. The following paper describes how a calibration program calculates the calibration factors. The calibration factors serve to determine the size of an unknown reference body.
An automated 3D reconstruction method of UAV images
NASA Astrophysics Data System (ADS)
Liu, Jun; Wang, He; Liu, Xiaoyang; Li, Feng; Sun, Guangtong; Song, Ping
2015-10-01
In this paper a novel fully automated 3D reconstruction approach based on low-altitude unmanned aerial vehicle system (UAVs) images will be presented, which does not require previous camera calibration or any other external prior knowledge. Dense 3D point clouds are generated by integrating orderly feature extraction, image matching, structure from motion (SfM) and multi-view stereo (MVS) algorithms, overcoming many of the cost, time limitations of rigorous photogrammetry techniques. An image topology analysis strategy is introduced to speed up large scene reconstruction by taking advantage of the flight-control data acquired by UAV. Image topology map can significantly reduce the running time of feature matching by limiting the combination of images. A high-resolution digital surface model of the study area is produced base on UAV point clouds by constructing the triangular irregular network. Experimental results show that the proposed approach is robust and feasible for automatic 3D reconstruction of low-altitude UAV images, and has great potential for the acquisition of spatial information at large scales mapping, especially suitable for rapid response and precise modelling in disaster emergency.
NASA Astrophysics Data System (ADS)
Hostache, Renaud; Krein, Andreas; Barrière, Julien
2014-05-01
Coupling the 3D hydro-morphodynamic model Telemac-3D-sisyphe and seismic measurements to estimate bedload transport rates in a small gravel-bed river. Renaud Hostache, Andreas Krein, Julien Barrière During flood events, amounts of river bed material are transported via bedload. This causes problems, like the silting of reservoirs or the disturbance of biological habitats. Some current bedload measuring techniques have limited possibilities for studies in high temporal resolutions. Optical systems are usually not applicable because of high turbidity due to concentrated suspended sediment transported. Sediment traps or bedload samplers yield only summative information on bedload transport with low temporal resolution. An alternative bedload measuring technique is the use of seismological systems installed next to the rivers. The potential advantages are observations in real time and under undisturbed conditions. The study area is a 120 m long reach of River Colpach (21.5 km2), a small gravel bed river in Northern Luxembourg. A combined approach of hydro-climatological observations, hydraulic measurements, sediment sampling, and seismological measurements is used in order to investigate bedload transport phenomena. Information derived from seismic measurements and results from a 3-dimensional hydro-morphodynamic model are exemplarily discussed for a November 2013 flood event. The 3-dimensional hydro-morphodynamic model is based on the Telemac hydroinformatic system. This allows for dynamically coupling a 3D hydrodynamic model (Telemac-3D) and a morphodynamic model (Sisyphe). The coupling is dynamic as these models exchange their information during simulations. This is a main advantage as it allows for taking into account the effects of the morphologic changes of the riverbed on the water hydrodynamic and the bedload processes. The coupled model has been calibrated using time series of gauged water depths and time series of bed material collected sequentially (after
A method of multi-view intraoral 3D measurement
NASA Astrophysics Data System (ADS)
Zhao, Huijie; Wang, Zhen; Jiang, Hongzhi; Xu, Yang; Lv, Peijun; Sun, Yunchun
2015-02-01
In dental restoration, its important to achieve a high-accuracy digital impression. Most of the existing intraoral measurement systems can only measure the tooth from a single view. Therfore - if we are wilng to acquire the whole data of a tooth, the scans of the tooth from multi-direction ad the data stitching based on the features of the surface are needed, which increases the measurement duration and influence the measurement accuracy. In this paper, we introduce a fringe-projection based on multi-view intraoral measurement system. It can acquire 3D data of the occlusal surface, the buccal surface and the lingual surface of a tooth synchronously, by using a senor with three mirrors, which aim at the three surfaces respectively and thus expand the measuring area. The constant relationship of the three mirrors is calibrated before measurement and can help stitch the data clouds acquired through different mirrors accurately. Therefore the system can obtain the 3D data of a tooth without the need to measure it from different directions for many times. Experiments proved the availability and reliability of this miniaturized measurement system.
MOM3D method of moments code theory manual
NASA Astrophysics Data System (ADS)
Shaeffer, John F.
1992-03-01
MOM3D is a FORTRAN algorithm that solves Maxwell's equations as expressed via the electric field integral equation for the electromagnetic response of open or closed three dimensional surfaces modeled with triangle patches. Two joined triangles (couples) form the vector current unknowns for the surface. Boundary conditions are for perfectly conducting or resistive surfaces. The impedance matrix represents the fundamental electromagnetic interaction of the body with itself. A variety of electromagnetic analysis options are possible once the impedance matrix is computed including backscatter radar cross section (RCS), bistatic RCS, antenna pattern prediction for user specified body voltage excitation ports, RCS image projection showing RCS scattering center locations, surface currents excited on the body as induced by specified plane wave excitation, and near field computation for the electric field on or near the body.
3D sensitivity of 6-electrode Focused Impedance Method (FIM)
NASA Astrophysics Data System (ADS)
Masum Iquebal, A. H.; Siddique-e Rabbani, K.
2010-04-01
The present work was taken up to have an understanding of the depth sensitivity of the 6 electrode FIM developed by our laboratory earlier, so that it may be applied judiciously for the measurement of organs in 3D, with electrodes on the skin surface. For a fixed electrode geometry sensitivity is expected to depend on the depth, size and conductivity of the target object. With current electrodes 18 cm apart and potential electrodes 5 cm apart, depth sensitivity of spherical conductors, insulators and of pieces of potato of different diameters were measured. The sensitivity dropped sharply with depth gradually leveling off to background, and objects could be sensed down to a depth of about twice their diameters. The sensitivity at a certain depth increases almost linearly with volume for objects with the same conductivity. Thus these results increase confidence in the use of FIM for studying organs at depths of the body.
MOM3D method of moments code theory manual
NASA Technical Reports Server (NTRS)
Shaeffer, John F.
1992-01-01
MOM3D is a FORTRAN algorithm that solves Maxwell's equations as expressed via the electric field integral equation for the electromagnetic response of open or closed three dimensional surfaces modeled with triangle patches. Two joined triangles (couples) form the vector current unknowns for the surface. Boundary conditions are for perfectly conducting or resistive surfaces. The impedance matrix represents the fundamental electromagnetic interaction of the body with itself. A variety of electromagnetic analysis options are possible once the impedance matrix is computed including backscatter radar cross section (RCS), bistatic RCS, antenna pattern prediction for user specified body voltage excitation ports, RCS image projection showing RCS scattering center locations, surface currents excited on the body as induced by specified plane wave excitation, and near field computation for the electric field on or near the body.
NASA Astrophysics Data System (ADS)
Wang, Brian; Goldstein, Moshe; Xu, X. George; Sahoo, Narayan
2005-03-01
Recently, the theoretical framework of the adjoint Monte Carlo (AMC) method has been developed using a simplified patient geometry. In this study, we extended our previous work by applying the AMC framework to a 3D anatomical model called VIP-Man constructed from the Visible Human images. First, the adjoint fluxes for the prostate (PTV) and rectum and bladder (organs at risk (OARs)) were calculated on a spherical surface of 1 m radius, centred at the centre of gravity of PTV. An importance ratio, defined as the PTV dose divided by the weighted OAR doses, was calculated for each of the available beamlets to select the beam angles. Finally, the detailed doses in PTV and OAR were calculated using a forward Monte Carlo simulation to include the electron transport. The dose information was then used to generate dose volume histograms (DVHs). The Pinnacle treatment planning system was also used to generate DVHs for the 3D plans with beam angles obtained from the AMC (3D-AMC) and a standard six-field conformal radiation therapy plan (3D-CRT). Results show that the DVHs for prostate from 3D-AMC and the standard 3D-CRT are very similar, showing that both methods can deliver prescribed dose to the PTV. A substantial improvement in the DVHs for bladder and rectum was found for the 3D-AMC method in comparison to those obtained from 3D-CRT. However, the 3D-AMC plan is less conformal than the 3D-CRT plan because only bladder, rectum and PTV are considered for calculating the importance ratios. Nevertheless, this study clearly demonstrated the feasibility of the AMC in selecting the beam directions as a part of a treatment planning based on the anatomical information in a 3D and realistic patient anatomy.
Reactor Dosimetry Applications Using RAPTOR-M3G:. a New Parallel 3-D Radiation Transport Code
NASA Astrophysics Data System (ADS)
Longoni, Gianluca; Anderson, Stanwood L.
2009-08-01
The numerical solution of the Linearized Boltzmann Equation (LBE) via the Discrete Ordinates method (SN) requires extensive computational resources for large 3-D neutron and gamma transport applications due to the concurrent discretization of the angular, spatial, and energy domains. This paper will discuss the development RAPTOR-M3G (RApid Parallel Transport Of Radiation - Multiple 3D Geometries), a new 3-D parallel radiation transport code, and its application to the calculation of ex-vessel neutron dosimetry responses in the cavity of a commercial 2-loop Pressurized Water Reactor (PWR). RAPTOR-M3G is based domain decomposition algorithms, where the spatial and angular domains are allocated and processed on multi-processor computer architectures. As compared to traditional single-processor applications, this approach reduces the computational load as well as the memory requirement per processor, yielding an efficient solution methodology for large 3-D problems. Measured neutron dosimetry responses in the reactor cavity air gap will be compared to the RAPTOR-M3G predictions. This paper is organized as follows: Section 1 discusses the RAPTOR-M3G methodology; Section 2 describes the 2-loop PWR model and the numerical results obtained. Section 3 addresses the parallel performance of the code, and Section 4 concludes this paper with final remarks and future work.
NASA Astrophysics Data System (ADS)
Druzgalski, Clara; Mani, Ali
2014-11-01
We have investigated the transport dynamics of an electrokinetic instability that occurs when ions are driven from bulk fluids to ion-selective membranes due to externally applied electric fields. This phenomenon is relevant to a wide range of electrochemical applications including electrodialysis for fresh water production. Using data from our 3D DNS, we show how electroconvective instability, arising from concentration polarization, results in a chaotic flow that significantly alters the net ion transport rate across the membrane surface. The 3D DNS results, which fully resolve the spatiotemporal scales including the electric double layers, enable visualization of instantaneous snapshots of current density directly on the membrane surface, as well as analysis of transport statistics such as concentration variance and fluctuating advective fluxes. Furthermore, we present a full spectral analysis revealing broadband spectra in both concentration and flow fields and deduce the key parameter controlling the range of contributing scales.
NASA Astrophysics Data System (ADS)
da Câmara Santa Clara Gomes, Tristan; De La Torre Medina, Joaquín; Velázquez-Galván, Yenni G.; Martínez-Huerta, Juan Manuel; Encinas, Armando; Piraux, Luc
2016-07-01
We have explored the interplay between the magnetic and magneto-transport properties of 3D interconnected nanowire networks made of various magnetic metals by electrodeposition into nanoporous membranes with crossed channels and controlled topology. The close relationship between their magnetic and structural properties has a direct impact on their magneto-transport behavior. In order to accurately and reliably describe the effective magnetic anisotropy and anisotropic magnetoresistance, an analytical model inherent to the topology of 3D nanowire networks is proposed and validated. The feasibility to obtain magneto-transport responses in nanowire network films based on interconnected nanowires makes them very attractive for the development of mechanically stable superstructures that are suitable for potential technological applications.
3D In Vitro Model for Breast Cancer Research Using Magnetic Levitation and Bioprinting Method.
Leonard, Fransisca; Godin, Biana
2016-01-01
Tumor microenvironment composition and architecture are known as a major factor in orchestrating the tumor growth and its response to various therapies. In this context, in vivo studies are necessary to evaluate the responses. However, while tumor cells can be of human origin, tumor microenvironment in the in vivo models is host-based. On the other hand, in vitro studies in a flat monoculture of tumor cells (the most frequently used in vitro tumor model) are unable to recapitulate the complexity of tumor microenvironment. Three-dimensional (3D) in vitro cell cultures of tumor cells have been proven to be an important experimental tool in understanding mechanisms of tumor growth, response to therapeutics, and transport of nutrients/drugs. We have recently described a novel tool to create 3D co-cultures of tumor cells and cells in the tumor microenvironment. Our method utilizes magnetic manipulation/levitation of the specific ratios of tumor cells and cells in the tumor microenvironment (from human or animal origin) aiding in the formation of tumor spheres with defined cellular composition and density, as quickly as within 24 h. This chapter describes the experimental protocols developed to model the 3D structure of the cancer environment using the above method. PMID:26820961
Development of 3-D Ice Accretion Measurement Method
NASA Technical Reports Server (NTRS)
Lee, Sam; Broeren, Andy P.; Addy, Harold E., Jr.; Sills, Robert; Pifer, Ellen M.
2012-01-01
A research plan is currently being implemented by NASA to develop and validate the use of a commercial laser scanner to record and archive fully three-dimensional (3-D) ice shapes from an icing wind tunnel. The plan focused specifically upon measuring ice accreted in the NASA Icing Research Tunnel (IRT). The plan was divided into two phases. The first phase was the identification and selection of the laser scanning system and the post-processing software to purchase and develop further. The second phase was the implementation and validation of the selected system through a series of icing and aerodynamic tests. Phase I of the research plan has been completed. It consisted of evaluating several scanning hardware and software systems against an established selection criteria through demonstrations in the IRT. The results of Phase I showed that all of the scanning systems that were evaluated were equally capable of scanning ice shapes. The factors that differentiated the scanners were ease of use and the ability to operate in a wide range of IRT environmental conditions.
Edge Transport Modeling using the 3D EMC3-Eirene code on Tokamaks and Stellarators
NASA Astrophysics Data System (ADS)
Lore, J. D.; Ahn, J. W.; Briesemeister, A.; Ferraro, N.; Labombard, B.; McLean, A.; Reinke, M.; Shafer, M.; Terry, J.
2015-11-01
The fluid plasma edge transport code EMC3-Eirene has been applied to aid data interpretation and understanding the results of experiments with 3D effects on several tokamaks. These include applied and intrinsic 3D magnetic fields, 3D plasma facing components, and toroidally and poloidally localized heat and particle sources. On Alcator C-Mod, a series of experiments explored the impact of toroidally and poloidally localized impurity gas injection on core confinement and asymmetries in the divertor fluxes, with the differences between the asymmetry in L-mode and H-mode qualitatively reproduced in the simulations due to changes in the impurity ionization in the private flux region. Modeling of NSTX experiments on the effect of 3D fields on detachment matched the trend of a higher density at which the detachment occurs when 3D fields are applied. On DIII-D, different magnetic field models were used in the simulation and compared against the 2D Thomson scattering diagnostic. In simulating each device different aspects of the code model are tested pointing to areas where the model must be further developed. The application to stellarator experiments will also be discussed. Work supported by U.S. DOE: DE-AC05-00OR22725, DE AC02-09CH11466, DE-FC02-99ER54512, and DE-FC02-04ER54698.
Ozone formation during an episode over Europe: A 3-D chemical/transport model simulation
NASA Technical Reports Server (NTRS)
Berntsen, Terje; Isaksen, Ivar S. A.
1994-01-01
A 3-D regional photochemical tracer/transport model for Europe and the Eastern Atlantic has been developed based on the NASA/GISS CTM. The model resolution is 4x5 degrees latitude and longitude with 9 layers in the vertical (7 in the troposphere). Advective winds, convection statistics and other meteorological data from the NASA/GISS GCM are used. An extensive gas-phase chemical scheme based on the scheme used in our global 2D model has been incorporated in the 3D model. In this work ozone formation in the troposphere is studied with the 3D model during a 5 day period starting June 30. Extensive local ozone production is found and the relationship between the source regions and the downwind areas are discussed. Variations in local ozone formation as a function of total emission rate, as well as the composition of the emissions (HC/NO(x)) ratio and isoprene emissions) are elucidated. An important vertical transport process in the troposphere is by convective clouds. The 3D model includes an explicit parameterization of this process. It is shown that this process has significant influence on the calculated surface ozone concentrations.
Advanced quadratures and periodic boundary conditions in parallel 3D S{sub n} transport
Manalo, K.; Yi, C.; Huang, M.; Sjoden, G.
2013-07-01
Significant updates in numerical quadratures have warranted investigation with 3D Sn discrete ordinates transport. We show new applications of quadrature departing from level symmetric (S{sub 2}o). investigating 3 recently developed quadratures: Even-Odd (EO), Linear-Discontinuous Finite Element - Surface Area (LDFE-SA), and the non-symmetric Icosahedral Quadrature (IC). We discuss implementation changes to 3D Sn codes (applied to Hybrid MOC-Sn TITAN and 3D parallel PENTRAN) that can be performed to accommodate Icosahedral Quadrature, as this quadrature is not 90-degree rotation invariant. In particular, as demonstrated using PENTRAN, the properties of Icosahedral Quadrature are suitable for trivial application using periodic BCs versus that of reflective BCs. In addition to implementing periodic BCs for 3D Sn PENTRAN, we implemented a technique termed 'angular re-sweep' which properly conditions periodic BCs for outer eigenvalue iterative loop convergence. As demonstrated by two simple transport problems (3-group fixed source and 3-group reflected/periodic eigenvalue pin cell), we remark that all of the quadratures we investigated are generally superior to level symmetric quadrature, with Icosahedral Quadrature performing the most efficiently for problems tested. (authors)
Optically directed molecular transport and 3D isoelectric positioning of amphoteric biomolecules
Hafeman, Dean G.; Harkins, James B.; WitkowskiII, Charles E.; Lewis, Nathan S.; Brown, Gilbert M; Warmack, Robert J Bruce; Thundat, Thomas George
2006-01-01
We demonstrate the formation of charged molecular packets and their transport within optically created electrical force-field traps in a pH-buffered electrolyte. We call this process photoelectrophoretic localization and transport (PELT). The electrolyte is in contact with a photoconductive semiconductor electrode and a counterelectrode that are connected through an external circuit. A light beam directed to coordinates on the photoconductive electrode surface produces a photocurrent within the circuit and electrolyte. Within the electrolyte, the photocurrent creates localized force-field traps centered at the illuminated coordinates. Charged molecules, including polypeptides and proteins, electrophoretically accumulate into the traps and subsequently can be transported in the electrolyte by moving the traps over the photoconductive electrode in response to movement of the light beam. The molecules in a single trap can be divided into aliquots, and the aliquots can be directed along multiple routes simultaneously by using multiple light beams. This photoelectrophoretic transport of charged molecules by PELT resembles the electrostatic transport of electrons within force-field wells of solid-state charge-coupled devices. The molecules, however, travel in a liquid electrolyte rather than a solid. Furthermore, we have used PELT to position amphoteric biomolecules in three dimensions. A 3D pH gradient was created in an electrolyte medium by controlling the illumination position on a photoconductive anode where protons were generated electrolytically. Photoelectrophoretic transport of amphoteric molecules through the pH gradient resulted in accumulation of the molecules at their apparent 3D isoelectric coordinates in the medium.
2D/1D approximations to the 3D neutron transport equation. II: Numerical comparisons
Kelley, B. W.; Collins, B.; Larsen, E. W.
2013-07-01
In a companion paper [1], (i) several new '2D/1D equations' are introduced as accurate approximations to the 3D Boltzmann transport equation, (ii) the simplest of these approximate equations is systematically discretized, and (iii) a theoretically stable iteration scheme is developed to solve the discrete equations. In this paper, numerical results are presented that confirm the theoretical predictions made in [1]. (authors)
3-D Particl-in-Cell Simulations of Transport Driven Currents
NASA Astrophysics Data System (ADS)
Tsung, F. S.; Dawson, J. M.
1997-11-01
In the advanced tokamak regime, transport phenomena can account for a signficant fraction of the toroidal current, possibly over that driven directly by the ohmic heating electric fields. Although bootstrap theory accounts for contributions of the collisional modification of banana orbits on the toroidal currents, the corresponding transport theory does not accurately predict the transport of particles and heat in present-day tokamak experiments. Furthermore, in our previous simulations in 21/2-D, currents were spontaneously generated in both the cylindrical and the toroidal geometries, contrary to neoclassical predictions. In these calculations, it was believed that the driving mechanism is the preferential loss of particles whose initial velocity is opposite to that of the plasma current. Because the preferential loss mechanism assumes the conservation of toroidal angular momentum, we have extended these simulations to three dimensions to study the effects of toroidal assymetries. A parallel, 3-D electromagnetic PIC code running on the IBM SP, with a localized field-solver has been developed to investigate the effects of perturbations parallel to the field lines, and direct comparisons has been made between the 21/2-D and 3-D simulations, and we have found good agreements between the 2 1/2-D calculations and the 3-D results. We will present these results at the meeting.
Impact of 3D root uptake on solute transport: a numerical study
NASA Astrophysics Data System (ADS)
Schröder, N.; Javaux, M.; Vanderborght, J.; Steffen, B.; Vereecken, H.
2011-12-01
Plant transpiration is an important component of the hydrological cycle. Through root water uptake, plants do not only affect the 3D soil water flow velocity distribution, but also solute movement in soil. This numerical study aims at investigating how solute fate is impacted by root uptake using the 3D biophysical model R-SWMS (Javaux et al., 2008). This model solves the Richards equation in 3D in the soil and the flow equation within the plant root xylem vessels. Furthermore, for solute transport simulations, the 3D particle tracker PARTRACE (Bechtold et al., 2011) was used. . We generated 3D virtual steady-state breakthrough curves (BTC) experiments in soils with transpiring plants. The averaged BTCs were then fitted with a 1D numerical flow model under steady-state conditions to obtain apparent CDE parameters. Two types of root architecture, a fibrous and a taprooted structure, were compared in virtual 3D experiments. The solute uptake type or the transpiration rate were also modified and we analyzed how these parameters affected apparent disperisivity and velocity profiles. Our simulation results show, that both, apparent velocity and dispersivity length are affected by water and solute root uptake. In addition, under high exclusion processes (slight or no active uptake), solute accumulates around roots and generates a long tailing to the breakthrough curves, which cannot be reproduced by 1D models that simulate root water uptake with solute exclusion. This observation may have an important impact on how to model pollutant mass transfer to groundwater at larger scales. Javaux, M., T. Schröder, J. Vanderborght, and H. Vereecken. 2008. Use of a three-dimensional detailed modeling approach for predicting root water uptake. Vadose Zone J. 7:1079-1088.doi: 10.2136/vzj2007.0115. Bechtold, M., S. Haber-Pohlmeier, J. Vanderborght, A. Pohlmeier, P.A. Ferre, and H. Vereecken. 2011. Near-surface solute redistribution during evaporation. Submitted to Geophys. Res. Lett
Teo, B G; Sarinder, K K S; Lim, L H S
2010-08-01
Three-dimensional (3D) models of the marginal hooks, dorsal and ventral anchors, bars and haptoral reservoirs of a parasite, Sundatrema langkawiense Lim & Gibson, 2009 (Monogenea) were developed using the polygonal modelling method in Autodesk 3ds Max (Version 9) based on two-dimensional (2D) illustrations. Maxscripts were written to rotate the modelled 3D structures. Appropriately orientated 3D haptoral hard-parts were then selected and positioned within the transparent 3D outline of the haptor and grouped together to form a complete 3D haptoral entity. This technique is an inexpensive tool for constructing 3D models from 2D illustrations for 3D visualisation of the spatial relationships between the different structural parts within organisms. PMID:20962723
A parametric study of mucociliary transport by numerical simulations of 3D non-homogeneous mucus.
Chatelin, Robin; Poncet, Philippe
2016-06-14
Mucociliary clearance is the natural flow of the mucus which covers and protects the lung from the outer world. Pathologies, like cystic fibrosis, highly change the biological parameters of the mucus flow leading to stagnation situations and pathogens proliferation. As the lung exhibits a complex dyadic structure, in-vivo experimental study of mucociliary clearance is almost impossible and numerical simulations can bring important knowledge about this biological flow. This paper brings a detailed study of the biological parameters influence on the mucociliary clearance, in particular for pathological situations such as cystic fibrosis. Using recent suitable numerical methods, a non-homogeneous mucus flow (including non-linearities) can be simulated efficiently in 3D, allowing the identification of the meaningful parameters involved in this biological flow. Among these parameters, it is shown that the mucus viscosity, the stiffness transition between pericilliary fluid and mucus, the pericilliary fluid height as well as both cilia length and beating frequency have a great influence on the mucociliary transport. PMID:27126985
Method for 3D fibre reconstruction on a microrobotic platform.
Hirvonen, J; Myllys, M; Kallio, P
2016-07-01
Automated handling of a natural fibrous object requires a method for acquiring the three-dimensional geometry of the object, because its dimensions cannot be known beforehand. This paper presents a method for calculating the three-dimensional reconstruction of a paper fibre on a microrobotic platform that contains two microscope cameras. The method is based on detecting curvature changes in the fibre centreline, and using them as the corresponding points between the different views of the images. We test the developed method with four fibre samples and compare the results with the references measured with an X-ray microtomography device. We rotate the samples through 16 different orientations on the platform and calculate the three-dimensional reconstruction to test the repeatability of the algorithm and its sensitivity to the orientation of the sample. We also test the noise sensitivity of the algorithm, and record the mismatch rate of the correspondences provided. We use the iterative closest point algorithm to align the measured three-dimensional reconstructions with the references. The average point-to-point distances between the reconstructed fibre centrelines and the references are 20-30 μm, and the mismatch rate is low. Given the manipulation tolerance, this shows that the method is well suited to automated fibre grasping. This has also been demonstrated with actual grasping experiments. PMID:26695385
New 3-D flow interpolation method on moving ADCP data
NASA Astrophysics Data System (ADS)
Tsubaki, R.; Kawahara, Y.; Muto, Y.; Fujita, I.
2012-05-01
A simple but accurate interpolation procedure for obtaining the three-dimensional distribution of three-component velocity data, from moving acoustic doppler current profiler (ADCP) observation data, is proposed. For understanding actual flow structure within a river with complex bathymetry, the three-dimensional mean velocity field provides a basic picture of the flow. For obtaining the three-dimensional distribution of three-component velocity data, in this work, anisotropic gridding was introduced in order to remove the random component of measured velocity data caused by the turbulence of the flow and measurement error. A continuity correction based on the pressure equation was used to reduce both random and systematic errors. The accuracy of the developed method was evaluated using three-dimensional flow simulation data from a detached-eddy simulation (DES). By using the procedure developed, the complex flow structure surrounding the spur dikes section in the Uji River was successfully visualized and explored. The proposed method shows superiorities in both accuracy and consistency for the interpolated velocity field, as compared to the kriging and inverse-distance weighted (IDW) methods.
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
Filtering method for 3D laser scanning point cloud
NASA Astrophysics Data System (ADS)
Liu, Da; Wang, Li; Hao, Yuncai; Zhang, Jun
2015-10-01
In recent years, with the rapid development of the hardware and software of the three-dimensional model acquisition, three-dimensional laser scanning technology is utilized in various aspects, especially in space exploration. The point cloud filter is very important before using the data. In the paper, considering both the processing quality and computing speed, an improved mean-shift point cloud filter method is proposed. Firstly, by analyze the relevance of the normal vector between the upcoming processing point and the near points, the iterative neighborhood of the mean-shift is selected dynamically, then the high frequency noise is constrained. Secondly, considering the normal vector of the processing point, the normal vector is updated. Finally, updated position is calculated for each point, then each point is moved in the normal vector according to the updated position. The experimental results show that the large features are retained, at the same time, the small sharp features are also existed for different size and shape of objects, so the target feature information is protected precisely. The computational complexity of the proposed method is not high, it can bring high precision results with fast speed, so it is very suitable for space application. It can also be utilized in civil, such as large object measurement, industrial measurement, car navigation etc. In the future, filter with the help of point strength will be further exploited.
Jacob's Interpretation Method Revisited: Accounting for 3-D Spatial Heterogeneity
NASA Astrophysics Data System (ADS)
Sanchez-Vila, X.; Riva, M.; Guadagnini, A.; Carrera, J.
2005-12-01
Traditional approaches to hydraulic test interpretation provide typically individual aquifer parameters, such as hydraulic conductivity (K) and storativity (S) values. The values obtained somehow incorporate some averaging values of aquifer heterogeneity, while the averaging functions are a direct consequence of the method of analysis employed. In recent years most work, casted in a stochastic framework, focused on the relationship between pumping rate and ensemble mean or variance of drawdown, thus having to pre-specify the parameters characterizing the underlying random spatial function. On the contrary, we contend that additional highly relevant information about heterogeneity can be obtained by looking to the spatial distribution of drawdown in individual realizations of the heterogeneous K field, without the need for invoking ergodic arguments. We present an analysis of the spatial distribution of time-dependent drawdown in a tridimensional aquifer produced by constant rate pumping in a fully penetrating well. The aquifer is considered of infinite extension in the x, y directions, and we assume no-flow boundaries in the aquifer top and bottom. The observation point is a fully penetrating piezometer. We consider an unknown spatial distribution of K(x,y,z), and using a perturbation expansion up to second order, we look at the late-time behavior of drawdown at any given observation vertical line. We conclude that: (1) at any given observation line the late-time behavior of drawdown would display a straight line in a drawdown versus log time plot, thus allowing the use of Jacob's method for test interpretation; (2) the slope of the straight line is the same for each observation line, thus providing a global average of K(x,y,z) through the aquifer; (3) the intercept point of the line in the same plot depends on location and is related to connectivity issues between the pumping and observation locations; (4) the intercept value is a weighted function of the local
Multi-crosswell profile 3D imaging and method
Washbourne, John K.; Rector, III, James W.; Bube, Kenneth P.
2002-01-01
Characterizing the value of a particular property, for example, seismic velocity, of a subsurface region of ground is described. In one aspect, the value of the particular property is represented using at least one continuous analytic function such as a Chebychev polynomial. The seismic data may include data derived from at least one crosswell dataset for the subsurface region of interest and may also include other data. In either instance, data may simultaneously be used from a first crosswell dataset in conjunction with one or more other crosswell datasets and/or with the other data. In another aspect, the value of the property is characterized in three dimensions throughout the region of interest using crosswell and/or other data. In still another aspect, crosswell datasets for highly deviated or horizontal boreholes are inherently useful. The method is performed, in part, by fitting a set of vertically spaced layer boundaries, represented by an analytic function such as a Chebychev polynomial, within and across the region encompassing the boreholes such that a series of layers is defined between the layer boundaries. Initial values of the particular property are then established between the layer boundaries and across the subterranean region using a series of continuous analytic functions. The continuous analytic functions are then adjusted to more closely match the value of the particular property across the subterranean region of ground to determine the value of the particular property for any selected point within the region.
Improvement of advanced nodal method used in 3D core design system
Rauck, S.; Dall'Osso, A.
2006-07-01
This paper deals with AREVA NP progress in the modelling of neutronic phenomena, evaluated through 3D determinist core codes and using 2-group diffusion theory. Our report highlights the advantages of taking into account the assembly environment in the process used for the building of the 2-group collapsed neutronic parameters, such as cross sections or discontinuity factors. The interest of the present method, developed in order to account for the impact of the environment on the above mentioned parameters, resides (i) in the very definition of a global correlation between collapsed neutronic data calculated in an infinite medium and those calculated in a 3D-geometry, and (ii) in the use of a re-homogenization method. Using this approach, computations match better with actual measurements on control rod worth. They also present smaller differences on pin by pin power values compared to the ones computed with another code considered as a reference since it relies on multigroup transport theory. (authors)
A method of 3-D data information storage with virtual holography
NASA Astrophysics Data System (ADS)
Huang, Zhen; Liu, Guodong; Ren, Zhong; Zeng, Lüming
2008-12-01
In this paper, a new method of 3-D data cube based on virtual holographic storage is presented. Firstly, the data information is encoded in the form of 3-D data cube with a certain algorithm, in which the interval along coordinates between every data is d. Using the plane-scanning method, the 3-D cube can be described as a assembly of slices which are parallel planes along the coordinates at an interval of d. The dot on the slice represents a bit. The bright one means "1", while the dark one means "0". Secondly, a hologram of the 3-D cube is obtained by computer with virtual optics technology. All the information of a 3-D cube can be described by a 2-D hologram. At last, the hologram is inputted in the SLM, and recorded in the recording material by intersecting two coherent laser beams. When the 3-D data is exported, a reference light illuminates the hologram, and a CCD is used to get the object image which is a hologram of the 3-D data. Then the 3-D data is computed with virtual optical technology. Compared with 2-D data page storage, the 3-D data cube storage has outstanding performance in larger capacity of data storage and higher security of data.
Quasi 3D modeling of water flow and solute transport in vadose zone and groundwater
NASA Astrophysics Data System (ADS)
Yakirevich, A.; Kuznetsov, M.; Weisbrod, N.; Pachepsky, Y. A.
2013-12-01
The complexity of subsurface flow systems calls for a variety of concepts leading to the multiplicity of simplified flow models. One commonly used simplification is based on the assumption that lateral flow and transport in unsaturated zone is insignificant unless the capillary fringe is involved. In such cases the flow and transport in the unsaturated zone above groundwater level can be simulated as a 1D phenomenon, whereas through groundwater they are viewed as 2D or 3D phenomena. A new approach for a numerical scheme for 3D variably saturated flow and transport is presented. A Quasi-3D approach allows representing flow in the 'vadose zone - aquifer' system by a series of 1D Richards' equations solved in variably-saturated zone and by 3D-saturated flow equation in groundwater (modified MODFLOW code). The 1D and 3D equations are coupled at the phreatic surface in a way that aquifer replenishment is calculated using the Richards' equation, and solving for the moving water table does not require definition of the specific yield parameter. The 3D advection-dispersion equation is solved in the entire domain by the MT3D code. Using implicit finite differences approximation to couple processes in the vadose zone and groundwater provides mass conservation and increase of computational efficiency. The above model was applied to simulate the impact of irrigation on groundwater salinity in the Alto Piura aquifer (Northern Peru). Studies on changing groundwater quality in arid and semi-arid lands show that irrigation return flow is one of the major factors contributing to aquifer salinization. Existing mathematical models do not account explicitly for the solute recycling during irrigation on a daily scale. Recycling occurs throughout the unsaturated and saturated zones, as function of the solute mass extracted from pumping wells. Salt concentration in irrigation water is calculated at each time step as a function of concentration of both surface water and groundwater
Optically directed molecular transport and 3D isoelectric positioning of amphoteric biomolecules
Hafeman, Dean G.; Harkins, James B.; Witkowski, Charles E.; Lewis, Nathan S.; Warmack, Robert J.; Brown, Gilbert M.; Thundat, Thomas
2006-01-01
We demonstrate the formation of charged molecular packets and their transport within optically created electrical force-field traps in a pH-buffered electrolyte. We call this process photoelectrophoretic localization and transport (PELT). The electrolyte is in contact with a photoconductive semiconductor electrode and a counterelectrode that are connected through an external circuit. A light beam directed to coordinates on the photoconductive electrode surface produces a photocurrent within the circuit and electrolyte. Within the electrolyte, the photocurrent creates localized force-field traps centered at the illuminated coordinates. Charged molecules, including polypeptides and proteins, electrophoretically accumulate into the traps and subsequently can be transported in the electrolyte by moving the traps over the photoconductive electrode in response to movement of the light beam. The molecules in a single trap can be divided into aliquots, and the aliquots can be directed along multiple routes simultaneously by using multiple light beams. This photoelectrophoretic transport of charged molecules by PELT resembles the electrostatic transport of electrons within force-field wells of solid-state charge-coupled devices. The molecules, however, travel in a liquid electrolyte rather than a solid. Furthermore, we have used PELT to position amphoteric biomolecules in three dimensions. A 3D pH gradient was created in an electrolyte medium by controlling the illumination position on a photoconductive anode where protons were generated electrolytically. Photoelectrophoretic transport of amphoteric molecules through the pH gradient resulted in accumulation of the molecules at their apparent 3D isoelectric coordinates in the medium. PMID:16618926
Accurate, finite-volume methods for 3D MHD on unstructured Lagrangian meshes
Barnes, D.C.; Rousculp, C.L.
1998-10-01
Previous 2D methods for magnetohydrodynamics (MHD) have contributed both to development of core code capability and to physics applications relevant to AGEX pulsed-power experiments. This strategy is being extended to 3D by development of a modular extension of an ASCI code. Extension to 3D not only increases complexity by problem size, but also introduces new physics, such as magnetic helicity transport. The authors have developed a method which incorporates all known conservation properties into the difference scheme on a Lagrangian unstructured mesh. Because the method does not depend on the mesh structure, mesh refinement is possible during a calculation to prevent the well known problem of mesh tangling. Arbitrary polyhedral cells are decomposed into tetrahedrons. The action of the magnetic vector potential, A {center_dot} {delta}l, is centered on the edges of this extended mesh. For ideal flow, this maintains {del} {center_dot} B = 0 to round-off error. Vertex forces are derived by the variation of magnetic energy with respect to vertex positions, F = {minus}{partial_derivative}W{sub B}/{partial_derivative}r. This assures symmetry as well as magnetic flux, momentum, and energy conservation. The method is local so that parallelization by domain decomposition is natural for large meshes. In addition, a simple, ideal-gas, finite pressure term has been included. The resistive diffusion part is calculated using the support operator method, to obtain an energy conservative, symmetric method on an arbitrary mesh. Implicit time difference equations are solved by preconditioned, conjugate gradient methods. Results of convergence tests are presented. Initial results of an annular Z-pinch implosion problem illustrate the application of these methods to multi-material problems.
Evaluation of a new method for stenosis quantification from 3D x-ray angiography images
NASA Astrophysics Data System (ADS)
Betting, Fabienne; Moris, Gilles; Knoplioch, Jerome; Trousset, Yves L.; Sureda, Francisco; Launay, Laurent
2001-05-01
A new method for stenosis quantification from 3D X-ray angiography images has been evaluated on both phantom and clinical data. On phantoms, for the parts larger or equal to 3 mm, the standard deviation of the measurement error has always found to be less or equal to 0.4 mm, and the maximum measurement error less than 0.17 mm. No clear relationship has been observed between the performances of the quantification method and the acquisition FoV. On clinical data, the 3D quantification method proved to be more robust to vessel bifurcations than its 3D equivalent. On a total of 15 clinical cases, the differences between 2D and 3D quantification were always less than 0.7 mm. The conclusion is that stenosis quantification from 3D X-4ay angiography images is an attractive alternative to quantification from 2D X-ray images.
Reconstructing photorealistic 3D models from image sequence using domain decomposition method
NASA Astrophysics Data System (ADS)
Xiong, Hanwei; Pan, Ming; Zhang, Xiangwei
2009-11-01
In the fields of industrial design, artistic design and heritage conservation, physical objects are usually digitalized by reverse engineering through some 3D scanning methods. Structured light and photogrammetry are two main methods to acquire 3D information, and both are expensive. Even if these expensive instruments are used, photorealistic 3D models are seldom available. In this paper, a new method to reconstruction photorealistic 3D models using a single camera is proposed. A square plate glued with coded marks is used to place the objects, and a sequence of about 20 images is taken. From the coded marks, the images are calibrated, and a snake algorithm is used to segment object from the background. A rough 3d model is obtained using shape from silhouettes algorithm. The silhouettes are decomposed into a combination of convex curves, which are used to partition the rough 3d model into some convex mesh patches. For each patch, the multi-view photo consistency constraints and smooth regulations are expressed as a finite element formulation, which can be resolved locally, and the information can be exchanged along the patches boundaries. The rough model is deformed into a fine 3d model through such a domain decomposition finite element method. The textures are assigned to each element mesh, and a photorealistic 3D model is got finally. A toy pig is used to verify the algorithm, and the result is exciting.
NASA Astrophysics Data System (ADS)
Daughney, C.; Toews, M. W.; Morgenstern, U.; Cornaton, F. J.; Jackson, B. M.
2013-12-01
Lake Rotorua is a focus of culture and tourism in New Zealand. The lake's water quality has declined since the 1970s, partly due to nutrient inputs that reach the lake via the groundwater system. Improved land use management within the catchment requires prediction of the spatial variations of groundwater transit time from land surface to the lake, and from this the prediction of current and future nutrient inflows to the lake. This study combines the two main methods currently available for determination of water age: numerical groundwater models and hydrological tracers. A steady-state 3D finite element model was constructed to simulate groundwater flow and transport of tritium and age at the catchment scale (555 km2). The model materials were defined using a 3D geologic model and included ignimbrites, rhyolites, alluvial and lake bottom sediments. The steady-state saturated groundwater flow model was calibrated using observed groundwater levels in boreholes (111 locations) and stream flow measurements from groundwater-fed streams and springs (61 locations). Hydraulic conductivities and Cauchy boundary conditions associated with the streams, springs and lake were parameterized. The transport parameters for the model were calibrated using 191 tritium samples from 105 locations (springs, streams and boreholes), with most locations having two sample dates. The transport model used steady-state flow, but simulated the transient transport and decay of tritium from rainfall recharge between 1945 and 2012. An additional 1D unsaturated sub-model was added to account for tritium decay from the ground surface to the water table. The sub-model is linked on top of the 3D model, and uses the water table depths and material properties from the 3D model. The adjustable calibration parameters for the transport model were porosity and van Genuchten parameters related to the unsaturated sub-models. Calibration of the flow model was achieved using a combination of automated least
3D-2D registration of cerebral angiograms: a method and evaluation on clinical images.
Mitrovic, Uroš; Špiclin, Žiga; Likar, Boštjan; Pernuš, Franjo
2013-08-01
Endovascular image-guided interventions (EIGI) involve navigation of a catheter through the vasculature followed by application of treatment at the site of anomaly using live 2D projection images for guidance. 3D images acquired prior to EIGI are used to quantify the vascular anomaly and plan the intervention. If fused with the information of live 2D images they can also facilitate navigation and treatment. For this purpose 3D-2D image registration is required. Although several 3D-2D registration methods for EIGI achieve registration accuracy below 1 mm, their clinical application is still limited by insufficient robustness or reliability. In this paper, we propose a 3D-2D registration method based on matching a 3D vasculature model to intensity gradients of live 2D images. To objectively validate 3D-2D registration methods, we acquired a clinical image database of 10 patients undergoing cerebral EIGI and established "gold standard" registrations by aligning fiducial markers in 3D and 2D images. The proposed method had mean registration accuracy below 0.65 mm, which was comparable to tested state-of-the-art methods, and execution time below 1 s. With the highest rate of successful registrations and the highest capture range the proposed method was the most robust and thus a good candidate for application in EIGI. PMID:23649179
3D Space Radiation Transport in a Shielded ICRU Tissue Sphere
NASA Technical Reports Server (NTRS)
Wilson, John W.; Slaba, Tony C.; Badavi, Francis F.; Reddell, Brandon D.; Bahadori, Amir A.
2014-01-01
A computationally efficient 3DHZETRN code capable of simulating High Charge (Z) and Energy (HZE) and light ions (including neutrons) under space-like boundary conditions with enhanced neutron and light ion propagation was recently developed for a simple homogeneous shield object. Monte Carlo benchmarks were used to verify the methodology in slab and spherical geometry, and the 3D corrections were shown to provide significant improvement over the straight-ahead approximation in some cases. In the present report, the new algorithms with well-defined convergence criteria are extended to inhomogeneous media within a shielded tissue slab and a shielded tissue sphere and tested against Monte Carlo simulation to verify the solution methods. The 3D corrections are again found to more accurately describe the neutron and light ion fluence spectra as compared to the straight-ahead approximation. These computationally efficient methods provide a basis for software capable of space shield analysis and optimization.
Momentum Transport: 2D and 3D Cloud Resolving Model Simulations
NASA Technical Reports Server (NTRS)
Tao, Wei-Kuo
2001-01-01
The major objective of this study is to investigate the momentum budgets associated with several convective systems that developed during the TOGA COARE IOP (west Pacific warm pool region) and GATE (east Atlantic region). The tool for this study is the improved Goddard Cumulas Ensemble (GCE) model which includes a 3-class ice-phase microphysical scheme, explicit cloud radiative interactive processes and air-sea interactive surface processes. The model domain contains 256 x 256 grid points (with 2 km resolution) in the horizontal and 38 grid points (to a depth of 22 km) in the vertical. The 2D domain has 1024 grid points. The simulations were performed over a 7-day time period (December 19-26, 1992, for TOGA COARE and September 1-7, 1994 for GATE). Cyclic literal boundary conditions are required for this type of long-term integration. Two well organized squall systems (TOGA, COARE February 22, 1993, and GATE September 12, 1994) were also simulated using the 3D GCE model. Only 9 h simulations were required to cover the life time of the squall systems. the lateral boundary conditions were open for these two squall systems simulations. the following will be examined: (1) the momentum budgets in the convective and stratiform regions, (2) the relationship between momentum transport and cloud organization (i.e., well organized squall lines versus less organized convective), (3) the differences and similarities in momentum transport between 2D and 3D simulated convective systems, and (4) the differences and similarities in momentum budgets between cloud systems simulated with open and cyclic lateral boundary conditions. Preliminary results indicate that there are only small differences between 2D and 3D simulated momentum budgets. Major differences occur, however, between momentum budgets associated with squall systems simulated using different lateral boundary conditions.
M3D-K simulations of sawteeth and energetic particle transport in tokamak plasmas
NASA Astrophysics Data System (ADS)
Shen, Wei; Fu, G. Y.; Sheng, Zheng-Mao; Breslau, J. A.; Wang, Feng
2014-09-01
Nonlinear simulations of sawteeth and related energetic particle transport are carried out using the kinetic/magnetohydrodynamic (MHD) hybrid code M3D-K. MHD simulations show repeated sawtooth cycles for a model tokamak equilibrium. Furthermore, test particle simulations are carried out to study the energetic particle transport due to a sawtooth crash. The results show that energetic particles are redistributed radially in the plasma core, depending on pitch angle and energy. For trapped particles, the redistribution occurs for particle energy below a critical value in agreement with existing theories. For co-passing particles, the redistribution is strong with little dependence on particle energy. In contrast, the redistribution level of counter-passing particles decreases with increasing particle energy.
M3D-K Simulations of Sawteeth and Energetic Particle Transport in Tokamak Plasmas
NASA Astrophysics Data System (ADS)
Shen, Wei; Fu, Guoyong; Sheng, Zhengmao; Breslau, Joshua; Wang, Feng
2013-10-01
Nonlinear simulations of Sawteeth and energetic particle transport are carried out using the kinetic/MHD hybrid code M3D-K. MHD simulations show repeated sawtooth cycles due to a resistive (1,1) internal kink mode for a model tokamak equilibrium. Furthermore, test particle simulations are carried out to study the energetic particle transport due to a sawtooth crash. The results show that energetic particles are redistributed radially in plasma core depending on pitch angle and energy. For trapped particles, the redistribution occurs for particle energy below a critical value in agreement with previous theory. For co-passing particles, the redistribution is strong with little dependence on particle energy. In contrast, the redistribution level of counter-passing particles decreases as particle energy becomes large.
M3D-K simulations of sawteeth and energetic particle transport in tokamak plasmas
Shen, Wei; Sheng, Zheng-Mao; Fu, G. Y.; Breslau, J. A.; Wang, Feng
2014-09-15
Nonlinear simulations of sawteeth and related energetic particle transport are carried out using the kinetic/magnetohydrodynamic (MHD) hybrid code M3D-K. MHD simulations show repeated sawtooth cycles for a model tokamak equilibrium. Furthermore, test particle simulations are carried out to study the energetic particle transport due to a sawtooth crash. The results show that energetic particles are redistributed radially in the plasma core, depending on pitch angle and energy. For trapped particles, the redistribution occurs for particle energy below a critical value in agreement with existing theories. For co-passing particles, the redistribution is strong with little dependence on particle energy. In contrast, the redistribution level of counter-passing particles decreases with increasing particle energy.
3D numerical simulation of the transport of chemical signature compounds from buried landmines
NASA Astrophysics Data System (ADS)
Irrazabal, Maik; Borrero, Ernesto; Briano, Julio G.; Castro, Miguel; Hernandez, Samuel P.
2005-06-01
The transport of the chemical signature compounds from buried landmines in a three-dimensional (3D) array has been numerically modeled using the finite-volume technique. Compounds such as trinitrotoluene, dinitrotoluene, and their degradation products, are semi volatile and somewhat soluble in water. Furthermore, they can strongly adsorb to the soil and undergo chemical and biological degradation. Consequently, the spatial and temporal concentration distributions of such chemicals depend on the mobility of the water and gaseous phases, their molecular and mechanical diffusion, adsorption characteristics, soil water content, compaction, and environmental factors. A 3D framework is required since two-dimensional (2D) symmetry may easily fade due to terrain topography: non-flat surfaces, soil heterogeneity, or underground fractures. The spatial and temporal distribution of the chemical-signature-compounds, in an inclined grid has been obtained. The fact that the chemicals may migrate horizontally, giving higher surface concentrations at positions not directly on top of the objects, emphasizes the need for understanding the transport mechanism when a chemical detector is used. Deformation in the concentration contours after rainfall is observed in the inclined surface and is attributed to both: the advective flux, and to the water flux at the surface caused by the slope. The analysis of the displacements in the position of the maximum concentrations at the surface, respect to the actual location of the mine, in an inclined system, is presented.
Analysis of transport connectivity in karstic aquifers spanned by 3D conduit networks
NASA Astrophysics Data System (ADS)
Ronayne, M. J.
2013-12-01
Karst aquifers are characterized by interconnected conduits that behave as structural pathways for groundwater and solutes. This modeling study assesses the influence of conduit network geometry on solute transport behavior within karst systems. Synthetic karst aquifers containing 3D conduit networks were considered. Networks of varying complexity were generated using a directed percolation model. Flow and transport simulations were conducted for each synthetic aquifer by modeling the conduits as discretized high-permeability features within a uniform matrix material. Transport connectivity and dispersive properties were evaluated using statistical moments of the solute arrival time distribution at the downgradient conduit outlet (karst spring). In addition, a new connectivity metric that quantifies solute residence time within conduits was considered. Results show that a more complex network leads to enhanced mixing between the conduit and matrix domains, which has the effect of reducing transport connectivity. This modeling study illustrates how typically available transport data (e.g., solute breakthrough curves at the conduit outlet) may reveal information about the internal network structure, thus providing guidance for future inverse modeling.
A Monte Carlo method for 3D thermal infrared radiative transfer
NASA Astrophysics Data System (ADS)
Chen, Y.; Liou, K. N.
2006-09-01
A 3D Monte Carlo model for specific application to the broadband thermal radiative transfer has been developed in which the emissivities for gases and cloud particles are parameterized by using a single cubic element as the building block in 3D space. For spectral integration in the thermal infrared, the correlated k-distribution method has been used for the sorting of gaseous absorption lines in multiple-scattering atmospheres involving 3D clouds. To check the Monte-Carlo simulation, we compare a variety of 1D broadband atmospheric fluxes and heating rates to those computed from the conventional plane-parallel (PP) model and demonstrate excellent agreement between the two. Comparisons of the Monte Carlo results for broadband thermal cooling rates in 3D clouds to those computed from the delta-diffusion approximation for 3D radiative transfer and the independent pixel-by-pixel approximation are subsequently carried out to understand the relative merits of these approaches.
MODIS volcanic ash retrievals vs FALL3D transport model: a quantitative comparison
NASA Astrophysics Data System (ADS)
Corradini, S.; Merucci, L.; Folch, A.
2010-12-01
Satellite retrievals and transport models represents the key tools to monitor the volcanic clouds evolution. Because of the harming effects of fine ash particles on aircrafts, the real-time tracking and forecasting of volcanic clouds is key for aviation safety. Together with the security reasons also the economical consequences of a disruption of airports must be taken into account. The airport closures due to the recent Icelandic Eyjafjöll eruption caused millions of passengers to be stranded not only in Europe, but across the world. IATA (the International Air Transport Association) estimates that the worldwide airline industry has lost a total of about 2.5 billion of Euro during the disruption. Both security and economical issues require reliable and robust ash cloud retrievals and trajectory forecasting. The intercomparison between remote sensing and modeling is required to assure precise and reliable volcanic ash products. In this work we perform a quantitative comparison between Moderate Resolution Imaging Spectroradiometer (MODIS) retrievals of volcanic ash cloud mass and Aerosol Optical Depth (AOD) with the FALL3D ash dispersal model. MODIS, aboard the NASA-Terra and NASA-Aqua polar satellites, is a multispectral instrument with 36 spectral bands operating in the VIS-TIR spectral range and spatial resolution varying between 250 and 1000 m at nadir. The MODIS channels centered around 11 and 12 micron have been used for the ash retrievals through the Brightness Temperature Difference algorithm and MODTRAN simulations. FALL3D is a 3-D time-dependent Eulerian model for the transport and deposition of volcanic particles that outputs, among other variables, cloud column mass and AOD. Three MODIS images collected the October 28, 29 and 30 on Mt. Etna volcano during the 2002 eruption have been considered as test cases. The results show a general good agreement between the retrieved and the modeled volcanic clouds in the first 300 km from the vents. Even if the
Age, double porosity, and simple reaction modifications for the MOC3D ground-water transport model
Goode, Daniel J.
1999-01-01
This report documents modifications for the MOC3D ground-water transport model to simulate (a) ground-water age transport; (b) double-porosity exchange; and (c) simple but flexible retardation, decay, and zero-order growth reactions. These modifications are incorporated in MOC3D version 3.0. MOC3D simulates the transport of a single solute using the method-ofcharacteristics numerical procedure. The age of ground water, that is the time since recharge to the saturated zone, can be simulated using the transport model with an additional source term of unit strength, corresponding to the rate of aging. The output concentrations of the model are in this case the ages at all locations in the model. Double porosity generally refers to a separate immobilewater phase within the aquifer that does not contribute to ground-water flow but can affect solute transport through diffusive exchange. The solute mass exchange rate between the flowing water in the aquifer and the immobile-water phase is the product of the concentration difference between the two phases and a linear exchange coefficient. Conceptually, double porosity can approximate the effects of dead-end pores in a granular porous media, or matrix diffusion in a fractured-rock aquifer. Options are provided for decay and zero-order growth reactions within the immobilewater phase. The simple reaction terms here extend the original model, which included decay and retardation. With these extensions, (a) the retardation factor can vary spatially within each model layer, (b) the decay rate coefficient can vary spatially within each model layer and can be different for the dissolved and sorbed phases, and (c) a zero-order growth reaction is added that can vary spatially and can be different in the dissolved and sorbed phases. The decay and growth reaction terms also can change in time to account for changing geochemical conditions during transport. The report includes a description of the theoretical basis of the model, a
A simple method for producing freestanding 3D microstructures by integrated photomask micromolding
NASA Astrophysics Data System (ADS)
Li, Hui
2015-12-01
Freestanding three-dimensional (3D) microstructures are widely used in micro-electro-mechanical system (MEMS) applications or can function as microdevices themselves. However, microfabrication methods for freestanding 3D microstructures have limitations in shape, size, cost, and mass production, etc. In this work, integrated photomask micromolding is demonstrated, which uses a portable UV light source and chrome glass micromolding to fabricate 3D microstructures without alignment. Specifically, a chrome layer on one side of the glass micromold shields the excess filling SU-8 photoresist from UV exposure and only the SU-8 photoresist in mold cavities is crosslinked. The 3D microstructures produced using this method have very high dimensional accuracy and the profile error is approximately 1.5%. This method can be used with features of virtually any size and shape and can be integrated into highly-parallel micromolding processes and has potential for MEMS applications.
NASA Astrophysics Data System (ADS)
Wei, Bing; Wang, Lidong; Wang, Yang; Yuan, Yinan; Miao, Qinghua; Yang, Ziyue; Fei, Weidong
2016-03-01
In this study, a new, effective strategy is reported for the fabrication of composites using manganese oxide (MnO2) grown in situ on three-dimensional (3D) graphene by the reverse microemulsion (water-in-oil) method. A uniform coating of nanoscale MnO2 layers can be observed on the internal surface of 3D graphene, which could benefit rapid ionic and electronic transport. The electrochemical performance of the MnO2/3D graphene composites is optimized by the control of the composite structure and mass loading of MnO2. The MnO2/3D graphene composite thus prepared exhibits a significantly high specific capacitance of 659.7 F g-1 at 0.3 A g-1 and an excellent retention life of 106% after 1000 cycles. The facile synthesis and excellent electrochemical performance of the MnO2/3D graphene composites indicate that the developed method demonstrates potential applications for the fabrication of novel electrode materials for use in energy storage devices.
3D modeling method for computer animate based on modified weak structured light method
NASA Astrophysics Data System (ADS)
Xiong, Hanwei; Pan, Ming; Zhang, Xiangwei
2010-11-01
A simple and affordable 3D scanner is designed in this paper. Three-dimensional digital models are playing an increasingly important role in many fields, such as computer animate, industrial design, artistic design and heritage conservation. For many complex shapes, optical measurement systems are indispensable to acquiring the 3D information. In the field of computer animate, such an optical measurement device is too expensive to be widely adopted, and on the other hand, the precision is not as critical a factor in that situation. In this paper, a new cheap 3D measurement system is implemented based on modified weak structured light, using only a video camera, a light source and a straight stick rotating on a fixed axis. For an ordinary weak structured light configuration, one or two reference planes are required, and the shadows on these planes must be tracked in the scanning process, which destroy the convenience of this method. In the modified system, reference planes are unnecessary, and size range of the scanned objects is expanded widely. A new calibration procedure is also realized for the proposed method, and points cloud is obtained by analyzing the shadow strips on the object. A two-stage ICP algorithm is used to merge the points cloud from different viewpoints to get a full description of the object, and after a series of operations, a NURBS surface model is generated in the end. A complex toy bear is used to verify the efficiency of the method, and errors range from 0.7783mm to 1.4326mm comparing with the ground truth measurement.
BioFVM: an efficient, parallelized diffusive transport solver for 3-D biological simulations
Ghaffarizadeh, Ahmadreza; Friedman, Samuel H.; Macklin, Paul
2016-01-01
Motivation: Computational models of multicellular systems require solving systems of PDEs for release, uptake, decay and diffusion of multiple substrates in 3D, particularly when incorporating the impact of drugs, growth substrates and signaling factors on cell receptors and subcellular systems biology. Results: We introduce BioFVM, a diffusive transport solver tailored to biological problems. BioFVM can simulate release and uptake of many substrates by cell and bulk sources, diffusion and decay in large 3D domains. It has been parallelized with OpenMP, allowing efficient simulations on desktop workstations or single supercomputer nodes. The code is stable even for large time steps, with linear computational cost scalings. Solutions are first-order accurate in time and second-order accurate in space. The code can be run by itself or as part of a larger simulator. Availability and implementation: BioFVM is written in C ++ with parallelization in OpenMP. It is maintained and available for download at http://BioFVM.MathCancer.org and http://BioFVM.sf.net under the Apache License (v2.0). Contact: paul.macklin@usc.edu. Supplementary information: Supplementary data are available at Bioinformatics online. PMID:26656933
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.
Accurate compressed look up table method for CGH in 3D holographic display.
Gao, Chuan; Liu, Juan; Li, Xin; Xue, Gaolei; Jia, Jia; Wang, Yongtian
2015-12-28
Computer generated hologram (CGH) should be obtained with high accuracy and high speed in 3D holographic display, and most researches focus on the high speed. In this paper, a simple and effective computation method for CGH is proposed based on Fresnel diffraction theory and look up table. Numerical simulations and optical experiments are performed to demonstrate its feasibility. The proposed method can obtain more accurate reconstructed images with lower memory usage compared with split look up table method and compressed look up table method without sacrificing the computational speed in holograms generation, so it is called accurate compressed look up table method (AC-LUT). It is believed that AC-LUT method is an effective method to calculate the CGH of 3D objects for real-time 3D holographic display where the huge information data is required, and it could provide fast and accurate digital transmission in various dynamic optical fields in the future. PMID:26831987
Time-dependent 3-D dterministic transport on parallel architectures using Dantsys/MPI
Baker, R.S.; Alcouffe, R.E.
1996-12-31
In addition to the ability to solve the static transport equation, we have also incorporated time dependence into our parallel 3-D S{sub {ital N}} code DANTSYS/MPI. Using a semi-implicit scheme, DANTSYS/MPI is capable of performing time-dependent calculations for both fissioning and pure source driven problems. We have applied this to various types of problems such as nuclear well logging and prompt fission experiments. This paper describes the form of the time- dependent equations implemented, their solution strategies in DANTSYS/MPI including iteration acceleration, and the strategies used for time-step control. Results are presented for a model nuclear well logging calculation.
Distributed network of integrated 3D sensors for transportation security applications
NASA Astrophysics Data System (ADS)
Hejmadi, Vic; Garcia, Fred
2009-05-01
The US Port Security Agency has strongly emphasized the needs for tighter control at transportation hubs. Distributed arrays of miniature CMOS cameras are providing some solutions today. However, due to the high bandwidth required and the low valued content of such cameras (simple video feed), large computing power and analysis algorithms as well as control software are needed, which makes such an architecture cumbersome, heavy, slow and expensive. We present a novel technique by integrating cheap and mass replicable stealth 3D sensing micro-devices in a distributed network. These micro-sensors are based on conventional structures illumination via successive fringe patterns on the object to be sensed. The communication bandwidth between each sensor remains very small, but is of very high valued content. Key technologies to integrate such a sensor are digital optics and structured laser illumination.
Development of a 3D to 1D Particle Transport Model to Predict Deposition in the Lungs
NASA Astrophysics Data System (ADS)
Oakes, Jessica M.; Grandmont, Celine; Shadden, Shawn C.; Vignon-Clementel, Irene E.
2014-11-01
Aerosolized particles are commonly used for therapeutic drug delivery as they can be delivered to the body systemically or be used to treat lung diseases. Recent advances in computational resources have allowed for sophisticated pulmonary simulations, however it is currently impossible to solve for airflow and particle transport for all length and time scales of the lung. Instead, multi-scale methods must be used. In our recent work, where computational methods were employed to solve for airflow and particle transport in the rat airways (Oakes et al. (2014), Annals of Biomedical Engineering 42, 899), the number of particles to exit downstream of the 3D domain was determined. In this current work, the time-dependent Lagrangian description of particles was used to numerically solve a 1D convection-diffusion model (trumpet model, Taulbee and Yu (1975), Journal of Applied Physiology, 38, 77) parameterized specifically for the lung. The expansion of the airway dimensions was determined based on data collected from our aerosol exposure experiments (Oakes et al. (2014), Journal of Applied Physiology, 116, 1561). This 3D-1D framework enables us to predict the fate of particles in the whole lung. This work was supported by the Whitaker Foundation at the IIE, a INRIA Associated Team Postdoc Grant, and a UC Presidential Fellowship.
A multimaterial bioink method for 3D printing tunable, cell-compatible hydrogels.
Rutz, Alexandra L; Hyland, Kelly E; Jakus, Adam E; Burghardt, Wesley R; Shah, Ramille N
2015-03-01
A multimaterial bio-ink method using polyethylene glycol crosslinking is presented for expanding the biomaterial palette required for 3D bioprinting of more mimetic and customizable tissue and organ constructs. Lightly crosslinked, soft hydrogels are produced from precursor solutions of various materials and 3D printed. Rheological and biological characterizations are presented, and the promise of this new bio-ink synthesis strategy is discussed. PMID:25641220
Geometrically-compatible 3-D Monte Carlo and discrete-ordinates methods
Morel, J.E.; Wareing, T.A.; McGhee, J.M.; Evans, T.M.
1998-12-31
This is the final report of a three-year, Laboratory Directed Research and Development (LDRD) project at the Los Alamos National Laboratory (LANL). The purpose of this project was two-fold. The first purpose was to develop a deterministic discrete-ordinates neutral-particle transport scheme for unstructured tetrahedral spatial meshes, and implement it in a computer code. The second purpose was to modify the MCNP Monte Carlo radiation transport code to use adjoint solutions from the tetrahedral-mesh discrete-ordinates code to reduce the statistical variance of Monte Carlo solutions via a weight-window approach. The first task has resulted in a deterministic transport code that is much more efficient for modeling complex 3-D geometries than any previously existing deterministic code. The second task has resulted in a powerful new capability for dramatically reducing the cost of difficult 3-D Monte Carlo calculations.
Multi-scale self-organisation of edge plasma turbulent transport in 3D global simulations
NASA Astrophysics Data System (ADS)
Tamain, P.; Ghendrih, Ph; Bufferand, H.; Ciraolo, G.; Colin, C.; Fedorczak, N.; Nace, N.; Schwander, F.; Serre, E.
2015-05-01
The 3D global edge turbulence code TOKAM3X is used to study the properties of edge particle turbulent transport in circular limited plasmas, including both closed and open flux surfaces. Turbulence is driven by an incoming particle flux from the core plasma and no scale separation between the equilibrium and the fluctuations is assumed. Simulations show the existence of a complex self-organization of turbulence transport coupling scales ranging from a few Larmor radii up to the machine scale. Particle transport is largely dominated by small scale turbulence with fluctuations forming quasi field-aligned filaments. Radial particle transport is intermittent and associated with the propagation of coherent structures on long distances via avalanches. Long range correlations are also found in the poloidal and toroidal direction. The statistical properties of fluctuations vary with the radial and poloidal directions, with larger fluctuation levels and intermittency found in the outboard scrape-off layer (SOL). Radial turbulent transport is strongly ballooned, with 90% of the flux at the separatrix flowing through the low-field side. One of the main consequences is the existence of quasi-sonic asymmetric parallel flows driving a net rotation of the plasma. Simulations also show the spontaneous onset of an intermittent E × B rotation characterized by a larger shear at the separatrix. Strong correlation is found between the turbulent particle flux and the E × B flow shear in a phenomenology reminiscent of H-mode physics. The poloidal position of the limiter is a key player in the observed dynamics.
A practical salient region feature based 3D multi-modality registration method for medical images
NASA Astrophysics Data System (ADS)
Hahn, Dieter A.; Wolz, Gabriele; Sun, Yiyong; Hornegger, Joachim; Sauer, Frank; Kuwert, Torsten; Xu, Chenyang
2006-03-01
We present a novel representation of 3D salient region features and its integration into a hybrid rigid-body registration framework. We adopt scale, translation and rotation invariance properties of those intrinsic 3D features to estimate a transform between underlying mono- or multi-modal 3D medical images. Our method combines advantageous aspects of both feature- and intensity-based approaches and consists of three steps: an automatic extraction of a set of 3D salient region features on each image, a robust estimation of correspondences and their sub-pixel accurate refinement with outliers elimination. We propose a region-growing based approach for the extraction of 3D salient region features, a solution to the problem of feature clustering and a reduction of the correspondence search space complexity. Results of the developed algorithm are presented for both mono- and multi-modal intra-patient 3D image pairs (CT, PET and SPECT) that have been acquired for change detection, tumor localization, and time based intra-person studies. The accuracy of the method is clinically evaluated by a medical expert with an approach that measures the distance between a set of selected corresponding points consisting of both anatomical and functional structures or lesion sites. This demonstrates the robustness of the proposed method to image overlap, missing information and artefacts. We conclude by discussing potential medical applications and possibilities for integration into a non-rigid registration framework.
NASA Astrophysics Data System (ADS)
Bai, Lu; Wang, Hongbo; Zhou, Jiangfan; Yang, Rong; Zhang, Hui
2014-11-01
In this paper, the aperture change of the aluminium alloy aerospace structure under real load is researched. Static experiments are carried on which is simulated the load environment of flight course. Compared with the traditional methods, through experiments results, it's proved that 3D digital speckle correlation method has good adaptability and precision on testing aperture change, and it can satisfy measurement on non-contact,real-time 3D deformation or stress concentration. The test results of new method is compared with the traditional method.
High efficient methods of content-based 3D model retrieval
NASA Astrophysics Data System (ADS)
Wu, Yuanhao; Tian, Ling; Li, Chenggang
2013-03-01
Content-based 3D model retrieval is of great help to facilitate the reuse of existing designs and to inspire designers during conceptual design. However, there is still a gap to apply it in industry due to the low time efficiency. This paper presents two new methods with high efficiency to build a Content-based 3D model retrieval system. First, an improvement is made on the "Shape Distribution (D2)" algorithm, and a new algorithm named "Quick D2" is proposed. Four sample 3D mechanical models are used in an experiment to compare the time cost of the two algorithms. The result indicates that the time cost of Quick D2 is much lower than that of D2, while the descriptors extracted by the two algorithms are almost the same. Second, an expandable 3D model repository index method with high performance, namely, RBK index, is presented. On the basis of RBK index, the search space is pruned effectively during the search process, leading to a speed up of the whole system. The factors that influence the values of the key parameters of RBK index are discussed and an experimental method to find the optimal values of the key parameters is given. Finally, "3D Searcher", a content-based 3D model retrieval system is developed. By using the methods proposed, the time cost for the system to respond one query online is reduced by 75% on average. The system has been implemented in a manufacturing enterprise, and practical query examples during a case of the automobile rear axle design are also shown. The research method presented shows a new research perspective and can effectively improve the content-based 3D model retrieval efficiency.
Modified Anderson Method for Accelerating 3D-RISM Calculations Using Graphics Processing Unit.
Maruyama, Yutaka; Hirata, Fumio
2012-09-11
A fast algorithm is proposed to solve the three-dimensional reference interaction site model (3D-RISM) theory on a graphics processing unit (GPU). 3D-RISM theory is a powerful tool for investigating biomolecular processes in solution; however, such calculations are often both memory-intensive and time-consuming. We sought to accelerate these calculations using GPUs, but to work around the problem of limited memory size in GPUs, we modified the less memory-intensive "Anderson method" to give faster convergence to 3D-RISM calculations. Using this method on a Tesla C2070 GPU, we reduced the total computational time by a factor of 8, 1.4 times by the modified Andersen method and 5.7 times by GPU, compared to calculations on an Intel Xeon machine (eight cores, 3.33 GHz) with the conventional method. PMID:26605714
Sippl, Wolfgang
2002-12-01
We have recently reported the development of a 3-D QSAR model for estrogen receptor ligands showing a significant correlation between calculated molecular interaction fields and experimentally measured binding affinity. The ligand alignment obtained from docking simulations was taken as basis for a comparative field analysis applying the GRID/GOLPE program. Using the interaction field derived with a water probe and applying the smart region definition (SRD) variable selection procedure, a significant and robust model was obtained (q(2)(LOO)=0.921, SDEP=0.345). To further analyze the robustness and the predictivity of the established model several recently developed estrogen receptor ligands were selected as external test set. An excellent agreement between predicted and experimental binding data was obtained indicated by an external SDEP of 0.531. Two other traditionally used prediction techniques were applied in order to check the performance of the receptor-based 3-D QSAR procedure. The interaction energies calculated on the basis of receptor-ligand complexes were correlated with experimentally observed affinities. Also ligand-based 3-D QSAR models were generated using program FlexS. The interaction energy-based model, as well as the ligand-based 3-D QSAR models yielded models with lower predictivity. The comparison with the interaction energy-based model and with the ligand-based 3-D QSAR models, respectively, indicates that the combination of receptor-based and 3-D QSAR methods is able to improve the quality of prediction. PMID:12413831
Implementation of algebraic stress models in a general 3-D Navier-Stokes method (PAB3D)
NASA Technical Reports Server (NTRS)
Abdol-Hamid, Khaled S.
1995-01-01
A three-dimensional multiblock Navier-Stokes code, PAB3D, which was developed for propulsion integration and general aerodynamic analysis, has been used extensively by NASA Langley and other organizations to perform both internal (exhaust) and external flow analysis of complex aircraft configurations. This code was designed to solve the simplified Reynolds Averaged Navier-Stokes equations. A two-equation k-epsilon turbulence model has been used with considerable success, especially for attached flows. Accurate predicting of transonic shock wave location and pressure recovery in separated flow regions has been more difficult. Two algebraic Reynolds stress models (ASM) have been recently implemented in the code that greatly improved the code's ability to predict these difficult flow conditions. Good agreement with Direct Numerical Simulation (DNS) for a subsonic flat plate was achieved with ASM's developed by Shih, Zhu, and Lumley and Gatski and Speziale. Good predictions were also achieved at subsonic and transonic Mach numbers for shock location and trailing edge boattail pressure recovery on a single-engine afterbody/nozzle model.
Simulation of bacteria transport processes in a river with Flow3D
NASA Astrophysics Data System (ADS)
Schwarzwälder, Kordula; Bui, Minh Duc; Rutschmann, Peter
2014-05-01
Water quality aspects are getting more and more important due to the European water Framework directive (WFD). One problem related to this topic is the inflow of untreated wastewater due to combined sewer overflows into a river. The wastewater mixture contains even bacteria like E. coli and Enterococci which are markers for water quality. In our work we investigated the transport of these bacteria in river Isar by using a large-scale flume in the outside area of our lab (Oskar von Miller Institute). Therefor we could collect basic data and knowledge about the processes which occur during bacteria sedimentation and remobilisation. In our flume we could use the real grain with the exact size distribution curve as in the river Isar which we want to simulate and we had the chance to nurture a biofilm which is realistic for the analysed situation. This biofilm plays an important role in the remobilisation processes, because the bacteria are hindered to be washed out back into the bulk phase as fast and in such an amount as this would happen without biofilm. The results of our experiments are now used for a module in the 3D software Flow3D to simulate the effects of a point source inlet of raw wastewater on the water quality. Therefor we have to implement the bacteria not as a problem of concentration with advection and diffusion but as single particles which can be inactivated during the process of settling and need to be hindered from remobilisation by the biofilm. This biofilm has special characteristic, it is slippery and has a special thickness which influences the chance of bacteria being removed. To achieve realistic results we have to include the biofilm with more than a probabilistic-tool to make sure that our module is transferable. The module should be as flexible as possible to be improved step by step with increasing quality of dataset.
Efficient fabrication method of nano-grating for 3D holographic display with full parallax views.
Wan, Wenqiang; Qiao, Wen; Huang, Wenbin; Zhu, Ming; Fang, Zongbao; Pu, Donglin; Ye, Yan; Liu, Yanhua; Chen, Linsen
2016-03-21
Without any special glasses, multiview 3D displays based on the diffractive optics can present high resolution, full-parallax 3D images in an ultra-wide viewing angle. The enabling optical component, namely the phase plate, can produce arbitrarily distributed view zones by carefully designing the orientation and the period of each nano-grating pixel. However, such 3D display screen is restricted to a limited size due to the time-consuming fabricating process of nano-gratings on the phase plate. In this paper, we proposed and developed a lithography system that can fabricate the phase plate efficiently. Here we made two phase plates with full nano-grating pixel coverage at a speed of 20 mm^{2}/mins, a 500 fold increment in the efficiency when compared to the method of E-beam lithography. One 2.5-inch phase plate generated 9-view 3D images with horizontal-parallax, while the other 6-inch phase plate produced 64-view 3D images with full-parallax. The angular divergence in horizontal axis and vertical axis was 1.5 degrees, and 1.25 degrees, respectively, slightly larger than the simulated value of 1.2 degrees by Finite Difference Time Domain (FDTD). The intensity variation was less than 10% for each viewpoint, in consistency with the simulation results. On top of each phase plate, a high-resolution binary masking pattern containing amplitude information of all viewing zone was well aligned. We achieved a resolution of 400 pixels/inch and a viewing angle of 40 degrees for 9-view 3D images with horizontal parallax. In another prototype, the resolution of each view was 160 pixels/inch and the view angle was 50 degrees for 64-view 3D images with full parallax. As demonstrated in the experiments, the homemade lithography system provided the key fabricating technology for multiview 3D holographic display. PMID:27136814
Method for making a single-step etch mask for 3D monolithic nanostructures.
Grishina, D A; Harteveld, C A M; Woldering, L A; Vos, W L
2015-12-18
Current nanostructure fabrication by etching is usually limited to planar structures as they are defined by a planar mask. The realization of three-dimensional (3D) nanostructures by etching requires technologies beyond planar masks. We present a method for fabricating a 3D mask that allows one to etch three-dimensional monolithic nanostructures using only CMOS-compatible processes. The mask is written in a hard-mask layer that is deposited on two adjacent inclined surfaces of a Si wafer. By projecting in a single step two different 2D patterns within one 3D mask on the two inclined surfaces, the mutual alignment between the patterns is ensured. Thereby after the mask pattern is defined, the etching of deep pores in two oblique directions yields a three-dimensional structure in Si. As a proof of concept we demonstrate 3D mask fabrication for three-dimensional diamond-like photonic band gap crystals in silicon. The fabricated crystals reveal a broad stop gap in optical reflectivity measurements. We propose how 3D nanostructures with five different Bravais lattices can be realized, namely cubic, tetragonal, orthorhombic, monoclinic and hexagonal, and demonstrate a mask for a 3D hexagonal crystal. We also demonstrate the mask for a diamond-structure crystal with a 3D array of cavities. In general, the 2D patterns on the different surfaces can be completely independently structured and still be in perfect mutual alignment. Indeed, we observe an alignment accuracy of better than 3.0 nm between the 2D mask patterns on the inclined surfaces, which permits one to etch well-defined monolithic 3D nanostructures. PMID:26581317
Method for making a single-step etch mask for 3D monolithic nanostructures
NASA Astrophysics Data System (ADS)
Grishina, D. A.; Harteveld, C. A. M.; Woldering, L. A.; Vos, W. L.
2015-12-01
Current nanostructure fabrication by etching is usually limited to planar structures as they are defined by a planar mask. The realization of three-dimensional (3D) nanostructures by etching requires technologies beyond planar masks. We present a method for fabricating a 3D mask that allows one to etch three-dimensional monolithic nanostructures using only CMOS-compatible processes. The mask is written in a hard-mask layer that is deposited on two adjacent inclined surfaces of a Si wafer. By projecting in a single step two different 2D patterns within one 3D mask on the two inclined surfaces, the mutual alignment between the patterns is ensured. Thereby after the mask pattern is defined, the etching of deep pores in two oblique directions yields a three-dimensional structure in Si. As a proof of concept we demonstrate 3D mask fabrication for three-dimensional diamond-like photonic band gap crystals in silicon. The fabricated crystals reveal a broad stop gap in optical reflectivity measurements. We propose how 3D nanostructures with five different Bravais lattices can be realized, namely cubic, tetragonal, orthorhombic, monoclinic and hexagonal, and demonstrate a mask for a 3D hexagonal crystal. We also demonstrate the mask for a diamond-structure crystal with a 3D array of cavities. In general, the 2D patterns on the different surfaces can be completely independently structured and still be in perfect mutual alignment. Indeed, we observe an alignment accuracy of better than 3.0 nm between the 2D mask patterns on the inclined surfaces, which permits one to etch well-defined monolithic 3D nanostructures.
2D-3D hybrid stabilized finite element method for tsunami runup simulations
NASA Astrophysics Data System (ADS)
Takase, S.; Moriguchi, S.; Terada, K.; Kato, J.; Kyoya, T.; Kashiyama, K.; Kotani, T.
2016-05-01
This paper presents a two-dimensional (2D)-three-dimensional (3D) hybrid stabilized finite element method that enables us to predict a propagation process of tsunami generated in a hypocentral region, which ranges from offshore propagation to runup to urban areas, with high accuracy and relatively low computational costs. To be more specific, the 2D shallow water equation is employed to simulate the propagation of offshore waves, while the 3D Navier-Stokes equation is employed for the runup in urban areas. The stabilized finite element method is utilized for numerical simulations for both of the 2D and 3D domains that are independently discretized with unstructured meshes. The multi-point constraint and transmission methods are applied to satisfy the continuity of flow velocities and pressures at the interface between the resulting 2D and 3D meshes, since neither their spatial dimensions nor node arrangements are consistent. Numerical examples are presented to demonstrate the performance of the proposed hybrid method to simulate tsunami behavior, including offshore propagation and runup to urban areas, with substantially lower computation costs in comparison with full 3D computations.
2D-3D hybrid stabilized finite element method for tsunami runup simulations
NASA Astrophysics Data System (ADS)
Takase, S.; Moriguchi, S.; Terada, K.; Kato, J.; Kyoya, T.; Kashiyama, K.; Kotani, T.
2016-09-01
This paper presents a two-dimensional (2D)-three-dimensional (3D) hybrid stabilized finite element method that enables us to predict a propagation process of tsunami generated in a hypocentral region, which ranges from offshore propagation to runup to urban areas, with high accuracy and relatively low computational costs. To be more specific, the 2D shallow water equation is employed to simulate the propagation of offshore waves, while the 3D Navier-Stokes equation is employed for the runup in urban areas. The stabilized finite element method is utilized for numerical simulations for both of the 2D and 3D domains that are independently discretized with unstructured meshes. The multi-point constraint and transmission methods are applied to satisfy the continuity of flow velocities and pressures at the interface between the resulting 2D and 3D meshes, since neither their spatial dimensions nor node arrangements are consistent. Numerical examples are presented to demonstrate the performance of the proposed hybrid method to simulate tsunami behavior, including offshore propagation and runup to urban areas, with substantially lower computation costs in comparison with full 3D computations.
NASA Technical Reports Server (NTRS)
Zhao, W.; Newman, J. C., Jr.; Sutton, M. A.; Wu, X. R.; Shivakumar, K. N.
1995-01-01
Stress intensity factors for quarter-elliptical corner cracks emanating from a circular hole are determined using a 3-D weight function method combined with a 3-D finite element method. The 3-D finite element method is used to analyze uncracked configuration and provide stress distribution in the region where crack is to occur. Using this stress distribution as input, the 3-D weight function method is used to determine stress intensity factors. Three different loading conditions, i.e. remote tension, remote bending and wedge loading, are considered for a wide range in geometrical parameters. The significance in using 3-D uncracked stress distribution and the difference between single and double corner cracks are studied. Typical crack opening displacements are also provided. Comparisons are made with solutions available in the literature.
NASA Astrophysics Data System (ADS)
Lamorski, Krzysztof; Sławiński, Cezary; Barna, Gyöngyi
2014-05-01
There are some important macroscopic properties of the soil porous media such as: saturated permeability and water retention characteristics. These soil characteristics are very important as they determine soil transport processes and are commonly used as a parameters of general models of soil transport processes used extensively for scientific developments and engineering practise. These characteristics are usually measured or estimated using some statistical or phenomenological modelling, i.e. pedotransfer functions. On the physical basis, saturated soil permeability arises from physical transport processes occurring at the pore level. Current progress in modelling techniques, computational methods and X-ray micro-tomographic technology gives opportunity to use direct methods of physical modelling for pore level transport processes. Physically valid description of transport processes at micro-scale based on Navier-Stokes type modelling approach gives chance to recover macroscopic porous medium characteristics from micro-flow modelling. Water microflow transport processes occurring at the pore level are dependent on the microstructure of porous body and interactions between the fluid and the medium. In case of soils, i.e. the medium there exist relatively big pores in which water can move easily but also finer pores are present in which water transport processes are dominated by strong interactions between the medium and the fluid - full physical description of these phenomena is a challenge. Ten samples of different soils were scanned using X-ray computational microtomograph. The diameter of samples was 5 mm. The voxel resolution of CT scan was 2.5 µm. Resulting 3D soil samples images were used for reconstruction of the pore space for further modelling. 3D image threshholding was made to determine the soil grain surface. This surface was triangulated and used for computational mesh construction for the pore space. Numerical modelling of water flow through the
Simulation of light transport in scintillators based on 3D characterization of crystal surfaces
Cherry, Simon R.
2013-01-01
In the development of positron emission tomography (PET) detectors, understanding and optimizing scintillator light collection is critical for achieving high performance, particularly when the design incorporates depth-of-interaction (DOI) encoding or time-of-flight information. Monte-Carlo simulations play an important role in guiding research in detector designs and popular software such as GATE now include models of light transport in scintillators. Although current simulation toolkits are able to provide accurate models of perfectly polished surfaces, they do not successfully predict light output for other surface finishes, for example those often used in DOI-encoding detectors. The lack of accuracy of those models mainly originates from a simplified description of rough surfaces as an ensemble of micro-facets determined by the distribution of their normal, typically a Gaussian distribution. The user can specify the standard deviation of this distribution, but this parameter does not provide a full description of the surface reflectance properties. We propose a different approach based on 3D measurements of the surface using atomic force microscopy (AFM). Polished and rough (unpolished) crystals were scanned to compute the surface reflectance properties. The angular distributions of reflectance and reflected rays were computed and stored in look-up tables (LUTs). The LUTs account for the effect of incidence angle and were integrated in a light transport model. Crystals of different sizes were simulated with and without reflector. The simulated maximum light output and the light output as a function of DOI showed very good agreement with experimental characterization of the crystals, indicating that our approach provides an accurate model of polished and rough surfaces and could be used to predict light collection in scintillators. This model is based on a true 3D representation of the surface, makes no assumption about the surface and provides insight on the
Detecting and estimating errors in 3D restoration methods using analog models.
NASA Astrophysics Data System (ADS)
José Ramón, Ma; Pueyo, Emilio L.; Briz, José Luis
2015-04-01
Some geological scenarios may be important for a number of socio-economic reasons, such as water or energy resources, but the available underground information is often limited, scarce and heterogeneous. A truly 3D reconstruction, which is still necessary during the decision-making process, may have important social and economic implications. For this reason, restoration methods were developed. By honoring some geometric or mechanical laws, they help build a reliable image of the subsurface. Pioneer methods were firstly applied in 2D (balanced and restored cross-sections) during the sixties and seventies. Later on, and due to the improvements of computational capabilities, they were extended to 3D. Currently, there are some academic and commercial restoration solutions; Unfold by the Université de Grenoble, Move by Midland Valley Exploration, Kine3D (on gOcad code) by Paradigm, Dynel3D by igeoss-Schlumberger. We have developed our own restoration method, Pmag3Drest (IGME-Universidad de Zaragoza), which is designed to tackle complex geometrical scenarios using paleomagnetic vectors as a pseudo-3D indicator of deformation. However, all these methods have limitations based on the assumptions they need to establish. For this reason, detecting and estimating uncertainty in 3D restoration methods is of key importance to trust the reconstructions. Checking the reliability and the internal consistency of every method, as well as to compare the results among restoration tools, is a critical issue never tackled so far because of the impossibility to test out the results in Nature. To overcome this problem we have developed a technique using analog models. We built complex geometric models inspired in real cases of superposed and/or conical folding at laboratory scale. The stratigraphic volumes were modeled using EVA sheets (ethylene vinyl acetate). Their rheology (tensile and tear strength, elongation, density etc) and thickness can be chosen among a large number of values
A Multiscale Constraints Method Localization of 3D Facial Feature Points
Li, Hong-an; Zhang, Yongxin; Li, Zhanli; Li, Huilin
2015-01-01
It is an important task to locate facial feature points due to the widespread application of 3D human face models in medical fields. In this paper, we propose a 3D facial feature point localization method that combines the relative angle histograms with multiscale constraints. Firstly, the relative angle histogram of each vertex in a 3D point distribution model is calculated; then the cluster set of the facial feature points is determined using the cluster algorithm. Finally, the feature points are located precisely according to multiscale integral features. The experimental results show that the feature point localization accuracy of this algorithm is better than that of the localization method using the relative angle histograms. PMID:26539244
A simple method for the production of anti-C3d monoclonal antibody.
Cruz, Carlos; León, Graciela
2007-12-01
Production of monoclonal antibodies to C3d usually involves the purification of protein. Our method does not require C3 purification; it relies on attachment of C3b to mouse erythrocytes by activation of alternative pathways and further conversion in C3d. We prepared human complement-coated mouse red cells and sensitized mice of the same strain with our own schedule of immunization and applied the classical methods to obtain a mouse monoclonal antibody. We obtained a clone called BMS-11 which produces a monoclonal antibody of IgM class, to C3d with a title of 1:500000. The monoclonal antibody obtained has shown that it is suitable for use as an antiglobulin reagent. PMID:18158789
Simulations of Coulomb systems with slab geometry using an efficient 3D Ewald summation method
NASA Astrophysics Data System (ADS)
dos Santos, Alexandre P.; Girotto, Matheus; Levin, Yan
2016-04-01
We present a new approach to efficiently simulate electrolytes confined between infinite charged walls using a 3d Ewald summation method. The optimal performance is achieved by separating the electrostatic potential produced by the charged walls from the electrostatic potential of electrolyte. The electric field produced by the 3d periodic images of the walls is constant inside the simulation cell, with the field produced by the transverse images of the charged plates canceling out. The non-neutral confined electrolyte in an external potential can be simulated using 3d Ewald summation with a suitable renormalization of the electrostatic energy, to remove a divergence, and a correction that accounts for the conditional convergence of the resulting lattice sum. The new algorithm is at least an order of magnitude more rapid than the usual simulation methods for the slab geometry and can be further sped up by adopting a particle-particle particle-mesh approach.
A Review of Failure Analysis Methods for Advanced 3D Microelectronic Packages
NASA Astrophysics Data System (ADS)
Li, Yan; Srinath, Purushotham Kaushik Muthur; Goyal, Deepak
2016-01-01
Advanced three dimensional (3D) packaging is a key enabler in driving form factor reduction, performance benefits, and package cost reduction, especially in the fast paced mobility and ultraportable consumer electronics segments. The high level of functional integration and the complex package architecture pose a significant challenge for conventional fault isolation (FI) and failure analysis (FA) methods. Innovative FI/FA tools and techniques are required to tackle the technical and throughput challenges. In this paper, the applications of FI and FA techniques such as Electro Optic Terahertz Pulse Reflectometry, 3D x-ray computed tomography, lock-in thermography, and novel physical sample preparation methods to 3D packages with package on package and stacked die with through silicon via configurations are reviewed, along with the key FI and FA challenges.
Uncovering the true nature of deformation microstructures using 3D analysis methods
NASA Astrophysics Data System (ADS)
Ferry, M.; Quadir, M. Z.; Afrin, N.; Xu, W.; Loeb, A.; Soe, B.; McMahon, C.; George, C.; Bassman, L.
2015-08-01
Three-dimensional electron backscatter diffraction (3D EBSD) has emerged as a powerful technique for generating 3D crystallographic information in reasonably large volumes of a microstructure. The technique uses a focused ion beam (FIB) as a high precision serial sectioning device for generating consecutive ion milled surfaces of a material, with each milled surface subsequently mapped by EBSD. The successive EBSD maps are combined using a suitable post-processing method to generate a crystallographic volume of the microstructure. The first part of this paper shows the usefulness of 3D EBSD for understanding the origin of various structural features associated with the plastic deformation of metals. The second part describes a new method for automatically identifying the various types of low and high angle boundaries found in deformed and annealed metals, particularly those associated with grains exhibiting subtle and gradual variations in orientation. We have adapted a 2D image segmentation technique, fast multiscale clustering, to 3D EBSD data using a novel variance function to accommodate quaternion data. This adaptation is capable of segmenting based on subtle and gradual variation as well as on sharp boundaries within the data. We demonstrate the excellent capabilities of this technique with application to 3D EBSD data sets generated from a range of cold rolled and annealed metals described in the paper.
Device and methods for "gold standard" registration of clinical 3D and 2D cerebral angiograms
NASA Astrophysics Data System (ADS)
Madan, Hennadii; Likar, Boštjan; Pernuš, Franjo; Å piclin, Žiga
2015-03-01
Translation of any novel and existing 3D-2D image registration methods into clinical image-guidance systems is limited due to lack of their objective validation on clinical image datasets. The main reason is that, besides the calibration of the 2D imaging system, a reference or "gold standard" registration is very difficult to obtain on clinical image datasets. In the context of cerebral endovascular image-guided interventions (EIGIs), we present a calibration device in the form of a headband with integrated fiducial markers and, secondly, propose an automated pipeline comprising 3D and 2D image processing, analysis and annotation steps, the result of which is a retrospective calibration of the 2D imaging system and an optimal, i.e., "gold standard" registration of 3D and 2D images. The device and methods were used to create the "gold standard" on 15 datasets of 3D and 2D cerebral angiograms, whereas each dataset was acquired on a patient undergoing EIGI for either aneurysm coiling or embolization of arteriovenous malformation. The use of the device integrated seamlessly in the clinical workflow of EIGI. While the automated pipeline eliminated all manual input or interactive image processing, analysis or annotation. In this way, the time to obtain the "gold standard" was reduced from 30 to less than one minute and the "gold standard" of 3D-2D registration on all 15 datasets of cerebral angiograms was obtained with a sub-0.1 mm accuracy.
Methods of constructing a 3D geological model from scatter data
Horsman, J.; Bethel, W.
1995-04-01
Most geoscience applications, such as assessment of an oil reservoir or hazardous waste site, require geological characterization of the site. Geological characterization involves analysis of spatial distributions of lithology, porosity, etc. Because of the complexity of the spatial relationships, the authors find that a 3-D model of geology is better suited for integration of many different types of data and provides a better representation of a site than a 2-D one. A 3-D model of geology is constructed from sample data obtained from field measurements, which are usually scattered. To create a volume model from scattered data, interpolation between points is required. The interpolation can be computed using one of several computational algorithms. Alternatively, a manual method may be employed, in which an interactive graphics device is used to input by hand the information that lies between the data points. For example, a mouse can be used to draw lines connecting data points with equal values. The combination of these two methods presents yet another approach. In this study, the authors will compare selected methods of 3-D geological modeling, They used a flow-based, modular visualization environment (AVS) to construct the geological models computationally. Within this system, they used three modules, scat{_}3d, trivar and scatter{_}to{_}ucd, as examples of computational methods. They compare these methods to the combined manual and computational approach. Because there are no tools readily available in AVS for this type of construction, they used a geological modeling system to demonstrate this method.
Flatbed-type 3D display systems using integral imaging method
NASA Astrophysics Data System (ADS)
Hirayama, Yuzo; Nagatani, Hiroyuki; Saishu, Tatsuo; Fukushima, Rieko; Taira, Kazuki
2006-10-01
We have developed prototypes of flatbed-type autostereoscopic display systems using one-dimensional integral imaging method. The integral imaging system reproduces light beams similar of those produced by a real object. Our display architecture is suitable for flatbed configurations because it has a large margin for viewing distance and angle and has continuous motion parallax. We have applied our technology to 15.4-inch displays. We realized horizontal resolution of 480 with 12 parallaxes due to adoption of mosaic pixel arrangement of the display panel. It allows viewers to see high quality autostereoscopic images. Viewing the display from angle allows the viewer to experience 3-D images that stand out several centimeters from the surface of the display. Mixed reality of virtual 3-D objects and real objects are also realized on a flatbed display. In seeking reproduction of natural 3-D images on the flatbed display, we developed proprietary software. The fast playback of the CG movie contents and real-time interaction are realized with the aid of a graphics card. Realization of the safety 3-D images to the human beings is very important. Therefore, we have measured the effects on the visual function and evaluated the biological effects. For example, the accommodation and convergence were measured at the same time. The various biological effects are also measured before and after the task of watching 3-D images. We have found that our displays show better results than those to a conventional stereoscopic display. The new technology opens up new areas of application for 3-D displays, including arcade games, e-learning, simulations of buildings and landscapes, and even 3-D menus in restaurants.
NASA Astrophysics Data System (ADS)
Crockett, Ethan Van
The need for clinically intuitive metrics for patient-specific quality assurance in radiation therapy has been well-documented (Zhen, Nelms et al. 2011). A novel transform method has shown to be effective at converting full-density 3D dose measurements made in a phantom to dose values in the patient geometry, enabling comparisons using clinically intuitive metrics such as dose-volume histograms (Oldham et al. 2011). This work investigates the transform method and compares its calculated dose-volume histograms (DVHs) to DVH values calculated by a Delta4 QA device (Scandidos), marking the first comparison of a true 3D system to a semi-3D device using clinical metrics. Measurements were made using Presage 3D dosimeters, which were readout by an in-house optical-CT scanner. Three patient cases were chosen for the study: one head-and-neck VMAT treatment and two spine IMRT treatments. The transform method showed good agreement with the planned dose values for all three cases. Furthermore, the transformed DVHs adhered to the planned dose with more accuracy than the Delta4 DVHs. The similarity between the Delta4 DVHs and the transformed DVHs, however, was greater for one of the spine cases than it was for the head-and-neck case, implying that the accuracy of the Delta4 Anatomy software may vary from one treatment site to another. Overall, the transform method, which incorporates data from full-density 3D dose measurements, provides clinically intuitive results that are more accurate and consistent than the corresponding results from a semi-3D Delta 4 system.
NASA Astrophysics Data System (ADS)
Oh, J.; Jiang, X.; Waliser, D. E.; Moncrieff, M. W.; Johnson, R. H.
2013-12-01
As one of the most prominent tropical atmospheric variability modes, the Madden-Julian Oscillation (MJO) exerts profound influences on global weather and climate, and serves as a critical predictability source for extend-range forecast. While credible representation of the MJO still represents a great challenge for current general circulation models (GCMs), previous studies on the vertical structure of the MJO have largely focused on collective impacts from multi-scale convective systems on thermodynamic properties of the MJO. Most recently, limited observational studies and idealized modeling work suggested that convective momentum transport (CMT) could also play an important role in interpreting the observed MJO features. In this study, the 3D CMT structure associated with the MJO is examined by analyzing model output from three recent high-quality reanalysis systems, including NOAA's Climate Forecast System Reanalysis (CFSR), NASA's Modern Era Retrospective-analysis for Research and Applications (MERRA), and ECMWF-the Year of Tropical Convection (YOTC) reanalysis. Consistent with previous cloud-resolving model study, a well-organized three-layer vertical structure in the CMT associated with the MJO is also discerned based on reanalyses. The result suggests that CMT tends to intensify the MJO circulation, particularly in the lower troposphere. Relative roles of meso-scale systems (MCS) and synoptic waves in contributing the total CMT profiles of the MJO will also be explored. Differences in CMT profiles in these several reanalysis models will be discussed.
3D imaging of amplitude objects embedded in phase objects using transport of intensity
NASA Astrophysics Data System (ADS)
Banerjee, Partha; Basunia, Mahmudunnabi
2015-09-01
The amplitude and phase of the complex optical field in the Helmholtz equation obey a pair of coupled equations, arising from equating the real and imaginary parts. The imaginary part yields the transport of intensity equation (TIE), which can be used to derive the phase distribution at the observation plane. If a phase object is approximately imaged on the recording plane(s), TIE yields the phase without the need for phase unwrapping. In our experiment, the 3D image of a phase object and an amplitude object embedded in a phase object is recovered. The phase object is created by heating a liquid, comprising a solution of red dye in alcohol, using a focused 514 nm laser beam to the point where self-phase modulation of the beam is observed. The optical intensities are recorded at various planes during propagation of a low power 633 nm laser beam through the liquid. In the process of applying TIE to derive the phase at the observation plane, the real part of the complex equation is also examined as a cross-check of our calculations. For pure phase objects, it is shown that the real part of the complex equation is best satisfied around the image plane. Alternatively, it is proposed that this information can be used to determine the optimum image plane.
3-D simulation of gases transport under condition of inert gas injection into goaf
NASA Astrophysics Data System (ADS)
Liu, Mao-Xi; Shi, Guo-Qing; Guo, Zhixiong; Wang, Yan-Ming; Ma, Li-Yang
2016-02-01
To prevent coal spontaneous combustion in mines, it is paramount to understand O2 gas distribution under condition of inert gas injection into goaf. In this study, the goaf was modeled as a 3-D porous medium based on stress distribution. The variation of O2 distribution influenced by CO2 or N2 injection was simulated based on the multi-component gases transport and the Navier-Stokes equations using Fluent. The numerical results without inert gas injection were compared with field measurements to validate the simulation model. Simulations with inert gas injection show that CO2 gas mainly accumulates at the goaf floor level; however, a notable portion of N2 gas moves upward. The evolution of the spontaneous combustion risky zone with continuous inert gas injection can be classified into three phases: slow inerting phase, rapid accelerating inerting phase, and stable inerting phase. The asphyxia zone with CO2 injection is about 1.25-2.4 times larger than that with N2 injection. The efficacy of preventing and putting out mine fires is strongly related with the inert gas injecting position. Ideal injections are located in the oxidation zone or the transitional zone between oxidation zone and heat dissipation zone.
A 3-D aerodynamic method for the analysis of isolated horizontal-axis wind turbines
Ammara, I.; Masson, C.; Paraschivoiu, I.
1997-12-31
In most existing performance-analysis methods, wind turbines are considered isolated so that interference effects caused by other rotors or by the site topography are neglected. The main objective of this paper is to propose a practical 3-D method suitable for the study of these effects, in order to optimize the arrangement and the positioning of Horizontal-Axis Wind Turbines (HAWTs) in a wind farm. In the proposed methodology, the flow field around isolated HAWTs is predicted by solving the 3-D, time-averaged, steady-state, incompressible, Navier-Stokes equations in which the turbines are represented by distributions of momentum sources. The resulting governing equations are solved using a Control-Volume Finite Element Method (CVFEM). The fundamental aspects related to the development of a practical 3-D method are discussed in this paper, with an emphasis on some of the challenges that arose during its implementation. The current implementation is limited to the analysis of isolated HAWTs. Preliminary results have indicated that, the proposed 3-D method reaches the same level of accuracy, in terms of performance predictions, that the previously developed 2-D axisymmetric model and the well-known momentum-strip theory, while still using reasonable computers resources. It can be considered as a useful tool for the design of HAWTs. Its main advantages, however, are its intrinsic capacity to predict the details of the flow in the wake, and its capabilities of modelling arbitrary wind-turbine arrangements and including ground effects.
High fidelity digital inline holographic method for 3D flow measurements.
Toloui, Mostafa; Hong, Jiarong
2015-10-19
Among all the 3D optical flow diagnostic techniques, digital inline holographic particle tracking velocimetry (DIH-PTV) provides the highest spatial resolution with low cost, simple and compact optical setups. Despite these advantages, DIH-PTV suffers from major limitations including poor longitudinal resolution, human intervention (i.e. requirement for manually determined tuning parameters during tracer field reconstruction and extraction), limited tracer concentration, and expensive computations. These limitations prevent this technique from being widely used for high resolution 3D flow measurements. In this study, we present a novel holographic particle extraction method with the goal of overcoming all the major limitations of DIH-PTV. The proposed method consists of multiple steps involving 3D deconvolution, automatic signal-to-noise ratio enhancement and thresholding, and inverse iterative particle extraction. The entire method is implemented using GPU-based algorithm to increase the computational speed significantly. Validated with synthetic particle holograms, the proposed method can achieve particle extraction rate above 95% with fake particles less than 3% and maximum position error below 1.6 particle diameter for holograms with particle concentration above 3000 particles/mm^{3}. The applicability of the proposed method for DIH-PTV has been further validated using the experiment of laminar flow in a microchannel and the synthetic tracer flow fields generated using a DNS turbulent channel flow database. Such improvements will substantially enhance the implementation of DIH-PTV for 3D flow measurements and enable the potential commercialization of this technique. PMID:26480377
Accident or homicide--virtual crime scene reconstruction using 3D methods.
Buck, Ursula; Naether, Silvio; Räss, Beat; Jackowski, Christian; Thali, Michael J
2013-02-10
The analysis and reconstruction of forensically relevant events, such as traffic accidents, criminal assaults and homicides are based on external and internal morphological findings of the injured or deceased person. For this approach high-tech methods are gaining increasing importance in forensic investigations. The non-contact optical 3D digitising system GOM ATOS is applied as a suitable tool for whole body surface and wound documentation and analysis in order to identify injury-causing instruments and to reconstruct the course of event. In addition to the surface documentation, cross-sectional imaging methods deliver medical internal findings of the body. These 3D data are fused into a whole body model of the deceased. Additional to the findings of the bodies, the injury inflicting instruments and incident scene is documented in 3D. The 3D data of the incident scene, generated by 3D laser scanning and photogrammetry, is also included into the reconstruction. Two cases illustrate the methods. In the fist case a man was shot in his bedroom and the main question was, if the offender shot the man intentionally or accidentally, as he declared. In the second case a woman was hit by a car, driving backwards into a garage. It was unclear if the driver drove backwards once or twice, which would indicate that he willingly injured and killed the woman. With this work, we demonstrate how 3D documentation, data merging and animation enable to answer reconstructive questions regarding the dynamic development of patterned injuries, and how this leads to a real data based reconstruction of the course of event. PMID:22727689
A fast and accurate method to predict 2D and 3D aerodynamic boundary layer flows
NASA Astrophysics Data System (ADS)
Bijleveld, H. A.; Veldman, A. E. P.
2014-12-01
A quasi-simultaneous interaction method is applied to predict 2D and 3D aerodynamic flows. This method is suitable for offshore wind turbine design software as it is a very accurate and computationally reasonably cheap method. This study shows the results for a NACA 0012 airfoil. The two applied solvers converge to the experimental values when the grid is refined. We also show that in separation the eigenvalues remain positive thus avoiding the Goldstein singularity at separation. In 3D we show a flow over a dent in which separation occurs. A rotating flat plat is used to show the applicability of the method for rotating flows. The shown capabilities of the method indicate that the quasi-simultaneous interaction method is suitable for design methods for offshore wind turbine blades.
Estimating the complexity of 3D structural models using machine learning methods
NASA Astrophysics Data System (ADS)
Mejía-Herrera, Pablo; Kakurina, Maria; Royer, Jean-Jacques
2016-04-01
Quantifying the complexity of 3D geological structural models can play a major role in natural resources exploration surveys, for predicting environmental hazards or for forecasting fossil resources. This paper proposes a structural complexity index which can be used to help in defining the degree of effort necessary to build a 3D model for a given degree of confidence, and also to identify locations where addition efforts are required to meet a given acceptable risk of uncertainty. In this work, it is considered that the structural complexity index can be estimated using machine learning methods on raw geo-data. More precisely, the metrics for measuring the complexity can be approximated as the difficulty degree associated to the prediction of the geological objects distribution calculated based on partial information on the actual structural distribution of materials. The proposed methodology is tested on a set of 3D synthetic structural models for which the degree of effort during their building is assessed using various parameters (such as number of faults, number of part in a surface object, number of borders, ...), the rank of geological elements contained in each model, and, finally, their level of deformation (folding and faulting). The results show how the estimated complexity in a 3D model can be approximated by the quantity of partial data necessaries to simulated at a given precision the actual 3D model without error using machine learning algorithms.
NASA Astrophysics Data System (ADS)
Toschi, I.; Capra, A.; De Luca, L.; Beraldin, J.-A.; Cournoyer, L.
2014-05-01
This paper discusses a methodology to evaluate the accuracy of recently developed image-based 3D modelling techniques. So far, the emergence of these novel methods has not been supported by the definition of an internationally recognized standard which is fundamental for user confidence and market growth. In order to provide an element of reflection and solution to the different communities involved in 3D imaging, a promising approach is presented in this paper for the assessment of both metric quality and limitations of an open-source suite of tools (Apero/MicMac), developed for the extraction of dense 3D point clouds from a set of unordered 2D images. The proposed procedural workflow is performed within a metrological context, through inter-comparisons with "reference" data acquired with two hemispherical laser scanners, one total station, and one laser tracker. The methodology is applied to two case studies, designed in order to analyse the software performances in dealing with both outdoor and environmentally controlled conditions, i.e. the main entrance of Cathédrale de la Major (Marseille, France) and a custom-made scene located at National Research Council of Canada 3D imaging Metrology Laboratory (Ottawa). Comparative data and accuracy evidence produced for both tests allow the study of some key factors affecting 3D model accuracy.
NASA Astrophysics Data System (ADS)
Bhavanam, Sharada
The aim of this thesis is to numerically evaluate the mixed-mode Stress Intensity Factors (SIFs) of complex 3D structural geometries with arbitrary 3D cracks using the Symmetric Galerkin Boundary Element Method-Finite Element Method (SGBEM-FEM) Alternating Method. Various structural geometries with different loading scenarios and crack configurations were examined in this thesis to understand the behavior and trends of the mixed-mode SIFs as well as the fatigue life for these complex structural geometries. Although some 3D structures have empirical and numerical solutions that are readily available in the open literature, some do not; therefore this thesis presents the results of fracture and fatigue analyses of these 3D complex structures using the SGBEM-FEM Alternating Method to serve as reference for future studies. Furthermore, there are advantages of using the SGBEM-FEM Alternating Method compared to traditional FEM methods. For example, the fatigue-crack-growth and fatigue life can be better estimated for a structure because different fatigue models (i.e. Walker, Paris, and NASGRO) can be used within the same framework of the SGBEM-FEM Alternating Method. The FEM (un-cracked structure)/BEM(crack model) meshes are modeled independently, which speeds up the computation process and reduces the cost of human labor. A simple coarse mesh can be used for all fracture and fatigue analyses of complex structures. In this thesis, simple coarse meshes were used for 3D complex structures, which were below 5000 elements as compared to traditional FEM, which require meshes where the elements range on the order of ˜250,000 to ˜106 and sometimes even more than that.
2D and 3D Method of Characteristic Tools for Complex Nozzle Development
NASA Technical Reports Server (NTRS)
Rice, Tharen
2003-01-01
This report details the development of a 2D and 3D Method of Characteristic (MOC) tool for the design of complex nozzle geometries. These tools are GUI driven and can be run on most Windows-based platforms. The report provides a user's manual for these tools as well as explains the mathematical algorithms used in the MOC solutions.
Simulation of surface tension in 2D and 3D with smoothed particle hydrodynamics method
NASA Astrophysics Data System (ADS)
Zhang, Mingyu
2010-09-01
The methods for simulating surface tension with smoothed particle hydrodynamics (SPH) method in two dimensions and three dimensions are developed. In 2D surface tension model, the SPH particle on the boundary in 2D is detected dynamically according to the algorithm developed by Dilts [G.A. Dilts, Moving least-squares particle hydrodynamics II: conservation and boundaries, International Journal for Numerical Methods in Engineering 48 (2000) 1503-1524]. The boundary curve in 2D is reconstructed locally with Lagrangian interpolation polynomial. In 3D surface tension model, the SPH particle on the boundary in 3D is detected dynamically according to the algorithm developed by Haque and Dilts [A. Haque, G.A. Dilts, Three-dimensional boundary detection for particle methods, Journal of Computational Physics 226 (2007) 1710-1730]. The boundary surface in 3D is reconstructed locally with moving least squares (MLS) method. By transforming the coordinate system, it is guaranteed that the interface function is one-valued in the local coordinate system. The normal vector and curvature of the boundary surface are calculated according to the reconstructed boundary surface and then surface tension force can be calculated. Surface tension force acts only on the boundary particle. Density correction is applied to the boundary particle in order to remove the boundary inconsistency. The surface tension models in 2D and 3D have been applied to benchmark tests for surface tension. The ability of the current method applying to the simulation of surface tension in 2D and 3D is proved.
3D Spectral Element Method Simulations Of The Seismic Response of Caracas (Venezuela) Basin
NASA Astrophysics Data System (ADS)
Delavaud, E.; Vilotte, J.; Festa, G.; Cupillard, P.
2007-12-01
We present here 3D numerical simulations of the response of the Caracas (Venezuela) valley up to 5 Hz for different scenarios of plane wave excitation based on the regional seismicity. Attention is focused on the effects of the 3D basin geometry and of the adjacent regional topography. The simulations are performed using Spectral Element method (SEM) together with an unstructured hexahedral mesh discretization and perfectly matched layers (PML). These simulations show 3D amplification phenomena associated with complex wave reflexion, diffraction and focalisation patterns linked to the geometry of the basin. Time and frequency analysis reveal some interesting features both in terms of amplification and energy residence in the basin. The low frequency amplification pattern is mainly controlled by the early response of the basin to the incident plane wave while the high frequency amplification patterns result mainly from late arrivals where complex 3D wave diffraction phenomena are dominating and the memory of the initial excitation is lost. Interestingly enough, it is shown that H/V method correctly predict the low frequency amplification pattern when apply to the late part of the recorded seismograms. The complex high frequency amplification pattern is shown to be associated with surface wave generation at, and propagation from, sharp edges of the basin. Importance of 3D phenomena is assessed by comparison with simple 2D simulations. Significant differences in terms of time of residence, energy and amplification levels point out the interest of complete 3D modeling. In conclusions some of the limitations associated with the use of unstructured hexahedral meshes will be adressed. Despite the use of unstructured meshing tool, modeling the geometry of geological basins remain a complex and time consuming task. Possible extensions using more elaborate techniques like non conforming domain decomposition will be also discussed in conclusion.
3D electro-thermal Monte Carlo study of transport in confined silicon devices
NASA Astrophysics Data System (ADS)
Mohamed, Mohamed Y.
The simultaneous explosion of portable microelectronics devices and the rapid shrinking of microprocessor size have provided a tremendous motivation to scientists and engineers to continue the down-scaling of these devices. For several decades, innovations have allowed components such as transistors to be physically reduced in size, allowing the famous Moore's law to hold true. As these transistors approach the atomic scale, however, further reduction becomes less probable and practical. As new technologies overcome these limitations, they face new, unexpected problems, including the ability to accurately simulate and predict the behavior of these devices, and to manage the heat they generate. This work uses a 3D Monte Carlo (MC) simulator to investigate the electro-thermal behavior of quasi-one-dimensional electron gas (1DEG) multigate MOSFETs. In order to study these highly confined architectures, the inclusion of quantum correction becomes essential. To better capture the influence of carrier confinement, the electrostatically quantum-corrected full-band MC model has the added feature of being able to incorporate subband scattering. The scattering rate selection introduces quantum correction into carrier movement. In addition to the quantum effects, scaling introduces thermal management issues due to the surge in power dissipation. Solving these problems will continue to bring improvements in battery life, performance, and size constraints of future devices. We have coupled our electron transport Monte Carlo simulation to Aksamija's phonon transport so that we may accurately and efficiently study carrier transport, heat generation, and other effects at the transistor level. This coupling utilizes anharmonic phonon decay and temperature dependent scattering rates. One immediate advantage of our coupled electro-thermal Monte Carlo simulator is its ability to provide an accurate description of the spatial variation of self-heating and its effect on non
A new method of 3D scene recognition from still images
NASA Astrophysics Data System (ADS)
Zheng, Li-ming; Wang, Xing-song
2014-04-01
Most methods of monocular visual three dimensional (3D) scene recognition involve supervised machine learning. However, these methods often rely on prior knowledge. Specifically, they learn the image scene as part of a training dataset. For this reason, when the sampling equipment or scene is changed, monocular visual 3D scene recognition may fail. To cope with this problem, a new method of unsupervised learning for monocular visual 3D scene recognition is here proposed. First, the image is made using superpixel segmentation based on the CIELAB color space values L, a, and b and on the coordinate values x and y of pixels, forming a superpixel image with a specific density. Second, a spectral clustering algorithm based on the superpixels' color characteristics and neighboring relationships was used to reduce the dimensions of the superpixel image. Third, the fuzzy distribution density functions representing sky, ground, and façade are multiplied with the segment pixels, where the expectations of these segments are obtained. A preliminary classification of sky, ground, and façade is generated in this way. Fourth, the most accurate classification images of sky, ground, and façade were extracted through the tier-1 wavelet sampling and Manhattan direction feature. Finally, a depth perception map is generated based on the pinhole imaging model and the linear perspective information of ground surface. Here, 400 images of Make3D Image data from the Cornell University website were used to test the algorithm. The experimental results showed that this unsupervised learning method provides a more effective monocular visual 3D scene recognition model than other methods.
NASA Astrophysics Data System (ADS)
Wen, D. H.; Jiang, Y. S.; Zhang, Y. Z.; Gao, Q.
2014-03-01
The theoretical and experimental investigations on the polarization imagery system of speckle statistical characteristics and speckle removing method are researched. A method to obtain two images encoded by polarization degree with a single measurement process is proposed. A theoretical model for polarization imagery system on Müller matrix is proposed. According to modern charge coupled device (CCD) imaging characteristics, speckles are divided into two kinds, namely small speckle and big speckle. Based on this model, a speckle reduction algorithm based on a dual-tree complex wavelet transform (DTCWT) and blockmatching 3D filter (BM3D) is proposed (DTBM3D). Original laser image data transformed by logarithmic compression is decomposed by DTCWT into approximation and detail subbands. Bilateral filtering is applied to the approximation subbands, and a suited BM3D filter is applied to the detail subbands. The despeckling results show that contrast improvement index and edge preserve index outperform those of traditional methods. The researches have important reference value in research of speckle noise level and removing speckle noise.
Gap-filling methods for 3D PlanTIS data.
Loukiala, A; Tuna, U; Beer, S; Jahnke, S; Ruotsalainen, U
2010-10-21
The range of positron emitters and their labeled compounds have led to high-resolution PET scanners becoming widely used, not only in clinical and pre-clinical studies but also in plant studies. A high-resolution PET scanner, plant tomographic imaging system (PlanTIS), was designed to study metabolic and physiological functions of plants noninvasively. The gantry of the PlanTIS scanner has detector-free regions. Even when the gantry of the PlanTIS is rotated during the scan, these regions result in missing sinogram bins in the acquired data. Missing data need to be estimated prior to the analytical image reconstructions in order to avoid artifacts in the final reconstructed images. In this study, we propose three gap-filling methods for estimation of the unique gaps existing in the 3D PlanTIS sinogram data. The 3D sinogram data were gap-filled either by linear interpolation in the transaxial planes or by the bicubic interpolation method (proposed for the ECAT high-resolution research tomograph) in the transradial planes or by the inpainting method in the transangular planes. Each gap-filling method independently compensates for slices in one of three orthogonal sinogram planes (transaxial, transradial and transangular planes). A 3D numerical Shepp-Logan phantom and the NEMA image quality phantom were used to evaluate the methods. The gap-filled sinograms were reconstructed using the analytical 3D reprojection (3DRP) method. The NEMA phantom sinograms were also reconstructed by the iterative reconstruction method, ordered subsets maximum a posteriori one step late (OSMAPOSL), to compare the results of gap filling followed by 3DRP with the results of OSMAPOSL reconstruction without gap filling. The three methods were evaluated quantitatively (by mean square error and coefficients of variation) over the selected regions of the 3D numerical Shepp-Logan phantom at eight different Poisson noise levels. Moreover, the NEMA phantom scan data were used in visual assessments
Gap-filling methods for 3D PlanTIS data
NASA Astrophysics Data System (ADS)
Loukiala, A.; Tuna, U.; Beer, S.; Jahnke, S.; Ruotsalainen, U.
2010-10-01
The range of positron emitters and their labeled compounds have led to high-resolution PET scanners becoming widely used, not only in clinical and pre-clinical studies but also in plant studies. A high-resolution PET scanner, plant tomographic imaging system (PlanTIS), was designed to study metabolic and physiological functions of plants noninvasively. The gantry of the PlanTIS scanner has detector-free regions. Even when the gantry of the PlanTIS is rotated during the scan, these regions result in missing sinogram bins in the acquired data. Missing data need to be estimated prior to the analytical image reconstructions in order to avoid artifacts in the final reconstructed images. In this study, we propose three gap-filling methods for estimation of the unique gaps existing in the 3D PlanTIS sinogram data. The 3D sinogram data were gap-filled either by linear interpolation in the transaxial planes or by the bicubic interpolation method (proposed for the ECAT high-resolution research tomograph) in the transradial planes or by the inpainting method in the transangular planes. Each gap-filling method independently compensates for slices in one of three orthogonal sinogram planes (transaxial, transradial and transangular planes). A 3D numerical Shepp-Logan phantom and the NEMA image quality phantom were used to evaluate the methods. The gap-filled sinograms were reconstructed using the analytical 3D reprojection (3DRP) method. The NEMA phantom sinograms were also reconstructed by the iterative reconstruction method, ordered subsets maximum a posteriori one step late (OSMAPOSL), to compare the results of gap filling followed by 3DRP with the results of OSMAPOSL reconstruction without gap filling. The three methods were evaluated quantitatively (by mean square error and coefficients of variation) over the selected regions of the 3D numerical Shepp-Logan phantom at eight different Poisson noise levels. Moreover, the NEMA phantom scan data were used in visual assessments
Efficient calculation method for realistic deep 3D scene hologram using orthographic projection
NASA Astrophysics Data System (ADS)
Igarashi, Shunsuke; Nakamura, Tomoya; Matsushima, Kyoji; Yamaguchi, Masahiro
2016-03-01
We propose a fast calculation method to synthesize a computer-generated hologram (CGH) of realistic deep three-dimensional (3D) scene. In our previous study, we have proposed a calculation method of CGH for reproducing such scene called ray-sampling-plane (RSP) method, in which light-ray information of a scene is converted to wavefront, and the wavefront is numerically propagated based on diffraction theory. In this paper, we introduce orthographic projection to the RSP method for accelerating calculation time. By numerical experiments, we verified the accelerated calculation with the ratio of 28-times compared to the conventional RSP method. The calculated CGH was fabricated by the printing system using laser lithography and demonstrated deep 3D image reconstruction in 52mm×52mm with realistic appearance effect such as gloss and translucent effect.
Computational methods for constructing protein structure models from 3D electron microscopy maps
Esquivel-Rodríguez, Juan; Kihara, Daisuke
2013-01-01
Protein structure determination by cryo-electron microscopy (EM) has made significant progress in the past decades. Resolutions of EM maps have been improving as evidenced by recently reported structures that are solved at high resolutions close to 3 Å. Computational methods play a key role in interpreting EM data. Among many computational procedures applied to an EM map to obtain protein structure information, in this article we focus on reviewing computational methods that model protein three-dimensional (3D) structures from a 3D EM density map that is constructed from two-dimensional (2D) maps. The computational methods we discuss range from de novo methods, which identify structural elements in an EM map, to structure fitting methods, where known high resolution structures are fit into a low-resolution EM map. A list of available computational tools is also provided. PMID:23796504
A new 3D tracking method exploiting the capabilities of digital holography in microscopy
NASA Astrophysics Data System (ADS)
Miccio, L.; Memmolo, P.; Merola, F.; Fusco, S.; Embrione, V.; Netti, P. A.; Ferraro, P.
2013-04-01
A method for 3D tracking has been developed exploiting Digital Holographic Microscopy (DHM) features. In the framework of self-consistent platform for manipulation and measurement of biological specimen we use DHM for quantitative and completely label free analysis of specimen with low amplitude contrast. Tracking capability extend the potentiality of DHM allowing to monitor the motion of appropriate probes and correlate it with sample properties. Complete 3D tracking has been obtained for the probes avoiding the issue of amplitude refocusing in traditional tracking processing. Our technique belongs to the video tracking methods that, conversely from Quadrant Photo-Diode method, opens the possibility to track multiples probes. All the common used video tracking algorithms are based on the numerical analysis of amplitude images in the focus plane and the shift of the maxima in the image plane are measured after the application of an appropriate threshold. Our approach for video tracking uses different theoretical basis. A set of interferograms is recorded and the complex wavefields are managed numerically to obtain three dimensional displacements of the probes. The procedure works properly on an higher number of probes and independently from their size. This method overcomes the traditional video tracking issues as the inability to measure the axial movement and the choice of suitable threshold mask. The novel configuration allows 3D tracking of micro-particles and simultaneously can furnish Quantitative Phase-contrast maps of tracked micro-objects by interference microscopy, without changing the configuration. In this paper, we show a new concept for a compact interferometric microscope that can ensure the multifunctionality, accomplishing accurate 3D tracking and quantitative phase-contrast analysis. Experimental results are presented and discussed for in vitro cells. Through a very simple and compact optical arrangement we show how two different functionalities
Importance of a 3D forward modeling tool for surface wave analysis methods
NASA Astrophysics Data System (ADS)
Pageot, Damien; Le Feuvre, Mathieu; Donatienne, Leparoux; Philippe, Côte; Yann, Capdeville
2016-04-01
Since a few years, seismic surface waves analysis methods (SWM) have been widely developed and tested in the context of subsurface characterization and have demonstrated their effectiveness for sounding and monitoring purposes, e.g., high-resolution tomography of the principal geological units of California or real time monitoring of the Piton de la Fournaise volcano. Historically, these methods are mostly developed under the assumption of semi-infinite 1D layered medium without topography. The forward modeling is generally based on Thomson-Haskell matrix based modeling algorithm and the inversion is driven by Monte-Carlo sampling. Given their efficiency, SWM have been transfered to several scale of which civil engineering structures in order to, e.g., determine the so-called V s30 parameter or assess other critical constructional parameters in pavement engineering. However, at this scale, many structures may often exhibit 3D surface variations which drastically limit the efficiency of SWM application. Indeed, even in the case of an homogeneous structure, 3D geometry can bias the dispersion diagram of Rayleigh waves up to obtain discontinuous phase velocity curves which drastically impact the 1D mean velocity model obtained from dispersion inversion. Taking advantages of high-performance computing center accessibility and wave propagation modeling algorithm development, it is now possible to consider the use of a 3D elastic forward modeling algorithm instead of Thomson-Haskell method in the SWM inversion process. We use a parallelized 3D elastic modeling code based on the spectral element method which allows to obtain accurate synthetic data with very low numerical dispersion and a reasonable numerical cost. In this study, we choose dike embankments as an illustrative example. We first show that their longitudinal geometry may have a significant effect on dispersion diagrams of Rayleigh waves. Then, we demonstrate the necessity of 3D elastic modeling as a forward
Charged-particle Gun Design with 3D Finite-element Methods
NASA Astrophysics Data System (ADS)
Humphries, Stanley
2002-04-01
The DARHT second-axis injector poses a major challenge for computer simulation. The relativistic electrons are subject to strong beam-generated electric and magnetic forces. The beam and applied fields are fully three-dimensional. Furthermore, accurate field calculations at surfaces are critical to model Child-law emission. Although several 2D relativistic beam codes are available, there is presently no 3D tool that can address all important processes in the DARHT injector. As a result, we created the OmniTrak 3D finite-element code suite. This talk gives a basic tutorial on finite-element methods with emphasis on electron gun design via the ray-tracing technique. Four main areas are covered: 1) the mesh as a tool to organize space, 2) transformation of the Poisson equation through the minimum residual principle, 3) orbit tracking in a complex environment and 4) handling self-consistent beam-generated fields. The components of a volume mesh (elements, nodes and facets) are reviewed. We consider motivations for choosing a 3D mesh style: structured versus unstructured, tetrahedrons versus hexahedrons. We discuss methods for taking volume integrals over arbitrary hexahedrons through normal coordinates and shape functions, leading to the fundamental field equations. The special problems of 3D magnetic field solutions and the advantages of the reduced potential method are outlined. Accurate field interpolations for orbit calculations require fast identification of occupied elements. A method for fast element identification that also yields the orbit penetration point on the element surface is described. The final topics are the assignment of charge and current to meshes from calculated orbits and techniques for space-charge-limited emission from multiple arbitrary 3D surfaces.
The arctic winter stratosphere: simulated with a 3-D chemistry transport model
NASA Astrophysics Data System (ADS)
van den Broek, Martina Maria Petronella
2004-10-01
During the past two decades, the ozone layer has developed a “hole” each winter and spring above the Antarctic continent. Also in cold Arctic winters substantial stratospheric ozone depletion has been measured, although less than in the Antarctic stratosphere. In the Arctic winter stratosphere, the amount of ozone depletion varies interannually and within one winter, depending on polar vortex stability and temperature. The simulation of transport and chemical conversion of ozone and related species requires a three-dimensional (3D) chemistry transport model (CTM), because of the non-symmetric behaviour of the Arctic polar vortex. This thesis reports on several studies of the Arctic winter stratosphere carried out with such a CTM, using off-line meteorological fields. In Chapter II, chlorine activation and ozone depletion in the Arctic winter stratosphere of 1996-1997 are modelled with the newly developed stratospheric version of our CTM. Comparisons have been made with total O3 columns and ClO concentrations observed by satellites, and with ozone loss rates derived from observations during February and March 1997. ClO concentrations and ozone depletion are somewhat underestimated by the model. Key model parameters have been varied to explain this underestimation. Next to temperature, the formation mechanism of solid and/or liquid PSC particles constitutes the main model uncertainty. The representation of tracer transport is a third uncertain parameter, influencing both ozone and inorganic chlorine. In Chapter III, we have used the CTM with different horizontal resolutions to evaluate this stratospheric transport by simulating the long-lived tracers HF and CH4 during the Arctic winter of 1999/2000. Outside the vortex the model results agree well with the observations, but inside, the model underestimates the observed vertical gradient in HF and CH4. Too strong mixing through the vortex edge could be a cause for these model discrepancies, e.g. associated with the
Flexible 3D reconstruction method based on phase-matching in multi-sensor system.
Wu, Qingyang; Zhang, Baichun; Huang, Jinhui; Wu, Zejun; Zeng, Zeng
2016-04-01
Considering the measuring range limitation of a single sensor system, multi-sensor system has become essential in obtaining complete image information of the object in the field of 3D image reconstruction. However, for the traditional multi-sensors worked independently in its system, there was some point in calibrating each sensor system separately. And the calibration between all single sensor systems was complicated and required a long time. In this paper, we present a flexible 3D reconstruction method based on phase-matching in multi-sensor system. While calibrating each sensor, it realizes the data registration of multi-sensor system in a unified coordinate system simultaneously. After all sensors are calibrated, the whole 3D image data directly exist in the unified coordinate system, and there is no need to calibrate the positions between sensors any more. Experimental results prove that the method is simple in operation, accurate in measurement, and fast in 3D image reconstruction. PMID:27137020
A support-operator method for viscoelastic wave modelling in 3-D heterogeneous media
NASA Astrophysics Data System (ADS)
Ely, Geoffrey P.; Day, Steven M.; Minster, Jean-Bernard
2008-01-01
We apply the method of support operators (SOM) to solve the 3-D, viscoelastic equations of motion for use in earthquake simulations. SOM is a generalized finite-difference method that can utilize meshes of arbitrary structure and incorporate irregular geometry. Our implementation uses a 3-D, logically rectangular, hexahedral mesh. Calculations are second-order in space and time. A correction term is employed for suppression of spurious zero-energy modes (hourglass oscillations). We develop a free surface boundary condition, and an absorbing boundary condition using the method of perfectly matched layers (PML). Numerical tests using a layered material model in a highly deformed mesh show good agreement with the frequency-wavenumber method, for resolutions greater than 10 nodes per wavelength. We also test a vertically incident P wave on a semi-circular canyon, for which results match boundary integral solutions at resolutions greater that 20 nodes per wavelength. We also demonstrate excellent parallel scalability of our code.
A correction method of color projection fringes in 3D contour measurement
NASA Astrophysics Data System (ADS)
Song, Li-mei; Li, Zong-yan; Chen, Chang-man; Xi, Jiang-tao; Guo, Qing-hua; Li, Xiao-jie
2015-07-01
In the three-dimensional (3D) contour measurement, the phase shift profilometry (PSP) method is the most widely used one. However, the measurement speed of PSP is very low because of the multiple projections. In order to improve the measurement speed, color grating stripes are used for measurement in this paper. During the measurement, only one color sinusoidal fringe is projected on the measured object. Therefore, the measurement speed is greatly improved. Since there is coupling or interference phenomenon between the adjacent color grating stripes, a color correction method is used to improve the measurement results. A method for correcting nonlinear error of measurement system is proposed in this paper, and the sinusoidal property of acquired image after correction is better than that before correction. Experimental results show that with these correction methods, the measurement errors can be reduced. Therefore, it can support a good foundation for the high-precision 3D reconstruction.
Numerical solution of 3-D magnetotelluric using vector finite element method
NASA Astrophysics Data System (ADS)
Prihantoro, Rudy; Sutarno, Doddy; Nurhasan
2015-09-01
Magnetotelluric (MT) is a passive electromagnetic (EM) method which measure natural variations of electric and magnetic vector fields at the Earth surface to map subsurface electrical conductivity/resistivity structure. In this study, we obtained numerical solution of three-dimensional (3-D) MT using vector finite element method by solving second order Maxwell differential equation describing diffusion of plane wave through the conductive earth. Rather than the nodes of the element, the edges of the element is used as a vector basis to overcome the occurrence of nonphysical solutions that usually faced by scalar (node based) finite element method. Electric vector fields formulation was used and the resulting system of equation was solved using direct solution method to obtain the electric vector field distribution throughout the earth resistivity model structure. The resulting MT response functions was verified with 1-D layered Earth and 3-D2 COMMEMI outcropping structure. Good agreement is achieved for both structure models.
An improved 3D shape context registration method for non-rigid surface registration
NASA Astrophysics Data System (ADS)
Xiao, Di; Zahra, David; Bourgeat, Pierrick; Berghofer, Paula; Acosta Tamayo, Oscar; Wimberley, Catriona; Gregoire, Marie-Claude; Salvado, Olivier
2010-03-01
3D shape context is a method to define matching points between similar shapes as a pre-processing step to non-rigid registration. The main limitation of the approach is point mismatching, which includes long geodesic distance mismatch and neighbors crossing mismatch. In this paper, we propose a topological structure verification method to correct the long geodesic distance mismatch and a correspondence field smoothing method to correct the neighbors crossing mismatch. A robust 3D shape context model is proposed and further combined with thin-plate spline model for non-rigid surface registration. The method was tested on phantoms and rat hind limb skeletons from micro CT images. The results from experiments on mouse hind limb skeletons indicate that the approach is robust.
Robust method for extracting the pulmonary vascular trees from 3D MDCT images
NASA Astrophysics Data System (ADS)
Taeprasartsit, Pinyo; Higgins, William E.
2011-03-01
Segmentation of pulmonary blood vessels from three-dimensional (3D) multi-detector CT (MDCT) images is important for pulmonary applications. This work presents a method for extracting the vascular trees of the pulmonary arteries and veins, applicable to both contrast-enhanced and unenhanced 3D MDCT image data. The method finds 2D elliptical cross-sections and evaluates agreement of these cross-sections in consecutive slices to find likely cross-sections. It next employs morphological multiscale analysis to separate vessels from adjoining airway walls. The method then tracks the center of the likely cross-sections to connect them to the pulmonary vessels in the mediastinum and forms connected vascular trees spanning both lungs. A ground-truth study indicates that the method was able to detect on the order of 98% of the vessel branches having diameter >= 3.0 mm. The extracted vascular trees can be utilized for the guidance of safe bronchoscopic biopsy.
Small pitch fringe projection method with multiple linear fiber arrays for 3D shape measurement
NASA Astrophysics Data System (ADS)
Hayashi, Takumi; Fujigaki, Motoharu; Murata, Yorinobu
2014-07-01
3-D shape measurement systems by contactless method are required in the quality inspections of metal molds and electronic parts in industrial fields. A grating projection method with phase-shifting method has advantages of high precision and high speed. Recently, the size of a BGA (ball grid array) becomes smaller. So the pitch of a grating pattern projected onto the specimen should be smaller. In conventional method, fringe pattern is projected using an imaging lens. The focal depth becomes smaller in the case of reduced projection. It is therefore difficult to project a grating pattern with small pitch onto an object with large incident angles. Authors recently proposed a light source stepping method using a linear LED device. It is easy to shrink the projected grating pitch with a lens because this projection method does not use an imaging lens. The pitch of the projected grating depends on the width of the light source. There is a limit to shrink the projected grating pitch according to the size of the LED chip. In this paper, a small pitch fringe projection method with multiple linear fiber arrays for 3D shape measurement is proposed. The width of the fiber array is 30μm. It is one digit smaller than the width of the LED chip. The experimental result of 3-D shape measurement with small pitch projection with large incident angles is shown.
NASA Astrophysics Data System (ADS)
Hosseininaveh Ahmadabadian, Ali; Robson, Stuart; Boehm, Jan; Shortis, Mark
2013-04-01
Multi-View Stereo (MVS) as a low cost technique for precise 3D reconstruction can be a rival for laser scanners if the scale of the model is resolved. A fusion of stereo imaging equipment with photogrammetric bundle adjustment and MVS methods, known as photogrammetric MVS, can generate correctly scaled 3D models without using any known object distances. Although a huge number of stereo images (e.g. 200 high resolution images from a small object) captured of the object contains redundant data that allows detailed and accurate 3D reconstruction, the capture and processing time is increased when a vast amount of high resolution images are employed. Moreover, some parts of the object are often missing due to the lack of coverage of all areas. These problems demand a logical selection of the most suitable stereo camera views from the large image dataset. This paper presents a method for clustering and choosing optimal stereo or optionally single images from a large image dataset. The approach focusses on the two key steps of image clustering and iterative image selection. The method is developed within a software application called Imaging Network Designer (IND) and tested by the 3D recording of a gearbox and three metric reference objects. A comparison is made between IND and CMVS, which is a free package for selecting vantage images. The final 3D models obtained from the IND and CMVS approaches are compared with datasets generated with an MMDx Nikon Laser scanner. Results demonstrate that IND can provide a better image selection for MVS than CMVS in terms of surface coordinate uncertainty and completeness.
Development of direct-inverse 3-D methods for applied aerodynamic design and analysis
NASA Technical Reports Server (NTRS)
Carlson, Leland A.
1988-01-01
Several inverse methods have been compared and initial results indicate that differences in results are primarily due to coordinate systems and fuselage representations and not to design procedures. Further, results from a direct-inverse method that includes 3-D wing boundary layer effects, wake curvature, and wake displacement are presented. These results show that boundary layer displacements must be included in the design process for accurate results.
Performance and sensitivity evaluation of 3D spot detection methods in confocal microscopy.
Štěpka, Karel; Matula, Pavel; Matula, Petr; Wörz, Stefan; Rohr, Karl; Kozubek, Michal
2015-08-01
Reliable 3D detection of diffraction-limited spots in fluorescence microscopy images is an important task in subcellular observation. Generally, fluorescence microscopy images are heavily degraded by noise and non-specifically stained background, making reliable detection a challenging task. In this work, we have studied the performance and parameter sensitivity of eight recent methods for 3D spot detection. The study is based on both 3D synthetic image data and 3D real confocal microscopy images. The synthetic images were generated using a simulator modeling the complete imaging setup, including the optical path as well as the image acquisition process. We studied the detection performance and parameter sensitivity under different noise levels and under the influence of uneven background signal. To evaluate the parameter sensitivity, we propose a novel measure based on the gradient magnitude of the F1 score. We measured the success rate of the individual methods for different types of the image data and found that the type of image degradation is an important factor. Using the F1 score and the newly proposed sensitivity measure, we found that the parameter sensitivity is not necessarily proportional to the success rate of a method. This also provided an explanation why the best performing method for synthetic data was outperformed by other methods when applied to the real microscopy images. On the basis of the results obtained, we conclude with the recommendation of the HDome method for data with relatively low variations in quality, or the Sorokin method for image sets in which the quality varies more. We also provide alternative recommendations for high-quality images, and for situations in which detailed parameter tuning might be deemed expensive. PMID:26033916
Crossover from 3D to 2D Quantum Transport in Bi2Se3/In2Se3 Superlattices
NASA Astrophysics Data System (ADS)
Yanfei, Zhao; Haiwen, Liu; Xin, Guo; Ying, Jiang; Yi, Sun; Huichao, Wang; Yong, Wang; Handong, Li; Maohai, Xie; Xincheng, Xie; Jian, Wang
2015-03-01
The topological insulator/normal insulator (TI/NI) superlattices (SLs) with multiple Dirac channels are predicted to offer great opportunity to design novel materials and investigate new quantum phenomena. Here, we report first transport studies on the SLs composed of TI Bi2Se3 layers sandwiched by NI In2Se3 layers artificially grown by molecular beam epitaxy (MBE). The transport properties of two kinds of SL samples show convincing evidence that the transport dimensionality changes from three-dimensional (3D) to two-dimensional (2D) when decreasing the thickness of building block Bi2Se3 layers, corresponding to the crossover from coherent TI transport to separated TI channels. Our findings provide the possibility to realizing 3D surface states in TI/NI SLs.
Crossover from 3D to 2D quantum transport in Bi2Se3/In2Se3 superlattices.
Zhao, Yanfei; Liu, Haiwen; Guo, Xin; Jiang, Ying; Sun, Yi; Wang, Huichao; Wang, Yong; Li, Han-Dong; Xie, Mao-Hai; Xie, Xin-Cheng; Wang, Jian
2014-09-10
The topological insulator/normal insulator (TI/NI) superlattices (SLs) with multiple Dirac channels are predicted to offer great opportunity to design novel materials and investigate new quantum phenomena. Here, we report first transport studies on the SLs composed of TI Bi2Se3 layers sandwiched by NI In2Se3 layers artificially grown by molecular beam epitaxy (MBE). The transport properties of two kinds of SL samples show convincing evidence that the transport dimensionality changes from three-dimensional (3D) to two-dimensional (2D) when decreasing the thickness of building block Bi2Se3 layers, corresponding to the crossover from coherent TI transport to separated TI channels. Our findings provide the possibility to realizing "3D surface states" in TI/NI SLs. PMID:25102289
A highly heterogeneous 3D PWR core benchmark: deterministic and Monte Carlo method comparison
NASA Astrophysics Data System (ADS)
Jaboulay, J.-C.; Damian, F.; Douce, S.; Lopez, F.; Guenaut, C.; Aggery, A.; Poinot-Salanon, C.
2014-06-01
Physical analyses of the LWR potential performances with regards to the fuel utilization require an important part of the work dedicated to the validation of the deterministic models used for theses analyses. Advances in both codes and computer technology give the opportunity to perform the validation of these models on complex 3D core configurations closed to the physical situations encountered (both steady-state and transient configurations). In this paper, we used the Monte Carlo Transport code TRIPOLI-4®; to describe a whole 3D large-scale and highly-heterogeneous LWR core. The aim of this study is to validate the deterministic CRONOS2 code to Monte Carlo code TRIPOLI-4®; in a relevant PWR core configuration. As a consequence, a 3D pin by pin model with a consistent number of volumes (4.3 millions) and media (around 23,000) is established to precisely characterize the core at equilibrium cycle, namely using a refined burn-up and moderator density maps. The configuration selected for this analysis is a very heterogeneous PWR high conversion core with fissile (MOX fuel) and fertile zones (depleted uranium). Furthermore, a tight pitch lattice is selcted (to increase conversion of 238U in 239Pu) that leads to harder neutron spectrum compared to standard PWR assembly. In these conditions two main subjects will be discussed: the Monte Carlo variance calculation and the assessment of the diffusion operator with two energy groups for the core calculation.
NASA Astrophysics Data System (ADS)
Fargier, Yannick; Dore, Ludovic; Antoine, Raphael; Palma Lopes, Sérgio; Fauchard, Cyrille
2016-04-01
The extraction of subsurface materials is a key element for the economy of a nation. However, natural degradation of underground quarries is a major issue from an economic and public safety point of view. Consequently, the quarries stakeholders require relevant tools to define hazards associated to these structures. Safety assessment methods of underground quarries are recent and mainly based on rock physical properties. This kind of method leads to a certain homogeneity assumption of pillar internal properties that can cause an underestimation of the risk. Electrical Resistivity Imaging (ERI) is a widely used method that possesses two advantages to overcome this limitation. The first is to provide a qualitative understanding for the detection and monitoring of anomalies in the pillar body (e.g. faults). The second is to provide a quantitative description of the electrical resistivity distribution inside the pillar. This quantitative description can be interpreted with constitutive laws to help decision support (water content decreases the mechanical resistance of a chalk). However, conventional 2D and 3D Imaging techniques are usually applied to flat surface surveys or to surfaces with moderate topography. A 3D inversion of more complex media (case of the pillar) requires a full consideration of the geometry that was never taken into account before. The Photogrammetric technique presents a cost effective solution to obtain an accurate description of the external geometry of a complex media. However, this method has never been fully coupled with a geophysical method to enhance/improve the inversion process. Consequently we developed a complete procedure showing that photogrammetric and ERI tools can be efficiently combined to assess a complex 3D structure. This procedure includes in a first part a photogrammetric survey, a processing stage with an open source software and a post-processing stage finalizing a 3D surface model. The second part necessitates the
Finite volume and finite element methods applied to 3D laminar and turbulent channel flows
Louda, Petr; Příhoda, Jaromír; Sváček, Petr; Kozel, Karel
2014-12-10
The work deals with numerical simulations of incompressible flow in channels with rectangular cross section. The rectangular cross section itself leads to development of various secondary flow patterns, where accuracy of simulation is influenced by numerical viscosity of the scheme and by turbulence modeling. In this work some developments of stabilized finite element method are presented. Its results are compared with those of an implicit finite volume method also described, in laminar and turbulent flows. It is shown that numerical viscosity can cause errors of same magnitude as different turbulence models. The finite volume method is also applied to 3D turbulent flow around backward facing step and good agreement with 3D experimental results is obtained.
An Application of the Method of Arbitrary Lines to 3D Elastic Stress Analysis
NASA Astrophysics Data System (ADS)
Kaminishi, Ken; Ando, Ryuma
The MAL (Method of Arbitrary Lines) is a technique of reducing a partial differential equation to a system of ordinary differential equations. It is known that relevant use of this procedure yields high accuracy in some problems of two-dimensional elasticity and elastoplasticity. Since the basic concept of MAL is simple and based on generality, it is expected that many problems in other fields will be effectively solvable by this method. In this study, we consider the application of MAL to 3D (three-dimensional) elasticity analysis. We first give a MAL formulation of 3D elasticity problems, and demonstrate its effectiveness and accuracy for a typical problem. The reported numerical results are compared with the exact solution or that of the finite element method (FEM).
A fast method to measure the 3D surface of the human heart
NASA Astrophysics Data System (ADS)
Cao, Yiping; Su, Xianyu; Xiang, Liqun; Chen, Wenjing; Zhang, Qican
2003-12-01
Three-dimensional (3-D) automatic measurement of an object is widely used in many fields. In Biology and Medicine society, it can be applicable for surgery, orthopedics, viscera disease analysis and diagnosis etc. Here a new fast method to measure the 3D surface of human heart is proposed which can provide doctors a lot of information, such as the size of heart profile, the sizes of the left or right heart ventricle, and the curvature center and radius of heart ventricle, to fully analyze and diagnose pathobiology of human heart. The new fast method is optically and noncontacted and based upon the Phase Measurement Profilometry (PMP), which has higher measuring precision. A human heart specimen experiment has verified our method.
Analysis of method of 3D shape reconstruction using scanning deflectometry
NASA Astrophysics Data System (ADS)
Novák, Jiří; Novák, Pavel; Mikš, Antonín.
2013-04-01
This work presents a scanning deflectometric approach to solving a 3D surface reconstruction problem, which is based on measurements of a surface gradient of optically smooth surfaces. It is shown that a description of this problem leads to a nonlinear partial differential equation (PDE) of the first order, from which the surface shape can be reconstructed numerically. The method for effective finding of the solution of this differential equation is proposed, which is based on the transform of the problem of PDE solving to the optimization problem. We describe different types of surface description for the shape reconstruction and a numerical simulation of the presented method is performed. The reconstruction process is analyzed by computer simulations and presented on examples. The performed analysis confirms a robustness of the reconstruction method and a good possibility for measurements and reconstruction of the 3D shape of specular surfaces.
Intrathoracic tumour motion estimation from CT imaging using the 3D optical flow method
NASA Astrophysics Data System (ADS)
Guerrero, Thomas; Zhang, Geoffrey; Huang, Tzung-Chi; Lin, Kang-Ping
2004-09-01
The purpose of this work was to develop and validate an automated method for intrathoracic tumour motion estimation from breath-hold computed tomography (BH CT) imaging using the three-dimensional optical flow method (3D OFM). A modified 3D OFM algorithm provided 3D displacement vectors for each voxel which were used to map tumour voxels on expiration BH CT onto inspiration BH CT images. A thoracic phantom and simulated expiration/inspiration BH CT pairs were used for validation. The 3D OFM was applied to the measured inspiration and expiration BH CT images from one lung cancer and one oesophageal cancer patient. The resulting displacements were plotted in histogram format and analysed to provide insight regarding the tumour motion. The phantom tumour displacement was measured as 1.20 and 2.40 cm with full-width at tenth maximum (FWTM) for the distribution of displacement estimates of 0.008 and 0.006 cm, respectively. The maximum error of any single voxel's motion estimate was 1.1 mm along the z-dimension or approximately one-third of the z-dimension voxel size. The simulated BH CT pairs revealed an rms error of less than 0.25 mm. The displacement of the oesophageal tumours was nonuniform and up to 1.4 cm, this was a new finding. A lung tumour maximum displacement of 2.4 cm was found in the case evaluated. In conclusion, 3D OFM provided an accurate estimation of intrathoracic tumour motion, with estimated errors less than the voxel dimension in a simulated motion phantom study. Surprisingly, oesophageal tumour motion was large and nonuniform, with greatest motion occurring at the gastro-oesophageal junction. Presented at The IASTED Second International Conference on Biomedical Engineering (BioMED 2004), Innsbruck, Austria, 16-18 February 2004.
Quantum transport theory of 3D time-reversal invariant topological insulators
NASA Astrophysics Data System (ADS)
Dellabetta, Brian Jon
We consider the potential technological role of a recently predicted and discovered phase of quantum matter - topological insulators (TIs), which are characterized by an insulating bulk and topologically protected, gapless, spin-momentum locked surface modes. Precise engineering of these gapless modes may yield new potential materials for novel electronic devices, but many materials issues and open questions in application remain in the nascent field. The quasiparticle dynamics of TI systems can be elegantly written in terms of a low-energy effective momentum-space Hamiltonian, but analytic methods quickly become intractable in multifarious systems and disordered heterostructures which in general lack translational invariance, as momentum is no longer a good quantum number. Computational methods possess a clear advantage in this regime, for understanding systems in which geometry, contact layout, and disorder play a dominant role. We employ computationally intensive methods to calculate observable, non-equilibrium transport dynamics of real-space topological systems, to propose and identify experimental signatures of topological behavior, and to connect interesting experimental observations to the underlying topological properties in normal, disordered, and superconducting systems. The customizability of these computational methods allows us to determine the salient underlying physics involved in a number of different scenarios, including surface transport corrugated TI channels, the Aharonov-Bohm effect in TI nanowires, supercurrent in TI Josephson junctions, and the superconducting proximity effect and resulting transport in TI-superconductor heterostructures. In doing so, we expand the understanding of quantum and mesoscopic transport in heterostructured TI systems as a first step in exploring their long-term place in novel device applications.
NASA Astrophysics Data System (ADS)
Uehara, Shin-ichi; Ujike, Hiroyasu; Hamagishi, Goro; Taira, Kazuki; Koike, Takafumi; Kato, Chiaki; Nomura, Toshio; Horikoshi, Tsutomu; Mashitani, Ken; Yuuki, Akimasa; Izumi, Kuniaki; Hisatake, Yuzo; Watanabe, Naoko; Umezu, Naoaki; Nakano, Yoshihiko
2010-02-01
We are engaged in international standardization activities for 3D displays. We consider that for a sound development of 3D displays' market, the standards should be based on not only mechanism of 3D displays, but also human factors for stereopsis. However, we think that there is no common understanding on what the 3D display should be and that the situation makes developing the standards difficult. In this paper, to understand the mechanism and human factors, we focus on a double image, which occurs in some conditions on an autostereoscopic display. Although the double image is generally considered as an unwanted effect, we consider that whether the double image is unwanted or not depends on the situation and that there are some allowable double images. We tried to classify the double images into the unwanted and the allowable in terms of the display mechanism and visual ergonomics for stereopsis. The issues associated with the double image are closely related to performance characteristics for the autostereoscopic display. We also propose performance characteristics, measurement and analysis methods to represent interocular crosstalk and motion parallax.
Comparison of Parallel MRI Reconstruction Methods for Accelerated 3D Fast Spin-Echo Imaging
Xiao, Zhikui; Hoge, W. Scott; Mulkern, R.V.; Zhao, Lei; Hu, Guangshu; Kyriakos, Walid E.
2014-01-01
Parallel MRI (pMRI) achieves imaging acceleration by partially substituting gradient-encoding steps with spatial information contained in the component coils of the acquisition array. Variable-density subsampling in pMRI was previously shown to yield improved two-dimensional (2D) imaging in comparison to uniform subsampling, but has yet to be used routinely in clinical practice. In an effort to reduce acquisition time for 3D fast spin-echo (3D-FSE) sequences, this work explores a specific nonuniform sampling scheme for 3D imaging, subsampling along two phase-encoding (PE) directions on a rectilinear grid. We use two reconstruction methods—2D-GRAPPA-Operator and 2D-SPACE RIP—and present a comparison between them. We show that high-quality images can be reconstructed using both techniques. To evaluate the proposed sampling method and reconstruction schemes, results via simulation, phantom study, and in vivo 3D human data are shown. We find that fewer artifacts can be seen in the 2D-SPACE RIP reconstructions than in 2D-GRAPPA-Operator reconstructions, with comparable reconstruction times. PMID:18727083
Novel high speed method using gray level vector modulation for 3D shape measurement
NASA Astrophysics Data System (ADS)
Lin, Gui-Wen; Li, Dong; Tian, Jin-Dong
2014-11-01
Binocular Vision Technique is widely used in three-dimensional (3-D) measurement. Matching of pictures captured from two cameras is the most critical and difficult step in 3-D shape reconstruction. The method combines codedstructured light and spatial phase is usually adopted. However, being time consuming in matching, this method could not meet the requirements of real-time 3-D vision. In order to satisfy the high speed characteristic of real-time measurement, a novel method using gray level vector modulation is introduced. Combining binary code with gray coding principle, new coding patterns using gray level vector method is designed and projected onto the object surface. Each pixel corresponds to the designed sequence of gray values as a feature vector. The unique gray level vector is then dimensionally reduced to a resulting value which could be used as characteristic information for binocular matching. Experimental results further demonstrated the correctness and feasibility of the proposed method with fewer component patterns and less computational time.
NASA Technical Reports Server (NTRS)
Strahan, Susan E.; Douglass, Anne R.; Einaudi, Franco (Technical Monitor)
2001-01-01
The Global Modeling Initiative (GMI) Team developed objective criteria for model evaluation in order to identify the best representation of the stratosphere. This work created a method to quantitatively and objectively discriminate between different models. In the original GMI study, 3 different meteorological data sets were used to run an offline chemistry and transport model (CTM). Observationally-based grading criteria were derived and applied to these simulations and various aspects of stratospheric transport were evaluated; grades were assigned. Here we report on the application of the GMI evaluation criteria to CTM simulations integrated with a new assimilated wind data set and a new general circulation model (GCM) wind data set. The Finite Volume Community Climate Model (FV-CCM) is a new GCM developed at Goddard which uses the NCAR CCM physics and the Lin and Rood advection scheme. The FV-Data Assimilation System (FV-DAS) is a new data assimilation system which uses the FV-CCM as its core model. One year CTM simulations of 2.5 degrees longitude by 2 degrees latitude resolution were run for each wind data set. We present the evaluation of temperature and annual transport cycles in the lower and middle stratosphere in the two new CTM simulations. We include an evaluation of high latitude transport which was not part of the original GMI criteria. Grades for the new simulations will be compared with those assigned during the original GMT evaluations and areas of improvement will be identified.
3D modelling of the transport and fate of riverine fine sediment exported to a semi-enclosed system
NASA Astrophysics Data System (ADS)
Delandmeter, Philippe; Lambrechts, Jonathan; Lewis, Stephen; Legat, Vincent; Deleersnijder, Eric; Wolanski, Eric
2015-04-01
Understanding the transport and fate of suspended sediment exported by rivers is crucial for the management of sensitive marine ecosystems. Sediment transport and fate can vary considerably depending on the geophysical characteristics of the offshore environment (i.e. open, semi-enclosed and enclosed systems and the nature of the continental shelf). In this presentation, we focus on a semi-enclosed setting in the Great Barrier Reef, NE Australia. In this system, the large tropical Burdekin River discharges to a long and narrow continental shelf containing numerous headlands and embayments. Using a new 3D sediment model we developed and SLIM 3D, a Finite Element 3D model for coastal flows, we highlight the key processes of sediment transport for such a system. We validate the model with available measured data from the region. Wind direction and speed during the high river flows are showed to largely control the dynamics and final fate of the sediments. Most (71%) of the sediment load delivered by the river is deposited and retained near the river mouth. The remaining sediment is transported further afield in riverine freshwater plumes. The suspended sediment transported longer distances in the freshwater plumes can reach sensitive marine ecosystems. These results are compared to previous studies on the Burdekin River sediment fate and differences are analysed. The model suggests that wind-driven resuspension events will redistribute sediments within an embayment but have little influence on transporting sediments from bay to bay.
Earthquake source tensor inversion with the gCAP method and 3D Green's functions
NASA Astrophysics Data System (ADS)
Zheng, J.; Ben-Zion, Y.; Zhu, L.; Ross, Z.
2013-12-01
We develop and apply a method to invert earthquake seismograms for source properties using a general tensor representation and 3D Green's functions. The method employs (i) a general representation of earthquake potency/moment tensors with double couple (DC), compensated linear vector dipole (CLVD), and isotropic (ISO) components, and (ii) a corresponding generalized CAP (gCap) scheme where the continuous wave trains are broken into Pnl and surface waves (Zhu & Ben-Zion, 2013). For comparison, we also use the waveform inversion method of Zheng & Chen (2012) and Ammon et al. (1998). Sets of 3D Green's functions are calculated on a grid of 1 km3 using the 3-D community velocity model CVM-4 (Kohler et al. 2003). A bootstrap technique is adopted to establish robustness of the inversion results using the gCap method (Ross & Ben-Zion, 2013). Synthetic tests with 1-D and 3-D waveform calculations show that the source tensor inversion procedure is reasonably reliable and robust. As initial application, the method is used to investigate source properties of the March 11, 2013, Mw=4.7 earthquake on the San Jacinto fault using recordings of ~45 stations up to ~0.2Hz. Both the best fitting and most probable solutions include ISO component of ~1% and CLVD component of ~0%. The obtained ISO component, while small, is found to be a non-negligible positive value that can have significant implications for the physics of the failure process. Work on using higher frequency data for this and other earthquakes is in progress.
A Quality Assurance Method that Utilizes 3D Dosimetry and Facilitates Clinical Interpretation
Oldham, Mark; Thomas, Andrew; O'Daniel, Jennifer; Juang, Titania; Ibbott, Geoffrey; Adamovics, John; Kirkpatrick, John P.
2012-10-01
Purpose: To demonstrate a new three-dimensional (3D) quality assurance (QA) method that provides comprehensive dosimetry verification and facilitates evaluation of the clinical significance of QA data acquired in a phantom. Also to apply the method to investigate the dosimetric efficacy of base-of-skull (BOS) intensity-modulated radiotherapy (IMRT) treatment. Methods and Materials: Two types of IMRT QA verification plans were created for 6 patients who received BOS IMRT. The first plan enabled conventional 2D planar IMRT QA using the Varian portal dosimetry system. The second plan enabled 3D verification using an anthropomorphic head phantom. In the latter, the 3D dose distribution was measured using the DLOS/Presage dosimetry system (DLOS = Duke Large-field-of-view Optical-CT System, Presage Heuris Pharma, Skillman, NJ), which yielded isotropic 2-mm data throughout the treated volume. In a novel step, measured 3D dose distributions were transformed back to the patient's CT to enable calculation of dose-volume histograms (DVH) and dose overlays. Measured and planned patient DVHs were compared to investigate clinical significance. Results: Close agreement between measured and calculated dose distributions was observed for all 6 cases. For gamma criteria of 3%, 2 mm, the mean passing rate for portal dosimetry was 96.8% (range, 92.0%-98.9%), compared to 94.9% (range, 90.1%-98.9%) for 3D. There was no clear correlation between 2D and 3D passing rates. Planned and measured dose distributions were evaluated on the patient's anatomy, using DVH and dose overlays. Minor deviations were detected, and the clinical significance of these are presented and discussed. Conclusions: Two advantages accrue to the methods presented here. First, treatment accuracy is evaluated throughout the whole treated volume, yielding comprehensive verification. Second, the clinical significance of any deviations can be assessed through the generation of DVH curves and dose overlays on the patient
Enhanced Rgb-D Mapping Method for Detailed 3d Modeling of Large Indoor Environments
NASA Astrophysics Data System (ADS)
Tang, Shengjun; Zhu, Qing; Chen, Wu; Darwish, Walid; Wu, Bo; Hu, Han; Chen, Min
2016-06-01
RGB-D sensors are novel sensing systems that capture RGB images along with pixel-wise depth information. Although they are widely used in various applications, RGB-D sensors have significant drawbacks with respect to 3D dense mapping of indoor environments. First, they only allow a measurement range with a limited distance (e.g., within 3 m) and a limited field of view. Second, the error of the depth measurement increases with increasing distance to the sensor. In this paper, we propose an enhanced RGB-D mapping method for detailed 3D modeling of large indoor environments by combining RGB image-based modeling and depth-based modeling. The scale ambiguity problem during the pose estimation with RGB image sequences can be resolved by integrating the information from the depth and visual information provided by the proposed system. A robust rigid-transformation recovery method is developed to register the RGB image-based and depth-based 3D models together. The proposed method is examined with two datasets collected in indoor environments for which the experimental results demonstrate the feasibility and robustness of the proposed method
Parlea, Lorena G; Sweeney, Blake A; Hosseini-Asanjan, Maryam; Zirbel, Craig L; Leontis, Neocles B
2016-07-01
RNA 3D motifs occupy places in structured RNA molecules that correspond to the hairpin, internal and multi-helix junction "loops" of their secondary structure representations. As many as 40% of the nucleotides of an RNA molecule can belong to these structural elements, which are distinct from the regular double helical regions formed by contiguous AU, GC, and GU Watson-Crick basepairs. With the large number of atomic- or near atomic-resolution 3D structures appearing in a steady stream in the PDB/NDB structure databases, the automated identification, extraction, comparison, clustering and visualization of these structural elements presents an opportunity to enhance RNA science. Three broad applications are: (1) identification of modular, autonomous structural units for RNA nanotechnology, nanobiology and synthetic biology applications; (2) bioinformatic analysis to improve RNA 3D structure prediction from sequence; and (3) creation of searchable databases for exploring the binding specificities, structural flexibility, and dynamics of these RNA elements. In this contribution, we review methods developed for computational extraction of hairpin and internal loop motifs from a non-redundant set of high-quality RNA 3D structures. We provide a statistical summary of the extracted hairpin and internal loop motifs in the most recent version of the RNA 3D Motif Atlas. We also explore the reliability and accuracy of the extraction process by examining its performance in clustering recurrent motifs from homologous ribosomal RNA (rRNA) structures. We conclude with a summary of remaining challenges, especially with regard to extraction of multi-helix junction motifs. PMID:27125735
Equivalent Body Force Finite Elements Method and 3-D Earth Model Applied In 2004 Sumatra Earthquake
NASA Astrophysics Data System (ADS)
Qu, W.; Cheng, H.; Shi, Y.
2015-12-01
The 26 December 2004 Sumatra-Andaman earthquake with moment magnitude (Mw) of 9.1 to 9.3 is the first great earthquake recorded by digital broadband, high-dynamic-range seismometers and global positioning system (GPS) equipment, which recorded many high-quality geophysical data sets. The spherical curvature is not negligible in far field especially for large event and the real Earth is laterally inhomogeneity and the analytical results still are difficult to explain the geodetic measurements. We use equivalent body force finite elements method Zhang et al. (2015) and mesh the whole earth, to compute global co-seismic displacements using four fault slip models of the 2004 Sumatra earthquake provided by different authors. Comparisons of calculated co-seismic displacements and GPS show that the confidences are well in near field for four models, and the confidences are according to different models. In the whole four models, the Chlieh model (Chlieh et al., 2007) is the best as this slip model not only accord well with near field data but also far field data. And then we use the best slip model, Chlieh model to explore influence of three dimensional lateral earth structure on both layered spherically symmetric (PREM) and real 3-D heterogeneous earth model (Crust 1.0 model and GyPSuM). Results show that the effects of 3-D heterogeneous earth model are not negligible and decrease concomitantly with increasing distance from the epicenter. The relative effects of 3-D crust model are 23% and 40% for horizontal and vertical displacements, respectively. The effects of the 3-D mantle model are much smaller than that of 3-D crust model but with wider impacting area.
A novel 3D constellation-masked method for physical security in hierarchical OFDMA system.
Zhang, Lijia; Liu, Bo; Xin, Xiangjun; Liu, Deming
2013-07-01
This paper proposes a novel 3D constellation-masked method to ensure the physical security in hierarchical optical orthogonal frequency division multiplexing access (OFDMA) system. The 3D constellation masking is executed on the two levels of hierarchical modulation and among different OFDM subcarriers, which is realized by the masking vectors. The Lorenz chaotic model is adopted for the generation of masking vectors in the proposed scheme. A 9.85 Gb/s encrypted hierarchical QAM OFDM signal is successfully demonstrated in the experiment. The performance of illegal optical network unit (ONU) with different masking vectors is also investigated. The proposed method is demonstrated to be secure and efficient against the commonly known attacks in the experiment. PMID:23842348
Methods for Measuring the Orientation and Rotation Rate of 3D-printed Particles in Turbulence.
Cole, Brendan C; Marcus, Guy G; Parsa, Shima; Kramel, Stefan; Ni, Rui; Voth, Greg A
2016-01-01
Experimental methods are presented for measuring the rotational and translational motion of anisotropic particles in turbulent fluid flows. 3D printing technology is used to fabricate particles with slender arms connected at a common center. Shapes explored are crosses (two perpendicular rods), jacks (three perpendicular rods), triads (three rods in triangular planar symmetry), and tetrads (four arms in tetrahedral symmetry). Methods for producing on the order of 10,000 fluorescently dyed particles are described. Time-resolved measurements of their orientation and solid-body rotation rate are obtained from four synchronized videos of their motion in a turbulent flow between oscillating grids with Rλ = 91. In this relatively low-Reynolds number flow, the advected particles are small enough that they approximate ellipsoidal tracer particles. We present results of time-resolved 3D trajectories of position and orientation of the particles as well as measurements of their rotation rates. PMID:27404898
Optoranger: A 3D pattern matching method for bin picking applications
NASA Astrophysics Data System (ADS)
Sansoni, Giovanna; Bellandi, Paolo; Leoni, Fabio; Docchio, Franco
2014-03-01
This paper presents a new method, based on 3D vision, for the recognition of free-form objects in the presence of clutters and occlusions, ideal for robotic bin picking tasks. The method can be considered as a compromise between complexity and effectiveness. A 3D point cloud representing the scene is generated by a triangulation-based scanning system, where a fast camera acquires a blade projected by a laser source. Image segmentation is based on 2D images, and on the estimation of the distances between point pairs, to search for empty areas. Object recognition is performed using commercial software libraries integrated with custom-developed segmentation algorithms, and a database of model clouds created by means of the same scanning system.
Efficient solution on solving 3D Maxwell equations using stable semi-implicit splitting method
NASA Astrophysics Data System (ADS)
Cen, Wei; Gu, Ning
2016-05-01
In this paper, we propose an efficient solution on solving 3-dimensional (3D) time-domain Maxwell equations using the semi-implicit Crank-Nicholson (CN) method for time domain discretization with advantage of unconditional time stability. By applying the idea of fractional steps method (FSM) to the CN scheme, the proposed method provides a much simpler and efficient implementation than a direct implementation of the CN scheme. Compared with the alternating-direction implicit (ADI) method and explicit finite-difference time-domain approach (FDTD), it significantly saves the computational resource like memory and CPU time while remains similar numerical accuracy.
Perfetti, Christopher M; Rearden, Bradley T
2014-01-01
This work introduces a new approach for calculating sensitivity coefficients for generalized neutronic responses to nuclear data uncertainties using continuous-energy Monte Carlo methods. The approach presented in this paper, known as the GEAR-MC method, allows for the calculation of generalized sensitivity coefficients for multiple responses in a single Monte Carlo calculation with no nuclear data perturbations or knowledge of nuclear covariance data. The theory behind the GEAR-MC method is presented here, and proof of principle is demonstrated by using the GEAR-MC method to calculate sensitivity coefficients for responses in several 3D, continuous-energy Monte Carlo applications.
Gras, Christophe; Smith, Nikaïa; Sengmanivong, Lucie; Gandini, Mariana; Kubelka, Claire Fernandes; Herbeuval, Jean-Philippe
2013-01-31
The apoptotic ligand TNF-related apoptosis ligand (TRAIL) is expressed on the membrane of immune cells during HIV infection. The intracellular stockade of TRAIL in human primary CD4(+) T cells is not known. Here we investigated whether primary CD4(+) T cells expressed TRAIL in their intracellular compartment and whether TRAIL is relocalized on the plasma membrane under HIV activation. We found that TRAIL protein was stocked in intracellular compartment in non activated CD4(+) T cells and that the total level of TRAIL protein was not increased under HIV-1 stimulation. However, TRAIL was massively relocalized on plasma membrane when cells were cultured with HIV. Using three dimensional (3D) microscopy we localized TRAIL protein in human T cells and developed a new method to visualize plasma membrane without the need of a membrane marker. This method used the 3D interactive surface plot and bright light acquired images. PMID:23085529