Development and application of a 3D Cartesian grid Euler method
NASA Technical Reports Server (NTRS)
Melton, John E.; Aftosmis, Michael J.; Berger, Marsha J.; Wong, Michael D.
1995-01-01
This report describes recent progress in the development and application of 3D Cartesian grid generation and Euler flow solution techniques. Improvements to flow field grid generation algorithms, geometry representations, and geometry refinement criteria are presented, including details of a procedure for correctly identifying and resolving extremely thin surface features. An initial implementation of automatic flow field refinement is also presented. Results for several 3D multi-component configurations are provided and discussed.
Development and Applications of 3D Cartesian CFD Technology
NASA Technical Reports Server (NTRS)
Melton, John E.; Berger, Marsha J.; VanDalsem, William (Technical Monitor)
1994-01-01
The urgent need for dramatic reductions in aircraft design cycle time is focusing scrutiny upon all aspects of computational fluid dynamics (CFD). These reductions will most likely come not from increased reliance upon user-interactive (and therefore time-expensive) methods, but instead from methods that can be fully automated and incorporated into 'black box' solutions. In comparison with tetrahedral methods, three-dimensional Cartesian grid approaches are in relative infancy, but initial experiences with automated Cartesian techniques are quite promising. Our research is targeted at furthering the development of Cartesian methods so that they can become key elements of a completely automatic grid generation/flow solution procedure applicable to the Euler analysis of complex aircraft geometries.
The 3D Euler solutions using automated Cartesian grid generation
NASA Technical Reports Server (NTRS)
Melton, John E.; Enomoto, Francis Y.; Berger, Marsha J.
1993-01-01
Viewgraphs on 3-dimensional Euler solutions using automated Cartesian grid generation are presented. Topics covered include: computational fluid dynamics (CFD) and the design cycle; Cartesian grid strategy; structured body fit; grid generation; prolate spheroid; and ONERA M6 wing.
Shared Memory Parallelism for 3D Cartesian Discrete Ordinates Solver
NASA Astrophysics Data System (ADS)
Moustafa, Salli; Dutka-Malen, Ivan; Plagne, Laurent; Ponçot, Angélique; Ramet, Pierre
2014-06-01
This paper describes the design and the performance of DOMINO, a 3D Cartesian SN solver that implements two nested levels of parallelism (multicore+SIMD) on shared memory computation nodes. DOMINO is written in C++, a multi-paradigm programming language that enables the use of powerful and generic parallel programming tools such as Intel TBB and Eigen. These two libraries allow us to combine multi-thread parallelism with vector operations in an efficient and yet portable way. As a result, DOMINO can exploit the full power of modern multi-core processors and is able to tackle very large simulations, that usually require large HPC clusters, using a single computing node. For example, DOMINO solves a 3D full core PWR eigenvalue problem involving 26 energy groups, 288 angular directions (S16), 46 × 106 spatial cells and 1 × 1012 DoFs within 11 hours on a single 32-core SMP node. This represents a sustained performance of 235 GFlops and 40:74% of the SMP node peak performance for the DOMINO sweep implementation. The very high Flops/Watt ratio of DOMINO makes it a very interesting building block for a future many-nodes nuclear simulation tool.
Triangle geometry processing for surface modeling and cartesian grid generation
Aftosmis, Michael J [San Mateo, CA; Melton, John E [Hollister, CA; Berger, Marsha J [New York, NY
2002-09-03
Cartesian mesh generation is accomplished for component based geometries, by intersecting components subject to mesh generation to extract wetted surfaces with a geometry engine using adaptive precision arithmetic in a system which automatically breaks ties with respect to geometric degeneracies. During volume mesh generation, intersected surface triangulations are received to enable mesh generation with cell division of an initially coarse grid. The hexagonal cells are resolved, preserving the ability to directionally divide cells which are locally well aligned.
Triangle Geometry Processing for Surface Modeling and Cartesian Grid Generation
NASA Technical Reports Server (NTRS)
Aftosmis, Michael J. (Inventor); Melton, John E. (Inventor); Berger, Marsha J. (Inventor)
2002-01-01
Cartesian mesh generation is accomplished for component based geometries, by intersecting components subject to mesh generation to extract wetted surfaces with a geometry engine using adaptive precision arithmetic in a system which automatically breaks ties with respect to geometric degeneracies. During volume mesh generation, intersected surface triangulations are received to enable mesh generation with cell division of an initially coarse grid. The hexagonal cells are resolved, preserving the ability to directionally divide cells which are locally well aligned.
Courau, T.; Moustafa, S.; Plagne, L.; Poncot, A.
2013-07-01
As part of its activity, EDF R and D is developing a new nuclear core simulation code named COCAGNE. This code relies on DIABOLO, a Simplified PN (SPN) method to compute the neutron flux inside the core for eigenvalue calculations. In order to assess the accuracy of SPN calculations, we have developed DOMINO, a new 3D Cartesian SN solver. The parallel implementation of DOMINO is very efficient and allows to complete an eigenvalue calculation involving around 300 x 10{sup 9} degrees of freedom within a few hours on a single shared-memory supercomputing node. This computation corresponds to a 26-group S{sub 8} 3D PWR core model used to assess the SPN accuracy. At the pin level, the maximal error for the SP{sub 5} DIABOLO fission production rate is lower than 0.2% compared to the S{sub 8} DOMINO reference for this 3D PWR core model. (authors)
Discovering Structural Regularity in 3D Geometry
Pauly, Mark; Mitra, Niloy J.; Wallner, Johannes; Pottmann, Helmut; Guibas, Leonidas J.
2010-01-01
We introduce a computational framework for discovering regular or repeated geometric structures in 3D shapes. We describe and classify possible regular structures and present an effective algorithm for detecting such repeated geometric patterns in point- or mesh-based models. Our method assumes no prior knowledge of the geometry or spatial location of the individual elements that define the pattern. Structure discovery is made possible by a careful analysis of pairwise similarity transformations that reveals prominent lattice structures in a suitable model of transformation space. We introduce an optimization method for detecting such uniform grids specifically designed to deal with outliers and missing elements. This yields a robust algorithm that successfully discovers complex regular structures amidst clutter, noise, and missing geometry. The accuracy of the extracted generating transformations is further improved using a novel simultaneous registration method in the spatial domain. We demonstrate the effectiveness of our algorithm on a variety of examples and show applications to compression, model repair, and geometry synthesis. PMID:21170292
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.
Tensor decomposition in electronic structure calculations on 3D Cartesian grids
Khoromskij, B.N. Khoromskaia, V.; Chinnamsetty, S.R.; Flad, H.-J.
2009-09-01
In this paper, we investigate a novel approach based on the combination of Tucker-type and canonical tensor decomposition techniques for the efficient numerical approximation of functions and operators in electronic structure calculations. In particular, we study applicability of tensor approximations for the numerical solution of Hartree-Fock and Kohn-Sham equations on 3D Cartesian grids. We show that the orthogonal Tucker-type tensor approximation of electron density and Hartree potential of simple molecules leads to low tensor rank representations. This enables an efficient tensor-product convolution scheme for the computation of the Hartree potential using a collocation-type approximation via piecewise constant basis functions on a uniform nxnxn grid. Combined with the Richardson extrapolation, our approach exhibits O(h{sup 3}) convergence in the grid-size h=O(n{sup -1}). Moreover, this requires O(3rn+r{sup 3}) storage, where r denotes the Tucker rank of the electron density with r=O(logn), almost uniformly in n. For example, calculations of the Coulomb matrix and the Hartree-Fock energy for the CH{sub 4} molecule, with a pseudopotential on the C atom, achieved accuracies of the order of 10{sup -6} hartree with a grid-size n of several hundreds. Since the tensor-product convolution in 3D is performed via 1D convolution transforms, our scheme markedly outperforms the 3D-FFT in both the computing time and storage requirements.
Using Cabri3D Diagrams for Teaching Geometry
ERIC Educational Resources Information Center
Accascina, Giuseppe; Rogora, Enrico
2006-01-01
Cabri3D is a potentially very useful software for learning and teaching 3D geometry. The dynamic nature of the digital diagrams produced with it provides a useful aid for helping students to better develop concept images of geometric concepts. However, since any Cabri3D diagram represents three-dimensional objects on the two dimensional screen of…
Unit cell geometry of 3-D braided structures
NASA Technical Reports Server (NTRS)
Du, Guang-Wu; Ko, Frank K.
1993-01-01
The traditional approach used in modeling of composites reinforced by three-dimensional (3-D) braids is to assume a simple unit cell geometry of a 3-D braided structure with known fiber volume fraction and orientation. In this article, we first examine 3-D braiding methods in the light of braid structures, followed by the development of geometric models for 3-D braids using a unit cell approach. The unit cell geometry of 3-D braids is identified and the relationship of structural parameters such as yarn orientation angle and fiber volume fraction with the key processing parameters established. The limiting geometry has been computed by establishing the point at which yarns jam against each other. Using this factor makes it possible to identify the complete range of allowable geometric arrangements for 3-D braided preforms. This identified unit cell geometry can be translated to mechanical models which relate the geometrical properties of fabric preforms to the mechanical responses of composite systems.
3D Euler flow solutions using unstructured Cartesian and prismatic grids
NASA Technical Reports Server (NTRS)
Melton, John E.; Pandya, Shishir A.; Steger, Joseph L.
1993-01-01
A hyperbolic prismatic grid generation technique is combined with a background Cartesian grid for the study of inviscid three-dimensional flows. The mathematics of the hyperbolic prismatic grid generation algorithm are described, and some simple inviscid demonstration cases are presented. By combining the simplicity of the Cartesian background grid with the geometric flexibility and computational efficiencies inherent to prismatic grids, this approach shows promise for improving computational aerodynamic simulations.
A 3D Geometry Model Search Engine to Support Learning
ERIC Educational Resources Information Center
Tam, Gary K. L.; Lau, Rynson W. H.; Zhao, Jianmin
2009-01-01
Due to the popularity of 3D graphics in animation and games, usage of 3D geometry deformable models increases dramatically. Despite their growing importance, these models are difficult and time consuming to build. A distance learning system for the construction of these models could greatly facilitate students to learn and practice at different…
Using 3D Geometric Models to Teach Spatial Geometry Concepts.
ERIC Educational Resources Information Center
Bertoline, Gary R.
1991-01-01
An explanation of 3-D Computer Aided Design (CAD) usage to teach spatial geometry concepts using nontraditional techniques is presented. The software packages CADKEY and AutoCAD are described as well as their usefulness in solving space geometry problems. (KR)
Automatic visualization of 3D geometry contained in online databases
NASA Astrophysics Data System (ADS)
Zhang, Jie; John, Nigel W.
2003-04-01
In this paper, the application of the Virtual Reality Modeling Language (VRML) for efficient database visualization is analyzed. With the help of JAVA programming, three examples of automatic visualization from a database containing 3-D Geometry are given. The first example is used to create basic geometries. The second example is used to create cylinders with a defined start point and end point. The third example is used to processs data from an old copper mine complex in Cheshire, United Kingdom. Interactive 3-D visualization of all geometric data in an online database is achieved with JSP technology.
Solwnd: A 3D Compressible MHD Code for Solar Wind Studies. Version 1.0: Cartesian Coordinates
NASA Technical Reports Server (NTRS)
Deane, Anil E.
1996-01-01
Solwnd 1.0 is a three-dimensional compressible MHD code written in Fortran for studying the solar wind. Time-dependent boundary conditions are available. The computational algorithm is based on Flux Corrected Transport and the code is based on the existing code of Zalesak and Spicer. The flow considered is that of shear flow with incoming flow that perturbs this base flow. Several test cases corresponding to pressure balanced magnetic structures with velocity shear flow and various inflows including Alfven waves are presented. Version 1.0 of solwnd considers a rectangular Cartesian geometry. Future versions of solwnd will consider a spherical geometry. Some discussions of this issue is presented.
A fast nested dissection solver for Cartesian 3D elliptic problems using hierarchical matrices
NASA Astrophysics Data System (ADS)
Schmitz, Phillip G.; Ying, Lexing
2014-02-01
We present a fast algorithm for solutions to linear systems arising from three dimensional elliptic problems on a regular Cartesian mesh. We follow the approach of Schmitz and Ying (2012) on combining the nested dissection matrix factorization method with hierarchical matrices in two dimensions and extend it to the three dimensional case. A theoretical linear time complexity is derived and a more practical variant with slightly worse scaling is demonstrated.
A Parallel Cartesian Approach for External Aerodynamics of Vehicles with Complex Geometry
NASA Technical Reports Server (NTRS)
Aftosmis, M. J.; Berger, M. J.; Adomavicius, G.
2001-01-01
This workshop paper presents the current status in the development of a new approach for the solution of the Euler equations on Cartesian meshes with embedded boundaries in three dimensions on distributed and shared memory architectures. The approach uses adaptively refined Cartesian hexahedra to fill the computational domain. Where these cells intersect the geometry, they are cut by the boundary into arbitrarily shaped polyhedra which receive special treatment by the solver. The presentation documents a newly developed multilevel upwind solver based on a flexible domain-decomposition strategy. One novel aspect of the work is its use of space-filling curves (SFC) for memory efficient on-the-fly parallelization, dynamic re-partitioning and automatic coarse mesh generation. Within each subdomain the approach employs a variety reordering techniques so that relevant data are on the same page in memory permitting high-performance on cache-based processors. Details of the on-the-fly SFC based partitioning are presented as are construction rules for the automatic coarse mesh generation. After describing the approach, the paper uses model problems and 3- D configurations to both verify and validate the solver. The model problems demonstrate that second-order accuracy is maintained despite the presence of the irregular cut-cells in the mesh. In addition, it examines both parallel efficiency and convergence behavior. These investigations demonstrate a parallel speed-up in excess of 28 on 32 processors of an SGI Origin 2000 system and confirm that mesh partitioning has no effect on convergence behavior.
Onset of buoyancy-driven convection in Cartesian and cylindrical geometries
NASA Astrophysics Data System (ADS)
Myint, Philip C.; Firoozabadi, Abbas
2013-04-01
We perform a linear stability analysis to examine the onset of buoyancy-driven convection relevant to subsurface carbon dioxide sequestration in confined, porous Cartesian and cylindrical domains. Our work amends the analysis in an earlier study on cylindrical geometries. We consider Cartesian geometries where the aspect ratio between the two horizontal dimensions is not necessarily equal to one. Two key elements of the stability analysis are: (1) the critical time and (2) the critical wavenumber. Lateral boundaries have a much greater influence on the critical wavenumber than on the critical time. The confinement due to these boundaries impedes the onset of convection to the extent that convection cannot even occur in domains that are smaller than a certain size. Large aspect ratios can significantly reduce boundary effects. Patterns of the earliest-growing perturbation mode in the horizontal plane reveal many interesting dynamics which have not been examined in previous stability analyses. We illustrate several differences between patterns in Cartesian geometries and patterns in cylindrical geometries. Based on observations from earlier papers, we hypothesize that the contrasts between the Cartesian and cylindrical patterns may lead to significantly different behavior in the two geometries after the onset of convection. Our results may guide future numerical studies that can investigate this hypothesis and may help with understanding the onset of buoyancy-driven convection in real systems where lateral boundary effects are significant.
Implicit Approaches for Moving Boundaries in a 3-D Cartesian Method
NASA Technical Reports Server (NTRS)
Murman, Scott M.; Aftosmis, Michael J.; Berger, Marsha J.; Kwak, Dochan
2003-01-01
This work considers numerical simulation of three-dimensional flows with time-evolving boundaries. Such problems pose a variety of challenges for numerical schemes, and have received a substantial amount of attention in the recent literature. Since such simulations are unsteady, time-accurate solution of the governing equations is required. In special cases, the body motion can be treated by a uniform rigid motion of the computational domain. For the more general situation of relative-body motion, however, this simplification is unavailable and the simulations require a mechanism for ensuring that the mesh evolves with the moving boundaries. This involves a "remeshing" of the computational domain (either localized or global) at each physical timestep, and places a premium on both the speed and robustness of the remeshing algorithms. This work presents a method which includes unsteady flow simulation, rigid domain motion, and relative body motion using a time-evolving Cartesian grid system in three dimensions.
Novel 3D Compression Methods for Geometry, Connectivity and Texture
NASA Astrophysics Data System (ADS)
Siddeq, M. M.; Rodrigues, M. A.
2016-06-01
A large number of applications in medical visualization, games, engineering design, entertainment, heritage, e-commerce and so on require the transmission of 3D models over the Internet or over local networks. 3D data compression is an important requirement for fast data storage, access and transmission within bandwidth limitations. The Wavefront OBJ (object) file format is commonly used to share models due to its clear simple design. Normally each OBJ file contains a large amount of data (e.g. vertices and triangulated faces, normals, texture coordinates and other parameters) describing the mesh surface. In this paper we introduce a new method to compress geometry, connectivity and texture coordinates by a novel Geometry Minimization Algorithm (GM-Algorithm) in connection with arithmetic coding. First, each vertex ( x, y, z) coordinates are encoded to a single value by the GM-Algorithm. Second, triangle faces are encoded by computing the differences between two adjacent vertex locations, which are compressed by arithmetic coding together with texture coordinates. We demonstrate the method on large data sets achieving compression ratios between 87 and 99 % without reduction in the number of reconstructed vertices and triangle faces. The decompression step is based on a Parallel Fast Matching Search Algorithm (Parallel-FMS) to recover the structure of the 3D mesh. A comparative analysis of compression ratios is provided with a number of commonly used 3D file formats such as VRML, OpenCTM and STL highlighting the performance and effectiveness of the proposed method.
NASA Astrophysics Data System (ADS)
Gilmanov, Anvar; Sotiropoulos, Fotis
2005-08-01
A numerical method is developed for solving the 3D, unsteady, incompressible Navier-Stokes equations in Cartesian domains containing immersed boundaries of arbitrary geometrical complexity moving with prescribed kinematics. The governing equations are discretized on a hybrid staggered/non-staggered grid layout using second-order accurate finite-difference formulas. The discrete equations are integrated in time via a second-order accurate dual-time-stepping, artificial compressibility iteration scheme. Unstructured, triangular meshes are employed to discretize complex immersed boundaries. The nodes of the surface mesh constitute a set of Lagrangian control points used to track the motion of the flexible body. At every instant in time, the influence of the body on the flow is accounted for by applying boundary conditions at Cartesian grid nodes located in the exterior but in the immediate vicinity of the body by reconstructing the solution along the local normal to the body surface. Grid convergence tests are carried out for the flow induced by an oscillating sphere in a cubic cavity, which show that the method is second-order accurate. The method is validated by applying it to calculate flow in a Cartesian domain containing a rigid sphere rotating at constant angular velocity as well as flow induced by a flapping wing. The ability of the method to simulate flows in domains with arbitrarily complex moving bodies is demonstrated by applying to simulate flow past an undulating fish-like body and flow past an anatomically realistic planktonic copepod performing an escape-like maneuver.
A linguistic geometry for 3D strategic planning
NASA Technical Reports Server (NTRS)
Stilman, Boris
1995-01-01
This paper is a new step in the development and application of the Linguistic Geometry. This formal theory is intended to discover the inner properties of human expert heuristics, which have been successful in a certain class of complex control systems, and apply them to different systems. In this paper we investigate heuristics extracted in the form of hierarchical networks of planning paths of autonomous agents. Employing Linguistic Geometry tools the dynamic hierarchy of networks is represented as a hierarchy of formal attribute languages. The main ideas of this methodology are shown in this paper on the new pilot example of the solution of the extremely complex 3D optimization problem of strategic planning for the space combat of autonomous vehicles. This example demonstrates deep and highly selective search in comparison with conventional search algorithms.
Cartesian grid simulations of gas-solids flow systems with complex geometry
Dietiker, Jean-Francois; Li, Tingwen; Garg, Rahul; Shahnam, Mehrdad
2013-02-01
Complex geometries encountered in many applications of gas–solids flow need special treatment in most legacy multiphase flow solvers with Cartesian numerical grid. This paper briefly outlines the implementation of a cut cell technique in the open-source multiphase flow solver—MFIX for accurate representation of complex geometries. Specifically, applications of the Cartesian cut cell method to different gas–solids fluidization systems including a small scale bubbling fluidized bed with submerged tube bundle and a complete pilot-scale circulating fluidized bed will be presented. In addition to qualitative predictions on the general flow behaviors inside each system, quantitative comparison with the available experimental data will be presented. Furthermore, some results on extending the current cut-cell technique to Lagrangian–Eulerian simulations will be presented.
Stochastic Modeling of Calcium in 3D Geometry
Mazel, Tomáš; Raymond, Rebecca; Raymond-Stintz, Mary; Jett, Stephen; Wilson, Bridget S.
2009-01-01
Release of inflammatory mediators by mast cells in type 1 immediate-hypersensitivity allergic reactions relies on antigen-dependent increases in cytosolic calcium. Here, we used a series of electron microscopy images to build a 3D reconstruction representing a slice through a rat tumor mast cell, which then served as a basis for stochastic modeling of inositol-trisphosphate-mediated calcium responses. The stochastic approach was verified by reaction-diffusion modeling within the same geometry. Local proximity of the endoplasmic reticulum to either the plasma membrane or mitochondria is predicted to differentially impact local inositol trisphosphate receptor transport. The explicit consideration of organelle spatial relationships represents an important step toward building a comprehensive, realistic model of cellular calcium dynamics. PMID:19254531
Aerodynamic Design of Complex Configurations Using Cartesian Methods and CAD Geometry
NASA Technical Reports Server (NTRS)
Nemec, Marian; Aftosmis, Michael J.; Pulliam, Thomas H.
2003-01-01
The objective for this paper is to present the development of an optimization capability for the Cartesian inviscid-flow analysis package of Aftosmis et al. We evaluate and characterize the following modules within the new optimization framework: (1) A component-based geometry parameterization approach using a CAD solid representation and the CAPRI interface. (2) The use of Cartesian methods in the development Optimization techniques using a genetic algorithm. The discussion and investigations focus on several real world problems of the optimization process. We examine the architectural issues associated with the deployment of a CAD-based design approach in a heterogeneous parallel computing environment that contains both CAD workstations and dedicated compute nodes. In addition, we study the influence of noise on the performance of optimization techniques, and the overall efficiency of the optimization process for aerodynamic design of complex three-dimensional configurations. of automated optimization tools. rithm and a gradient-based algorithm.
Geometry optimization for peptides and proteins: comparison of Cartesian and internal coordinates.
Koslover, Elena F; Wales, David J
2007-12-21
We present the results of several benchmarks comparing the relative efficiency of different coordinate systems in optimizing polypeptide geometries. Cartesian, natural internal, and primitive internal coordinates are employed in quasi-Newton minimization for a variety of biomolecules. The peptides and proteins used in these benchmarks range in size from 16 to 999 residues. They vary in complexity from polyalanine helices to a beta-barrel enzyme. We find that the relative performance of the different coordinate systems depends on the parameters of the optimization method, the starting point for the optimization, and the size of the system studied. In general, internal coordinates were found to be advantageous for small peptides. For larger structures, Cartesians appear to be more efficient for empirical potentials where the energy and gradient can be evaluated relatively quickly compared to the cost of the coordinate transformations. PMID:18154373
Geometry optimization for peptides and proteins: Comparison of Cartesian and internal coordinates
NASA Astrophysics Data System (ADS)
Koslover, Elena F.; Wales, David J.
2007-12-01
We present the results of several benchmarks comparing the relative efficiency of different coordinate systems in optimizing polypeptide geometries. Cartesian, natural internal, and primitive internal coordinates are employed in quasi-Newton minimization for a variety of biomolecules. The peptides and proteins used in these benchmarks range in size from 16 to 999 residues. They vary in complexity from polyalanine helices to a β-barrel enzyme. We find that the relative performance of the different coordinate systems depends on the parameters of the optimization method, the starting point for the optimization, and the size of the system studied. In general, internal coordinates were found to be advantageous for small peptides. For larger structures, Cartesians appear to be more efficient for empirical potentials where the energy and gradient can be evaluated relatively quickly compared to the cost of the coordinate transformations.
3D stochastic geophysical inversion for contact surface geometry
NASA Astrophysics Data System (ADS)
Lelièvre, Peter; Farquharson, Colin; Bijani, Rodrigo
2015-04-01
Geologists' interpretations about the Earth typically involve distinct rock units with contacts (interfaces) between them. As such, 3D geological Earth models typically comprise wireframe contact surfaces of tessellated triangles or other polygonal planar facets. In contrast, standard minimum-structure geophysical inversions are performed on meshes of space-filling cells (typically prisms or tetrahedra) and recover smoothly varying physical property distributions that are inconsistent with typical geological interpretations. There are several approaches through which mesh-based geophysical inversion can help recover models with some of the desired characteristics. However, a more effective strategy is to consider a fundamentally different type of inversion that works directly with models that comprise surfaces representing contacts between rock units. We are researching such an approach, our goal being to perform geophysical forward and inverse modelling directly with 3D geological models of any complexity. Geological and geophysical models should be specified using the same parameterization such that they are, in essence, the same Earth model. We parameterize the wireframe contact surfaces in a 3D model as the coordinates of the nodes (facet vertices). The physical properties of each rock unit in a model remain fixed while the geophysical inversion controls the position of the contact surfaces via the control nodes, perturbing the surfaces as required to fit the geophysical data responses. This is essentially a "geometry inversion", which can be used to recover the unknown geometry of a target body or to investigate the viability of a proposed Earth model. We apply global optimization strategies to solve the inverse problem, including stochastic sampling to obtain statistical information regarding the likelihood of particular features in the model, helping to assess the viability of a proposed model. Jointly inverting multiple types of geophysical data is simple
NASA Technical Reports Server (NTRS)
Aftosmis, M. J.; Berger, M. J.; Murman, S. M.; Kwak, Dochan (Technical Monitor)
2002-01-01
The proposed paper will present recent extensions in the development of an efficient Euler solver for adaptively-refined Cartesian meshes with embedded boundaries. The paper will focus on extensions of the basic method to include solution adaptation, time-dependent flow simulation, and arbitrary rigid domain motion. The parallel multilevel method makes use of on-the-fly parallel domain decomposition to achieve extremely good scalability on large numbers of processors, and is coupled with an automatic coarse mesh generation algorithm for efficient processing by a multigrid smoother. Numerical results are presented demonstrating parallel speed-ups of up to 435 on 512 processors. Solution-based adaptation may be keyed off truncation error estimates using tau-extrapolation or a variety of feature detection based refinement parameters. The multigrid method is extended to for time-dependent flows through the use of a dual-time approach. The extension to rigid domain motion uses an Arbitrary Lagrangian-Eulerlarian (ALE) formulation, and results will be presented for a variety of two- and three-dimensional example problems with both simple and complex geometry.
Alj, Domenico; Caputo, Roberto; Umeton, Cesare
2014-11-01
We report on the realization of a liquid crystal (LC)-based optical diffraction grating showing a polar symmetry of the director alignment. This has been obtained as a natural evolution of the POLICRYPS technique, which enables the realization of highly efficient, switchable, planar diffraction gratings. Performances exhibited in the Cartesian geometry are extended to the polar one by exploiting the spherical aberration produced by simple optical elements. This enables producing the required highly stable polar pattern that allows fabricating a circular optical diffraction grating. Results are promising for their possible application in fields in which a rotational structure of the optical beam is needed. PMID:25361314
The COMET method in 3-D hexagonal geometry
Connolly, K. J.; Rahnema, F.
2012-07-01
The hybrid stochastic-deterministic coarse mesh radiation transport (COMET) method developed at Georgia Tech now solves reactor core problems in 3-D hexagonal geometry. In this paper, the method is used to solve three preliminary test problems designed to challenge the method with steep flux gradients, high leakage, and strong asymmetry and heterogeneity in the core. The test problems are composed of blocks taken from a high temperature test reactor benchmark problem. As the method is still in development, these problems and their results are strictly preliminary. Results are compared to whole core Monte Carlo reference solutions in order to verify the method. Relative errors are on the order of 50 pcm in core eigenvalue, and mean relative error in pin fission density calculations is less than 1% in these difficult test cores. The method requires the one-time pre-computation of a response expansion coefficient library, which may be compiled in a comparable amount of time to a single whole core Monte Carlo calculation. After the library has been computed, COMET may solve any number of core configurations on the order of an hour, representing a significant gain in efficiency over other methods for whole core transport calculations. (authors)
Indoor Modelling Benchmark for 3D Geometry Extraction
NASA Astrophysics Data System (ADS)
Thomson, C.; Boehm, J.
2014-06-01
A combination of faster, cheaper and more accurate hardware, more sophisticated software, and greater industry acceptance have all laid the foundations for an increased desire for accurate 3D parametric models of buildings. Pointclouds are the data source of choice currently with static terrestrial laser scanning the predominant tool for large, dense volume measurement. The current importance of pointclouds as the primary source of real world representation is endorsed by CAD software vendor acquisitions of pointcloud engines in 2011. Both the capture and modelling of indoor environments require great effort in time by the operator (and therefore cost). Automation is seen as a way to aid this by reducing the workload of the user and some commercial packages have appeared that provide automation to some degree. In the data capture phase, advances in indoor mobile mapping systems are speeding up the process, albeit currently with a reduction in accuracy. As a result this paper presents freely accessible pointcloud datasets of two typical areas of a building each captured with two different capture methods and each with an accurate wholly manually created model. These datasets are provided as a benchmark for the research community to gauge the performance and improvements of various techniques for indoor geometry extraction. With this in mind, non-proprietary, interoperable formats are provided such as E57 for the scans and IFC for the reference model. The datasets can be found at: http://indoor-bench.github.io/indoor-bench.
Adjoint Sensitivity Computations for an Embedded-Boundary Cartesian Mesh Method and CAD Geometry
NASA Technical Reports Server (NTRS)
Nemec, Marian; Aftosmis,Michael J.
2006-01-01
Cartesian-mesh methods are perhaps the most promising approach for addressing the issues of flow solution automation for aerodynamic design problems. In these methods, the discretization of the wetted surface is decoupled from that of the volume mesh. This not only enables fast and robust mesh generation for geometry of arbitrary complexity, but also facilitates access to geometry modeling and manipulation using parametric Computer-Aided Design (CAD) tools. Our goal is to combine the automation capabilities of Cartesian methods with an eficient computation of design sensitivities. We address this issue using the adjoint method, where the computational cost of the design sensitivities, or objective function gradients, is esseutially indepeudent of the number of design variables. In previous work, we presented an accurate and efficient algorithm for the solution of the adjoint Euler equations discretized on Cartesian meshes with embedded, cut-cell boundaries. Novel aspects of the algorithm included the computation of surface shape sensitivities for triangulations based on parametric-CAD models and the linearization of the coupling between the surface triangulation and the cut-cells. The objective of the present work is to extend our adjoint formulation to problems involving general shape changes. Central to this development is the computation of volume-mesh sensitivities to obtain a reliable approximation of the objective finction gradient. Motivated by the success of mesh-perturbation schemes commonly used in body-fitted unstructured formulations, we propose an approach based on a local linearization of a mesh-perturbation scheme similar to the spring analogy. This approach circumvents most of the difficulties that arise due to non-smooth changes in the cut-cell layer as the boundary shape evolves and provides a consistent approximation tot he exact gradient of the discretized abjective function. A detailed gradient accurace study is presented to verify our approach
Types of Reasoning in 3D Geometry Thinking and Their Relation with Spatial Ability
ERIC Educational Resources Information Center
Pittalis, Marios; Christou, Constantinos
2010-01-01
The aim of this study is to describe and analyse the structure of 3D geometry thinking by identifying different types of reasoning and to examine their relation with spatial ability. To achieve this goal, two tests were administered to students in grades 5 to 9. The results of the study showed that 3D geometry thinking could be described by four…
Software-based geometry operations for 3D computer graphics
NASA Astrophysics Data System (ADS)
Sima, Mihai; Iancu, Daniel; Glossner, John; Schulte, Michael; Mamidi, Suman
2006-02-01
In order to support a broad dynamic range and a high degree of precision, many of 3D renderings fundamental algorithms have been traditionally performed in floating-point. However, fixed-point data representation is preferable over floating-point representation in graphics applications on embedded devices where performance is of paramount importance, while the dynamic range and precision requirements are limited due to the small display sizes (current PDA's are 640 × 480 (VGA), while cell-phones are even smaller). In this paper we analyze the efficiency of a CORDIC-augmented Sandbridge processor when implementing a vertex processor in software using fixed-point arithmetic. A CORDIC-based solution for vertex processing exhibits a number of advantages over classical Multiply-and-Acumulate solutions. First, since a single primitive is used to describe the computation, the code can easily be vectorized and multithreaded, and thus fits the major Sandbridge architectural features. Second, since a CORDIC iteration consists of only a shift operation followed by an addition, the computation may be deeply pipelined. Initially, we outline the Sandbridge architecture extension which encompasses a CORDIC functional unit and the associated instructions. Then, we consider rigid-body rotation, lighting, exponentiation, vector normalization, and perspective division (which are some of the most important data-intensive 3D graphics kernels) and propose a scheme to implement them on the CORDIC-augmented Sandbridge processor. Preliminary results indicate that the performance improvement within the extended instruction set ranges from 3× to 10× (with the exception of rigid body rotation).
Aerodynamic Optimization of Rocket Control Surface Geometry Using Cartesian Methods and CAD Geometry
NASA Technical Reports Server (NTRS)
Nelson, Andrea; Aftosmis, Michael J.; Nemec, Marian; Pulliam, Thomas H.
2004-01-01
Aerodynamic design is an iterative process involving geometry manipulation and complex computational analysis subject to physical constraints and aerodynamic objectives. A design cycle consists of first establishing the performance of a baseline design, which is usually created with low-fidelity engineering tools, and then progressively optimizing the design to maximize its performance. Optimization techniques have evolved from relying exclusively on designer intuition and insight in traditional trial and error methods, to sophisticated local and global search methods. Recent attempts at automating the search through a large design space with formal optimization methods include both database driven and direct evaluation schemes. Databases are being used in conjunction with surrogate and neural network models as a basis on which to run optimization algorithms. Optimization algorithms are also being driven by the direct evaluation of objectives and constraints using high-fidelity simulations. Surrogate methods use data points obtained from simulations, and possibly gradients evaluated at the data points, to create mathematical approximations of a database. Neural network models work in a similar fashion, using a number of high-fidelity database calculations as training iterations to create a database model. Optimal designs are obtained by coupling an optimization algorithm to the database model. Evaluation of the current best design then gives either a new local optima and/or increases the fidelity of the approximation model for the next iteration. Surrogate methods have also been developed that iterate on the selection of data points to decrease the uncertainty of the approximation model prior to searching for an optimal design. The database approximation models for each of these cases, however, become computationally expensive with increase in dimensionality. Thus the method of using optimization algorithms to search a database model becomes problematic as the
Reentry-Vehicle Shape Optimization Using a Cartesian Adjoint Method and CAD Geometry
NASA Technical Reports Server (NTRS)
Nemec, Marian; Aftosmis, Michael J.
2006-01-01
A DJOINT solutions of the governing flow equations are becoming increasingly important for the development of efficient analysis and optimization algorithms. A well-known use of the adjoint method is gradient-based shape. Given an objective function that defines some measure of performance, such as the lift and drag functionals, its gradient is computed at a cost that is essentially independent of the number of design variables (e.g., geometric parameters that control the shape). Classic aerodynamic applications of gradient-based optimization include the design of cruise configurations for transonic and supersonic flow, as well as the design of high-lift systems. are perhaps the most promising approach for addressing the issues of flow solution automation for aerodynamic design problems. In these methods, the discretization of the wetted surface is decoupled from that of the volume mesh. This not only enables fast and robust mesh generation for geometry of arbitrary complexity, but also facilitates access to geometry modeling and manipulation using parametric computer-aided design (CAD). In previous work on Cartesian adjoint solvers, Melvin et al. developed an adjoint formulation for the TRANAIR code, which is based on the full-potential equation with viscous corrections. More recently, Dadone and Grossman presented an adjoint formulation for the two-dimensional Euler equations using a ghost-cell method to enforce the wall boundary conditions. In Refs. 18 and 19, we presented an accurate and efficient algorithm for the solution of the adjoint Euler equations discretized on Cartesian meshes with embedded, cut-cell boundaries. Novel aspects of the algorithm were the computation of surface shape sensitivities for triangulations based on parametric-CAD models and the linearization of the coupling between the surface triangulation and the cut-cells. The accuracy of the gradient computation was verified using several three-dimensional test cases, which included design
ERIC Educational Resources Information Center
Goodson-Espy, Tracy; Lynch-Davis, Kathleen; Schram, Pamela; Quickenton, Art
2010-01-01
This paper describes the genesis and purpose of our geometry methods course, focusing on a geometry-teaching technology we created using NVIDIA[R] Chameleon demonstration. This article presents examples from a sequence of lessons centered about a 3D computer graphics demonstration of the chameleon and its geometry. In addition, we present data…
Verification of internal flow analyses in complex 3-D geometries
NASA Astrophysics Data System (ADS)
Choi, S. K.; Buggeln, R. C.
1992-11-01
Analysis of internal flow in advanced rocket propulsion systems is complicated by hardware geometry, high Reynolds numbers, rotation, high frequency phenomena, and near incompressibility. Typical of such a problem is the Space Shuttle Main Engine (SSME) hot gas manifold (HGM). Previous analyses of flow in the SSME HGM have been compared to air flow data and found to be inaccurate with respect to system losses, outer wall static pressures, and transfer duct environments. Such discrepancies could arise from flow measurement methodology, low order algorithms, turbulence modeling, and/or inadequate grid resolution. The objective of this work is to compare internal flow computational analyses to LDV flow measurements for the MSFC HGM pilot model configuration using two grids of different node density in the near wall region. Grids were generated with the EAGLE grid generator and calculations were made with the SRA MINT code. The calculated results were compared with HGM experimental data obtained in the MSFC water flow facility.
Development of a new two-dimensional Cartesian geometry nodal multigroup discrete-ordinates method
Pevey, R.E.
1982-07-01
The purpose of this work is the development and testing of a new family of methods for calculating the spatial dependence of the neutron density in nuclear systems described in two-dimensional Cartesian geometry. The energy and angular dependence of the neutron density is approximated using the multigroup and discrete ordinates techniques, respectively. The resulting FORTRAN computer code is designed to handle an arbitrary number of spatial, energy, and angle subdivisions. Any degree of scattering anisotropy can be handled by the code for either external source or fission systems. The basic approach is to (1) approximate the spatial variation of the neutron source across each spatial subdivision as an expansion in terms of a user-supplied set of exponential basis functions; (2) solve analytically for the resulting neutron density inside each region; and (3) approximate this density in the basis function space in order to calculate the next iteration flux-dependent source terms. In the general case the calculation is iterative due to neutron sources which depend on the neutron density itself, such as scattering interactions.
Super Cooled Large Droplet Analysis of Several Geometries Using LEWICE3D Version 3
NASA Technical Reports Server (NTRS)
Bidwell, Colin S.
2011-01-01
Super Cooled Large Droplet (SLD) collection efficiency calculations were performed for several geometries using the LEWICE3D Version 3 software. The computations were performed using the NASA Glenn Research Center SLD splashing model which has been incorporated into the LEWICE3D Version 3 software. Comparisons to experiment were made where available. The geometries included two straight wings, a swept 64A008 wing tip, two high lift geometries, and the generic commercial transport DLR-F4 wing body configuration. In general the LEWICE3D Version 3 computations compared well with the 2D LEWICE 3.2.2 results and with experimental data where available.
High-resolution, real-time simultaneous 3D surface geometry and temperature measurement.
An, Yatong; Zhang, Song
2016-06-27
This paper presents a method to simultaneously measure three-dimensional (3D) surface geometry and temperature in real time. Specifically, we developed 1) a holistic approach to calibrate both a structured light system and a thermal camera under exactly the same world coordinate system even though these two sensors do not share the same wavelength; and 2) a computational framework to determine the sub-pixel corresponding temperature for each 3D point as well as discard those occluded points. Since the thermal 2D imaging and 3D visible imaging systems do not share the same spectrum of light, they can perform sensing simultaneously in real time: we developed a hardware system that can achieve real-time 3D geometry and temperature measurement at 26 Hz with 768 × 960 points per frame. PMID:27410608
KENO3D Visualization Tool for KENO V.a and KENO-VI Geometry Models
Horwedel, J.E.; Bowman, S.M.
2000-06-01
Criticality safety analyses often require detailed modeling of complex geometries. Effective visualization tools can enhance checking the accuracy of these models. This report describes the KENO3D visualization tool developed at the Oak Ridge National Laboratory (ORNL) to provide visualization of KENO V.a and KENO-VI criticality safety models. The development of KENO3D is part of the current efforts to enhance the SCALE (Standardized Computer Analyses for Licensing Evaluations) computer software system.
NASA Astrophysics Data System (ADS)
Polkowski, Marcin
2016-04-01
Seismic wave travel time calculation is the most common numerical operation in seismology. The most efficient is travel time calculation in 1D velocity model - for given source, receiver depths and angular distance time is calculated within fraction of a second. Unfortunately, in most cases 1D is not enough to encounter differentiating local and regional structures. Whenever possible travel time through 3D velocity model has to be calculated. It can be achieved using ray calculation or time propagation in space. While single ray path calculation is quick it is complicated to find the ray path that connects source with the receiver. Time propagation in space using Fast Marching Method seems more efficient in most cases, especially when there are multiple receivers. In this presentation a Python module pySeismicFMM is presented - simple and very efficient tool for calculating travel time from sources to receivers. Calculation requires regular 2D or 3D velocity grid either in Cartesian or geographic coordinates. On desktop class computer calculation speed is 200k grid cells per second. Calculation has to be performed once for every source location and provides travel time to all receivers. pySeismicFMM is free and open source. Development of this tool is a part of authors PhD thesis. National Science Centre Poland provided financial support for this work via NCN grant DEC-2011/02/A/ST10/00284.
Conformal geometry and its applications on 3D shape matching, recognition, and stitching.
Wang, Sen; Wang, Yang; Jin, Miao; Gu, Xianfeng David; Samaras, Dimitris
2007-07-01
Three-dimensional shape matching is a fundamental issue in computer vision with many applications such as shape registration, 3D object recognition, and classification. However, shape matching with noise, occlusion, and clutter is a challenging problem. In this paper, we analyze a family of quasi-conformal maps including harmonic maps, conformal maps, and least-squares conformal maps with regards to 3D shape matching. As a result, we propose a novel and computationally efficient shape matching framework by using least-squares conformal maps. According to conformal geometry theory, each 3D surface with disk topology can be mapped to a 2D domain through a global optimization and the resulting map is a diffeomorphism, i.e., one-to-one and onto. This allows us to simplify the 3D shape-matching problem to a 2D image-matching problem, by comparing the resulting 2D parametric maps, which are stable, insensitive to resolution changes and robust to occlusion, and noise. Therefore, highly accurate and efficient 3D shape matching algorithms can be achieved by using the above three parametric maps. Finally, the robustness of least-squares conformal maps is evaluated and analyzed comprehensively in 3D shape matching with occlusion, noise, and resolution variation. In order to further demonstrate the performance of our proposed method, we also conduct a series of experiments on two computer vision applications, i.e., 3D face recognition and 3D nonrigid surface alignment and stitching. PMID:17496378
Spray Coating of Photoresist for 3D Microstructures with Different Geometries
NASA Astrophysics Data System (ADS)
Yu, Liming; Yeow Lee, Yong; Tay, Francis E. H.; Iliescu, Ciprian
2006-04-01
This paper presents the advantages of spray coating technique as compared to the conventional spin coating method for photoresist coating of 3D microstructures. An optimized mix of photoresist AZ4620: MEK: PGMEA (1:1.5:0.5) was used to achieve good coverage and uniformity of photoresist not only on planar surface, but also along the trenches' sidewall. In order to achieve the ideal coverage of photoresist layer, the effects of the geometries of the microstructures were also considered. Then, we implement this technique for our application in a MEMS device to prove the viability and potentiality of spray coating of photoresist for fabrication of 3D microstructures.
NASA Astrophysics Data System (ADS)
Zapiór, Maciej; Martínez-Gómez, David
2016-02-01
Based on the data collected by the Vacuum Tower Telescope located in the Teide Observatory in the Canary Islands, we analyzed the three-dimensional (3D) motion of so-called knots in a solar prominence of 2014 June 9. Trajectories of seven knots were reconstructed, giving information of the 3D geometry of the magnetic field. Helical motion was detected. From the equipartition principle, we estimated the lower limit of the magnetic field in the prominence to ≈1-3 G and from the Ampère’s law the lower limit of the electric current to ≈1.2 × 109 A.
Experiments with Uas Imagery for Automatic Modeling of Power Line 3d Geometry
NASA Astrophysics Data System (ADS)
Jóźków, G.; Vander Jagt, B.; Toth, C.
2015-08-01
The ideal mapping technology for transmission line inspection is the airborne LiDAR executed from helicopter platforms. It allows for full 3D geometry extraction in highly automated manner. Large scale aerial images can be also used for this purpose, however, automation is possible only for finding transmission line positions (2D geometry), and the sag needs to be estimated manually. For longer lines, these techniques are less expensive than ground surveys, yet they are still expensive. UAS technology has the potential to reduce these costs, especially if using inexpensive platforms with consumer grade cameras. This study investigates the potential of using high resolution UAS imagery for automatic modeling of transmission line 3D geometry. The key point of this experiment was to employ dense matching algorithms to appropriately acquired UAS images to have points created also on wires. This allowed to model the 3D geometry of transmission lines similarly to LiDAR acquired point clouds. Results showed that the transmission line modeling is possible with a high internal accuracy for both, horizontal and vertical directions, even when wires were represented by a partial (sparse) point cloud.
3D FEM Geometry and Material Flow Optimization of Porthole-Die Extrusion
Ceretti, Elisabetta; Mazzoni, Luca; Giardini, Claudio
2007-05-17
The aim of this work is to design and to improve the geometry of a porthole-die for the production of aluminum components by means of 3D FEM simulations. In fact, the use of finite element models will allow to investigate the effects of the die geometry (webs, extrusion cavity) on the material flow and on the stresses acting on the die so to reduce the die wear and to improve the tool life. The software used to perform the simulations was a commercial FEM code, Deform 3D. The technological data introduced in the FE model have been furnished by METRA S.p.A. Company, partner in this research. The results obtained have been considered valid and helpful by the Company for building a new optimized extrusion porthole-die.
Effect of geometry on drug release from 3D printed tablets.
Goyanes, Alvaro; Robles Martinez, Pamela; Buanz, Asma; Basit, Abdul W; Gaisford, Simon
2015-10-30
The aim of this work was to explore the feasibility of combining hot melt extrusion (HME) with 3D printing (3DP) technology, with a view to producing different shaped tablets which would be otherwise difficult to produce using traditional methods. A filament extruder was used to obtain approx. 4% paracetamol loaded filaments of polyvinyl alcohol with characteristics suitable for use in fused-deposition modelling 3DP. Five different tablet geometries were successfully 3D-printed-cube, pyramid, cylinder, sphere and torus. The printing process did not affect the stability of the drug. Drug release from the tablets was not dependent on the surface area but instead on surface area to volume ratio, indicating the influence that geometrical shape has on drug release. An erosion-mediated process controlled drug release. This work has demonstrated the potential of 3DP to manufacture tablet shapes of different geometries, many of which would be challenging to manufacture by powder compaction. PMID:25934428
Simulations of Coulomb systems with slab geometry using an efficient 3D Ewald summation method
NASA Astrophysics Data System (ADS)
dos Santos, Alexandre P.; Girotto, Matheus; Levin, Yan
2016-04-01
We present a new approach to efficiently simulate electrolytes confined between infinite charged walls using a 3d Ewald summation method. The optimal performance is achieved by separating the electrostatic potential produced by the charged walls from the electrostatic potential of electrolyte. The electric field produced by the 3d periodic images of the walls is constant inside the simulation cell, with the field produced by the transverse images of the charged plates canceling out. The non-neutral confined electrolyte in an external potential can be simulated using 3d Ewald summation with a suitable renormalization of the electrostatic energy, to remove a divergence, and a correction that accounts for the conditional convergence of the resulting lattice sum. The new algorithm is at least an order of magnitude more rapid than the usual simulation methods for the slab geometry and can be further sped up by adopting a particle-particle particle-mesh approach.
Laboratory Study of Magnetic Reconnection in 3D Geometry Relevant to Magnetopause and Magnetotail
NASA Astrophysics Data System (ADS)
Ren, Y.; Lu, Q.; Ji, H.; Mao, A.; Wang, X.; E, P.; Wang, Z.; Xiao, Q.; Ding, W.; Zheng, J.
2015-12-01
Laboratory Study of Magnetic Reconnection in 3D Geometry Relevant to Magnetopause and Magnetotail Y. Ren1,2, Quaming Lu3, Hantao Ji1,2, Aohua Mao1, Xiaogang Wang1, Peng E1, Zhibin Wang1, Qingmei Xiao1, Weixing Ding4, Jinxing Zheng51 Harbin Institute of Technology, Harbin, China2 Princeton Plasma Physics Laboratory, Princeton University, Princeton, NJ 08543 3University of Science and Technology of China, Hefei, China 4University of California at Los Angeles, Los Angeles, CA, 90095 5ASIPP, Hefei, China A new magnetic reconnection experiment, Harbin reconnection eXperiment (HRX), is currently being designed as a key part of Space Plasma Environment Research Facility (SPERF) at Harbin Institute of Technology in Harbin, China. HRX aims to provide a unique experimental platform for studying reconnections in 3D geometry relevant to magnetopause and magnetotail to address: the role of electron and ion-scale dynamics in the current sheet; particle and energy transfer from magnetosheath to magnetosphere; particle energization/heating mechanisms during magnetic reconnection; 3D effects in fast reconnection, e.g. the role of 3D magnetic null point. HRX employs a unique set of coils to generate the required 3D magnetic geometry and provides a wide range of plasma parameters. Here, important motivating scientific problems are reviewed and the physics design of HRX is presented, including plasma parameters determined from Vlasov scaling law, reconnection scenarios explored using vacuum magnetic field calculations and numerical simulations of HRX using hybrid and MHD codes. Plasma diagnostics plan and engineering design of important coils will also be briefly presented.
NASA Astrophysics Data System (ADS)
Mekuria, Rufael; Cesar, Pablo; Doumanis, Ioannis; Frisiello, Antonella
2015-09-01
Compression of 3D object based video is relevant for 3D Immersive applications. Nevertheless, the perceptual aspects of the degradation introduced by codecs for meshes and point clouds are not well understood. In this paper we evaluate the subjective and objective degradations introduced by such codecs in a state of art 3D immersive virtual room. In the 3D immersive virtual room, users are captured with multiple cameras, and their surfaces are reconstructed as photorealistic colored/textured 3D meshes or point clouds. To test the perceptual effect of compression and transmission, we render degraded versions with different frame rates in different contexts (near/far) in the scene. A quantitative subjective study with 16 users shows that negligible distortion of decoded surfaces compared to the original reconstructions can be achieved in the 3D virtual room. In addition, a qualitative task based analysis in a full prototype field trial shows increased presence, emotion, user and state recognition of the reconstructed 3D Human representation compared to animated computer avatars.
The 3-D world modeling with updating capability based on combinatorial geometry
NASA Technical Reports Server (NTRS)
Goldstein, M.; Pin, F. G.; Desaussure, G.; Weisbin, C. R.
1987-01-01
A 3-D world modeling technique using range data is discribed. Range data quantify the distances from the sensor focal plane to the object surface, i.e., the 3-D coordinates of discrete points on the object surface are known. The approach proposed herein for 3-D world modeling is based on the Combinatorial Geometry (CG) method which is widely used in Monte Carlo particle transport calculations. First, each measured point on the object surface is surrounded by a small sphere with a radius determined by the range to that point. Then, the 3-D shapes of the visible surfaces are obtained by taking the (Boolean) union of all the spheres. The result is an unambiguous representation of the object's boundary surfaces. The pre-learned partial knowledge of the environment can be also represented using the CG Method with a relatively small amount of data. Using the CG type of representation, distances in desired directions to boundary surfaces of various objects are efficiently calculated. This feature is particularly useful for continuously verifying the world model against the data provided by a range finder, and for integrating range data from successive locations of the robot during motion. The efficiency of the proposed approach is illustrated by simulations of a spherical robot in a 3-D room in the presence of moving obstacles and inadequate prelearned partial knowledge of the environment.
Xu, Gang; Xing, Mengdao; Xia, Xiang-Gen; Zhang, Lei; Chen, Qianqian; Bao, Zheng
2016-05-01
In the current scenario of high-resolution inverse synthetic aperture radar (ISAR) imaging, the non-cooperative targets may have strong maneuverability, which tends to cause time-variant Doppler modulation and imaging plane in the echoed data. Furthermore, it is still a challenge to realize ISAR imaging of maneuvering targets from sparse aperture (SA) data. In this paper, we focus on the problem of 3D geometry and motion estimations of maneuvering targets for interferometric ISAR (InISAR) with SA. For a target of uniformly accelerated rotation, the rotational modulation in echo is formulated as chirp sensing code under a chirp-Fourier dictionary to represent the maneuverability. In particular, a joint multi-channel imaging approach is developed to incorporate the multi-channel data and treat the multi-channel ISAR image formation as a joint-sparsity constraint optimization. Then, a modified orthogonal matching pursuit (OMP) algorithm is employed to solve the optimization problem to produce high-resolution range-Doppler (RD) images and chirp parameter estimation. The 3D target geometry and the motion estimations are followed by using the acquired RD images and chirp parameters. Herein, a joint estimation approach of 3D geometry and rotation motion is presented to realize outlier removing and error reduction. In comparison with independent single-channel processing, the proposed joint multi-channel imaging approach performs better in 2D imaging, 3D imaging, and motion estimation. Finally, experiments using both simulated and measured data are performed to confirm the effectiveness of the proposed algorithm. PMID:26930684
Mapping 3-D functional capillary geometry in rat skeletal muscle in vivo
Milkovich, Stephanie; Goldman, Daniel; Ellis, Christopher G.
2012-01-01
We have developed a novel mapping software package to reconstruct microvascular networks in three dimensions (3-D) from in vivo video images for use in blood flow and O2 transport modeling. An intravital optical imaging system was used to collect video sequences of blood flow in microvessels at different depths in the tissue. Functional images of vessels were produced from the video sequences and were processed using automated edge tracking software to yield location and geometry data for construction of the 3-D network. The same video sequences were analyzed for hemodynamic and O2 saturation data from individual capillaries in the network. Simple user-driven commands allowed the connection of vessel segments at bifurcations, and semiautomated registration enabled the tracking of vessels across multiple focal planes and fields of view. The reconstructed networks can be rotated and manipulated in 3-D to verify vessel connections and continuity. Hemodynamic and O2 saturation measurements made in vivo can be indexed to corresponding vessels and visualized using colorized maps of the vascular geometry. Vessels in each reconstruction are saved as text-based files that can be easily imported into flow or O2 transport models with complete geometry, hemodynamic, and O2 transport conditions. The results of digital morphometric analysis of seven microvascular networks showed mean capillary diameters and overall capillary density consistent with previous findings using histology and corrosion cast techniques. The described mapping software is a valuable tool for the quantification of in vivo microvascular geometry, hemodynamics, and oxygenation, thus providing rich data sets for experiment-based computational models. PMID:22140042
3-D geometry and physical property of the Mega-Splay Fault in Nankai trough
NASA Astrophysics Data System (ADS)
Masui, R.; Tsuji, T.; Yamada, Y.; Environmental Resource; System Engineering laboratory
2011-12-01
The Nankai trough is a subduction zone, where the Philippine Sea plate is being subducted beneath southwest Japan at a rate of ~4-6.5 cm/y at an azimuth of ~300°-315°. A lot of operations have been done in Nankai, such as three-dimensional seismic reflection surveys and Deep Sea Drilling Project (DSDP), Ocean Drilling Program (ODP), Integrated Ocean Drilling Program (IODP). They revealed that there is a large splay fault, referred to as 'Mega-Splay'. The Mega-Splay Fault has caused a series of catastrophic earthquakes and submarine landslides, which may have led to TSUNAMI. Since fault development history may have affected the geometry of the Mega-Splay Fault and physical property within the fault zone, they need to be examined in detail. In this research, we used 3-D pre-stack depth migration (PSDM), 3-D pre-stack time migration (PSTM) and P-wave velocity in C0004B well (Logging data), in order to interpret 3-D structure of Mega-Splay Fault. The analysis in this research is basically divided into two parts. One is structural interpretation of Splay Fault, based on the high amplitude reflection surface on seismic profiles. The other part is acoustic impedance inversion (AI inversion), in which we inverted seismic waveform into physical property (in this study, acoustic impedance), with the P-wave velocity data at C0004B near Mega-Splay Fault. The 3-D PSDM (or PSTM) clearly images details of Splay Fault, with good continuity of reflections along the fault. It is possible on each seismic profile to trace the high amplitude lines, where rock-properties significantly change. Since Mega-Splay Fault has 45-59m width along the wells, we interpreted the upper limit and the lower limit of the Mega-Splay Fault, based on the high amplitude surfaces along 3-D PSDM. Our interpretation shows that the width of Mega-Splay Fault has variation along the fault, and the plan geometry of the fault toe has a salient at the middle of the 3D box area, suggesting the fault could be
Li, Yong Gang; Yang, Yang; Short, Michael P.; Ding, Ze Jun; Zeng, Zhi; Li, Ju
2015-01-01
SRIM-like codes have limitations in describing general 3D geometries, for modeling radiation displacements and damage in nanostructured materials. A universal, computationally efficient and massively parallel 3D Monte Carlo code, IM3D, has been developed with excellent parallel scaling performance. IM3D is based on fast indexing of scattering integrals and the SRIM stopping power database, and allows the user a choice of Constructive Solid Geometry (CSG) or Finite Element Triangle Mesh (FETM) method for constructing 3D shapes and microstructures. For 2D films and multilayers, IM3D perfectly reproduces SRIM results, and can be ∼102 times faster in serial execution and > 104 times faster using parallel computation. For 3D problems, it provides a fast approach for analyzing the spatial distributions of primary displacements and defect generation under ion irradiation. Herein we also provide a detailed discussion of our open-source collision cascade physics engine, revealing the true meaning and limitations of the “Quick Kinchin-Pease” and “Full Cascades” options. The issues of femtosecond to picosecond timescales in defining displacement versus damage, the limitation of the displacements per atom (DPA) unit in quantifying radiation damage (such as inadequacy in quantifying degree of chemical mixing), are discussed. PMID:26658477
NASA Astrophysics Data System (ADS)
Li, Yong Gang; Yang, Yang; Short, Michael P.; Ding, Ze Jun; Zeng, Zhi; Li, Ju
2015-12-01
SRIM-like codes have limitations in describing general 3D geometries, for modeling radiation displacements and damage in nanostructured materials. A universal, computationally efficient and massively parallel 3D Monte Carlo code, IM3D, has been developed with excellent parallel scaling performance. IM3D is based on fast indexing of scattering integrals and the SRIM stopping power database, and allows the user a choice of Constructive Solid Geometry (CSG) or Finite Element Triangle Mesh (FETM) method for constructing 3D shapes and microstructures. For 2D films and multilayers, IM3D perfectly reproduces SRIM results, and can be ∼102 times faster in serial execution and > 104 times faster using parallel computation. For 3D problems, it provides a fast approach for analyzing the spatial distributions of primary displacements and defect generation under ion irradiation. Herein we also provide a detailed discussion of our open-source collision cascade physics engine, revealing the true meaning and limitations of the “Quick Kinchin-Pease” and “Full Cascades” options. The issues of femtosecond to picosecond timescales in defining displacement versus damage, the limitation of the displacements per atom (DPA) unit in quantifying radiation damage (such as inadequacy in quantifying degree of chemical mixing), are discussed.
RV functional imaging: 3-D echo-derived dynamic geometry and flow field simulations.
Pasipoularides, Ares D; Shu, Ming; Womack, Michael S; Shah, Ashish; Von Ramm, Olaf; Glower, Donald D
2003-01-01
We describe a novel functional imaging approach for quantitative analysis of right ventricular (RV) blood flow patterns in specific experimental animals (or humans) using real-time, three-dimensional (3-D) echocardiography (RT3D). The method is independent of the digital imaging modality used. It comprises three parts. First, a semiautomated segmentation aided by intraluminal contrast medium locates the RV endocardial surface. Second, a geometric scheme for dynamic RV chamber reconstruction applies a time interpolation procedure to the RT3D data to quantify wall geometry and motion at 400 Hz. A volumetric prism method validated the dynamic geometric reconstruction against simultaneous sonomicrometric canine measurements. Finally, the RV endocardial border motion information is used for mesh generation on a computational fluid dynamics solver to simulate development of the early RV diastolic inflow field. Boundary conditions (tessellated endocardial surface nodal velocities) for the solver are directly derived from the endocardial geometry and motion information. The new functional imaging approach may yield important kinematic information on the distribution of instantaneous velocities in the RV diastolic flow field of specific normal or diseased hearts. PMID:12388220
NASA Astrophysics Data System (ADS)
Giri, Abhra; Tarafdar, Sujata; Gouze, Philippe; Dutta, Tapati
2013-03-01
Several studies, both theoretical and experimental, show that sedimentary rocks have a fractal pore-grain interface. In this paper a computer simulated 3-D sedimentary rock structure generated by the Relaxed Ballistic Bidisperse Deposition Model (RBBDM), is investigated to characterize the micro structure of its pores. The pore volume and the rock-pore interface show the same fractal dimension indicating that the pore volume is a fractal. The two point density correlation is computed for the pore space and the results compare favourably with the range reported from experiments. An array of 2-D X-ray tomography micrograph sections of a real sedimentary rock, an oolitic limestone (pure calcite) from the Mondeville formation of Middle Jurassic age (Paris Basin, France), was used to generate a 3-D bitmap. The 3-D real rock sample generated in this manner, was analysed for similar studies as the simulated structure. The results were compared with those obtained from simulation. The simulation results agree qualitatively with the real rock sample. Diffusion through the connected pore space of the simulated structure was studied using a random walk algorithm and the results compared with the similar simulation study done on the 3-D oolitic limestone specimen. In both cases diffusion was found to be anomalous indicating that the sedimentary rock has a fractal geometry. The favourable comparability of results between the simulated and real rock supports the usefulness of the model of sedimentary rock generation which can be applicable to transport phenomena.
Color constancy using 3D scene geometry derived from a single image.
Elfiky, Noha; Gevers, Theo; Gijsenij, Arjan; Gonzalez, Jordi
2014-09-01
The aim of color constancy is to remove the effect of the color of the light source. As color constancy is inherently an ill-posed problem, most of the existing color constancy algorithms are based on specific imaging assumptions (e.g., gray-world and white patch assumption). In this paper, 3D geometry models are used to determine which color constancy method to use for the different geometrical regions (depth/layer) found in images. The aim is to classify images into stages (rough 3D geometry models). According to stage models, images are divided into stage regions using hard and soft segmentation. After that, the best color constancy methods are selected for each geometry depth. To this end, we propose a method to combine color constancy algorithms by investigating the relation between depth, local image statistics, and color constancy. Image statistics are then exploited per depth to select the proper color constancy method. Our approach opens the possibility to estimate multiple illuminations by distinguishing nearby light source from distant illuminations. Experiments on state-of-the-art data sets show that the proposed algorithm outperforms state-of-the-art single color constancy algorithms with an improvement of almost 50% of median angular error. When using a perfect classifier (i.e, all of the test images are correctly classified into stages); the performance of the proposed method achieves an improvement of 52% of the median angular error compared with the best-performing single color constancy algorithm. PMID:25051548
Disentangling Fault Scarp Geometry and Slip-Distribution in 3D
NASA Astrophysics Data System (ADS)
Mackenzie, D.; Walker, R. T.
2015-12-01
We present a new and inherently 3D approach to the analysis of fault scarp geometry using high resolution topography. Recent advance in topographic measurement techniques (LiDAR and Structure from Motion) has allowed the extensive measurement of single earthquake scarps and multiple event cumulative scarps to draw conclusions about along-strike slip variation and characteristic slip. Present analysis of the resulting point clouds and digital elevation models is generally achieved by taking vertical or map view profiles of geomorphic markers across the scarp. Profiles are done at numerous locations along strike carefully chosen to avoid regions degraded by erosion/deposition. The resulting slip distributions are almost always extremely variable and "noisy", both for strike-slip and dip-slip faulting scarps and it is often unclear whether this reflects slip variation, noise/erosion, site effects or geometric variation. When observing palaeo-earthquake and even modern event scarps, the full geometry, such as the degree of oblique slip or the fault dip, is often poorly constrained. We first present the results of synthetic tests to demonstrate the introduction of significant apparent noise by simply varying terrain, fault and measurement geometry (slope angle, slope azimuth, fault dip and slip obliquity). Considering fully 3-dimensional marker surfaces (e.g. Planar or conical) we use the variation in apparent offset with terrain and measurement geometry, to constrain the slip geometry in 3D. Combining measurements windowed along strike, we show that determining the slip vector is reduced to a simple linear problem. We conclude that for scarps in regions of significant topography or with oblique slip, our method will give enhanced slip resolution while standard methods will give poor slip resolution. We test our method using a Structure from Motion pointcloud and digital elevation model covering a ~25 km stretch of a thrust fault scarp in the Kazakh Tien Shan.
3D Digital Surveying and Modelling of Cave Geometry: Application to Paleolithic Rock Art
González-Aguilera, Diego; Muñoz-Nieto, Angel; Gómez-Lahoz, Javier; Herrero-Pascual, Jesus; Gutierrez-Alonso, Gabriel
2009-01-01
3D digital surveying and modelling of cave geometry represents a relevant approach for research, management and preservation of our cultural and geological legacy. In this paper, a multi-sensor approach based on a terrestrial laser scanner, a high-resolution digital camera and a total station is presented. Two emblematic caves of Paleolithic human occupation and situated in northern Spain, “Las Caldas” and “Peña de Candamo”, have been chosen to put in practise this approach. As a result, an integral and multi-scalable 3D model is generated which may allow other scientists, pre-historians, geologists…, to work on two different levels, integrating different Paleolithic Art datasets: (1) a basic level based on the accurate and metric support provided by the laser scanner; and (2) a advanced level using the range and image-based modelling. PMID:22399958
Subramanian, K R; Thubrikar, M J; Fowler, B; Mostafavi, M T; Funk, M W
2000-01-01
We present a technique that accurately reconstructs complex three dimensional blood vessel geometry from 2D intravascular ultrasound (IVUS) images. Biplane x-ray fluoroscopy is used to image the ultrasound catheter tip at a few key points along its path as the catheter is pulled through the blood vessel. An interpolating spline describes the continuous catheter path. The IVUS images are located orthogonal to the path, resulting in a non-uniform structured scalar volume of echo densities. Isocontour surfaces are used to view the vessel geometry, while transparency and clipping enable interactive exploration of interior structures. The two geometries studied are a bovine artery vascular graft having U-shape and a constriction, and a canine carotid artery having multiple branches and a constriction. Accuracy of the reconstructions is established by comparing the reconstructions to (1) silicone moulds of the vessel interior, (2) biplane x-ray images, and (3) the original echo images. Excellent shape and geometry correspondence was observed in both geometries. Quantitative measurements made at key locations of the 3D reconstructions also were in good agreement with those made in silicone moulds. The proposed technique is easily adoptable in clinical practice, since it uses x-rays with minimal exposure and existing IVUS technology. PMID:11105284
CasimirSim - A Tool to Compute Casimir Polder Forces for Nontrivial 3D Geometries
Sedmik, Rene; Tajmar, Martin
2007-01-30
The so-called Casimir effect is one of the most interesting macro-quantum effects. Being negligible on the macro-scale it becomes a governing factor below structure sizes of 1 {mu}m where it accounts for typically 100 kN m-2. The force does not depend on gravity, or electric charge but solely on the materials properties, and geometrical shape. This makes the effect a strong candidate for micro(nano)-mechanical devices M(N)EMS. Despite a long history of research the theory lacks a uniform description valid for arbitrary geometries which retards technical application. We present an advanced state-of-the-art numerical tool overcoming all the usual geometrical restrictions, capable of calculating arbitrary 3D geometries by utilizing the Casimir Polder approximation for the Casimir force.
Waveform tomography in 2.5-D to appropriately handle 3-D geometry
NASA Astrophysics Data System (ADS)
Smithyman, B.; Clowes, R. M.
2011-12-01
In order to improve the tractability of waveform tomography when applied to field data acquired along a crooked-line, we implement 2.5-D forward modeling and inversion. Waveform tomography combines conventional velocity-model building (i.e. tomography) with full-waveform inversion to reconstruct an image of subsurface acoustic velocity. For reasons of computational efficiency, it is desirable to use 2-D full-waveform inversion when processing data acquired with 2-D seismic survey geometry. However, crooked-line acquisition results in a cross-line component of the source-receiver offset that cannot be accounted for by 2-D forward modeling. If the cross-line geometry components are significant, full-waveform inversion may be intractable. To address the latter difficulty, we first apply 3-D traveltime tomography to generate a 2-D cross-sectional initial velocity model by taking a representative average slice through the 3-D model. Then this initial model from traveltime inversion is iteratively updated by 2.5-D full-waveform inversion using a frequency-domain viscoacoustic implementation. The 2.5-D method generates waveform data by combining the solutions of multiple 2-D wave equation components. Each wavefield represents the solution of a modified wave equation in which the cross-line wavenumber takes a value between zero and ~ω/c. The results are combined by inverse Fourier transform in the cross-line coordinate. This produces a synthetic wavefield that is a solution to the 3-D viscoacoustic wave equation in a 2-D velocity model. Consequently, the 2.5-D synthetic wavefield better approximates seismic field data (including crooked-line geometry), when compared to a 2-D synthetic result. Cross-line source-receiver offsets can be accounted for by reconstructing the wavefield out-of-plane with respect to the source. The amplitude and phase of the wavefield are consistent with a 3-D solution in a model that is homogeneous in one direction. The 2-D model is ideally
Development of an algorithm to measure defect geometry using a 3D laser scanner
NASA Astrophysics Data System (ADS)
Kilambi, S.; Tipton, S. M.
2012-08-01
Current fatigue life prediction models for coiled tubing (CT) require accurate measurements of the defect geometry. Three-dimensional (3D) laser imaging has shown promise toward becoming a nondestructive, non-contacting method of surface defect characterization. Laser imaging provides a detailed photographic image of a flaw, in addition to a detailed 3D surface map from which its critical dimensions can be measured. This paper describes algorithms to determine defect characteristics, specifically depth, width, length and projected cross-sectional area. Curve-fitting methods were compared and implicit algebraic fits have higher probability of convergence compared to explicit geometric fits. Among the algebraic fits, the Taubin circle fit has the least error. The algorithm was able to extract the dimensions of the flaw geometry from the scanned data of CT to within a tolerance of about 0.127 mm, close to the tolerance specified for the laser scanner itself, compared to measurements made using traveling microscopes. The algorithm computes the projected surface area of the flaw, which could previously only be estimated from the dimension measurements and the assumptions made about cutter shape. Although shadows compromised the accuracy of the shape characterization, especially for deep and narrow flaws, the results indicate that the algorithm with laser scanner can be used for non-destructive evaluation of CT in the oil field industry. Further work is needed to improve accuracy, to eliminate shadow effects and to reduce radial deviation.
Image-based reconstruction of 3D myocardial infarct geometry for patient specific applications
NASA Astrophysics Data System (ADS)
Ukwatta, Eranga; Rajchl, Martin; White, James; Pashakhanloo, Farhad; Herzka, Daniel A.; McVeigh, Elliot; Lardo, Albert C.; Trayanova, Natalia; Vadakkumpadan, Fijoy
2015-03-01
Accurate reconstruction of the three-dimensional (3D) geometry of a myocardial infarct from two-dimensional (2D) multi-slice image sequences has important applications in the clinical evaluation and treatment of patients with ischemic cardiomyopathy. However, this reconstruction is challenging because the resolution of common clinical scans used to acquire infarct structure, such as short-axis, late-gadolinium enhanced cardiac magnetic resonance (LGE-CMR) images, is low, especially in the out-of-plane direction. In this study, we propose a novel technique to reconstruct the 3D infarct geometry from low resolution clinical images. Our methodology is based on a function called logarithm of odds (LogOdds), which allows the broader class of linear combinations in the LogOdds vector space as opposed to being limited to only a convex combination in the binary label space. To assess the efficacy of the method, we used high-resolution LGE-CMR images of 36 human hearts in vivo, and 3 canine hearts ex vivo. The infarct was manually segmented in each slice of the acquired images, and the manually segmented data were downsampled to clinical resolution. The developed method was then applied to the downsampled image slices, and the resulting reconstructions were compared with the manually segmented data. Several existing reconstruction techniques were also implemented, and compared with the proposed method. The results show that the LogOdds method significantly outperforms all the other tested methods in terms of region overlap.
3D Simulation of Velocity Profile of Turbulent Flow in Open Channel with Complex Geometry
NASA Astrophysics Data System (ADS)
Kamel, Benoumessad; Ilhem, Kriba; Ali, Fourar; Abdelbaki, Djebaili
Simulation of open channel flow or river flow presents unique challenge to numerical simulators, which is widely used in the applications of computational fluid dynamics. The prediction is extremely difficult because the flow in open channel is usually transient and turbulent, the geometry is irregular and curved, and the free-surface elevation is varying with time. The results from a 3D non-linear k- ɛ turbulence model are presented to investigate the flow structure, the velocity distribution and mass transport process in a meandering compound open channel and a straight open channel. The 3D numerical model for calculating flow is set up in cylinder coordinates in order to calculate the complex boundary channel. The finite volume method is used to disperse the governing equations and the SIMPLE algorithm is applied to acquire the coupling of velocity and pressure. The non-linear k- ɛ turbulent model has good useful value because of taking into account the anisotropy and not increasing the computational time. The main contributions of this study are developing a numerical method that can be applied to predict the flow in river bends with various bend curvatures and different width-depth ratios. This study demonstrates that the 3D non-linear k- ɛ turbulence model can be used for analyzing flow structures, the velocity distribution and pollutant transport in the complex boundary open channel, this model is applicable for real river and wetland problem.
SU-E-J-128: 3D Surface Reconstruction of a Patient Using Epipolar Geometry
Kotoku, J; Nakabayashi, S; Kumagai, S; Ishibashi, T; Kobayashi, T; Haga, A; Saotome, N; Arai, N
2014-06-01
Purpose: To obtain a 3D surface data of a patient in a non-invasive way can substantially reduce the effort for the registration of patient in radiation therapy. To achieve this goal, we introduced the multiple view stereo technique, which is known to be used in a 'photo tourism' on the internet. Methods: 70 Images were taken with a digital single-lens reflex camera from different angles and positions. The camera positions and angles were inferred later in the reconstruction step. A sparse 3D reconstruction model was locating by SIFT features, which is robust for rotation and shift variance, in each image. We then found a set of correspondences between pairs of images by computing the fundamental matrix using the eight-point algorithm with RANSAC. After the pair matching, we optimized the parameter including camera positions to minimize the reprojection error by use of bundle adjustment technique (non-linear optimization). As a final step, we performed dense reconstruction and associate a color with each point using the library of PMVS. Results: Surface data were reconstructed well by visual inspection. The human skin is reconstructed well, althogh the reconstruction was time-consuming for direct use in daily clinical practice. Conclusion: 3D reconstruction using multi view stereo geometry is a promising tool for reducing the effort of patient setup. This work was supported by JSPS KAKENHI(25861128)
MHD Modeling in Complex 3D Geometries: Towards Predictive Simulation of SIHI Current Drive
NASA Astrophysics Data System (ADS)
Hansen, Christopher James
The HIT-SI experiment studies Steady Inductive Helicity Injection (SIHI) for the purpose of forming and sustaining a spheromak plasma. A spheromak is formed in a nearly axisymmetric flux conserver, with a bow tie cross section, by means of two semi-toroidal injectors. The plasma-facing surfaces of the device, which are made of copper for its low resistivity, are covered in an insulating coating in order to operate in a purely inductive manner. Following formation, the spheromak flux and current are increased during a quiescent period marked by a decrease in the global mode activity. A proposed mechanism, Imposed Dynamo Current Drive (IDCD), is expected to be responsible for this phase of quiescent current drive. Due to the geometric complexity of the experiment, previous numerical modeling efforts have used a simplified geometry that excludes the injector volumes from the simulated domain. The effect of helicity injection is then modeled by boundary conditions on this reduced plasma volume. The work presented here has explored and developed more complete computational models of the HIT-SI device. This work is separated into 3 distinct but complementary areas: 1) Development of a 3D MHD equilibrium code that can incorporate the non-axisymmetric injector fields present in HIT-SI and investigation of equilibria of interest during spheromak sustainment. 2) A 2D axisymmetric MHD equilibrium code that was used to explore reduced order models for mean-field evolution using equations derived from IDCD theory including coupling to 3D equilibria. 3) A 3D time-dependent non-linear MHD code that is capable of modeling the entire plasma volume including dynamics within the injectors. Although HIT-SI was the motivation for, and experiment studied in this research, the tools and methods developed are general --- allowing their application to a broad range of magnetic confinement experiments. These tools constitute a significant advance for modeling plasma dynamics in devices with
Amoeboid migration mode adaption in quasi-3D spatial density gradients of varying lattice geometry
NASA Astrophysics Data System (ADS)
Gorelashvili, Mari; Emmert, Martin; Hodeck, Kai F.; Heinrich, Doris
2014-07-01
Cell migration processes are controlled by sensitive interaction with external cues such as topographic structures of the cell’s environment. Here, we present systematically controlled assays to investigate the specific effects of spatial density and local geometry of topographic structure on amoeboid migration of Dictyostelium discoideum cells. This is realized by well-controlled fabrication of quasi-3D pillar fields exhibiting a systematic variation of inter-pillar distance and pillar lattice geometry. By time-resolved local mean-squared displacement analysis of amoeboid migration, we can extract motility parameters in order to elucidate the details of amoeboid migration mechanisms and consolidate them in a two-state contact-controlled motility model, distinguishing directed and random phases. Specifically, we find that directed pillar-to-pillar runs are found preferably in high pillar density regions, and cells in directed motion states sense pillars as attractive topographic stimuli. In contrast, cell motion in random probing states is inhibited by high pillar density, where pillars act as obstacles for cell motion. In a gradient spatial density, these mechanisms lead to topographic guidance of cells, with a general trend towards a regime of inter-pillar spacing close to the cell diameter. In locally anisotropic pillar environments, cell migration is often found to be damped due to competing attraction by different pillars in close proximity and due to lack of other potential stimuli in the vicinity of the cell. Further, we demonstrate topographic cell guidance reflecting the lattice geometry of the quasi-3D environment by distinct preferences in migration direction. Our findings allow to specifically control amoeboid cell migration by purely topographic effects and thus, to induce active cell guidance. These tools hold prospects for medical applications like improved wound treatment, or invasion assays for immune cells.
NASA Astrophysics Data System (ADS)
Wada, I.; Wang, K.; He, J.
2013-12-01
In this study, we revisit the effects of along-strike variation in slab geometry and oblique subduction on subduction zone thermal structures. Along-strike variations in slab dip cause changes in the descending rate of the slab and generate trench-parallel pressure gradients that drive trench-parallel mantle flow (e.g., Kneller and van Keken, 2007). Oblique subduction also drives trench-parallel mantle flow. In this study, we use a finite element code PGCtherm3D and examine a range of generic subduction geometries and parameters to investigate the effects of the above two factors. This exercise is part of foundational work towards developing detailed 3-D thermal models for NE Japan, Nankai, and Cascadia to better constrain their 3-D thermal structures and to understand the role of temperature in controlling metamorphic, seismogenic, and volcanic processes. The 3-D geometry of the subducting slabs in the forearc and arc regions are well delineated at these three subduction zones. Further, relatively large compilations of surface heat flow data at these subduction zones make them excellent candidates for this study. At NE Japan, a megathrust earthquake occurred on March 11, 2011; at Nankai and Cascadia, there has been a great effort to constrain the scale of the next subduction thrust earthquake for the purpose of disaster prevention. Temperature influences the slip behavior of subduction faults by (1) affecting the rheology of the interface material and (2) controlling dehydration reactions, which can lead to elevated pore fluid pressure. Beyond the depths of subduction thrust earthquakes, the thermal structure is affected strongly by the pattern of mantle wedge flow. This flow is driven by viscous coupling between the subducting slab and the overriding mantle, and it brings in hot flowing mantle into the wedge. The trench-ward (up-dip) extent of the slab-mantle coupling is thus a key factor that controls the thermal structure. Slab-mantle decoupling at shallow
Relative Scale Estimation and 3D Registration of Multi-Modal Geometry Using Growing Least Squares.
Mellado, Nicolas; Dellepiane, Matteo; Scopigno, Roberto
2016-09-01
The advent of low cost scanning devices and the improvement of multi-view stereo techniques have made the acquisition of 3D geometry ubiquitous. Data gathered from different devices, however, result in large variations in detail, scale, and coverage. Registration of such data is essential before visualizing, comparing and archiving them. However, state-of-the-art methods for geometry registration cannot be directly applied due to intrinsic differences between the models, e.g., sampling, scale, noise. In this paper we present a method for the automatic registration of multi-modal geometric data, i.e., acquired by devices with different properties (e.g., resolution, noise, data scaling). The method uses a descriptor based on Growing Least Squares, and is robust to noise, variation in sampling density, details, and enables scale-invariant matching. It allows not only the measurement of the similarity between the geometry surrounding two points, but also the estimation of their relative scale. As it is computed locally, it can be used to analyze large point clouds composed of millions of points. We implemented our approach in two registration procedures (assisted and automatic) and applied them successfully on a number of synthetic and real cases. We show that using our method, multi-modal models can be automatically registered, regardless of their differences in noise, detail, scale, and unknown relative coverage. PMID:26672045
A formal classification of 3D medial axis points and their local geometry.
Giblin, Peter; Kimia, Benjamin B
2004-02-01
This paper proposes a novel hypergraph skeletal representation for 3D shape based on a formal derivation of the generic structure of its medial axis. By classifying each skeletal point by its order of contact, we show that, generically, the medial axis consists of five types of points, which are then organized into sheets, curves, and points: 1) sheets (manifolds with boundary) which are the locus of bitangent spheres with regular tangency A1(2) (Ak(n) notation means n distinct k-fold tangencies of the sphere of contact, as explained in the text); two types of curves, 2) the intersection curve of three sheets and the locus of centers of tritangent spheres, A1(3), and 3) the boundary of sheets, which are the locus of centers of spheres whose radius equals the larger principal curvature, i.e., higher order contact A3 points; and two types of points, 4) centers of quad-tangent spheres, A1(4), and 5) centers of spheres with one regular tangency and one higher order tangency, A1A3. The geometry of the 3D medial axis thus consists of sheets (A1(2)) bounded by one type of curve (A3) on their free end, which corresponds to ridges on the surface, and attached to two other sheets at another type of curve (A1(3)), which support a generalized cylinder description. The A3 curves can only end in A1A3 points where they must meet an A1(3) curve. The A1(3) curves meet together in fours at an A1(4) point. This formal result leads to a compact representation for 3D shape, referred to as the medial axis hypergraph representation consisting of nodes (A1(4) and A1A3 points), links between pairs of nodes (A1(3) and A3 curves) and hyperlinks between groups of links (A1(2) sheets). The description of the local geometry at nodes by itself is sufficient to capture qualitative aspects of shapes, in analogy to 2D. We derive a pointwise reconstruction formula to reconstruct a surface from this medial axis hypergraph together with the radius function. Thus, this information completely
NASA Astrophysics Data System (ADS)
Cai, Tao
2016-04-01
In this paper, we have described a 'stratified' semi-implicit spectral method to study compressible convection in Cartesian geometry. The full set of compressible hydrodynamic equations are solved in conservative forms. The numerical scheme is accurate and efficient, based on fast Fourier/sin/cos spectral transforms in the horizontal directions, Chebyshev spectral transform or second-order finite difference scheme in the vertical direction, and second order semi-implicit scheme in time marching of linear terms. We have checked the validity of both the fully pseudo-spectral scheme and the mixed finite-difference pseudo-spectral scheme by studying the onset of compressible convection. The difference of the critical Rayleigh number between our numerical result and the linear stability analysis is within two percent. Besides, we have computed the Mach numbers with different Rayleigh numbers in compressible convection. It shows good agreement with the numerical results of finite difference methods and finite volume method. This model has wide application in studying laminar and turbulent flow. Illustrative examples of application on horizontal convection, gravity waves, and long-lived vortex are given in this paper.
PDE constrained optimization of electrical defibrillation in a 3D ventricular slice geometry.
Chamakuri, Nagaiah; Kunisch, Karl; Plank, Gernot
2016-04-01
A computational study of an optimal control approach for cardiac defibrillation in a 3D geometry is presented. The cardiac bioelectric activity at the tissue and bath volumes is modeled by the bidomain model equations. The model includes intramural fiber rotation, axially symmetric around the fiber direction, and anisotropic conductivity coefficients, which are extracted from a histological image. The dynamics of the ionic currents are based on the regularized Mitchell-Schaeffer model. The controls enter in the form of electrodes, which are placed at the boundary of the bath volume with the goal of dampening undesired arrhythmias. The numerical optimization is based on Newton techniques. We demonstrated the parallel architecture environment for the computation of potentials on multidomains and for the higher order optimization techniques. Copyright © 2015 John Wiley & Sons, Ltd. PMID:26249168
Effects of 3D geometries on cellular gradient sensing and polarization
NASA Astrophysics Data System (ADS)
Spill, Fabian; Andasari, Vivi; Mak, Michael; Kamm, Roger D.; Zaman, Muhammad H.
2016-06-01
During cell migration, cells become polarized, change their shape, and move in response to various internal and external cues. Cell polarization is defined through the spatio-temporal organization of molecules such as PI3K or small GTPases, and is determined by intracellular signaling networks. It results in directional forces through actin polymerization and myosin contractions. Many existing mathematical models of cell polarization are formulated in terms of reaction–diffusion systems of interacting molecules, and are often defined in one or two spatial dimensions. In this paper, we introduce a 3D reaction–diffusion model of interacting molecules in a single cell, and find that cell geometry has an important role affecting the capability of a cell to polarize, or change polarization when an external signal changes direction. Our results suggest a geometrical argument why more roundish cells can repolarize more effectively than cells which are elongated along the direction of the original stimulus, and thus enable roundish cells to turn faster, as has been observed in experiments. On the other hand, elongated cells preferentially polarize along their main axis even when a gradient stimulus appears from another direction. Furthermore, our 3D model can accurately capture the effect of binding and unbinding of important regulators of cell polarization to and from the cell membrane. This spatial separation of membrane and cytosol, not possible to capture in 1D or 2D models, leads to marked differences of our model from comparable lower-dimensional models.
NASA Astrophysics Data System (ADS)
Zhang, Y.; Chen, C.; Du, J.; Sun, S.; Liang, Q.
2015-12-01
In the study of the inversion of gravity and magnetic data, the discretization of underground space is usually achieved by the use of structured grids. For instance, using the regular block as the module unit to divide model space in Cartesian coordinate system and the tesseroid in spherical coordinate system. Structured grids show clear spatial structures and mathematical properties. However, the block can only provide a rough approximation to the given terrain and using the tesseroid to approximate the terrain even seems impracticable. These shape determining errors cause the reduction of forward modeling precision. Moreover, the precision decreases again while using the tesseroid as no analytical algorithm has been acquired. On the other hand, since most terrain data has a limited resolution, unstructured grids, based on the polyhedron or tetrahedron, could fill the space completely, which allows us to reduce errors in shape determination to the minima. In addition, the analytical algorithms for polyhedron have been proposed. In our study, we use the tetrahedron as the module unit to divide the underground space. Moreover, based on the former researches, we supplement new analytical algorithms for tetrahedron to forward modeling gravity and magnetic fields and their gradient tensors in both Cartesian and spherical coordinate systems. The algorithm is testified by comparing the forward gravity and magnetic data of a block with the data obtained using the existed algorithms. The absolute difference between these two data is under 10e-9 mGal. Our approach is suitable for the inversion of gravity and magnetic data in both Cartesian and spherical coordinate systems.This study is supported by Natural Science Fund of Hubei Province (Grant No.: 2015CFB361) and International Cooperation Project in Science and Technology of China (Grant No.: 2010DFA24580).
Lei Liu; Feng Zhou; Xue-Ru Bai; Ming-Liang Tao; Zi-Jing Zhang
2016-04-01
Traditionally, the factorization method is applied to reconstruct the 3D geometry of a target from its sequential inverse synthetic aperture radar images. However, this method requires performing cross-range scaling to all the sub-images and thus has a large computational burden. To tackle this problem, this paper proposes a novel method for joint cross-range scaling and 3D geometry reconstruction of steadily moving targets. In this method, we model the equivalent rotational angular velocity (RAV) by a linear polynomial with time, and set its coefficients randomly to perform sub-image cross-range scaling. Then, we generate the initial trajectory matrix of the scattering centers, and solve the 3D geometry and projection vectors by the factorization method with relaxed constraints. After that, the coefficients of the polynomial are estimated from the projection vectors to obtain the RAV. Finally, the trajectory matrix is re-scaled using the estimated rotational angle, and accurate 3D geometry is reconstructed. The two major steps, i.e., the cross-range scaling and the factorization, are performed repeatedly to achieve precise 3D geometry reconstruction. Simulation results have proved the effectiveness and robustness of the proposed method. PMID:26886991
Amundsen, Morten; Linder, Jacob
2016-01-01
An extension of quasiclassical Keldysh-Usadel theory to higher spatial dimensions than one is crucial in order to describe physical phenomena like charge/spin Hall effects and topological excitations like vortices and skyrmions, none of which are captured in one-dimensional models. We here present a numerical finite element method which solves the non-linearized 2D and 3D quasiclassical Usadel equation relevant for the diffusive regime. We show the application of this on three model systems with non-trivial geometries: (i) a bottlenecked Josephson junction with external flux, (ii) a nanodisk ferromagnet deposited on top of a superconductor and (iii) superconducting islands in contact with a ferromagnet. In case (i), we demonstrate that one may control externally not only the geometrical array in which superconducting vortices arrange themselves, but also to cause coalescence and tune the number of vortices. In case (iii), we show that the supercurrent path can be tailored by incorporating magnetic elements in planar Josephson junctions which also lead to a strong modulation of the density of states. The finite element method presented herein paves the way for gaining insight in physical phenomena which have remained largely unexplored due to the complexity of solving the full quasiclassical equations in higher dimensions. PMID:26961921
Amundsen, Morten; Linder, Jacob
2016-01-01
An extension of quasiclassical Keldysh-Usadel theory to higher spatial dimensions than one is crucial in order to describe physical phenomena like charge/spin Hall effects and topological excitations like vortices and skyrmions, none of which are captured in one-dimensional models. We here present a numerical finite element method which solves the non-linearized 2D and 3D quasiclassical Usadel equation relevant for the diffusive regime. We show the application of this on three model systems with non-trivial geometries: (i) a bottlenecked Josephson junction with external flux, (ii) a nanodisk ferromagnet deposited on top of a superconductor and (iii) superconducting islands in contact with a ferromagnet. In case (i), we demonstrate that one may control externally not only the geometrical array in which superconducting vortices arrange themselves, but also to cause coalescence and tune the number of vortices. In case (iii), we show that the supercurrent path can be tailored by incorporating magnetic elements in planar Josephson junctions which also lead to a strong modulation of the density of states. The finite element method presented herein paves the way for gaining insight in physical phenomena which have remained largely unexplored due to the complexity of solving the full quasiclassical equations in higher dimensions. PMID:26961921
3D CFD modeling of subsonic and transonic flowing-gas DPALs with different pumping geometries
NASA Astrophysics Data System (ADS)
Yacoby, Eyal; Sadot, Oren; Barmashenko, Boris D.; Rosenwaks, Salman
2015-10-01
Three-dimensional computational fluid dynamics (3D CFD) modeling of subsonic (Mach number M ~ 0.2) and transonic (M ~ 0.9) diode pumped alkali lasers (DPALs), taking into account fluid dynamics and kinetic processes in the lasing medium is reported. The performance of these lasers is compared with that of supersonic (M ~ 2.7 for Cs and M ~ 2.4 for K) DPALs. The motivation for this study stems from the fact that subsonic and transonic DPALs require much simpler hardware than supersonic ones where supersonic nozzle, diffuser and high power mechanical pump (due to a drop in the gas total pressure in the nozzle) are required for continuous closed cycle operation. For Cs DPALs with 5 x 5 cm2 flow cross section pumped by large cross section (5 x 2 cm2) beam the maximum achievable power of supersonic devices is higher than that of the transonic and subsonic devices by only ~ 3% and ~ 10%, respectively. Thus in this case the supersonic operation mode has no substantial advantage over the transonic one. The main processes limiting the power of Cs supersonic DPALs are saturation of the D2 transition and large ~ 60% losses of alkali atoms due to ionization, whereas the influence of gas heating is negligible. For K transonic DPALs both the gas heating and ionization effects are shown to be unimportant. The maximum values of the power are higher than those in Cs transonic laser by ~ 11%. The power achieved in the supersonic and transonic K DPAL is higher than for the subsonic version, with the same resonator and K density at the inlet, by ~ 84% and ~ 27%, respectively, showing a considerable advantaged of the supersonic device over the transonic one. For pumping by rectangular beams of the same (5 x 2 cm2) cross section, comparison between end-pumping - where the laser beam and pump beam both propagate at along the same axis, and transverse-pumping - where they propagate perpendicularly to each other, shows that the output power and optical-to-optical efficiency are not
Yu, Dongliang; Yin, Min; Lu, Linfeng; Zhang, Hanzhong; Chen, Xiaoyuan; Zhu, Xufei; Che, Jianfei; Li, Dongdong
2015-11-01
High-performance thin-film hydrogenated amorphous silicon solar cells are achieved by combining macroscale 3D tubular substrates and nanoscaled 3D cone-like antireflective films. The tubular geometry delivers a series of advantages for large-scale deployment of photovoltaics, such as omnidirectional performance, easier encapsulation, decreased wind resistance, and easy integration with a second device inside the glass tube. PMID:26418573
NASA Astrophysics Data System (ADS)
Yu, Peicheng; Xu, Xinlu; Davidson, Asher; Tableman, Adam; Dalichaouch, Thamine; Li, Fei; Meyers, Michael D.; An, Weiming; Tsung, Frank S.; Decyk, Viktor K.; Fiuza, Frederico; Vieira, Jorge; Fonseca, Ricardo A.; Lu, Wei; Silva, Luis O.; Mori, Warren B.
2016-07-01
When modeling laser wakefield acceleration (LWFA) using the particle-in-cell (PIC) algorithm in a Lorentz boosted frame, the plasma is drifting relativistically at βb c towards the laser, which can lead to a computational speedup of ∼ γb2 = (1 - βb2)-1. Meanwhile, when LWFA is modeled in the quasi-3D geometry in which the electromagnetic fields and current are decomposed into a limited number of azimuthal harmonics, speedups are achieved by modeling three dimensional (3D) problems with the computational loads on the order of two dimensional r - z simulations. Here, we describe a method to combine the speedups from the Lorentz boosted frame and quasi-3D algorithms. The key to the combination is the use of a hybrid Yee-FFT solver in the quasi-3D geometry that significantly mitigates the Numerical Cerenkov Instability (NCI) which inevitably arises in a Lorentz boosted frame due to the unphysical coupling of Langmuir modes and EM modes of the relativistically drifting plasma in these simulations. In addition, based on the space-time distribution of the LWFA data in the lab and boosted frame, we propose to use a moving window to follow the drifting plasma, instead of following the laser driver as is done in the LWFA lab frame simulations, in order to further reduce the computational loads. We describe the details of how the NCI is mitigated for the quasi-3D geometry, the setups for simulations which combine the Lorentz boosted frame, quasi-3D geometry, and the use of a moving window, and compare the results from these simulations against their corresponding lab frame cases. Good agreement is obtained among these sample simulations, particularly when there is no self-trapping, which demonstrates it is possible to combine the Lorentz boosted frame and the quasi-3D algorithms when modeling LWFA. We also discuss the preliminary speedups achieved in these sample simulations.
NASA Astrophysics Data System (ADS)
Capone, M.; Campi, M.; Catuogno, R.
2015-02-01
This paper is part of a research about ribbed vaults systems in French Gothic Cathedrals. Our goal is to compare some different gothic cathedrals to understand the complex geometry of the ribbed vaults. The survey isn't the main objective but it is the way to verify the theoretical hypotheses about geometric configuration of the flamboyant churches in Paris. The survey method's choice generally depends on the goal; in this case we had to study many churches in a short time, so we chose 3D reconstruction method based on image dense stereo matching. This method allowed us to obtain the necessary information to our study without bringing special equipment, such as the laser scanner. The goal of this paper is to test image matching 3D reconstruction method in relation to some particular study cases and to show the benefits and the troubles. From a methodological point of view this is our workflow: - theoretical study about geometrical configuration of rib vault systems; - 3D model based on theoretical hypothesis about geometric definition of the vaults' form; - 3D model based on image matching 3D reconstruction methods; - comparison between 3D theoretical model and 3D model based on image matching;
Energy Science and Technology Software Center (ESTSC)
2013-06-24
Version 07 TART2012 is a coupled neutron-photon Monte Carlo transport code designed to use three-dimensional (3-D) combinatorial geometry. Neutron and/or photon sources as well as neutron induced photon production can be tracked. It is a complete system to assist you with input preparation, running Monte Carlo calculations, and analysis of output results. TART2012 is also incredibly FAST; if you have used similar codes, you will be amazed at how fast this code is compared tomore » other similar codes. Use of the entire system can save you a great deal of time and energy. TART2012 extends the general utility of the code to even more areas of application than available in previous releases by concentrating on improving the physics, particularly with regard to improved treatment of neutron fission, resonance self-shielding, molecular binding, and extending input options used by the code. Several utilities are included for creating input files and displaying TART results and data. TART2012 uses the latest ENDF/B-VI, Release 8, data. New for TART2012 is the use of continuous energy neutron cross sections, in addition to its traditional multigroup cross sections. For neutron interaction, the data are derived using ENDF-ENDL2005 and include both continuous energy cross sections and 700 group neutron data derived using a combination of ENDF/B-VI, Release 8, and ENDL data. The 700 group structure extends from 10-5 eV up to 1 GeV. Presently nuclear data are only available up to 20 MeV, so that only 616 of the groups are currently used. For photon interaction, 701 point photon data were derived using the Livermore EPDL97 file. The new 701 point structure extends from 100 eV up to 1 GeV, and is currently used over this entire energy range. TART2012 completely supersedes all older versions of TART, and it is strongly recommended that one use only the most recent version of TART2012 and its data files. Check authors homepage for related information: http
Self-calibration of cone-beam CT geometry using 3D-2D image registration.
Ouadah, S; Stayman, J W; Gang, G J; Ehtiati, T; Siewerdsen, J H
2016-04-01
Robotic C-arms are capable of complex orbits that can increase field of view, reduce artifacts, improve image quality, and/or reduce dose; however, it can be challenging to obtain accurate, reproducible geometric calibration required for image reconstruction for such complex orbits. This work presents a method for geometric calibration for an arbitrary source-detector orbit by registering 2D projection data to a previously acquired 3D image. It also yields a method by which calibration of simple circular orbits can be improved. The registration uses a normalized gradient information similarity metric and the covariance matrix adaptation-evolution strategy optimizer for robustness against local minima and changes in image content. The resulting transformation provides a 'self-calibration' of system geometry. The algorithm was tested in phantom studies using both a cone-beam CT (CBCT) test-bench and a robotic C-arm (Artis Zeego, Siemens Healthcare) for circular and non-circular orbits. Self-calibration performance was evaluated in terms of the full-width at half-maximum (FWHM) of the point spread function in CBCT reconstructions, the reprojection error (RPE) of steel ball bearings placed on each phantom, and the overall quality and presence of artifacts in CBCT images. In all cases, self-calibration improved the FWHM-e.g. on the CBCT bench, FWHM = 0.86 mm for conventional calibration compared to 0.65 mm for self-calibration (p < 0.001). Similar improvements were measured in RPE-e.g. on the robotic C-arm, RPE = 0.73 mm for conventional calibration compared to 0.55 mm for self-calibration (p < 0.001). Visible improvement was evident in CBCT reconstructions using self-calibration, particularly about high-contrast, high-frequency objects (e.g. temporal bone air cells and a surgical needle). The results indicate that self-calibration can improve even upon systems with presumably accurate geometric calibration and is
Self-calibration of cone-beam CT geometry using 3D-2D image registration
NASA Astrophysics Data System (ADS)
Ouadah, S.; Stayman, J. W.; Gang, G. J.; Ehtiati, T.; Siewerdsen, J. H.
2016-04-01
Robotic C-arms are capable of complex orbits that can increase field of view, reduce artifacts, improve image quality, and/or reduce dose; however, it can be challenging to obtain accurate, reproducible geometric calibration required for image reconstruction for such complex orbits. This work presents a method for geometric calibration for an arbitrary source-detector orbit by registering 2D projection data to a previously acquired 3D image. It also yields a method by which calibration of simple circular orbits can be improved. The registration uses a normalized gradient information similarity metric and the covariance matrix adaptation-evolution strategy optimizer for robustness against local minima and changes in image content. The resulting transformation provides a ‘self-calibration’ of system geometry. The algorithm was tested in phantom studies using both a cone-beam CT (CBCT) test-bench and a robotic C-arm (Artis Zeego, Siemens Healthcare) for circular and non-circular orbits. Self-calibration performance was evaluated in terms of the full-width at half-maximum (FWHM) of the point spread function in CBCT reconstructions, the reprojection error (RPE) of steel ball bearings placed on each phantom, and the overall quality and presence of artifacts in CBCT images. In all cases, self-calibration improved the FWHM—e.g. on the CBCT bench, FWHM = 0.86 mm for conventional calibration compared to 0.65 mm for self-calibration (p < 0.001). Similar improvements were measured in RPE—e.g. on the robotic C-arm, RPE = 0.73 mm for conventional calibration compared to 0.55 mm for self-calibration (p < 0.001). Visible improvement was evident in CBCT reconstructions using self-calibration, particularly about high-contrast, high-frequency objects (e.g. temporal bone air cells and a surgical needle). The results indicate that self-calibration can improve even upon systems with presumably accurate geometric calibration and is
NASA Astrophysics Data System (ADS)
Wada, Ikuko; He, Jiangheng; Hasegawa, Akira; Nakajima, Junichi
2015-09-01
We develop a 3-D thermal model for the Northeast Japan subduction margin, using a realistic slab geometry for the subducting Pacific plate, and investigate the effects of oblique subduction and 3-D slab geometry on the mantle wedge flow pattern and the thermal structure. In the Tohoku region, the mantle wedge flow pattern is nearly two-dimensional resulting in a thermal structure similar to those obtained by a 2-D model, owing to the simple slab geometry and subduction nearly perpendicular to the margin. However, in Hokkaido, oblique subduction leads to 3-D mantle wedge flow with northerly inflow and west-northwestward outflow and also results in lower temperatures in the shallow part of the mantle wedge than in Tohoku due to lower sinking rate of the slab. Between Hokkaido and Tohoku, the slab has a hinge-like shape due to a relatively sharp change in the dip direction. In this hinge zone, northerly mantle inflow from Hokkaido and westerly mantle inflow from Tohoku converge, discouraging inflow from northwest and resulting in a cooler mantle wedge. The model-predicted mantle wedge flow patterns are consistent with observed seismic anisotropy and may explain the orientations of volcanic cross-chains. The predicted 3-D thermal structure correlates well with the along-arc variations in the location of the frontal arc volcanoes and help to provide new insights into the surface heat flow pattern and the down-dip extent of interplate earthquakes.
NASA Technical Reports Server (NTRS)
Colborn, B. L.; Armstrong, T. W.
1992-01-01
A computer model of the three dimensional geometry and material distributions for the LDEF spacecraft, experiment trays, and, for selected trays, the components of experiments within a tray was developed for use in ionizing radiation assessments. The model is being applied to provide 3-D shielding distributions around radiation dosimeters to aid in data interpretation, particularly in assessing the directional properties of the radiation exposure. Also, the model has been interfaced with radiation transport codes for 3-D dosimetry response predictions and for calculations related to determining the accuracy of trapped proton and cosmic ray environment models. The methodology is described used in developing the 3-D LDEF model and the level of detail incorporated. Currently, the trays modeled in detail are F2, F8, and H12 and H3. Applications of the model which are discussed include the 3-D shielding distributions around various dosimeters, the influence of shielding on dosimetry responses, and comparisons of dose predictions based on the present 3-D model vs those from 1-D geometry model approximations used in initial estimates.
Gai, Jiading; Obeid, Nady; Holtrop, Joseph L.; Wu, Xiao-Long; Lam, Fan; Fu, Maojing; Haldar, Justin P.; Hwu, Wen-mei W.; Liang, Zhi-Pei; Sutton, Bradley P.
2013-01-01
Several recent methods have been proposed to obtain significant speed-ups in MRI image reconstruction by leveraging the computational power of GPUs. Previously, we implemented a GPU-based image reconstruction technique called the Illinois Massively Parallel Acquisition Toolkit for Image reconstruction with ENhanced Throughput in MRI (IMPATIENT MRI) for reconstructing data collected along arbitrary 3D trajectories. In this paper, we improve IMPATIENT by removing computational bottlenecks by using a gridding approach to accelerate the computation of various data structures needed by the previous routine. Further, we enhance the routine with capabilities for off-resonance correction and multi-sensor parallel imaging reconstruction. Through implementation of optimized gridding into our iterative reconstruction scheme, speed-ups of more than a factor of 200 are provided in the improved GPU implementation compared to the previous accelerated GPU code. PMID:23682203
geomIO: A tool for geodynamicists to turn 2D cross-sections into 3D geometries
NASA Astrophysics Data System (ADS)
Baumann, Tobias; Bauville, Arthur
2016-04-01
In numerical deformation models, material properties are usually defined on elements (e.g., in body-fitted finite elements), or on a set of Lagrangian markers (Eulerian, ALE or mesh-free methods). In any case, geometrical constraints are needed to assign different material properties to the model domain. Whereas simple geometries such as spheres, layers or cuboids can easily be programmed, it quickly gets complex and time-consuming to create more complicated geometries for numerical model setups, especially in three dimensions. geomIO (geometry I/O, http://geomio.bitbucket.org/) is a MATLAB-based library that has two main functionalities. First, it can be used to create 3D volumes based on series of 2D vector drawings similar to a CAD program; and second, it uses these 3D volumes to assign material properties to the numerical model domain. The drawings can conveniently be created using the open-source vector graphics software Inkscape. Adobe Illustrator is also partially supported. The drawings represent a series of cross-sections in the 3D model domain, for example, cross-sectional interpretations of seismic tomography. geomIO is then used to read the drawings and to create 3D volumes by interpolating between the cross-sections. In the second part, the volumes are used to assign material phases to markers inside the volumes. Multiple volumes can be created at the same time and, depending on the order of assignment, unions or intersections can be built to assign additional material phases. geomIO also offers the possibility to create 3D temperature structures for geodynamic models based on depth dependent parameterisations, for example the half space cooling model. In particular, this can be applied to geometries of subducting slabs of arbitrary shape. Yet, geomIO is held very general, and can be used for a variety of applications. We present examples of setup generation from pictures of micro-scale tectonics and lithospheric scale setups of 3D present-day model
Ardekani, Siamak; Jain, Saurabh; Sanzi, Alianna; Corona-Villalobos, Celia P; Abraham, Theodore P; Abraham, M Roselle; Zimmerman, Stefan L; Wu, Katherine C; Winslow, Raimond L; Miller, Michael I; Younes, Laurent
2016-04-01
The focus of this study was to develop advanced mathematical tools to construct high-resolution 3D models of left-ventricular (LV) geometry to evaluate focal geometric differences between patients with hypertrophic cardiomyopathy (HCM) and hypertensive heart disease (HHD) using cardiac magnetic resonance (MR) cross-sectional images. A limiting factor in 3D analysis of cardiac MR cross-sections is the low out-of-plane resolution of the acquired images. To overcome this problem, we have developed a mathematical framework to construct a population-based high-resolution 3D LV triangulated surface (template) in which an iterative matching algorithm maps a surface mesh of a normal heart to a set of cross-sectional contours that were extracted from short-axis cine cardiac MR images of patients who were diagnosed with either HCM or HHD. A statistical analysis was conducted on deformations that were estimated at each surface node to identify shape differences at end-diastole (ED), end-systole (ES), and motion-related shape variation from ED to ES. Some significant shape difference in radial thickness was detected at ES. Differences of LV 3D surface geometry were identified focally on the basal anterior septum wall. Further research is needed to relate these findings to the HCM morphological substrate and to design a classifier to discriminate among different etiologies of LV hypertrophy. PMID:26766206
Effects of Training Method and Gender on Learning 2D/3D Geometry
ERIC Educational Resources Information Center
Khairulanuar, Samsudin; Nazre, Abd Rashid; Jamilah, H.; Sairabanu, Omar Khan; Norasikin, Fabil
2010-01-01
This article reports the findings of an experimental study involving 36 primary school students (16 girls, 20 boys, Mean age = 9.5 years, age range: 8-10 years) in geometrical understanding of 2D and 3D objects. Students were assigned into two experimental groups and one control group based on a stratified random sampling procedure. The first…
GSRP/David Marshall: Fully Automated Cartesian Grid CFD Application for MDO in High Speed Flows
NASA Technical Reports Server (NTRS)
2003-01-01
With the renewed interest in Cartesian gridding methodologies for the ease and speed of gridding complex geometries in addition to the simplicity of the control volumes used in the computations, it has become important to investigate ways of extending the existing Cartesian grid solver functionalities. This includes developing methods of modeling the viscous effects in order to utilize Cartesian grids solvers for accurate drag predictions and addressing the issues related to the distributed memory parallelization of Cartesian solvers. This research presents advances in two areas of interest in Cartesian grid solvers, viscous effects modeling and MPI parallelization. The development of viscous effects modeling using solely Cartesian grids has been hampered by the widely varying control volume sizes associated with the mesh refinement and the cut cells associated with the solid surface. This problem is being addressed by using physically based modeling techniques to update the state vectors of the cut cells and removing them from the finite volume integration scheme. This work is performed on a new Cartesian grid solver, NASCART-GT, with modifications to its cut cell functionality. The development of MPI parallelization addresses issues associated with utilizing Cartesian solvers on distributed memory parallel environments. This work is performed on an existing Cartesian grid solver, CART3D, with modifications to its parallelization methodology.
Geometry-based vs. intensity-based medical image registration: A comparative study on 3D CT data.
Savva, Antonis D; Economopoulos, Theodore L; Matsopoulos, George K
2016-02-01
Spatial alignment of Computed Tomography (CT) data sets is often required in numerous medical applications and it is usually achieved by applying conventional exhaustive registration techniques, which are mainly based on the intensity of the subject data sets. Those techniques consider the full range of data points composing the data, thus negatively affecting the required processing time. Alternatively, alignment can be performed using the correspondence of extracted data points from both sets. Moreover, various geometrical characteristics of those data points can be used, instead of their chromatic properties, for uniquely characterizing each point, by forming a specific geometrical descriptor. This paper presents a comparative study reviewing variations of geometry-based, descriptor-oriented registration techniques, as well as conventional, exhaustive, intensity-based methods for aligning three-dimensional (3D) CT data pairs. In this context, three general image registration frameworks were examined: a geometry-based methodology featuring three distinct geometrical descriptors, an intensity-based methodology using three different similarity metrics, as well as the commonly used Iterative Closest Point algorithm. All techniques were applied on a total of thirty 3D CT data pairs with both known and unknown initial spatial differences. After an extensive qualitative and quantitative assessment, it was concluded that the proposed geometry-based registration framework performed similarly to the examined exhaustive registration techniques. In addition, geometry-based methods dramatically improved processing time over conventional exhaustive registration. PMID:26771247
Mesoscopic hydrogel molding to control the 3D geometry of bioartificial muscle tissues
Bian, Weining; Liau, Brian; Badie, Nima
2010-01-01
This protocol describes a cell/hydrogel molding method for precise and reproducible biomimetic fabrication of three-dimensional (3D) muscle tissue architectures in vitro. Using a high aspect ratio soft lithography technique, we fabricate polydimethylsiloxane (PDMS) molds containing arrays of mesoscopic posts with defined size, elongation and spacing. On cell/hydrogel molding, these posts serve to enhance the diffusion of nutrients to cells by introducing elliptical pores in the cell-laden hydrogels and to guide local 3D cell alignment by governing the spatial pattern of mechanical tension. Instead of ultraviolet or chemical cross-linking, this method utilizes natural hydrogel polymerization and topographically constrained cell-mediated gel compaction to create the desired 3D tissue structures. We apply this method to fabricate several square centimeter large, few hundred micron-thick bioartificial muscle tissues composed of viable, dense, uniformly aligned and highly differentiated cardiac or skeletal muscle fibers. The protocol takes 4–5 d to fabricate PDMS molds followed by 2 weeks of cell culture. PMID:19798085
The Need (?) for Descriptive Geometry in a World of 3D Modeling.
ERIC Educational Resources Information Center
Croft, Frank M. Jr.
1998-01-01
Evaluates the use of modern CAD methods to solve geometric problems. Solves descriptive geometry problems using the layout and position of the successive auxiliary views from the projection of three-dimensional figures onto a two-dimensional plane of paper. (CCM)
Math in Action. 3-D Geometry Comes in All Shapes and Sizes.
ERIC Educational Resources Information Center
Bresser, Rusty; Sheffield, Stephanie
1996-01-01
Two mathematics activities introduce elementary students to geometry, using depth, width, and height. A primary level activity involves students with sorting and classifying geometric solids and with comparing numerical quantities. The activity requires building blocks or geometric solids. An intermediate level activity teaches students about…
The Emergence of 3D Geometry from Children's (Teacher-Guided) Classification Tasks
ERIC Educational Resources Information Center
Roth, Wolff-Michael; Thom, Jennifer S.
2009-01-01
Geometry, classification, and the classification of geometrical objects are integral aspects of recent curriculum documents in mathematics education. Such curriculum documents, however, leave open how the "work" of classifying objects according to geometrical properties can be accomplished given that the knowledge of these properties is the…
Kloxin, April M.; Lewis, Katherine J. R.; DeForest, Cole A.; Seedorf, Gregory; Tibbitt, Mark W.; Balasubramaniam, Vivek
2012-01-01
We describe the development of a well-based cell culture platform that enables experimenters to control the geometry and connectivity of cellular microenvironments spatiotemporally. The base material is a hydrogel comprised of photolabile and enzyme-labile crosslinks and pendant cell adhesion sequences, enabling spatially-specific, in situ patterning with light and cell-dictated microenvironment remodeling through enzyme secretion. Arrays of culture wells of varying shape and size were patterned into the hydrogel surface using photolithography, where well depth was correlated with irradiation dose. The geometry of these devices can be subsequently modified through sequential patterning, while simultaneously monitoring changes in cell geometry and connectivity. Towards establishing the utility of these devices for dynamic evaluation of the influence of physical cues on tissue morphogenesis, the effect of well shape on lung epithelial cell differentiation (i.e., primary mouse alveolar type II cells, ATII cells) was assessed. Shapes inspired by alveoli were degraded into hydrogel surfaces. ATII cells were seeded within the well-based arrays and encapsulated by the addition of a top hydrogel layer. Cell differentiation in response to these geometries was characterized over 7 days of culture with immunocytochemistry (surfactant protein C, ATII; T1α protein, alveolar type I (ATI) differentiated epithelial cells) and confocal image analysis. Individual cell clusters were further connected by eroding channels between wells during culture via controlled two-photon irradiation. Collectively, these studies demonstrate the development and utility of responsive hydrogel culture devices to study how a range of microenvironment geometries of evolving shape and connectivity might influence or direct cell function. PMID:23138879
NASA Astrophysics Data System (ADS)
Kincaid, C. R.; MacDougall, J. G.; Druken, K. A.; Fischer, K. M.
2010-12-01
Understanding patterns in plate scale mantle flow in subduction zones is key to models of thermal structure, dehydration reactions, volatile distributions and magma generation and transport in convergent margins. Different patterns of flow in the mantle wedge can generate distinct signatures in seismological observables. Observed shear wave fast polarization directions in several subduction zones are inconsistent with predictions of simple 2-D wedge corner flow. Geochemical signatures in a number of subduction zones also indicate 3-D flow and entrainment patterns in the wedge. We report on a series of laboratory experiments on subduction driven flow to characterize spatial and temporal variability in 3-D patterns in flow and shear-induced finite strain. Cases focus on how rollback subduction, along-strike dip changes in subducting plates and evolving gaps or tears in subduction zones control temporal-spatial patterns in 3-D wedge flow. Models utilize a glucose working fluid with a temperature dependent viscosity to represent the upper 2000 km of the mantle. Subducting lithosphere is modeled with two rubber-reinforced continuous belts. Belts pass around trench and upper/lower mantle rollers. The deeper rollers can move laterally to allow for time varying dip angle. Each belt has independent speed control and dip adjustment, allowing for along-strike changes in convergence rate and the evolution of slab gaps. Rollback is modeled using a translation system to produce either uniform and asymmetric lateral trench motion. Neutral density finite strain markers are distributed throughout the fluid and used as proxies for tracking the evolution of anisotropy through space and time in the evolving flow fields. Particle image velocimetry methods are also used to track time varying 3-D velocity fields for directly calculating anisotropy patterns. Results show that complex plate motions (rollback, steepening) and morphologies (gaps) in convergent margins produce flows with
3-D in vivo brain tumor geometry study by scaling analysis
NASA Astrophysics Data System (ADS)
Torres Hoyos, F.; Martín-Landrove, M.
2012-02-01
A new method, based on scaling analysis, is used to calculate fractal dimension and local roughness exponents to characterize in vivo 3-D tumor growth in the brain. Image acquisition was made according to the standard protocol used for brain radiotherapy and radiosurgery, i.e., axial, coronal and sagittal magnetic resonance T1-weighted images, and comprising the brain volume for image registration. Image segmentation was performed by the application of the k-means procedure upon contrasted images. We analyzed glioblastomas, astrocytomas, metastases and benign brain tumors. The results show significant variations of the parameters depending on the tumor stage and histological origin.
NASA Astrophysics Data System (ADS)
Schneeberger, Raphael; de la Varga, Miguel; Florian Wellmann, J.; Kober, Florian; Berger, Alfons; Herwegh, Marco
2016-04-01
Fluid circulation in crystalline rocks is of key importance when exploring crystalline basement in light of, for example, deep-seated geothermal energy projects or selection of sites for nuclear waste repositories. Due to their enhanced permeability, fluid circulation within crystalline bedrock is mainly controlled by fault zones, which may originate from ductile mylonites but show a strong brittle overprint. In order to better constrain 3D flow paths, a well-founded knowledge on the 3D nature of the fault zone pattern is indispensable. We attempt to constrain the geometry of a complex 3D fault zone pattern in a case study of the Grimsel Test Site (GTS, central Switzerland). The constraints are based on mapping of both the surface as well as the GTS underground tunnel system, offering a unique opportunity to test the 3D model and associated uncertainties. We investigate the effect of increasing geoinformation on the quality and accuracy of the 3D model by using: (i) remote sensing surface data only, (ii) field surface mapping in combination with (i), and (iii) underground data combined with (i) and (ii). This approach allows for defining different steps in 3D geological modelling of a specific area, including a measure of the remaining uncertainty after each step. We obtain a best-estimate model by fitting results between surface and underground data by using a combination of field data and orientation obtained by Delaunay triangulation. We incorporate novel approaches to uncertainty analysis of fault orientations and investigate different fault planes showing the possible variation range of the structures investigated.
The Monte Carlo SRNA-VOX code for 3D proton dose distribution in voxelized geometry using CT data
NASA Astrophysics Data System (ADS)
Ilic, Radovan D.; Spasic-Jokic, Vesna; Belicev, Petar; Dragovic, Milos
2005-03-01
This paper describes the application of the SRNA Monte Carlo package for proton transport simulations in complex geometry and different material compositions. The SRNA package was developed for 3D dose distribution calculation in proton therapy and dosimetry and it was based on the theory of multiple scattering. The decay of proton induced compound nuclei was simulated by the Russian MSDM model and our own using ICRU 63 data. The developed package consists of two codes: the SRNA-2KG, which simulates proton transport in combinatorial geometry and the SRNA-VOX, which uses the voxelized geometry using the CT data and conversion of the Hounsfield's data to tissue elemental composition. Transition probabilities for both codes are prepared by the SRNADAT code. The simulation of the proton beam characterization by multi-layer Faraday cup, spatial distribution of positron emitters obtained by the SRNA-2KG code and intercomparison of computational codes in radiation dosimetry, indicate immediate application of the Monte Carlo techniques in clinical practice. In this paper, we briefly present the physical model implemented in the SRNA package, the ISTAR proton dose planning software, as well as the results of the numerical experiments with proton beams to obtain 3D dose distribution in the eye and breast tumour.
The Monte Carlo SRNA-VOX code for 3D proton dose distribution in voxelized geometry using CT data.
Ilić, Radovan D; Spasić-Jokić, Vesna; Belicev, Petar; Dragović, Milos
2005-03-01
This paper describes the application of the SRNA Monte Carlo package for proton transport simulations in complex geometry and different material compositions. The SRNA package was developed for 3D dose distribution calculation in proton therapy and dosimetry and it was based on the theory of multiple scattering. The decay of proton induced compound nuclei was simulated by the Russian MSDM model and our own using ICRU 63 data. The developed package consists of two codes: the SRNA-2KG, which simulates proton transport in combinatorial geometry and the SRNA-VOX, which uses the voxelized geometry using the CT data and conversion of the Hounsfield's data to tissue elemental composition. Transition probabilities for both codes are prepared by the SRNADAT code. The simulation of the proton beam characterization by multi-layer Faraday cup, spatial distribution of positron emitters obtained by the SRNA-2KG code and intercomparison of computational codes in radiation dosimetry, indicate immediate application of the Monte Carlo techniques in clinical practice. In this paper, we briefly present the physical model implemented in the SRNA package, the ISTAR proton dose planning software, as well as the results of the numerical experiments with proton beams to obtain 3D dose distribution in the eye and breast tumour. PMID:15798273
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)
Facsko, Gabor; Reshetnyk, Volodymyr; Agapitov, Oleksiy; Opitz, Andrea; Szabo, Adam; McComas, David
Tangential discontinuties (TDs) are usually considered as thin planar current sheets frozen in the solar wind flow. Previous studies based on the magnetic field measurements onboard of ACE, Wind, and STEREO A, and B proved that this hypotesis is not valid. The curvature of the TDs were determined in several cases. After applying minimum variance and the cross product methods for Ulysses, ACE and STEREO A and B magnetometer measurements, numerous TDs are identified in 2008 and 2009. The time shift of the TD observations is determinated by correlation analysis of the solar wind speed and the magnetic field variations. The 3D topology of the TD is then determinated in some special cases when the four spacecraft are on the same side of the Sun. After fitting a simple model, the location of the TD formation region can be outlined.
The Geometry of the Subducting Slabs Beneath the PRVI Microplate Based on 3D Tomography
NASA Astrophysics Data System (ADS)
Xu, X.; Keller, G. R.; Holland, A. A.; Keranen, K. M.; Li, H.
2011-12-01
The Puerto Rico and the Virgin Islands (PRVI) microplate is located between two subduction zones, with the Puerto Rico trench to the north and the Muertos trough to the south. The Puerto Rico trench is caused by southward-directed subduction of the North American Plate, and the Muertos trough is the northern boundary of the Caribbean Plate. There is no active volcanism on Puerto Rico; however, earthquake depths and seismic tomography imply that the slab of Caribbean plate continues northward beneath Puerto Rico. Puerto Rico overlies these two slabs with extension to both the west (Mona Passage) and southeast (Anacapa Passage). The cause of the extension is unknown, but GPS measurements show that Puerto Rico is experiencing anti-clockwise rotation, and the extension associated with the Anacapa Passage may be produced by this rotation. To the west, it is debated whether the Mona Passage is a boundary between two micro-plates or simple a local rift basin. To address the sources of the extension and the cause of the rotation, we are investigating if the deep structures can be the dynamic source for the observed kinematic movements. We collected data on earthquakes occurring between 2009-2011 in the PRVI region and relocated them using the SEISAN code provided by the Institute of Solid Earth Physics, University of Bergen. The FMTOMO code from Australian National University was used for 3D tomography from P and S wave arrival times. By comparing the relocated epicenters and the 3D tomography results, the subducting slabs were identified. When integrated with the results of previous studies, the geometric model of the slabs is a critical key to understanding the evolution of the PRVI microplate in the past and the future.
Irregular geometries in normal unmyelinated axons: a 3D serial EM analysis.
Greenberg, M M; Leitao, C; Trogadis, J; Stevens, J K
1990-12-01
Axons have generally been represented as straight cylinders. It is not at all uncommon for anatomists to take single cross-sections of an axonal bundle, and from the axonal diameter compute expected conduction velocities. This assumes that each cross-section represents a slice through a perfect cylinder. We have examined the three-dimensional geometry of 98 central and peripheral unmyelinated axons, using computer-assisted serial electron microscopy. These reconstructions reveal that virtually all unmyelinated axons have highly irregular axial shapes consisting of periodic varicosities. The varicosities were, without exception, filled with membranous organelles frequently including mitochondria, and have obligatory volumes similar to that described in other neurites. The mitochondria make contact with microtubules, while the other membraneous organelles were frequently found free floating in the cytoplasm. We conclude that unmyelinated axons are fundamentally varicose structures created by the presence of organelles, and that an axon's calibre is dynamic in both space and time. These irregular axonal geometries raise serious doubts about standard two dimensional morphometric analysis and suggest that electrical properties may be more heterogeneous than expected from single section data. These results also suggest that the total number of microtubules contained in an axon, rather than its single section diameter, may prove to be a more accurate predictor of properties such as conduction velocity. Finally, these results offer an explanation for a number of pathological changes that have been described in unmyelinated axons. PMID:2292722
NASA Astrophysics Data System (ADS)
Watt, J. T.; Ponce, D. A.; Johnson, S. Y.
2014-12-01
Integrated geological and geophysical studies of three fault intersections in central California reveal an abundance of serpentinite and related ultramafic ophiolitic rocks in the subsurface; this is a significant discovery that not only constrains the geometry of structures, but also helps to explain fault development and fault behavior. We present and compare fault framework models from the Hayward-Calaveras, Calaveras-San Andreas, and Hosgri-Shoreline fault intersections and discuss the importance of ultramafic ophiolitic rocks and the implications for seismic hazard. Fault framework models were created based on integrated analysis of geology, potential-field, seismic-reflection, multibeam bathymetry, and seismicity data. In each case, ultramafic ophiolitic rocks have very limited surface expression, but are prevalent in the subsurface from ~1 to ~8 km depth. High-resolution aero- and marine-magnetic data are particularly important for determining the geometry of these strongly magnetic rock bodies and the faults along which they occur. At each fault intersection, the ultramafic ophiolitic rocks characterize a pre-existing structure that either promotes or represents a barrier to faulting, and significantly influences local deformation patterns. Although the faults at each intersection do not connect directly at the surface, combined geologic and geophysical data suggest that these faults connect in the subsurface either along a single through-going fault plane or along multiple splay faults. The continuity of these fault intersections at seismogenic depths suggests the possibility of a combined rupture that should be evaluated in future seismic hazard studies.
Sebastian Schunert; Yousry Y. Azmy; Damien Fournier
2011-05-01
We present a comprehensive error estimation of four spatial discretization schemes of the two-dimensional Discrete Ordinates (SN) equations on Cartesian grids utilizing a Method of Manufactured Solution (MMS) benchmark suite based on variants of Larsen’s benchmark featuring different orders of smoothness of the underlying exact solution. The considered spatial discretization schemes include the arbitrarily high order transport methods of the nodal (AHOTN) and characteristic (AHOTC) types, the discontinuous Galerkin Finite Element method (DGFEM) and the recently proposed higher order diamond difference method (HODD) of spatial expansion orders 0 through 3. While AHOTN and AHOTC rely on approximate analytical solutions of the transport equation within a mesh cell, DGFEM and HODD utilize a polynomial expansion to mimick the angular flux profile across each mesh cell. Intuitively, due to the higher degree of analyticity, we expect AHOTN and AHOTC to feature superior accuracy compared with DGFEM and HODD, but at the price of potentially longer grind times and numerical instabilities. The latter disadvantages can result from the presence of exponential terms evaluated at the cell optical thickness that arise from the semianalytical solution process. This work quantifies the order of accuracy and the magnitude of the error of all four discretization methods for different optical thicknesses, scattering ratios and degrees of smoothness of the underlying exact solutions in order to verify or contradict the aforementioned intuitive expectation.
FARGO3D: Hydrodynamics/magnetohydrodynamics code
NASA Astrophysics Data System (ADS)
Benítez Llambay, Pablo; Masset, Frédéric
2015-09-01
A successor of FARGO (ascl:1102.017), FARGO3D is a versatile HD/MHD code that runs on clusters of CPUs or GPUs, with special emphasis on protoplanetary disks. FARGO3D offers Cartesian, cylindrical or spherical geometry; 1-, 2- or 3-dimensional calculations; and orbital advection (aka FARGO) for HD and MHD calculations. As in FARGO, a simple Runge-Kutta N-body solver may be used to describe the orbital evolution of embedded point-like objects. There is no need to know CUDA; users can develop new functions in C and have them translated to CUDA automatically to run on GPUs.
Full-field drift Hamiltonian particle orbits in 3D geometry
NASA Astrophysics Data System (ADS)
Cooper, W. A.; Graves, J. P.; Brunner, S.; Isaev, M. Yu
2011-02-01
A Hamiltonian/Lagrangian theory to describe guiding centre orbit drift motion which is canonical in the Boozer coordinate frame has been extended to include full electromagnetic perturbed fields in anisotropic pressure 3D equilibria with nested magnetic flux surfaces. A redefinition of the guiding centre velocity to eliminate the motion due to finite equilibrium radial magnetic fields and the choice of a gauge condition that sets the radial component of the electromagnetic vector potential to zero are invoked to guarantee that the Boozer angular coordinates retain the canonical structure. The canonical momenta are identified and the guiding centre particle radial drift motion and parallel gyroradius evolution are derived. The particle coordinate position is linearly modified by wave-particle interactions. All the nonlinear wave-wave interactions appear explicitly only in the evolution of the parallel gyroradius. The radial variation of the electrostatic potential is related to the binormal component of the displacement vector for MHD-type perturbations. The electromagnetic vector potential projections can then be determined from the electrostatic potential and the radial component of the MHD displacement vector.
Non-invasive 3D geometry extraction of a Sea lion foreflipper
NASA Astrophysics Data System (ADS)
Friedman, Chen; Watson, Martha; Zhang, Pamela; Leftwich, Megan
2015-11-01
We are interested in underwater propulsion that leaves little traceable wake structure while producing high levels of thrust. A potential biological model is the California sea lion, a highly maneuverable aquatic mammal that produces thrust primarily with its foreflippers without a characteristic flapping frequency. The foreflippers are used for thrust, stability, and control during swimming motions. Recently, the flipper's kinematics during the thrust phase was extracted using 2D video tracking. This work extends the tracking ability to 3D using a non-invasive Direct Linear Transformation technique employed on non-research sea lions. marker-less flipper tracking is carried out manually for complete dorsal-ventral flipper motions. Two cameras are used (3840 × 2160 pixels resolution), calibrated in space using a calibration target inserted into the sea lion habitat, and synchronized in time using a simple light flash. The repeatability and objectivity of the tracked data is assessed by having two people tracking the same clap and comparing the results. The number of points required to track a flipper with sufficient detail is also discussed. Changes in the flipper pitch angle during the clap, an important feature for fluid dynamics modeling, will also be presented.
Flow properties along field lines in a 3-D tilted-dipole geometry
NASA Technical Reports Server (NTRS)
Pizzo, V. J.
1995-01-01
A 3-D MHD simulation of a global, tilted-dipole solar wind flow pattern is analyzed to determine flow properties along individual magnetic field lines. In the model, flow conditions near the Sun are chosen to provide a reasonable match to the interplanetary configuration prevailing during the recent south polar passage by Ulysses, i.e., a streamer belt inclined approximately 30 deg to the solar equator and speeds ranging from 325-800 km/s. Field lines all across the stream pattern are traced from 1 to 10 AU by following the motion of marker particles embedded in the flow. It is found that those field lines threading the core of the interaction region are subject to significant latitudinal and relative longitudinal displacement over this range of heliocentric distance. Thus, observations taken at a fixed latitude in the inner solar system sample, over the course of a solar rotation, field lines which connect to a range of latitudes in the outer heliosphere. Maps of the field line displacements are presented to help visualize these connections. In addition, it is found that depending upon the location relative to the CIR structure, the radial evolution of fluid elements frozen to different field lines can deviate considerably from that of the canonical solar wind. That is, for selected subsets of field lines, large speed changes (not just at shocks) can be experienced; the density variation can be far from 1/r(exp 2), and the magnetic field intensity need not decay monotonically with distance.
High order spatial expansion for the method of characteristics applied to 3-D geometries
Naymeh, L.; Masiello, E.; Sanchez, R.
2013-07-01
The method of characteristics is an efficient and flexible technique to solve the neutron transport equation and has been extensively used in two-dimensional calculations because it permits to deal with complex geometries. However, because of a very fast increase in storage requirements and number of floating operations, its direct application to three-dimensional routine transport calculations it is not still possible. In this work we introduce and analyze several modifications aimed to reduce memory requirements and to diminish the computing burden. We explore high-order spatial approximation, the use of intermediary trajectory-dependent flux expansions and the possibility of dynamic trajectory reconstruction from local tracking for typed subdomains. (authors)
Characterization of 3D filament dynamics in a MAST SOL flux tube geometry
NASA Astrophysics Data System (ADS)
Walkden, N. R.; Dudson, B. D.; Fishpool, G.
2013-10-01
Non-linear simulations of filament propagation in a realistic MAST SOL flux tube geometry using the BOUT++ fluid modelling framework show an isolation of the dynamics of the filament in the divertor region from the midplane region due to three features of the magnetic geometry; the variation of magnetic curvature along the field line, the expansion of the flux tube and strong magnetic shear. Of the three effects, the latter two lead to a midplane ballooning feature of the filament, whilst the former leads to a ballooning around the X-points. In simulations containing all three effects the filament is observed to balloon at the midplane, suggesting that the role of curvature variation is sub-dominant to the flux expansion and magnetic shear. The magnitudes of these effects are all strongest near the X-point which leads to the formation of parallel density gradients. The filaments simulated, which represent filaments in MAST, are identified as resistive ballooning, meaning that their motion is inertially limited, not sheath limited. Parallel density gradients can drive the filament towards a Boltzmann response when the collisionalityof the plasma is low. The results here show that the formation of parallel density gradients is a natural and inevitable consequence of a realistic magnetic geometry and therefore the transition to the Boltzmann response is a consequence of the use of realistic magnetic geometry and does not require initializing specifically varying background profiles as in slab simulations. The filaments studied here are stable to the linear resistive drift-wave instability but are subject to the non-linear effects associated with the Boltzmann response, particularly Boltzmann spinning. The Boltzmann response causes the filament to spin on an axis. In later stages of its evolution a non-linear turbulent state develops where the vorticity evolves into a turbulent eddy field on the same length scale as the parallel current. The transition from interchange
TART97 a coupled neutron-photon 3-D, combinatorial geometry Monte Carlo transport code
Cullen, D.E.
1997-11-22
TART97 is a coupled neutron-photon, 3 Dimensional, combinatorial geometry, time dependent Monte Carlo transport code. This code can on any modern computer. It is a complete system to assist you with input preparation, running Monte Carlo calculations, and analysis of output results. TART97 is also incredibly FAST; if you have used similar codes, you will be amazed at how fast this code is compared to other similar codes. Use of the entire system can save you a great deal of time and energy. TART97 is distributed on CD. This CD contains on- line documentation for all codes included in the system, the codes configured to run on a variety of computers, and many example problems that you can use to familiarize yourself with the system. TART97 completely supersedes all older versions of TART, and it is strongly recommended that users only use the most recent version of TART97 and its data riles.
THE EFFECTS OF APONEUROSIS GEOMETRY ON STRAIN INJURY SUSCEPTIBILITY EXPLORED WITH A 3D MUSCLE MODEL
Rehorn, Michael R.; Blemker, Silvia S.
2010-01-01
In the musculoskeletal system, some muscles are injured more frequently than others. For example, the biceps femoris longhead (BFLH) is the most commonly injured hamstring muscle. It is thought that acute injuries result from large strains within the muscle tissue, but the mechanism behind this type of strain injury is still poorly understood. The purpose of this study was to build computational models to analyze the stretch distributions within the BFLH muscle and to explore the effects of aponeurosis geometry on the magnitude and location of peak stretches within the model. We created a three-dimensional finite element (FE) model of the BFLH based on magnetic resonance (MR) images. We also created a series of simplified models with a similar geometry to the MR-based model. We analyzed the stretches predicted by the MR-based model during lengthening contractions to determine the region of peak local fiber stretch. The peak along-fiber stretch was 1.64 and was located adjacent to the proximal myotendinous junction (MTJ). In contrast, the average along-fiber stretch across all the muscle tissue was 0.95. By analyzing the simple models, we found that varying the dimensions of the aponeuroses (width, length, and thickness) had a substantial impact on the location and magnitude of peak stretches within the muscle. Specifically, the difference in widths between the proximal and distal aponeurosis in the BFLH contributed most to the location and magnitude of peak stretch, as decreasing the proximal aponeurosis width by 80% increased peak average stretches along the proximal MTJ by greater than 60% while slightly decreasing stretches along the distal MTJ. These results suggest that the aponeurosis morphology of the BFLH plays a significant role in determining stretch distributions throughout the muscle. Furthermore, this study introduces the new hypothesis that aponeurosis widths may be important in determining muscle injury susceptibility. PMID:20541207
NASA Technical Reports Server (NTRS)
Temmer, M.; Veronig, A. M.; Gopalswamy, N.; Yashiro, S.
2011-01-01
We study the kinematical characteristics and 3D geometry of a large-scale coronal wave that occurred in association with the 26 April 2008 flare-CME event. The wave was observed with the EUVI instruments aboard both STEREO spacecraft (STEREO-A and STEREO-B) with a mean speed of approx 240 km/s. The wave is more pronounced in the eastern propagation direction, and is thus, better observable in STEREO-B images. From STEREO-B observations we derive two separate initiation centers for the wave, and their locations fit with the coronal dimming regions. Assuming a simple geometry of the wave we reconstruct its 3D nature from combined STEREO-A and STEREO-B observations. We find that the wave structure is asymmetric with an inclination toward East. The associated CME has a deprojected speed of approx 750 +/- 50 km/s, and it shows a non-radial outward motion toward the East with respect to the underlying source region location. Applying the forward fitting model developed by Thernisien, Howard, and Vourlidas we derive the CME flux rope position on the solar surface to be close to the dimming regions. We conclude that the expanding flanks of the CME most likely drive and shape the coronal wave.
NASA Astrophysics Data System (ADS)
Keller, Lukas M.; Holzer, Lorenz; Wepf, Roger; Gasser, Philippe; Münch, Beat; Marschall, Paul
The evaluation and optimization of radioactive disposal systems requires a comprehensive understanding of mass transport processes. Among others, mass transport in porous geomaterials depends crucially on the topology and geometry of the pore space. Thus, understanding the mechanism of mass transport processes ultimately requires a 3D characterization of the pore structure. Here, we demonstrate the potential of focused ion beam nanotomography (FIB-nT) in characterizing the 3D geometry of pore space in clay rocks, i.e. Opalinus clay. In order to preserve the microstructure and to reduce sample preparation artefacts we used high pressure freezing and subsequent freeze drying to prepare the samples. Resolution limitations placed the lower limit in pore radii that can be analyzed by FIB-nT to about 10-15 nm. Image analysis and the calculation of pore size distribution revealed that pores with radii larger than 15 nm are related to a porosity of about 3 vol.%. To validate the method, we compared the pores size distribution obtained by FIB-nT with the one obtained by N 2 adsorption analysis. The latter yielded a porosity of about 13 vol.%. This means that FIB-nT can describe around 20-30% of the total pore space. For pore radii larger than 15 nm the pore size distribution obtained by FIB-nT and N 2 adsorption analysis were in good agreement. This suggests that FIB-nT can provide representative data on the spatial distribution of pores for pore sizes in the range of about 10-100 nm. Based on the spatial analysis of 3D data we extracted information on the spatial distribution of pore space geometrical properties.
A Coupled Neutron-Photon 3-D Combinatorial Geometry Monte Carlo Transport Code
Energy Science and Technology Software Center (ESTSC)
1998-06-12
TART97 is a coupled neutron-photon, 3 dimensional, combinatorial geometry, time dependent Monte Carlo transport code. This code can run on any modern computer. It is a complete system to assist you with input preparation, running Monte Carlo calculations, and analysis of output results. TART97 is also incredibly fast: if you have used similar codes, you will be amazed at how fast this code is compared to other similar codes. Use of the entire system canmore » save you a great deal of time and energy. TART 97 is distributed on CD. This CD contains on-line documentation for all codes included in the system, the codes configured to run on a variety of computers, and many example problems that you can use to familiarize yourself with the system. TART97 completely supersedes all older versions of TART, and it is strongly recommended that users only use the most recent version of TART97 and ist data files.« less
A 3D Computational fluid dynamics model validation for candidate molybdenum-99 target geometry
NASA Astrophysics Data System (ADS)
Zheng, Lin; Dale, Greg; Vorobieff, Peter
2014-11-01
Molybdenum-99 (99Mo) is the parent product of technetium-99m (99mTc), a radioisotope used in approximately 50,000 medical diagnostic tests per day in the U.S. The primary uses of this product include detection of heart disease, cancer, study of organ structure and function, and other applications. The US Department of Energy seeks new methods for generating 99Mo without the use of highly enriched uranium, to eliminate proliferation issues and provide a domestic supply of 99mTc for medical imaging. For this project, electron accelerating technology is used by sending an electron beam through a series of 100Mo targets. During this process a large amount of heat is created, which directly affects the operating temperature dictated by the tensile stress limit of the wall material. To maintain the required temperature range, helium gas is used as a cooling agent that flows through narrow channels between the target disks. In our numerical study, we investigate the cooling performance on a series of new geometry designs of the cooling channel. This research is supported by Los Alamos National Laboratory.
Validation and Analysis of Forward Osmosis CFD Model in Complex 3D Geometries
Gruber, Mathias F.; Johnson, Carl J.; Tang, Chuyang; Jensen, Mogens H.; Yde, Lars; Hélix-Nielsen, Claus
2012-01-01
In forward osmosis (FO), an osmotic pressure gradient generated across a semi-permeable membrane is used to generate water transport from a dilute feed solution into a concentrated draw solution. This principle has shown great promise in the areas of water purification, wastewater treatment, seawater desalination and power generation. To ease optimization and increase understanding of membrane systems, it is desirable to have a comprehensive model that allows for easy investigation of all the major parameters in the separation process. Here we present experimental validation of a computational fluid dynamics (CFD) model developed to simulate FO experiments with asymmetric membranes. Simulations are compared with experimental results obtained from using two distinctly different complex three-dimensional membrane chambers. It is found that the CFD model accurately describes the solute separation process and water permeation through membranes under various flow conditions. It is furthermore demonstrated how the CFD model can be used to optimize membrane geometry in such as way as to promote the mass transfer. PMID:24958428
Validation and Analysis of Forward Osmosis CFD Model in Complex 3D Geometries.
Gruber, Mathias F; Johnson, Carl J; Tang, Chuyang; Jensen, Mogens H; Yde, Lars; Hélix-Nielsen, Claus
2012-01-01
In forward osmosis (FO), an osmotic pressure gradient generated across a semi-permeable membrane is used to generate water transport from a dilute feed solution into a concentrated draw solution. This principle has shown great promise in the areas of water purification, wastewater treatment, seawater desalination and power generation. To ease optimization and increase understanding of membrane systems, it is desirable to have a comprehensive model that allows for easy investigation of all the major parameters in the separation process. Here we present experimental validation of a computational fluid dynamics (CFD) model developed to simulate FO experiments with asymmetric membranes. Simulations are compared with experimental results obtained from using two distinctly different complex three-dimensional membrane chambers. It is found that the CFD model accurately describes the solute separation process and water permeation through membranes under various flow conditions. It is furthermore demonstrated how the CFD model can be used to optimize membrane geometry in such as way as to promote the mass transfer. PMID:24958428
3D geometry and evolutionary sequence of fold-thrust systems in NW Taiwan
NASA Astrophysics Data System (ADS)
Chang, Hao-Yun; Yang, Kenn-Ming; Hsieh, Ching-Yun; Yang, Tzu-Ruei; Chuang, Hui-Ju; Chen, Yi-Ju
2015-04-01
During the arc-continental collision from the Pliocene to the Pleistocene, two sets of fold-and-thrust system developed in NW Taiwan, a series of NNE-SSW striking low-angle thrust faults and their related folds (set A) and the other series of NEE-SWW striking high-angle thrust faults and their related folds (set B). The latter one cuts the former one and extends forelandward. The geometry of intersection and development sequence of both sets of structures are still in debate. In this study, we utilized a grid of seismic profiles to constrain our interpretation on the subsurface structural geometry of the two structural sets, which then was tested by structural restoration. We also made some simulations on the formation of fault-related folds by trishear modeling. The influence of normal fault reactivation on and the transitional relationships among the structures were investigated to establish an evolutionary sequence for the fold-and thrust systems of NW Taiwan. The strike of set A is NNE-SSW in the northern part of the study area but becomes N-S to the south. The location of the strike change is cut by a NEE-SWW high-angle fault of set B. According to the seismic interpretation, shallower anticline is asymmetric whereas deeper anticline is symmetric. The low-angle thrust of set A extends to the south and transfers into high-angle where it is cut by the high-angle fault of set B. The trishear model suggests that the shallower anticline resulted from low-angle fault thrusting in early period, whereas the deeper one was caused by basal detachment faulting in the late stage. Seismic interpretation also reveals an asymmetric and gentle fold cut by a high-angle thrust fault of set B. The result of trishear modeling indicates that the anticline was formed by slip along a high angle thrust, which is a low-angle fault in the deep but turns into high angle along a pre-existing normal fault up to the surface. In summary, the development of the shallower anticline of set A
3D Geometry and Kinematics of the Taiwan Arc-continent Collision
NASA Astrophysics Data System (ADS)
Carena, S.; Suppe, J.; Wu, Y. M.
2015-12-01
In Taiwan two subduction zones (Manila trench and Ryukyu trench) come together in a quasi-orthogonal, kinematically stable configuration. Subduction is ongoing in both trenches, even though the tectonic setting in the Manila trench is that of an arc-continent collision. The upper crust of Eurasia (EU) is decoupled from the rest of the lithosphere by a detachment horizon, which is the main subduction interface between EU and Philippine Sea plate (PSP). The interface is visible in both seismicity and crustal tomography at shallow depths, and it can be followed into the mantle to 450-500 km depth with global tomography. Shortening across the plate boundary is accomplished by a combination of subduction of EU lithosphere, folding and thrusting in the Eurasian upper crust, and a secondary subduction zone within the PSP. We hypothesize that: (1) once arc-continent collision occurs, subduction of Eurasian continental lower crust and upper mantle can continue by lithospheric delamination and by continuity with the much larger Eurasian slab to the south; (2) the upper crust of EU deforms by faulting and folding; (3) the present convergence rate of about 90 mm/yr is limited at most to the last 2 Ma, whereas the long-term rate is about 30 mm/yr and in Taiwan the difference is being taken up by secondary subduction within the PSP margin; (4) a margin-parallel STEP (Subduction-Transform-Edge-Propagator) fault forms the northern limit of Eurasian subduction, which allows the whole system to propagate self-similarly southwestward. No slab breakoff is required for the kinematics of the margin, and none is observed in geophysical or geological data either. This kinematics is consistent with geologic observations: from timing of opening of the southern Okinawa trough, to geometry of geologic boundaries within the Taiwan mountain belt, to geographic distribution, geochemical character, and timing of Quaternary volcanism in the northern Taiwan volcanic zone. We constrained the long
NASA Astrophysics Data System (ADS)
Coetzee, André; Kisters, Alexander
2016-05-01
Dolerites in the Karoo Basin of South Africa commonly represent kilometre-scale, interconnected saucer-shaped structures that consist of inner sills, bounded by inclined sheets connected to stratigraphically higher outer sills. Based on information from over 3000 boreholes and mining operations extending over an area of ca. 500 km2 and covering a > 3 km vertical section from Karoo strata into underlying basement rocks, this paper presents the results of a 3D modelling exercise that describes the geometry and spatial relationships of a regional-scale saucer complex, locally referred to as the number 8 sill, from the Secunda (coal mine) Complex in the northern parts of the Karoo Basin. The composite number 8 sill complex consists of three main dolerite saucers (dolerites A to C). These dolerite saucers are hosted by the Karoo Supergroup and the connectivity and geometry of the saucers support a lateral, sill-feeding-sill relationship between dolerite saucers A, B and C. The saucers are underlain and fed by a shallowly-dipping sheet (dolerite D) in the basement rocks below the Karoo sequence. The 3D geometric strata model agrees well with experimental results of saucer formation from underlying feeders in sedimentary basins, but demonstrates a more intricate relationship where a single feeder can give rise to several split level saucers in one regionally extensive saucer complex. More localised dome- or ridge-shape protrusions are common in the flat lying sill parts of the regional-scale saucers. We suggest a mode of emplacement for these kilometre-scale dome- and ridge structures having formed as a result of lobate magma flow processes. Magma lobes, propagating in different directions ahead of the main magma sheet, undergo successive episodes of lobe arrest and inflation. The inflation of lobes initiates failure of the overlying strata and the formation of curved faults. Magma exploiting these faults transgresses the stratigraphy and coalesces to form a ring
NASA Astrophysics Data System (ADS)
Van Gent, Heijn; Urai, Janos L.; de Keijzer, Martin
2011-03-01
We present a first look at the large-scale, complexly folded and faulted internal structure of Zechstein salt bodies in NW Europe using 3D reflection seismic reflection data from two surveys on the Groningen High and the Cleaver Bank High. We focus on a relatively brittle, folded and boudinaged, claystone-carbonate-anhydrite layer (the Z3 stringer) enclosed in ductile salt. A first classification of the structures is presented and compared with observations from salt mines and analogue and numerical models. Z3 stringers not only are reservoirs for hydrocarbons but can also present a serious drilling problem in some areas. Results of this study could provide the basis for better prediction of zones of drilling problems. More generally, the techniques presented here can be used to predict the internal structure of salt bodies, to estimate the geometry of economic deposits of all kinds and locate zones suitable for storage caverns. Structures observed include an extensive network of zones with increased thickness of the stringer. These we infer to have formed by early diagenesis, karstification, gravitational sliding and associated local sedimentation. Later, this template was deformed into large-scale folds and boudins during salt tectonics. Salt flow was rarely plane strain, producing complex fold and boudin geometries. Deformation was further complicated by the stronger zones of increased thickness, which led to strongly non-cylindrical structures. We present some indications that the thicker zones also influence the locations of later suprasalt structures, suggesting a feedback between the early internal evolution of this salt giant and later salt tectonics. This study opens the possibility to study the internal structure of the Zechstein and other salt giants in 3D using this technique, exposing a previously poorly known structure which is comparable in size and complexity to the internal parts of some orogens.
Freels, James D; Jain, Prashant K
2011-01-01
A research and development project is ongoing to convert the currently operating High Flux Isotope Reactor (HFIR) of Oak Ridge National Laboratory (ORNL) from highly-enriched Uranium (HEU U3O8) fuel to low-enriched Uranium (LEU U-10Mo) fuel. Because LEU HFIR-specific testing and experiments will be limited, COMSOL is chosen to provide the needed multiphysics simulation capability to validate against the HEU design data and calculations, and predict the performance of the LEU fuel for design and safety analyses. The focus of this paper is on the unique issues associated with COMSOL modeling of the 3D geometry, meshing, and solution of the HFIR fuel plate and assembled fuel elements. Two parallel paths of 3D model development are underway. The first path follows the traditional route through examination of all flow and heat transfer details using the Low-Reynolds number k-e turbulence model provided by COMSOL v4.2. The second path simplifies the fluid channel modeling by taking advantage of the wealth of knowledge provided by decades of design and safety analyses, data from experiments and tests, and HFIR operation. By simplifying the fluid channel, a significant level of complexity and computer resource requirements are reduced, while also expanding the level and type of analysis that can be performed with COMSOL. Comparison and confirmation of validity of the first (detailed) and second (simplified) 3D modeling paths with each other, and with available data, will enable an expanded level of analysis. The detailed model will be used to analyze hot-spots and other micro fuel behavior events. The simplified model will be used to analyze events such as routine heat-up and expansion of the entire fuel element, and flow blockage. Preliminary, coarse-mesh model results of the detailed individual fuel plate are presented. Examples of the solution for an entire fuel element consisting of multiple individual fuel plates produced by the simplified model are also presented.
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.
Ramanathan, Muruganathan; Lokitz, Bradley S.; Messman, Jamie M.; Stafford, Christopher M.; Kilbey II, S. Michael
2013-01-01
We report a simple, one step process for developing wrinkling patterns in azlactone-based polymer thin films and brushes in 2D and 3D surfaces. The polymer used in this work wrinkles spontaneously upon deposition and solidification on a substrate without applying any external strain to the substrate, with the mode of deposition defining the direction of the wrinkles. Wrinkle formation is shown to occur on a variety of substrates over large areas. We also find that a very thin brush-like layer of an azlactone-containing block copolymer also exhibits wrinkled topology. Given the spontaneity and versatility of wrinkle formation, we further demonstrate two proofs-of-concept, i) that these periodic wrinkled structures are not limited to planar surfaces, but are also developed in complex geometries including tubes, cones and other 3D structures; and ii) that this one-step wrinkling process can be used to guide the deposition of metal nanoparticles and quantum dots, creating a periodic, nanopatterned film.
NASA Astrophysics Data System (ADS)
Jensen, E. A.; Mulligan, T.; Jackson, B. V.; Tokumaru, M.
2006-12-01
Multiple reconstructions of the October 28-29, 2003 CME/ICME using white-light observations, ground-based cosmic-ray and in situ magnetic field flux rope modeling show two possible flux-rope configurations that pass Earth on opposite sides of the central symmetry axis of the disturbance. An analysis of flux rope model geometries initiated over a wide range in parameter space to test the uniqueness of the single spacecraft inversion reveals the fit is degenerate over a range of impact parameters such that two solutions are obtained. In one case (fit A) the disturbance passes Earth to the west of the rope center with the rope axis at a low inclination of 20 deg to the ecliptic, similar to the ground-based flux rope analysis by Kuwabara et al.~(2004). In the second case (fit B) the disturbance passes Earth to the east of the flux rope axis, with the rope axis more highly inclined at 42 deg from the ecliptic, consistent with the SMEI white-light analysis of Jackson et al.~(2006). The current densities in both solutions indicate a nearly force-free structure. Multipoint studies of ICMEs show the radius of curvature in the plane of the rope is between that of a dipole field line connected to the Sun and that of a circular field line connected to the Sun. Assuming a dipole field geometry for the large- scale axial field curvature of the rope results in a 3-D reconstruction for case B that is consistent with the loop structure and observed speed in the white-light LASCO images and SMEI density reconstruction, but not for case A. Multipoint measurements of large-scale solar wind transients is one of the key objectives of the Stereo mission, allowing more accurate 3-D reconstructions of in situ data for comparison with white-light observations. Until they become available, the large-scale axial field orientation and loop geometry of these rope reconstructions provides another tool to constrain magnetic flux rope fits of ICMEs using single spacecraft measurements.
NASA Astrophysics Data System (ADS)
Maginot, Peter G.; Morel, Jim E.; Ragusa, Jean C.
2012-08-01
We present a new nonlinear spatial finite-element method for the linearized Boltzmann transport equation with Sn angular discretization in 1-D and 2-D Cartesian geometries. This method has two central characteristics. First, it is equivalent to the linear-discontinuous (LD) Galerkin method whenever that method yields a strictly non-negative solution. Second, it always satisfies both the zeroth and first spatial moment equations. Because it yields the LD solution when that solution is non-negative, one might interpret our method as a classical fix-up to the LD scheme. However, fix-up schemes for the LD equations derived in the past have given up solution of the first moment equations when the LD solution is negative in order to satisfy positivity in a simple manner. We present computational results comparing our method in 1-D to the strictly non-negative linear exponential-discontinuous method and to the LD method. We present computational results in 2-D comparing our method to a recently developed LD fix-up scheme and to the LD scheme. It is demonstrated that our method is a valuable alternative to existing methods.
NASA Astrophysics Data System (ADS)
Lücke, O. H.; Arroyo, I. G.
2015-10-01
The eastern part of the oceanic Cocos Plate presents a heterogeneous crustal structure due to diverse origins and ages as well as plate-hot spot interactions which originated the Cocos Ridge, a structure that converges with the Caribbean Plate in southeastern Costa Rica. The complex structure of the oceanic plate directly influences the dynamics and geometry of the subduction zone along the Middle American Trench. In this paper an integrated interpretation of the slab geometry in Costa Rica is presented based on 3-D density modeling of combined satellite and surface gravity data, constrained by available geophysical and geological data and seismological information obtained from local networks. The results show the continuation of steep subduction geometry from the Nicaraguan margin into northwestern Costa Rica, followed by a moderate dipping slab under the Central Cordillera toward the end of the Central American Volcanic Arc. Contrary to commonly assumed, to the southeast end of the volcanic arc, our preferred model shows a steep, coherent slab that extends up to the landward projection of the Panama Fracture Zone. Overall, a gradual change in the depth of the intraplate seismicity is observed, reaching 220 km in the northwestern part, and becoming progressively shallower toward the southeast, where it reaches a maximum depth of 75 km. The changes in the terminal depth of the observed seismicity correlate with the increased density in the modeled slab. The absence of intermediate depth (> 75 km) intraplate seismicity in the southeastern section and the higher densities for the subducted slab in this area, support a model in which dehydration reactions in the subducted slab cease at a shallower depth, originating an anhydrous and thus aseismic slab.
Brosed, Francisco Javier; Aguilar, Juan José; Guillomía, David; Santolaria, Jorge
2011-01-01
This article discusses different non contact 3D measuring strategies and presents a model for measuring complex geometry parts, manipulated through a robot arm, using a novel vision system consisting of a laser triangulation sensor and a motorized linear stage. First, the geometric model incorporating an automatic simple module for long term stability improvement will be outlined in the article. The new method used in the automatic module allows the sensor set up, including the motorized linear stage, for the scanning avoiding external measurement devices. In the measurement model the robot is just a positioning of parts with high repeatability. Its position and orientation data are not used for the measurement and therefore it is not directly "coupled" as an active component in the model. The function of the robot is to present the various surfaces of the workpiece along the measurement range of the vision system, which is responsible for the measurement. Thus, the whole system is not affected by the robot own errors following a trajectory, except those due to the lack of static repeatability. For the indirect link between the vision system and the robot, the original model developed needs only one first piece measuring as a "zero" or master piece, known by its accurate measurement using, for example, a Coordinate Measurement Machine. The strategy proposed presents a different approach to traditional laser triangulation systems on board the robot in order to improve the measurement accuracy, and several important cues for self-recalibration are explored using only a master piece. Experimental results are also presented to demonstrate the technique and the final 3D measurement accuracy. PMID:22346569
Brosed, Francisco Javier; Aguilar, Juan José; Guillomía, David; Santolaria, Jorge
2011-01-01
This article discusses different non contact 3D measuring strategies and presents a model for measuring complex geometry parts, manipulated through a robot arm, using a novel vision system consisting of a laser triangulation sensor and a motorized linear stage. First, the geometric model incorporating an automatic simple module for long term stability improvement will be outlined in the article. The new method used in the automatic module allows the sensor set up, including the motorized linear stage, for the scanning avoiding external measurement devices. In the measurement model the robot is just a positioning of parts with high repeatability. Its position and orientation data are not used for the measurement and therefore it is not directly “coupled” as an active component in the model. The function of the robot is to present the various surfaces of the workpiece along the measurement range of the vision system, which is responsible for the measurement. Thus, the whole system is not affected by the robot own errors following a trajectory, except those due to the lack of static repeatability. For the indirect link between the vision system and the robot, the original model developed needs only one first piece measuring as a “zero” or master piece, known by its accurate measurement using, for example, a Coordinate Measurement Machine. The strategy proposed presents a different approach to traditional laser triangulation systems on board the robot in order to improve the measurement accuracy, and several important cues for self-recalibration are explored using only a master piece. Experimental results are also presented to demonstrate the technique and the final 3D measurement accuracy. PMID:22346569
NASA Astrophysics Data System (ADS)
Hartl, Darren J.; Lagoudas, Dimitris C.
2007-04-01
This work describes the thermomechanical characterization and FEA modeling of commercial jet engine chevrons incorporating active Shape Memory Alloy (SMA) beam components. The reduction of community noise at airports generated during aircraft take-off has become a major research goal. Serrated aerodynamic devices along the trailing edge of a jet engine primary and secondary exhaust nozzle, known as chevrons, have been shown to greatly reduce jet noise by encouraging advantageous mixing of the streams. To achieve the noise reduction, the secondary exhaust nozzle chevrons are typically immersed into the fan flow which results in drag, or thrust losses during cruise. SMA materials have been applied to this problem of jet engine noise. Active chevrons, utilizing SMA components, have been developed and tested to create maximum deflection during takeoff and landing while minimizing deflection into the flow during the remainder of flight, increasing efficiency. Boeing has flight tested one Variable Geometry Chevron (VGC) system which includes active SMA beams encased in a composite structure with a complex 3-D configuration. The SMA beams, when activated, induce the necessary bending forces on the chevron structure to deflect it into the fan flow and reduce noise. The SMA composition chosen for the fabrication of these beams is a Ni60Ti40 (wt%) alloy. In order to calibrate the material parameters of the constitutive SMA model, various thermomechanical experiments are performed on trained (stabilized) standard SMA tensile specimens. Primary among these tests are thermal cycles at various constant stress levels. Material properties for the shape memory alloy components are derived from this tensile experimentation. Using this data, a 3-D FEA implementation of a phenomenological SMA model is calibrated and used to analyze the response of the chevron. The primary focus of this work is the full 3-D modeling of the active chevron system behavior by considering the SMA beams as
Simulations of 3D LPI's relevant to IFE using the PIC code OSIRIS
NASA Astrophysics Data System (ADS)
Tsung, F. S.; Mori, W. B.; Winjum, B. J.
2014-10-01
We will study three dimensional effects of laser plasma instabilities, including backward raman scattering, the high frequency hybrid instability, and the two plasmon instability using OSIRIS in 3D Cartesian geometry and cylindrical 2D OSIRIS with azimuthal mode decompositions. With our new capabilities we hope to demonstrate that we are capable of studying single speckle physics relevant to IFE in an efficent manner.
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.
NASA Astrophysics Data System (ADS)
Ghisetti, F.; Sibson, R. H.; Hamling, I. J.
2014-12-01
The Alpine Fault is the principal component of the transform boundary between the Australian and Pacific plates across the South Island of New Zealand, linking the opposite dipping Hikurangi and Puysegur subduction zones. In the northern South Island, the transition from the subducted W-dipping Pacific slab of the Hikurangi margin to the intra-continental transform margin is defined by earthquake foci from 350 to 100 km deep. West of the Alpine Fault the Australian crust above the slab has been incorporated into the collisional plate boundary and uplifted in a compressional belt up to 100 km wide. Retro-deformation and back-stripping of 10 regional transects utilising surface geology, seismic reflection lines and exploration wells define the progressive deformation of the Australian crust since 35 Ma along the collisional margin. The reconstructed geometry of faulted basement blocks is tied to localisation and evolution of overlying sedimentary basins, coeval with displacement on the Alpine Fault. Amounts of shortening, uplift and subsidence and fault activity are heterogeneous in space and time across the margin, and are controlled by compressional reactivation of inherited high-angle, N-S Paleogene normal faults oblique to the margin. However, significant differences also occur along the strike of the collisional margin, with major contrasts in uplift and subsidence north and south of lat. 41°.7, i.e. the region overlying the southern termination of the Hikurangi slab. These differences are highlighted by present day hydrographic anomalies in the Buller region, and by the pattern of filtered topography at > 75 km wavelength. Our data show that the 3D geometry of the Australian plate cannot be entirely attributed to inherited crustal heterogeneity of a flexured "retro-foreland" domain in the footwall of the Alpine Fault, and suggest the need of deeper dynamic interaction between the Pacific and Australian lithosphere along the subduction-collision margin.
Modelling of 3D fields due to ferritic inserts and test blanket modules in toroidal geometry at ITER
NASA Astrophysics Data System (ADS)
Liu, Yueqiang; Äkäslompolo, Simppa; Cavinato, Mario; Koechl, Florian; Kurki-Suonio, Taina; Li, Li; Parail, Vassili; Saibene, Gabriella; Särkimäki, Konsta; Sipilä, Seppo; Varje, Jari
2016-06-01
Computations in toroidal geometry are systematically performed for the plasma response to 3D magnetic perturbations produced by ferritic inserts (FIs) and test blanket modules (TBMs) for four ITER plasma scenarios: the 15 MA baseline, the 12.5 MA hybrid, the 9 MA steady state, and the 7.5 MA half-field helium plasma. Due to the broad toroidal spectrum of the FI and TBM fields, the plasma response for all the n = 1–6 field components are computed and compared. The plasma response is found to be weak for the high-n (n > 4) components. The response is not globally sensitive to the toroidal plasma flow speed, as long as the latter is not reduced by an order of magnitude. This is essentially due to the strong screening effect occurring at a finite flow, as predicted for ITER plasmas. The ITER error field correction coils (EFCC) are used to compensate the n = 1 field errors produced by FIs and TBMs for the baseline scenario for the purpose of avoiding mode locking. It is found that the middle row of the EFCC, with a suitable toroidal phase for the coil current, can provide the best correction of these field errors, according to various optimisation criteria. On the other hand, even without correction, it is predicted that these n = 1 field errors will not cause substantial flow damping for the 15 MA baseline scenario.
NASA Astrophysics Data System (ADS)
Moustafa, Salli; Févotte, François; Lathuilière, Bruno; Plagne, Laurent
2014-06-01
The past few years have been marked by a noticeable increase in the interest in 3D whole-core heterogeneous deterministic neutron transport solvers for reference calculations. Due to the extremely large problem sizes tackled by such solvers, they need to use adapted numerical methods and need to be efficiently implemented to take advantage of the full computing power of modern systems. As for numerical methods, one possible approach consists in iterating over resolutions of 2D and 1D MOC problems by taking advantage of prismatic geometries. The MICADO solver, developed at EDF R&D, is a parallel implementation of such a method in distributed and shared memory systems. However it is currently unable to use SIMD vectorization to leverage the full computing power of modern CPUs. In this paper, we describe our first effort to support vectorization in MICADO, typically targeting Intel© SSE CPUs. Both the 2D and 1D algorithms are vectorized, allowing for high expected speedups for the whole spatial solver. We present benchmark computations, which show nearly optimal speedups for our vectorized implementation on the TAKEDA case.
NASA Astrophysics Data System (ADS)
Bergmann, Ryan
general 3D geometries on GPUs, but compared to production codes like Serpent and MCNP, WARP ha
NASA Astrophysics Data System (ADS)
Kanaujia, Jyotima; Kumar, Ashwani; Gupta, S. C.
2016-02-01
We investigate the upper crustal velocity structure beneath the Tehri region of the Garhwal Himalaya. The investigated region is situated within the 700-km-long central seismic gap of the Himalaya that has experienced three gap-filling earthquakes since 1991 including the recent 2015 Nepal earthquake (Mw 7.8). The local tomographic inversion is based on a dataset of 1365 events collected from January 2008 to December 2012 by a 12-station local network that covers an area of about 100 × 80 km around Tehri Dam. We perform a simultaneous inversion for P- and S-wave velocity anomalies. Tomograms are interpreted in the backdrop of the regional geological and tectonic framework of the region. The spatial distribution of relocated events from the 3- D velocity model has shed new light on the pattern of seismicity in the vicinity of the Main Central thrust (MCT), and has elucidated the structure of the underthrusting Indian plate. Our model exhibits a significant negative velocity anomaly up to ˜5 per cent beneath the central part of the Garhwal Inner Lesser Himalaya, and a P-wave low velocity anomaly near the Chamoli region. The seismicity zone around the Chamoli region may be attributed to the presence of fluid filled rocks. Furthermore, an area with˜3-4 per cent positive velocity anomaly is delineated to the northwest of the Uttarkashi thrust in the vicinity of the MCT. Significant findings of the study include: a flat-ramp-flat type sub-surface geometry of the underthrusting Indian plate below the Garhwal Himalaya, high velocity images representing the trend and configuration of Delhi-Haridwar-ridge below the Sub Himalaya and Lesser Himalaya, and a seismically active zone representing geometrical asperity on the basement thrust in the vicinity of the MCT.
NASA Astrophysics Data System (ADS)
Bergmann, Ryan
general 3D geometries on GPUs, but compared to production codes like Serpent and MCNP, WARP has limited capabilities. Despite WARP's lack of features, its novel algorithm implementations show that high performance can be achieved on a GPU despite the inherently divergent program flow and sparse data access patterns. WARP is not ready for everyday nuclear reactor calculations, but is a good platform for further development of GPU-accelerated Monte Carlo neutron transport. In it's current state, it may be a useful tool for multiplication factor searches, i.e. determining reactivity coefficients by perturbing material densities or temperatures, since these types of calculations typically do not require many flux tallies. (Abstract shortened by UMI.)
NASA Astrophysics Data System (ADS)
José Ramón, Mª; Pueyo, Emilio L.; Briz, José Luis; Caumon, Guillaume; Fernández, Óscar; Ciria, José Carlos; Pocovi, Andrés; Ros, Luis H.
2013-04-01
Three-dimensional reconstructions of the underground involve the integration of discrete and heterogeneous datasets and have significant socio- economic implications. The problem arises when there are limited data to build 3D models or when deformation processes are complex; these reasons inspired the development of restoration methods to validate subsurface reconstructions. The restoration is based on the application of simple geometric (or mechanic) laws that help reduce the uncertainty and increase geomodel accuracy. Apart from mechanical approaches, geometric methods are based on the initial assumption of global conservation of volume during deformation in addition to the paleo-horizontality of the stratigraphic horizons in the undeformed stage. The problem is that the bedding plane cannot be used as a three-dimensional reference system, because a single vector defines it and additional constraints are required. This is particularly important when dealing with complex structures, such as non-cylindrical structures and the superposition of non-coaxial geometries. In this context, paleomagnetism (known in both the deformed and undeformed stages) can contribute to building a more complete reference system and to reducing the uncertainty in restoration processes. The use of paleomagnetism in restoration tools was suggested in the early 1990's and only a few quantitative map-view applications have been developed since then. In this contribution, we introduce the two first surface restoration methods that use paleomagnetic vectors as a primary reference. The first one is a simple geometric approach based on the piecewise restoration of a triangulated surface into which paleomagnetic variables can be easily incorporated. It is valid for complexly folded structures. The surface is modelled by a mesh and the method starts from a pin-element. Triangles are laid flat and then fitted together to minimize distances between common vertices and paleomagnetic error. However
NASA Astrophysics Data System (ADS)
Munn, Jonathan; Parker, Beth
2013-04-01
In bedrock aquifers where matrix permeability is low, the nature and distribution of the fracture network has a strong impact on the transport and fate of contaminants. Accurate fracture characterization is therefore essential to fully understand the flow system and to predict contaminant migration. Powerful DFN models exist, yet the limitation is often on obtaining field data of sufficient quality to use as input parameters. One major contributing factor is the common practice of using only vertical coreholes to characterize bedrock aquifers. This can lead to datasets that are significantly biased toward fractures perpendicular to the corehole and are therefore not well suited for three-dimensional (3-D) fracture geometry characterization. This bias is particularly pronounced in flat-lying sedimentary strata where fracture networks are typically comprised of flat-lying bedding parallel fractures and vertical, or near vertical joints. An examination of such bias was conducted at a contaminated site in Guelph, Ontario, Canada, in a Silurian dolostone aquifer. Three inclined coreholes plunging 60 degrees with varying azimuths were drilled between 2010 and 2012 to supplement existing data from eleven vertical coreholes from previous investigations. Depth discrete datasets were collected in the coreholes including lithological and fracture logs from rock core, downhole geophysical surveys (e.g, acoustic televiewer, formation conductivity, temperature, natural gamma), and hydraulic testing including the first use of flexible liner profiling in inclined coreholes. These datasets were integrated to provide estimates of fracture frequency, orientation and aperture distributions and to estimate values of bulk effective fracture porosity. Orientation analysis revealed three dominant fracture sets on site that vary in intensity through mechanical layers. These sets consist of a horizontal, bedding-plane set with an average spacing of 0.3m, and two high-angle sets, NE-SW and
ERIC Educational Resources Information Center
Widder, Mirela; Gorsky, Paul
2013-01-01
In schools, learning spatial geometry is usually dependent upon a student's ability to visualize three dimensional geometric configurations from two dimensional drawings. Such a process, however, often creates visual obstacles which are unique to spatial geometry. Useful software programs which realistically depict three dimensional geometric…
Olubamiji, Adeola D; Izadifar, Zohreh; Si, Jennifer L; Cooper, David M L; Eames, B Frank; Chen, Daniel X B
2016-06-01
Three-dimensional (3D)-printed poly(ε)-caprolactone (PCL)-based scaffolds are increasingly being explored for cartilage tissue engineering (CTE) applications. However, ensuring that the mechanical properties of these PCL-based constructs are comparable to that of articular cartilage that they are meant to regenerate is an area that has been under-explored. This paper presents the effects of PCL's molecular weight (MW) and scaffold's pore geometric configurations; strand size (SZ), strand spacing (SS), and strand orientation (SO), on mechanical properties of 3D-printed PCL scaffolds. The results illustrate that MW has significant effect on compressive moduli and yield strength of 3D-printed PCL scaffolds. Specifically, PCL with MW of 45 K was a more feasible choice for fabrication of visco-elastic, flexible and load-bearing PCL scaffolds. Furthermore, pore geometric configurations; SZ, SS, and SO, all significantly affect on tensile moduli of scaffolds. However, only SZ and SS have statistically significant effects on compressive moduli and porosity of these scaffolds. That said, inverse linear relationship was observed between porosity and mechanical properties of 3D-printed PCL scaffolds in Pearson's correlation test. Altogether, this study illustrates that modulating MW of PCL and pore geometrical configurations of the scaffolds enabled design and fabrication of PCL scaffolds with mechanical and biomimetic properties that better mimic mechanical behaviour of human articular cartilage. Thus, the modulated PCL scaffold proposed in this study is a framework that offers great potentials for CTE applications. PMID:27328736
NASA Technical Reports Server (NTRS)
Colborn, B. L.; Armstong, T. W.
1993-01-01
A three-dimensional geometry and mass model of the Long Duration Exposure Facility (LDEF) spacecraft and experiment trays was developed for use in predictions and data interpretation related to ionizing radiation measurements. The modeling approach, level of detail incorporated, example models for specific experiments and radiation dosimeters, and example applications of the model are described.
Ambosta, Althea H.; Reichert, James F.; Kelly, Debbie M.
2013-01-01
Studies have shown that animals, including humans, use the geometric properties of environments to orient. It has been proposed that orientation is accomplished primarily by encoding the principal axes (i.e., global geometry) of an environment. However, recent research has shown that animals use local information such as wall length and corner angles as well as local shape parameters (i.e., medial axes) to orient. The goal of the current study was to determine whether adult humans reorient according to global geometry based on principal axes or whether reliance is on local geometry such as wall length and sense information or medial axes. Using a virtual environment task, participants were trained to select a response box located at one of two geometrically identical corners within a featureless rectangular-shaped environment. Participants were subsequently tested in a transformed L-shaped environment that allowed for a dissociation of strategies based on principal axes, medial axes and local geometry. Results showed that participants relied primarily on a medial axes strategy to reorient in the L-shaped test environment. Importantly, the search behaviour of participants could not be explained by a principal axes-based strategy. PMID:24223869
Dual anode contact geometries for x-ray and gamma-ray spectroscopy and 3D localization
NASA Astrophysics Data System (ADS)
Martin, J. W.; Garson, A. B., III; Li, Q.; Lee, K.; Groza, M.; Buliga, V.; Burger, A.; Krawczynski, H.
2009-08-01
We report on the continued development and testing of unique types of Cadmium Zinc Telluride (CZT) detectors. Using large volume (10×20×20 mm3) CZT crystals, we contact various "dual anode" detector designs. We incorporate a segmented cathode with five regions so that the charge on all seven contacts can be used to determine the energy and the 3-D interaction location of detected X-ray and gamma-ray photons. We describe the status of the detector development program, emphasize strengths and weaknesses of the different contact configurations, and discuss possible applications of Dual Anode Detectors in radiation detection applications.
Manukyan, Liana; Milinkovitch, Michel C.
2015-01-01
While recent imaging techniques provide insights into biological processes from the molecular to the cellular scale, phenotypes at larger scales remain poorly amenable to quantitative analyses. For example, investigations of the biophysical mechanisms generating skin morphological complexity and diversity would greatly benefit from 3D geometry and colour-texture reconstructions. Here, we report on R2OBBIE-3D, an integrated system that combines a robotic arm, a high-resolution digital colour camera, an illumination basket of high-intensity light-emitting diodes and state-of-the-art 3D-reconstruction approaches. We demonstrate that R2OBBIE generates accurate 3D models of biological objects between 1 and 100 cm, makes multiview photometric stereo scanning possible in practical processing times, and enables the capture of colour-texture and geometric resolutions better than 15 μm without the use of magnifying lenses. R2OBBIE has the potential to greatly improve quantitative analyses of phenotypes in addition to providing multiple new applications in, e.g., biomedical science. PMID:26039509
Martins, António F; Bessant, Michel; Manukyan, Liana; Milinkovitch, Michel C
2015-01-01
While recent imaging techniques provide insights into biological processes from the molecular to the cellular scale, phenotypes at larger scales remain poorly amenable to quantitative analyses. For example, investigations of the biophysical mechanisms generating skin morphological complexity and diversity would greatly benefit from 3D geometry and colour-texture reconstructions. Here, we report on R(2)OBBIE-3D, an integrated system that combines a robotic arm, a high-resolution digital colour camera, an illumination basket of high-intensity light-emitting diodes and state-of-the-art 3D-reconstruction approaches. We demonstrate that R(2)OBBIE generates accurate 3D models of biological objects between 1 and 100 cm, makes multiview photometric stereo scanning possible in practical processing times, and enables the capture of colour-texture and geometric resolutions better than 15 μm without the use of magnifying lenses. R(2)OBBIE has the potential to greatly improve quantitative analyses of phenotypes in addition to providing multiple new applications in, e.g., biomedical science. PMID:26039509
Shapiro, A.B.
1983-08-01
The computer code FACET calculates the radiation geometric view factor (alternatively called shape factor, angle factor, or configuration factor) between surfaces for axisymmetric, two-dimensional planar and three-dimensional geometries with interposed third surface obstructions. FACET was developed to calculate view factors for input to finite-element heat-transfer analysis codes. The first section of this report is a brief review of previous radiation-view-factor computer codes. The second section presents the defining integral equation for the geometric view factor between two surfaces and the assumptions made in its derivation. Also in this section are the numerical algorithms used to integrate this equation for the various geometries. The third section presents the algorithms used to detect self-shadowing and third-surface shadowing between the two surfaces for which a view factor is being calculated. The fourth section provides a user's input guide followed by several example problems.
NASA Astrophysics Data System (ADS)
Bai, Chao-ying; Huang, Guo-jiao; Li, Xing-wang; Greenhalgh, Stewart
2015-10-01
Traditionally, traveltime tomography entails inversion of either the velocity field and the reflector geometry sequentially, or the velocity field and the hypocenter locations simultaneously or in a cascaded fashion, but seldom are all three types (velocities, geometry of reflectors, and source locations) updated simultaneously because of the compromise between the different classes of model variable and the lack of different seismic phases to constrain these variables. By using a state-of-the-art ray-tracing algorithm for the first and later arrivals combined with a popular linearized inversion solver, it is possible to simultaneously recover the three classes of model variables. In the work discussed in this paper we combined the multistage irregular shortest-path ray-tracing algorithm with a subspace inversion solver to achieve simultaneous inversion of multi-class variables, using arrival times for different phases to concurrently obtain the velocity field, the reflector shapes, and the hypocenter locations. Simulation and comparison tests for two sets of source-receiver arrangements (one the ideal case and the other an approximated real case) indicate that the combined triple-class inversion algorithm is capable of obtaining nearly the same results as the double-class affect inversion scheme (velocity and reflector geometry, or velocity and source locations) even if a lower ray density and irregular source-receiver geometry are used to simulate the real situation. In addition, the new simultaneous inversion method is not sensitive to a modest amount of picking error in the traveltime data and reasonable uncertainty in earthquake hypocenter locations, which shows it to be a feasible and promising approach in real applications.
An IPOT meshless method using DC PSE approximation for fluid flow equations in 2D and 3D geometries
NASA Astrophysics Data System (ADS)
Bourantas, G. C.; Loukopoulos, V. C.; Skouras, E. D.; Burganos, V. N.; Nikiforidis, G. C.
2016-06-01
Navier-Stokes (N-S) equations, in their primitive variable (u-v-p) formulation, are numerically solved using the Implicit Potential (IPOT) numerical scheme in the context of strong form Meshless Point Collocation (MPC) method. The unknown field functions are computed using the Discretization Correction Particle Strength Exchange (DC PSE) approximation method. The latter makes use of discrete moment conditions to derive the operator kernels, which leads to low condition number for the moment matrix compared to other meshless interpolation methods and increased stability for the numerical solution. The proposed meshless scheme is applied on 2D and 3D spatial domains, using uniform or irregular set of nodes to represent the domain. The numerical results obtained are compared against those obtained using well-established methods.
NASA Astrophysics Data System (ADS)
Suleiman, Adamu; Jackson, Christopher; Magee, Craig; Fraser, Alastair
2016-04-01
Recent studies of regional unconformities in the circum-South Atlantic tectonic plates have linked unconformity age to the timing of changes in the azimuth of oceanic fracture zones, caused by plate interactions during opening of the South Atlantic. This observation is significant, proposing that a plate boundary geodynamic processes are transmitted into and expressed in plate interiors. However, it is not yet clear if and how other geologic events, such as intra-plate magmatism, may be linked to changes in the oceanic fracture azimuthal geometry. Here we use 2D and 3D seismic reflection, geochemical, borehole datasets and outcrop observations from the Bornu Basin, one of several intra-continental rift basins located in NE Nigeria to constrain the 3D geometry of igneous bodies and magmatic emplacement processes. This allows us to link South Atlantic plate boundary geodynamics and magmatism in the surrounding continental rift basins. Seismic attributes, reflection intensity, relative acoustic impedance, were used to identify and map igneous intrusions. Saucer-shaped sills are the most common type of intrusion, although en-echelon sills, up to 1.4 km in length, were also identified. The 3D geometry of the sills reveals the detailed structural components like inner sill, inclined sheets and outer sill. A mapped bifurcating network of the sills suggests magma emplacement process through upward and outward propagation. Seismic-stratigraphic observations indicate that igneous activity occurred in the Early Cretaceous, Late Cretaceous and Paleogene corresponding to the timing of major azimuth changes observed in the Kane Oceanic fracture zone in the South Atlantic Ocean. Overall, our study, suggests a possible influence of plate boundary geodynamics on intra-plate magmatism as reflected in the link between the time of changes in the azimuth of oceanic fracture zones and magmatic emplacement observed in the tectono-stratigraphy of the intra-continental rift basins.
NASA Astrophysics Data System (ADS)
Suleiman, A. A.; Jackson, C. A. L.; Magee, C.; Fraser, A.
2015-12-01
Recent studies of regional unconformities in the circum-South Atlantic tectonic plates have linked unconformity age to the timing of changes in the azimuth of oceanic fracture zones, caused by plate interactions during opening of the South Atlantic. This observation is significant, proposing that a plate boundary geodynamic processes are transmitted into and expressed in plate interiors. However, it is not yet clear if and how other geologic events, such as intra-plate magmatism, may be linked to changes in the oceanic fracture azimuthal geometry. Here we use 2D and 3D seismic reflection, geochemical, borehole datasets and outcrop observations from the Bornu Basin, one of several intra-continental rift basins located in NE Nigeria to constrain the 3D geometry of igneous bodies and magmatic emplacement processes. This allows us to link South Atlantic plate boundary geodynamics and magmatism in the surrounding continental rift basins. Seismic attributes, reflection intensity, relative acoustic impedance, were used to identify and map igneous intrusions. Saucer-shaped sills are the most common type of intrusion, although en-echelon sills, up to 1.4 km in length, were also identified. The 3D geometry of the sills reveals the detailed structural components like inner sill, inclined sheets and outer sill (Fig.1). A mapped bifurcating network of the sills suggests magma emplacement process through upward and outward propagation. Seismic-stratigraphic observations indicate that igneous activity occurred in the Early Cretaceous, Late Cretaceous and Paleogene corresponding to the timing of major azimuth changes observed in the Kane Oceanic fracture zone in the South Atlantic Ocean. Overall, our study, suggests a possible influence of plate boundary geodynamics on intra-plate magmatism as reflected in the link between the time of changes in the azimuth of oceanic fracture zones and magmatic emplacement observed in the tectono-stratigraphy of the intra-continental rift basins.
NASA Astrophysics Data System (ADS)
Lücke, O. H.; Arroyo, I. G.
2015-07-01
The eastern part of the oceanic Cocos Plate presents a heterogeneous crustal structure due to diverse origins and ages as well as plate-hot spot interactions which originated the Cocos Ridge, a structure that converges with the Caribbean Plate in southeastern Costa Rica. The complex structure of the oceanic plate directly influences the dynamics and geometry of the subduction zone along the Middle American Trench. In this paper an integrated interpretation of the slab geometry is presented based on three-dimensional density modeling of combined satellite and surface gravity data, constrained by available geophysical and geological data and seismological information obtained from local networks. The results show the continuation of steep subduction geometry from the Nicaraguan margin into Northwestern Costa Rica, followed by a moderate dipping slab under the Central Cordillera toward the end of the Central American Volcanic Arc. To the southeast end of the volcanic arc, our preferred model shows a steep, coherent slab that extends up to the landward projection of the Panama Fracture Zone. Overall, a gradual change in the depth of the intraplate seismicity is observed, reaching 220 km in the northwestern part, and becoming progressively shallower toward the southeast, where it reaches a terminal depth of 75 km. The changes in the terminal depth of the observed seismicity correlate with the increased density in the modeled slab. The absence of intermediate depth intraplate seismicity in the southeastern section and the higher densities for the subducted slab in this area, support a model in which dehydration reactions in the subducted slab cease at a shallower depth, originating an anhydrous and thus aseismic slab.
2009-01-01
Background Discovery of new bioactive molecules that could enter drug discovery programs or that could serve as chemical probes is a very complex and costly endeavor. Structure-based and ligand-based in silico screening approaches are nowadays extensively used to complement experimental screening approaches in order to increase the effectiveness of the process and facilitating the screening of thousands or millions of small molecules against a biomolecular target. Both in silico screening methods require as input a suitable chemical compound collection and most often the 3D structure of the small molecules has to be generated since compounds are usually delivered in 1D SMILES, CANSMILES or in 2D SDF formats. Results Here, we describe the new open source program DG-AMMOS which allows the generation of the 3D conformation of small molecules using Distance Geometry and their energy minimization via Automated Molecular Mechanics Optimization. The program is validated on the Astex dataset, the ChemBridge Diversity database and on a number of small molecules with known crystal structures extracted from the Cambridge Structural Database. A comparison with the free program Balloon and the well-known commercial program Omega generating the 3D of small molecules is carried out. The results show that the new free program DG-AMMOS is a very efficient 3D structure generator engine. Conclusion DG-AMMOS provides fast, automated and reliable access to the generation of 3D conformation of small molecules and facilitates the preparation of a compound collection prior to high-throughput virtual screening computations. The validation of DG-AMMOS on several different datasets proves that generated structures are generally of equal quality or sometimes better than structures obtained by other tested methods. PMID:19912625
NASA Astrophysics Data System (ADS)
Shahbazmohamadi, Sina; Jordan, Eric H.
2012-12-01
Creation of three-dimensional representations of surfaces from images taken at two or more view angles is a well-established technique applied to optical images and is frequently used in combination with scanning electron microscopy (SEM). The present work describes specific steps taken to optimize and enhance the repeatability of three-dimensional surfaces reconstructed from SEM images. The presented steps result in an approximately tenfold improvement in the repeatability of the surface reconstruction compared to more standard techniques. The enhanced techniques presented can be used with any SEM friendly samples. In this work the modified technique was developed in order to accurately quantify surface geometry changes in metallic bond coats used with thermal barrier coatings (TBCs) to provide improved turbine hot part durability. Bond coat surfaces are quite rough, and accurate determination of surface geometry change (rumpling) requires excellent repeatability. Rumpling is an important contributor to TBC failure, and accurate quantification of rumpling is important to better understanding of the failure behavior of TBCs.
NASA Astrophysics Data System (ADS)
Silva, D. J.; Wahl, U.; Correia, J. G.; Augustyns, V.; Lima, T. A. L.; Costa, A.; Bosne, E.; da Silva, M. R.; Araújo, J. P.; Pereira, L. M. C.
2016-03-01
Although the formation of transition metal-boron pairs is currently well established in silicon processing, the geometry of these complexes is still not completely understood. We investigated the lattice location of the transition metals manganese, iron, cobalt and nickel in n- and p+ -type silicon by means of electron emission channeling. For manganese, iron and cobalt, we observed an increase of sites near the ideal tetrahedral interstitial position by changing the doping from n- to p+ -type Si. Such increase was not observed for Ni. We ascribe this increase to the formation of pairs with boron, driven by Coulomb interactions, since the majority of iron, manganese and cobalt is positively charged in p+ -type silicon while Ni is neutral. We propose that breathing mode relaxation around the boron ion within the pair causes the observed displacement from the ideal tetrahedral interstitial site. We discuss the application of the emission channeling technique in this system and, in particular, how it provides insight on the geometry of such pairs.
NASA Astrophysics Data System (ADS)
Aquè, R.; Tavarnelli, E.
2012-04-01
The three-dimensional (3D) reconstruction of complex geological settings and of original, pre-thrusting basin geometry is one of the challenges for modern structural geology. It has indeed a critical role in many industrial applications, such as in the hydrocarbon exploration. By using commercial specific softwares to produce balanced cross-sections and inferred 3D reconstructions (2DMove™, Gocad™), we modelled a portion of the Umbria-Marche fold-and-thrust belt, in the outer zones of the Northern Apennines of Italy, in order to infer the pre-thrusting geometry of the Mesozoic-Cenozoic extensional basins and to test the applicability of existing computer tools in areas that have experienced the effects of positive tectonic inversion. In the study area, the accurate reconstruction of the structural setting, cross-cut relationships and timing of the deformation, was inferred by using field data, map analysis and cross-section balancing techniques. The structural overprinting relationships among the investigated thrusts made it possible to infer a general piggy-back thrusting sequence, with new thrust faults to the East, developed in the footwall of formerly emplaced thrust sheets, in the West. This allowed to sequentially remove the effects of the deformation for progressively older structures, and to back-strip the thrust sheets in sequential evolutionary steps, in order to reconstruct a viable pre-thrusting template. Four balanced cross-sections have been drawn, providing the initial skeleton for 3D modelling, together with the map trace of the major tectonic features. The cross-sections and the geological map have been digitized and geo-referred in 2D-Move™. Starting from the inferred geometries, a coherent 3D model was built in Gocad™. The surfaces represent post-thrust normal faults, thrust planes, and pre-thrust normal faults, and five key stratigraphic surfaces, from bottom; the base and top of the Calcare Massiccio fm. (Lower Liassic), the base of the
NASA Astrophysics Data System (ADS)
Dalichaouch, Thamine; Yu, Peicheng; Davidson, Asher; Mori, Warren; Vieira, Jorge; Fonseca, Ricardo
2015-11-01
Laser wakefield acceleration (LWFA) has attracted a lot of interest as a possible compact particle accelerator. However, 3D simulations of plasma-based accelerators are computationally intensive, sometimes taking millions of core hours on today's computers. A quasi-3D particle-In-cell (PIC) approach has been developed to take advantage of azimuthal symmetry in LWFA (and PWFA) simulations by using a particle-in-cell description in r-z and a Fourier description in φ. Quasi-3D simulations of LWFA are computationally more efficient and faster than Full-3D simulations because only first few azimuthal harmonics are needed to capture the physics of the problem. We have developed a cylindrical mode decomposition diagnostic for 3D Cartesian geometry simulations to analyze the agreement between full-3D and quasi-3D PIC simulations of laser and beam-plasma interactions. The diagnostic interpolates field data from Full-3D PIC simulations onto an irregular cylindrical grid (r , φ , z). A Fourier decomposition is then performed on the interpolated 3D simulation data along the azimuthal direction. This diagnostic has the added advantage of separating out the wakefields from the laser field. Preliminary results for this diagnostic of LWFA and PWFA simulations with symmetric and nearly symmetric spot sizes as well as of laser-plasma interactions using lasers with orbital angular momentum (higher order Laguerre-Gaussian modes) will be presented.
Speidel, M; Hatt, C; Tomkowiak, M; Raval, A; Funk, T
2014-06-15
Purpose: To develop a method for the fusion of 3D echocardiography and Scanning-Beam Digital X-ray (SBDX) fluoroscopy to assist with catheter device and soft tissue visualization during interventional procedures. Methods: SBDX is a technology for low-dose inverse geometry x-ray fluoroscopy that performs digital tomosynthesis at multiple planes in real time. In this study, transesophageal echocardiography (TEE) images were fused with SBDX images by estimating the 3D position and orientation (the “pose”) of the TEE probe within the x-ray coordinate system and then spatially transforming the TEE image data to match this pose. An initial pose estimate was obtained through tomosynthesis-based 3D localization of points along the probe perimeter. Position and angle estimates were then iteratively refined by comparing simulated projections of a 3D probe model against SBDX x-ray images. Algorithm performance was quantified by imaging a TEE probe in different known orientations and locations within the x-ray field (0-30 degree tilt angle, up to 50 mm translation). Fused 3D TEE/SBDX imaging was demonstrated by imaging a tissue-mimicking polyvinyl alcohol cylindrical cavity as a catheter was navigated along the cavity axis. Results: Detected changes in probe tilt angle agreed with the known changes to within 1.2 degrees. For a 50 mm translation along the source-detector axis, the detected translation was 50.3 mm. Errors for in-plane translations ranged from 0.1 to 0.9 mm. In a fused 3D TEE/SBDX display, the catheter device was well visualized and coincident with the device shadow in the TEE images. The TEE images portrayed phantom boundaries that were not evident under x-ray. Conclusion: Registration of soft tissue anatomy derived from TEE imaging and device imaging from SBDX x-ray fluoroscopy is feasible. The simultaneous 3D visualization of these two modalities may be useful in interventional procedures involving the navigation of devices to soft tissue anatomy.
Cullen, D E
1998-11-22
TART98 is a coupled neutron-photon, 3 Dimensional, combinatorial geometry, time dependent Monte Carlo radiation transport code. This code can run on any modern computer. It is a complete system to assist you with input preparation, running Monte Carlo calculations, and analysis of output results. TART98 is also incredibly FAST; if you have used similar codes, you will be amazed at how fast this code is compared to other similar codes. Use of the entire system can save you a great deal of time and energy. TART98 is distributed on CD. This CD contains on-line documentation for all codes included in the system, the codes configured to run on a variety of computers, and many example problems that you can use to familiarize yourself with the system. TART98 completely supersedes all older versions of TART, and it is strongly recommended that users only use the most recent version of TART98 and its data files.
High Energy Boundary Conditions for a Cartesian Mesh Euler Solver
NASA Technical Reports Server (NTRS)
Pandya, Shishir; Murman, Scott; Aftosmis, Michael
2003-01-01
Inlets and exhaust nozzles are common place in the world of flight. Yet, many aerodynamic simulation packages do not provide a method of modelling such high energy boundaries in the flow field. For the purposes of aerodynamic simulation, inlets and exhausts are often fared over and it is assumed that the flow differences resulting from this assumption are minimal. While this is an adequate assumption for the prediction of lift, the lack of a plume behind the aircraft creates an evacuated base region thus effecting both drag and pitching moment values. In addition, the flow in the base region is often mis-predicted resulting in incorrect base drag. In order to accurately predict these quantities, a method for specifying inlet and exhaust conditions needs to be available in aerodynamic simulation packages. A method for a first approximation of a plume without accounting for chemical reactions is added to the Cartesian mesh based aerodynamic simulation package CART3D. The method consists of 3 steps. In the first step, a components approach where each triangle is assigned a component number is used. Here, a method for marking the inlet or exhaust plane triangles as separate components is discussed. In step two, the flow solver is modified to accept a reference state for the components marked inlet or exhaust. In the third step, the flow solver uses these separated components and the reference state to compute the correct flow condition at that triangle. The present method is implemented in the CART3D package which consists of a set of tools for generating a Cartesian volume mesh from a set of component triangulations. The Euler equations are solved on the resulting unstructured Cartesian mesh. The present methods is implemented in this package and its usefulness is demonstrated with two validation cases. A generic missile body is also presented to show the usefulness of the method on a real world geometry.
NASA Astrophysics Data System (ADS)
Dutta, Rishabh; Jónsson, Sigurjón
2016-04-01
Earthquake faults are generally considered planar (or of other simple geometry) in earthquake source parameter estimations. However, simplistic fault geometries likely result in biases in estimated slip distributions and increased fault slip uncertainties. In case of large subduction zone earthquakes, these biases and uncertainties propagate into tsunami waveform modeling and other calculations related to postseismic studies, Coulomb failure stresses, etc. In this research, we parameterize 3D non-planar fault geometry for the 2011 Tohoku-Oki earthquake (Mw 9.1) and estimate these geometrical parameters along with fault slip parameters from onland and offshore GPS using Bayesian inference. This non-planar fault is formed using several 3rd degree polynomials in along-strike (X-Y plane) and along-dip (X-Z plane) directions that are tied together using a triangular mesh. The coefficients of these polynomials constitute the fault geometrical parameters. We use the trench and locations of past seismicity as a priori information to constrain these fault geometrical parameters and the Laplacian to characterize the fault slip smoothness. Hyper-parameters associated to these a priori constraints are estimated empirically and the posterior probability distribution of the model (fault geometry and slip) parameters is sampled using an adaptive Metropolis Hastings algorithm. The across-strike uncertainties in the fault geometry (effectively the local fault location) around high-slip patches increases from 6 km at 10km depth to about 35 km at 50km depth, whereas around low-slip patches the uncertainties are larger (from 7 km to 70 km). Uncertainties in reverse slip are found to be higher at high slip patches than at low slip patches. In addition, there appears to be high correlation between adjacent patches of high slip. Our results demonstrate that we can constrain complex non-planar fault geometry together with fault slip from GPS data using past seismicity as a priori
NASA Astrophysics Data System (ADS)
Daehne, A.; Travelletti, J.; Malet, J.-P.; Corsini, A.; Ronchetti, F.
2009-04-01
Bedrock geometry drastically influences the kinematic deformation pattern of slow-moving landslides exhibiting some flow characteristics. The development of extension and compression zones within the landslide body is largely controlled by the geometry (crests, bumps, hollows) and roughness of the topography covered by the moving mass. A challenge to progress in the forecast of such type of landslides is to precisely define 3D geometrical and geomechanical models. The objective of this work is to present a methodology for 3D geometrical modelling of the landslide structure, and to discuss the main possible errors in integrating multi-source and multi-resolution data in the modelling. The methodology is presented through the analysis of three landslides with similar geomorphological features (e.g. flow-like geomorphology) and development patterns (retrogression of the crown, roto-translational failures of the upper part, and translational movements in the lower part), and for which an extensive dataset of geophysical, geotechnical and geomorphological information is available. The three cases studies are the complex Valoria earth-slide-flow located in the Northern Apennines, the Super-Sauze and La Valette mudslides in the French South Alps. All three landslides are predominantly developed in a clay-shale soil formation. First, interpretation of the multi-data information, their resolution and accuracy is presented for the landslides. Second, a procedure to construct 3D geometrical models of the landslides is proposed (by using the Rockware's Rockworks geological modeller) and the influence of the interpolation algorithms is discussed. It is demonstrated that the model uncertainty is strongly depending on the density and distribution of the input data which vary for the three landslides. The quality of several geometrical models is then compared; a best-fit is achieved by using available geological and geomorphological site interpretation.
NASA Astrophysics Data System (ADS)
Bailey, B.; Stoll, R., II; Miller, N. E.; Pardyjak, E.; Mahaffee, W.
2014-12-01
Plants cover the majority of Earth's land surface, and thus play a critical role in the surface energy balance. Within individual plant communities, the leaf energy balance is a fundamental component of most biophysical processes. Absorbed radiation drives the energy balance and provides the means by which plants produce food. Available energy is partitioned into sensible and latent heat fluxes to determine surface temperature, which strongly influences rates of metabolic activity and growth. The energy balance of an individual leaf is coupled with other leaves in the community through longwave radiation emission and advection through the air. This complex coupling can make scaling models from leaves to whole-canopies difficult, specifically in canopies with complex, heterogeneous geometries. We present a new three-dimensional canopy model that simultaneously resolves sub-tree to whole-canopy scales. The model provides spatially explicit predictions of net radiation exchange, boundary-layer and stomatal conductances, evapotranspiration rates, and ultimately leaf surface temperature. The radiation model includes complex physics such as anisotropic emission and scattering. Radiation calculations are accelerated by leveraging graphics processing unit (GPU) technology, which allows canopy-scale problems to be performed on a standard desktop workstation. Since validating the three-dimensional distribution of leaf temperature can be extremely challenging, we used several independent measurement techniques to quantify errors in measured and modeled values. When compared with measured leaf temperatures, the model gave a mean error of about 2°C, which was close to the estimated measurement uncertainty.
Luo, Jiajia; Betschart, Cornelia; Chen, Luyun; Ashton-Miller, James A.; DeLancey, John O. L.
2014-01-01
Introduction and hypothesis A method was developed using 3D stress magnetic resonance imaging (MRI) and was piloted to test hypotheses concerning changes in apical ligament lengths and lines of action from rest to maximal Valsalva. Methods Ten women with (cases) and ten without (controls) pelvic organ prolapse (POP) were selected from an ongoing case-control study. Supine, multiplanar stress MRI was performed at rest and at maximal Valsalva and was imported into 3D Slicer v. 3.4.1 and aligned. The 3D reconstructions of the uterus and vagina, cardinal ligament (CL), deep uterosacral ligament (USLd), and pelvic bones were created. Ligament length and orientation were then measured. Results Adequate ligament representations were possible in all 20 study participants. When cases were compared with controls, the curve length of the CL at rest was 71 ±16 mm vs. 59±9 mm (p =0.051), and the USLd was 38±16 mm vs. 36±11 mm (p =0.797). Similarly, the increase in CL length from rest to strain was 30±16 mm vs. 15±9 mm (p =0.033), and USLd was 15±12 mm vs. 7±4 mm (p =0.094). Likewise, the change in USLd angle was significantly different from CL (p <0.001). Conclusions This technique allows quantification of 3D geometry at rest and at strain. In our pilot sample, at maximal Valsalva, CL elongation was greater in cases than controls, whereas USLd was not; CL also exhibited greater changes in ligament length, and USLd exhibited greater changes in ligament inclination angle. PMID:24008367
Modelling and 3D optimisation of CdTe pixels detector array geometry - Extension to small pixels
NASA Astrophysics Data System (ADS)
Zumbiehl, A.; Hage-Ali, M.; Fougeres, P.; Koebel, J. M.; Regal, R.; Rit, C.; Ayoub, M.; Siffert, P.
2001-08-01
CdTe and CdZnTe pixel detectors offer great interest for many applications, especially for medical and industrial imaging. Up to now, the material, generally, used and investigated for pixel arrays was CZT (Hamel et al., IEEE Trans. Nucl. Sci. 43 (3) (1996) 1422; Barrett et al., Phys. Rev. Lett. 75 (1) (1995) 156; Bennett et al., Nucl. Instr. and Meth. A 392 (1997) 260; Eskin et al., J. Appl. Phys. 85 (2) (1999) 647; Brunett et al., J. Appl. Phys. 86 (7) (1999) 3926; Luke, Nucl. Instr. and Meth. A 380 (1996) 232), but cadmium telluride can also be an appropriate choice, as shown here. However, we clearly demonstrate here that the optimal pixel configuration is highly dependent on the electrical transport properties of the material. Depending on the field of primary interest, either energy resolution or counting rate efficiency in the photopeak, the geometry for each case has to be optimised. For that purpose, we have developed a calculation of the signal induced onto the pixel. Two distinct parts are used: after showing our approach for the weighting potential calculation, we present our results performed by a "pseudo-Monte Carlo" simulation. Results are supported by a few experimental comparisons. We argue about the optimum sizes with clarifying the problems caused by too small and too large pixel sizes. The study field is chosen to be vast, i.e. pixel size to detector thickness ratios ( W/ L) of 1/8-1, and detector thickness of 1.0-8.0 mm. In addition, several electrical transport properties are used. Since efficiency is often of primary interest, thick detectors could be very attractive, which are shown to be really feasible even on CdTe.
Static Aeroelastic Analysis with an Inviscid Cartesian Method
NASA Technical Reports Server (NTRS)
Rodriguez, David L.; Aftosmis, Michael J.; Nemec, Marian; Smith, Stephen C.
2014-01-01
An embedded-boundary Cartesian-mesh flow solver is coupled with a three degree-offreedom structural model to perform static, aeroelastic analysis of complex aircraft geometries. The approach solves the complete system of aero-structural equations using a modular, loosely-coupled strategy which allows the lower-fidelity structural model to deform the highfidelity CFD model. The approach uses an open-source, 3-D discrete-geometry engine to deform a triangulated surface geometry according to the shape predicted by the structural model under the computed aerodynamic loads. The deformation scheme is capable of modeling large deflections and is applicable to the design of modern, very-flexible transport wings. The interface is modular so that aerodynamic or structural analysis methods can be easily swapped or enhanced. This extended abstract includes a brief description of the architecture, along with some preliminary validation of underlying assumptions and early results on a generic 3D transport model. The final paper will present more concrete cases and validation of the approach. Preliminary results demonstrate convergence of the complete aero-structural system and investigate the accuracy of the approximations used in the formulation of the structural model.
Stable boundary conditions for Cartesian grid calculations
NASA Technical Reports Server (NTRS)
Berger, M. J.; Leveque, R. J.
1990-01-01
The inviscid Euler equations in complicated geometries are solved using a Cartesian grid. This requires solid wall boundary conditions in the irregular grid cells near the boundary. Since these cells may be orders of magnitude smaller than the regular grid cells, stability is a primary concern. An approach to this problem is presented and its use is illustrated.
NASA Astrophysics Data System (ADS)
Hong, Daeki; Cho, Heemoon; Cho, Hyosung; Choi, Sungil; Je, Uikyu; Park, Yeonok; Park, Chulkyu; Lim, Hyunwoo; Park, Soyoung; Woo, Taeho
2015-11-01
In this work, we performed a feasibility study on the three-dimensional (3D) image reconstruction in a truncated Archimedean-like spiral geometry with a long-rectangular detector for application to high-accurate, cost-effective dental x-ray imaging. Here an x-ray tube and a detector rotate together around the rotational axis several times and, concurrently, the detector moves horizontally in the detector coordinate at a constant speed to cover the whole imaging volume during the projection data acquisition. We established a table-top setup which mainly consists of an x-ray tube (60 kVp, 5 mA), a narrow CMOS-type detector (198-μm pixel resolution, 184 (W)×1176 (H) pixel dimension), and a rotational stage for sample mounting and performed a systematic experiment to demonstrate the viability of the proposed approach to volumetric dental imaging. For the image reconstruction, we employed a compressed-sensing (CS)-based algorithm, rather than a common filtered-backprojection (FBP) one, for more accurate reconstruction. We successfully reconstructed 3D images of considerably high quality and investigated the image characteristics in terms of the image value profile, the contrast-to-noise ratio (CNR), and the spatial resolution.
NASA Astrophysics Data System (ADS)
Rapoport, Diego L.
2011-01-01
In this transdisciplinary article which stems from philosophical considerations (that depart from phenomenology—after Merleau-Ponty, Heidegger and Rosen—and Hegelian dialectics), we develop a conception based on topological (the Moebius surface and the Klein bottle) and geometrical considerations (based on torsion and non-orientability of manifolds), and multivalued logics which we develop into a unified world conception that surmounts the Cartesian cut and Aristotelian logic. The role of torsion appears in a self-referential construction of space and time, which will be further related to the commutator of the True and False operators of matrix logic, still with a quantum superposed state related to a Moebius surface, and as the physical field at the basis of Spencer-Brown's primitive distinction in the protologic of the calculus of distinction. In this setting, paradox, self-reference, depth, time and space, higher-order non-dual logic, perception, spin and a time operator, the Klein bottle, hypernumbers due to Musès which include non-trivial square roots of ±1 and in particular non-trivial nilpotents, quantum field operators, the transformation of cognition to spin for two-state quantum systems, are found to be keenly interwoven in a world conception compatible with the philosophical approach taken for basis of this article. The Klein bottle is found not only to be the topological in-formation for self-reference and paradox whose logical counterpart in the calculus of indications are the paradoxical imaginary time waves, but also a classical-quantum transformer (Hadamard's gate in quantum computation) which is indispensable to be able to obtain a complete multivalued logical system, and still to generate the matrix extension of classical connective Boolean logic. We further find that the multivalued logic that stems from considering the paradoxical equation in the calculus of distinctions, and in particular, the imaginary solutions to this equation
Lyras, Dimitrios N.; Loucks, Craig; Greenhow, Robert
2016-01-01
Background: The aim of this study is to evaluate the geometry of the distal femur and the proximal tibia in the osteoarthritic knee using 3D reconstructive CT scan imaging. Methods: 449 patients with knee osteoarthritis were treated surgically in our center with patient-specific technology total knee arthroplasty. Preoperatively, all the patients underwent a CT scan according to a standard protocol. Using this database, the Hip-Knee-Angle (HKA), the Femur Valgus Angle (FVA), the Tibia Varus Angle (TVA), the Posterior Tibia Slope (PTS), and the angle between the posterior condylar axis and the anatomical transepicondylar axis (PCA) for each patient were recorded and statistically evaluated. Results: In overall, the mean HKA angle was 177.3±5.55, the mean FVA angle was 3.19±2.08, the mean TVA was 3.28±2.35, the PTS angle was 9.02±3.46, and the PCA angle was 2.86±0.78. Evaluation of the correlations between HKA and PCA (r=0.035), HKA and PTS (r=-0.047), and PCA and PTS (r=0.05) showed non-significant relationships (P=0.46, P=0.32, and P=0.29 respectively). No significant differences were revealed from the comparison of male patients with female patients, regarding the mean HKA, FVA, TVA, PTS, and PCA. Conclusion: The posterior condylar axis is a well-defined but not a reliable axis, while the transepicondylar and the anteroposterior are reliable, but not easily defined axes. Given the large ranges and standard deviations of the location of posterior condylar axis, and the important inter- and intraobserver variability in the intraoperative location of the transepicondylar and the anteroposterior axes, the use of a preoperative 3D CT scan is recommended. PMID:27200387
NASA Astrophysics Data System (ADS)
Phillips, Thomas; Jackson, Christopher; Bell, Rebecca; Duffy, Oliver; Fossen, Haakon
2015-04-01
Rift basins form within lithosphere containing a range of heterogeneities, such as thin-skinned thrust belts and larger scale structures such as thick-skinned shear zones or crustal sutures. How the presence and reactivation of these structures during later rift events affect the geometry and evolution of rifts remains poorly understood as they are not typically well imaged on seismic data. The main reasons for this are that crystalline basement is often buried beneath thick sedimentary successions and contains small impedance contrasts. Furthermore, larger, crustal-scale, lineaments and sutures may not be imaged at all on seismic data due to their large scale and depth. In this study, we use borehole-constrained 2D and 3D seismic reflection data located around the Egersund and Farsund Basins, offshore south Norway. In both areas, crystalline basement is exceptionally well-imaged on typical 2D and 3D reflection data due to large impedance contrasts within a highly heterogeneous, shallow basement. This allows us to map a series of intrabasement reflections and overlying rift systems. Within the Egersund area, two main types of intrabasement structure are identified and mapped: i) thin (100 m), shallowly dipping (0-10°W) reflections showing a ramp-flat geometry; and ii) thick (1-1.5 km), low angle (c. 30°W) structures comprising of packages of reflections. These structures correlate along-strike northwards to Caledonian orogeny related structures mapped onshore Norway. The thin structures are interpreted as thin-skinned Caledonian thrusts, whereas the thicker structures represent thick-skinned Devonian shear zones formed through orogenic collapse of the Caledonides. Through seismic-stratigraphic analysis of the cover, we document multiple stages of extensional reactivation along these structures during Devonian, Permian-Triassic and Late Jurassic-Early Cretaceous extension followed by reverse reactivation during Late Cretaceous compression. The Farsund Basin is
NASA Astrophysics Data System (ADS)
Gerke, Kirill; Vasilyev, Roman; Khirevich, Siarhei; Karsanina, Marina; Collins, Daniel; Korost, Dmitry; Mallants, Dirk
2015-04-01
In this contribution we introduce a novel free software which solves the Stokes equation to obtain velocity fields for low Reynolds-number flows within externally generated 3D pore geometries. Provided with velocity fields, one can calculate permeability for known pressure gradient boundary conditions via Darcy's equation. Finite-difference schemes of 2nd and 4th order of accuracy are used together with an artificial compressibility method to iteratively converge to a steady-state solution of Stokes' equation. This numerical approach is much faster and less computationally demanding than the majority of open-source or commercial softwares employing other algorithms (finite elements/volumes, lattice Boltzmann, etc.) The software consists of two parts: 1) a pre and post-processing graphical interface, and 2) a solver. The latter is efficiently parallelized to use any number of available cores (the speedup on 16 threads was up to 10-12 depending on hardware). Due to parallelization and memory optimization our software can be used to obtain solutions for 300x300x300 voxels geometries on modern desktop PCs. The software was successfully verified by testing it against lattice Boltzmann simulations and analytical solutions. To illustrate the software's applicability for numerous problems in Earth Sciences, a number of case studies have been developed: 1) identifying the representative elementary volume for permeability determination within a sandstone sample, 2) derivation of permeability/hydraulic conductivity values for rock and soil samples and comparing those with experimentally obtained values, 3) revealing the influence of the amount of fine-textured material such as clay on filtration properties of sandy soil. This work was partially supported by RSF grant 14-17-00658 (pore-scale modelling) and RFBR grants 13-04-00409-a and 13-05-01176-a.
Fevotte, F.; Lathuiliere, B.
2013-07-01
The large increase in computing power over the past few years now makes it possible to consider developing 3D full-core heterogeneous deterministic neutron transport solvers for reference calculations. Among all approaches presented in the literature, the method first introduced in [1] seems very promising. It consists in iterating over resolutions of 2D and ID MOC problems by taking advantage of prismatic geometries without introducing approximations of a low order operator such as diffusion. However, before developing a solver with all industrial options at EDF, several points needed to be clarified. In this work, we first prove the convergence of this iterative process, under some assumptions. We then present our high-performance, parallel implementation of this algorithm in the MICADO solver. Benchmarking the solver against the Takeda case shows that the 2D-1D coupling algorithm does not seem to affect the spatial convergence order of the MOC solver. As for performance issues, our study shows that even though the data distribution is suited to the 2D solver part, the efficiency of the ID part is sufficient to ensure a good parallel efficiency of the global algorithm. After this study, the main remaining difficulty implementation-wise is about the memory requirement of a vector used for initialization. An efficient acceleration operator will also need to be developed. (authors)
A Cartesian scheme for compressible multimaterial models in 3D
NASA Astrophysics Data System (ADS)
de Brauer, Alexia; Iollo, Angelo; Milcent, Thomas
2016-05-01
We model the three-dimensional interaction of compressible materials separated by sharp interfaces. We simulate fluid and hyperelastic solid flows in a fully Eulerian framework. The scheme is the same for all materials and can handle large deformations and frictionless contacts. Necessary conditions for hyperbolicity of the hyperelastic neohookean model in three dimensions are proved thanks to an explicit computation of the characteristic speeds. We present stiff multimaterial interactions including air-helium and water-air shock interactions, projectile-shield impacts in air and rebounds.
Static Aeroelastic Analysis with an Inviscid Cartesian Method
NASA Technical Reports Server (NTRS)
Rodriguez, David L.; Aftosmis, Michael J.; Nemec, Marian; Smith, Stephen C.
2014-01-01
An embedded-boundary, Cartesian-mesh flow solver is coupled with a three degree-of-freedom structural model to perform static, aeroelastic analysis of complex aircraft geometries. The approach solves a nonlinear, aerostructural system of equations using a loosely-coupled strategy. An open-source, 3-D discrete-geometry engine is utilized to deform a triangulated surface geometry according to the shape predicted by the structural model under the computed aerodynamic loads. The deformation scheme is capable of modeling large deflections and is applicable to the design of modern, very-flexible transport wings. The coupling interface is modular so that aerodynamic or structural analysis methods can be easily swapped or enhanced. After verifying the structural model with comparisons to Euler beam theory, two applications of the analysis method are presented as validation. The first is a relatively stiff, transport wing model which was a subject of a recent workshop on aeroelasticity. The second is a very flexible model recently tested in a low speed wind tunnel. Both cases show that the aeroelastic analysis method produces results in excellent agreement with experimental data.
NASA Astrophysics Data System (ADS)
Koehl, M.; Fabre, Ph.; Schlussel, B.
2014-06-01
Turckheim is a small town located in Alsace, north-east of France. In the heart of the Alsatian vineyard, this city has many historical monuments including its old church. To understand the effectiveness of the project described in this paper, it is important to have a look at the history of this church. Indeed there are many historical events that explain its renovation and even its partial reconstruction. The first mention of a christian sanctuary in Turckheim dates back to 898. It will be replaced in the 12th century by a roman church (chapel), which subsists today as the bell tower. Touched by a lightning in 1661, the tower then was enhanced. In 1736, it was repaired following damage sustained in a tornado. In 1791, the town installs an organ to the church. Last milestone, the church is destroyed by fire in 1978. The organ, like the heart of the church will then have to be again restored (1983) with a simplified architecture. From this heavy and rich past, it unfortunately and as it is often the case, remains only very few documents and information available apart from facts stated in some sporadic writings. And with regard to the geometry, the positioning, the physical characteristics of the initial building, there are very little indication. Some assumptions of positions and right-of-way were well issued by different historians or archaeologists. The acquisition and 3D modeling project must therefore provide the current state of the edifice to serve as the basis of new investigations and for the generation of new hypotheses on the locations and historical shapes of this church and its original chapel (Fig. 1)
Blodwell, J.F.
1987-10-01
It is argued that the point structure of space and time must be constructed from the primitive extensional character of space and time. A procedure for doing this is laid down and applied to one-dimensional and two-dimensional systems of abstract extensions. Topological and metrical properties of the constructed point systems, which differ nontrivially from the usual R and R/sup 2/ models, are examined. Briefly, constructed points are associated with directions and the Cartesian point is split. In one-dimension each point splits into a point pair compatible with the linear ordering. An application to one-dimensional particle motion is given, with the result that natural topological assumptions force the number of left point, right point transitions to remain locally finite in a continuous motion. In general, Cartesian points are seen to correspond to certain filters on a suitable Boolean algebra. Constructed points correspond to ultrafilters. Thus, point construction gives a natural refinement of the Cartesian systems.
Software for Automated Generation of Cartesian Meshes
NASA Technical Reports Server (NTRS)
Aftosmis, Michael J.; Melton, John E.; Berger, Marshal J.
2006-01-01
Cart3D is a collection of computer programs for generating Cartesian meshes [for computational fluid dynamics (CFD) and other applications] in volumes bounded by solid objects. Aspects of Cart3D at earlier stages of development were reported in "Robust and Efficient Generation of Cartesian Meshes for CFD" (ARC-14275), NASA Tech Briefs, Vol. 23, No. 8 (August 1999), page 30. The geometric input to Cart3D comprises surface triangulations like those commonly generated by computer-aided-design programs. Complexly shaped objects can be represented as assemblies of simpler ones. Cart3D deletes all portions of such an assembled object that are not on the exterior surface. Intersections between components are preserved in the resulting triangulation. A tie-breaking routine unambiguously resolves geometric degeneracies. Then taking the intersected surface triangulation as input, the volume mesh is generated through division of cells of an initially coarse hexahedral grid. Cells are subdivided to refine the grid in regions of increased surface curvature and/or increased flow gradients. Cells that become split into multiple unconnected regions by thin pieces of surface are identified.
Electronic Absolute Cartesian Autocollimator
NASA Technical Reports Server (NTRS)
Leviton, Douglas B.
2006-01-01
An electronic absolute Cartesian autocollimator performs the same basic optical function as does a conventional all-optical or a conventional electronic autocollimator but differs in the nature of its optical target and the manner in which the position of the image of the target is measured. The term absolute in the name of this apparatus reflects the nature of the position measurement, which, unlike in a conventional electronic autocollimator, is based absolutely on the position of the image rather than on an assumed proportionality between the position and the levels of processed analog electronic signals. The term Cartesian in the name of this apparatus reflects the nature of its optical target. Figure 1 depicts the electronic functional blocks of an electronic absolute Cartesian autocollimator along with its basic optical layout, which is the same as that of a conventional autocollimator. Referring first to the optical layout and functions only, this or any autocollimator is used to measure the compound angular deviation of a flat datum mirror with respect to the optical axis of the autocollimator itself. The optical components include an illuminated target, a beam splitter, an objective or collimating lens, and a viewer or detector (described in more detail below) at a viewing plane. The target and the viewing planes are focal planes of the lens. Target light reflected by the datum mirror is imaged on the viewing plane at unit magnification by the collimating lens. If the normal to the datum mirror is parallel to the optical axis of the autocollimator, then the target image is centered on the viewing plane. Any angular deviation of the normal from the optical axis manifests itself as a lateral displacement of the target image from the center. The magnitude of the displacement is proportional to the focal length and to the magnitude (assumed to be small) of the angular deviation. The direction of the displacement is perpendicular to the axis about which the
The program FANS-3D (finite analytic numerical simulation 3-dimensional) and its applications
NASA Technical Reports Server (NTRS)
Bravo, Ramiro H.; Chen, Ching-Jen
1992-01-01
In this study, the program named FANS-3D (Finite Analytic Numerical Simulation-3 Dimensional) is presented. FANS-3D was designed to solve problems of incompressible fluid flow and combined modes of heat transfer. It solves problems with conduction and convection modes of heat transfer in laminar flow, with provisions for radiation and turbulent flows. It can solve singular or conjugate modes of heat transfer. It also solves problems in natural convection, using the Boussinesq approximation. FANS-3D was designed to solve heat transfer problems inside one, two and three dimensional geometries that can be represented by orthogonal planes in a Cartesian coordinate system. It can solve internal and external flows using appropriate boundary conditions such as symmetric, periodic and user specified.
NASA Astrophysics Data System (ADS)
Ouadah, S.; Stayman, J. W.; Gang, G.; Uneri, A.; Ehtiati, T.; Siewerdsen, J. H.
2015-03-01
Purpose: Robotic C-arm systems are capable of general noncircular orbits whose trajectories can be driven by the particular imaging task. However obtaining accurate calibrations for reconstruction in such geometries can be a challenging problem. This work proposes a method to perform a unique geometric calibration of an arbitrary C-arm orbit by registering 2D projections to a previously acquired 3D image to determine the transformation parameters representing the system geometry. Methods: Experiments involved a cone-beam CT (CBCT) bench system, a robotic C-arm, and three phantoms. A robust 3D-2D registration process was used to compute the 9 degree of freedom (DOF) transformation between each projection and an existing 3D image by maximizing normalized gradient information with a digitally reconstructed radiograph (DRR) of the 3D volume. The quality of the resulting "self-calibration" was evaluated in terms of the agreement with an established calibration method using a BB phantom as well as image quality in the resulting CBCT reconstruction. Results: The self-calibration yielded CBCT images without significant difference in spatial resolution from the standard ("true") calibration methods (p-value >0.05 for all three phantoms), and the differences between CBCT images reconstructed using the "self" and "true" calibration methods were on the order of 10-3 mm-1. Maximum error in magnification was 3.2%, and back-projection ray placement was within 0.5 mm. Conclusion: The proposed geometric "self" calibration provides a means for 3D imaging on general noncircular orbits in CBCT systems for which a geometric calibration is either not available or not reproducible. The method forms the basis of advanced "task-based" 3D imaging methods now in development for robotic C-arms.
A Cartesian grid approach with hierarchical refinement for compressible flows
NASA Technical Reports Server (NTRS)
Quirk, James J.
1994-01-01
Many numerical studies of flows that involve complex geometries are limited by the difficulties in generating suitable grids. We present a Cartesian boundary scheme for two-dimensional, compressible flows that is unfettered by the need to generate a computational grid and so it may be used, routinely, even for the most awkward of geometries. In essence, an arbitrary-shaped body is allowed to blank out some region of a background Cartesian mesh and the resultant cut-cells are singled out for special treatment. This is done within a finite-volume framework and so, in principle, any explicit flux-based integration scheme can take advantage of this method for enforcing solid boundary conditions. For best effect, the present Cartesian boundary scheme has been combined with a sophisticated, local mesh refinement scheme, and a number of examples are shown in order to demonstrate the efficacy of the combined algorithm for simulations of shock interaction phenomena.
NASA Astrophysics Data System (ADS)
Hassen, Imen; Gibson, Helen; Hamzaoui-Azaza, Fadoua; Negro, François; Rachid, Khanfir; Bouhlila, Rachida
2016-08-01
The challenge of this study was to create a 3D geological and structural model of the Kasserine Aquifer System (KAS) in central Tunisia and its natural extension into north-east Algeria. This was achieved using an implicit 3D method, which honors prior geological data for both formation boundaries and faults. A current model is presented which provides defendable predictions for the spatial distribution of geology and water resources in aquifers throughout the model-domain. This work has allowed validation of regional scale geology and fault networks in the KAS, and has facilitated the first-ever estimations of groundwater resources in this region by a 3D method. The model enables a preliminary assessment of the hydraulic significance of the major faults by evaluating their influence and role on groundwater flow within and between four compartments of the multi-layered, KAS hydrogeological system. Thus a representative hydrogeological model of the study area is constructed. The possible dual nature of faults in the KAS is discussed in the context that some faults appear to be acting both as barriers to horizontal groundwater flow, and simultaneously as conduits for vertical flow. Also discussed is the possibility that two flow directions occur within the KAS, at a small syncline area of near Feriana. In summary, this work evaluates the influence of aquifer connectivity and the role of faults and geology in groundwater flow within the KAS aquifer system. The current KAS geological model can now be used to guide groundwater managers on the best placement for drilling to test and further refine the understanding of the groundwater system, including the faults connectivity. As more geological data become available, the current model can be easily edited and re-computed to provide an updated model ready for the next stage of investigation by numerical flow modeling.
S2PLOT: Three-dimensional (3D) Plotting Library
NASA Astrophysics Data System (ADS)
Barnes, D. G.; Fluke, C. J.; Bourke, P. D.; Parry, O. T.
2011-03-01
We present a new, three-dimensional (3D) plotting library with advanced features, and support for standard and enhanced display devices. The library - S2PLOT - is written in C and can be used by C, C++ and FORTRAN programs on GNU/Linux and Apple/OSX systems. S2PLOT draws objects in a 3D (x,y,z) Cartesian space and the user interactively controls how this space is rendered at run time. With a PGPLOT inspired interface, S2PLOT provides astronomers with elegant techniques for displaying and exploring 3D data sets directly from their program code, and the potential to use stereoscopic and dome display devices. The S2PLOT architecture supports dynamic geometry and can be used to plot time-evolving data sets, such as might be produced by simulation codes. In this paper, we introduce S2PLOT to the astronomical community, describe its potential applications, and present some example uses of the library.
NASA Astrophysics Data System (ADS)
Jin, BoCheng
2011-12-01
Organic and inorganic fiber reinforced composites with innumerable fiber orientation distributions and fiber geometries are abundantly available in several natural and synthetic structures. Inorganic glass fiber composites have been introduced to numerous applications due to their economical fabrication and tailored structural properties. Numerical characterization of such composite material systems is necessitated due to their intrinsic statistical nature, which renders extensive experimentation prohibitively time consuming and costly. To predict various mechanical behavior and characterizations of Uni-Directional Fiber Composites (UDFC) and Random Fiber Composites (RaFC), we numerically developed Representative Volume Elements (RVE) with high accuracy and efficiency and with complex fiber geometric representations encountered in uni-directional and random fiber networks. In this thesis, the numerical simulations of unidirectional RaFC fiber strand RVE models (VF>70%) are first presented by programming in ABAQUS PYTHON. Secondly, when the cross sectional aspect ratios (AR) of the second phase fiber inclusions are not necessarily one, various types of RVE models with different cross sectional shape fibers are simulated and discussed. A modified random sequential absorption algorithm is applied to enhance the volume fraction number (VF) of the RVE, which the mechanical properties represents the composite material. Thirdly, based on a Spatial Segment Shortest Distance (SSSD) algorithm, a 3-Dimentional RaFC material RVE model is simulated in ABAQUS PYTHON with randomly oriented and distributed straight fibers of high fiber aspect ratio (AR=100:1) and volume fraction (VF=31.8%). Fourthly, the piecewise multi-segments fiber geometry is obtained in MATLAB environment by a modified SSSD algorithm. Finally, numerical methods including the polynomial curve fitting and piecewise quadratic and cubic B-spline interpolation are applied to optimize the RaFC fiber geometries
Noo, Frédéric; Clackdoyle, Rolf; Wagner, Jean-Marc
2002-08-01
This work presents new mathematical results on the inversion of the exponential x-ray transform. It is shown that a reconstruction formula can be obtained for any dataset whose projection directions consist of a union of half great circles on the unit sphere. A basic example of such a dataset is the semi-equatorial band. The discussion in the paper is mostly focused on this example. The reconstruction formula takes the form of a Neumann (geometric) series and is both exact and stable. The exponential x-ray transform has been mainly studied in SPECT imaging. In this context, our results demonstrate mathematically that fully 3D image reconstruction in SPECT with non-zero attenuation does not always require symmetric datasets (opposing views). PMID:12200935
NASA Astrophysics Data System (ADS)
Zhou, Yu; Walker, Richard T.; Elliott, John R.; Parsons, Barry
2016-04-01
Fault dips are usually measured from outcrops in the field or inferred through geodetic or seismological modeling. Here we apply the classic structural geology approach of calculating dip from a fault's 3-D surface trace using recent, high-resolution topography. A test study applied to the 2010 El Mayor-Cucapah earthquake shows very good agreement between our results and those previously determined from field measurements. To obtain a reliable estimate, a fault segment ≥120 m long with a topographic variation ≥15 m is suggested. We then applied this method to the 2013 Balochistan earthquake, getting dips similar to previous estimates. Our dip estimates show a switch from north to south dipping at the southern end of the main trace, which appears to be a response to local extension within a stepover. We suggest that this previously unidentified geometrical complexity may act as the endpoint of earthquake ruptures for the southern end of the Hoshab fault.
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)
Mack, Tobias; Cierpka, Christian; Kähler, Christian J.
2012-11-01
Astigmatism-PTV is a method that allows to measure with a single camera the fully three-dimensional, three-component velocity field. The technique is ideally suited for microfluidic velocity measurements without errors due to in-plane and out-of-plane averaging (Cierpka et al. Meas Scie Tech 21, 2010). Recently it was shown, that the interface between two fluids or the surrounding fluid and droplets or bubbles can be estimated as well with the technique (Rossi et al., Meas Scie Tech 22, 2010). In this contribution the advantages of both techniques are combined to measure the shape of a droplet inside a micro channel along with the internal 3D flow field of the droplet induced by the surrounding fluid. For the current investigation, particles were only distributed within oil-droplets. Therefore the shape of the droplet could be later reconstructed by the volumetric particle positions and the velocity can be estimated tracking the same particles in consecutive frames of the same dataset. The procedure allows the simultaneous determination of the shape and the droplet velocity as well as the inner flow field and offers a great potential for current research.
NASA Astrophysics Data System (ADS)
Kulikov, G. M.; Mamontov, A. A.; Plotnikova, S. V.; Mamontov, S. A.
2015-11-01
A hybrid-mixed ANS four-node shell element by using the sampling surfaces (SaS) technique is developed. The SaS formulation is based on choosing inside the nth layer In not equally spaced SaS parallel to the middle surface of the shell in order to introduce the displacements of these surfaces as basic shell variables. Such choice of unknowns with the consequent use of Lagrange polynomials of degree In - 1 in the thickness direction for each layer permits the presentation of the layered shell formulation in a very compact form. The SaS are located inside each layer at Chebyshev polynomial nodes that allows one to minimize uniformly the error due to the Lagrange interpolation. To implement the efficient analytical integration throughout the element, the enhanced ANS method is employed. The proposed hybrid-mixed four-node shell element is based on the Hu-Washizu variational equation and exhibits a superior performance in the case of coarse meshes. It could be useful for the 3D stress analysis of thick and thin doubly-curved shells since the SaS formulation gives the possibility to obtain numerical solutions with a prescribed accuracy, which asymptotically approach the exact solutions of elasticity as the number of SaS tends to infinity.
NASA Astrophysics Data System (ADS)
Frener, Gernot; Thum, Katharina; Hirz, Mario; Harrich, Alexander
2012-06-01
State of the art automotive development processes are based on virtual product models, which include a digital representation of complete vehicle geometry and structures. Increasing computation performance and continuously growing demands on virtual development processes lead to the application of precise product representation within common CAD software packages. A specific challenge represents the creation of PVC-seams, which are used for corrosion protection of sheet metal components in automotive body in white design. Besides the high requirements of accurate geometric modeling in digital representation, modern development processes call for an increasing level of design automation. To fit both, the present approach introduces a method for an automatic generation of PVC-seams using the functionalities of commercial CAD software. [Figure not available: see fulltext.
Dumane, V; Knoll, M; Green, S; Bakst, R; Hunt, M; Steinberger, E
2014-06-01
Purpose: To predict the dosimetric gain of VMAT over 3D for the treatment ofchestwall/IMN/supraclavicular nodes using geometric parameters acquired during simulation Methods: CT scans for 20 left and 20 right sided patients were retrospectively analyzed toobtain percent ipsilateral lung volume included in the PWT and supraclavicular fields, central lung depth (CLD), maximum lung depth (MLD), separation, chestwall concavity (defined here as the product of CLD and separation) and the maximum heart depth (MHD). VMAT, PWT and P/E plans were done for each case. The ipsilateral lung V20 Gy and mean, total lung V20 Gy and mean, heart V25 Gy and mean were noted for each plan. Correlation coefficients were obtained and linear regression models were built using data from the above training set of patients and then tested on 4 new patients. Results: The decrease in ipsilateral lung V20 Gy, total lung V20 Gy, ipsilateral lung mean and total lung mean with VMAT over PWT significantly (p<0.05) correlated with the percent volume of ipsilateral lung included in the PWT and supraclavicular fields with correlation coefficient values of r = 0.83, r = 0.77, r = 0.78 and r = 0.75 respectively. Significant correlations were also found between MHD and the decrease in heart V25 Gy and mean of r = 0.77 and r = 0.67 respectively. Dosimetric improvement with VMAT over P/E plans showed no correlation to any of the geometric parameters investigated in this study. The dosimetric gain predicted for the 4 test cases by the linear regression models given their respective percent ipsilateral lung volumes fell within the 95% confidence intervals around the best regression fit. Conclusion: The percent ipsilateral lung volume appears to be a strong predictor of the dosimetric gain on using VMAT over PWT apriori.
Price, Sedona H.; Kriek, Mariska; Brammer, Gabriel B.; Conroy, Charlie; Schreiber, Natascha M. Förster; Wuyts, Stijn; Franx, Marijn; Fumagalli, Mattia; Lundgren, Britt; Momcheva, Ivelina; Nelson, Erica J.; Van Dokkum, Pieter G.; Skelton, Rosalind E.; Whitaker, Katherine E.
2014-06-10
The nature of dust in distant galaxies is not well understood, and until recently few direct dust measurements have been possible. We investigate dust in distant star-forming galaxies using near-infrared grism spectra of the 3D-HST survey combined with archival multi-wavelength photometry. These data allow us to make a direct comparison between dust around star-forming regions (A {sub V,} {sub H} {sub II}) and the integrated dust content (A {sub V,} {sub star}). We select a sample of 163 galaxies between 1.36 ≤ z ≤ 1.5 with Hα signal-to-noise ratio ≥5 and measure Balmer decrements from stacked spectra to calculate A {sub V,} {sub H} {sub II}. First, we stack spectra in bins of A {sub V,} {sub star}, and find that A {sub V,} {sub H} {sub II} = 1.86 A {sub V,} {sub star}, with a significance of σ = 1.7. Our result is consistent with the two-component dust model, in which galaxies contain both diffuse and stellar birth cloud dust. Next, we stack spectra in bins of specific star formation rate (log SSFR), star formation rate (log SFR), and stellar mass (log M {sub *}). We find that on average A {sub V,} {sub H} {sub II} increases with SFR and mass, but decreases with increasing SSFR. Interestingly, the data hint that the amount of extra attenuation decreases with increasing SSFR. This trend is expected from the two-component model, as the extra attenuation will increase once older stars outside the star-forming regions become more dominant in the galaxy spectrum. Finally, using Balmer decrements we derive dust-corrected Hα SFRs, and find that stellar population modeling produces incorrect SFRs if rapidly declining star formation histories are included in the explored parameter space.
NASA Technical Reports Server (NTRS)
Price, Sedona H.; Kriek, Mariska; Brammer, Gabriel B; Conroy, Charlie; Schreiber, Natascha M. Foerster; Franx, Marijn; Fumagalli, Mattia; Lundren, Britt; Momcheva, Ivelina; Nelson, Erica J.; Rix, Hans-Walter; Skelton, Rosalind E.; VanDokkum, Pieter G.; Tease, Katherine Whitaker; Wuyts, Stijn
2013-01-01
The nature of dust in distant galaxies is not well understood, and until recently few direct dust measurements have been possible. We investigate dust in distant star-forming galaxies using near-infrared grism spectra of the 3D-HST survey combined with archival multi-wavelength photometry. These data allow us to make a direct comparison between dust towards star-forming regions (measured using Balmer decrements) and the integrated dust properties (derived by comparing spectral energy distributions [SEDs] with stellar population and dust models) for a statistically significant sample of distant galaxies. We select a sample of 163 galaxies between 1.36< or = z< or = 1.5 with H(alpha) SNR > or = 5 and measure Balmer decrements from stacked spectra. First, we stack spectra in bins of integrated stellar dust attenuation, and find that there is extra dust extinction towards star-forming regions (AV,HII is 1.81 times the integrated AV, star), though slightly lower than found for low-redshift starburst galaxies. Next, we stack spectra in bins of specific star formation rate (log sSFR), star formation rate (log SFR), and stellar mass (logM*). We find that on average AV,HII increases with SFR and mass, but decreases with increasing sSFR. The amount of extra extinction also decreases with increasing sSFR and decreasing stellar mass. Our results are consistent with the two-phase dust model - in which galaxies contain both a diffuse and a stellar birth cloud dust component - as the extra extinction will increase once older stars outside the star-forming regions become more dominant. Finally, using our Balmer decrements we derive dust-corrected H(alpha) SFRs, and find evidence that SED fitting produces incorrect SFRs if very rapidly declining SFHs are included in the explored parameter space. Subject headings: dust, extinction- galaxies: evolution- galaxies: high-redshift
NASA Astrophysics Data System (ADS)
McEvilly, A.; Abimbola, A.; Chan, J. H.; Strayer, L. M.
2015-12-01
California State University, East Bay (CSUEB), located in Hayward, California, lies atop the San Leandro block (SLB) in the Hayward fault zone. The SLB is a J-K aged lithotectonic assemblage dominated by gabbro and intercalated with minor volcanics and sediments. It is bound by the subparallel northwest-trending western Hayward and eastern Chabot (CF) faults and pervasively cut by anastomosing secondary faults. The block itself is ~30 km along strike and 2-3 km wide. Previous studies suggest the block dips steeply to the northeast and extends to a depth of at least 7 km. In May of 2015, as part of an ongoing collaborative effort led by the USGS to create a 3D velocity model of the San Francisco Bay Area, researchers from CSUEB and the USGS conducted a seismic survey on the CSUEB campus. The primary goal of this pilot study was to locate the trace of the CF on the CSUEB campus and to determine bedrock depth. We deployed a 60-channel, 300m profile using 4.5Hz sensors spaced at 5m intervals. Active seismic sources were used at each geophone location. A 226kg accelerated weight-drop was used to generate P and Rayleigh waves for P-wave tomography and multichannel analysis of surface waves (MASW), and a 3.5kg sledgehammer and block were used to generate S and Love waves for S-wave tomography and multichannel analysis of Love waves (MALW). Preliminary P-wave tomography, MASW, and MALW results from this pilot study suggest the location of an eastward-dipping CF as well as the presence of a high-velocity unit at about 20m depth, presumably an unmapped sliver of bedrock from the San Leandro block. Further studies planned for the fall of 2015 include additional seismic lines and surface mapping along the Chabot fault on and near the CSUEB campus. These new geophysical, GPS, and field geological data will be integrated with LiDAR imagery and existing geological, gravity and magnetic maps to create a 3-dimensional model of the portion of the SLB that contains the CSUEB campus.
NASA Astrophysics Data System (ADS)
Liu, Y.; Luhmann, J.; Wang, L.; Lynch, B.; Huttunen, E.; Lin, R.; Bale, S.; Russell, C.; Galvin, T.
2008-05-01
Coronal mass ejections (CMEs) are often assumed to be magnetic flux ropes, but direct proof has been lacking. A key feature, resulting from the translational symmetry of a flux rope, is that the total transverse pressure as well as the axial magnetic field has the same functional form over the vector potential along any crossing of the flux rope. We test this feature (and hence the flux-rope structure) by reconstructing the 22 May 2007 magnetic cloud (MC) observed at STEREO B, Wind/ACE and possibly STEREO A with the Grad-Shafranov (GS) method. The model output from reconstruction at STEREO B agrees fairly well with the magnetic field and thermal pressure observed at ACE/Wind; the separation between STEREO B and ACE/Wind is about 0.06 AU, almost half of the MC radial width. For the first time, we reproduce observations at one spacecraft with data from another well- separated spacecraft, which provides compelling evidence for the flux-rope geometry and is of importance for understanding CME initiation and propagation. The magnetic field line length, calculated from the velocity dispersion of energetic electrons within the MC, is substantially longer than the nominal Parker spiral field, which is consistent with a flux-rope structure. The reconstruction gives a global configuration of the MC at different spacecraft with a flattened cross section owing to the solar wind radial expansion. A similar approach is also performed with a force-free flux-rope fitting model in an effort to compare with the GS reconstruction results.
NASA Astrophysics Data System (ADS)
Ando, R.
2014-12-01
The boundary integral equation method formulated in the real space and time domain (BIEM-ST) has been used as a powerful tool to analyze the earthquake rupture dynamics on non-planar faults. Generally, BIEM is more accurate than volumetric methods such as the finite difference method and the finite difference method. With the recent development of the high performance computing environment, the earthquake rupture simulation studies have been conducted considering three dimensional realistic fault geometry models. However, the utility of BIEM-ST has been limited due to its heavy computational demanding increased depending on square of time steps (N2), which was needed to evaluate the historic integration. While BIEM can be efficient with the spectral domain formulation, the applications of such a method are limited to planar fault cases. In this study, we propose a new method to reduce the calculation time of BIEM-ST to linear of time step (N) without degrading the accuracy in the 3 dimensional modeling space. We extends the method proposed earlier for the case of the 2 dimensional framework, applying the asymptotic expressions of the elasto-dynamic Green's functions. This method uses the physical nature of the stress Green's function as dividing the causality cone according to the distances from the wave-fronts. The scalability of this method is shown on the parallel computing environment of the distributed memory. We demonstrate the applicability to analyses of subduction earthquake cases, suffering long time from the numerical limitations of previously available BIEMs. We analyze the dynamic rupture processes on dipping reverse faults embed in a three dimensional elastic half space.
ERIC Educational Resources Information Center
Planinsic, G.; Kos, M.; Jerman, R.
2004-01-01
It is quite easy to make a version of the well known Cartesian diver experiment that uses two immiscible liquids. This allows students to test their knowledge of density and pressure in explaining the diver's behaviour. Construction details are presented here together with a mathematical model to explain the observations.
On the Use of CAD and Cartesian Methods for Aerodynamic Optimization
NASA Technical Reports Server (NTRS)
Nemec, M.; Aftosmis, M. J.; Pulliam, T. H.
2004-01-01
The objective for this paper is to present the development of an optimization capability for Curt3D, a Cartesian inviscid-flow analysis package. We present the construction of a new optimization framework and we focus on the following issues: 1) Component-based geometry parameterization approach using parametric-CAD models and CAPRI. A novel geometry server is introduced that addresses the issue of parallel efficiency while only sparingly consuming CAD resources; 2) The use of genetic and gradient-based algorithms for three-dimensional aerodynamic design problems. The influence of noise on the optimization methods is studied. Our goal is to create a responsive and automated framework that efficiently identifies design modifications that result in substantial performance improvements. In addition, we examine the architectural issues associated with the deployment of a CAD-based approach in a heterogeneous parallel computing environment that contains both CAD workstations and dedicated compute engines. We demonstrate the effectiveness of the framework for a design problem that features topology changes and complex geometry.
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.
Cartesian control of redundant robots
NASA Technical Reports Server (NTRS)
Colbaugh, R.; Glass, K.
1989-01-01
A Cartesian-space position/force controller is presented for redundant robots. The proposed control structure partitions the control problem into a nonredundant position/force trajectory tracking problem and a redundant mapping problem between Cartesian control input F is a set member of the set R(sup m) and robot actuator torque T is a set member of the set R(sup n) (for redundant robots, m is less than n). The underdetermined nature of the F yields T map is exploited so that the robot redundancy is utilized to improve the dynamic response of the robot. This dynamically optimal F yields T map is implemented locally (in time) so that it is computationally efficient for on-line control; however, it is shown that the map possesses globally optimal characteristics. Additionally, it is demonstrated that the dynamically optimal F yields T map can be modified so that the robot redundancy is used to simultaneously improve the dynamic response and realize any specified kinematic performance objective (e.g., manipulability maximization or obstacle avoidance). Computer simulation results are given for a four degree of freedom planar redundant robot under Cartesian control, and demonstrate that position/force trajectory tracking and effective redundancy utilization can be achieved simultaneously with the proposed controller.
NASA Astrophysics Data System (ADS)
Crameri, Fabio; Tackley, Paul J.
2014-07-01
We present temporally evolving 3-D global mantle convection models with single-sided subduction and a free surface in both 3-D Cartesian and fully spherical geometry. Special focus is given to the spontaneous development of three-dimensional structures at the surface and in the upper mantle. We find that an arcuate shape is the natural form for trenches and slabs. Cartesian models are used first to study the dynamic evolution of subduction zones, spreading ridges, and interconnected transform features. These experiments highlight the strong variation of spontaneously developing, arcuate slab curvature and subduction polarity along the trench strike. The spontaneous development of spreading ridges leads to lateral offsets between separated segments that are characterized by normal transform motion. Spherical models then allow insights into the evolution of plate tectonics on a sphere. Investigated are the spontaneous evolution of slab geometry, trench motion, and subduction-induced mantle flow. Two new dynamical features are discovered: "back-slab spiral flow" and "slab tunneling." 2014. American Geophysical Union. All Rights Reserved.
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
A Cartesian cut cell method for rarefied flow simulations around moving obstacles
NASA Astrophysics Data System (ADS)
Dechristé, G.; Mieussens, L.
2016-06-01
For accurate simulations of rarefied gas flows around moving obstacles, we propose a cut cell method on Cartesian grids: it allows exact conservation and accurate treatment of boundary conditions. Our approach is designed to treat Cartesian cells and various kinds of cut cells by the same algorithm, with no need to identify the specific shape of each cut cell. This makes the implementation quite simple, and allows a direct extension to 3D problems. Such simulations are also made possible by using an adaptive mesh refinement technique and a hybrid parallel implementation. This is illustrated by several test cases, including a 3D unsteady simulation of the Crookes radiometer.
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.
SU-E-I-41: Non-Cartesian MR Image Reconstruction with Integrated Gradient Non-Linearity Correction
Tao, S; Trzasko, JD; Polley, TW; Shu, Y; Bernstein, MA
2014-06-01
Purpose: Nonlinearities in the spatial encoding gradients of MRI systems cause geometric distortion in images. Typically, this is retrospectively corrected via image-domain interpolation (a.k.a., “gradwarp”) albeit with a loss of spatial resolution. For non-Cartesian MRI, the latter problem is exaggerated by noise and undersampling artifact. In this study, we describe a novel correction strategy that accounts for gradient nonlinearities during — rather than after — non-Cartesian MRI reconstruction, and demonstrate that this approach mitigates the resolution loss that can occur with standard methods. Methods: To test the proposed method, the American College of Radiology (ACR) quality control phantom was scanned on at 1.5 T (General Electric, v16.0, “zoom” gradient) using a 1.6x undersampled 3D non- Cartesian Shells trajectory (GRE, FOV=24 cm3, 120 shells, 16552 shots, 512 readout, matrix=2403). Image reconstruction was first performed via standard k-space density-compensated gridding and retrospectively corrected via cubic spline interpolation. Image reconstruction was then separately performed using a k-space and image-domain densitycompensated type-3 non-uniform fast Fourier transform (NUFFT), which provides a direct mapping between non-Cartesian k-space samples and warped image space voxel locations. Thus, no separate distortion correction procedure is needed for the proposed approach. The gradient distortion field was determined using vendor provided calibration data. Results: Phantom scan results show that both processing approaches successfully correct geometric distortion. However, visual inspection of the ACR phantom spatial resolution inserts shows that the proposed strategy preserves the resolution of the nominal (uncorrected) reconstruction while “gradwarp” imparts marked spatial blurring (especially for the 1.0 and 1.1 mm inserts) and thus resolution loss. Conclusion: We've presented a novel reconstruction strategy for non-Cartesian MRI
A Cartesian embedded boundary method for hyperbolic conservation laws
Sjogreen, B; Petersson, N A
2006-12-04
The authors develop an embedded boundary finite difference technique for solving the compressible two- or three-dimensional Euler equations in complex geometries on a Cartesian grid. The method is second order accurate with an explicit time step determined by the grid size away from the boundary. Slope limiters are used on the embedded boundary to avoid non-physical oscillations near shock waves. They show computed examples of supersonic flow past a cylinder and compare with results computed on a body fitted grid. Furthermore, they discuss the implementation of the method for thin geometries, and show computed examples of transonic flow past an airfoil.
NASA Astrophysics Data System (ADS)
Pereira, J. P.; Duarte, C. A.; Jiao, X.; Guoy, D.
2009-06-01
This paper presents a study of generalized enrichment functions for 3D curved crack fronts. Two coordinate systems used in the definition of singular curved crack front enrichment functions are analyzed. In the first one, a set of Cartesian coordinate systems defined along the crack front is used. In the second case, the geometry of the crack front is approximated by a set of curvilinear coordinate systems. A description of the computation of derivatives of enrichment functions and curvilinear base vectors is presented. The coordinate systems are automatically defined using geometrical information provided by an explicit representation of the crack surface. A detailed procedure to accurately evaluate the surface normal, conormal and tangent vectors along curvilinear crack fronts in explicit crack surface representations is also presented. An accurate and robust definition of orthonormal vectors along crack fronts is crucial for the proper definition of enrichment functions. Numerical experiments illustrate the accuracy and robustness of the proposed approaches.
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.
3-D Spherical Mantle Convection with Radial Basis Functions
NASA Astrophysics Data System (ADS)
Flyer, N.; Wright, G. B.; Yuen, D.
2008-12-01
In the past 25 years a wide variety of numerical methods, such as finite-difference, finite-volume , finite- elements, and pseudospectral methods have been employed to study the problem of 3-D mantle convection. All have specialized strengths but also serious weaknesses. The first three methods are generally considered low-order and can involve high algorithmic complexity (as in triangular elements). Spectrally accurate methods do not practically allow for local mesh refinement and often involve cumbersome algebra. Here, we introduce a new grid/mesh-free approach using radial basis functions (RBFs). It has the advantage of being spectrally accurate for arbitrary node layouts in multi-dimensions with extreme algorithmic simplicity, and naturally permits local node refinement. It has been shown for shallow-water equations and vortex flows that RBFs outperform other numerical methods in the sense that they obtain a much higher accuracy for the same spatial resolution while being able to take unusually large time steps. One virtue of the RBF scheme is the ability to use a simple Cartesian geometry while implementing the required boundary conditions for the temperature, velocity and stresses on a spherical surface of both the outer( planetary surface ) and inner shell ( core-mantle boundary ). The velocity and stress components are expressed in terms of the scalar potential approach (Zebib and Schubert, 1982) and the other remaining variable is the perturbed temperature field. We have studied the problem from the onset of convection to a modest nonlinear regime.
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?
ATHENA 3D: A finite element code for ultrasonic wave propagation
NASA Astrophysics Data System (ADS)
Rose, C.; Rupin, F.; Fouquet, T.; Chassignole, B.
2014-04-01
The understanding of wave propagation phenomena requires use of robust numerical models. 3D finite element (FE) models are generally prohibitively time consuming. However, advances in computing processor speed and memory allow them to be more and more competitive. In this context, EDF R&D developed the 3D version of the well-validated FE code ATHENA2D. The code is dedicated to the simulation of wave propagation in all kinds of elastic media and in particular, heterogeneous and anisotropic materials like welds. It is based on solving elastodynamic equations in the calculation zone expressed in terms of stress and particle velocities. The particularity of the code relies on the fact that the discretization of the calculation domain uses a Cartesian regular 3D mesh while the defect of complex geometry can be described using a separate (2D) mesh using the fictitious domains method. This allows combining the rapidity of regular meshes computation with the capability of modelling arbitrary shaped defects. Furthermore, the calculation domain is discretized with a quasi-explicit time evolution scheme. Thereby only local linear systems of small size have to be solved. The final step to reduce the computation time relies on the fact that ATHENA3D has been parallelized and adapted to the use of HPC resources. In this paper, the validation of the 3D FE model is discussed. A cross-validation of ATHENA 3D and CIVA is proposed for several inspection configurations. The performances in terms of calculation time are also presented in the cases of both local computer and computation cluster use.
NASA Astrophysics Data System (ADS)
Ge, Wenjun; Modest, Michael F.; Roy, Somesh P.
2016-03-01
The high-order spherical harmonics (PN) method for 2-D Cartesian domains is extracted from the 3-D formulation. The number of equations and intensity coefficients reduces to (N + 1)2 / 4 in the 2-D Cartesian formulation compared with N(N + 1) / 2 for the general 3-D PN formulation. The Marshak boundary conditions are extended to solve problems with nonblack and mixed diffuse-specular surfaces. Additional boundary conditions for specified radiative wall flux, for symmetry/specular reflection boundaries have also been developed. The mathematical details of the formulations and their implementation in the OpenFOAM finite volume based CFD software platform are presented. The accuracy and computational cost of the 2-D Cartesian PN are compared with that of the 3-D PN solver and a Photon Monte Carlo solver for a square enclosure, as well as a 45° wedge geometry with variable radiative properties. The new boundary conditions have been applied for both test cases, and the boundary condition for mixed diffuse-specular surfaces is further illustrated by numerical examples of a rectangular geometry enclosed by walls with different surface characteristics.
GPU accelerated simulations of 3D deterministic particle transport using discrete ordinates method
NASA Astrophysics Data System (ADS)
Gong, Chunye; Liu, Jie; Chi, Lihua; Huang, Haowei; Fang, Jingyue; Gong, Zhenghu
2011-07-01
Graphics Processing Unit (GPU), originally developed for real-time, high-definition 3D graphics in computer games, now provides great faculty in solving scientific applications. The basis of particle transport simulation is the time-dependent, multi-group, inhomogeneous Boltzmann transport equation. The numerical solution to the Boltzmann equation involves the discrete ordinates ( Sn) method and the procedure of source iteration. In this paper, we present a GPU accelerated simulation of one energy group time-independent deterministic discrete ordinates particle transport in 3D Cartesian geometry (Sweep3D). The performance of the GPU simulations are reported with the simulations of vacuum boundary condition. The discussion of the relative advantages and disadvantages of the GPU implementation, the simulation on multi GPUs, the programming effort and code portability are also reported. The results show that the overall performance speedup of one NVIDIA Tesla M2050 GPU ranges from 2.56 compared with one Intel Xeon X5670 chip to 8.14 compared with one Intel Core Q6600 chip for no flux fixup. The simulation with flux fixup on one M2050 is 1.23 times faster than on one X5670.
GPU accelerated simulations of 3D deterministic particle transport using discrete ordinates method
Gong Chunye; Liu Jie; Chi Lihua; Huang Haowei; Fang Jingyue; Gong Zhenghu
2011-07-01
Graphics Processing Unit (GPU), originally developed for real-time, high-definition 3D graphics in computer games, now provides great faculty in solving scientific applications. The basis of particle transport simulation is the time-dependent, multi-group, inhomogeneous Boltzmann transport equation. The numerical solution to the Boltzmann equation involves the discrete ordinates (S{sub n}) method and the procedure of source iteration. In this paper, we present a GPU accelerated simulation of one energy group time-independent deterministic discrete ordinates particle transport in 3D Cartesian geometry (Sweep3D). The performance of the GPU simulations are reported with the simulations of vacuum boundary condition. The discussion of the relative advantages and disadvantages of the GPU implementation, the simulation on multi GPUs, the programming effort and code portability are also reported. The results show that the overall performance speedup of one NVIDIA Tesla M2050 GPU ranges from 2.56 compared with one Intel Xeon X5670 chip to 8.14 compared with one Intel Core Q6600 chip for no flux fixup. The simulation with flux fixup on one M2050 is 1.23 times faster than on one X5670.
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)
3D whiteboard: collaborative sketching with 3D-tracked smart phones
NASA Astrophysics Data System (ADS)
Lue, James; Schulze, Jürgen P.
2014-02-01
We present the results of our investigation of the feasibility of a new approach for collaborative drawing in 3D, based on Android smart phones. Our approach utilizes a number of fiduciary markers, placed in the working area where they can be seen by the smart phones' cameras, in order to estimate the pose of each phone in the room. Our prototype allows two users to draw 3D objects with their smart phones by moving their phones around in 3D space. For example, 3D lines are drawn by recording the path of the phone as it is moved around in 3D space, drawing line segments on the screen along the way. Each user can see the virtual drawing space on their smart phones' displays, as if the display was a window into this space. Besides lines, our prototype application also supports 3D geometry creation, geometry transformation operations, and it shows the location of the other user's phone.
A Deformable Generic 3D Model of Haptoral Anchor of Monogenean
Teo, Bee Guan; Dhillon, Sarinder Kaur; Lim, Lee Hong Susan
2013-01-01
In this paper, a digital 3D model which allows for visualisation in three dimensions and interactive manipulation is explored as a tool to help us understand the structural morphology and elucidate the functions of morphological structures of fragile microorganisms which defy live studies. We developed a deformable generic 3D model of haptoral anchor of dactylogyridean monogeneans that can subsequently be deformed into different desired anchor shapes by using direct manipulation deformation technique. We used point primitives to construct the rectangular building blocks to develop our deformable 3D model. Point primitives are manually marked on a 2D illustration of an anchor on a Cartesian graph paper and a set of Cartesian coordinates for each point primitive is manually extracted from the graph paper. A Python script is then written in Blender to construct 3D rectangular building blocks based on the Cartesian coordinates. The rectangular building blocks are stacked on top or by the side of each other following their respective Cartesian coordinates of point primitive. More point primitives are added at the sites in the 3D model where more structural variations are likely to occur, in order to generate complex anchor structures. We used Catmull-Clark subdivision surface modifier to smoothen the surface and edge of the generic 3D model to obtain a smoother and more natural 3D shape and antialiasing option to reduce the jagged edges of the 3D model. This deformable generic 3D model can be deformed into different desired 3D anchor shapes through direct manipulation deformation technique by aligning the vertices (pilot points) of the newly developed deformable generic 3D model onto the 2D illustrations of the desired shapes and moving the vertices until the desire 3D shapes are formed. In this generic 3D model all the vertices present are deployed for displacement during deformation. PMID:24204903
A deformable generic 3D model of haptoral anchor of Monogenean.
Teo, Bee Guan; Dhillon, Sarinder Kaur; Lim, Lee Hong Susan
2013-01-01
In this paper, a digital 3D model which allows for visualisation in three dimensions and interactive manipulation is explored as a tool to help us understand the structural morphology and elucidate the functions of morphological structures of fragile microorganisms which defy live studies. We developed a deformable generic 3D model of haptoral anchor of dactylogyridean monogeneans that can subsequently be deformed into different desired anchor shapes by using direct manipulation deformation technique. We used point primitives to construct the rectangular building blocks to develop our deformable 3D model. Point primitives are manually marked on a 2D illustration of an anchor on a Cartesian graph paper and a set of Cartesian coordinates for each point primitive is manually extracted from the graph paper. A Python script is then written in Blender to construct 3D rectangular building blocks based on the Cartesian coordinates. The rectangular building blocks are stacked on top or by the side of each other following their respective Cartesian coordinates of point primitive. More point primitives are added at the sites in the 3D model where more structural variations are likely to occur, in order to generate complex anchor structures. We used Catmull-Clark subdivision surface modifier to smoothen the surface and edge of the generic 3D model to obtain a smoother and more natural 3D shape and antialiasing option to reduce the jagged edges of the 3D model. This deformable generic 3D model can be deformed into different desired 3D anchor shapes through direct manipulation deformation technique by aligning the vertices (pilot points) of the newly developed deformable generic 3D model onto the 2D illustrations of the desired shapes and moving the vertices until the desire 3D shapes are formed. In this generic 3D model all the vertices present are deployed for displacement during deformation. PMID:24204903
A Cartesian Adaptive Level Set Method for Two-Phase Flows
NASA Technical Reports Server (NTRS)
Ham, F.; Young, Y.-N.
2003-01-01
In the present contribution we develop a level set method based on local anisotropic Cartesian adaptation as described in Ham et al. (2002). Such an approach should allow for the smallest possible Cartesian grid capable of resolving a given flow. The remainder of the paper is organized as follows. In section 2 the level set formulation for free surface calculations is presented and its strengths and weaknesses relative to the other free surface methods reviewed. In section 3 the collocated numerical method is described. In section 4 the method is validated by solving the 2D and 3D drop oscilation problem. In section 5 we present some results from more complex cases including the 3D drop breakup in an impulsively accelerated free stream, and the 3D immiscible Rayleigh-Taylor instability. Conclusions are given in section 6.
NASA Technical Reports Server (NTRS)
Aftosmis, M. J.; Berger, M. J.; Adomavicius, G.; Nixon, David (Technical Monitor)
1998-01-01
The work presents a new method for on-the-fly domain decomposition technique for mapping grids and solution algorithms to parallel machines, and is applicable to both shared-memory and message-passing architectures. It will be demonstrated on the Cray T3E, HP Exemplar, and SGI Origin 2000. Computing time has been secured on all these platforms. The decomposition technique is an outgrowth of techniques used in computational physics for simulations of N-body problems and the event horizons of black holes, and has not been previously used by the CFD community. Since the technique offers on-the-fly partitioning, it offers a substantial increase in flexibility for computing in heterogeneous environments, where the number of available processors may not be known at the time of job submission. In addition, since it is dynamic it permits the job to be repartitioned without global communication in cases where additional processors become available after the simulation has begun, or in cases where dynamic mesh adaptation changes the mesh size during the course of a simulation. The platform for this partitioning strategy is a completely new Cartesian Euler solver tarcreted at parallel machines which may be used in conjunction with Ames' "Cart3D" arbitrary geometry simulation package.
Extending a CAD-Based Cartesian Mesh Generator for the Lattice Boltzmann Method
Cantrell, J Nathan; Inclan, Eric J; Joshi, Abhijit S; Popov, Emilian L; Jain, Prashant K
2012-01-01
This paper describes the development of a custom preprocessor for the PaRAllel Thermal Hydraulics simulations using Advanced Mesoscopic methods (PRATHAM) code based on an open-source mesh generator, CartGen [1]. PRATHAM is a three-dimensional (3D) lattice Boltzmann method (LBM) based parallel flow simulation software currently under development at the Oak Ridge National Laboratory. The LBM algorithm in PRATHAM requires a uniform, coordinate system-aligned, non-body-fitted structured mesh for its computational domain. CartGen [1], which is a GNU-licensed open source code, already comes with some of the above needed functionalities. However, it needs to be further extended to fully support the LBM specific preprocessing requirements. Therefore, CartGen is being modified to (i) be compiler independent while converting a neutral-format STL (Stereolithography) CAD geometry to a uniform structured Cartesian mesh, (ii) provide a mechanism for PRATHAM to import the mesh and identify the fluid/solid domains, and (iii) provide a mechanism to visually identify and tag the domain boundaries on which to apply different boundary conditions.
NASA Astrophysics Data System (ADS)
Ronchin, Erika; Masterlark, Timothy; Dawson, John; Saunders, Steve; Martí Molist, Joan
2015-04-01
In this study, we present a method to fully integrate a family of finite element models (FEMs) into the regularized linear inversion of InSAR data collected at Rabaul caldera (PNG) between February 2007 and December 2010. During this period the caldera experienced a long-term steady subsidence that characterized surface movement both inside the caldera and outside, on its western side. The inversion is based on an array of FEM sources in the sense that the Green's function matrix is a library of forward numerical displacement solutions generated by the sources of an array common to all FEMs. Each entry of the library is the LOS surface displacement generated by injecting a unity mass of fluid, of known density and bulk modulus, into a different source cavity of the array for each FEM. By using FEMs, we are taking advantage of their capability of including topography and heterogeneous distribution of elastic material properties. All FEMs of the family share the same mesh in which only one source is activated at the time by removing the corresponding elements and applying the unity fluid flux. The domain therefore only needs to be discretized once. This precludes remeshing for each activated source, thus reducing computational requirements, often a downside of FEM-based inversions. Without imposing an a-priori source, the method allows us to identify, from a least-squares standpoint, a complex distribution of fluid flux (or change in pressure) with a 3D free geometry within the source array, as dictated by the data. The results of applying the proposed inversion to Rabaul InSAR data show a shallow magmatic system under the caldera made of two interconnected lobes located at the two opposite sides of the caldera. These lobes could be consistent with feeding reservoirs of the ongoing Tavuvur volcano eruption of andesitic products, on the eastern side, and of the past Vulcan volcano eruptions of more evolved materials, on the western side. The interconnection and
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.
JAR3D Webserver: Scoring and aligning RNA loop sequences to known 3D motifs.
Roll, James; Zirbel, Craig L; Sweeney, Blake; Petrov, Anton I; Leontis, Neocles
2016-07-01
Many non-coding RNAs have been identified and may function by forming 2D and 3D structures. RNA hairpin and internal loops are often represented as unstructured on secondary structure diagrams, but RNA 3D structures show that most such loops are structured by non-Watson-Crick basepairs and base stacking. Moreover, different RNA sequences can form the same RNA 3D motif. JAR3D finds possible 3D geometries for hairpin and internal loops by matching loop sequences to motif groups from the RNA 3D Motif Atlas, by exact sequence match when possible, and by probabilistic scoring and edit distance for novel sequences. The scoring gauges the ability of the sequences to form the same pattern of interactions observed in 3D structures of the motif. The JAR3D webserver at http://rna.bgsu.edu/jar3d/ takes one or many sequences of a single loop as input, or else one or many sequences of longer RNAs with multiple loops. Each sequence is scored against all current motif groups. The output shows the ten best-matching motif groups. Users can align input sequences to each of the motif groups found by JAR3D. JAR3D will be updated with every release of the RNA 3D Motif Atlas, and so its performance is expected to improve over time. PMID:27235417
An Improved Version of TOPAZ 3D
Krasnykh, Anatoly
2003-07-29
An improved version of the TOPAZ 3D gun code is presented as a powerful tool for beam optics simulation. In contrast to the previous version of TOPAZ 3D, the geometry of the device under test is introduced into TOPAZ 3D directly from a CAD program, such as Solid Edge or AutoCAD. In order to have this new feature, an interface was developed, using the GiD software package as a meshing code. The article describes this method with two models to illustrate the results.
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.
Coolant side heat transfer with rotation: User manual for 3D-TEACH with rotation
NASA Technical Reports Server (NTRS)
Syed, S. A.; James, R. H.
1989-01-01
This program solves the governing transport equations in Reynolds average form for the flow of a 3-D, steady state, viscous, heat conducting, multiple species, single phase, Newtonian fluid with combustion. The governing partial differential equations are solved in physical variables in either a Cartesian or cylindrical coordinate system. The effects of rotation on the momentum and enthalpy calculations modeled in Cartesian coordinates are examined. The flow of the fluid should be confined and subsonic with a maximum Mach number no larger than 0.5. This manual describes the operating procedures and input details for executing a 3D-TEACH computation.
3D barcodes: theoretical aspects and practical implementation
NASA Astrophysics Data System (ADS)
Gladstein, David; Kakarala, Ramakrishna; Baharav, Zachi
2015-02-01
This paper introduces the concept of three dimensional (3D) barcodes. A 3D barcode is composed of an array of 3D cells, called modules, and each can be either filled or empty, corresponding to two possible values of a bit. These barcodes have great theoretical promise thanks to their very large information capacity, which grows as the cube of the linear size of the barcode, and in addition are becoming practically manufacturable thanks to the ubiquitous use of 3D printers. In order to make these 3D barcodes practical for consumers, it is important to keep the decoding simple using commonly available means like smartphones. We therefore limit ourselves to decoding mechanisms based only on three projections of the barcode, which imply specific constraints on the barcode itself. The three projections produce the marginal sums of the 3D cube, which are the counts of filled-in modules along each Cartesian axis. In this paper we present some of the theoretical aspects of the 2D and 3D cases, and describe the resulting complexity of the 3D case. We then describe a method to reduce these complexities into a practical application. The method features an asymmetric coding scheme, where the decoder is much simpler than the encoder. We close by demonstrating 3D barcodes we created and their usability.
Modeling ICF Spherical Implosion Instabilities in 3D with Exact Energy Conservation
NASA Astrophysics Data System (ADS)
Fatenejad, Milad; Moses, Gregory
2009-11-01
We will present the results of 3D instability simulations performed on spherically convergent geometries with a new 3D Lagrangian hydrodynamics code, cooper. The code uses a compatible discretization of the conservation equations to ensure that energy is conserved to within machine round off error [Caramana JCP 146, 227 (1998)]. Modifications are made to the discrete equations to ensure that spherically symmetric implosions can be performed on non-orthogonal Cartesian grids [Caramana JCP 157, 89 (2000)]. Subzonal restoring forces counteract anomalous grid distortions [Carmana JCP 142, 521 (1998)] and an edge-centered viscosity is used to capture shocks [Caramana JCP 215, 385 (2006)]. Cooper is parallelized using domain decomposition. This is necessary due to the large processor and memory requirements associated with simulations in three dimensions. Advanced computational libraries are used to reduce the complexity of the code without sacrificing features. One example is the MOAB library [Tautges Engr. Comput. 20, 286 (2004)] which manages the mesh and is responsible for communicating information between processes.
Magnetic Properties of 3D Printed Toroids
NASA Astrophysics Data System (ADS)
Bollig, Lindsey; Otto, Austin; Hilpisch, Peter; Mowry, Greg; Nelson-Cheeseman, Brittany; Renewable Energy; Alternatives Lab (REAL) Team
Transformers are ubiquitous in electronics today. Although toroidal geometries perform most efficiently, transformers are traditionally made with rectangular cross-sections due to the lower manufacturing costs. Additive manufacturing techniques (3D printing) can easily achieve toroidal geometries by building up a part through a series of 2D layers. To get strong magnetic properties in a 3D printed transformer, a composite filament is used containing Fe dispersed in a polymer matrix. How the resulting 3D printed toroid responds to a magnetic field depends on two structural factors of the printed 2D layers: fill factor (planar density) and fill pattern. In this work, we investigate how the fill factor and fill pattern affect the magnetic properties of 3D printed toroids. The magnetic properties of the printed toroids are measured by a custom circuit that produces a hysteresis loop for each toroid. Toroids with various fill factors and fill patterns are compared to determine how these two factors can affect the magnetic field the toroid can produce. These 3D printed toroids can be used for numerous applications in order to increase the efficiency of transformers by making it possible for manufacturers to make a toroidal geometry.
3D Magnetron simulation with CST STUDIO SUITE
Balk, Monika C.
2011-07-01
The modeling of magnetrons compared to other tubes is more difficult since it requires 3D modeling rather than a 2D investigation. This is not only due to the geometry which can include complicated details to be modeled in 3D but also due to the interaction process itself. The electric field, magnetic field and particle movement span a 3D space. In this paper 3D simulations of a strapped magnetron with CSTSTUDIO SUITE{sup TM} are presented. (author)
Turing instabilities on Cartesian product networks
Asllani, Malbor; Busiello, Daniel M.; Carletti, Timoteo; Fanelli, Duccio; Planchon, Gwendoline
2015-01-01
The problem of Turing instabilities for a reaction-diffusion system defined on a complex Cartesian product network is considered. To this end we operate in the linear regime and expand the time dependent perturbation on a basis formed by the tensor product of the eigenvectors of the discrete Laplacian operators, associated to each of the individual networks that build the Cartesian product. The dispersion relation which controls the onset of the instability depends on a set of discrete wavelengths, the eigenvalues of the aforementioned Laplacians. Patterns can develop on the Cartesian network, if they are supported on at least one of its constitutive sub-graphs. Multiplex networks are also obtained under specific prescriptions. In this case, the criteria for the instability reduce to compact explicit formulae. Numerical simulations carried out for the Mimura-Murray reaction kinetics confirm the adequacy of the proposed theory. PMID:26245138
Turing instabilities on Cartesian product networks
NASA Astrophysics Data System (ADS)
Asllani, Malbor; Busiello, Daniel M.; Carletti, Timoteo; Fanelli, Duccio; Planchon, Gwendoline
2015-08-01
The problem of Turing instabilities for a reaction-diffusion system defined on a complex Cartesian product network is considered. To this end we operate in the linear regime and expand the time dependent perturbation on a basis formed by the tensor product of the eigenvectors of the discrete Laplacian operators, associated to each of the individual networks that build the Cartesian product. The dispersion relation which controls the onset of the instability depends on a set of discrete wavelengths, the eigenvalues of the aforementioned Laplacians. Patterns can develop on the Cartesian network, if they are supported on at least one of its constitutive sub-graphs. Multiplex networks are also obtained under specific prescriptions. In this case, the criteria for the instability reduce to compact explicit formulae. Numerical simulations carried out for the Mimura-Murray reaction kinetics confirm the adequacy of the proposed theory.
Turing instabilities on Cartesian product networks.
Asllani, Malbor; Busiello, Daniel M; Carletti, Timoteo; Fanelli, Duccio; Planchon, Gwendoline
2015-01-01
The problem of Turing instabilities for a reaction-diffusion system defined on a complex Cartesian product network is considered. To this end we operate in the linear regime and expand the time dependent perturbation on a basis formed by the tensor product of the eigenvectors of the discrete Laplacian operators, associated to each of the individual networks that build the Cartesian product. The dispersion relation which controls the onset of the instability depends on a set of discrete wavelengths, the eigenvalues of the aforementioned Laplacians. Patterns can develop on the Cartesian network, if they are supported on at least one of its constitutive sub-graphs. Multiplex networks are also obtained under specific prescriptions. In this case, the criteria for the instability reduce to compact explicit formulae. Numerical simulations carried out for the Mimura-Murray reaction kinetics confirm the adequacy of the proposed theory. PMID:26245138