1,2,3-D Diffusion Depletion Multi-Group
Energy Science and Technology Software Center (ESTSC)
1992-04-20
CITATION is designed to solve problems using the finite difference representation of neutron diffusion theory, treating up to three space dimensions with arbitrary group to group scattering. X-y-z, theta-r-z, hexagonal z, and triagonal z geometries may be treated. Depletion problems may be solved and fuel managed for multi-cycle analysis. Extensive first order perturbation results may be obtained given microscopic data and nuclide concentrations. Statics problems may be solved and perturbation results obtained with microscopic data.
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.
Solves the Multigroup Neutron Diffusion Equation
Energy Science and Technology Software Center (ESTSC)
1995-06-23
GNOMER is a program which solves the multigroup neutron diffusion equation in 1D, 2D and 3D cartesian geometry. The program is designed to calculate the global core power distributions (with thermohydraulic feedbacks), as well as power distribution and homogenized cross sections over a fuel assembly.
Multigroup Complex Geometry Neutron Diffusion Code System.
Energy Science and Technology Software Center (ESTSC)
2002-12-18
Version 01 SNAP-3D is based on SNAP2 and is a one- two- or three-dimensional multigroup diffusion code system. It is primarily intended for neutron diffusion calculations, but it can also carry out gamma-ray calculations if the diffusion approximation is accurate enough. It is suitable for fast and thermal reactor core calculations and for shield calculations. SNAP-3D can solve the multi-group neutron diffusion equations using finite difference methods in (x,y,z), (r,theta,z), (TRI,z), (HEX,z) or (spherical) coordinates.more » The one-dimensional slab and cylindrical geometries and the two-dimensional (x,y), (r,z), (r,theta), (HEX) and (TRI) are all treated as simple special cases of three-dimensional geometries. Numerous reflective and periodic symmetry options are available and may be used to reduce the number of mesh points necessary to represent the system. Extrapolation lengths can be specified at internal and external boundaries. The problem classes are: 1) eigenvalue search for critical k-effective, 2) eigenvalue search for critical buckling, 3) eigenvalue search for critical time-constant, 4) fixed source problems in which the sources are functions of regions, 5) fixed source problems in which the sources are provided, on disc, for every mesh point and group.« less
Multigroup 3-Dimensional Neutron Diffusion Nodal Code System with Thermohydraulic Feedbacks.
Energy Science and Technology Software Center (ESTSC)
1994-02-07
Version 01 GNOMER is a program which solves the multigroup neutron diffusion equation on coarse mesh in 1D, 2D, and 3D Cartesian geometry. The program is designed to calculate the global core power distributions (with thermohydraulic feedbacks) as well as power distributions and homogenized cross sections over a fuel assembly.
Radiosity diffusion model in 3D
NASA Astrophysics Data System (ADS)
Riley, Jason D.; Arridge, Simon R.; Chrysanthou, Yiorgos; Dehghani, Hamid; Hillman, Elizabeth M. C.; Schweiger, Martin
2001-11-01
We present the Radiosity-Diffusion model in three dimensions(3D), as an extension to previous work in 2D. It is a method for handling non-scattering spaces in optically participating media. We present the extension of the model to 3D including an extension to the model to cope with increased complexity of the 3D domain. We show that in 3D more careful consideration must be given to the issues of meshing and visibility to model the transport of light within reasonable computational bounds. We demonstrate the model to be comparable to Monte-Carlo simulations for selected geometries, and show preliminary results of comparisons to measured time-resolved data acquired on resin phantoms.
950809 Charged particle transport updated multi-group diffusion
Corman, E.G.; Perkins, S.T.; Dairiki, N.T.
1995-09-01
In 1974, a charged particle transport scheme was introduced which utilized a multi-group diffusion method for the spatial transport and slowing down of energetic ions in a hot plasma. In this treatment a diffusion coefficient was used which was flux-limited to provide, hopefully, some degree of accuracy when the slowing down of an energetic charged particle is dominated by Coulomb collisions with thermal ions and electrons in a plasma medium. An advantage of this method was a very fast, memory-contained program for calculating the behavior of energetic charged particles which resulted in smoothly varying particle number densities and energy depositions. The main limitation of the original multi-group charged particle diffusion scheme is its constraint to a basic ten group structure; the same ten group structure for each of the five energetic ions tracked. This is regarded as a severe limitation, inasmuch as more groups would be desired to simulate more accurately the corresponding Monte Carlo results of energies deposited over spatial zones from a charged particle source. More generally, it seems preferable to have a different group structure for each particle type since they are created at inherently different energies. In this paper, the basic theory and multi-group description will be given. This is followed by the specific techniques that were used to solve the resultant equations. Finally, the modifications that were made to the cross section data as well as the methods used for energy and momentum deposition are described.
A multigroup radiation diffusion test problem: Comparison of code results with analytic solution
Shestakov, A I; Harte, J A; Bolstad, J H; Offner, S R
2006-12-21
We consider a 1D, slab-symmetric test problem for the multigroup radiation diffusion and matter energy balance equations. The test simulates diffusion of energy from a hot central region. Opacities vary with the cube of the frequency and radiation emission is given by a Wien spectrum. We compare results from two LLNL codes, Raptor and Lasnex, with tabular data that define the analytic solution.
Diffusion approximation for modeling of 3-D radiation distributions
Zardecki, A.; Gerstl, S.A.W.; De Kinder, R.E. Jr.
1985-01-01
A three-dimensional transport code DIF3D, based on the diffusion approximation, is used to model the spatial distribution of radiation energy arising from volumetric isotropic sources. Future work will be concerned with the determination of irradiances and modeling of realistic scenarios, relevant to the battlefield conditions. 8 refs., 4 figs.
Energy Science and Technology Software Center (ESTSC)
1985-10-10
MARCOPOLO calculates the radial and axial diffusion coefficients in one-group and multi-group theory for a cylinderized cell (Wigner-Seitz theory) with several concentric zones according to the isotropic shock or linear anisotropic shock hypotheses.
GPU-accelerated 3D neutron diffusion code based on finite difference method
Xu, Q.; Yu, G.; Wang, K.
2012-07-01
Finite difference method, as a traditional numerical solution to neutron diffusion equation, although considered simpler and more precise than the coarse mesh nodal methods, has a bottle neck to be widely applied caused by the huge memory and unendurable computation time it requires. In recent years, the concept of General-Purpose computation on GPUs has provided us with a powerful computational engine for scientific research. In this study, a GPU-Accelerated multi-group 3D neutron diffusion code based on finite difference method was developed. First, a clean-sheet neutron diffusion code (3DFD-CPU) was written in C++ on the CPU architecture, and later ported to GPUs under NVIDIA's CUDA platform (3DFD-GPU). The IAEA 3D PWR benchmark problem was calculated in the numerical test, where three different codes, including the original CPU-based sequential code, the HYPRE (High Performance Pre-conditioners)-based diffusion code and CITATION, were used as counterpoints to test the efficiency and accuracy of the GPU-based program. The results demonstrate both high efficiency and adequate accuracy of the GPU implementation for neutron diffusion equation. A speedup factor of about 46 times was obtained, using NVIDIA's Geforce GTX470 GPU card against a 2.50 GHz Intel Quad Q9300 CPU processor. Compared with the HYPRE-based code performing in parallel on an 8-core tower server, the speedup of about 2 still could be observed. More encouragingly, without any mathematical acceleration technology, the GPU implementation ran about 5 times faster than CITATION which was speeded up by using the SOR method and Chebyshev extrapolation technique. (authors)
Simulate-HEX - The multi-group diffusion equation in hexagonal-z geometry
Lindahl, S. O.
2013-07-01
The multigroup diffusion equation is solved for the hexagonal-z geometry by dividing each hexagon into 6 triangles. In each triangle, the Fourier solution of the wave equation is approximated by 8 plane waves to describe the intra-nodal flux accurately. In the end an efficient Finite Difference like equation is obtained. The coefficients of this equation depend on the flux solution itself and they are updated once per power/void iteration. A numerical example demonstrates the high accuracy of the method. (authors)
Gifford, Kent A; Wareing, Todd A; Failla, Gregory; Horton, John L; Eifel, Patricia J; Mourtada, Firas
2010-01-01
A patient dose distribution was calculated by a 3D multi-group S N particle transport code for intracavitary brachytherapy of the cervix uteri and compared to previously published Monte Carlo results. A Cs-137 LDR intracavitary brachytherapy CT data set was chosen from our clinical database. MCNPX version 2.5.c, was used to calculate the dose distribution. A 3D multi-group S N particle transport code, Attila version 6.1.1 was used to simulate the same patient. Each patient applicator was built in SolidWorks, a mechanical design package, and then assembled with a coordinate transformation and rotation for the patient. The SolidWorks exported applicator geometry was imported into Attila for calculation. Dose matrices were overlaid on the patient CT data set. Dose volume histograms and point doses were compared. The MCNPX calculation required 14.8 hours, whereas the Attila calculation required 22.2 minutes on a 1.8 GHz AMD Opteron CPU. Agreement between Attila and MCNPX dose calculations at the ICRU 38 points was within +/- 3%. Calculated doses to the 2 cc and 5 cc volumes of highest dose differed by not more than +/- 1.1% between the two codes. Dose and DVH overlays agreed well qualitatively. Attila can calculate dose accurately and efficiently for this Cs-137 CT-based patient geometry. Our data showed that a three-group cross-section set is adequate for Cs-137 computations. Future work is aimed at implementing an optimized version of Attila for radiotherapy calculations. PMID:20160682
Converged accelerated finite difference scheme for the multigroup neutron diffusion equation
Terranova, N.; Mostacci, D.; Ganapol, B. D.
2013-07-01
Computer codes involving neutron transport theory for nuclear engineering applications always require verification to assess improvement. Generally, analytical and semi-analytical benchmarks are desirable, since they are capable of high precision solutions to provide accurate standards of comparison. However, these benchmarks often involve relatively simple problems, usually assuming a certain degree of abstract modeling. In the present work, we show how semi-analytical equivalent benchmarks can be numerically generated using convergence acceleration. Specifically, we investigate the error behavior of a 1D spatial finite difference scheme for the multigroup (MG) steady-state neutron diffusion equation in plane geometry. Since solutions depending on subsequent discretization can be envisioned as terms of an infinite sequence converging to the true solution, extrapolation methods can accelerate an iterative process to obtain the limit before numerical instability sets in. The obtained results have been compared to the analytical solution to the 1D multigroup diffusion equation when available, using FORTRAN as the computational language. Finally, a slowing down problem has been solved using a cascading source update, showing how a finite difference scheme performs for ultra-fine groups (104 groups) in a reasonable computational time using convergence acceleration. (authors)
NASA Technical Reports Server (NTRS)
Collier, G.
1967-01-01
Computer program VARI-QUIR 3 provides Gauss-Seidel type of solution with inner and outer iterations for steady-state, multigroup, two-dimensional neutron diffusion equations. The program has no restrictions on any of the input parameters such as the number of groups, regions, or materials.
Olson, Gordon L.
2011-08-20
Highlights: {yields} An existing multigroup transport algorithm is extended to be second-order in time. {yields} A new algorithm is presented that does not require a grey acceleration solution. {yields} The two algorithms are tested with 2D, multi-material problems. {yields} The two algorithms have comparable computational requirements. - Abstract: An existing solution method for solving the multigroup radiation equations, linear multifrequency-grey acceleration, is here extended to be second order in time. This method works for simple diffusion and for flux-limited diffusion, with or without material conduction. A new method is developed that does not require the solution of an averaged grey transport equation. It is effective solving both the diffusion and P{sub 1} forms of the transport equation. Two dimensional, multi-material test problems are used to compare the solution methods.
Hasan, M.Z.
1986-07-01
FENAT solves the two-dimensional energy dependent diffusion equation in Cartesian (X-Y) and cylindrical/toroidal (R-Z) coordinates. The boundary conditions allowed are: vacuum, reflection, albedo and surface source. The energy variable is treated by multigroup method. The resulting multigroup diffusion equation is solved by finite element Galerkin's method with triangular element discretization of the spatial domain. The algebraic matrix equation is solved by the direct method of Crout variation of Gauss' elimination. Dynamic memory allocation has been used so that the maximum problem size is limited by the size of active core storage of the machine. When necessary, the global matrix is stored in a binary disk file. FENAT is particularly suitable for the transport of neutral atoms in fusion plasmas.
Modeling the diffusion of phosphorus in silicon in 3-D
Baker, K.R.
1994-12-31
The use of matrix preconditioning in semiconductor process simulation is examined. The simplified nonlinear single-species model for the diffusion of phosphorus into silicon is considered. The experimental three-dimensional simulator, PEPPER3, which uses finite differences and the numerical method of lines to implement the reaction-diffusion equation is modified to allow NSPCG to be called to solve the linear system in the inner Newton loop. Use of NSPCG allowed various accelerators such as Generalized Minimal Residual (GMRES) and Conjugate Gradient (CG) to be used in conjunction with preconditioners such as Richardson, Jacobi, and Incomplete Cholesky.
VENTURE/PC manual: A multidimensional multigroup neutron diffusion code system. Version 3
Shapiro, A.; Huria, H.C.; Cho, K.W.
1991-12-01
VENTURE/PC is a recompilation of part of the Oak Ridge BOLD VENTURE code system, which will operate on an IBM PC or compatible computer. Neutron diffusion theory solutions are obtained for multidimensional, multigroup problems. This manual contains information associated with operating the code system. The purpose of the various modules used in the code system, and the input for these modules are discussed. The PC code structure is also given. Version 2 included several enhancements not given in the original version of the code. In particular, flux iterations can be done in core rather than by reading and writing to disk, for problems which allow sufficient memory for such in-core iterations. This speeds up the iteration process. Version 3 does not include any of the special processors used in the previous versions. These special processors utilized formatted input for various elements of the code system. All such input data is now entered through the Input Processor, which produces standard interface files for the various modules in the code system. In addition, a Standard Interface File Handbook is included in the documentation which is distributed with the code, to assist in developing the input for the Input Processor.
VENTURE/PC manual: A multidimensional multigroup neutron diffusion code system
Shapiro, A.; Huria, H.C.; Cho, K.W. )
1991-12-01
VENTURE/PC is a recompilation of part of the Oak Ridge BOLD VENTURE code system, which will operate on an IBM PC or compatible computer. Neutron diffusion theory solutions are obtained for multidimensional, multigroup problems. This manual contains information associated with operating the code system. The purpose of the various modules used in the code system, and the input for these modules are discussed. The PC code structure is also given. Version 2 included several enhancements not given in the original version of the code. In particular, flux iterations can be done in core rather than by reading and writing to disk, for problems which allow sufficient memory for such in-core iterations. This speeds up the iteration process. Version 3 does not include any of the special processors used in the previous versions. These special processors utilized formatted input for various elements of the code system. All such input data is now entered through the Input Processor, which produces standard interface files for the various modules in the code system. In addition, a Standard Interface File Handbook is included in the documentation which is distributed with the code, to assist in developing the input for the Input Processor.
3D parameter reconstruction in hyperspectral diffuse optical tomography
NASA Astrophysics Data System (ADS)
Saibaba, Arvind K.; Krishnamurthy, Nishanth; Anderson, Pamela G.; Kainerstorfer, Jana M.; Sassaroli, Angelo; Miller, Eric L.; Fantini, Sergio; Kilmer, Misha E.
2015-03-01
The imaging of shape perturbation and chromophore concentration using Diffuse Optical Tomography (DOT) data can be mathematically described as an ill-posed and non-linear inverse problem. The reconstruction algorithm for hyperspectral data using a linearized Born model is prohibitively expensive, both in terms of computation and memory. We model the shape of the perturbation using parametric level-set approach (PaLS). We discuss novel computational strategies for reducing the computational cost based on a Krylov subspace approach for parameteric linear systems and a compression strategy for the parameter-to-observation map. We will demonstrate the validity of our approach by comparison with experiments.
3D shape measurements for non-diffusive objects using fringe projection techniques
NASA Astrophysics Data System (ADS)
Su, Wei-Hung; Tseng, Bae-Heng; Cheng, Nai-Jen
2013-09-01
A scanning approach using holographic techniques to perform the 3D shape measurement for a non-diffusive object is proposed. Even though the depth discontinuity on the inspected surface is pretty high, the proposed method can retrieve the 3D shape precisely.
Intelligent speckle reducing anisotropic diffusion algorithm for automated 3-D ultrasound images.
Wu, Jun; Wang, Yuanyuan; Yu, Jinhua; Shi, Xinling; Zhang, Junhua; Chen, Yue; Pang, Yun
2015-02-01
A novel 3-D filtering method is presented for speckle reduction and detail preservation in automated 3-D ultrasound images. First, texture features of an image are analyzed by using the improved quadtree (QT) decomposition. Then, the optimal homogeneous and the obvious heterogeneous regions are selected from QT decomposition results. Finally, diffusion parameters and diffusion process are automatically decided based on the properties of these two selected regions. The computing time needed for 2-D speckle reduction is very short. However, the computing time required for 3-D speckle reduction is often hundreds of times longer than 2-D speckle reduction. This may limit its potential application in practice. Because this new filter can adaptively adjust the time step of iteration, the computation time is reduced effectively. Both synthetic and real 3-D ultrasound images are used to evaluate the proposed filter. It is shown that this filter is superior to other methods in both practicality and efficiency. PMID:26366596
FPGA-based real-time anisotropic diffusion filtering of 3D ultrasound images
NASA Astrophysics Data System (ADS)
Castro-Pareja, Carlos R.; Dandekar, Omkar S.; Shekhar, Raj
2005-02-01
Three-dimensional ultrasonic imaging, especially the emerging real-time version of it, is particularly valuable in medical applications such as echocardiography, obstetrics and surgical navigation. A known problem with ultrasound images is their high level of speckle noise. Anisotropic diffusion filtering has been shown to be effective in enhancing the visual quality of 3D ultrasound images and as preprocessing prior to advanced image processing. However, due to its arithmetic complexity and the sheer size of 3D ultrasound images, it is not possible to perform online, real-time anisotropic diffusion filtering using standard software implementations. We present an FPGA-based architecture that allows performing anisotropic diffusion filtering of 3D images at acquisition rates, thus enabling the use of this filtering technique in real-time applications, such as visualization, registration and volume rendering.
Wang, Y.
2013-07-01
Nonlinear diffusion acceleration (NDA) can improve the performance of a neutron transport solver significantly especially for the multigroup eigenvalue problems. The high-order transport equation and the transport-corrected low-order diffusion equation form a nonlinear system in NDA, which can be solved via a Picard iteration. The consistency of the correction of the low-order equation is important to ensure the stabilization and effectiveness of the iteration. It also makes the low-order equation preserve the scalar flux of the high-order equation. In this paper, the consistent correction for a particular discretization scheme, self-adjoint angular flux (SAAF) formulation with discrete ordinates method (S{sub N}) and continuous finite element method (CFEM) is proposed for the multigroup neutron transport equation. Equations with the anisotropic scatterings and a void treatment are included. The Picard iteration with this scheme has been implemented and tested with RattleS{sub N}ake, a MOOSE-based application at INL. Convergence results are presented. (authors)
Shestakov, Aleksei I. Offner, Stella S.R.
2008-01-10
We present a scheme to solve the nonlinear multigroup radiation diffusion (MGD) equations. The method is incorporated into a massively parallel, multidimensional, Eulerian radiation-hydrodynamic code with Adaptive Mesh Refinement (AMR). The patch-based AMR algorithm refines in both space and time creating a hierarchy of levels, coarsest to finest. The physics modules are time-advanced using operator splitting. On each level, separate 'level-solve' packages advance the modules. Our multigroup level-solve adapts an implicit procedure which leads to a two-step iterative scheme that alternates between elliptic solves for each group with intra-cell group coupling. For robustness, we introduce pseudo transient continuation ({psi}tc). We analyze the magnitude of the {psi}tc parameter to ensure positivity of the resulting linear system, diagonal dominance and convergence of the two-step scheme. For AMR, a level defines a subdomain for refinement. For diffusive processes such as MGD, the refined level uses Dirichlet boundary data at the coarse-fine interface and the data is derived from the coarse level solution. After advancing on the fine level, an additional procedure, the sync-solve (SS), is required in order to enforce conservation. The MGD SS reduces to an elliptic solve on a combined grid for a system of G equations, where G is the number of groups. We adapt the 'partial temperature' scheme for the SS; hence, we reuse the infrastructure developed for scalar equations. Results are presented. We consider a multigroup test problem with a known analytic solution. We demonstrate utility of {psi}tc by running with increasingly larger timesteps. Lastly, we simulate the sudden release of energy Y inside an Al sphere (r = 15 cm) suspended in air at STP. For Y = 11 kT, we find that gray radiation diffusion and MGD produce similar results. However, if Y = 1 MT, the two packages yield different results. Our large Y simulation contradicts a long-standing theory and demonstrates
Shestakov, A I; Offner, S R
2007-03-02
We present a scheme to solve the nonlinear multigroup radiation diffusion (MGD) equations. The method is incorporated into a massively parallel, multidimensional, Eulerian radiation-hydrodynamic code with adaptive mesh refinement (AMR). The patch-based AMR algorithm refines in both space and time creating a hierarchy of levels, coarsest to finest. The physics modules are time-advanced using operator splitting. On each level, separate 'level-solve' packages advance the modules. Our multigroup level-solve adapts an implicit procedure which leads to a two-step iterative scheme that alternates between elliptic solves for each group with intra-cell group coupling. For robustness, we introduce pseudo transient continuation ({Psi}tc). We analyze the magnitude of the {Psi}tc parameter to ensure positivity of the resulting linear system, diagonal dominance and convergence of the two-step scheme. For AMR, a level defines a subdomain for refinement. For diffusive processes such as MGD, the refined level uses Dirichet boundary data at the coarse-fine interface and the data is derived from the coarse level solution. After advancing on the fine level, an additional procedure, the sync-solve (SS), is required in order to enforce conservation. The MGD SS reduces to an elliptic solve on a combined grid for a system of G equations, where G is the number of groups. We adapt the 'partial temperature' scheme for the SS; hence, we reuse the infrastructure developed for scalar equations. Results are presented. We consider a multigroup test problem with a known analytic solution. We demonstrate utility of {Psi}tc by running with increasingly larger timesteps. Lastly, we simulate the sudden release of energy Y inside an Al sphere (r = 15 cm) suspended in air at STP. For Y = 11 kT, we find that gray radiation diffusion and MGD produce similar results. However, if Y = 1 MT, the two packages yield different results. Our large Y simulation contradicts a long-standing theory and demonstrates
NASA Astrophysics Data System (ADS)
Shestakov, Aleksei I.; Offner, Stella S. R.
2008-01-01
We present a scheme to solve the nonlinear multigroup radiation diffusion (MGD) equations. The method is incorporated into a massively parallel, multidimensional, Eulerian radiation-hydrodynamic code with Adaptive Mesh Refinement (AMR). The patch-based AMR algorithm refines in both space and time creating a hierarchy of levels, coarsest to finest. The physics modules are time-advanced using operator splitting. On each level, separate "level-solve" packages advance the modules. Our multigroup level-solve adapts an implicit procedure which leads to a two-step iterative scheme that alternates between elliptic solves for each group with intra-cell group coupling. For robustness, we introduce pseudo transient continuation (Ψtc). We analyze the magnitude of the Ψtc parameter to ensure positivity of the resulting linear system, diagonal dominance and convergence of the two-step scheme. For AMR, a level defines a subdomain for refinement. For diffusive processes such as MGD, the refined level uses Dirichlet boundary data at the coarse-fine interface and the data is derived from the coarse level solution. After advancing on the fine level, an additional procedure, the sync-solve (SS), is required in order to enforce conservation. The MGD SS reduces to an elliptic solve on a combined grid for a system of G equations, where G is the number of groups. We adapt the "partial temperature" scheme for the SS; hence, we reuse the infrastructure developed for scalar equations. Results are presented. We consider a multigroup test problem with a known analytic solution. We demonstrate utility of Ψtc by running with increasingly larger timesteps. Lastly, we simulate the sudden release of energy Y inside an Al sphere (r = 15 cm) suspended in air at STP. For Y = 11 kT, we find that gray radiation diffusion and MGD produce similar results. However, if Y = 1 MT, the two packages yield different results. Our large Y simulation contradicts a long-standing theory and demonstrates the
Shestakov, A I; Offner, S R
2006-09-21
We present a scheme to solve the nonlinear multigroup radiation diffusion (MGD) equations. The method is incorporated into a massively parallel, multidimensional, Eulerian radiation-hydrodynamic code with adaptive mesh refinement (AMR). The patch-based AMR algorithm refines in both space and time creating a hierarchy of levels, coarsest to finest. The physics modules are time-advanced using operator splitting. On each level, separate 'level-solve' packages advance the modules. Our multigroup level-solve adapts an implicit procedure which leads to a two-step iterative scheme that alternates between elliptic solves for each group with intra-cell group coupling. For robustness, we introduce pseudo transient continuation ({Psi}tc). We analyze the magnitude of the {Psi}tc parameter to ensure positivity of the resulting linear system, diagonal dominance and convergence of the two-step scheme. For AMR, a level defines a subdomain for refinement. For diffusive processes such as MGD, the refined level uses Dirichet boundary data at the coarse-fine interface and the data is derived from the coarse level solution. After advancing on the fine level, an additional procedure, the sync-solve (SS), is required in order to enforce conservation. The MGD SS reduces to an elliptic solve on a combined grid for a system of G equations, where G is the number of groups. We adapt the 'partial temperature' scheme for the SS; hence, we reuse the infrastructure developed for scalar equations. Results are presented. We consider a multigroup test problem with a known analytic solution. We demonstrate utility of {Psi}tc by running with increasingly larger timesteps. Lastly, we simulate the sudden release of energy Y inside an Al sphere (r = 15 cm) suspended in air at STP. For Y = 11 kT, we find that gray radiation diffusion and MGD produce similar results. However, if Y = 1 MT, the two packages yield different results. Our large Y simulation contradicts a long-standing theory and demonstrates
Howison, Mark
2010-05-06
We compare the performance of hand-tuned CUDA implementations of bilateral and anisotropic diffusion filters for denoising 3D MRI datasets. Our tests sweep comparable parameters for the two filters and measure total runtime, memory bandwidth, computational throughput, and mean squared errors relative to a noiseless reference dataset.
3D hydrodynamic interactions lead to divergences in 2D diffusion.
Bleibel, Johannes; Domínguez, Alvaro; Oettel, Martin
2015-05-20
We investigate the influence of 3D hydrodynamic interactions on confined colloidal suspensions, where only the colloids are restricted to one or two dimensions. In the absence of static interactions among the colloids, i.e., an ideal gas of colloidal particles with a finite hydrodynamic radius, we find a divergent collective diffusion coefficient. The origin of the divergence is traced back to the dimensional mismatch of 3D hydrodynamic interactions and the colloidal particles moving only in 1D or 2D. Our results from theory are confirmed by Stokesian dynamics simulations and supported by light scattering observational data for particles at a fluid interface. PMID:25923320
3D hydrodynamic interactions lead to divergences in 2D diffusion
NASA Astrophysics Data System (ADS)
Bleibel, Johannes; Domínguez, Alvaro; Oettel, Martin
2015-05-01
We investigate the influence of 3D hydrodynamic interactions on confined colloidal suspensions, where only the colloids are restricted to one or two dimensions. In the absence of static interactions among the colloids, i.e., an ideal gas of colloidal particles with a finite hydrodynamic radius, we find a divergent collective diffusion coefficient. The origin of the divergence is traced back to the dimensional mismatch of 3D hydrodynamic interactions and the colloidal particles moving only in 1D or 2D. Our results from theory are confirmed by Stokesian dynamics simulations and supported by light scattering observational data for particles at a fluid interface.
Motion-induced phase error estimation and correction in 3D diffusion tensor imaging.
Van, Anh T; Hernando, Diego; Sutton, Bradley P
2011-11-01
A multishot data acquisition strategy is one way to mitigate B0 distortion and T2∗ blurring for high-resolution diffusion-weighted magnetic resonance imaging experiments. However, different object motions that take place during different shots cause phase inconsistencies in the data, leading to significant image artifacts. This work proposes a maximum likelihood estimation and k-space correction of motion-induced phase errors in 3D multishot diffusion tensor imaging. The proposed error estimation is robust, unbiased, and approaches the Cramer-Rao lower bound. For rigid body motion, the proposed correction effectively removes motion-induced phase errors regardless of the k-space trajectory used and gives comparable performance to the more computationally expensive 3D iterative nonlinear phase error correction method. The method has been extended to handle multichannel data collected using phased-array coils. Simulation and in vivo data are shown to demonstrate the performance of the method. PMID:21652284
Multilayer Spheroids To Quantify Drug Uptake and Diffusion in 3D
2015-01-01
There is a need for new quantitative in vitro models of drug uptake and diffusion to help assess drug toxicity/efficacy as well as new more predictive models for drug discovery. We report a three-dimensional (3D) multilayer spheroid model and a new algorithm to quantitatively study uptake and inward diffusion of fluorescent calcein via gap junction intercellular communication (GJIC). When incubated with calcein-AM, a substrate of the efflux transporter P-glycoprotein (Pgp), spheroids from a variety of cell types accumulated calcein over time. Accumulation decreased in spheroids overexpressing Pgp (HEK-MDR) and was increased in the presence of Pgp inhibitors (verapamil, loperamide, cyclosporin A). Inward diffusion of calcein was negligible in spheroids that lacked GJIC (OVCAR-3, SK-OV-3) and was reduced in the presence of an inhibitor of GJIC (carbenoxolone). In addition to inhibiting Pgp, verapamil and loperamide, but not cyclosporin A, inhibited inward diffusion of calcein, suggesting that they also inhibit GJIC. The dose response curves of verapamil’s inhibition of Pgp and GJIC were similar (IC50: 8 μM). The method is amenable to many different cell types and may serve as a quantitative 3D model that more accurately replicates in vivo barriers to drug uptake and diffusion. PMID:24641346
Planar Gradient Diffusion System to Investigate Chemotaxis in a 3D Collagen Matrix.
Stout, David A; Toyjanova, Jennet; Franck, Christian
2015-01-01
The importance of cell migration can be seen through the development of human life. When cells migrate, they generate forces and transfer these forces to their surrounding area, leading to cell movement and migration. In order to understand the mechanisms that can alter and/or affect cell migration, one can study these forces. In theory, understanding the fundamental mechanisms and forces underlying cell migration holds the promise of effective approaches for treating diseases and promoting cellular transplantation. Unfortunately, modern chemotaxis chambers that have been developed are usually restricted to two dimensions (2D) and have complex diffusion gradients that make the experiment difficult to interpret. To this end, we have developed, and describe in this paper, a direct-viewing chamber for chemotaxis studies, which allows one to overcome modern chemotaxis chamber obstacles able to measure cell forces and specific concentration within the chamber in a 3D environment to study cell 3D migration. More compelling, this approach allows one to successfully model diffusion through 3D collagen matrices and calculate the coefficient of diffusion of a chemoattractant through multiple different concentrations of collagen, while keeping the system simple and user friendly for traction force microscopy (TFM) and digital volume correlation (DVC) analysis. PMID:26131645
3D MRI of non-Gaussian 3He gas diffusion in the rat lung
NASA Astrophysics Data System (ADS)
Jacob, Richard E.; Laicher, Gernot; Minard, Kevin R.
2007-10-01
In 3He magnetic resonance images of pulmonary air spaces, the confining architecture of the parenchymal tissue results in a non-Gaussian distribution of signal phase that non-exponentially attenuates image intensity as diffusion weighting is increased. Here, two approaches previously used for the analysis of non-Gaussian effects in the lung are compared and related using diffusion-weighted 3He MR images of mechanically ventilated rats. One approach is model-based and was presented by Yablonskiy et al., while the other approach utilizes the second order decay contribution that is predicted from the cumulant expansion theorem. Total lung coverage is achieved using a hybrid 3D pulse sequence that combines conventional phase encoding with sparse radial sampling for efficient gas usage. This enables the acquisition of nine 3D images using a total of only ˜1 L of hyperpolarized 3He gas. Diffusion weighting ranges from 0 s/cm 2 to 40 s/cm 2. Results show that the non-Gaussian effects of 3He gas diffusion in healthy rat lungs are directly attributed to the anisotropic geometry of lung microstructure as predicted by the Yablonskiy model, and that quantitative analysis over the entire lung can be reliably repeated in time-course studies of the same animal.
Proton exchange membrane micro fuel cells on 3D porous silicon gas diffusion layers
NASA Astrophysics Data System (ADS)
Kouassi, S.; Gautier, G.; Thery, J.; Desplobain, S.; Borella, M.; Ventura, L.; Laurent, J.-Y.
2012-10-01
Since the 90's, porous silicon has been studied and implemented in many devices, especially in MEMS technology. In this article, we present a new approach to build miniaturized proton exchange membrane micro-fuel cells using porous silicon as a hydrogen diffusion layer. In particular, we propose an innovative process to build micro fuel cells from a “corrugated iron like” 3D structured porous silicon substrates. This structure is able to increase up to 40% the cell area keeping a constant footprint on the silicon wafer. We propose here a process route to perform electrochemically 3D porous gas diffusion layers and to deposit fuel cell active layers on such substrates. The prototype peak power performance was measured to be 90 mW cm-2 in a “breathing configuration” at room temperature. These performances are less than expected if we compare with a reference 2D micro fuel cell. Actually, the active layer deposition processes are not fully optimized but this prototype demonstrates the feasibility of these 3D devices.
Fractality in the neuron axonal topography of the human brain based on 3-D diffusion MRI
NASA Astrophysics Data System (ADS)
Katsaloulis, P.; Ghosh, A.; Philippe, A. C.; Provata, A.; Deriche, R.
2012-05-01
In this work the fractal architecture of the neuron axonal topography of the human brain is evaluated, as derived from 3-D diffusion MRI (dMRI) acquisitions. This is a 3D extension of work performed previously in 2D regions of interest (ROIs), where the fractal dimension of the neuron axonal topography was computed from dMRI data. A group study with 18 subjects is here conducted and the fractal dimensions D f of the entire 3-D volume of the brains is estimated via the box counting, the correlation dimension and the fractal mass dimension methods. The neuron axon data is obtained using tractography algorithms on diffusion tensor imaging of the brain. We find that all three calculations of D f give consistent results across subjects, namely, they demonstrate fractal characteristics in the short and medium length scales: different fractal exponents prevail at different length scales, an indication of multifractality. We surmise that this complexity stems as a collective property emerging when many local brain units, performing different functional tasks and having different local topologies, are recorded together.
Double depth-enhanced 3D integral imaging in projection-type system without diffuser
NASA Astrophysics Data System (ADS)
Zhang, Lei; Jiao, Xiao-xue; Sun, Yu; Xie, Yan; Liu, Shao-peng
2015-05-01
Integral imaging is a three dimensional (3D) display technology without any additional equipment. A new system is proposed in this paper which consists of the elemental images of real images in real mode (RIRM) and the ones of virtual images in real mode (VIRM). The real images in real mode are the same as the conventional integral images. The virtual images in real mode are obtained by changing the coordinates of the corresponding points in elemental images which can be reconstructed by the lens array in virtual space. In order to reduce the spot size of the reconstructed images, the diffuser in conventional integral imaging is given up in the proposed method. Then the spot size is nearly 1/20 of that in the conventional system. And an optical integral imaging system is constructed to confirm that our proposed method opens a new way for the application of the passive 3D display technology.
Parametric estimation of 3D tubular structures for diffuse optical tomography
Larusson, Fridrik; Anderson, Pamela G.; Rosenberg, Elizabeth; Kilmer, Misha E.; Sassaroli, Angelo; Fantini, Sergio; Miller, Eric L.
2013-01-01
We explore the use of diffuse optical tomography (DOT) for the recovery of 3D tubular shapes representing vascular structures in breast tissue. Using a parametric level set method (PaLS) our method incorporates the connectedness of vascular structures in breast tissue to reconstruct shape and absorption values from severely limited data sets. The approach is based on a decomposition of the unknown structure into a series of two dimensional slices. Using a simplified physical model that ignores 3D effects of the complete structure, we develop a novel inter-slice regularization strategy to obtain global regularity. We report on simulated and experimental reconstructions using realistic optical contrasts where our method provides a more accurate estimate compared to an unregularized approach and a pixel based reconstruction. PMID:23411913
Schilling, Kurt; Janve, Vaibhav; Gao, Yurui; Stepniewska, Iwona; Landman, Bennett A; Anderson, Adam W
2016-04-01
The ability of diffusion MRI (dMRI) fiber tractography to non-invasively map three-dimensional (3D) anatomical networks in the human brain has made it a valuable tool in both clinical and research settings. However, there are many assumptions inherent to any tractography algorithm that can limit the accuracy of the reconstructed fiber tracts. Among them is the assumption that the diffusion-weighted images accurately reflect the underlying fiber orientation distribution (FOD) in the MRI voxel. Consequently, validating dMRI's ability to assess the underlying fiber orientation in each voxel is critical for its use as a biomedical tool. Here, using post-mortem histology and confocal microscopy, we present a method to perform histological validation of orientation functions in 3D, which has previously been limited to two-dimensional analysis of tissue sections. We demonstrate the ability to extract the 3D FOD from confocal z-stacks, and quantify the agreement between the MRI estimates of orientation information obtained using constrained spherical deconvolution (CSD) and the true geometry of the fibers. We find an orientation error of approximately 6° in voxels containing nearly parallel fibers, and 10-11° in crossing fiber regions, and note that CSD was unable to resolve fibers crossing at angles below 60° in our dataset. This is the first time that the 3D white matter orientation distribution is calculated from histology and compared to dMRI. Thus, this technique serves as a gold standard for dMRI validation studies - providing the ability to determine the extent to which the dMRI signal is consistent with the histological FOD, and to establish how well different dMRI models can predict the ground truth FOD. PMID:26804781
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
Frank, Lawrence R.; Jung, Youngkyoo; Inati, Souheil; Tyszka, J. Michael; Wong, Eric C.
2009-01-01
We present an acquisition and reconstruction method designed to acquire high resolution 3D fast spin echo diffusion tensor images while mitigating the major sources of artifacts in DTI - field distortions, eddy currents and motion. The resulting images, being 3D, are of high SNR, and being fast spin echoes, exhibit greatly reduced field distortions. This sequence utilizes variable density spiral acquisition gradients, which allow for the implementation of a self-navigation scheme by which both eddy current and motion artifacts are removed. The result is that high resolution 3D DTI images are produced without the need for eddy current compensating gradients or B0 field correction. In addition, a novel method for fast and accurate reconstruction of the non-Cartesian data is employed. Results are demonstrated in the brains of normal human volunteers. PMID:19778618
Diffuse reflectance optical topography: location of inclusions in 3D and detectability limits
Carbone, N. A.; Baez, G. R.; García, H. A.; Waks Serra, M. V.; Di Rocco, H. O.; Iriarte, D. I.; Pomarico, J. A.; Grosenick, D.; Macdonald, R.
2014-01-01
In the present contribution we investigate the images of CW diffusely reflected light for a point-like source, registered by a CCD camera imaging a turbid medium containing an absorbing lesion. We show that detection of μa variations (absorption anomalies) is achieved if images are normalized to background intensity. A theoretical analysis based on the diffusion approximation is presented to investigate the sensitivity and the limitations of our proposal and a novel procedure to find the location of the inclusions in 3D is given and tested. An analysis of the noise and its influence on the detection capabilities of our proposal is provided. Experimental results on phantoms are also given, supporting the proposed approach. PMID:24876999
Fast and Robust Sixth Order Multigrid Computation for 3D Convection Diffusion Equation.
Wang, Yin; Zhang, Jun
2010-10-15
We present a sixth order explicit compact finite difference scheme to solve the three dimensional (3D) convection diffusion equation. We first use multiscale multigrid method to solve the linear systems arising from a 19-point fourth order discretization scheme to compute the fourth order solutions on both the coarse grid and the fine grid. Then an operator based interpolation scheme combined with an extrapolation technique is used to approximate the sixth order accurate solution on the fine grid. Since the multigrid method using a standard point relaxation smoother may fail to achieve the optimal grid independent convergence rate for solving convection diffusion equation with a high Reynolds number, we implement the plane relaxation smoother in the multigrid solver to achieve better grid independency. Supporting numerical results are presented to demonstrate the efficiency and accuracy of the sixth order compact scheme (SOC), compared with the previously published fourth order compact scheme (FOC). PMID:21151737
Ex Vivo 3D Diffusion Tensor Imaging and Quantification of Cardiac Laminar Structure
Helm, Patrick A.; Tseng, Hsiang-Jer; Younes, Laurent; McVeigh, Elliot R.; Winslow, Raimond L.
2007-01-01
A three-dimensional (3D) diffusion-weighted imaging (DWI) method for measuring cardiac fiber structure at high spatial resolution is presented. The method was applied to the ex vivo reconstruction of the fiber architecture of seven canine hearts. A novel hypothesis-testing method was developed and used to show that distinct populations of secondary and tertiary eigenvalues may be distinguished at reasonable confidence levels (P ≤ 0.01) within the canine ventricle. Fiber inclination and sheet angles are reported as a function of transmural depth through the anterior, lateral, and posterior left ventricle (LV) free wall. Within anisotropic regions, two consistent and dominant orientations were identified, supporting published results from histological studies and providing strong evidence that the tertiary eigenvector of the diffusion tensor (DT) defines the sheet normal. PMID:16149057
Kafieh, Raheleh; Rabbani, Hossein; Abramoff, Michael D.; Sonka, Milan
2013-01-01
Optical coherence tomography (OCT) is a powerful and noninvasive method for retinal imaging. In this paper, we introduce a fast segmentation method based on a new variant of spectral graph theory named diffusion maps. The research is performed on spectral domain (SD) OCT images depicting macular and optic nerve head appearance. The presented approach does not require edge-based image information in localizing most of boundaries and relies on regional image texture. Consequently, the proposed method demonstrates robustness in situations of low image contrast or poor layer-to-layer image gradients. Diffusion mapping applied to 2D and 3D OCT datasets is composed of two steps, one for partitioning the data into important and less important sections, and another one for localization of internal layers. In the first step, the pixels/voxels are grouped in rectangular/cubic sets to form a graph node. The weights of the graph are calculated based on geometric distances between pixels/voxels and differences of their mean intensity. The first diffusion map clusters the data into three parts, the second of which is the area of interest. The other two sections are eliminated from the remaining calculations. In the second step, the remaining area is subjected to another diffusion map assessment and the internal layers are localized based on their textural similarities. The proposed method was tested on 23 datasets from two patient groups (glaucoma and normals). The mean unsigned border positioning errors (mean ± SD) was 8.52 ± 3.13 and 7.56 ± 2.95 μm for the 2D and 3D methods, respectively. PMID:23837966
3D modeling for solving forward model of no-contact fluorescence diffuse optical tomography method
NASA Astrophysics Data System (ADS)
Nouizi, F.; Chabrier, R.; Torregrossa, M.; Poulet, P.
2009-07-01
This paper presents detailed computational aspects of a new 3D modeling for solving the direct problem in a no-contact time-resolved Fluorescent Diffuse Optical Tomography (FDOT) method that rely on near-infrared scattered and fluorescent photons to image the optical properties and distribution of fluorescent probes in small laboratory animals. An optical scanner allowing performing in-vivo measurements in no-contact scheme was built in our laboratory and is presented. We use the three-dimensional Finite Element Method (FEM) to solve the coupled diffusion equations of excitation and fluorescence photons in highly scattering objects. The computed results allowed yielding photon density maps and the temporal profiles of photons on the surface of the small animal. Our 3D modeling of propagation of photons in the void space between the surface of the object and the detectors allows calculating the quantity of photons reaching the optodes. Simulations were carried-out on two test objects: a resin cylinder and a mouse phantom. The results demonstrate the potential applications of the method to pre-clinical imaging.
3D choroid neovascularization growth prediction based on reaction-diffusion model
NASA Astrophysics Data System (ADS)
Zhu, Shuxia; Chen, Xinjian; Shi, Fei; Xiang, Dehui; Zhu, Weifang; Chen, Haoyu
2016-03-01
Choroid neovascularization (CNV) is a kind of pathology from the choroid and CNV-related disease is one important cause of vision loss. It is desirable to predict the CNV growth rate so that appropriate treatment can be planned. In this paper, we seek to find a method to predict the growth of CNV based on 3D longitudinal Optical Coherence Tomography (OCT) images. A reaction-diffusion model is proposed for prediction. The method consists of four phases: pre-processing, meshing, CNV growth modeling and prediction. We not only apply the reaction-diffusion model to the disease region, but also take the surrounding tissues into consideration including outer retinal layer, inner retinal layer and choroid layer. The diffusion in these tissues is considered as isotropic. The finite-element-method (FEM) is used to solve the partial differential equations (PDE) in the diffusion model. The curve of CNV growth with treatment are fitted and then we can predict the CNV status in a future time point. The preliminary results demonstrated that our proposed method is accurate and the validity and feasibility of our model is obvious.
Confocal fluorometer for diffusion tracking in 3D engineered tissue constructs
NASA Astrophysics Data System (ADS)
Daly, D.; Zilioli, A.; Tan, N.; Buttenschoen, K.; Chikkanna, B.; Reynolds, J.; Marsden, B.; Hughes, C.
2016-03-01
We present results of the development of a non-contacting instrument, called fScan, based on scanning confocal fluorometry for assessing the diffusion of materials through a tissue matrix. There are many areas in healthcare diagnostics and screening where it is now widely accepted that the need for new quantitative monitoring technologies is a major pinch point in patient diagnostics and in vitro testing. With the increasing need to interpret 3D responses this commonly involves the need to track the diffusion of compounds, pharma-active species and cells through a 3D matrix of tissue. Methods are available but to support the advances that are currently only promised, this monitoring needs to be real-time, non-invasive, and economical. At the moment commercial meters tend to be invasive and usually require a sample of the medium to be removed and processed prior to testing. This methodology clearly has a number of significant disadvantages. fScan combines a fiber based optical arrangement with a compact, free space optical front end that has been integrated so that the sample's diffusion can be measured without interference. This architecture is particularly important due to the "wet" nature of the samples. fScan is designed to measure constructs located within standard well plates and a 2-D motion stage locates the required sample with respect to the measurement system. Results are presented that show how the meter has been used to evaluate movements of samples through collagen constructs in situ without disturbing their kinetic characteristics. These kinetics were little understood prior to these measurements.
Sensitivity of an asymmetric 3D diffuser to vortex-generator induced inlet condition perturbations
NASA Astrophysics Data System (ADS)
Grundmann, S.; Sayles, E. L.; Elkins, Christopher J.; Eaton, J. K.
2012-01-01
Modifications of the turbulent separated flow in an asymmetric three-dimensional diffuser due to inlet condition perturbations were investigated using conventional static pressure measurements and velocity data acquired using magnetic resonance velocimetry (MRV). Previous experiments and simulations revealed a strong sensitivity of the diffuser performance to weak secondary flows in the inlet. The present, more detailed experiments were conducted to obtain a better understanding of this sensitivity. Pressure data were acquired in an airflow apparatus at an inlet Reynolds number of 10,000. The diffuser pressure recovery was strongly affected by a pair of longitudinal vortices injected along one wall of the inlet channel using either dielectric barrier discharge plasma actuators or conventional half-delta wing vortex generators. MRV measurements were obtained in a water flow apparatus at matched Reynolds number for two different cases with passive vortex generators. The first case had a pair of counter-rotating longitudinal vortices embedded in the boundary layer near the center of the expanding wall of the diffuser such that the flow on the outsides of the vortices was directed toward the wall. The MRV data showed that the three-dimensional separation bubble initially grew much slower causing a rapid early reduction in the core flow velocity and a consequent reduction of total pressure losses due to turbulent mixing. This produced a 13% increase in the overall pressure recovery. For the second case, the vortices rotated in the opposite sense, and the image vortices pushed them into the corners. This led to a very rapid initial growth of the separation bubble and formation of strong swirl at the diffuser exit. These changes resulted in a 17% reduction in the overall pressure recovery for this case. The results emphasize the extreme sensitivity of 3D separated flows to weak perturbations.
NASA Astrophysics Data System (ADS)
Yu, Jhao-Ming; Chen, Liang-Yu; Pan, Min-Cheng; Hsu, Ya-Fen; Pan, Min-Chun
2015-03-01
Diffuse optical imaging (DOI) providing functional information of tissues has drawn great attention for the last two decades. Near infrared (NIR) DOI systems composed of scanning bench, opt-electrical measurement module, system control, and data processing and image reconstruction schemes are developed for the screening and diagnosis of breast tumors. Mostly, the scanning bench belonging to fixed source-and-detector configuration limits computed image resolution to an extent. To cope with the issue, we propose, design and implement a 3D prostrate ring-scanning equipment for NIR DOI with flexible combinations of illumination and detection, and with the function of radial, circular and vertical movement without hard compression of breast tissue like the imaging system using or incorporating with X-ray mammographic bench. Especially, a rotation-sliding-and-moving mechanism was designed for the guidance of source- and detection-channel movement. Following the previous justification for synthesized image reconstruction, in the paper the validation using varied phantoms is further conducted and 3D image reconstruction for their absorption and scattering coefficients is illustrated through the computation of our in-house coded schemes. The source and detection NIR data are acquired to reconstruct the 3D images through the operation of scanning bench in the movement of vertical, radial and circular directions. Rather than the fixed configuration, the addressed screening/diagnosing equipment has the flexibility for optical-channel expansion with a compromise among construction cost, operation time, and spatial resolution of reconstructed μa and μs' images.
Description of a parallel, 3D, finite element, hydrodynamics-diffusion code
Milovich, J L; Prasad, M K; Shestakov, A I
1999-04-11
We describe a parallel, 3D, unstructured grid finite element, hydrodynamic diffusion code for inertial confinement fusion (ICF) applications and the ancillary software used to run it. The code system is divided into two entities, a controller and a stand-alone physics code. The code system may reside on different computers; the controller on the user's workstation and the physics code on a supercomputer. The physics code is composed of separate hydrodynamic, equation-of-state, laser energy deposition, heat conduction, and radiation transport packages and is parallelized for distributed memory architectures. For parallelization, a SPMD model is adopted; the domain is decomposed into a disjoint collection of subdomains, one per processing element (PE). The PEs communicate using MPI. The code is used to simulate the hydrodynamic implosion of a spherical bubble.
Dynamic Implicit 3D Adaptive Mesh Refinement for Non-Equilibrium Radiation Diffusion
Philip, Bobby; Wang, Zhen; Berrill, Mark A; Rodriguez Rodriguez, Manuel; Pernice, Michael
2014-01-01
The time dependent non-equilibrium radiation diffusion equations are important for solving the transport of energy through radiation in optically thick regimes and find applications in several fields including astrophysics and inertial confinement fusion. The associated initial boundary value problems that are encountered exhibit a wide range of scales in space and time and are extremely challenging to solve. To efficiently and accurately simulate these systems we describe our research on combining techniques that will also find use more broadly for long term time integration of nonlinear multiphysics systems: implicit time integration for efficient long term time integration of stiff multiphysics systems, local control theory based step size control to minimize the required global number of time steps while controlling accuracy, dynamic 3D adaptive mesh refinement (AMR) to minimize memory and computational costs, Jacobian Free Newton Krylov methods on AMR grids for efficient nonlinear solution, and optimal multilevel preconditioner components that provide level independent linear solver convergence.
A CNN-based approach to integrate the 3-D turbolent diffusion equation
NASA Astrophysics Data System (ADS)
Nunnari, G.
2003-04-01
The paper deals with the integration of the 3-D turbulent diffusion equation. This problem is relevant in several application fields including fluid dynamics, air/water pollution, volcanic ash emissions and industrial hazard assessment. As it is well known numerical solution of such a kind of equation is very time consuming even by using modern digital computers and this represents a short-coming for on-line applications. To overcome this drawback a Cellular Neural Network Approach is proposed in this paper. CNN's proposed by Chua and Yang in 1988 are massive parallel analog non-linear circuits with local interconnections between the computing elements that allow very fast distributed computations. Nowadays several producers of semiconductors such as SGS-Thomson are producing on chip CNN's so that their massive use for heavy computing applications is expected in the near future. In the paper the methodological background of the proposed approach will be outlined. Further some results both in terms of accuracy and computation time will be presented also in comparison with traditional three-dimensional computation schemes. Some results obtained to model 3-D pollution problems in the industrial area of Siracusa (Italy), characterised by a large concentration of petrol-chemical plants, will be presented.
PFLOW: A 3-D Numerical Modeling Tool for Calculating Fluid-Pressure Diffusion from Coulomb Strain
NASA Astrophysics Data System (ADS)
Wolf, L. W.; Lee, M.; Meir, A.; Dyer, G.; Ma, K.; Chan, C.
2009-12-01
A new 3D time-dependent pore-pressure diffusion model PFLOW is developed to investigate the response of pore fluids to the crustal deformation generated by strong earthquakes in heterogeneous geologic media. Given crustal strain generated by changes in Coulomb stress, this MATLAB-based code uses Skempton's coefficient to calculate resulting changes fluid pressure. Pore-pressure diffusion can be tracked over time in a user-defined model space with user-prescribed Neumann or Dirchilet boundary conditions and with spatially variable values of permeability. PFLOW employs linear or quadratic finite elements for spatial discretization and first order or second order, explicit or implicit finite difference discretization in time. PFLOW is easily interfaced with output from deformation modeling programs such as Coulomb (Toda et al., 2007) or 3D-DEF (Gomberg and Ellis, 1994). The code is useful for investigating to first-order the evolution of pore pressure changes induced by changes in Coulomb stress and their possible relation to water-level changes in wells or changes in stream discharge. It can also be used for student research and classroom instruction. As an example application, we calculate the coseismic pore pressure changes and diffusion induced by volumetric strain associated with the 1999 Chi-Chi earthquake (Mw = 7.6) in Taiwan. The Chi-Chi earthquake provides an unique opportunity to investigate the spatial and time-dependent poroelastic response of near-field rocks and sediments because there exist extensive observational data of water-level changes and crustal deformation. The integrated model allows us to explore whether changes in Coulomb stress can adequately explain hydrologic anomalies observed in areas such as Taiwan’s western foothills and the Choshui River alluvial plain. To calculate coseismic strain, we use the carefully calibrated finite fault-rupture model of Ma et al. (2005) and the deformation modeling code Coulomb 3.1 (Toda et al., 2007
An Extension of Implicit Monte Carlo Diffusion: Multigroup and The Difference Formulation
Cleveland, M A; Gentile, N; Palmer, T S
2010-04-19
Implicit Monte Carlo (IMC) and Implicit Monte Carlo Diffusion (IMD) are approaches to the numerical solution of the equations of radiative transfer. IMD was previously derived and numerically tested on grey, or frequency-integrated problems. In this research, we extend Implicit Monte Carlo Diffusion (IMD) to account for frequency dependence, and we implement the difference formulation as a source manipulation variance reduction technique. We derive the relevant probability distributions and present the frequency dependent IMD algorithm, with and without the difference formulation. The IMD code with and without the difference formulation was tested using both grey and frequency dependent benchmark problems. The Su and Olson semi-analytic Marshak wave benchmark was used to demonstrate the validity of the code for grey problems. The Su and Olson semi-analytic picket fence benchmark was used for the frequency dependent problems. The frequency dependent IMD algorithm reproduces the results of both Su and Olson benchmark problems. Frequency group refinement studies indicate that the computational cost of refining the group structure is likely less than that of group refinement in deterministic solutions of the radiation diffusion methods. Our results show that applying the difference formulation to the IMD algorithm can result in an overall increase in the figure of merit for frequency dependent problems. However, the creation of negatively weighted particles from the difference formulation can cause significant numerical instabilities in regions of the problem with sharp spatial gradients in the solution. An adaptive implementation of the difference formulation may be necessary to focus its use in regions that are at or near thermal equilibrium.
An extension of implicit Monte Carlo diffusion: Multigroup and the difference formulation
NASA Astrophysics Data System (ADS)
Cleveland, Mathew A.; Gentile, Nick A.; Palmer, Todd S.
2010-08-01
Implicit Monte Carlo (IMC) and Implicit Monte Carlo Diffusion (IMD) are approaches to the numerical solution of the equations of radiative transfer. IMD was previously derived and numerically tested on grey, or frequency-integrated problems [1]. In this research, we extend Implicit Monte Carlo Diffusion (IMD) to account for frequency dependence, and we implement the difference formulation[2] as a source manipulation variance reduction technique. We derive the relevant probability distributions and present the frequency dependent IMD algorithm, with and without the difference formulation. The IMD code with and without the difference formulation was tested using both grey and frequency dependent benchmark problems. The Su and Olson semi-analytic Marshak wave benchmark was used to demonstrate the validity of the code for grey problems [3]. The Su and Olson semi-analytic picket fence benchmark was used for the frequency dependent problems [4]. The frequency dependent IMD algorithm reproduces the results of both Su and Olson benchmark problems. Frequency group refinement studies indicate that the computational cost of refining the group structure is likely less than that of group refinement in deterministic solutions of the radiation diffusion methods. Our results show that applying the difference formulation to the IMD algorithm can result in an overall increase in the figure of merit for frequency dependent problems. However, the creation of negatively weighted particles from the difference formulation can cause significant numerical instabilities in regions of the problem with sharp spatial gradients in the solution. An adaptive implementation of the difference formulation may be necessary to focus its use in regions that are at or near thermal equilibrium.
Mapping the holes: 3D ISM maps and diffuse X-ray background
NASA Astrophysics Data System (ADS)
Lallement, R.; Vergely, J.-L.; Puspitarini, L.; Snowden, S.; Galeazzi, M.; Koutroumpa, D.
3D maps of Galactic interstellar dust and gas reveal empty regions, including cavities carved by stellar winds and supernovae. Such cavities are often filled with hot gas and are sources of soft X-ray background emission. We discuss the combined analysis of the diffuse soft (0.25 keV) X-ray background and the 3D distribution of nearby (<1 kpc) dust, including studies of shadows cast by nearby clouds in the background. This analysis benefits from recent progress in the estimate of the foreground X-ray emission from the heliosphere. New and past X-ray data are found to be consistent with the maps if the ≃ 100-150 pc wide Local Bubble surrounding the Sun is filled with 106K gas with a pressure 2nT ≃ 10,000 K cm-3. On the other hand, the giant cavity found in the 3rd Galactic quadrant has a weaker volume emission than the LB and is very likely filled to a large extent with warm ionized gas. Its geometry suggests a link with the tilted Gould belt, and a potential mechanism for the formation of the whole structure has been recently proposed. According to it, the local inclination of gas and stars, the velocity pattern and enhanced star formation could have been initiated 60-70 Myr ago when a massive globular cluster crossed the Galactic Plane in the vicinity of the Sun. The destabilization of stellar orbits around the Sun may have generated enhanced asteroid falls of the Cretaceous-Tertiary (KT) extinction events. Additionally, a short gamma ray burst may have occurred in the cluster during the crossing, producing intense ionization and subsequent shock waves leading to the star formations seen today in the form of the giant ionized region and OB associations at its periphery. Gaia measurements of nearby stars and clusters should help shedding light on the local history.
Parareal in time 3D numerical solver for the LWR Benchmark neutron diffusion transient model
NASA Astrophysics Data System (ADS)
Baudron, Anne-Marie; Lautard, Jean-Jacques; Maday, Yvon; Riahi, Mohamed Kamel; Salomon, Julien
2014-12-01
In this paper we present a time-parallel algorithm for the 3D neutrons calculation of a transient model in a nuclear reactor core. The neutrons calculation consists in numerically solving the time dependent diffusion approximation equation, which is a simplified transport equation. The numerical resolution is done with finite elements method based on a tetrahedral meshing of the computational domain, representing the reactor core, and time discretization is achieved using a θ-scheme. The transient model presents moving control rods during the time of the reaction. Therefore, cross-sections (piecewise constants) are taken into account by interpolations with respect to the velocity of the control rods. The parallelism across the time is achieved by an adequate use of the parareal in time algorithm to the handled problem. This parallel method is a predictor corrector scheme that iteratively combines the use of two kinds of numerical propagators, one coarse and one fine. Our method is made efficient by means of a coarse solver defined with large time step and fixed position control rods model, while the fine propagator is assumed to be a high order numerical approximation of the full model. The parallel implementation of our method provides a good scalability of the algorithm. Numerical results show the efficiency of the parareal method on large light water reactor transient model corresponding to the Langenbuch-Maurer-Werner benchmark.
Dynamic implicit 3D adaptive mesh refinement for non-equilibrium radiation diffusion
B. Philip; Z. Wang; M.A. Berrill; M. Birke; M. Pernice
2014-04-01
The time dependent non-equilibrium radiation diffusion equations are important for solving the transport of energy through radiation in optically thick regimes and find applications in several fields including astrophysics and inertial confinement fusion. The associated initial boundary value problems that are encountered often exhibit a wide range of scales in space and time and are extremely challenging to solve. To efficiently and accurately simulate these systems we describe our research on combining techniques that will also find use more broadly for long term time integration of nonlinear multi-physics systems: implicit time integration for efficient long term time integration of stiff multi-physics systems, local control theory based step size control to minimize the required global number of time steps while controlling accuracy, dynamic 3D adaptive mesh refinement (AMR) to minimize memory and computational costs, Jacobian Free Newton–Krylov methods on AMR grids for efficient nonlinear solution, and optimal multilevel preconditioner components that provide level independent solver convergence.
Parareal in time 3D numerical solver for the LWR Benchmark neutron diffusion transient model
Baudron, Anne-Marie; Riahi, Mohamed Kamel; Salomon, Julien
2014-12-15
In this paper we present a time-parallel algorithm for the 3D neutrons calculation of a transient model in a nuclear reactor core. The neutrons calculation consists in numerically solving the time dependent diffusion approximation equation, which is a simplified transport equation. The numerical resolution is done with finite elements method based on a tetrahedral meshing of the computational domain, representing the reactor core, and time discretization is achieved using a θ-scheme. The transient model presents moving control rods during the time of the reaction. Therefore, cross-sections (piecewise constants) are taken into account by interpolations with respect to the velocity of the control rods. The parallelism across the time is achieved by an adequate use of the parareal in time algorithm to the handled problem. This parallel method is a predictor corrector scheme that iteratively combines the use of two kinds of numerical propagators, one coarse and one fine. Our method is made efficient by means of a coarse solver defined with large time step and fixed position control rods model, while the fine propagator is assumed to be a high order numerical approximation of the full model. The parallel implementation of our method provides a good scalability of the algorithm. Numerical results show the efficiency of the parareal method on large light water reactor transient model corresponding to the Langenbuch–Maurer–Werner benchmark.
Gerhard Strydom; Cristian Rabiti; Andrea Alfonsi
2012-10-01
PHISICS is a neutronics code system currently under development at the Idaho National Laboratory (INL). Its goal is to provide state of the art simulation capability to reactor designers. The different modules for PHISICS currently under development are a nodal and semi-structured transport core solver (INSTANT), a depletion module (MRTAU) and a cross section interpolation (MIXER) module. The INSTANT module is the most developed of the mentioned above. Basic functionalities are ready to use, but the code is still in continuous development to extend its capabilities. This paper reports on the effort of coupling the nodal kinetics code package PHISICS (INSTANT/MRTAU/MIXER) to the thermal hydraulics system code RELAP5-3D, to enable full core and system modeling. This will enable the possibility to model coupled (thermal-hydraulics and neutronics) problems with more options for 3D neutron kinetics, compared to the existing diffusion theory neutron kinetics module in RELAP5-3D (NESTLE). In the second part of the paper, an overview of the OECD/NEA MHTGR-350 MW benchmark is given. This benchmark has been approved by the OECD, and is based on the General Atomics 350 MW Modular High Temperature Gas Reactor (MHTGR) design. The benchmark includes coupled neutronics thermal hydraulics exercises that require more capabilities than RELAP5-3D with NESTLE offers. Therefore, the MHTGR benchmark makes extensive use of the new PHISICS/RELAP5-3D coupling capabilities. The paper presents the preliminary results of the three steady state exercises specified in Phase I of the benchmark using PHISICS/RELAP5-3D.
Doblhoff-Dier, Katharina; Meyer, Jörg; Hoggan, Philip E; Kroes, Geert-Jan; Wagner, Lucas K
2016-06-14
Transition metals and transition metal compounds are important to catalysis, photochemistry, and many superconducting systems. We study the performance of diffusion Monte Carlo (DMC) applied to transition metal containing dimers (TMCDs) using single-determinant Slater-Jastrow trial wavefunctions and investigate the possible influence of the locality and pseudopotential errors. We find that the locality approximation can introduce nonsystematic errors of up to several tens of kilocalories per mole in the absolute energy of Cu and CuH if Ar or Mg core pseudopotentials (PPs) are used for the 3d transition metal atoms. Even for energy differences such as binding energies, errors due to the locality approximation can be problematic if chemical accuracy is sought. The use of the Ne core PPs developed by Burkatzki et al. (J. Chem. Phys. 2008, 129, 164115), the use of linear energy minimization rather than unreweighted variance minimization for the optimization of the Jastrow function, and the use of large Jastrow parametrizations reduce the locality errors. In the second section of this article, we study the general performance of DMC for 3d TMCDs using a database of binding energies of 20 TMCDs, for which comparatively accurate experimental data is available. Comparing our DMC results to these data for our results that compare best with experiment, we find a mean unsigned error (MUE) of 4.5 kcal/mol. This compares well with the achievable accuracy in CCSDT(2)Q (MUE = 4.6 kcal/mol) and the best all-electron DFT results (MUE = 4.5 kcal/mol) for the same set of systems (Truhlar et al. J. Chem. Theory Comput. 2015, 11, 2036-2052). The mean errors in DMC depend less on the exchange-correlation functionals used to generate the trial wavefunction than the corresponding mean errors in the underlying DFT calculations. Furthermore, the QMC results obtained for each molecule individually vary less with the functionals used. These observations are relevant for systems such as
Study of a non-diffusing radiochromic gel dosimeter for 3D radiation dose imaging
NASA Astrophysics Data System (ADS)
Marsden, Craig Michael
2000-12-01
This thesis investigates the potential of a new radiation gel dosimeter, based on nitro-blue tetrazolium (NBTZ) suspended in a gelatin mold. Unlike all Fricke based gel dosimeters this dosimeter does not suffer from diffusive loss of image stability. Images are obtained by an optical tomography method. Nitro blue tetrazolium is a common biological indicator that when irradiated in an aqueous medium undergoes reduction to a highly colored formazan, which has an absorbance maximum at 525nm. Tetrazolium is water soluble while the formazan product is insoluble. The formazan product sticks to the gelatin matrix and the dose image is maintained for three months. Methods to maximize the sensitivity of the system were evaluated. It was found that a chemical detergent, Triton X-100, in combination with sodium formate, increased the dosimeter sensitivity significantly. An initial G-value of formazan production for a dosimeter composed of 1mM NBTZ, gelatin, and water was on the order of 0.2. The addition of Triton and formate produced a G-value in excess of 5.0. The effects of NBTZ, triton, formate, and gel concentration were all investigated. All the gels provided linear dose vs. absorbance plots for doses from 0 to >100 Gy. It was determined that gel concentration had minimal if any effect on sensitivity. Sensitivity increased slightly with increasing NBTZ concentration. Triton and formate individually and together provided moderate to large increases in dosimeter sensitivity. The dosimeter described in this work can provide stable 3D radiation dose images for all modalities of radiation therapy equipment. Methods to increase sensitivity are developed and discussed.
A remark on the Beale-Kato-Majda criterion for the 3D MHD equations with zero magnetic diffusivity
NASA Astrophysics Data System (ADS)
Gala, Sadek; Ragusa, Maria Alessandra
2016-06-01
In this work, we show that a smooth solution of the 3D MHD equations with zero magnetic diffusivity in the whole space ℝ3 breaks down if and only if a certain norm of the magnetic field blows up at the same time.
Fogtmann, Mads; Seshamani, Sharmishtaa; Kroenke, Christopher; Cheng, Xi; Chapman, Teresa; Wilm, Jakob; Rousseau, François
2014-01-01
This paper presents an approach to 3-D diffusion tensor image (DTI) reconstruction from multi-slice diffusion weighted (DW) magnetic resonance imaging acquisitions of the moving fetal brain. Motion scatters the slice measurements in the spatial and spherical diffusion domain with respect to the underlying anatomy. Previous image registration techniques have been described to estimate the between slice fetal head motion, allowing the reconstruction of 3-D a diffusion estimate on a regular grid using interpolation. We propose Approach to Unified Diffusion Sensitive Slice Alignment and Reconstruction (AUDiSSAR) that explicitly formulates a process for diffusion direction sensitive DW-slice-to-DTI-volume alignment. This also incorporates image resolution modeling to iteratively deconvolve the effects of the imaging point spread function using the multiple views provided by thick slices acquired in different anatomical planes. The algorithm is implemented using a multi-resolution iterative scheme and multiple real and synthetic data are used to evaluate the performance of the technique. An accuracy experiment using synthetically created motion data of an adult head and a experiment using synthetic motion added to sedated fetal monkey dataset show a significant improvement in motion-trajectory estimation compared to a state-of-the-art approaches. The performance of the method is then evaluated on challenging but clinically typical in utero fetal scans of four different human cases, showing improved rendition of cortical anatomy and extraction of white matter tracts. While the experimental work focuses on DTI reconstruction (second-order tensor model), the proposed reconstruction framework can employ any 5-D diffusion volume model that can be represented by the spatial parameterizations of an orientation distribution function. PMID:24108711
A fully implicit method for 3D quasi-steady state magnetic advection-diffusion.
Siefert, Christopher; Robinson, Allen Conrad
2009-09-01
We describe the implementation of a prototype fully implicit method for solving three-dimensional quasi-steady state magnetic advection-diffusion problems. This method allows us to solve the magnetic advection diffusion equations in an Eulerian frame with a fixed, user-prescribed velocity field. We have verified the correctness of method and implementation on two standard verification problems, the Solberg-White magnetic shear problem and the Perry-Jones-White rotating cylinder problem.
NASA Technical Reports Server (NTRS)
Fletcher, Michael J.; Won, Mark J.; Cosentino, Gary B.; Te, Alexander
1993-01-01
Subsonic inlet ducts for advanced, high-performance aircraft are evolving towards complex three-dimensional shapes for reasons of overall integration and weight. These factors lead to diffuser geometries that may sacrifice inlet performance, unless careful attention to design details and boundary layer management techniques are employed. The ability of viscous computational fluid dynamic (CFD) analysis of such geometries to aid the aircraft configurator in this complex design problem is herein examined. The RANS-3D Reynolds-Averaged Navier-Stokes solver is applied to model the complex flowfield occurring in a representative diffuser geometry and the solutions are compared to experimental results from a static test of the inlet duct. The computational results are shown to compare very favorably with experimental results over a range of mass flow rates, including those involving large amounts of separation in the diffuser. In addition, a novel grid topology is presented, and two turbulence models are evaluated in this study as part of the RANS-3D code.
NASA Astrophysics Data System (ADS)
Xiu, Ziyang; Deng, Zongquan; Wang, Xiaofeng; Wu, Gaohui
A dense and uniform Sip/LG5 composite were fabricated by squeeze casting technology, and high temperature diffusion treatment was adapted to the composite. Microstructure observation indicated that Si transformed from irregular particles to 3D-structure. Fine dispersive precipitates Si were also observed on Si-Al interface and within Al matrix, smoothing and improving the interface. Based on the microstructure observation, three transformation stages were designated: melting, dissolution and precipitation, solidification. Thermodynamics and kinetics of the transformation can be explained by Gibbs-Thomson effect.
Nodal Diffusion & Transport Theory
Energy Science and Technology Software Center (ESTSC)
1992-02-19
DIF3D solves multigroup diffusion theory eigenvalue, adjoint, fixed source, and criticality (concentration, buckling, and dimension search) problems in 1, 2, and 3-space dimensions for orthogonal (rectangular or cylindrical), triangular, and hexagonal geometries. Anisotropic diffusion theory coefficients are permitted. Flux and power density maps by mesh cell and regionwise balance integrals are provided. Although primarily designed for fast reactor problems, upscattering and internal black boundary conditions are also treated.
NASA Astrophysics Data System (ADS)
Peng, Hongyun; Wen, Huanyao; Zhu, Changjiang
2014-12-01
In this paper, we study the existence of global classical solutions and the vanishing diffusion limit of a 3D conservation laws derived from the well-known Keller-Segel model. First, we establish the global well-posedness of classical solutions to the Cauchy problem for the model with smooth initial data which is of small L 2 norm, together with some a priori estimates uniform for t and . Then, we investigate the zero diffusion limit and get the global well-posedness of classical solutions to the Cauchy problem for the non-diffusive model. Finally, we derive the convergence rate of the model toward the non-diffusive model. It is shown that the convergence rate in L ∞ norm is of the order . It should be noted that the initial data are small in L 2-norm but can be of large oscillations with constant state at far field. As a byproduct, we improve the corresponding result on the well-posedness of the non-diffusive model which requires small oscillations.
3-D Deep Penetration Neutron Imaging of Thick Absorgin and Diffusive Objects Using Transport Theory
Ragusa, Jean; Bangerth, Wolfgang
2011-08-01
here explores the inverse problem of optical tomography applied to heterogeneous domains. The neutral particle transport equation was used as the forward model for how neutral particles stream through and interact within these heterogeneous domains. A constrained optimization technique that uses Newtons method served as the basis of the inverse problem. Optical tomography aims at reconstructing the material properties using (a) illuminating sources and (b) detector readings. However, accurate simulations for radiation transport require that the particle (gamma and/or neutron) energy be appropriate discretize in the multigroup approximation. This, in turns, yields optical tomography problems where the number of unknowns grows (1) about quadratically with respect to the number of energy groups, G, (notably to reconstruct the scattering matrix) and (2) linearly with respect to the number of unknown material regions. As pointed out, a promising approach could rely on algorithms to appropriately select a material type per material zone rather than G2 values. This approach, though promising, still requires further investigation: (a) when switching from cross-section values unknowns to material type indices (discrete integer unknowns), integer programming techniques are needed since derivative information is no longer available; and (b) the issue of selecting the initial material zoning remains. The work reported here proposes an approach to solve the latter item, whereby a material zoning is proposed using one-group or few-groups transport approximations. The capabilities and limitations of the presented method were explored; they are briefly summarized next and later described in fuller details in the Appendices. The major factors that influenced the ability of the optimization method to reconstruct the cross sections of these domains included the locations of the sources used to illuminate the domains, the number of separate experiments used in the reconstruction, the
NASA Astrophysics Data System (ADS)
Hu, Y.; Ji, Y.; Egbert, G. D.
2015-12-01
The fictitious time domain method (FTD), based on the correspondence principle for wave and diffusion fields, has been developed and used over the past few years primarily for marine electromagnetic (EM) modeling. Here we present results of our efforts to apply the FTD approach to land and airborne TEM problems which can reduce the computer time several orders of magnitude and preserve high accuracy. In contrast to the marine case, where sources are in the conductive sea water, we must model the EM fields in the air; to allow for topography air layers must be explicitly included in the computational domain. Furthermore, because sources for most TEM applications generally must be modeled as finite loops, it is useful to solve directly for the impulse response appropriate to the problem geometry, instead of the point-source Green functions typically used for marine problems. Our approach can be summarized as follows: (1) The EM diffusion equation is transformed to a fictitious wave equation. (2) The FTD wave equation is solved with an explicit finite difference time-stepping scheme, with CPML (Convolutional PML) boundary conditions for the whole computational domain including the air and earth , with FTD domain source corresponding to the actual transmitter geometry. Resistivity of the air layers is kept as low as possible, to compromise between efficiency (longer fictitious time step) and accuracy. We have generally found a host/air resistivity contrast of 10-3 is sufficient. (3)A "Modified" Fourier Transform (MFT) allow us recover system's impulse response from the fictitious time domain to the diffusion (frequency) domain. (4) The result is multiplied by the Fourier transformation （FT） of the real source current avoiding time consuming convolutions in the time domain. (5) The inverse FT is employed to get the final full waveform and full time response of the system in the time domain. In general, this method can be used to efficiently solve most time-domain EM
FAST SIMULATION OF SOLID TUMORS THERMAL ABLATION TREATMENTS WITH A 3D REACTION DIFFUSION MODEL *
BERTACCINI, DANIELE; CALVETTI, DANIELA
2007-01-01
An efficient computational method for near real-time simulation of thermal ablation of tumors via radio frequencies is proposed. Model simulations of the temperature field in a 3D portion of tissue containing the tumoral mass for different patterns of source heating can be used to design the ablation procedure. The availability of a very efficient computational scheme makes it possible update the predicted outcome of the procedure in real time. In the algorithms proposed here a discretization in space of the governing equations is followed by an adaptive time integration based on implicit multistep formulas. A modification of the ode15s MATLAB function which uses Krylov space iterative methods for the solution of for the linear systems arising at each integration step makes it possible to perform the simulations on standard desktop for much finer grids than using the built-in ode15s. The proposed algorithm can be applied to a wide class of nonlinear parabolic differential equations. PMID:17173888
Fast simulation of solid tumors thermal ablation treatments with a 3D reaction diffusion model.
Bertaccini, Daniele; Calvetti, Daniela
2007-08-01
An efficient computational method for near real-time simulation of thermal ablation of tumors via radio frequencies is proposed. Model simulations of the temperature field in a 3D portion of tissue containing the tumoral mass for different patterns of source heating can be used to design the ablation procedure. The availability of a very efficient computational scheme makes it possible to update the predicted outcome of the procedure in real time. In the algorithms proposed here a discretization in space of the governing equations is followed by an adaptive time integration based on implicit multistep formulas. A modification of the ode15s MATLAB function which uses Krylov space iterative methods for the solution of the linear systems arising at each integration step makes it possible to perform the simulations on standard desktop for much finer grids than using the built-in ode15s. The proposed algorithm can be applied to a wide class of nonlinear parabolic differential equations. PMID:17173888
1,2 or 3-D Few-Group Diffusion Depletion
Energy Science and Technology Software Center (ESTSC)
1982-12-01
The PDQ series of programs is designed to solve the neutron diffusion-depletion problem in one, two or three dimensions. The three-dimensional spatial calculation may be either explicit or discontinuous trial function synthesis. Up to five lethargy groups are permitted. Adjoint, fixed source, one iteration, additive fixed source, eigenvalue, and boundary value calculations may be performed. The programs utilize the HARMONY system for time-dependent representation of cross section variation and generalized depletion chain solutions. Geometries availablemore » include rectangular, cylindrical, spherical, and hexagonal. All allow variable mesh in all dimensions. Various control searches as well as temperature and xenon feedback are provided.« less
MO-G-BRF-07: Anomalously Fast Diffusion of Carbon Nanotubes Carriers in 3D Tissue Model
Wang, Y; Bahng, J; Kotov, N
2014-06-15
Purpose: We aim to investigate and understand diffusion process of carbon nanotubes (CNTs) and other nanoscale particles in tissue and organs. Methods: In this research, we utilized a 3D model tissue of hepatocellular carcinoma (HCC)cultured in inverted colloidal crystal (ICC) scaffolds to compare the diffusivity of CNTs with small molecules such as Rhodamine and FITC in vitro, and further investigated the transportation of CNTs with and without targeting ligand, TGFβ1. The real-time permeation profiles of CNTs in HCC tissue model with high temporal and spatial resolution was demonstrated by using standard confocal microscopy. Quantitative analysis of the diffusion process in 3D was carried out using luminescence intensity in a series of Z-stack images obtained for different time points of the diffusion process after initial addition of CNTs or small molecules to the cell culture and the image data was analyzed by software ImageJ and Mathematica. Results: CNTs display diffusion rate in model tissues substantially faster than small molecules of the similar charge such as FITC, and the diffusion rate of CNTs are significantly enhanced with targeting ligand, TGFβ1. Conclusion: In terms of the advantages of in-vitro model, we were able to have access to measuring the rate of CNT penetration at designed conditions with variable parameters. And the findings by using this model, changed our understanding about advantages of CNTs as nanoscale drug carriers and provides design principles for making new drug carriers for both treatment and diagnostics. Additionally the fast diffusion opens the discussion of the best possible drug carriers to reach deep parts of cancerous tissues, which is often a prerequisite for successful cancer treatment. This work was supported by the Center for Photonic and Multiscale Nanomaterials funded by National Science Foundation Materials Research Science and Engineering Center program DMR 1120923. The work was also partially supported by NSF
Time-resolved diffusion tomographic 2D and 3D imaging in highly scattering turbid media
NASA Technical Reports Server (NTRS)
Alfano, Robert R. (Inventor); Cai, Wei (Inventor); Gayen, Swapan K. (Inventor)
2000-01-01
A method for imaging objects in highly scattering turbid media. According to one embodiment of the invention, the method involves using a plurality of intersecting source/detectors sets and time-resolving equipment to generate a plurality of time-resolved intensity curves for the diffusive component of light emergent from the medium. For each of the curves, the intensities at a plurality of times are then inputted into the following inverse reconstruction algorithm to form an image of the medium: wherein W is a matrix relating output at source and detector positions r.sub.s and r.sub.d, at time t, to position r, .LAMBDA. is a regularization matrix, chosen for convenience to be diagonal, but selected in a way related to the ratio of the noise,
Time-resolved diffusion tomographic 2D and 3D imaging in highly scattering turbid media
NASA Technical Reports Server (NTRS)
Alfano, Robert R. (Inventor); Cai, Wei (Inventor); Liu, Feng (Inventor); Lax, Melvin (Inventor); Das, Bidyut B. (Inventor)
1999-01-01
A method for imaging objects in highly scattering turbid media. According to one embodiment of the invention, the method involves using a plurality of intersecting source/detectors sets and time-resolving equipment to generate a plurality of time-resolved intensity curves for the diffusive component of light emergent from the medium. For each of the curves, the intensities at a plurality of times are then inputted into the following inverse reconstruction algorithm to form an image of the medium: ##EQU1## wherein W is a matrix relating output at source and detector positions r.sub.s and r.sub.d, at time t, to position r, .LAMBDA. is a regularization matrix, chosen for convenience to be diagonal, but selected in a way related to the ratio of the noise,
Diffusion of 3D-transition elements in Ti-54at%Al
Lee, C.G.; Iijima, Y.; Kim, S.E.; Lee, Y.T.; Kim, H.M.
1995-12-31
The diffusion coefficients of iron, cobalt, manganese and chromium in Ti-54at%Al alloy were measured over the temperature range from 1,000 to 1,540 K. This study used the radioactive tracers {sup 59}Fe, {sup 57}Co, {sup 54}Mn and {sup 51}Cr and employed the serial the radio-frequency sputter-microsectioning method to measure the penetration profiles of the radioisotope into the specimen. The temperature dependence of the diffusion coefficients was analyzed to be expressed by the following Arrhenius equations: D{sub Fe/TiAl} = (1.3 {times} 10{sup {minus}3}) exp ({minus}363 kJ mol{sup {minus}1}/RT) m{sup 2}/s, D{sub Co/TiAl} = (1.1 {times} 10{sup {minus}3}) exp ({minus}318 kJ mol{sup {minus}1}/RT) m{sup 2}/s, D{sub Mn/TiAl} = (1.2 {times} 10{sup {minus}3}) exp ({minus}326 kJ mol{sup {minus}1}/RT) m{sup 2}/s, D{sub Cr/TiAl} = (4.4 {times} 10{sup {minus}3}) exp ({minus}350 kJ mol{sup {minus}1}/RT) m{sup 2}/s.
Schütze, Friedrich; Röhrig, Florian; Vorlová, Sandra; Gätzner, Sabine; Kuhn, Anja; Ergün, Süleyman; Henke, Erik
2015-01-01
Tumors are characterized by a rigid, highly cross-linked extracellular matrix (ECM), which impedes homogeneous drug distribution and potentially protects malignant cells from exposure to therapeutics. Lysyl oxidases are major contributors to tissue stiffness and the elevated expression of these enzymes observed in most cancers might influence drug distribution and efficacy. We examined the effect of lysyl oxidases on drug distribution and efficacy in 3D in vitro assay systems. In our experiments elevated lysyl oxidase activity was responsible for reduced drug diffusion under hypoxic conditions and consequently impaired cytotoxicity of various chemotherapeutics. This effect was only observed in 3D settings but not in 2D-cell culture, confirming that lysyl oxidases affect drug efficacy by modification of the ECM and do not confer a direct desensitizing effect. Both drug diffusion and efficacy were strongly enhanced by inhibition of lysyl oxidases. The results from the in vitro experiments correlated with tumor drug distribution in vivo, and predicted response to therapeutics in murine tumor models. Our results demonstrate that lysyl oxidase activity modulates the physical barrier function of ECM for small molecule drugs influencing their therapeutic efficacy. Targeting this process has the potential to significantly enhance therapeutic efficacy in the treatment of malignant diseases. PMID:26620400
Van, Anh T; Aksoy, Murat; Holdsworth, Samantha J; Kopeinigg, Daniel; Vos, Sjoerd B; Bammer, Roland
2014-01-01
Purpose To propose a method for mitigating slab boundary artifacts in 3D multislab diffusion imaging with no or minimal increases in scan time. Methods The multislab acquisition was treated as parallel imaging acquisition where the slab profiles acted as the traditional receiver sensitivity profiles. All the slabs were then reconstructed simultaneously along the slab direction using Cartesian-based sensitivity encoding (SENSE) reconstruction. The slab profile estimation was performed using either a Bloch simulation or a calibration scan. Results Both phantom and in vivo results showed negligible slab boundary artifacts after reconstruction using the proposed method. The performance of the proposed method is comparable to the state-of-the-art slab combination method without the scan time penalty that depends on the number of acquired volumes. The obtained g-factor map of the SENSE reconstruction problem showed a maximum g-factor of 1.7 in the region of interest. Conclusion We proposed a novel method for mitigating slab boundary artifacts in 3D diffusion imaging by treating the multislab acquisition as a parallel imaging acquisition and reconstructing all slabs simultaneously using Cartesian SENSE. Unlike existing methods, the scan time increase, if any, does not scale with the number of image volumes acquired. PMID:24691843
Lawrence, R.D.
1983-03-01
A nodal method is developed for the solution of the neutron-diffusion equation in two- and three-dimensional hexagonal geometries. The nodal scheme has been incorporated as an option in the finite-difference diffusion-theory code DIF3D, and is intended for use in the analysis of current LMFBR designs. The nodal equations are derived using higher-order polynomial approximations to the spatial dependence of the flux within the hexagonal-z node. The final equations, which are cast in the form of inhomogeneous response-matrix equations for each energy group, involved spatial moments of the node-interior flux distribution plus surface-averaged partial currents across the faces of the node. These equations are solved using a conventional fission-source iteration accelerated by coarse-mesh rebalance and asymptotic source extrapolation. This report describes the mathematical development and numerical solution of the nodal equations, as well as the use of the nodal option and details concerning its programming structure. This latter information is intended to supplement the information provided in the separate documentation of the DIF3D code.
Cai, Yu; McMurray, Matthew S.; Oguz, Ipek; Yuan, Hong; Styner, Martin A.; Lin, Weili; Johns, Josephine M.; An, Hongyu
2011-01-01
High resolution diffusion tensor imaging (DTI) can provide important information on brain development, yet it is challenging in live neonatal rats due to the small size of neonatal brain and motion-sensitive nature of DTI. Imaging in live neonatal rats has clear advantages over fixed brain scans, as longitudinal and functional studies would be feasible to understand neuro-developmental abnormalities. In this study, we developed imaging strategies that can be used to obtain high resolution 3D DTI images in live neonatal rats at postnatal day 5 (PND5) and PND14, using only 3 h of imaging acquisition time. An optimized 3D DTI pulse sequence and appropriate animal setup to minimize physiological motion artifacts are the keys to successful high resolution 3D DTI imaging. Thus, a 3D rapid acquisition relaxation enhancement DTI sequence with twin navigator echoes was implemented to accelerate imaging acquisition time and minimize motion artifacts. It has been suggested that neonatal mammals possess a unique ability to tolerate mild-to-moderate hypothermia and hypoxia without long term impact. Thus, we additionally utilized this ability to minimize motion artifacts in magnetic resonance images by carefully suppressing the respiratory rate to around 15/min for PND5 and 30/min for PND14 using mild-to-moderate hypothermia. These imaging strategies have been successfully implemented to study how the effect of cocaine exposure in dams might affect brain development in their rat pups. Image quality resulting from this in vivo DTI study was comparable to ex vivo scans. fractional anisotropy values were also similar between the live and fixed brain scans. The capability of acquiring high quality in vivo DTI imaging offers a valuable opportunity to study many neurological disorders in brain development in an authentic living environment. PMID:22013426
Lee, Shiu-Hang; Kamae, Tuneyoshi; Ellison, Donald C.
2008-07-02
We present a 3-dimensional model of supernova remnants (SNRs) where the hydrodynamical evolution of the remnant is modeled consistently with nonlinear diffusive shock acceleration occurring at the outer blast wave. The model includes particle escape and diffusion outside of the forward shock, and particle interactions with arbitrary distributions of external ambient material, such as molecular clouds. We include synchrotron emission and cooling, bremsstrahlung radiation, neutral pion production, inverse-Compton (IC), and Coulomb energy-loss. Boardband spectra have been calculated for typical parameters including dense regions of gas external to a 1000 year old SNR. In this paper, we describe the details of our model but do not attempt a detailed fit to any specific remnant. We also do not include magnetic field amplification (MFA), even though this effect may be important in some young remnants. In this first presentation of the model we don't attempt a detailed fit to any specific remnant. Our aim is to develop a flexible platform, which can be generalized to include effects such as MFA, and which can be easily adapted to various SNR environments, including Type Ia SNRs, which explode in a constant density medium, and Type II SNRs, which explode in a pre-supernova wind. When applied to a specific SNR, our model will predict cosmic-ray spectra and multi-wavelength morphology in projected images for instruments with varying spatial and spectral resolutions. We show examples of these spectra and images and emphasize the importance of measurements in the hard X-ray, GeV, and TeV gamma-ray bands for investigating key ingredients in the acceleration mechanism, and for deducing whether or not TeV emission is produced by IC from electrons or pion-decay from protons.
NASA Astrophysics Data System (ADS)
Bhowmik, Tanmoy; Liu, Hanli; Ye, Zhou; Oraintara, Soontorn
2016-03-01
Diffuse optical tomography (DOT) is a relatively low cost and portable imaging modality for reconstruction of optical properties in a highly scattering medium, such as human tissue. The inverse problem in DOT is highly ill-posed, making reconstruction of high-quality image a critical challenge. Because of the nature of sparsity in DOT, sparsity regularization has been utilized to achieve high-quality DOT reconstruction. However, conventional approaches using sparse optimization are computationally expensive and have no selection criteria to optimize the regularization parameter. In this paper, a novel algorithm, Dimensionality Reduction based Optimization for DOT (DRO-DOT), is proposed. It reduces the dimensionality of the inverse DOT problem by reducing the number of unknowns in two steps and thereby makes the overall process fast. First, it constructs a low resolution voxel basis based on the sensing-matrix properties to find an image support. Second, it reconstructs the sparse image inside this support. To compensate for the reduced sensitivity with increasing depth, depth compensation is incorporated in DRO-DOT. An efficient method to optimally select the regularization parameter is proposed for obtaining a high-quality DOT image. DRO-DOT is also able to reconstruct high-resolution images even with a limited number of optodes in a spatially limited imaging set-up.
Bhowmik, Tanmoy; Liu, Hanli; Ye, Zhou; Oraintara, Soontorn
2016-01-01
Diffuse optical tomography (DOT) is a relatively low cost and portable imaging modality for reconstruction of optical properties in a highly scattering medium, such as human tissue. The inverse problem in DOT is highly ill-posed, making reconstruction of high-quality image a critical challenge. Because of the nature of sparsity in DOT, sparsity regularization has been utilized to achieve high-quality DOT reconstruction. However, conventional approaches using sparse optimization are computationally expensive and have no selection criteria to optimize the regularization parameter. In this paper, a novel algorithm, Dimensionality Reduction based Optimization for DOT (DRO-DOT), is proposed. It reduces the dimensionality of the inverse DOT problem by reducing the number of unknowns in two steps and thereby makes the overall process fast. First, it constructs a low resolution voxel basis based on the sensing-matrix properties to find an image support. Second, it reconstructs the sparse image inside this support. To compensate for the reduced sensitivity with increasing depth, depth compensation is incorporated in DRO-DOT. An efficient method to optimally select the regularization parameter is proposed for obtaining a high-quality DOT image. DRO-DOT is also able to reconstruct high-resolution images even with a limited number of optodes in a spatially limited imaging set-up. PMID:26940661
Bhowmik, Tanmoy; Liu, Hanli; Ye, Zhou; Oraintara, Soontorn
2016-01-01
Diffuse optical tomography (DOT) is a relatively low cost and portable imaging modality for reconstruction of optical properties in a highly scattering medium, such as human tissue. The inverse problem in DOT is highly ill-posed, making reconstruction of high-quality image a critical challenge. Because of the nature of sparsity in DOT, sparsity regularization has been utilized to achieve high-quality DOT reconstruction. However, conventional approaches using sparse optimization are computationally expensive and have no selection criteria to optimize the regularization parameter. In this paper, a novel algorithm, Dimensionality Reduction based Optimization for DOT (DRO-DOT), is proposed. It reduces the dimensionality of the inverse DOT problem by reducing the number of unknowns in two steps and thereby makes the overall process fast. First, it constructs a low resolution voxel basis based on the sensing-matrix properties to find an image support. Second, it reconstructs the sparse image inside this support. To compensate for the reduced sensitivity with increasing depth, depth compensation is incorporated in DRO-DOT. An efficient method to optimally select the regularization parameter is proposed for obtaining a high-quality DOT image. DRO-DOT is also able to reconstruct high-resolution images even with a limited number of optodes in a spatially limited imaging set-up. PMID:26940661
Ye, Chunfang; Kamysiak, Keith T; Sullivan, Amy C; McLeod, Robert R
2012-03-12
We demonstrate single-mode uniform and parabolically tapered three-dimensional waveguides fabricated via direct-write lithography in diffusion-based photopolymers. Modulation of the writing power is shown to compensate Beer-Lambert absorption in the single-photon initiator and to provide precise control of modal tapers. A laminated sample preparation is introduced to enable full 3D characterization of these modal tapers without the need for sample polishing which is difficult for this class of polymer. The accuracy and repeatability of this modal characterization is shown to allow precise measurement of propagation loss from single samples. These testing procedures are used to demonstrate single-mode waveguides with 0.147 dB/cm excess propagation loss and symmetrical tapers up to 1:2.5 using 1.5 microwatts of continuous write power. PMID:22418540
Chervin, Christopher N; Parker, Joseph F; Nelson, Eric S; Rolison, Debra R; Long, Jeffrey W
2016-04-29
The ability to effectively screen and validate gas-diffusion electrodes is critical to the development of next-generation metal-air batteries and regenerative fuel cells. The limiting electrode in a classic two-terminal device such as a battery or fuel cell is difficult to discern without an internal reference electrode, but the flooded electrolyte characteristic of three-electrode electroanalytical cells negates the prime function of an air electrode-a void volume freely accessible to gases. The nanostructured catalysts that drive the energy-conversion reactions (e.g., oxygen reduction and evolution in the air electrode of metal-air batteries) are best evaluated in the electrode structure as-used in the practical device. We have designed, 3D-printed, and characterized an air-breathing, thermodynamically referenced electroanalytical cell that allows us to mimic the Janus arrangement of the gas-diffusion electrode in a metal-air cell: one face freely exposed to gases, the other wetted by electrolyte. PMID:26987282
NASA Astrophysics Data System (ADS)
Chervin, Christopher N.; Parker, Joseph F.; Nelson, Eric S.; Rolison, Debra R.; Long, Jeffrey W.
2016-04-01
The ability to effectively screen and validate gas-diffusion electrodes is critical to the development of next-generation metal-air batteries and regenerative fuel cells. The limiting electrode in a classic two-terminal device such as a battery or fuel cell is difficult to discern without an internal reference electrode, but the flooded electrolyte characteristic of three-electrode electroanalytical cells negates the prime function of an air electrode—a void volume freely accessible to gases. The nanostructured catalysts that drive the energy-conversion reactions (e.g., oxygen reduction and evolution in the air electrode of metal-air batteries) are best evaluated in the electrode structure as-used in the practical device. We have designed, 3D-printed, and characterized an air-breathing, thermodynamically referenced electroanalytical cell that allows us to mimic the Janus arrangement of the gas-diffusion electrode in a metal-air cell: one face freely exposed to gases, the other wetted by electrolyte.
Lockwood, Sarah Y; Meisel, Jayda E; Monsma, Frederick J; Spence, Dana M
2016-02-01
The process of bringing a drug to market involves many steps, including the preclinical stage, where various properties of the drug candidate molecule are determined. These properties, which include drug absorption, distribution, metabolism, and excretion, are often displayed in a pharmacokinetic (PK) profile. While PK profiles are determined in animal models, in vitro systems that model in vivo processes are available, although each possesses shortcomings. Here, we present a 3D-printed, diffusion-based, and dynamic in vitro PK device. The device contains six flow channels, each with integrated porous membrane-based insert wells. The pores of these membranes enable drugs to freely diffuse back and forth between the flow channels and the inserts, thus enabling both loading and clearance portions of a standard PK curve to be generated. The device is designed to work with 96-well plate technology and consumes single-digit milliliter volumes to generate multiple PK profiles, simultaneously. Generation of PK profiles by use of the device was initially performed with fluorescein as a test molecule. Effects of such parameters as flow rate, loading time, volume in the insert well, and initial concentration of the test molecule were investigated. A prediction model was generated from this data, enabling the user to predict the concentration of the test molecule at any point along the PK profile within a coefficient of variation of ∼ 5%. Depletion of the analyte from the well was characterized and was determined to follow first-order rate kinetics, indicated by statistically equivalent (p > 0.05) depletion half-lives that were independent of the starting concentration. A PK curve for an approved antibiotic, levofloxacin, was generated to show utility beyond the fluorescein test molecule. PMID:26727249
NASA Technical Reports Server (NTRS)
Asenov, Asen; Brown, A. R.; Slavcheva, G.; Davies, J. H.
2000-01-01
When MOSFETs are scaled to deep submicron dimensions the discreteness and randomness of the dopant charges in the channel region introduces significant fluctuations in the device characteristics. This effect, predicted 20 year ago, has been confirmed experimentally and in simulation studies. The impact of the fluctuations on the functionality, yield, and reliability of the corresponding systems shifts the paradigm of the numerical device simulation. It becomes insufficient to simulate only one device representing one macroscopical design in a continuous charge approximation. An ensemble of macroscopically identical but microscopically different devices has to be characterized by simulation of statistically significant samples. The aims of the numerical simulations shift from predicting the characteristics of a single device with continuous doping towards estimating the mean values and the standard deviations of basic design parameters such as threshold voltage, subthreshold slope, transconductance, drive current, etc. for the whole ensemble of 'atomistically' different devices in the system. It has to be pointed out that even the mean values obtained from 'atomistic' simulations are not identical to the values obtained from continuous doping simulations. In this paper we present a hierarchical approach to the 'atomistic' simulation of aggressively scaled decanano MOSFETs. A full scale 3D drift-diffusion'atomostic' simulation approach is first described and used for verification of the more economical, but also more restricted, options. To reduce the processor time and memory requirements at high drain voltage we have developed a self-consistent option based on a thin slab solution of the current continuity equation only in the channel region. This is coupled to the Poisson's equation solution in the whole simulation domain in the Gummel iteration cycles. The accuracy of this approach is investigated in comparison with the full self-consistent solution. At low drain
NASA Technical Reports Server (NTRS)
Agrawal, Ajay K.; Yang, Tah-Teh
1993-01-01
This paper describes the 3D computations of a flow field in the compressor/combustor diffusers of an industrial gas turbine. The geometry considered includes components such as the combustor support strut, the transition piece and the impingement sleeve with discrete cooling air holes on its surface. Because the geometry was complex and 3D, the airflow path was divided into two computational domains sharing an interface region. The body-fitted grid was generated independently in each of the two domains. The governing equations for incompressible Navier-Stokes equations were solved using the finite volume approach. The results show that the flow in the prediffuser is strongly coupled with the flow in the dump diffuser and vice versa. The computations also revealed that the flow in the dump diffuser is highly nonuniform.
MCNP: Multigroup/adjoint capabilities
Wagner, J.C.; Redmond, E.L. II; Palmtag, S.P.; Hendricks, J.S.
1994-04-01
This report discusses various aspects related to the use and validity of the general purpose Monte Carlo code MCNP for multigroup/adjoint calculations. The increased desire to perform comparisons between Monte Carlo and deterministic codes, along with the ever-present desire to increase the efficiency of large MCNP calculations has produced a greater user demand for the multigroup/adjoint capabilities. To more fully utilize these capabilities, we review the applications of the Monte Carlo multigroup/adjoint method, describe how to generate multigroup cross sections for MCNP with the auxiliary CRSRD code, describe how to use the multigroup/adjoint capability in MCNP, and provide examples and results indicating the effectiveness and validity of the MCNP multigroup/adjoint treatment. This information should assist users in taking advantage of the MCNP multigroup/adjoint capabilities.
NASA Astrophysics Data System (ADS)
Tiwari, Anand P.; Yoo, Heejoun; Lee, Jeongtaik; Kim, Doyoung; Park, Jong Hyeok; Lee, Hyoyoung
2015-07-01
We report new three-dimensional (3D)-nanostructured MoS2-carbonaceous materials in which MoS2 sheets are intercalated between the graphite layers that possess a multiply repeated graphite/MoS2/graphite structure which prevents the aggregation of MoS2 and diffusion of sulfur from carbonaceous materials, enhancing the cycling stability of Li-ion batteries. We developed an efficient and scalable process applicable to mass production for synthesizing non-aggregated MoS2-intercalated 3D hybrid-nanostructured graphite based on stress induced and microwave irradiation. X-ray diffraction, X-ray photospectroscopy, Raman spectroscopy, field emission scanning electron microscopy, and high-resolution transmission electron microscopy analyses demonstrated that the as-synthesized materials consisted of MoS2-intercalated 3D hybrid-nanostructured graphite platelets that had a multiply repeated graphite/MoS2/graphite structure. The obtained MoS2-graphite powder surpasses MoS2 as an anode material in terms of specific capacity, cyclic stability, and rate performances at high current densities for Li-ion batteries. The electrochemical impedance spectroscopy demonstrated that the graphite sheets not only reduced the contact resistance in the electrode but also facilitated electron transfer in the lithiation/delithiation processes. The superior electrochemical performances especially for the cycling stability of the Li-ion battery originate from prevention of the sulfur diffusion of the MoS2-intercalated 3D-nanostructured graphite.We report new three-dimensional (3D)-nanostructured MoS2-carbonaceous materials in which MoS2 sheets are intercalated between the graphite layers that possess a multiply repeated graphite/MoS2/graphite structure which prevents the aggregation of MoS2 and diffusion of sulfur from carbonaceous materials, enhancing the cycling stability of Li-ion batteries. We developed an efficient and scalable process applicable to mass production for synthesizing non
Adluru, Nagesh; Lee, Jee Eun; Lazar, Mariana; Lainhart, Janet E.; Alexander, Andrew L.
2011-01-01
We present a novel cosine series representation for encoding fiber bundles consisting of multiple 3D curves. The coordinates of curves are parameterized as coefficients of cosine series expansion. We address the issue of registration, averaging and statistical inference on curves in a unified Hilbert space framework. Unlike traditional splines, the proposed method does not have internal knots and explicitly represents curves as a linear combination of cosine basis. This simplicity in the representation enables us to design statistical models, register curves and perform subsequent analysis in a more unified statistical framework than splines. The proposed representation is applied in characterizing abnormal shape of white matter fiber tracts passing through the splenium of the corpus callosum in autistic subjects. For an arbitrary tract, a 19 degree expansion is usually found to be sufficient to reconstruct the tract with 60 parameters. PMID:23316267
Ivanov, Konstantin L. Lukzen, Nikita N.; Sadovsky, Vladimir M.
2015-08-28
In this work, we treat spin-selective recombination of a geminate radical pair (RP) in a spherical “microreactor,” i.e., of a RP confined in a micelle, vesicle, or liposome. We consider the microreactor model proposed earlier, in which one of the radicals is located at the center of the micelle and the other one undergoes three-dimensional diffusion inside the micelle. In addition, we suggest a two-dimensional model, in which one of the radicals is located at the “pole” of the sphere, while the other one diffuses on the spherical surface. For this model, we have obtained a general analytical expression for the RP recombination yield in terms of the free Green function of two-dimensional diffusion motion. In turn, this Green function is expressed via the Legendre functions and thus takes account of diffusion over a restricted spherical surface and its curvature. The obtained expression allows one to calculate the RP recombination efficiency at an arbitrary magnetic field strength. We performed a comparison of the two models taking the same geometric parameters (i.e., the microreactor radius and the closest approach distance of the radicals), chemical reactivity, magnetic interactions in the RP and diffusion coefficient. Significant difference between the predictions of the two models is found, which is thus originating solely from the dimensionality effect: for different dimensionality of space, the statistics of diffusional contacts of radicals becomes different altering the reaction yield. We have calculated the magnetic field dependence of the RP reaction yield and chemically induced dynamic nuclear polarization of the reaction products at different sizes of the microreactor, exchange interaction, and spin relaxation rates. Interestingly, due to the intricate interplay of diffusional contacts of reactants and spin dynamics, the dependence of the reaction yield on the microreactor radius is non-monotonous. Our results are of importance for (i) interpreting
NASA Astrophysics Data System (ADS)
Ivanov, Konstantin L.; Sadovsky, Vladimir M.; Lukzen, Nikita N.
2015-08-01
In this work, we treat spin-selective recombination of a geminate radical pair (RP) in a spherical "microreactor," i.e., of a RP confined in a micelle, vesicle, or liposome. We consider the microreactor model proposed earlier, in which one of the radicals is located at the center of the micelle and the other one undergoes three-dimensional diffusion inside the micelle. In addition, we suggest a two-dimensional model, in which one of the radicals is located at the "pole" of the sphere, while the other one diffuses on the spherical surface. For this model, we have obtained a general analytical expression for the RP recombination yield in terms of the free Green function of two-dimensional diffusion motion. In turn, this Green function is expressed via the Legendre functions and thus takes account of diffusion over a restricted spherical surface and its curvature. The obtained expression allows one to calculate the RP recombination efficiency at an arbitrary magnetic field strength. We performed a comparison of the two models taking the same geometric parameters (i.e., the microreactor radius and the closest approach distance of the radicals), chemical reactivity, magnetic interactions in the RP and diffusion coefficient. Significant difference between the predictions of the two models is found, which is thus originating solely from the dimensionality effect: for different dimensionality of space, the statistics of diffusional contacts of radicals becomes different altering the reaction yield. We have calculated the magnetic field dependence of the RP reaction yield and chemically induced dynamic nuclear polarization of the reaction products at different sizes of the microreactor, exchange interaction, and spin relaxation rates. Interestingly, due to the intricate interplay of diffusional contacts of reactants and spin dynamics, the dependence of the reaction yield on the microreactor radius is non-monotonous. Our results are of importance for (i) interpreting
Ivanov, Konstantin L; Sadovsky, Vladimir M; Lukzen, Nikita N
2015-08-28
In this work, we treat spin-selective recombination of a geminate radical pair (RP) in a spherical "microreactor," i.e., of a RP confined in a micelle, vesicle, or liposome. We consider the microreactor model proposed earlier, in which one of the radicals is located at the center of the micelle and the other one undergoes three-dimensional diffusion inside the micelle. In addition, we suggest a two-dimensional model, in which one of the radicals is located at the "pole" of the sphere, while the other one diffuses on the spherical surface. For this model, we have obtained a general analytical expression for the RP recombination yield in terms of the free Green function of two-dimensional diffusion motion. In turn, this Green function is expressed via the Legendre functions and thus takes account of diffusion over a restricted spherical surface and its curvature. The obtained expression allows one to calculate the RP recombination efficiency at an arbitrary magnetic field strength. We performed a comparison of the two models taking the same geometric parameters (i.e., the microreactor radius and the closest approach distance of the radicals), chemical reactivity, magnetic interactions in the RP and diffusion coefficient. Significant difference between the predictions of the two models is found, which is thus originating solely from the dimensionality effect: for different dimensionality of space, the statistics of diffusional contacts of radicals becomes different altering the reaction yield. We have calculated the magnetic field dependence of the RP reaction yield and chemically induced dynamic nuclear polarization of the reaction products at different sizes of the microreactor, exchange interaction, and spin relaxation rates. Interestingly, due to the intricate interplay of diffusional contacts of reactants and spin dynamics, the dependence of the reaction yield on the microreactor radius is non-monotonous. Our results are of importance for (i) interpreting
Terryn, Christine; Garnotel, Roselyne; Jeannesson, Pierre; Sockalingum, Ganesh D.; Manfait, Michel; Perraut, François; Dinten, Jean-Marc; Koenig, Anne; Piot, Olivier
2016-01-01
During aging, alterations of extracellular matrix proteins contribute to various pathological phenotypes. Among these alterations, type I collagen cross-linking and associated glycation products accumulation over time detrimentally affects its physico-chemical properties, leading to alterations of tissue biomechanical stability. Here, different-age collagen 3D matrices using non-destructive and label-free biophotonic techniques were analysed to highlight the impact of collagen I aging on 3D constructs, at macroscopic and microscopic levels. Matrices were prepared with collagens extracted from tail tendons of rats (newborns, young and old adults) to be within the physiological aging process. The data of diffuse reflectance spectroscopy reveal that aging leads to an inhibition of fibril assembly and a resulting decrease of gel density. Investigations by confocal reflectance microscopy highlight poor-fibrillar structures in oldest collagen networks most likely related to the glycation products accumulation. Complementarily, an infrared analysis brings out marked spectral variations in the Amide I profile, specific of the peptidic bond conformation and for carbohydrates vibrations as function of collagen-age. Interestingly, we also highlight an unexpected behavior for newborn collagen, exhibiting poorly-organized networks and microscopic features close to the oldest collagen. These results demonstrate that changes in collagen optical properties are relevant for investigating the incidence of aging in 3D matrix models. PMID:26885896
Zemskova, Varvara; Garaud, Pascale; Deal, Morgan; Vauclair, Sylvie
2014-11-10
Iron-rich layers are known to form in the stellar subsurface through a combination of gravitational settling and radiative levitation. Their presence, nature, and detailed structure can affect the excitation process of various stellar pulsation modes and must therefore be modeled carefully in order to better interpret Kepler asteroseismic data. In this paper, we study the interplay between atomic diffusion and fingering convection in A-type stars, as well as its role in the establishment and evolution of iron accumulation layers. To do so, we use a combination of three-dimensional idealized numerical simulations of fingering convection (which neglect radiative transfer and complex opacity effects) and one-dimensional realistic stellar models. Using the three-dimensional simulations, we first validate the mixing prescription for fingering convection recently proposed by Brown et al. (within the scope of the aforementioned approximation) and identify what system parameters (total mass of iron, iron diffusivity, thermal diffusivity, etc.) play a role in the overall evolution of the layer. We then implement the Brown et al. prescription in the Toulouse-Geneva Evolution Code to study the evolution of the iron abundance profile beneath the stellar surface. We find, as first discussed by Théado et al., that when the concurrent settling of helium is ignored, this accumulation rapidly causes an inversion in the mean molecular weight profile, which then drives fingering convection. The latter mixes iron with the surrounding material very efficiently, and the resulting iron layer is very weak. However, taking helium settling into account partially stabilizes the iron profile against fingering convection, and a large iron overabundance can accumulate. The opacity also increases significantly as a result, and in some cases it ultimately triggers dynamical convection. The direct effects of radiative acceleration on the dynamics of fingering convection (especially in the
NASA Technical Reports Server (NTRS)
Zhang, Jun; Ge, Lixin; Kouatchou, Jules
2000-01-01
A new fourth order compact difference scheme for the three dimensional convection diffusion equation with variable coefficients is presented. The novelty of this new difference scheme is that it Only requires 15 grid points and that it can be decoupled with two colors. The entire computational grid can be updated in two parallel subsweeps with the Gauss-Seidel type iterative method. This is compared with the known 19 point fourth order compact differenCe scheme which requires four colors to decouple the computational grid. Numerical results, with multigrid methods implemented on a shared memory parallel computer, are presented to compare the 15 point and the 19 point fourth order compact schemes.
NASA Astrophysics Data System (ADS)
Hudspeth, J. M.; Goossens, D. J.
2012-01-01
A new vial-in-vial vapour diffusion method for growing single crystals of fully deuterated triglycine sulphate (TGS) has been developed. Single crystals of hydrogenous TGS were also grown for comparison purposes. The crystals have been characterised using x-ray diffraction and differential scanning calorimetry. The phase transition temperature was 334.0±0.5 K for fully deuterated TGS compared to 322.3±0.3 K for hydrogenous TGS. These values compare well with the expected TC.
3D He-3 diffusion MRI as a local in vivo morphometric tool to evaluate emphysematous rat lungs
Jacob, Rick E.; Minard, Kevin R.; Laicher, Gernot J.; Timchalk, Charles
2008-08-21
In this work, we validate 3He magnetic resonance imaging as a non-invasive morphometric tool to assess emphysematous disease state on a local level. Emphysema was induced intratracheally in rats with 25U/100g body weight of porcine pancreatic elastase dissolved in 200 μL saline. Rats were then paired with saline-dosed controls. Nine three-dimensional 3He diffusion-weighted images were acquired at one-, two-, or three-weeks post-dose, after which the lungs were harvested and prepared for histological analysis. Recently introduced indices sensitive to the heterogeneity of the airspace size distribution were calculated. These indices, D1 and D2, were derived from the moments of the mean equivalent airway diameters. Averaged over the entire lung, it is shown that the 3He diffusivity (Dave) and anisotropy (Dan) both correlate with histology (R = 0.85, p < 0.0001 and R = 0.88, p < 0.0001, respectively). By matching small (0.046 cm2) regions in 3He images with corresponding regions in histological slices, Dave and Dan each correlate significantly with both D1 and D2 (R = 0.93, p < 0.0001). It is concluded that 3He MRI is a viable non-invasive morphometric tool for localized in vivo emphysema assessment.
NASA Astrophysics Data System (ADS)
Zhang, Weian; Wang, Long; Dong, Qixin
2011-06-01
The omni-directional laser warning equipment based on infrared fish-eye lens and short-wave infrared FPA has been used to protect large-scale targets, which can detect the threat laser scattered by the attacked targets or the objects surrounding them, and image the laser spot on FPA, then fix the position of spot. The application offsets the disadvantage of direct interception warner which need disposed largely. Before study of imaging mechanism about the scattered laser spot, the definition of geometry relationship is needed firstly. In this paper we developed a 3D geometry model by analyzing the position relationships in typical battlefield environment among the enemy's threat laser source, the laser spot radiated on one flat surface and our omni-directional laser warning fish-eye lens. The model including R, α, β, d, θ, φ, ψ, δ etc. 8 parameters and 4 coordinate systems was suitable for any general situations. After achievement of the model foundation, we obtained analytic expression of the laser spot contour on flat surface, then attained analytic expression of spot contour on image surface by calculating the object space half-field angle and the azimuth angle relative to fish-eye lens of an arbitrary point at the spot edge on flat surface. The attainment of the expression makes possible that we can analyze the spot energy distributions on image surface and the imaging characteristic of the scattered laser spot via fish-eye lens, then can compute the transmission direction of the threat laser. The foundation of the model in this paper has an importantly basic and guiding meaning to the latter research on this aspect.
Kim, Tae-Yoo; Son, Hwa-Jin; Lim, Seung-Kyu; Song, Young-Il; Park, Hwa-Sun; Suh, Su-Jeong
2014-12-01
Electroless Ni-P films were investigated with the aim of application as barrier and seed layers in 3D interconnect technology. Different shapes of blind-via holes were fabricated with a deep reactive ion etcher and SiO2 formed on these holes as an insulating layer. The surface of the substrate has been made hydrophilic by O2 plasma treatment with 100 W of power for 20 min. Electroless Ni-P films were deposited as both a diffusion barrier and a seed layer for Cu filling process. Prior to plating, substrates were activated in a palladium chloride solution after sensitization in a tin chloride solution with various conditions in order to deposit uniform films in TSV. After the formation of the electroless barrier layer, electro Cu was plated directly on the barrier layer. Ni-P films fabricated in blind-via holes were observed by scanning electron microscope. Energy dispersive spectroscopy line scanning was carried out for evaluating the diffusion barrier properties of the Ni-P films. The electroless Ni-P layer worked well as a Cu diffusion barrier until 300 degrees C. However, Cu ions diffused into barrier layer when the annealing temperature increases over 400 degrees C. PMID:25971093
The Suppression of Energy Discretization Errors in Multigroup Transport Calculations
Larsen, Edward
2013-06-17
The Objective of this project is to develop, implement, and test new deterministric methods to solve, as efficiently as possible, multigroup neutron transport problems having an extremely large number of groups. Our approach was to (i) use the standard CMFD method to "coarsen" the space-angle grid, yielding a multigroup diffusion equation, and (ii) use a new multigrid-in-space-and-energy technique to efficiently solve the multigroup diffusion problem. The overall strategy of (i) how to coarsen the spatial an energy grids, and (ii) how to navigate through the various grids, has the goal of minimizing the overall computational effort. This approach yields not only the fine-grid solution, but also coarse-group flux-weighted cross sections that can be used for other related problems.
NASA Astrophysics Data System (ADS)
Monga, O.; Garnier, P.; Pot, V.; Coucheney, E.; Nunan, N.; Otten, W.; Chenu, C.
2013-10-01
This paper deals with the simulation of microbial degradation in soil within pore space at microscopic scale. Pore space was described using sphere network coming from a geometrical modeling algorithm. The biological model was improved regarding previous work in order to include transformation of dissolved organic compounds and diffusion processes. Our model was tested using experimental results of a simple substrate decomposition (Fructose) within a simple media (the sand). Diverse microbial communities were inoculated. Separated incubations in microcosms were carried out using 5 different bacterial communities at 2 different water potentials of -10 cm and -100 cm of water. We calibrated the biological parameters by means of experimental data obtained at high water content and we tested the model without any parameters change at low water content. Same as for experimental data, our simulation results showed the decrease in water content involved the decrease of mineralisation. The model was able to simulate the decrease of connectivity between substrate and microorganism due the decrease of water content.
NASA Astrophysics Data System (ADS)
Monga, O.; Garnier, P.; Pot, V.; Coucheney, E.; Nunan, N.; Otten, W.; Chenu, C.
2014-04-01
This paper deals with the simulation of microbial degradation of organic matter in soil within the pore space at a microscopic scale. Pore space was analysed with micro-computed tomography and described using a sphere network coming from a geometrical modelling algorithm. The biological model was improved regarding previous work in order to include the transformation of dissolved organic compounds and diffusion processes. We tested our model using experimental results of a simple substrate decomposition experiment (fructose) within a simple medium (sand) in the presence of different bacterial strains. Separate incubations were carried out in microcosms using five different bacterial communities at two different water potentials of -10 and -100 cm of water. We calibrated the biological parameters by means of experimental data obtained at high water content, and we tested the model without changing any parameters at low water content. Same as for the experimental data, our simulation results showed that the decrease in water content caused a decrease of mineralization rate. The model was able to simulate the decrease of connectivity between substrate and microorganism due the decrease of water content.
Doll, A; Abu Eid, M; Kehrli, P; Esposito, P; Gillis, C; Bogorin, A; Jacques, C; Dietemann, J L
2000-06-01
We propose to assess the usefulness of diffusion-weighted MR Imaging (DWI), fluid-attenuated inversion recovery (FLAIR) and constructive interference in steady state (CISS) sequences in depicting epidermoid cysts (EC). FLAIR, CISS and DWI were obtained in 7 patients among 22. All patients were studied with T1 and T2 sequences. On Spin Echo images, EC demonstrate signal similar to LCS, which may lead to difficult differentiation between EC and arachnoid cyst (AC), specially for inexperienced radiologists. EC appear with a heterogeneous signal on T1 images (32%), irregular limits (91%) and with extension through foramen of Pacchioni in 18% of cases. On FLAIR sequence, the tumors were heterogeneous, different from void signal of CSF in 86% of cases. On CISS sequence, the tumors appear heterogeneous, hyperintense but less than LCS and with irregular limits in all cases. Some more, CISS images allowed to appreciate exact tumor extension and their relations with nerves and vessels. On DWI images, signal is hyperintense in all cases. Our study exhibited the great usefulness of DWI, CISS and FLAIR sequences in diagnosis of EC and in differentiating EC from AC. PMID:10970961
NASA Astrophysics Data System (ADS)
Cunningham, G.; Tu, W.; Morley, S.; Chen, Y.; Haidecuk, J.; De Pascuale, S.; Kletzing, C.
2014-12-01
Modeling the variation of the MeV electron phase space density in the inner magnetosphere during active times is sensitive to many parameters, including the initial and time-varying boundary conditions, VLF wave spectral properties, plasma density, and magnetic field. Historically, diffusion codes like LANL's DREAM3D have relied on the statistically-derived dependence of these parameters on geomagnetic indices, e.g. the wave intensity as a function of the AE index. However, the large number of satellites currently sampling the inner magnetosphere presents modelers with an unparalleled opportunity to create 'event-specific' models for many of these parameters. Toward this goal, we recently showed that using an event-specific model of the chorus wave intensity, built from proxy observations of low-energy electron precipitation observed by POES, along with a low-energy time-varying boundary condition informed by the Van Allen Probes, allows DREAM3D to reproduce the large enhancement of PSD for MeV electrons observed during the October 8-9, 2012, storm. One major limitation of this work is the fact that we used the static Sheeley plasma density model and a dipole magnetic field. Here we will discuss new results that use measurements of the plasma density inferred from the Van Allen Probes' EMFISIS instrument to build an event-specific, global, time-dependent model of the plasma density that we use in DREAM3D in combination with the Tsyganenko 2004 storm-time model of the magnetic field. We show that this combination of plasma density and magnetic field model reproduce the ratio of cyclotron frequency to plasma frequency reported by EMFISIS during the entirety of the October 8-9, 2012, storm at all L-shells of interest, whereas our earlier results did not use the correct ratio at most locations and times. Because this ratio is a key parameter governing the effectiveness of chorus waves in accelerating electrons to higher energy, our new DREAM3D results resolve several
NASA Astrophysics Data System (ADS)
Cunningham, G.; Tu, W.; Chen, Y.; Reeves, G. D.; Henderson, M. G.; Baker, D. N.; Blake, J. B.; Spence, H.
2013-12-01
During the interval October 8-9, 2012, the phase-space density (PSD) of high-energy electrons exhibited a dropout preceding an intense enhancement observed by the MagEIS and REPT instruments aboard the Van Allen Probes. The evolution of the PSD suggests heating by chorus waves, which were observed to have high intensities at the time of the enhancement [1]. Although intense chorus waves were also observed during the first Dst dip on October 8, no PSD enhancement was observed at this time. We demonstrate a quantitative reproduction of the entire event that makes use of three recent modifications to the LANL DREAM3D diffusion code: 1) incorporation of a time-dependent, low-energy, boundary condition from the MagEIS instrument, 2) use of a time-dependent estimate of the chorus wave intensity derived from observations of POES low-energy electron precipitation, and 3) use of an estimate of the last closed drift shell, beyond which electrons are assumed to have a lifetime that is proportional to their drift period around earth. The key features of the event are quantitatively reproduced by the simulation, including the dropout on October 8, and a rapid increase in PSD early on October 9, with a peak near L*=4.2. The DREAM3D code predicts the dropout on October 8 because this feature is dominated by magnetospheric compression and outward radial diffusion-the L* of the last closed drift-shell reaches a minimum value of 5.33 at 1026 UT on October 8. We find that a ';statistical' wave model based on historical CRRES measurements binned in AE* does not reproduce the enhancement because the peak wave amplitudes are only a few 10's of pT, whereas an ';event-specific' model reproduces both the magnitude and timing of the enhancement very well, a s shown in the Figure, because the peak wave amplitudes are 10x higher. [1] 'Electron Acceleration in the Heart of the Van Allen Radiation Belts', G. D. Reeves et al., Science 1237743, Published online 25 July 2013 [DOI:10.1126/science
NASA Astrophysics Data System (ADS)
He, Jie; Weersink, Robert; Veilleux, Israel; Mayo, Kenwrick; Zhang, Anqi; Piao, Daqing; Alam, Adeel; Trachtenberg, John; Wilson, Brian C.
2013-03-01
Interstitial near-infrared laser thermal therapy (LITT) is currently undergoing clinical trials as an alternative to watchful waiting or radical surgery in patients with low-risk focal prostate cancer. Currently, we use magnetic resonance image (MRI)-based thermography to monitor treatment delivery and determine indirectly the completeness of the target tissue destruction while avoiding damage to adjacent normal tissues, particularly the rectal wall. However, incomplete tumor destruction has occurred in a significant fraction of patients due to premature termination of treatment, since the photocoagulation zone is not directly observed. Hence, we are developing transrectal diffuse optical tomography (TRDOT), in combination with transrectal 3D ultrasound (3D-TRUS), to address his limitation. This is based on the large changes in optical scattering expected upon tissue coagulation. Here, we present forward simulations of a growing coagulated lesion with optical scattering contrast, using an established finite element analysis software platform (NIRFAST). The simulations were validated in tissue-simulating phantoms, with measurements acquired by a state-of-the-art continuous wave (CW) TRDOT system and a recently assembled bench-top CW-DOT system, with specific source-detector configurations. Two image reconstruction schemes were investigated and evaluated, specifically for the accurate delineation of the posterior boundary of the coagulation zone as the critical parameter for treatment guidance in this clinical application.
NASA Astrophysics Data System (ADS)
Bergeon, N.; Mota, F. L.; Chen, L.; Tourret, D.; Debierre, J. M.; Guérin, R.; Karma, A.; Billia, B.; Trivedi, R.
2015-06-01
To clarify and characterize the fundamental physical mechanisms active in the dynamical formation of three-dimensional (3D) arrays of cells and dendrites under diffusive growth conditions, in situ monitoring of series of experiments on transparent model alloy succinonitrile - 0.24 wt% camphor was carried out under low gravity in the DECLIC Directional Solidification Insert on-board the International Space Station. These experiments offered the very unique opportunity to in situ observe and characterize the whole development of the microstructure in extended 3D patterns. The experimental methods will be first briefly described, including in particular the observation modes and the image analysis procedures developed to quantitatively characterize the patterns. Microgravity environment provided the conditions to get quantitative benchmark data: homogeneous patterns corresponding to homogeneous values of control parameters along the whole interface were obtained. The sequence of microstructure formation will be presented as well as the evolution of the primary spacing which is one of the most important pattern characteristic. Time evolution of this primary spacing during the microstructure development will be analysed to identify the mechanisms of spacing selection and adjustment; the importance of the macroscopic interfacial curvature will be pointed out.
NASA Technical Reports Server (NTRS)
Bune, Andris V.; Gillies, Donald C.; Lehoczky, Sandor L.
1998-01-01
Numerical simulation of the HgCdTe growth by the vertical Bridgman method was performed using FIDAP finite element code. Double-diffusive melt convection is analyzed, as the primary factor at controls inhomogeneity of the solidified material. Temperature and concentration fields in the model are also coupled via material properties, such as thermal and solutal expansion coefficients with the dependence on both temperature and concentration, and melting temperature evaluation from pseudobinary CdTe-HgTe phase diagram. Experimental measurements were used to obtain temperature boundary conditions. Parametric study of the melt convection dependence on the gravity conditions was undertaken. It was found, that the maximum convection velocity in the melt can be reduced under certain conditions. Optimal conditions to obtain a near flat solidified interface are discussed. The predicted interface shape is in agreement with one obtained experimentally by quenching. The results of 3-D calculations are compared with previous 2- D findings. A video film featuring 3-D melt convection will be presented.
Gignac, Paul M; Kley, Nathan J; Clarke, Julia A; Colbert, Matthew W; Morhardt, Ashley C; Cerio, Donald; Cost, Ian N; Cox, Philip G; Daza, Juan D; Early, Catherine M; Echols, M Scott; Henkelman, R Mark; Herdina, A Nele; Holliday, Casey M; Li, Zhiheng; Mahlow, Kristin; Merchant, Samer; Müller, Johannes; Orsbon, Courtney P; Paluh, Daniel J; Thies, Monte L; Tsai, Henry P; Witmer, Lawrence M
2016-06-01
Morphologists have historically had to rely on destructive procedures to visualize the three-dimensional (3-D) anatomy of animals. More recently, however, non-destructive techniques have come to the forefront. These include X-ray computed tomography (CT), which has been used most commonly to examine the mineralized, hard-tissue anatomy of living and fossil metazoans. One relatively new and potentially transformative aspect of current CT-based research is the use of chemical agents to render visible, and differentiate between, soft-tissue structures in X-ray images. Specifically, iodine has emerged as one of the most widely used of these contrast agents among animal morphologists due to its ease of handling, cost effectiveness, and differential affinities for major types of soft tissues. The rapid adoption of iodine-based contrast agents has resulted in a proliferation of distinct specimen preparations and scanning parameter choices, as well as an increasing variety of imaging hardware and software preferences. Here we provide a critical review of the recent contributions to iodine-based, contrast-enhanced CT research to enable researchers just beginning to employ contrast enhancement to make sense of this complex new landscape of methodologies. We provide a detailed summary of recent case studies, assess factors that govern success at each step of the specimen storage, preparation, and imaging processes, and make recommendations for standardizing both techniques and reporting practices. Finally, we discuss potential cutting-edge applications of diffusible iodine-based contrast-enhanced computed tomography (diceCT) and the issues that must still be overcome to facilitate the broader adoption of diceCT going forward. PMID:26970556
Gao, Yurui; Parvathaneni, Prasanna; Schilling, Kurt G.; Wang, Feng; Stepniewska, Iwona; Xu, Zhoubing; Choe, Ann S.; Ding, Zhaohua; Gore, John C.; Chen, Li Min; Landman, Bennett A.; Anderson, Adam W.
2016-01-01
Modern magnetic resonance imaging (MRI) brain atlases are high quality 3-D volumes with specific structures labeled in the volume. Atlases are essential in providing a common space for interpretation of results across studies, for anatomical education, and providing quantitative image-based navigation. Extensive work has been devoted to atlas construction for humans, macaque, and several non-primate species (e.g., rat). One notable gap in the literature is the common squirrel monkey – for which the primary published atlases date from the 1960’s. The common squirrel monkey has been used extensively as surrogate for humans in biomedical studies, given its anatomical neuro-system similarities and practical considerations. This work describes the continued development of a multi-modal MRI atlas for the common squirrel monkey, for which a structural imaging space and gray matter parcels have been previously constructed. This study adds white matter tracts to the atlas. The new atlas includes 49 white matter (WM) tracts, defined using diffusion tensor imaging (DTI) in three animals and combines these data to define the anatomical locations of these tracks in a standardized coordinate system compatible with previous development. An anatomist reviewed the resulting tracts and the inter-animal reproducibility (i.e., the Dice index of each WM parcel across animals in common space) was assessed. The Dice indices range from 0.05 to 0.80 due to differences of local registration quality and the variation of WM tract position across individuals. However, the combined WM labels from the 3 animals represent the general locations of WM parcels, adding basic connectivity information to the atlas. PMID:27064328
NASA Astrophysics Data System (ADS)
Gao, Yurui; Parvathaneni, Prasanna; Schilling, Kurt G.; Wang, Feng; Stepniewska, Iwona; Xu, Zhoubing; Choe, Ann S.; Ding, Zhaohua; Gore, John C.; Chen, Li min; Landman, Bennett A.; Anderson, Adam W.
2016-03-01
Modern magnetic resonance imaging (MRI) brain atlases are high quality 3-D volumes with specific structures labeled in the volume. Atlases are essential in providing a common space for interpretation of results across studies, for anatomical education, and providing quantitative image-based navigation. Extensive work has been devoted to atlas construction for humans, macaque, and several non-primate species (e.g., rat). One notable gap in the literature is the common squirrel monkey - for which the primary published atlases date from the 1960's. The common squirrel monkey has been used extensively as surrogate for humans in biomedical studies, given its anatomical neuro-system similarities and practical considerations. This work describes the continued development of a multi-modal MRI atlas for the common squirrel monkey, for which a structural imaging space and gray matter parcels have been previously constructed. This study adds white matter tracts to the atlas. The new atlas includes 49 white matter (WM) tracts, defined using diffusion tensor imaging (DTI) in three animals and combines these data to define the anatomical locations of these tracks in a standardized coordinate system compatible with previous development. An anatomist reviewed the resulting tracts and the inter-animal reproducibility (i.e., the Dice index of each WM parcel across animals in common space) was assessed. The Dice indices range from 0.05 to 0.80 due to differences of local registration quality and the variation of WM tract position across individuals. However, the combined WM labels from the 3 animals represent the general locations of WM parcels, adding basic connectivity information to the atlas.
CASTRO: A NEW COMPRESSIBLE ASTROPHYSICAL SOLVER. III. MULTIGROUP RADIATION HYDRODYNAMICS
Zhang, W.; Almgren, A.; Bell, J.; Howell, L.; Burrows, A.; Dolence, J.
2013-01-15
We present a formulation for multigroup radiation hydrodynamics that is correct to order O(v/c) using the comoving-frame approach and the flux-limited diffusion approximation. We describe a numerical algorithm for solving the system, implemented in the compressible astrophysics code, CASTRO. CASTRO uses a Eulerian grid with block-structured adaptive mesh refinement based on a nested hierarchy of logically rectangular variable-sized grids with simultaneous refinement in both space and time. In our multigroup radiation solver, the system is split into three parts: one part that couples the radiation and fluid in a hyperbolic subsystem, another part that advects the radiation in frequency space, and a parabolic part that evolves radiation diffusion and source-sink terms. The hyperbolic subsystem and the frequency space advection are solved explicitly with high-order Godunov schemes, whereas the parabolic part is solved implicitly with a first-order backward Euler method. Our multigroup radiation solver works for both neutrino and photon radiation.
CASTRO: A New Compressible Astrophysical Solver. III. Multigroup Radiation Hydrodynamics
NASA Astrophysics Data System (ADS)
Zhang, W.; Howell, L.; Almgren, A.; Burrows, A.; Dolence, J.; Bell, J.
2013-01-01
We present a formulation for multigroup radiation hydrodynamics that is correct to order O(v/c) using the comoving-frame approach and the flux-limited diffusion approximation. We describe a numerical algorithm for solving the system, implemented in the compressible astrophysics code, CASTRO. CASTRO uses a Eulerian grid with block-structured adaptive mesh refinement based on a nested hierarchy of logically rectangular variable-sized grids with simultaneous refinement in both space and time. In our multigroup radiation solver, the system is split into three parts: one part that couples the radiation and fluid in a hyperbolic subsystem, another part that advects the radiation in frequency space, and a parabolic part that evolves radiation diffusion and source-sink terms. The hyperbolic subsystem and the frequency space advection are solved explicitly with high-order Godunov schemes, whereas the parabolic part is solved implicitly with a first-order backward Euler method. Our multigroup radiation solver works for both neutrino and photon radiation.
NASA Astrophysics Data System (ADS)
Pletinckx, D.
2011-09-01
The current 3D hype creates a lot of interest in 3D. People go to 3D movies, but are we ready to use 3D in our homes, in our offices, in our communication? Are we ready to deliver real 3D to a general public and use interactive 3D in a meaningful way to enjoy, learn, communicate? The CARARE project is realising this for the moment in the domain of monuments and archaeology, so that real 3D of archaeological sites and European monuments will be available to the general public by 2012. There are several aspects to this endeavour. First of all is the technical aspect of flawlessly delivering 3D content over all platforms and operating systems, without installing software. We have currently a working solution in PDF, but HTML5 will probably be the future. Secondly, there is still little knowledge on how to create 3D learning objects, 3D tourist information or 3D scholarly communication. We are still in a prototype phase when it comes to integrate 3D objects in physical or virtual museums. Nevertheless, Europeana has a tremendous potential as a multi-facetted virtual museum. Finally, 3D has a large potential to act as a hub of information, linking to related 2D imagery, texts, video, sound. We describe how to create such rich, explorable 3D objects that can be used intuitively by the generic Europeana user and what metadata is needed to support the semantic linking.
Procedure to Generate the MPACT Multigroup Library
Kim, Kang Seog
2015-12-17
The CASL neutronics simulator MPACT is under development for the neutronics and T-H coupled simulation for the light water reactor. The objective of this document is focused on reviewing the current procedure to generate the MPACT multigroup library. Detailed methodologies and procedures are included in this document for further discussion to improve the MPACT multigroup library.
Continuous-energy eigenvalue sensitivity coefficient calculations in TSUNAMI-3D
Perfetti, C. M.; Rearden, B. T.
2013-07-01
Two methods for calculating eigenvalue sensitivity coefficients in continuous-energy Monte Carlo applications were implemented in the KENO code within the SCALE code package. The methods were used to calculate sensitivity coefficients for several test problems and produced sensitivity coefficients that agreed well with both reference sensitivities and multigroup TSUNAMI-3D sensitivity coefficients. The newly developed CLUTCH method was observed to produce sensitivity coefficients with high figures of merit and a low memory footprint, and both continuous-energy sensitivity methods met or exceeded the accuracy of the multigroup TSUNAMI-3D calculations. (authors)
3d-3d correspondence revisited
NASA Astrophysics Data System (ADS)
Chung, Hee-Joong; Dimofte, Tudor; Gukov, Sergei; Sułkowski, Piotr
2016-04-01
In fivebrane compactifications on 3-manifolds, we point out the importance of all flat connections in the proper definition of the effective 3d {N}=2 theory. The Lagrangians of some theories with the desired properties can be constructed with the help of homological knot invariants that categorify colored Jones polynomials. Higgsing the full 3d theories constructed this way recovers theories found previously by Dimofte-Gaiotto-Gukov. We also consider the cutting and gluing of 3-manifolds along smooth boundaries and the role played by all flat connections in this operation.
NASA Astrophysics Data System (ADS)
Meulien Ohlmann, Odile
2013-02-01
Today the industry offers a chain of 3D products. Learning to "read" and to "create in 3D" becomes an issue of education of primary importance. 25 years professional experience in France, the United States and Germany, Odile Meulien set up a personal method of initiation to 3D creation that entails the spatial/temporal experience of the holographic visual. She will present some different tools and techniques used for this learning, their advantages and disadvantages, programs and issues of educational policies, constraints and expectations related to the development of new techniques for 3D imaging. Although the creation of display holograms is very much reduced compared to the creation of the 90ies, the holographic concept is spreading in all scientific, social, and artistic activities of our present time. She will also raise many questions: What means 3D? Is it communication? Is it perception? How the seeing and none seeing is interferes? What else has to be taken in consideration to communicate in 3D? How to handle the non visible relations of moving objects with subjects? Does this transform our model of exchange with others? What kind of interaction this has with our everyday life? Then come more practical questions: How to learn creating 3D visualization, to learn 3D grammar, 3D language, 3D thinking? What for? At what level? In which matter? for whom?
Multigroup Reactor Lattice Cell Calculation
Energy Science and Technology Software Center (ESTSC)
1990-03-01
The Winfrith Improved Multigroup Scheme (WIMS), is a general code for reactor lattice cell calculations on a wide range of reactor systems. In particular, the code will accept rod or plate fuel geometries in either regular arrays or in clusters, and the energy group structure has been chosen primarily for thermal calculations. The basic library has been compiled with 14 fast groups, 13 resonance groups and 42 thermal groups, but the user is offered themore » choice of accurate solutions in many groups or rapid calculations in few groups. Temperature dependent thermal scattering matrices for a variety of scattering laws are available in the library for the principal moderators which include hydrogen, deuterium, graphite, beryllium and oxygen. WIMSD5 is a succesor version of WIMS-D/4.« less
ERIC Educational Resources Information Center
Hastings, S. K.
2002-01-01
Discusses 3 D imaging as it relates to digital representations in virtual library collections. Highlights include X-ray computed tomography (X-ray CT); the National Science Foundation (NSF) Digital Library Initiatives; output peripherals; image retrieval systems, including metadata; and applications of 3 D imaging for libraries and museums. (LRW)
A discretization of the multigroup PN radiative transfer equation on general meshes
NASA Astrophysics Data System (ADS)
Hermeline, F.
2016-05-01
We propose and study a finite volume method of discrete duality type for discretizing the multigroup PN approximation of radiative transfer equation on general meshes. This method is second order-accurate on a very large variety of meshes, stable under a Courant-Friedrichs-Lewy condition and it preserves naturally the diffusion asymptotic limit.
NASA Astrophysics Data System (ADS)
Dima, M.; Farisato, G.; Bergomi, M.; Viotto, V.; Magrin, D.; Greggio, D.; Farinato, J.; Marafatto, L.; Ragazzoni, R.; Piazza, D.
2014-08-01
In the last few years 3D printing is getting more and more popular and used in many fields going from manufacturing to industrial design, architecture, medical support and aerospace. 3D printing is an evolution of bi-dimensional printing, which allows to obtain a solid object from a 3D model, realized with a 3D modelling software. The final product is obtained using an additive process, in which successive layers of material are laid down one over the other. A 3D printer allows to realize, in a simple way, very complex shapes, which would be quite difficult to be produced with dedicated conventional facilities. Thanks to the fact that the 3D printing is obtained superposing one layer to the others, it doesn't need any particular work flow and it is sufficient to simply draw the model and send it to print. Many different kinds of 3D printers exist based on the technology and material used for layer deposition. A common material used by the toner is ABS plastics, which is a light and rigid thermoplastic polymer, whose peculiar mechanical properties make it diffusely used in several fields, like pipes production and cars interiors manufacturing. I used this technology to create a 1:1 scale model of the telescope which is the hardware core of the space small mission CHEOPS (CHaracterising ExOPlanets Satellite) by ESA, which aims to characterize EXOplanets via transits observations. The telescope has a Ritchey-Chrétien configuration with a 30cm aperture and the launch is foreseen in 2017. In this paper, I present the different phases for the realization of such a model, focusing onto pros and cons of this kind of technology. For example, because of the finite printable volume (10×10×12 inches in the x, y and z directions respectively), it has been necessary to split the largest parts of the instrument in smaller components to be then reassembled and post-processed. A further issue is the resolution of the printed material, which is expressed in terms of layers
Asymptotic, multigroup flux reconstruction and consistent discontinuity factors
Trahan, Travis J.; Larsen, Edward W.
2015-05-12
Recent theoretical work has led to an asymptotically derived expression for reconstructing the neutron flux from lattice functions and multigroup diffusion solutions. The leading-order asymptotic term is the standard expression for flux reconstruction, i.e., it is the product of a shape function, obtained through a lattice calculation, and the multigroup diffusion solution. The first-order asymptotic correction term is significant only where the gradient of the diffusion solution is not small. Inclusion of this first-order correction term can significantly improve the accuracy of the reconstructed flux. One may define discontinuity factors (DFs) to make certain angular moments of the reconstructed fluxmore » continuous across interfaces between assemblies in 1-D. Indeed, the standard assembly discontinuity factors make the zeroth moment (scalar flux) of the reconstructed flux continuous. The inclusion of the correction term in the flux reconstruction provides an additional degree of freedom that can be used to make two angular moments of the reconstructed flux continuous across interfaces by using current DFs in addition to flux DFs. Thus, numerical results demonstrate that using flux and current DFs together can be more accurate than using only flux DFs, and that making the second angular moment continuous can be more accurate than making the zeroth moment continuous.« less
Asymptotic, multigroup flux reconstruction and consistent discontinuity factors
Trahan, Travis J.; Larsen, Edward W.
2015-05-12
Recent theoretical work has led to an asymptotically derived expression for reconstructing the neutron flux from lattice functions and multigroup diffusion solutions. The leading-order asymptotic term is the standard expression for flux reconstruction, i.e., it is the product of a shape function, obtained through a lattice calculation, and the multigroup diffusion solution. The first-order asymptotic correction term is significant only where the gradient of the diffusion solution is not small. Inclusion of this first-order correction term can significantly improve the accuracy of the reconstructed flux. One may define discontinuity factors (DFs) to make certain angular moments of the reconstructed flux continuous across interfaces between assemblies in 1-D. Indeed, the standard assembly discontinuity factors make the zeroth moment (scalar flux) of the reconstructed flux continuous. The inclusion of the correction term in the flux reconstruction provides an additional degree of freedom that can be used to make two angular moments of the reconstructed flux continuous across interfaces by using current DFs in addition to flux DFs. Thus, numerical results demonstrate that using flux and current DFs together can be more accurate than using only flux DFs, and that making the second angular moment continuous can be more accurate than making the zeroth moment continuous.
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)
GPU-Accelerated Denoising in 3D (GD3D)
Energy Science and Technology Software Center (ESTSC)
2013-10-01
The raw computational power GPU Accelerators enables fast denoising of 3D MR images using bilateral filtering, anisotropic diffusion, and non-local means. This software addresses two facets of this promising application: what tuning is necessary to achieve optimal performance on a modern GPU? And what parameters yield the best denoising results in practice? To answer the first question, the software performs an autotuning step to empirically determine optimal memory blocking on the GPU. To answer themore » second, it performs a sweep of algorithm parameters to determine the combination that best reduces the mean squared error relative to a noiseless reference image.« less
Crandall, K.R.
1987-08-01
TRACE 3-D is an interactive beam-dynamics program that calculates the envelopes of a bunched beam, including linear space-charge forces, through a user-defined transport system. TRACE 3-D provides an immediate graphics display of the envelopes and the phase-space ellipses and allows nine types of beam-matching options. This report describes the beam-dynamics calculations and gives detailed instruction for using the code. Several examples are described in detail.
NASA Astrophysics Data System (ADS)
Seidenberger, K.; Wilhelm, F.; Schmitt, T.; Lehnert, W.; Scholta, J.
Understanding of both water management in PEM fuel cells and degradation mechanisms of the gas diffusion layer (GDL) and their mutual impact is still at least incomplete. Different modelling approaches contribute to gain deeper insight into the processes occurring during fuel cell operation. Considering the GDL, the models can help to obtain information about the distribution of liquid water within the material. Especially, flooded regions can be identified, and the water distribution can be linked to the system geometry. Employed for material development, this information can help to increase the life time of the GDL as a fuel cell component and the fuel cell as the entire system. The Monte Carlo (MC) model presented here helps to simulate and analyse the water household in PEM fuel cell GDLs. This model comprises a three-dimensional, voxel-based representation of the GDL substrate, a section of the flowfield channel and the corresponding rib. Information on the water distribution within the substrate part of the GDL can be estimated.
Si, S.
2012-07-01
The Universal Algorithm of Stiffness Confinement Method (UASCM) for neutron kinetics model of multi-dimensional and multi-group transport equations or diffusion equations has been developed. The numerical experiments based on transport theory code MGSNM and diffusion theory code MGNEM have demonstrated that the algorithm has sufficient accuracy and stability. (authors)
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.
NASA Astrophysics Data System (ADS)
Oldham, Mark
2015-01-01
Radiochromic materials exhibit a colour change when exposed to ionising radiation. Radiochromic film has been used for clinical dosimetry for many years and increasingly so recently, as films of higher sensitivities have become available. The two principle advantages of radiochromic dosimetry include greater tissue equivalence (radiologically) and the lack of requirement for development of the colour change. In a radiochromic material, the colour change arises direct from ionising interactions affecting dye molecules, without requiring any latent chemical, optical or thermal development, with important implications for increased accuracy and convenience. It is only relatively recently however, that 3D radiochromic dosimetry has become possible. In this article we review recent developments and the current state-of-the-art of 3D radiochromic dosimetry, and the potential for a more comprehensive solution for the verification of complex radiation therapy treatments, and 3D dose measurement in general.
NASA Astrophysics Data System (ADS)
Iliesiu, Luca; Kos, Filip; Poland, David; Pufu, Silviu S.; Simmons-Duffin, David; Yacoby, Ran
2016-03-01
We study the conformal bootstrap for a 4-point function of fermions < ψψψψ> in 3D. We first introduce an embedding formalism for 3D spinors and compute the conformal blocks appearing in fermion 4-point functions. Using these results, we find general bounds on the dimensions of operators appearing in the ψ × ψ OPE, and also on the central charge C T . We observe features in our bounds that coincide with scaling dimensions in the GrossNeveu models at large N . We also speculate that other features could coincide with a fermionic CFT containing no relevant scalar operators.
NASA Astrophysics Data System (ADS)
Puszka, Agathe; Di Sieno, Laura; Dalla Mora, Alberto; Pifferi, Antonio; Contini, Davide; Boso, Gianluca; Tosi, Alberto; Hervé, Lionel; Planat-Chrétien, Anne; Koenig, Anne; Dinten, Jean-Marc
2014-02-01
Fiber optic probes with a width limited to a few centimeters can enable diffuse optical tomography (DOT) in intern organs like the prostate or facilitate the measurements on extern organs like the breast or the brain. We have recently shown on 2D tomographic images that time-resolved measurements with a large dynamic range obtained with fast-gated single-photon avalanche diodes (SPADs) could push forward the imaged depth range in a diffusive medium at short source-detector separation compared with conventional non-gated approaches. In this work, we confirm these performances with the first 3D tomographic images reconstructed with such a setup and processed with the Mellin- Laplace transform. More precisely, we investigate the performance of hand-held probes with short interfiber distances in terms of spatial resolution and specifically demonstrate the interest of having a compact probe design featuring small source-detector separations. We compare the spatial resolution obtained with two probes having the same design but different scale factors, the first one featuring only interfiber distances of 15 mm and the second one, 10 mm. We evaluate experimentally the spatial resolution obtained with each probe on the setup with fast-gated SPADs for optical phantoms featuring two absorbing inclusions positioned at different depths and conclude on the potential of short source-detector separations for DOT.
3-D transient analysis of pebble-bed HTGR by TORT-TD/ATTICA3D
Seubert, A.; Sureda, A.; Lapins, J.; Buck, M.; Bader, J.; Laurien, E.
2012-07-01
As most of the acceptance criteria are local core parameters, application of transient 3-D fine mesh neutron transport and thermal hydraulics coupled codes is mandatory for best estimate evaluations of safety margins. This also applies to high-temperature gas cooled reactors (HTGR). Application of 3-D fine-mesh transient transport codes using few energy groups coupled with 3-D thermal hydraulics codes becomes feasible in view of increasing computing power. This paper describes the discrete ordinates based coupled code system TORT-TD/ATTICA3D that has recently been extended by a fine-mesh diffusion solver. Based on transient analyses for the PBMR-400 design, the transport/diffusion capabilities are demonstrated and 3-D local flux and power redistribution effects during a partial control rod withdrawal are shown. (authors)
NASA Astrophysics Data System (ADS)
Iizuka, Keigo
2008-02-01
In order to circumvent the fact that only one observer can view the image from a stereoscopic microscope, an attachment was devised for displaying the 3D microscopic image on a large LCD monitor for viewing by multiple observers in real time. The principle of operation, design, fabrication, and performance are presented, along with tolerance measurements relating to the properties of the cellophane half-wave plate used in the design.
Feasibility of a Multigroup Deterministic Solution Method for 3D Radiotherapy Dose Calculations
Vassiliev, Oleg N.; Wareing, Todd A.; Davis, Ian M.; McGhee, John; Barnett, Douglas; Horton, John L.; Gifford, Kent; Failla, Gregory; Titt, Uwe; Mourtada, Firas
2008-01-01
Purpose To investigate the potential of a novel deterministic solver, Attila, for external photon beam radiotherapy dose calculations. Methods and Materials Two hypothetical cases for prostate and head and neck cancer photon beam treatment plans were calculated using Attila and EGSnrc Monte Carlo simulations. Open beams were modeled as isotropic photon point sources collimated to specified field sizes (100 cm SSD). The sources had a realistic energy spectrum calculated by Monte Carlo for a Varian Clinac 2100 operated in a 6MV photon mode. The Attila computational grids consisted of 106,000 elements, or 424,000 spatial degrees of freedom, for the prostate case, and 123,000 tetrahedral elements, or 492,000 spatial degrees of freedom, for the head and neck cases. Results For both cases, results demonstrate excellent agreement between Attila and EGSnrc in all areas, including the build-up regions, near heterogeneities, and at the beam penumbra. Dose agreement for 99% of the voxels was within 3% (relative point-wise difference) or 3mm distance-to-agreement criterion. Localized differences between the Attila and EGSnrc results were observed at bone and soft tissue interfaces, and are attributable to the effect of voxel material homogenization in calculating dose-to-medium in EGSnrc. For both cases, Attila calculation times were under 20 CPU minutes on a single 2.2 GHz AMD Opteron processor. Conclusions The methods in Attila have the potential to be the basis for an efficient dose engine for patient specific treatment planning, providing accuracy similar to that obtained by Monte Carlo. PMID:18722273
NASA Technical Reports Server (NTRS)
1992-01-01
Ames Research Center research into virtual reality led to the development of the Convolvotron, a high speed digital audio processing system that delivers three-dimensional sound over headphones. It consists of a two-card set designed for use with a personal computer. The Convolvotron's primary application is presentation of 3D audio signals over headphones. Four independent sound sources are filtered with large time-varying filters that compensate for motion. The perceived location of the sound remains constant. Possible applications are in air traffic control towers or airplane cockpits, hearing and perception research and virtual reality development.
NASA Astrophysics Data System (ADS)
Kostrzewski, Andrew A.; Aye, Tin M.; Kim, Dai Hyun; Esterkin, Vladimir; Savant, Gajendra D.
1998-09-01
Physical Optics Corporation has developed an advanced 3-D virtual reality system for use with simulation tools for training technical and military personnel. This system avoids such drawbacks of other virtual reality (VR) systems as eye fatigue, headaches, and alignment for each viewer, all of which are due to the need to wear special VR goggles. The new system is based on direct viewing of an interactive environment. This innovative holographic multiplexed screen technology makes it unnecessary for the viewer to wear special goggles.
Cevidanes, Lucia; Tucker, Scott; Styner, Martin; Kim, Hyungmin; Chapuis, Jonas; Reyes, Mauricio; Proffit, William; Turvey, Timothy; Jaskolka, Michael
2009-01-01
This paper discusses the development of methods for computer-aided jaw surgery. Computer-aided jaw surgery allows us to incorporate the high level of precision necessary for transferring virtual plans into the operating room. We also present a complete computer-aided surgery (CAS) system developed in close collaboration with surgeons. Surgery planning and simulation include construction of 3D surface models from Cone-beam CT (CBCT), dynamic cephalometry, semi-automatic mirroring, interactive cutting of bone and bony segment repositioning. A virtual setup can be used to manufacture positioning splints for intra-operative guidance. The system provides further intra-operative assistance with the help of a computer display showing jaw positions and 3D positioning guides updated in real-time during the surgical procedure. The CAS system aids in dealing with complex cases with benefits for the patient, with surgical practice, and for orthodontic finishing. Advanced software tools for diagnosis and treatment planning allow preparation of detailed operative plans, osteotomy repositioning, bone reconstructions, surgical resident training and assessing the difficulties of the surgical procedures prior to the surgery. CAS has the potential to make the elaboration of the surgical plan a more flexible process, increase the level of detail and accuracy of the plan, yield higher operative precision and control, and enhance documentation of cases. Supported by NIDCR DE017727, and DE018962 PMID:20816308
EXTENSION OF THE 1D FOUR-GROUP ANALYTIC NODAL METHOD TO FULL MULTIGROUP
B. D. Ganapol; D. W. Nigg
2008-09-01
In the mid 80’s, a four-group/two-region, entirely analytical 1D nodal benchmark appeared. It was readily acknowledged that this special case was as far as one could go in terms of group number and still achieve an analytical solution. In this work, we show that by decomposing the solution to the multigroup diffusion equation into homogeneous and particular solutions, extension to any number of groups is a relatively straightforward exercise using the mathematics of linear algebra.
NASA Technical Reports Server (NTRS)
Bune, Andris V.; Gillies, Donald C.; Lehoczky, Sandor L.
1998-01-01
A numerical calculation for a non-dilute alloy solidification was performed using the FIDAP finite element code. For low growth velocities plane front solidification occurs. The location and the shape of the interface was determined using melting temperatures from the HgCdTe liquidus curve. The low thermal conductivity of the solid HgCdTe causes thermal short circuit through the ampoule walls, resulting in curved isotherms in the vicinity of the interface. Double-diffusive convection in the melt is caused by radial temperature gradients and by material density inversion with temperature. Cooling from below and the rejection at the solid-melt interface of the heavier HgTe-rich solute each tend to reduce convection. Because of these complicating factors dimensional rather then non-dimensional modeling was performed. Estimates of convection contributions for various gravity conditions was performed parametrically. For gravity levels higher then 1 0 -7 of earth's gravity it was found that the maximum convection velocity is extremely sensitive to gravity vector orientation and can be reduced at least by factor of 50% for precise orientation of the ampoule in the microgravity environment. The predicted interface shape is in agreement with one obtained experimentally by quenching. The results of 3-D modeling are compared with previous 2-D finding. A video film featuring melt convection will be presented.
Parallel computation of multigroup reactivity coefficient using iterative method
NASA Astrophysics Data System (ADS)
Susmikanti, Mike; Dewayatna, Winter
2013-09-01
One of the research activities to support the commercial radioisotope production program is a safety research target irradiation FPM (Fission Product Molybdenum). FPM targets form a tube made of stainless steel in which the nuclear degrees of superimposed high-enriched uranium. FPM irradiation tube is intended to obtain fission. The fission material widely used in the form of kits in the world of nuclear medicine. Irradiation FPM tube reactor core would interfere with performance. One of the disorders comes from changes in flux or reactivity. It is necessary to study a method for calculating safety terrace ongoing configuration changes during the life of the reactor, making the code faster became an absolute necessity. Neutron safety margin for the research reactor can be reused without modification to the calculation of the reactivity of the reactor, so that is an advantage of using perturbation method. The criticality and flux in multigroup diffusion model was calculate at various irradiation positions in some uranium content. This model has a complex computation. Several parallel algorithms with iterative method have been developed for the sparse and big matrix solution. The Black-Red Gauss Seidel Iteration and the power iteration parallel method can be used to solve multigroup diffusion equation system and calculated the criticality and reactivity coeficient. This research was developed code for reactivity calculation which used one of safety analysis with parallel processing. It can be done more quickly and efficiently by utilizing the parallel processing in the multicore computer. This code was applied for the safety limits calculation of irradiated targets FPM with increment Uranium.
Parallel computation of multigroup reactivity coefficient using iterative method
Susmikanti, Mike; Dewayatna, Winter
2013-09-09
One of the research activities to support the commercial radioisotope production program is a safety research target irradiation FPM (Fission Product Molybdenum). FPM targets form a tube made of stainless steel in which the nuclear degrees of superimposed high-enriched uranium. FPM irradiation tube is intended to obtain fission. The fission material widely used in the form of kits in the world of nuclear medicine. Irradiation FPM tube reactor core would interfere with performance. One of the disorders comes from changes in flux or reactivity. It is necessary to study a method for calculating safety terrace ongoing configuration changes during the life of the reactor, making the code faster became an absolute necessity. Neutron safety margin for the research reactor can be reused without modification to the calculation of the reactivity of the reactor, so that is an advantage of using perturbation method. The criticality and flux in multigroup diffusion model was calculate at various irradiation positions in some uranium content. This model has a complex computation. Several parallel algorithms with iterative method have been developed for the sparse and big matrix solution. The Black-Red Gauss Seidel Iteration and the power iteration parallel method can be used to solve multigroup diffusion equation system and calculated the criticality and reactivity coeficient. This research was developed code for reactivity calculation which used one of safety analysis with parallel processing. It can be done more quickly and efficiently by utilizing the parallel processing in the multicore computer. This code was applied for the safety limits calculation of irradiated targets FPM with increment Uranium.
NASA Technical Reports Server (NTRS)
2004-01-01
The Mars Exploration Rover Spirit took this 3-D navigation camera mosaic of the crater called 'Bonneville' after driving approximately 13 meters (42.7 feet) to get a better vantage point. Spirit's current position is close enough to the edge to see the interior of the crater, but high enough and far enough back to get a view of all of the walls. Because scientists and rover controllers are so pleased with this location, they will stay here for at least two more martian days, or sols, to take high resolution panoramic camera images of 'Bonneville' in its entirety. Just above the far crater rim, on the left side, is the rover's heatshield, which is visible as a tiny reflective speck.
Caspi, S.; Helm, M.; Laslett, L.J.
1991-03-30
We have developed an harmonic representation for the three dimensional field components within the windings of accelerator magnets. The form by which the field is presented is suitable for interfacing with other codes that make use of the 3D field components (particle tracking and stability). The field components can be calculated with high precision and reduced cup time at any location (r,{theta},z) inside the magnet bore. The same conductor geometry which is used to simulate line currents is also used in CAD with modifications more readily available. It is our hope that the format used here for magnetic fields can be used not only as a means of delivering fields but also as a way by which beam dynamics can suggest correction to the conductor geometry. 5 refs., 70 figs.
NASA Astrophysics Data System (ADS)
Gil, José J.; San José, Ignacio
2010-11-01
From our previous definition of the indices of polarimetric purity for 3D light beams [J.J. Gil, J.M. Correas, P.A. Melero and C. Ferreira, Monogr. Semin. Mat. G. de Galdeano 31, 161 (2004)], an analysis of their geometric and physical interpretation is presented. It is found that, in agreement with previous results, the first parameter is a measure of the degree of polarization, whereas the second parameter (called the degree of directionality) is a measure of the mean angular aperture of the direction of propagation of the corresponding light beam. This pair of invariant, non-dimensional, indices of polarimetric purity contains complete information about the polarimetric purity of a light beam. The overall degree of polarimetric purity is obtained as a weighted quadratic average of the degree of polarization and the degree of directionality.
NASA Technical Reports Server (NTRS)
1997-01-01
Many prominent rocks near the Sagan Memorial Station are featured in this image, taken in stereo by the Imager for Mars Pathfinder (IMP) on Sol 3. 3D glasses are necessary to identify surface detail. Wedge is at lower left; Shark, Half-Dome, and Pumpkin are at center. Flat Top, about four inches high, is at lower right. The horizon in the distance is one to two kilometers away.
Mars Pathfinder is the second in NASA's Discovery program of low-cost spacecraft with highly focused science goals. The Jet Propulsion Laboratory, Pasadena, CA, developed and manages the Mars Pathfinder mission for NASA's Office of Space Science, Washington, D.C. JPL is an operating division of the California Institute of Technology (Caltech). The Imager for Mars Pathfinder (IMP) was developed by the University of Arizona Lunar and Planetary Laboratory under contract to JPL. Peter Smith is the Principal Investigator.
Click below to see the left and right views individually. [figure removed for brevity, see original site] Left [figure removed for brevity, see original site] Right
NASA Technical Reports Server (NTRS)
2004-01-01
This 3-D, microscopic imager mosaic of a target area on a rock called 'Diamond Jenness' was taken after NASA's Mars Exploration Rover Opportunity ground into the surface with its rock abrasion tool for a second time.
Opportunity has bored nearly a dozen holes into the inner walls of 'Endurance Crater.' On sols 177 and 178 (July 23 and July 24, 2004), the rover worked double-duty on Diamond Jenness. Surface debris and the bumpy shape of the rock resulted in a shallow and irregular hole, only about 2 millimeters (0.08 inch) deep. The final depth was not enough to remove all the bumps and leave a neat hole with a smooth floor. This extremely shallow depression was then examined by the rover's alpha particle X-ray spectrometer.
On Sol 178, Opportunity's 'robotic rodent' dined on Diamond Jenness once again, grinding almost an additional 5 millimeters (about 0.2 inch). The rover then applied its Moessbauer spectrometer to the deepened hole. This double dose of Diamond Jenness enabled the science team to examine the rock at varying layers. Results from those grindings are currently being analyzed.
The image mosaic is about 6 centimeters (2.4 inches) across.
Large-scale 3D simulations of ICF and HEDP targets
NASA Astrophysics Data System (ADS)
Marinak, Michael M.
2000-10-01
The radiation hydrodynamics code HYDRA continues to be developed and applied to 3D simulations of a variety of targets for both inertial confinement fusion (ICF) and high energy density physics. Several packages have been added enabling this code to perform ICF target simulations with similar accuracy as two-dimensional codes of long-time historical use. These include a laser ray trace and deposition package, a heavy ion deposition package, implicit Monte Carlo photonics, and non-LTE opacities, derived from XSN or the linearized response matrix approach.(R. More, T. Kato, Phys. Rev. Lett. 81, 814 (1998), S. Libby, F. Graziani, R. More, T. Kato, Proceedings of the 13th International Conference on Laser Interactions and Related Plasma Phenomena, (AIP, New York, 1997).) LTE opacities can also be calculated for arbitrary mixtures online by combining tabular values generated by different opacity codes. Thermonuclear burn, charged particle transport, neutron energy deposition, electron-ion coupling and conduction, and multigroup radiation diffusion packages are also installed. HYDRA can employ ALE hydrodynamics; a number of grid motion algorithms are available. Multi-material flows are resolved using material interface reconstruction. Results from large-scale simulations run on up to 1680 processors, using a combination of massively parallel processing and symmetric multiprocessing, will be described. A large solid angle simulation of Rayleigh-Taylor instability growth in a NIF ignition capsule has resolved simultaneously the full spectrum of the most dangerous modes that grow from surface roughness. Simulations of a NIF hohlraum illuminated with the initial 96 beam configuration have also been performed. The effect of the hohlraum’s 3D intrinsic drive asymmetry on the capsule implosion will be considered. We will also discuss results from a Nova experiment in which a copper sphere is crushed by a planar shock. Several interacting hydrodynamic instabilities, including
3-D Finite Element Heat Transfer
Energy Science and Technology Software Center (ESTSC)
1992-02-01
TOPAZ3D is a three-dimensional implicit finite element computer code for heat transfer analysis. TOPAZ3D can be used to solve for the steady-state or transient temperature field on three-dimensional geometries. Material properties may be temperature-dependent and either isotropic or orthotropic. A variety of time-dependent and temperature-dependent boundary conditions can be specified including temperature, flux, convection, and radiation. By implementing the user subroutine feature, users can model chemical reaction kinetics and allow for any type of functionalmore » representation of boundary conditions and internal heat generation. TOPAZ3D can solve problems of diffuse and specular band radiation in an enclosure coupled with conduction in the material surrounding the enclosure. Additional features include thermal contact resistance across an interface, bulk fluids, phase change, and energy balances.« less
A Note on Multigroup Comparisons Using SAS PROC CALIS
ERIC Educational Resources Information Center
Jones-Farmer, L. Allison; Pitts, Jennifer P.; Rainer, R. Kelly
2008-01-01
Although SAS PROC CALIS is not designed to perform multigroup comparisons, it is believed that SAS can be "tricked" into doing so for groups of equal size. At present, there are no comprehensive examples of the steps involved in performing a multigroup comparison in SAS. The purpose of this article is to illustrate these steps. We demonstrate…
Multigroup calculations using VIM: A user's guide to ISOVIM
Blomquist, R.N.
1992-09-01
Monte Carlo calculations have long been used to benchmark more a mate approximate solution methods for reactor physics problems. The power of VIM (ref 1) lies partly in the detailed geometrical representations incorporating the (generally) curved surfaces of combinatorial geometry, and partly in the fine energy detail of pointwise cross sections which are independent of the neutron spectrum. When differences arise between Monte Carlo and deterministic calculations, the question arises, is the error in the multigroup cross sections, in the treatment of transport effects, or in the mesh-based treatment of space in the deterministic calculation The answers may not be obvious, but may be identified by combining the exact geometry capability of VIM with the multigroup formalism. We can now run VIM in a multigroup mode by producing special VIM Material files which contain point-wise data describing multigroup data with histograms. This forces VIM to solve the multigroup problem with only three small code modifications. P[sub N] scattering is simulated with the usual tabulated angular distributions with 20 equally-sized scattering angle cosine meshes. This document describes the VIM multigroup capability, the procedures for generating multigroup cross sections for VIM, and their use. The multigroup cross section generating code, ISOVIM, is described, and benchmark testing is documented.
NASA Technical Reports Server (NTRS)
2009-01-01
wavelengths. Since the amount of the wavelength shift is related to the speed of motion, one can determine how fast the debris are moving in either direction. Because Cas A is the result of an explosion, the stellar debris is expanding radially outwards from the explosion center. Using simple geometry, the scientists were able to construct a 3-D model using all of this information. A program called 3-D Slicer modified for astronomical use by the Astronomical Medicine Project at Harvard University in Cambridge, Mass. was used to display and manipulate the 3-D model. Commercial software was then used to create the 3-D fly-through.
The blue filaments defining the blast wave were not mapped using the Doppler effect because they emit a different kind of light synchrotron radiation that does not emit light at discrete wavelengths, but rather in a broad continuum. The blue filaments are only a representation of the actual filaments observed at the blast wave.
This visualization shows that there are two main components to this supernova remnant: a spherical component in the outer parts of the remnant and a flattened (disk-like) component in the inner region. The spherical component consists of the outer layer of the star that exploded, probably made of helium and carbon. These layers drove a spherical blast wave into the diffuse gas surrounding the star. The flattened component that astronomers were unable to map into 3-D prior to these Spitzer observations consists of the inner layers of the star. It is made from various heavier elements, not all shown in the visualization, such as oxygen, neon, silicon, sulphur, argon and iron.
High-velocity plumes, or jets, of this material are shooting out from the explosion in the plane of the disk-like component mentioned above. Plumes of silicon appear in the northeast and southwest, while those of iron are seen in the southeast and north. These jets were already known and Doppler velocity measurements have been made for these
NASA Astrophysics Data System (ADS)
Mediavilla, Evencio; Arribas, Santiago; Roth, Martin; Cepa-Nogué, Jordi; Sánchez, Francisco
2011-09-01
Preface; Acknowledgements; 1. Introductory review and technical approaches Martin M. Roth; 2. Observational procedures and data reduction James E. H. Turner; 3. 3D Spectroscopy instrumentation M. A. Bershady; 4. Analysis of 3D data Pierre Ferruit; 5. Science motivation for IFS and galactic studies F. Eisenhauer; 6. Extragalactic studies and future IFS science Luis Colina; 7. Tutorials: how to handle 3D spectroscopy data Sebastian F. Sánchez, Begona García-Lorenzo and Arlette Pécontal-Rousset.
3D Elevation Program—Virtual USA in 3D
Lukas, Vicki; Stoker, J.M.
2016-01-01
The U.S. Geological Survey (USGS) 3D Elevation Program (3DEP) uses a laser system called ‘lidar’ (light detection and ranging) to create a virtual reality map of the Nation that is very accurate. 3D maps have many uses with new uses being discovered all the time.
A microfluidic device for 2D to 3D and 3D to 3D cell navigation
NASA Astrophysics Data System (ADS)
Shamloo, Amir; Amirifar, Leyla
2016-01-01
Microfluidic devices have received wide attention and shown great potential in the field of tissue engineering and regenerative medicine. Investigating cell response to various stimulations is much more accurate and comprehensive with the aid of microfluidic devices. In this study, we introduced a microfluidic device by which the matrix density as a mechanical property and the concentration profile of a biochemical factor as a chemical property could be altered. Our microfluidic device has a cell tank and a cell culture chamber to mimic both 2D to 3D and 3D to 3D migration of three types of cells. Fluid shear stress is negligible on the cells and a stable concentration gradient can be obtained by diffusion. The device was designed by a numerical simulation so that the uniformity of the concentration gradients throughout the cell culture chamber was obtained. Adult neural cells were cultured within this device and they showed different branching and axonal navigation phenotypes within varying nerve growth factor (NGF) concentration profiles. Neural stem cells were also cultured within varying collagen matrix densities while exposed to NGF concentrations and they experienced 3D to 3D collective migration. By generating vascular endothelial growth factor concentration gradients, adult human dermal microvascular endothelial cells also migrated in a 2D to 3D manner and formed a stable lumen within a specific collagen matrix density. It was observed that a minimum absolute concentration and concentration gradient were required to stimulate migration of all types of the cells. This device has the advantage of changing multiple parameters simultaneously and is expected to have wide applicability in cell studies.
NASA Technical Reports Server (NTRS)
2009-01-01
wavelengths. Since the amount of the wavelength shift is related to the speed of motion, one can determine how fast the debris are moving in either direction. Because Cas A is the result of an explosion, the stellar debris is expanding radially outwards from the explosion center. Using simple geometry, the scientists were able to construct a 3-D model using all of this information. A program called 3-D Slicer modified for astronomical use by the Astronomical Medicine Project at Harvard University in Cambridge, Mass. was used to display and manipulate the 3-D model. Commercial software was then used to create the 3-D fly-through.
The blue filaments defining the blast wave were not mapped using the Doppler effect because they emit a different kind of light synchrotron radiation that does not emit light at discrete wavelengths, but rather in a broad continuum. The blue filaments are only a representation of the actual filaments observed at the blast wave.
This visualization shows that there are two main components to this supernova remnant: a spherical component in the outer parts of the remnant and a flattened (disk-like) component in the inner region. The spherical component consists of the outer layer of the star that exploded, probably made of helium and carbon. These layers drove a spherical blast wave into the diffuse gas surrounding the star. The flattened component that astronomers were unable to map into 3-D prior to these Spitzer observations consists of the inner layers of the star. It is made from various heavier elements, not all shown in the visualization, such as oxygen, neon, silicon, sulphur, argon and iron.
High-velocity plumes, or jets, of this material are shooting out from the explosion in the plane of the disk-like component mentioned above. Plumes of silicon appear in the northeast and southwest, while those of iron are seen in the southeast and north. These jets were already known and Doppler velocity measurements have been made for these
Multigroup neutron dose calculations for proton therapy
Kelsey Iv, Charles T; Prinja, Anil K
2009-01-01
We have developed tools for the preparation of coupled multigroup proton/neutron cross section libraries. Our method is to use NJOY to process evaluated nuclear data files for incident particles below 150 MeV and MCNPX to produce data for higher energies. We modified the XSEX3 program of the MCNPX code system to produce Legendre expansions of scattering matrices generated by sampling the physics models that are comparable to the output of the GROUPR routine of NJOY. Our code combines the low and high energy scattering data with user input stopping powers and energy deposition cross sections that we also calculated using MCNPX. Our code also calculates momentum transfer coefficients for the library and optionally applies an energy straggling model to the scattering cross sections and stopping powers. The motivation was initially for deterministic solution of space radiation shielding calculations using Attila, but noting that proton therapy treatment planning may neglect secondary neutron dose assessments because of difficulty and expense, we have also investigated the feasibility of multi group methods for this application. We have shown that multigroup MCNPX solutions for secondary neutron dose compare well with continuous energy solutions and are obtainable with less than half computational cost. This efficiency comparison neglects the cost of preparing the library data, but this becomes negligible when distributed over many multi group calculations. Our deterministic calculations illustrate recognized obstacles that may have to be overcome before discrete ordinates methods can be efficient alternatives for proton therapy neutron dose calculations.
NASA Technical Reports Server (NTRS)
1977-01-01
A market study of a proposed version of a 3-D eyetracker for initial use at NASA's Ames Research Center was made. The commercialization potential of a simplified, less expensive 3-D eyetracker was ascertained. Primary focus on present and potential users of eyetrackers, as well as present and potential manufacturers has provided an effective means of analyzing the prospects for commercialization.
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...
[3-D ultrasound in gastroenterology].
Zoller, W G; Liess, H
1994-06-01
Three-dimensional (3D) sonography represents a development of noninvasive diagnostic imaging by real-time two-dimensional (2D) sonography. The use of transparent rotating scans, comparable to a block of glass, generates a 3D effect. The objective of the present study was to optimate 3D presentation of abdominal findings. Additional investigations were made with a new volumetric program to determine the volume of selected findings of the liver. The results were compared with the estimated volumes of 2D sonography and 2D computer tomography (CT). For the processing of 3D images, typical parameter constellations were found for the different findings, which facilitated processing of 3D images. In more than 75% of the cases examined we found an optimal 3D presentation of sonographic findings with respect to the evaluation criteria developed by us for the 3D imaging of processed data. Great differences were found for the estimated volumes of the findings of the liver concerning the three different techniques applied. 3D ultrasound represents a valuable method to judge morphological appearance in abdominal findings. The possibility of volumetric measurements enlarges its potential diagnostic significance. Further clinical investigations are necessary to find out if definite differentiation between benign and malign findings is possible. PMID:7919882
2013-10-30
This video provides an overview of the Sandia National Laboratories developed 3-D World Model Building capability that provides users with an immersive, texture rich 3-D model of their environment in minutes using a laptop and color and depth camera.
None
2014-02-26
This video provides an overview of the Sandia National Laboratories developed 3-D World Model Building capability that provides users with an immersive, texture rich 3-D model of their environment in minutes using a laptop and color and depth camera.
Energy Science and Technology Software Center (ESTSC)
2013-10-01
Earth3D is a computer code designed to allow fast calculation of seismic rays and travel times through a 3D model of the Earth. LLNL is using this for earthquake location and global tomography efforts and such codes are of great interest to the Earth Science community.
NASA Astrophysics Data System (ADS)
Walsh, J. R.
2004-02-01
The Euro3D RTN is an EU funded Research Training Network to foster the exploitation of 3D spectroscopy in Europe. 3D spectroscopy is a general term for spectroscopy of an area of the sky and derives its name from its two spatial + one spectral dimensions. There are an increasing number of instruments which use integral field devices to achieve spectroscopy of an area of the sky, either using lens arrays, optical fibres or image slicers, to pack spectra of multiple pixels on the sky (``spaxels'') onto a 2D detector. On account of the large volume of data and the special methods required to reduce and analyse 3D data, there are only a few centres of expertise and these are mostly involved with instrument developments. There is a perceived lack of expertise in 3D spectroscopy spread though the astronomical community and its use in the armoury of the observational astronomer is viewed as being highly specialised. For precisely this reason the Euro3D RTN was proposed to train young researchers in this area and develop user tools to widen the experience with this particular type of data in Europe. The Euro3D RTN is coordinated by Martin M. Roth (Astrophysikalisches Institut Potsdam) and has been running since July 2002. The first Euro3D science conference was held in Cambridge, UK from 22 to 23 May 2003. The main emphasis of the conference was, in keeping with the RTN, to expose the work of the young post-docs who are funded by the RTN. In addition the team members from the eleven European institutes involved in Euro3D also presented instrumental and observational developments. The conference was organized by Andy Bunker and held at the Institute of Astronomy. There were over thirty participants and 26 talks covered the whole range of application of 3D techniques. The science ranged from Galactic planetary nebulae and globular clusters to kinematics of nearby galaxies out to objects at high redshift. Several talks were devoted to reporting recent observations with newly
Dawood, A; Marti Marti, B; Sauret-Jackson, V; Darwood, A
2015-12-01
3D printing has been hailed as a disruptive technology which will change manufacturing. Used in aerospace, defence, art and design, 3D printing is becoming a subject of great interest in surgery. The technology has a particular resonance with dentistry, and with advances in 3D imaging and modelling technologies such as cone beam computed tomography and intraoral scanning, and with the relatively long history of the use of CAD CAM technologies in dentistry, it will become of increasing importance. Uses of 3D printing include the production of drill guides for dental implants, the production of physical models for prosthodontics, orthodontics and surgery, the manufacture of dental, craniomaxillofacial and orthopaedic implants, and the fabrication of copings and frameworks for implant and dental restorations. This paper reviews the types of 3D printing technologies available and their various applications in dentistry and in maxillofacial surgery. PMID:26657435
NASA Technical Reports Server (NTRS)
Walatka, Pamela P.; Buning, Pieter G.; Pierce, Larry; Elson, Patricia A.
1990-01-01
PLOT3D is a computer graphics program designed to visualize the grids and solutions of computational fluid dynamics. Seventy-four functions are available. Versions are available for many systems. PLOT3D can handle multiple grids with a million or more grid points, and can produce varieties of model renderings, such as wireframe or flat shaded. Output from PLOT3D can be used in animation programs. The first part of this manual is a tutorial that takes the reader, keystroke by keystroke, through a PLOT3D session. The second part of the manual contains reference chapters, including the helpfile, data file formats, advice on changing PLOT3D, and sample command files.
3D reconstruction of tensors and vectors
Defrise, Michel; Gullberg, Grant T.
2005-02-17
Here we have developed formulations for the reconstruction of 3D tensor fields from planar (Radon) and line-integral (X-ray) projections of 3D vector and tensor fields. Much of the motivation for this work is the potential application of MRI to perform diffusion tensor tomography. The goal is to develop a theory for the reconstruction of both Radon planar and X-ray or line-integral projections because of the flexibility of MRI to obtain both of these type of projections in 3D. The development presented here for the linear tensor tomography problem provides insight into the structure of the nonlinear MRI diffusion tensor inverse problem. A particular application of tensor imaging in MRI is the potential application of cardiac diffusion tensor tomography for determining in vivo cardiac fiber structure. One difficulty in the cardiac application is the motion of the heart. This presents a need for developing future theory for tensor tomography in a motion field. This means developing a better understanding of the MRI signal for diffusion processes in a deforming media. The techniques developed may allow the application of MRI tensor tomography for the study of structure of fiber tracts in the brain, atherosclerotic plaque, and spine in addition to fiber structure in the heart. However, the relations presented are also applicable to other fields in medical imaging such as diffraction tomography using ultrasound. The mathematics presented can also be extended to exponential Radon transform of tensor fields and to other geometric acquisitions such as cone beam tomography of tensor fields.
PLOT3D/AMES, APOLLO UNIX VERSION USING GMR3D (WITH TURB3D)
NASA Technical Reports Server (NTRS)
Buning, P.
1994-01-01
PLOT3D is an interactive graphics program designed to help scientists visualize computational fluid dynamics (CFD) grids and solutions. Today, supercomputers and CFD algorithms can provide scientists with simulations of such highly complex phenomena that obtaining an understanding of the simulations has become a major problem. Tools which help the scientist visualize the simulations can be of tremendous aid. PLOT3D/AMES offers more functions and features, and has been adapted for more types of computers than any other CFD graphics program. Version 3.6b+ is supported for five computers and graphic libraries. Using PLOT3D, CFD physicists can view their computational models from any angle, observing the physics of problems and the quality of solutions. As an aid in designing aircraft, for example, PLOT3D's interactive computer graphics can show vortices, temperature, reverse flow, pressure, and dozens of other characteristics of air flow during flight. As critical areas become obvious, they can easily be studied more closely using a finer grid. PLOT3D is part of a computational fluid dynamics software cycle. First, a program such as 3DGRAPE (ARC-12620) helps the scientist generate computational grids to model an object and its surrounding space. Once the grids have been designed and parameters such as the angle of attack, Mach number, and Reynolds number have been specified, a "flow-solver" program such as INS3D (ARC-11794 or COS-10019) solves the system of equations governing fluid flow, usually on a supercomputer. Grids sometimes have as many as two million points, and the "flow-solver" produces a solution file which contains density, x- y- and z-momentum, and stagnation energy for each grid point. With such a solution file and a grid file containing up to 50 grids as input, PLOT3D can calculate and graphically display any one of 74 functions, including shock waves, surface pressure, velocity vectors, and particle traces. PLOT3D's 74 functions are organized into
PLOT3D/AMES, APOLLO UNIX VERSION USING GMR3D (WITHOUT TURB3D)
NASA Technical Reports Server (NTRS)
Buning, P.
1994-01-01
PLOT3D is an interactive graphics program designed to help scientists visualize computational fluid dynamics (CFD) grids and solutions. Today, supercomputers and CFD algorithms can provide scientists with simulations of such highly complex phenomena that obtaining an understanding of the simulations has become a major problem. Tools which help the scientist visualize the simulations can be of tremendous aid. PLOT3D/AMES offers more functions and features, and has been adapted for more types of computers than any other CFD graphics program. Version 3.6b+ is supported for five computers and graphic libraries. Using PLOT3D, CFD physicists can view their computational models from any angle, observing the physics of problems and the quality of solutions. As an aid in designing aircraft, for example, PLOT3D's interactive computer graphics can show vortices, temperature, reverse flow, pressure, and dozens of other characteristics of air flow during flight. As critical areas become obvious, they can easily be studied more closely using a finer grid. PLOT3D is part of a computational fluid dynamics software cycle. First, a program such as 3DGRAPE (ARC-12620) helps the scientist generate computational grids to model an object and its surrounding space. Once the grids have been designed and parameters such as the angle of attack, Mach number, and Reynolds number have been specified, a "flow-solver" program such as INS3D (ARC-11794 or COS-10019) solves the system of equations governing fluid flow, usually on a supercomputer. Grids sometimes have as many as two million points, and the "flow-solver" produces a solution file which contains density, x- y- and z-momentum, and stagnation energy for each grid point. With such a solution file and a grid file containing up to 50 grids as input, PLOT3D can calculate and graphically display any one of 74 functions, including shock waves, surface pressure, velocity vectors, and particle traces. PLOT3D's 74 functions are organized into
Stanton, M M; Samitier, J; Sánchez, S
2015-08-01
Three-dimensional (3D) bioprinting has recently emerged as an extension of 3D material printing, by using biocompatible or cellular components to build structures in an additive, layer-by-layer methodology for encapsulation and culture of cells. These 3D systems allow for cell culture in a suspension for formation of highly organized tissue or controlled spatial orientation of cell environments. The in vitro 3D cellular environments simulate the complexity of an in vivo environment and natural extracellular matrices (ECM). This paper will focus on bioprinting utilizing hydrogels as 3D scaffolds. Hydrogels are advantageous for cell culture as they are highly permeable to cell culture media, nutrients, and waste products generated during metabolic cell processes. They have the ability to be fabricated in customized shapes with various material properties with dimensions at the micron scale. 3D hydrogels are a reliable method for biocompatible 3D printing and have applications in tissue engineering, drug screening, and organ on a chip models. PMID:26066320
Unassisted 3D camera calibration
NASA Astrophysics Data System (ADS)
Atanassov, Kalin; Ramachandra, Vikas; Nash, James; Goma, Sergio R.
2012-03-01
With the rapid growth of 3D technology, 3D image capture has become a critical part of the 3D feature set on mobile phones. 3D image quality is affected by the scene geometry as well as on-the-device processing. An automatic 3D system usually assumes known camera poses accomplished by factory calibration using a special chart. In real life settings, pose parameters estimated by factory calibration can be negatively impacted by movements of the lens barrel due to shaking, focusing, or camera drop. If any of these factors displaces the optical axes of either or both cameras, vertical disparity might exceed the maximum tolerable margin and the 3D user may experience eye strain or headaches. To make 3D capture more practical, one needs to consider unassisted (on arbitrary scenes) calibration. In this paper, we propose an algorithm that relies on detection and matching of keypoints between left and right images. Frames containing erroneous matches, along with frames with insufficiently rich keypoint constellations, are detected and discarded. Roll, pitch yaw , and scale differences between left and right frames are then estimated. The algorithm performance is evaluated in terms of the remaining vertical disparity as compared to the maximum tolerable vertical disparity.
Arena3D: visualization of biological networks in 3D
Pavlopoulos, Georgios A; O'Donoghue, Seán I; Satagopam, Venkata P; Soldatos, Theodoros G; Pafilis, Evangelos; Schneider, Reinhard
2008-01-01
Background Complexity is a key problem when visualizing biological networks; as the number of entities increases, most graphical views become incomprehensible. Our goal is to enable many thousands of entities to be visualized meaningfully and with high performance. Results We present a new visualization tool, Arena3D, which introduces a new concept of staggered layers in 3D space. Related data – such as proteins, chemicals, or pathways – can be grouped onto separate layers and arranged via layout algorithms, such as Fruchterman-Reingold, distance geometry, and a novel hierarchical layout. Data on a layer can be clustered via k-means, affinity propagation, Markov clustering, neighbor joining, tree clustering, or UPGMA ('unweighted pair-group method with arithmetic mean'). A simple input format defines the name and URL for each node, and defines connections or similarity scores between pairs of nodes. The use of Arena3D is illustrated with datasets related to Huntington's disease. Conclusion Arena3D is a user friendly visualization tool that is able to visualize biological or any other network in 3D space. It is free for academic use and runs on any platform. It can be downloaded or lunched directly from . Java3D library and Java 1.5 need to be pre-installed for the software to run. PMID:19040715
NASA Astrophysics Data System (ADS)
Otis, Collin; Ferrero, Pietro; Candler, Graham; Givi, Peyman
2013-11-01
The scalar filtered mass density function (SFMDF) methodology is implemented into the computer code US3D. This is an unstructured Eulerian finite volume hydrodynamic solver and has proven very effective for simulation of compressible turbulent flows. The resulting SFMDF-US3D code is employed for large eddy simulation (LES) on unstructured meshes. Simulations are conducted of subsonic and supersonic flows under non-reacting and reacting conditions. The consistency and the accuracy of the simulated results are assessed along with appraisal of the overall performance of the methodology. The SFMDF-US3D is now capable of simulating high speed flows in complex configurations.
Multigroup Free-atom Doppler-broadening Approximation. Theory
Gray, Mark Girard
2015-11-06
Multigroup cross sections at a one target temperature can be Doppler-broadened to multigroup cross sections at a higher target temperature by matrix multiplication if the group structure suf- ficiently resolves the original temperature continuous energy cross section. Matrix elements are the higher temperature group weighted averages of the integral over the lower temperature group boundaries of the free-atom Doppler-broadening kernel. The results match theory for constant and 1/v multigroup cross sections at 618 lanl group structure resolution.
Solves Multigroup Diffusion Equations in One-Dimensional Systems.
Energy Science and Technology Software Center (ESTSC)
1984-11-01
Version 00 UNIMUG3 determines the multiplication factor, fluxes, and adjoint fluxes. It permits variation of the volume of any one of the regions for determining the critical state and performs the adjoint flux calculation and critical size search.
3D optical tomography in the presence of void regions
NASA Astrophysics Data System (ADS)
Riley, J.; Dehghani, Hamid; Schweiger, Martin; Arridge, Simon R.; Ripoll, Jorge; Nieto-Vesperinas, Manuel
2000-12-01
We present an investigation of the effect of a 3D non-scattering gap region on image reconstruction in diffuse optical tomography. The void gap is modelled by the Radiosity-Diffusion method and the inverse problem is solved using the adjoint field method. The case of a sphere with concentric spherical gap is used as an example.
3D optical tomography in the presence of void regions.
Riley, J; Dehghani, H; Schweiger, M; Arridge, S; Ripoll, J; Nieto-Vesperinas, M
2000-12-18
We present an investigation of the effect of a 3D non-scattering gap region on image reconstruction in diffuse optical tomography. The void gap is modelled by the Radiosity-Diffusion method and the inverse problem is solved using the adjoint field method. The case of a sphere with concentric spherical gap is used as an example. PMID:19407898
NASA Astrophysics Data System (ADS)
Yang, Xu; Zhang, Yong; Yang, Chenghua; Xu, Lu; Wang, Qiang; Zhao, Yuan
2016-06-01
Conventional three dimensional (3D) ghost imaging measures range of target based on pulse fight time measurement method. Due to the limit of data acquisition system sampling rate, range resolution of the conventional 3D ghost imaging is usually low. In order to take off the effect of sampling rate to range resolution of 3D ghost imaging, a heterodyne 3D ghost imaging (HGI) system is presented in this study. The source of HGI is a continuous wave laser instead of pulse laser. Temporal correlation and spatial correlation of light are both utilized to obtain the range image of target. Through theory analysis and numerical simulations, it is demonstrated that HGI can obtain high range resolution image with low sampling rate.
Chilcoat, S.R. Hildebrand, S.T.
1995-12-31
Travel time computation in inhomogeneous media is essential for pre-stack Kirchhoff imaging in areas such as the sub-salt province in the Gulf of Mexico. The 2D algorithm published by Vinje, et al, has been extended to 3D to compute wavefronts in complicated inhomogeneous media. The 3D wavefront construction algorithm provides many advantages over conventional ray tracing and other methods of computing travel times in 3D. The algorithm dynamically maintains a reasonably consistent ray density without making a priori guesses at the number of rays to shoot. The determination of caustics in 3D is a straight forward geometric procedure. The wavefront algorithm also enables the computation of multi-valued travel time surfaces.
Combinatorial 3D Mechanical Metamaterials
NASA Astrophysics Data System (ADS)
Coulais, Corentin; Teomy, Eial; de Reus, Koen; Shokef, Yair; van Hecke, Martin
2015-03-01
We present a class of elastic structures which exhibit 3D-folding motion. Our structures consist of cubic lattices of anisotropic unit cells that can be tiled in a complex combinatorial fashion. We design and 3d-print this complex ordered mechanism, in which we combine elastic hinges and defects to tailor the mechanics of the material. Finally, we use this large design space to encode smart functionalities such as surface patterning and multistability.
A Bayesian Approach for Multigroup Nonlinear Factor Analysis.
ERIC Educational Resources Information Center
Song, Xin-Yuan; Lee, Sik-Yum
2002-01-01
Developed a Bayesian approach for a general multigroup nonlinear factor analysis model that simultaneously obtains joint Bayesian estimates of the factor scores and the structural parameters subjected to some constraints across different groups. (SLD)
Experimental Study of Electrothermal 3D Mixing using 3D microPIV
NASA Astrophysics Data System (ADS)
Kauffmann, Paul; Loire, Sophie; Meinhart, Carl; Mezic, Igor
2012-11-01
Mixing is a keystep which can greatly accelerate bio-reactions. For thirty years, dynamical system theory has predicted that chaotic mixing must involve at least 3 dimensions (either time dependent 2D flows or 3D flows). So far, 3D embedded chaotic mixing has been scarcely studied at microscale. In that regard, electrokinetics has emerged as an efficient embedded actuation to drive microflows. Physiological mediums can be driven by electrothermal flows generated by the interaction of an electric field with conductivity and permittivity gradients induced by Joule heating We present original electrothermal time dependant 3D (3D+1) mixing in microwells. The key point of our chaotic mixer is to generate overlapping asymmetric vortices, which switch periodically. When the two vortex configurations blink, flows stretch and fold, thereby generating chaotic advection. Each flow configuration is characterized by an original 3D PIV (3 Components / 3 Dimensions) based on the decomposition of the flows by Proper Orthogonal Decomposition. Velocity field distribution are then compared to COMSOL simulation and discussed. Mixing efficiency of low diffusive particles is studied using the mix-variance coefficient and shows a dramatic increase of mixing efficiency compared to steady flow.
YouDash3D: exploring stereoscopic 3D gaming for 3D movie theaters
NASA Astrophysics Data System (ADS)
Schild, Jonas; Seele, Sven; Masuch, Maic
2012-03-01
Along with the success of the digitally revived stereoscopic cinema, events beyond 3D movies become attractive for movie theater operators, i.e. interactive 3D games. In this paper, we present a case that explores possible challenges and solutions for interactive 3D games to be played by a movie theater audience. We analyze the setting and showcase current issues related to lighting and interaction. Our second focus is to provide gameplay mechanics that make special use of stereoscopy, especially depth-based game design. Based on these results, we present YouDash3D, a game prototype that explores public stereoscopic gameplay in a reduced kiosk setup. It features live 3D HD video stream of a professional stereo camera rig rendered in a real-time game scene. We use the effect to place the stereoscopic effigies of players into the digital game. The game showcases how stereoscopic vision can provide for a novel depth-based game mechanic. Projected trigger zones and distributed clusters of the audience video allow for easy adaptation to larger audiences and 3D movie theater gaming.
Remote 3D Medical Consultation
NASA Astrophysics Data System (ADS)
Welch, Greg; Sonnenwald, Diane H.; Fuchs, Henry; Cairns, Bruce; Mayer-Patel, Ketan; Yang, Ruigang; State, Andrei; Towles, Herman; Ilie, Adrian; Krishnan, Srinivas; Söderholm, Hanna M.
Two-dimensional (2D) video-based telemedical consultation has been explored widely in the past 15-20 years. Two issues that seem to arise in most relevant case studies are the difficulty associated with obtaining the desired 2D camera views, and poor depth perception. To address these problems we are exploring the use of a small array of cameras to synthesize a spatially continuous range of dynamic three-dimensional (3D) views of a remote environment and events. The 3D views can be sent across wired or wireless networks to remote viewers with fixed displays or mobile devices such as a personal digital assistant (PDA). The viewpoints could be specified manually or automatically via user head or PDA tracking, giving the remote viewer virtual head- or hand-slaved (PDA-based) remote cameras for mono or stereo viewing. We call this idea remote 3D medical consultation (3DMC). In this article we motivate and explain the vision for 3D medical consultation; we describe the relevant computer vision/graphics, display, and networking research; we present a proof-of-concept prototype system; and we present some early experimental results supporting the general hypothesis that 3D remote medical consultation could offer benefits over conventional 2D televideo.
NASA Technical Reports Server (NTRS)
2002-01-01
In 1999, Genex submitted a proposal to Stennis Space Center for a volumetric 3-D display technique that would provide multiple users with a 360-degree perspective to simultaneously view and analyze 3-D data. The futuristic capabilities of the VolumeViewer(R) have offered tremendous benefits to commercial users in the fields of medicine and surgery, air traffic control, pilot training and education, computer-aided design/computer-aided manufacturing, and military/battlefield management. The technology has also helped NASA to better analyze and assess the various data collected by its satellite and spacecraft sensors. Genex capitalized on its success with Stennis by introducing two separate products to the commercial market that incorporate key elements of the 3-D display technology designed under an SBIR contract. The company Rainbow 3D(R) imaging camera is a novel, three-dimensional surface profile measurement system that can obtain a full-frame 3-D image in less than 1 second. The third product is the 360-degree OmniEye(R) video system. Ideal for intrusion detection, surveillance, and situation management, this unique camera system offers a continuous, panoramic view of a scene in real time.
Au, Anthony K; Huynh, Wilson; Horowitz, Lisa F; Folch, Albert
2016-03-14
The advent of soft lithography allowed for an unprecedented expansion in the field of microfluidics. However, the vast majority of PDMS microfluidic devices are still made with extensive manual labor, are tethered to bulky control systems, and have cumbersome user interfaces, which all render commercialization difficult. On the other hand, 3D printing has begun to embrace the range of sizes and materials that appeal to the developers of microfluidic devices. Prior to fabrication, a design is digitally built as a detailed 3D CAD file. The design can be assembled in modules by remotely collaborating teams, and its mechanical and fluidic behavior can be simulated using finite-element modeling. As structures are created by adding materials without the need for etching or dissolution, processing is environmentally friendly and economically efficient. We predict that in the next few years, 3D printing will replace most PDMS and plastic molding techniques in academia. PMID:26854878
3D Computations and Experiments
Couch, R; Faux, D; Goto, D; Nikkel, D
2004-04-05
This project consists of two activities. Task A, Simulations and Measurements, combines all the material model development and associated numerical work with the materials-oriented experimental activities. The goal of this effort is to provide an improved understanding of dynamic material properties and to provide accurate numerical representations of those properties for use in analysis codes. Task B, ALE3D Development, involves general development activities in the ALE3D code with the focus of improving simulation capabilities for problems of mutual interest to DoD and DOE. Emphasis is on problems involving multi-phase flow, blast loading of structures and system safety/vulnerability studies.
3D Computations and Experiments
Couch, R; Faux, D; Goto, D; Nikkel, D
2003-05-12
This project is in its first full year after the combining of two previously funded projects: ''3D Code Development'' and ''Dynamic Material Properties''. The motivation behind this move was to emphasize and strengthen the ties between the experimental work and the computational model development in the materials area. The next year's activities will indicate the merging of the two efforts. The current activity is structured in two tasks. Task A, ''Simulations and Measurements'', combines all the material model development and associated numerical work with the materials-oriented experimental activities. Task B, ''ALE3D Development'', is a continuation of the non-materials related activities from the previous project.
Energy Science and Technology Software Center (ESTSC)
2007-07-20
This software distribution contains MATLAB and C++ code to enable identity verification using 3D images that may or may not contain a texture component. The code is organized to support system performance testing and system capability demonstration through the proper configuration of the available user interface. Using specific algorithm parameters the face recognition system has been demonstrated to achieve a 96.6% verification rate (Pd) at 0.001 false alarm rate. The system computes robust facial featuresmore » of a 3D normalized face using Principal Component Analysis (PCA) and Fisher Linear Discriminant Analysis (FLDA). A 3D normalized face is obtained by alighning each face, represented by a set of XYZ coordinated, to a scaled reference face using the Iterative Closest Point (ICP) algorithm. The scaled reference face is then deformed to the input face using an iterative framework with parameters that control the deformed surface regulation an rate of deformation. A variety of options are available to control the information that is encoded by the PCA. Such options include the XYZ coordinates, the difference of each XYZ coordinates from the reference, the Z coordinate, the intensity/texture values, etc. In addition to PCA/FLDA feature projection this software supports feature matching to obtain similarity matrices for performance analysis. In addition, this software supports visualization of the STL, MRD, 2D normalized, and PCA synthetic representations in a 3D environment.« less
3D Printing: Exploring Capabilities
ERIC Educational Resources Information Center
Samuels, Kyle; Flowers, Jim
2015-01-01
As 3D printers become more affordable, schools are using them in increasing numbers. They fit well with the emphasis on product design in technology and engineering education, allowing students to create high-fidelity physical models to see and test different iterations in their product designs. They may also help students to "think in three…
Russ, Trina; Koch, Mark; Koudelka, Melissa; Peters, Ralph; Little, Charles; Boehnen, Chris; Peters, Tanya
2007-07-20
This software distribution contains MATLAB and C++ code to enable identity verification using 3D images that may or may not contain a texture component. The code is organized to support system performance testing and system capability demonstration through the proper configuration of the available user interface. Using specific algorithm parameters the face recognition system has been demonstrated to achieve a 96.6% verification rate (Pd) at 0.001 false alarm rate. The system computes robust facial features of a 3D normalized face using Principal Component Analysis (PCA) and Fisher Linear Discriminant Analysis (FLDA). A 3D normalized face is obtained by alighning each face, represented by a set of XYZ coordinated, to a scaled reference face using the Iterative Closest Point (ICP) algorithm. The scaled reference face is then deformed to the input face using an iterative framework with parameters that control the deformed surface regulation an rate of deformation. A variety of options are available to control the information that is encoded by the PCA. Such options include the XYZ coordinates, the difference of each XYZ coordinates from the reference, the Z coordinate, the intensity/texture values, etc. In addition to PCA/FLDA feature projection this software supports feature matching to obtain similarity matrices for performance analysis. In addition, this software supports visualization of the STL, MRD, 2D normalized, and PCA synthetic representations in a 3D environment.
ERIC Educational Resources Information Center
Manos, Harry
2016-01-01
Visual aids are important to student learning, and they help make the teacher's job easier. Keeping with the "TPT" theme of "The Art, Craft, and Science of Physics Teaching," the purpose of this article is to show how teachers, lacking equipment and funds, can construct a durable 3-D model reference frame and a model gravity…
Consistent Multigroup Theory Enabling Accurate Course-Group Simulation of Gen IV Reactors
Rahnema, Farzad; Haghighat, Alireza; Ougouag, Abderrafi
2013-11-29
The objective of this proposal is the development of a consistent multi-group theory that accurately accounts for the energy-angle coupling associated with collapsed-group cross sections. This will allow for coarse-group transport and diffusion theory calculations that exhibit continuous energy accuracy and implicitly treat cross- section resonances. This is of particular importance when considering the highly heterogeneous and optically thin reactor designs within the Next Generation Nuclear Plant (NGNP) framework. In such reactors, ignoring the influence of anisotropy in the angular flux on the collapsed cross section, especially at the interface between core and reflector near which control rods are located, results in inaccurate estimates of the rod worth, a serious safety concern. The scope of this project will include the development and verification of a new multi-group theory enabling high-fidelity transport and diffusion calculations in coarse groups, as well as a methodology for the implementation of this method in existing codes. This will allow for a higher accuracy solution of reactor problems while using fewer groups and will reduce the computational expense. The proposed research represents a fundamental advancement in the understanding and improvement of multi- group theory for reactor analysis.
TACO3D. 3-D Finite Element Heat Transfer Code
Mason, W.E.
1992-03-04
TACO3D is a three-dimensional, finite-element program for heat transfer analysis. An extension of the two-dimensional TACO program, it can perform linear and nonlinear analyses and can be used to solve either transient or steady-state problems. The program accepts time-dependent or temperature-dependent material properties, and materials may be isotropic or orthotropic. A variety of time-dependent and temperature-dependent boundary conditions and loadings are available including temperature, flux, convection, and radiation boundary conditions and internal heat generation. Additional specialized features treat enclosure radiation, bulk nodes, and master/slave internal surface conditions (e.g., contact resistance). Data input via a free-field format is provided. A user subprogram feature allows for any type of functional representation of any independent variable. A profile (bandwidth) minimization option is available. The code is limited to implicit time integration for transient solutions. TACO3D has no general mesh generation capability. Rows of evenly-spaced nodes and rows of sequential elements may be generated, but the program relies on separate mesh generators for complex zoning. TACO3D does not have the ability to calculate view factors internally. Graphical representation of data in the form of time history and spatial plots is provided through links to the POSTACO and GRAPE postprocessor codes.
Embedding objects during 3D printing to add new functionalities.
Yuen, Po Ki
2016-07-01
A novel method for integrating and embedding objects to add new functionalities during 3D printing based on fused deposition modeling (FDM) (also known as fused filament fabrication or molten polymer deposition) is presented. Unlike typical 3D printing, FDM-based 3D printing could allow objects to be integrated and embedded during 3D printing and the FDM-based 3D printed devices do not typically require any post-processing and finishing. Thus, various fluidic devices with integrated glass cover slips or polystyrene films with and without an embedded porous membrane, and optical devices with embedded Corning(®) Fibrance™ Light-Diffusing Fiber were 3D printed to demonstrate the versatility of the FDM-based 3D printing and embedding method. Fluid perfusion flow experiments with a blue colored food dye solution were used to visually confirm fluid flow and/or fluid perfusion through the embedded porous membrane in the 3D printed fluidic devices. Similar to typical 3D printed devices, FDM-based 3D printed devices are translucent at best unless post-polishing is performed and optical transparency is highly desirable in any fluidic devices; integrated glass cover slips or polystyrene films would provide a perfect optical transparent window for observation and visualization. In addition, they also provide a compatible flat smooth surface for biological or biomolecular applications. The 3D printed fluidic devices with an embedded porous membrane are applicable to biological or chemical applications such as continuous perfusion cell culture or biocatalytic synthesis but without the need for any post-device assembly and finishing. The 3D printed devices with embedded Corning(®) Fibrance™ Light-Diffusing Fiber would have applications in display, illumination, or optical applications. Furthermore, the FDM-based 3D printing and embedding method could also be utilized to print casting molds with an integrated glass bottom for polydimethylsiloxane (PDMS) device replication
3D light scanning macrography.
Huber, D; Keller, M; Robert, D
2001-08-01
The technique of 3D light scanning macrography permits the non-invasive surface scanning of small specimens at magnifications up to 200x. Obviating both the problem of limited depth of field inherent to conventional close-up macrophotography and the metallic coating required by scanning electron microscopy, 3D light scanning macrography provides three-dimensional digital images of intact specimens without the loss of colour, texture and transparency information. This newly developed technique offers a versatile, portable and cost-efficient method for the non-invasive digital and photographic documentation of small objects. Computer controlled device operation and digital image acquisition facilitate fast and accurate quantitative morphometric investigations, and the technique offers a broad field of research and educational applications in biological, medical and materials sciences. PMID:11489078
3-D Relativistic MHD Simulations
NASA Astrophysics Data System (ADS)
Nishikawa, K.-I.; Frank, J.; Koide, S.; Sakai, J.-I.; Christodoulou, D. M.; Sol, H.; Mutel, R. L.
1998-12-01
We present 3-D numerical simulations of moderately hot, supersonic jets propagating initially along or obliquely to the field lines of a denser magnetized background medium with Lorentz factors of W = 4.56 and evolving in a four-dimensional spacetime. The new results are understood as follows: Relativistic simulations have consistently shown that these jets are effectively heavy and so they do not suffer substantial momentum losses and are not decelerated as efficiently as their nonrelativistic counterparts. In addition, the ambient magnetic field, however strong, can be pushed aside with relative ease by the beam, provided that the degrees of freedom associated with all three spatial dimensions are followed self-consistently in the simulations. This effect is analogous to pushing Japanese ``noren'' or vertical Venetian blinds out of the way while the slats are allowed to bend in 3-D space rather than as a 2-D slab structure.
Forensic 3D Scene Reconstruction
LITTLE,CHARLES Q.; PETERS,RALPH R.; RIGDON,J. BRIAN; SMALL,DANIEL E.
1999-10-12
Traditionally law enforcement agencies have relied on basic measurement and imaging tools, such as tape measures and cameras, in recording a crime scene. A disadvantage of these methods is that they are slow and cumbersome. The development of a portable system that can rapidly record a crime scene with current camera imaging, 3D geometric surface maps, and contribute quantitative measurements such as accurate relative positioning of crime scene objects, would be an asset to law enforcement agents in collecting and recording significant forensic data. The purpose of this project is to develop a feasible prototype of a fast, accurate, 3D measurement and imaging system that would support law enforcement agents to quickly document and accurately record a crime scene.
Forensic 3D scene reconstruction
NASA Astrophysics Data System (ADS)
Little, Charles Q.; Small, Daniel E.; Peters, Ralph R.; Rigdon, J. B.
2000-05-01
Traditionally law enforcement agencies have relied on basic measurement and imaging tools, such as tape measures and cameras, in recording a crime scene. A disadvantage of these methods is that they are slow and cumbersome. The development of a portable system that can rapidly record a crime scene with current camera imaging, 3D geometric surface maps, and contribute quantitative measurements such as accurate relative positioning of crime scene objects, would be an asset to law enforcement agents in collecting and recording significant forensic data. The purpose of this project is to develop a fieldable prototype of a fast, accurate, 3D measurement and imaging system that would support law enforcement agents to quickly document and accurately record a crime scene.
NASA Astrophysics Data System (ADS)
Song, Yuanhe; Zhao, Hong; Chen, Wenyi; Tan, Yushan
1997-12-01
A new method of 360 degree turning 3D shape measurement in which light sectioning and phase shifting techniques are both used is presented in this paper. A sine light field is applied in the projected light stripe, meanwhile phase shifting technique is used to calculate phases of the light slit. Thereafter wrapped phase distribution of the slit is formed and the unwrapping process is made by means of the height information based on the light sectioning method. Therefore phase measuring results with better precision can be obtained. At last the target 3D shape data can be produced according to geometric relationships between phases and the object heights. The principles of this method are discussed in detail and experimental results are shown in this paper.
Optoplasmonics: hybridization in 3D
NASA Astrophysics Data System (ADS)
Rosa, L.; Gervinskas, G.; Žukauskas, A.; Malinauskas, M.; Brasselet, E.; Juodkazis, S.
2013-12-01
Femtosecond laser fabrication has been used to make hybrid refractive and di ractive micro-optical elements in photo-polymer SZ2080. For applications in micro- uidics, axicon lenses were fabricated (both single and arrays), for generation of light intensity patterns extending through the entire depth of a typically tens-of-micrometers deep channel. Further hybridisation of an axicon with a plasmonic slot is fabricated and demonstrated nu- merically. Spiralling chiral grooves were inscribed into a 100-nm-thick gold coating sputtered over polymerized micro-axicon lenses, using a focused ion beam. This demonstrates possibility of hybridisation between optical and plasmonic 3D micro-optical elements. Numerical modelling of optical performance by 3D-FDTD method is presented.
NASA Technical Reports Server (NTRS)
Pizarro, Yaritzmar Rosario; Schuler, Jason M.; Lippitt, Thomas C.
2013-01-01
Dexterous robotic hands are changing the way robots and humans interact and use common tools. Unfortunately, the complexity of the joints and actuations drive up the manufacturing cost. Some cutting edge and commercially available rapid prototyping machines now have the ability to print multiple materials and even combine these materials in the same job. A 3D model of a robotic hand was designed using Creo Parametric 2.0. Combining "hard" and "soft" materials, the model was printed on the Object Connex350 3D printer with the purpose of resembling as much as possible the human appearance and mobility of a real hand while needing no assembly. After printing the prototype, strings where installed as actuators to test mobility. Based on printing materials, the manufacturing cost of the hand was $167, significantly lower than other robotic hands without the actuators since they have more complex assembly processes.
Belenkov, E. A. Ali-Pasha, V. A.
2011-01-15
The structure of clusters of some new carbon 3D-graphite phases have been calculated using the molecular-mechanics methods. It is established that 3D-graphite polytypes {alpha}{sub 1,1}, {alpha}{sub 1,3}, {alpha}{sub 1,5}, {alpha}{sub 2,1}, {alpha}{sub 2,3}, {alpha}{sub 3,1}, {beta}{sub 1,2}, {beta}{sub 1,4}, {beta}{sub 1,6}, {beta}{sub 2,1}, and {beta}{sub 3,2} consist of sp{sup 2}-hybridized atoms, have hexagonal unit cells, and differ in regards to the structure of layers and order of their alternation. A possible way to experimentally synthesize new carbon phases is proposed: the polymerization and carbonization of hydrocarbon molecules.
3D Printable Graphene Composite.
Wei, Xiaojun; Li, Dong; Jiang, Wei; Gu, Zheming; Wang, Xiaojuan; Zhang, Zengxing; Sun, Zhengzong
2015-01-01
In human being's history, both the Iron Age and Silicon Age thrived after a matured massive processing technology was developed. Graphene is the most recent superior material which could potentially initialize another new material Age. However, while being exploited to its full extent, conventional processing methods fail to provide a link to today's personalization tide. New technology should be ushered in. Three-dimensional (3D) printing fills the missing linkage between graphene materials and the digital mainstream. Their alliance could generate additional stream to push the graphene revolution into a new phase. Here we demonstrate for the first time, a graphene composite, with a graphene loading up to 5.6 wt%, can be 3D printable into computer-designed models. The composite's linear thermal coefficient is below 75 ppm·°C(-1) from room temperature to its glass transition temperature (Tg), which is crucial to build minute thermal stress during the printing process. PMID:26153673
Rapid high-fidelity visualisation of multispectral 3D mapping
NASA Astrophysics Data System (ADS)
Tudor, Philip M.; Christy, Mark
2011-06-01
Mobile LIDAR scanning typically provides captured 3D data in the form of 3D 'Point Clouds'. Combined with colour imagery these data produce coloured point clouds or, if further processed, polygon-based 3D models. The use of point clouds is simple and rapid, but visualisation can appear ghostly and diffuse. Textured 3D models provide high fidelity visualisation, but their creation is time consuming, difficult to automate and can modify key terrain details. This paper describes techniques for the visualisation of fused multispectral 3D data that approach the visual fidelity of polygon-based models with the rapid turnaround and detail of 3D point clouds. The general approaches to data capture and data fusion are identified as well as the central underlying mathematical transforms, data management and graphics processing techniques used to support rapid, interactive visualisation of very large multispectral 3D datasets. Performance data with respect to real-world 3D mapping as well as illustrations of visualisation outputs are included.
DREAM3D simulations of inner-belt dynamics
Cunningham, Gregory Scott
2015-05-26
A 1973 paper by Lyons and Thorne explains the two-belt structure for electrons in the inner magnetosphere as a balance between inward radial diffusion and loss to the atmosphere, where the loss to the atmosphere is enabled by pitch-angle scattering from Coulomb and wave-particle interactions. In the 1973 paper, equilibrium solutions to a decoupled set of 1D radial diffusion equations, one for each value of the first invariant of motion, μ, were computed to produce the equilibrium two-belt structure. Each 1D radial diffusion equation incorporated an L-and μ-dependent `lifetime' due to the Coulomb and wave-particle interactions. This decoupling of the problem is appropriate under the assumption that radial diffusion is slow in comparison to pitch-angle scattering. However, for some values of μ and L the lifetime associated with pitch-angle scattering is comparable to the timescale associated with radial diffusion, suggesting that the true equilibrium solutions might reflect `coupled modes' involving pitch-angle scattering and radial diffusion and thus requiring a 3D diffusion model. In the work we show here, we have computed the equilibrium solutions using our 3D diffusion model, DREAM3D, that allows for such coupling. We find that the 3D equilibrium solutions are quite similar to the solutions shown in the 1973 paper when we use the same physical models for radial diffusion and pitch-angle scattering from hiss. However, we show that the equilibrium solutions are quite sensitive to various aspects of the physics model employed in the 1973 paper that can be improved, suggesting that additional work needs to be done to understand the two-belt structure.
[Real time 3D echocardiography].
Bauer, F; Shiota, T; Thomas, J D
2001-07-01
Three-dimensional representation of the heart is an old concern. Usually, 3D reconstruction of the cardiac mass is made by successive acquisition of 2D sections, the spatial localisation and orientation of which require complex guiding systems. More recently, the concept of volumetric acquisition has been introduced. A matricial emitter-receiver probe complex with parallel data processing provides instantaneous of a pyramidal 64 degrees x 64 degrees volume. The image is restituted in real time and is composed of 3 planes (planes B and C) which can be displaced in all spatial directions at any time during acquisition. The flexibility of this system of acquisition allows volume and mass measurement with greater accuracy and reproducibility, limiting inter-observer variability. Free navigation of the planes of investigation allows reconstruction for qualitative and quantitative analysis of valvular heart disease and other pathologies. Although real time 3D echocardiography is ready for clinical usage, some improvements are still necessary to improve its conviviality. Then real time 3D echocardiography could be the essential tool for understanding, diagnosis and management of patients. PMID:11494630
[Real time 3D echocardiography
NASA Technical Reports Server (NTRS)
Bauer, F.; Shiota, T.; Thomas, J. D.
2001-01-01
Three-dimensional representation of the heart is an old concern. Usually, 3D reconstruction of the cardiac mass is made by successive acquisition of 2D sections, the spatial localisation and orientation of which require complex guiding systems. More recently, the concept of volumetric acquisition has been introduced. A matricial emitter-receiver probe complex with parallel data processing provides instantaneous of a pyramidal 64 degrees x 64 degrees volume. The image is restituted in real time and is composed of 3 planes (planes B and C) which can be displaced in all spatial directions at any time during acquisition. The flexibility of this system of acquisition allows volume and mass measurement with greater accuracy and reproducibility, limiting inter-observer variability. Free navigation of the planes of investigation allows reconstruction for qualitative and quantitative analysis of valvular heart disease and other pathologies. Although real time 3D echocardiography is ready for clinical usage, some improvements are still necessary to improve its conviviality. Then real time 3D echocardiography could be the essential tool for understanding, diagnosis and management of patients.
PHISICS multi-group transport neutronic capabilities for RELAP5
Epiney, A.; Rabiti, C.; Alfonsi, A.; Wang, Y.; Cogliati, J.; Strydom, G.
2012-07-01
PHISICS is a neutronic code system currently under development at INL. Its goal is to provide state of the art simulation capability to reactor designers. This paper reports on the effort of coupling this package to the thermal hydraulic system code RELAP5. This will enable full prismatic core and system modeling and the possibility to model coupled (thermal-hydraulics and neutronics) problems with more options for 3D neutron kinetics, compared to the existing diffusion theory neutron kinetics module in RELAP5 (NESTLE). The paper describes the capabilities of the coupling and illustrates them with a set of sample problems. (authors)
NASA Astrophysics Data System (ADS)
Kent, G. M.; Harding, A. J.; Babcock, J. M.; Orcutt, J. A.; Bazin, S.; Singh, S.; Detrick, R. S.; Canales, J. P.; Carbotte, S. M.; Diebold, J.
2002-12-01
Multichannel seismic (MCS) images of crustal magma chambers are ideal targets for advanced visualization techniques. In the mid-ocean ridge environment, reflections originating at the melt-lens are well separated from other reflection boundaries, such as the seafloor, layer 2A and Moho, which enables the effective use of transparency filters. 3-D visualization of seismic reflectivity falls into two broad categories: volume and surface rendering. Volumetric-based visualization is an extremely powerful approach for the rapid exploration of very dense 3-D datasets. These 3-D datasets are divided into volume elements or voxels, which are individually color coded depending on the assigned datum value; the user can define an opacity filter to reject plotting certain voxels. This transparency allows the user to peer into the data volume, enabling an easy identification of patterns or relationships that might have geologic merit. Multiple image volumes can be co-registered to look at correlations between two different data types (e.g., amplitude variation with offsets studies), in a manner analogous to draping attributes onto a surface. In contrast, surface visualization of seismic reflectivity usually involves producing "fence" diagrams of 2-D seismic profiles that are complemented with seafloor topography, along with point class data, draped lines and vectors (e.g. fault scarps, earthquake locations and plate-motions). The overlying seafloor can be made partially transparent or see-through, enabling 3-D correlations between seafloor structure and seismic reflectivity. Exploration of 3-D datasets requires additional thought when constructing and manipulating these complex objects. As numbers of visual objects grow in a particular scene, there is a tendency to mask overlapping objects; this clutter can be managed through the effective use of total or partial transparency (i.e., alpha-channel). In this way, the co-variation between different datasets can be investigated
Tomographic system for 3D temperature reconstruction
NASA Astrophysics Data System (ADS)
Antos, Martin; Malina, Radomir
2003-11-01
The novel laboratory system for the optical tomography is used to obtain three-dimensional temperature field around a heated element. The Mach-Zehnder holographic interferometers with diffusive illumination of the phase object provide the possibility to scan of multidirectional holographic interferograms in the range of viewing angles from 0 deg to 108 deg. These interferograms form the input data for the computer tomography of the 3D distribution of the refractive index variation, which characterizes the physical state of the studied medium. The configuration of the system allows automatic projection scanning of the studied phase object. The computer calculates the wavefront deformation for each projection, making use of different methods of Fourier-transform and phase-sampling evaluations. The experimental set-up together with experimental results is presented.
A 3-D microfluidic combinatorial cell array.
Liu, Mike C; Tai, Yu-Chong
2011-02-01
We present the development of a three-dimensional (3-D) combinatorial cell culture array device featured with integrated three-input, eight-output combinatorial mixer and cell culture chambers. The device is designed for cell-based screening of multiple compounds simultaneously on a microfluidic platform. The final assembled device is composed of a porous membrane integrated in between a Parylene 3-D microfluidic chip and a PDMS microfluidic chip. The membrane turned the cell culture chambers into two-level configuration to facilitate cell loading and to maintain cells in a diffusion dominated space during device operation. Experimentally, we first characterized the combined compound concentration profile at each chamber using a fluorescence method. We then successfully demonstrated the functionality of the quantitative cell-based assay by culturing B35 rat neuronal cells on this device and screening the ability of three compounds (1,5-dihydroxyisoquinoline, deferoxamine, and 3-aminobenzoic acid) to attenuate cell death caused by cytotoxic hydrogen peroxide. In another experiment, we assayed for the combinatorial effects of three chemotherapeutic compound exposures (vinorelbine, paclitaxel, and γ-linolenic acid) on human breast cancer cells, MDA-MB-231. The same technology will enable the construction of inexpensive lab-on-a-chip devices with high-input combinatorial mixer for performing high-throughput cell-based assay and highly parallel and combinatorial chemical or biochemical reactions. PMID:21063783
3D geometry applied to atmospheric layers
NASA Astrophysics Data System (ADS)
Nadjib Kouahla, Mohamed; Moreels, Guy; Faivre, Michael
Epipolar geometry is an efficient method for generating 3D representations of objects. Here we present an original application of this method to the case of atmospheric layers. Two synchronized simultaneous images of the same scene are taken in two sites at a distance D. The 36*36 fields of view are oriented face to face along the same line of sight, but in opposite directions. The elevation angle of the optical axis above the horizon is 17. The observed objects are airglow emissions or cirrus clouds or aircraft trails. In the case of clouds, the shape of the objects is diffuse. To obtain a superposition of the common observed zone, it is necessary to calculate a normalized cross-correlation coefficient (NCC) to identify pairs of matching points in both images. The perspective effect in the rectangular images is inverted to produce a satellite-type view of the atmospheric layer as could be seen from an overlying satellite. We developed a triangulation algorithm to retrieve the 3D surface of the observed layer. The stereoscopic method was used to retrieve the wavy structure of the OH emissive layer at the altitude of 87 km. The distance between the observing sites was 600 km. Results obtained in Peru from the sites of Cerro Cosmos and Cerro Verde will be presented. We are currently extending the stereoscopic procedure to the study of troposphere cirruses, of natural origin or induced by aircraft engines. In this case, the distance between observation sites is D 60 km.
Dissipation mechanism in 3D magnetic reconnection
Fujimoto, Keizo
2011-11-15
Dissipation processes responsible for fast magnetic reconnection in collisionless plasmas are investigated using 3D electromagnetic particle-in-cell simulations. The present study revisits the two simulation runs performed in the previous study (Fujimoto, Phys. Plasmas 16, 042103 (2009)); one with small system size in the current density direction, and the other with larger system size. In the case with small system size, the reconnection processes are almost the same as those in 2D reconnection, while in the other case a kink mode evolves along the current density and deforms the current sheet structure drastically. Although fast reconnection is achieved in both the cases, the dissipation mechanism is very different between them. In the case without kink mode, the electrons transit the electron diffusion region without thermalization, so that the magnetic dissipation is supported by the inertia resistivity alone. On the other hand, in the kinked current sheet, the electrons are not only accelerated in bulk, but they are also partly scattered and thermalized by the kink mode, which results in the anomalous resistivity in addition to the inertia resistivity. It is demonstrated that in 3D reconnection the thickness of the electron current sheet becomes larger than the local electron inertia length, consistent with the theoretical prediction in Fujimoto and Sydora (Phys. Plasmas 16, 112309 (2009)).
Interactive 3D Mars Visualization
NASA Technical Reports Server (NTRS)
Powell, Mark W.
2012-01-01
The Interactive 3D Mars Visualization system provides high-performance, immersive visualization of satellite and surface vehicle imagery of Mars. The software can be used in mission operations to provide the most accurate position information for the Mars rovers to date. When integrated into the mission data pipeline, this system allows mission planners to view the location of the rover on Mars to 0.01-meter accuracy with respect to satellite imagery, with dynamic updates to incorporate the latest position information. Given this information so early in the planning process, rover drivers are able to plan more accurate drive activities for the rover than ever before, increasing the execution of science activities significantly. Scientifically, this 3D mapping information puts all of the science analyses to date into geologic context on a daily basis instead of weeks or months, as was the norm prior to this contribution. This allows the science planners to judge the efficacy of their previously executed science observations much more efficiently, and achieve greater science return as a result. The Interactive 3D Mars surface view is a Mars terrain browsing software interface that encompasses the entire region of exploration for a Mars surface exploration mission. The view is interactive, allowing the user to pan in any direction by clicking and dragging, or to zoom in or out by scrolling the mouse or touchpad. This set currently includes tools for selecting a point of interest, and a ruler tool for displaying the distance between and positions of two points of interest. The mapping information can be harvested and shared through ubiquitous online mapping tools like Google Mars, NASA WorldWind, and Worldwide Telescope.
NASA Technical Reports Server (NTRS)
2004-01-01
This 3-D cylindrical-perspective mosaic taken by the navigation camera on the Mars Exploration Rover Spirit on sol 82 shows the view south of the large crater dubbed 'Bonneville.' The rover will travel toward the Columbia Hills, seen here at the upper left. The rock dubbed 'Mazatzal' and the hole the rover drilled in to it can be seen at the lower left. The rover's position is referred to as 'Site 22, Position 32.' This image was geometrically corrected to make the horizon appear flat.
NASA Technical Reports Server (NTRS)
2004-01-01
This is a 3-D anaglyph showing a microscopic image taken of an area measuring 3 centimeters (1.2 inches) across on the rock called Adirondack. The image was taken at Gusev Crater on the 33rd day of the Mars Exploration Rover Spirit's journey (Feb. 5, 2004), after the rover used its rock abrasion tool brush to clean the surface of the rock. Dust, which was pushed off to the side during cleaning, can still be seen to the left and in low areas of the rock.
NASA Technical Reports Server (NTRS)
2004-01-01
This 3-D image captured by the Mars Exploration Rover Opportunity's rear hazard-identification camera shows the now-empty lander that carried the rover 283 million miles to Meridiani Planum, Mars. Engineers received confirmation that Opportunity's six wheels successfully rolled off the lander and onto martian soil at 3:01 a.m. PST, January 31, 2004, on the seventh martian day, or sol, of the mission. The rover is approximately 1 meter (3 feet) in front of the lander, facing north.
NASA Astrophysics Data System (ADS)
Manos, Harry
2016-03-01
Visual aids are important to student learning, and they help make the teacher's job easier. Keeping with the TPT theme of "The Art, Craft, and Science of Physics Teaching," the purpose of this article is to show how teachers, lacking equipment and funds, can construct a durable 3-D model reference frame and a model gravity well tailored to specific class lessons. Most of the supplies are readily available in the home or at school: rubbing alcohol, a rag, two colors of spray paint, art brushes, and masking tape. The cost of these supplies, if you don't have them, is less than 20.
Positional Awareness Map 3D (PAM3D)
NASA Technical Reports Server (NTRS)
Hoffman, Monica; Allen, Earl L.; Yount, John W.; Norcross, April Louise
2012-01-01
The Western Aeronautical Test Range of the National Aeronautics and Space Administration s Dryden Flight Research Center needed to address the aging software and hardware of its current situational awareness display application, the Global Real-Time Interactive Map (GRIM). GRIM was initially developed in the late 1980s and executes on older PC architectures using a Linux operating system that is no longer supported. Additionally, the software is difficult to maintain due to its complexity and loss of developer knowledge. It was decided that a replacement application must be developed or acquired in the near future. The replacement must provide the functionality of the original system, the ability to monitor test flight vehicles in real-time, and add improvements such as high resolution imagery and true 3-dimensional capability. This paper will discuss the process of determining the best approach to replace GRIM, and the functionality and capabilities of the first release of the Positional Awareness Map 3D.
3D Printable Graphene Composite
Wei, Xiaojun; Li, Dong; Jiang, Wei; Gu, Zheming; Wang, Xiaojuan; Zhang, Zengxing; Sun, Zhengzong
2015-01-01
In human being’s history, both the Iron Age and Silicon Age thrived after a matured massive processing technology was developed. Graphene is the most recent superior material which could potentially initialize another new material Age. However, while being exploited to its full extent, conventional processing methods fail to provide a link to today’s personalization tide. New technology should be ushered in. Three-dimensional (3D) printing fills the missing linkage between graphene materials and the digital mainstream. Their alliance could generate additional stream to push the graphene revolution into a new phase. Here we demonstrate for the first time, a graphene composite, with a graphene loading up to 5.6 wt%, can be 3D printable into computer-designed models. The composite’s linear thermal coefficient is below 75 ppm·°C−1 from room temperature to its glass transition temperature (Tg), which is crucial to build minute thermal stress during the printing process. PMID:26153673
3-D Relativistic MHD Simulations
NASA Astrophysics Data System (ADS)
Nishikaw, K.-I.; Frank, J.; Christodoulou, D. M.; Koide, S.; Sakai, J.-I.; Sol, H.; Mutel, R. L.
1998-12-01
We present 3-D numerical simulations of moderately hot, supersonic jets propagating initially along or obliquely to the field lines of a denser magnetized background medium with Lorentz factors of W=4.56 and evolving in a four-dimensional spacetime. The new results are understood as follows: Relativistic simulations have consistently shown that these jets are effectively heavy and so they do not suffer substantial momentum losses and are not decelerated as efficiently as their nonrelativistic counterparts. In addition, the ambient magnetic field, however strong, can be pushed aside with relative ease by the beam, provided that the degrees of freedom associated with all three spatial dimensions are followed self-consistently in the simulations. This effect is analogous to pushing Japanese ``noren'' or vertical Venetian blinds out of the way while the slats are allowed to bend in 3-D space rather than as a 2-D slab structure. We also simulate jets with the more realistic initial conditions for injecting jets for helical mangetic field, perturbed density, velocity, and internal energy, which are supposed to be caused in the process of jet generation. Three possible explanations for the observed variability are (i) tidal disruption of a star falling into the black hole, (ii) instabilities in the relativistic accretion disk, and (iii) jet-related PRocesses. New results will be reported at the meeting.
3D Printable Graphene Composite
NASA Astrophysics Data System (ADS)
Wei, Xiaojun; Li, Dong; Jiang, Wei; Gu, Zheming; Wang, Xiaojuan; Zhang, Zengxing; Sun, Zhengzong
2015-07-01
In human being’s history, both the Iron Age and Silicon Age thrived after a matured massive processing technology was developed. Graphene is the most recent superior material which could potentially initialize another new material Age. However, while being exploited to its full extent, conventional processing methods fail to provide a link to today’s personalization tide. New technology should be ushered in. Three-dimensional (3D) printing fills the missing linkage between graphene materials and the digital mainstream. Their alliance could generate additional stream to push the graphene revolution into a new phase. Here we demonstrate for the first time, a graphene composite, with a graphene loading up to 5.6 wt%, can be 3D printable into computer-designed models. The composite’s linear thermal coefficient is below 75 ppm·°C-1 from room temperature to its glass transition temperature (Tg), which is crucial to build minute thermal stress during the printing process.
3D acoustic atmospheric tomography
NASA Astrophysics Data System (ADS)
Rogers, Kevin; Finn, Anthony
2014-10-01
This paper presents a method for tomographically reconstructing spatially varying 3D atmospheric temperature profiles and wind velocity fields based. Measurements of the acoustic signature measured onboard a small Unmanned Aerial Vehicle (UAV) are compared to ground-based observations of the same signals. The frequency-shifted signal variations are then used to estimate the acoustic propagation delay between the UAV and the ground microphones, which are also affected by atmospheric temperature and wind speed vectors along each sound ray path. The wind and temperature profiles are modelled as the weighted sum of Radial Basis Functions (RBFs), which also allow local meteorological measurements made at the UAV and ground receivers to supplement any acoustic observations. Tomography is used to provide a full 3D reconstruction/visualisation of the observed atmosphere. The technique offers observational mobility under direct user control and the capacity to monitor hazardous atmospheric environments, otherwise not justifiable on the basis of cost or risk. This paper summarises the tomographic technique and reports on the results of simulations and initial field trials. The technique has practical applications for atmospheric research, sound propagation studies, boundary layer meteorology, air pollution measurements, analysis of wind shear, and wind farm surveys.
Mannoor, Manu S.; Jiang, Ziwen; James, Teena; Kong, Yong Lin; Malatesta, Karen A.; Soboyejo, Winston O.; Verma, Naveen; Gracias, David H.; McAlpine, Michael C.
2013-01-01
The ability to three-dimensionally interweave biological tissue with functional electronics could enable the creation of bionic organs possessing enhanced functionalities over their human counterparts. Conventional electronic devices are inherently two-dimensional, preventing seamless multidimensional integration with synthetic biology, as the processes and materials are very different. Here, we present a novel strategy for overcoming these difficulties via additive manufacturing of biological cells with structural and nanoparticle derived electronic elements. As a proof of concept, we generated a bionic ear via 3D printing of a cell-seeded hydrogel matrix in the precise anatomic geometry of a human ear, along with an intertwined conducting polymer consisting of infused silver nanoparticles. This allowed for in vitro culturing of cartilage tissue around an inductive coil antenna in the ear, which subsequently enables readout of inductively-coupled signals from cochlea-shaped electrodes. The printed ear exhibits enhanced auditory sensing for radio frequency reception, and complementary left and right ears can listen to stereo audio music. Overall, our approach suggests a means to intricately merge biologic and nanoelectronic functionalities via 3D printing. PMID:23635097
3D Ion Temperature Reconstruction
NASA Astrophysics Data System (ADS)
Tanabe, Hiroshi; You, Setthivoine; Balandin, Alexander; Inomoto, Michiaki; Ono, Yasushi
2009-11-01
The TS-4 experiment at the University of Tokyo collides two spheromaks to form a single high-beta compact toroid. Magnetic reconnection during the merging process heats and accelerates the plasma in toroidal and poloidal directions. The reconnection region has a complex 3D topology determined by the pitch of the spheromak magnetic fields at the merging plane. A pair of multichord passive spectroscopic diagnostics have been established to measure the ion temperature and velocity in the reconnection volume. One setup measures spectral lines across a poloidal plane, retrieving velocity and temperature from Abel inversion. The other, novel setup records spectral lines across another section of the plasma and reconstructs velocity and temperature from 3D vector and 2D scalar tomography techniques. The magnetic field linking both measurement planes is determined from in situ magnetic probe arrays. The ion temperature is then estimated within the volume between the two measurement planes and at the reconnection region. The measurement is followed over several repeatable discharges to follow the heating and acceleration process during the merging reconnection.
3D medical thermography device
NASA Astrophysics Data System (ADS)
Moghadam, Peyman
2015-05-01
In this paper, a novel handheld 3D medical thermography system is introduced. The proposed system consists of a thermal-infrared camera, a color camera and a depth camera rigidly attached in close proximity and mounted on an ergonomic handle. As a practitioner holding the device smoothly moves it around the human body parts, the proposed system generates and builds up a precise 3D thermogram model by incorporating information from each new measurement in real-time. The data is acquired in motion, thus it provides multiple points of view. When processed, these multiple points of view are adaptively combined by taking into account the reliability of each individual measurement which can vary due to a variety of factors such as angle of incidence, distance between the device and the subject and environmental sensor data or other factors influencing a confidence of the thermal-infrared data when captured. Finally, several case studies are presented to support the usability and performance of the proposed system.
Mannoor, Manu S; Jiang, Ziwen; James, Teena; Kong, Yong Lin; Malatesta, Karen A; Soboyejo, Winston O; Verma, Naveen; Gracias, David H; McAlpine, Michael C
2013-06-12
The ability to three-dimensionally interweave biological tissue with functional electronics could enable the creation of bionic organs possessing enhanced functionalities over their human counterparts. Conventional electronic devices are inherently two-dimensional, preventing seamless multidimensional integration with synthetic biology, as the processes and materials are very different. Here, we present a novel strategy for overcoming these difficulties via additive manufacturing of biological cells with structural and nanoparticle derived electronic elements. As a proof of concept, we generated a bionic ear via 3D printing of a cell-seeded hydrogel matrix in the anatomic geometry of a human ear, along with an intertwined conducting polymer consisting of infused silver nanoparticles. This allowed for in vitro culturing of cartilage tissue around an inductive coil antenna in the ear, which subsequently enables readout of inductively-coupled signals from cochlea-shaped electrodes. The printed ear exhibits enhanced auditory sensing for radio frequency reception, and complementary left and right ears can listen to stereo audio music. Overall, our approach suggests a means to intricately merge biologic and nanoelectronic functionalities via 3D printing. PMID:23635097
Larry Lawrence; Bruce Miller
2004-09-01
The Lott Ranch 3D seismic prospect located in Garza County, Texas is a project initiated in September of 1991 by the J.M. Huber Corp., a petroleum exploration and production company. By today's standards the 126 square mile project does not seem monumental, however at the time it was conceived it was the most intensive land 3D project ever attempted. Acquisition began in September of 1991 utilizing GEO-SEISMIC, INC., a seismic data contractor. The field parameters were selected by J.M. Huber, and were of a radical design. The recording instruments used were GeoCor IV amplifiers designed by Geosystems Inc., which record the data in signed bit format. It would not have been practical, if not impossible, to have processed the entire raw volume with the tools available at that time. The end result was a dataset that was thought to have little utility due to difficulties in processing the field data. In 1997, Yates Energy Corp. located in Roswell, New Mexico, formed a partnership to further develop the project. Through discussions and meetings with Pinnacle Seismic, it was determined that the original Lott Ranch 3D volume could be vastly improved upon reprocessing. Pinnacle Seismic had shown the viability of improving field-summed signed bit data on smaller 2D and 3D projects. Yates contracted Pinnacle Seismic Ltd. to perform the reprocessing. This project was initiated with high resolution being a priority. Much of the potential resolution was lost through the initial summing of the field data. Modern computers that are now being utilized have tremendous speed and storage capacities that were cost prohibitive when this data was initially processed. Software updates and capabilities offer a variety of quality control and statics resolution, which are pertinent to the Lott Ranch project. The reprocessing effort was very successful. The resulting processed data-set was then interpreted using modern PC-based interpretation and mapping software. Production data, log data
Goyanes, Alvaro; Det-Amornrat, Usanee; Wang, Jie; Basit, Abdul W; Gaisford, Simon
2016-07-28
Acne is a multifactorial inflammatory skin disease with high prevalence. In this work, the potential of 3D printing to produce flexible personalised-shape anti-acne drug (salicylic acid) loaded devices was demonstrated by two different 3D printing (3DP) technologies: Fused Deposition Modelling (FDM) and stereolithography (SLA). 3D scanning technology was used to obtain a 3D model of a nose adapted to the morphology of an individual. In FDM 3DP, commercially produced Flex EcoPLA™ (FPLA) and polycaprolactone (PCL) filaments were loaded with salicylic acid by hot melt extrusion (HME) (theoretical drug loading - 2% w/w) and used as feedstock material for 3D printing. Drug loading in the FPLA-salicylic acid and PCL-salicylic acid 3D printed patches was 0.4% w/w and 1.2% w/w respectively, indicating significant thermal degradation of drug during HME and 3D printing. Diffusion testing in Franz cells using a synthetic membrane revealed that the drug loaded printed samples released <187μg/cm(2) within 3h. FPLA-salicylic acid filament was successfully printed as a nose-shape mask by FDM 3DP, but the PCL-salicylic acid filament was not. In the SLA printing process, the drug was dissolved in different mixtures of poly(ethylene glycol) diacrylate (PEGDA) and poly(ethylene glycol) (PEG) that were solidified by the action of a laser beam. SLA printing led to 3D printed devices (nose-shape) with higher resolution and higher drug loading (1.9% w/w) than FDM, with no drug degradation. The results of drug diffusion tests revealed that drug diffusion was faster than with the FDM devices, 229 and 291μg/cm(2) within 3h for the two formulations evaluated. In this study, SLA printing was the more appropriate 3D printing technology to manufacture anti-acne devices with salicylic acid. The combination of 3D scanning and 3D printing has the potential to offer solutions to produce personalised drug loaded devices, adapted in shape and size to individual patients. PMID:27189134
Smith, L.A.; Gallmeier, F.X.; Gehin, J.C.
1995-05-01
The FOEHN critical experiment was analyzed to validate the use of multigroup cross sections and Oak Ridge National Laboratory neutronics computer codes in the design of the Advanced Neutron Source. The ANSL-V 99-group master cross section library was used for all the calculations. Three different critical configurations were evaluated using the multigroup KENO Monte Carlo transport code, the multigroup DORT discrete ordinates transport code, and the multigroup diffusion theory code VENTURE. The simple configuration consists of only the fuel and control elements with the heavy water reflector. The intermediate configuration includes boron endplates at the upper and lower edges of the fuel element. The complex configuration includes both the boron endplates and components in the reflector. Cross sections were processed using modules from the AMPX system. Both 99-group and 20-group cross sections were created and used in two-dimensional models of the FOEHN experiment. KENO calculations were performed using both 99-group and 20-group cross sections. The DORT and VENTURE calculations were performed using 20-group cross sections. Because the simple and intermediate configurations are azimuthally symmetric, these configurations can be explicitly modeled in R-Z geometry. Since the reflector components cannot be modeled explicitly using the current versions of these codes, three reflector component homogenization schemes were developed and evaluated for the complex configuration. Power density distributions were calculated with KENO using 99-group cross sections and with DORT and VENTURE using 20-group cross sections. The average differences between the measured values and the values calculated with the different computer codes range from 2.45 to 5.74%. The maximum differences between the measured and calculated thermal flux values for the simple and intermediate configurations are {approx} 13%, while the average differences are < 8%.
3D Printing of Graphene Aerogels.
Zhang, Qiangqiang; Zhang, Feng; Medarametla, Sai Pradeep; Li, Hui; Zhou, Chi; Lin, Dong
2016-04-01
3D printing of a graphene aerogel with true 3D overhang structures is highlighted. The aerogel is fabricated by combining drop-on-demand 3D printing and freeze casting. The water-based GO ink is ejected and freeze-cast into designed 3D structures. The lightweight (<10 mg cm(-3) ) 3D printed graphene aerogel presents superelastic and high electrical conduction. PMID:26861680
Non-isothermal 3D SDPD Simulations
NASA Astrophysics Data System (ADS)
Yang, Jun; Potami, Raffaele; Gatsonis, Nikolaos
2012-11-01
The study of fluids at micro and nanoscale requires new modeling and computational approaches. Smooth Particle Dissipative Dynamics (SDPD) is a mesh-free method that provides a bridge between the continuum equations of hydrodynamics embedded in the Smooth Particle Hydrodynamics approach and the molecular nature embedded in the DPD approach. SDPD is thermodynamically consistent, does not rely on arbitrary coefficients for its thermostat, involves realistic transport coefficients, and includes fluctuation terms. SDPD is implemented in our work for arbitrary 3D geometries with a methodology to model solid wall boundary conditions. We present simulations for isothermal flows for verification of our approach. The entropy equation is implemented with a velocity-entropy Verlet integration algorithm Flows with heat transfer are simulated for verification of the SDPD. We present also the self-diffusion coefficient derived from SDPD simulations for gases and liquids. Results show the scale dependence of self-diffusion coefficient on SDPD particle size. Computational Mathematics Program of the Air Force Office of Scientific Research under grant/contract number FA9550-06-1-0236.
Sinclair, Michael B
2012-01-05
ShowMe3D is a data visualization graphical user interface specifically designed for use with hyperspectral image obtained from the Hyperspectral Confocal Microscope. The program allows the user to select and display any single image from a three dimensional hyperspectral image stack. By moving a slider control, the user can easily move between images of the stack. The user can zoom into any region of the image. The user can select any pixel or region from the displayed image and display the fluorescence spectrum associated with that pixel or region. The user can define up to 3 spectral filters to apply to the hyperspectral image and view the image as it would appear from a filter-based confocal microscope. The user can also obtain statistics such as intensity average and variance from selected regions.
3D Elastic Wavefield Tomography
NASA Astrophysics Data System (ADS)
Guasch, L.; Warner, M.; Stekl, I.; Umpleby, A.; Shah, N.
2010-12-01
Wavefield tomography, or waveform inversion, aims to extract the maximum information from seismic data by matching trace by trace the response of the solid earth to seismic waves using numerical modelling tools. Its first formulation dates from the early 80's, when Albert Tarantola developed a solid theoretical basis that is still used today with little change. Due to computational limitations, the application of the method to 3D problems has been unaffordable until a few years ago, and then only under the acoustic approximation. Although acoustic wavefield tomography is widely used, a complete solution of the seismic inversion problem requires that we account properly for the physics of wave propagation, and so must include elastic effects. We have developed a 3D tomographic wavefield inversion code that incorporates the full elastic wave equation. The bottle neck of the different implementations is the forward modelling algorithm that generates the synthetic data to be compared with the field seismograms as well as the backpropagation of the residuals needed to form the direction update of the model parameters. Furthermore, one or two extra modelling runs are needed in order to calculate the step-length. Our approach uses a FD scheme explicit time-stepping by finite differences that are 4th order in space and 2nd order in time, which is a 3D version of the one developed by Jean Virieux in 1986. We chose the time domain because an explicit time scheme is much less demanding in terms of memory than its frequency domain analogue, although the discussion of wich domain is more efficient still remains open. We calculate the parameter gradients for Vp and Vs by correlating the normal and shear stress wavefields respectively. A straightforward application would lead to the storage of the wavefield at all grid points at each time-step. We tackled this problem using two different approaches. The first one makes better use of resources for small models of dimension equal
Energy Science and Technology Software Center (ESTSC)
2012-01-05
ShowMe3D is a data visualization graphical user interface specifically designed for use with hyperspectral image obtained from the Hyperspectral Confocal Microscope. The program allows the user to select and display any single image from a three dimensional hyperspectral image stack. By moving a slider control, the user can easily move between images of the stack. The user can zoom into any region of the image. The user can select any pixel or region from themore » displayed image and display the fluorescence spectrum associated with that pixel or region. The user can define up to 3 spectral filters to apply to the hyperspectral image and view the image as it would appear from a filter-based confocal microscope. The user can also obtain statistics such as intensity average and variance from selected regions.« less
RZ calculations for self shielded multigroup cross sections
Li, M.; Sanchez, R.; Zmijarevic, I.; Stankovski, Z.
2006-07-01
A collision probability method has been implemented for RZ geometries. The method accounts for white albedo, specular and translation boundary condition on the top and bottom surfaces of the geometry and for a white albedo condition on the outer radial surface. We have applied the RZ CP method to the calculation of multigroup self shielded cross sections for Gadolinia absorbers in BWRs. (authors)
Multigroup Confirmatory Factor Analysis: Locating the Invariant Referent Sets
ERIC Educational Resources Information Center
French, Brian F.; Finch, W. Holmes
2008-01-01
Multigroup confirmatory factor analysis (MCFA) is a popular method for the examination of measurement invariance and specifically, factor invariance. Recent research has begun to focus on using MCFA to detect invariance for test items. MCFA requires certain parameters (e.g., factor loadings) to be constrained for model identification, which are…
3D multiplexed immunoplasmonics microscopy
NASA Astrophysics Data System (ADS)
Bergeron, Éric; Patskovsky, Sergiy; Rioux, David; Meunier, Michel
2016-07-01
Selective labelling, identification and spatial distribution of cell surface biomarkers can provide important clinical information, such as distinction between healthy and diseased cells, evolution of a disease and selection of the optimal patient-specific treatment. Immunofluorescence is the gold standard for efficient detection of biomarkers expressed by cells. However, antibodies (Abs) conjugated to fluorescent dyes remain limited by their photobleaching, high sensitivity to the environment, low light intensity, and wide absorption and emission spectra. Immunoplasmonics is a novel microscopy method based on the visualization of Abs-functionalized plasmonic nanoparticles (fNPs) targeting cell surface biomarkers. Tunable fNPs should provide higher multiplexing capacity than immunofluorescence since NPs are photostable over time, strongly scatter light at their plasmon peak wavelengths and can be easily functionalized. In this article, we experimentally demonstrate accurate multiplexed detection based on the immunoplasmonics approach. First, we achieve the selective labelling of three targeted cell surface biomarkers (cluster of differentiation 44 (CD44), epidermal growth factor receptor (EGFR) and voltage-gated K+ channel subunit KV1.1) on human cancer CD44+ EGFR+ KV1.1+ MDA-MB-231 cells and reference CD44- EGFR- KV1.1+ 661W cells. The labelling efficiency with three stable specific immunoplasmonics labels (functionalized silver nanospheres (CD44-AgNSs), gold (Au) NSs (EGFR-AuNSs) and Au nanorods (KV1.1-AuNRs)) detected by reflected light microscopy (RLM) is similar to the one with immunofluorescence. Second, we introduce an improved method for 3D localization and spectral identification of fNPs based on fast z-scanning by RLM with three spectral filters corresponding to the plasmon peak wavelengths of the immunoplasmonics labels in the cellular environment (500 nm for 80 nm AgNSs, 580 nm for 100 nm AuNSs and 700 nm for 40 nm × 92 nm AuNRs). Third, the developed